Deficiência intelectual (DI) é um distúrbio frequente do neurodesenvolvimento que afeta ~1% da população mundial e a causa genética é desconhecida em grande parte dos casos. O sequenciamento de exoma é uma ferramenta valiosa na elucidação da etiologia da DI e foi utilizada neste trabalho para investigar 25 famílias com dois ou mais indivíduos afetados. Nosso objetivo foi identificar novos genes ou variantes em genes já conhecidos que pudessem explicar a DI nessas famílias. Identificamos variantes que podem explicar a DI em nove das 25 famílias (36%); três variantes em genes já associados a DI de herança autossômica dominante (SETBP1, MED13L e MBD5), duas em um gene já associado a DI de herança autossômica recessiva (CDK5RAP2 - heterozigose composta) e cinco em genes de herança ligada ao X (UBE2A, MED12, SCL6A8, ARHGAP4 e PHF6). Nenhuma dessas variantes havia sido previamente descrita. O gene ARHGAP4 aparece como um novo candidato e novos estudos serão necessários para estabelecer seu papel na DI. No gene PHF6, a variante identificada mapeia na região 3´UTR e potencialmente interfere com as interações de miRNA reguladores. A variante no gene SETBP1, que apresenta padrão de herança dominante, foi detectada em duas irmãs, indicando que um genitor seria portador obrigatório; essa variante foi detectada em 0,13% das células de sangue periférico da mãe com a utilização da técnica de PCR digital. Essa variante em SETBP1, juntamente com as variantes em MED13L e no gene MBD5, salientam a contribuição de mosaicismo gonadal assim como genitor também com quadro clínico na DI de herança autossômica dominante. A repetição da DI nas irmandades foi essencial para determinar a patogenicidade das variantes identificadas do tipo missense. Dois irmãos portadores de variante missense ainda não descrita no gene UBE2A apresentavam quadro clínico atipicamente leve, trazendo dúvidas sobre a patogenicidade dessa variante; o estudo funcional do gene demonstrou que a mutação de fato prejudicava a função da proteína. Embora estudos funcionais sejam o padrão-ouro para determinar a patogenicidade de uma variante, é difícil sua implementação na rotina diagnóstica / Intellectual deficiency (ID) is a common neurodevelopmental disorder affecting ~1% of the world population, and in which the genetic underlying condition is not determined in many cases. Whole exome sequencing is a valuable tool to elucidate ID etiology and was used in the present work to investigate 25 families with two or more affected individuals. Our main goal was to identify either new genes or variants in already known genes that might explain ID in these families. We detected variants in nine out of the 25 families (36%), with three variants in autossomal dominant ID genes (SETBP1, MED13L e MBD5), two in an autossomal recessive ID gene (compound heterozygote - CDK5RAP2) and five in X-linked ID genes (UBE2A, MED12, SCL6A8, ARHGAP4 e PHF6). None of the variants had been previously described. ARHGAP4 emerged as a new candidate gene and further studies are needed to establish its role in ID. In PHF6, the identified variant mapped to the 3\'UTR region and potentially interferes with miRNA interactions. The variant in SETBP1, which segregates as an autossomal dominant, is present in two sisters, indicating that one parent is an obligate carrier; in fact, the variant was only identified in 0.13% of the mother\'s blood cells after using digital PCR. This case, together with the variants in MED13L and MBD5, highlight the contribution of gonadal mosaicism and affected parent in the autossomal dominant ID. The ID reccurrence in the sibships was essential to determine the putative pathogenicity of missense variants. Two affected brothers carrying a novel missense variant in UBE2A exhibited an atypically mild phenotype, and raised questions regarding the pathogenicity of this variant; associated functional studies showed that the variant impairs gene function. Although functional studies are the gold standard to determine a variant pathogenicity, its implementation as a diagnostic routine is still difficult

A fisiopatologia das síndromes de falência da medula óssea (FMO) está relacionada a mecanismos adquiridos de destruição das células-tronco hematopoeiticas na medula ou a defeitos constitucionais em genes fundamentais para o reparo do DNA e manutenção dos telômeros. A anemia aplástica (AA), o protótipo das doenças de FMO, pode ter etiologia adquirida ou constitucional. A avaliação genética de pacientes com AA adquirida tem como objetivo a detecção de mutações somáticas que possam ser usadas como marcadores de resposta ao tratamento imunossupressor. Diferentemente, em pacientes com AA constitucional, a avaliação genética é fundamental para detecção de mutações etiológicas na doença do paciente, sendo essencial para o tratamento e seleção de doadores de medula óssea. Contudo, o papel das mutações constitucionais na fisiopatologia e modulação imunológica da AA adquirida ainda não é conhecido. Neste estudo, nós sequenciamos pacientes com AA de duas coortes independentes utilizando diferentes painéis de sequenciamento de genes alvos. A primeira coorte, composta por 13 pacientes com AA adquirida, foi sequenciada utilizando um painel com 165 genes relacionados à FMO, neoplasias hematológicas, reparo de DNA, manutenção dos telômeros e vias de resposta imune. A segunda coorte, composta por 59 pacientes investigados para doença constitucional, foi sequenciada com um painel de sequenciamento comercial com 49 genes relacionados à FMO hereditária. Foram identificadas alterações potencialmente patogênicas em três dos cinco pacientes com AA adquirida que não responderam à imunossupressão: dois pacientes com variantes em TERT e um com uma variante em DHX36. Não foram identificadas variantes funcionalmente relevantes nos pacientes que responderam ao tratamento imunossupressor. Em contraste, foram identificadas variantes potencialmente patogênicas em RTEL1 em 8 pacientes com AA constitucional. Variantes em RTEL1 foram associadas tanto ao encurtamento telomérico quanto à erosão excessiva do 3\' overhang, independentemente do comprimento dos telômeros. Desse modo, apenas a medida do comprimento dos telômeros não foi suficiente para identificar todos ospacientes com disfunções teloméricas. As plataformas de sequenciamento de nova geração diminuíram o custo e o tempo para a avaliação genética dos pacientes com FMO. Em nosso estudo, os pacientes com AA adquirida não apresentaram um padrão genético associado à sua resposta ao tratamento com imunossupressores, no entanto, o sequenciamento da coorte com suspeita de AA constitucional foi capaz de identificar o defeito genético associado à doença do paciente em 40% dos casos. O uso de dados clínicos, investigação familiar, análises in silico e testes funcionais foram essenciais para uma correta interpretação da patogenicidade de novas variantes identificadas por sequenciamento de nova geração. / The pathophysiology of bone marrow failure (BMF) can be immune, as in acquired aplastic anemia (AA), or constitutional, due to germline mutations in genes critical for DNA repair and telomere maintenance. The genetic screening of patients with constitutional AA is performed to detect germline mutations that are etiologic in patients\' disease. That is critical for treatment decisions and to identify a donor for a bone marrow transplant. In acquired AA, the genetic screening has been used to detect somatic mutations that can predict patients\' outcomes after treatment, as the role of germline mutations in this disease is yet not clear. To investigate the role of germline variants in AA, we screened two independent cohorts with two different targeting sequencing panels; a first cohort composed by 13 patients with acquired AA that was screened using a panel with 165 genes related to BMF, hematologic malignancies, DNA repair, telomere maintenance, and immune response pathways. A second cohort composed of 59 patients suspected to have a constitutional disease screened by a commercial Inherited Bone Marrow Failure Sequencing panel. In our first cohort, while patients without functional relevant germline variants responded to immunosuppression treatment (n=8), three out of 5 nonresponder patients were identified with variants in telomere biology genes. We found patients carrying TERT and DHX36 variants. In our constitutional AA cohort, we identified 8 patients carrying variants in the RTEL1 gene, a helicase critical to telomere maintenance. RTEL1 variants associated with both patients\' overall telomere shortening and single-stranded 3\' overhang erosion independent of telomere length. Also, 3\' overhang erosion was associated with patients\' predisposition to clonal evolution. In this context, the variants identified in the helicases genes DHX36 and RTEL1 were both associated with patients\' normal telomere length and poor outcomes. Also, telomere length measurement alone was insufficient to identify all primary telomere defects. The platforms of next-generation sequencing decreased the cost and time for the genetic screening of patients with BMF. In our study, acquired AA patients did not display a clear genetic pattern associated with their immunosuppressive treatment response. In contrast, the sequencing of the cohort selected based on their suspicion to have an inherited diseaseidentified a molecular defect that might be pathogenic in up to 40% of patients, including the RTEL1 variants. Pathogenicity assessment of genetic variants requires a combination of clinical, in silico, and functional data required to avoid misinterpretation of common variants.

Dezoito galinhas de uma população experimental utilizada em um cruzamento recíproco entre as linhagens de frangos de corte (TT) e de postura (CC) foram sequenciadas pela tecnologia de nova geração na plataforma Illumina com uma cobertura média de 10X. A descoberta de variantes genéticas foi realizada em uma região de locos de característica quantitativa (Quantitative Trait Locus, QTL), associado anteriormente com peso e percentagem de gordura abdominal no cromossomo 3 da galinha (GGA3), entre os marcadores microssatélites LEI0161 e ADL0371 (33,595,706-42,632,651 pb). O programa SAMtools foi utilizado na identificação de 136.054 SNPs únicos e 15.496 INDELs únicas nos 18 animais sequenciados e após a filtragem das mutações, 92.518 SNPs únicos e 9.298 INDELs únicas foram mantidas. Uma lista de 77 genes foi analisada buscando genes relacionados ao metabolismo de lipídios. Variantes localizadas na região codificante (386 SNPs e 15 INDELs) foram identificadas e associadas com vias metabólicas importantes. Variantes nos genes LOC771163, EGLN1, GNPAT, FAM120B, THBS2 e GGPS1 foram identificadas e podem ser responsáveis pela associação do QTL com a deposição de gordura na carcaça em galinhas. / Eighteen chickens from a parental generation used in a reciprocal cross with broiler and layer lines were sequenced by new generation technology with an average of 10-fold coverage. The DNA sequencing was performed by Illumina next generation platform. The genetic variants discovery was performed in a quantitative trait loci (QTL) region which was previously associated with abdominal fat weight and percentage in chicken chromosome 3 (GGA3) between the microsatellite markers LEI0161 and ADL0371 (33,595,706-42,632,651 bp). SAMtools software was used to detect 136,054 unique SNPs and 15,496 unique INDELs for the 18 chickens, and after quality filtration 92,518 unique SNPs and 9,298 unique INDELs were retained. One list of 77 genes was analised and genes related to lipid metabolism were searched. Variants located in coding region (386 SNPs and 15 INDELs) were identified and associated with important metabolic pathways. Loss of functional variants in the genes LOC771163, EGLN1, GNPAT, FAM120B, THBS2 and GGPS1 may be responsible for the QTL associated with fat deposition in chicken.

A produção brasileira de carne de frango tem uma grande importância econômica no mundo todo devido principalmente aos avanços do melhoramento genético. O surgimento de novas tecnologias de sequenciamento (sequenciamento de nova geração) tem se tornado uma ferramenta poderosa, pois por meio da identificação de SNPs (polimorfismo de nucleotídeo único) e INDELs (deleções/inserções) possibilita a adição de novas informações ao melhoramento genético. A deposição de músculo, em especial o músculo de peito, é uma das características que mais merecem destaque por causa da sua importância nutricional e econômica. Sendo assim o objetivo deste trabalho foi ressequenciar o genoma de 18 aves de duas linhagens distintas experimentais e identificar SNPs e INDELs em uma região de QTL no cromossomo 2 da galinha associado anteriormente com deposição de músculo do peito, além de caracterizar variantes potencialemente funcionais e propor mutações candidatas para estudos futuros. Para isso, dezoito galinhas de duas diferentes linhagens experimentais (corte e postura), ambas desenvolvidas pela Embrapa Suíno e Aves, foram sequenciadas pela plataforma de nova geração da Illumina. SNPs e INDELs foram identificados por meio de ferramentas de bioinformática em uma região de QTL no cromossomo 2 da galinha (105.848.755-112.648.761 pb) que foi previamente associada com deposição de músculo de peito. O sequenciamento dos 18 animais gerou em torno 2,7 bilhões de reads e após a filtragem por qualidade foram mantidas 77% das reads. Em seguida, as reads foram alinhadas ao genoma referência (Gallus_gallus-4.0, NCBI) pela ferramenta Bowtie2 e gerou em média 10,6X de cobertura de sequenciamento na região-alvo. , Foram identificados 722.832 SNPs e 63.727 INDELs para os 18 animais por meio do programa SAMtools, e após uma filtragem rigorosa, foram mantidos 77% dos SNPs (n=558.767) e 60% das INDELs (n=38.402). Com base nas variantes únicas para os 18 animais (85.765 SNPs e 7.824 INDELs) foi realizada a anotação funcional por meio da ferramenta ANNOVAR. Dentre os SNPs não sinônimos (n=153) e stopgain (n=3), 15 foram classificados como deletérios. Um dos SNPs deletérios que já foi depositado em banco de dados foi identificado no gene RB1CC1, que tem sua função relacionada ao desenvolvimento do músculo de peito. Utilizando a ferramenta DAVID foi possível analisar 37 genes relacionados aos SNPs não sinônimos, stopgain, INDELs frameshift e não frameshift. Dentre estes genes, três (DTNA, RB1CC1 e C-MOS) foram selecionados por terem suas funções relacionadas ao desenvolvimento muscular e suas mutações foram analisadas. Sendo assim, futuros estudos podem ser realizados nestes genes candidatos e nas mutações identificadas, por meio de análises de associação e validação em populações comerciais, permitindo assim uma melhor explicação o efeito do QTL estudado. / The Brazilian chicken meat production has a great economic importance in worldwide mainly due to advances in breeding. The emergence of new techniques of sequencing (nextgeneration sequencing) becomes a powerful tool because through identification of SNPs (single nucleotide polymorphism) and INDELs (deletions/insertions) allows the addition of new information for genetic improvement. The muscle deposition, particularly the breast muscle, is one of the features that are most noteworthy because of its nutritional and economic importance. Therefore the aim of this study was to perform the genome resequencing of 18 chicken from two distinct experimental lines and identify SNPs and INDELs in a QTL region on chromosome 2 previously associated with breast muscle, and characterize the variants to identify potentially function ones and propose candidate mutations for future studies. To achieve these objectives, eighteen chickens of two different experimental lines (broiler and layer), both developed by Embrapa Swine and Poultry were sequenced by Illumina next-generation platform. SNPs and INDELs were identified by bioinformatic tools in a QTL region on chicken chromosome 2 (105,848,755-112,648,761 bp) which was previously associated with breast muscle deposition. Sequencing of the eighteen animals generated around 2.7 billion of reads, and 77% of the reads were retained after filtering. The reads were aligned against the chicken genome reference (Gallus_gallus-4.0, NCBI) by Bowtie2 tool resulting in a 10.6X coverage across the target region. Using SAMtools, 722,832 SNPs and 63,727 INDELs were identified in the all individuals, and after a stringent filtration, 77% of SNPs (n=558,767) and 60% of INDELs (n=38,402) were maintained. Based on unique variants for all the animal (85,765 SNPs and 7,828 INDELs) were performed the functional annotation by ANNOVAR tool. Among the non-synonymous SNPs (n=153) and stopgain (n=3), fifteen were predicted like a deleterious mutation. One of deleterious SNPs has already deposited in public database, and it was identified in RB1CC1 gene, which function is related to breast muscle development. Using the DAVID tool was possible to analyze the 37 genes related to the non-synonymous SNPs, stopgain, frameshift and non-frameshift INDELs. Among these genes, three (DTNA, RB1CC1 and C-MOS) were selected due their functions related to muscle development and their mutations were analyzed. Therefore, further association studies can be performed with these candidate genes and their mutations, and also validation in commercial populations, allowing a better explanation of QTL effects.

For over two decades, the study of genetics has been making significant progress towards understanding the causes of common disease. Across a wide range of complex disorders there have been hundreds of associated loci identified, largely driven by common genetic variation. Now, with the advent of next-generation sequencing technology, we are able to interrogate rare and low frequency variation in a high throughput manner for the first time. This provides an exciting opportunity to investigate the role of rarer variation in complex disease risk on a genome-wide scale, potentially o↵ering novel insights into the biological mechanisms underlying disease pathogenesis. In this thesis I will assess the potential of this technology to further our understanding of the genetics of complex disease, using inflammatory bowel disease (IBD) as an example. After first reviewing the history of genetic studies into IBD, I will describe the analytical challenges that can occur when using sequencing to perform case-control association testing at scale, and the methods that can be used to overcome these. I then test for novel IBD associations in a low coverage whole genome sequencing dataset, and uncover a significant burden of rare, damaging missense variation in the gene NOD2, as well as a more general burden of such variation amongst known inflammatory bowel disease risk genes. Through imputation into both new and existing genotyped cohorts, I also describe the discovery of 26 novel IBD-associated loci, including a low frequency missense variant in ADCY7 that approximately doubles the risk of ulcerative colitis. I resolve biological associations underlying several of these novel associations, including a number of signals associated with monocyte-specific changes in integrin gene expression following immune stimulation. These results reveal important insights into the genetic architecture of inflammatory bowel disease, and suggest that a combination of continued array-based genome- wide association studies, imputed using substantial new reference panels, and large scale deep sequencing projects will be required in order to fully understand the genetic basis of complex diseases like IBD.

"Recently, sequencing technologies have generated massive and heterogeneous data sets. However, interpretation of these data sets is a major barrier to understand genomic heterogeneity in complex diseases. In this dissertation, we develop a Bayesian statistical method for single nucleotide level analysis and a global optimization method for gene expression level analysis to characterize genomic heterogeneity in mixed samples. The detection of rare single nucleotide variants (SNVs) is important for understanding genetic heterogeneity using next-generation sequencing (NGS) data. Various computational algorithms have been proposed to detect variants at the single nucleotide level in mixed samples. Yet, the noise inherent in the biological processes involved in NGS technology necessitates the development of statistically accurate methods to identify true rare variants. At the single nucleotide level, we propose a Bayesian probabilistic model and a variational expectation maximization (EM) algorithm to estimate non-reference allele frequency (NRAF) and identify SNVs in heterogeneous cell populations. We demonstrate that our variational EM algorithm has comparable sensitivity and specificity compared with a Markov Chain Monte Carlo (MCMC) sampling inference algorithm, and is more computationally efficient on tests of relatively low coverage (27x and 298x) data. Furthermore, we show that our model with a variational EM inference algorithm has higher specificity than many state-of-the-art algorithms. In an analysis of a directed evolution longitudinal yeast data set, we are able to identify a time-series trend in non-reference allele frequency and detect novel variants that have not yet been reported. Our model also detects the emergence of a beneficial variant earlier than was previously shown, and a pair of concomitant variants. Characterization of heterogeneity in gene expression data is a critical challenge for personalized treatment and drug resistance due to intra-tumor heterogeneity. Mixed membership factorization has become popular for analyzing data sets that have within-sample heterogeneity. In recent years, several algorithms have been developed for mixed membership matrix factorization, but they only guarantee estimates from a local optimum. At the gene expression level, we derive a global optimization (GOP) algorithm that provides a guaranteed epsilon-global optimum for a sparse mixed membership matrix factorization problem for molecular subtype classification. We test the algorithm on simulated data and find the algorithm always bounds the global optimum across random initializations and explores multiple modes efficiently. The GOP algorithm is well-suited for parallel computations in the key optimization steps. "

Next-Generation Sequencing has opened the doors to nearly limitless amounts of genomic data, but the clinical utility of this data is not yet clear. From examining at sequencing data of known familial cancer genes in hereditary cancer patients, the NCGENES study found a clear molecular diagnosis in about 5% of patients and an uncertain molecular result in about 15% of patients. The remaining 80% of hereditary cancer patients received a negative result for the screening of known cancer genes. These latter patients were followed up by whole exome sequencing analysis, and the data was used to perform a research sweep to potentially identify mutation(s) in gene(s) that have yet to be clearly associated with their phenotype. Hereditary breast cancer has a relatively well-established set of susceptibility genes, yet a large percentage of the molecular etiology is still unknown. There are many genes that are good candidates for breast cancer genes based on their protein's function, but they may not actually contribute to breast cancer susceptibility. The ClinGen consortium is aiming to establish the clinical validity of gene-disease associations so that clinicians and patients can better interpret and utilize sequencing results. Six breast cancer susceptibility genes were evaluated using the ClinGen clinical validity framework with the goal of both evaluating the genes already on hereditary breast cancer panels and evaluating genes not yet widely tested to determine if there is enough evidence to support their role in disease to warrant widespread testing. These genes have varying levels of evidence supporting their role in breast cancer susceptibility. The variants in each of the six genes were compared between a cancer patient cohort and a non-cancer patient cohort enrolled in the NCGENES whole exome sequencing study. One likely pathogenic variant and several variants of unknown significance were identified in various genes, and the burden of variants in cancer cases versus controls was evaluated, although the controls were not matched to the cancer cohort in any way. Research sweeps were performed for patients with VUSs to ensure that there were no other mutations in genes that would better fit the phenotype. This thesis presents a method for evaluating gene-disease associations and for utilizing whole exome sequencing data to pinpoint a molecular diagnosis in hereditary breast cancer patients. Overall, it was found that the ClinGen method of evaluating clinical validity of gene-disease associations could be helpful when determining if variants are pathogenic or benign. A new gene, RINT1, was found to have enough evidence to be moderately associated with hereditary breast cancer and it was subsequently added to the diagnostic list so that all cancer patients will now be screened for RINT1 variants. In addition, it was found that two of the genes currently on the diagnostic list, RAD51C and RAD51D, have "disputed" evidence with respect to breast cancer susceptibility. Interestingly, they have much more evidence for an association with ovarian cancer, so if variants are found in these genes, the patient's phenotype should be considered when evaluating them. It was also shown that PALB2, an established breast cancer susceptibility gene, indeed is definitively associated with breast cancer, and the NCGENES cancer patients have more truncating variants than the controls, further validating the clinical validity assertion. Finally, an ovarian cancer patient with two interesting variants, one in SLX4 and one in GEN1, were evaluated. Studies showed that knocking out both of these genes' pathways was highly destructive to the cell. A VUS was found in each of these genes, and it was hypothesized that perhaps these two variants together may be sufficient to contribute to this patient's cancer susceptibility.

The repeated emergence of drug resistance in Plasmodium falciparum underscores the importance of understanding the genetic architecture of current resistance pathways, as well as any associated fitness costs. Why resistance emerges in particular regions of the world has been linked to particular genetic backgrounds that better tolerate resistance-associated polymorphisms; this is likely to play a key role in driving the epidemiology of drug resistance, however is infrequently studied at a large scale in a laboratory setting. The first results chapter establishes a barcoding approach for P. falciparum with the aim of tracking parasite growth in vitro. The strategy used was adapted for P. falciparum by using a pseudogene (PfRh3) as a safe harbour to insert unique molecular barcodes. These libraries of barcoded P. falciparum vectors were also used as a readout of transfection efficiency. The second chapter establishes a proof of principle for phenotyping by barcode sequencing, using a panel of barcoded parasites generated in different genetic backgrounds that comprise sufficient genetic diversity to pilot the method. These were grown in the presence and absence of antimalarial compounds, and growth phenotypes were measured in parallel using BarSeq. The third results chapter studies the contribution of mutations in Pfkelch13, a molecular marker of artemisinin resistance, to parasite fitness. Combining CRISPR/Cas9-based genome editing and high throughput sequencing, the impact of Pfkelch13 alleles on fitness in the context of particular strain backgrounds is revealed. In particular, the impact of genetic background in the emergence and spread of drug-resistant lineages (referred to as KEL1) in Southeast Asia carrying a Y580 Pfkelch13 allele. Overall, given the current pace of genome sequencing of pathogenic organisms such as P. falciparum, it will be important to increase the scale of experimental genetics, in order to tackle in real-time natural variation that might be under constant selection from drugs, thus anticipating the emergence of drug resistance in changing parasite populations. Through this work, tools were developed to facilitate parallel phenotyping by measuring in vitro growth using high-throughput sequencing. The work also develops novel approaches to address the importance of genetic background and a potential role for positive epistasis in a lineage responsible for the recent outbreak of drug-resistant malaria in Southeast Asia.

A produção brasileira de carne de frango tem uma grande importância econômica no mundo todo devido principalmente aos avanços do melhoramento genético. O surgimento de novas tecnologias de sequenciamento (sequenciamento de nova geração) tem se tornado uma ferramenta poderosa, pois por meio da identificação de SNPs (polimorfismo de nucleotídeo único) e INDELs (deleções/inserções) possibilita a adição de novas informações ao melhoramento genético. A deposição de músculo, em especial o músculo de peito, é uma das características que mais merecem destaque por causa da sua importância nutricional e econômica. Sendo assim o objetivo deste trabalho foi ressequenciar o genoma de 18 aves de duas linhagens distintas experimentais e identificar SNPs e INDELs em uma região de QTL no cromossomo 2 da galinha associado anteriormente com deposição de músculo do peito, além de caracterizar variantes potencialemente funcionais e propor mutações candidatas para estudos futuros. Para isso, dezoito galinhas de duas diferentes linhagens experimentais (corte e postura), ambas desenvolvidas pela Embrapa Suíno e Aves, foram sequenciadas pela plataforma de nova geração da Illumina. SNPs e INDELs foram identificados por meio de ferramentas de bioinformática em uma região de QTL no cromossomo 2 da galinha (105.848.755-112.648.761 pb) que foi previamente associada com deposição de músculo de peito. O sequenciamento dos 18 animais gerou em torno 2,7 bilhões de reads e após a filtragem por qualidade foram mantidas 77% das reads. Em seguida, as reads foram alinhadas ao genoma referência (Gallus_gallus-4.0, NCBI) pela ferramenta Bowtie2 e gerou em média 10,6X de cobertura de sequenciamento na região-alvo. , Foram identificados 722.832 SNPs e 63.727 INDELs para os 18 animais por meio do programa SAMtools, e após uma filtragem rigorosa, foram mantidos 77% dos SNPs (n=558.767) e 60% das INDELs (n=38.402). Com base nas variantes únicas para os 18 animais (85.765 SNPs e 7.824 INDELs) foi realizada a anotação funcional por meio da ferramenta ANNOVAR. Dentre os SNPs não sinônimos (n=153) e stopgain (n=3), 15 foram classificados como deletérios. Um dos SNPs deletérios que já foi depositado em banco de dados foi identificado no gene RB1CC1, que tem sua função relacionada ao desenvolvimento do músculo de peito. Utilizando a ferramenta DAVID foi possível analisar 37 genes relacionados aos SNPs não sinônimos, stopgain, INDELs frameshift e não frameshift. Dentre estes genes, três (DTNA, RB1CC1 e C-MOS) foram selecionados por terem suas funções relacionadas ao desenvolvimento muscular e suas mutações foram analisadas. Sendo assim, futuros estudos podem ser realizados nestes genes candidatos e nas mutações identificadas, por meio de análises de associação e validação em populações comerciais, permitindo assim uma melhor explicação o efeito do QTL estudado. / The Brazilian chicken meat production has a great economic importance in worldwide mainly due to advances in breeding. The emergence of new techniques of sequencing (nextgeneration sequencing) becomes a powerful tool because through identification of SNPs (single nucleotide polymorphism) and INDELs (deletions/insertions) allows the addition of new information for genetic improvement. The muscle deposition, particularly the breast muscle, is one of the features that are most noteworthy because of its nutritional and economic importance. Therefore the aim of this study was to perform the genome resequencing of 18 chicken from two distinct experimental lines and identify SNPs and INDELs in a QTL region on chromosome 2 previously associated with breast muscle, and characterize the variants to identify potentially function ones and propose candidate mutations for future studies. To achieve these objectives, eighteen chickens of two different experimental lines (broiler and layer), both developed by Embrapa Swine and Poultry were sequenced by Illumina next-generation platform. SNPs and INDELs were identified by bioinformatic tools in a QTL region on chicken chromosome 2 (105,848,755-112,648,761 bp) which was previously associated with breast muscle deposition. Sequencing of the eighteen animals generated around 2.7 billion of reads, and 77% of the reads were retained after filtering. The reads were aligned against the chicken genome reference (Gallus_gallus-4.0, NCBI) by Bowtie2 tool resulting in a 10.6X coverage across the target region. Using SAMtools, 722,832 SNPs and 63,727 INDELs were identified in the all individuals, and after a stringent filtration, 77% of SNPs (n=558,767) and 60% of INDELs (n=38,402) were maintained. Based on unique variants for all the animal (85,765 SNPs and 7,828 INDELs) were performed the functional annotation by ANNOVAR tool. Among the non-synonymous SNPs (n=153) and stopgain (n=3), fifteen were predicted like a deleterious mutation. One of deleterious SNPs has already deposited in public database, and it was identified in RB1CC1 gene, which function is related to breast muscle development. Using the DAVID tool was possible to analyze the 37 genes related to the non-synonymous SNPs, stopgain, frameshift and non-frameshift INDELs. Among these genes, three (DTNA, RB1CC1 and C-MOS) were selected due their functions related to muscle development and their mutations were analyzed. Therefore, further association studies can be performed with these candidate genes and their mutations, and also validation in commercial populations, allowing a better explanation of QTL effects.

Dentre os vários marcadores genéticos existentes, alguns SNPs (Single Nucleotide Polymorphisms) podem estar associados à determinação de uma série de características fenotípicas (cor de pele, olhos e cabelos, estatura, forma do rosto, espessura do fio de cabelo) tendo sua predição um grande valor nas investigações forenses. Dentre os principais genes conhecidos por controlarem a pigmentação humana, através da produção da melanina, o SLC24A5 (solute carrier family 24, member 5) e o MC1R (melanocortin 1-receptor) apresentam um papel fundamental na melanogênese. O advento do sequenciamento de nova geração permitiu o processamento de várias regiões genômicas e indivíduos simultaneamente, aumentando a disponibilidade e precisão de dados de genomas completos, como alcançado em estudos como o 1000 Genomes Project. Nossas análises preliminares demonstraram que os dados de Fase 3 do 1000 Genomes são muito mais confiáveis do que as versões anteriores. Com esses dados, obtidos para 26 populações globais, foram realizadas análises populacionais (diversidade haplotípica, desequilíbrio de ligação, redes de haplótipo e de variância molecular) para os genes MC1R e SLC24A5 a fim de se compreender melhor seus padrões de diversidade, correlacionando-os à prováveis eventos de seleção natural que tenham moldado sua história evolutiva, como por exemplo a intensidade da radiação UV nas diferentes regiões geográficas. Alguns padrões foram encontradas entre os grupos africano, europeu e asiático, de acordo com a história evolutiva já descrita para estes grupos. Apesar disso, foi observado um padrão claro de varredura seletiva para o SLC24A5 nos resultados obtidos, como o alto desequilíbrio de ligação e baixa diversidade haplotípica. As análises da diversidade do SLC24A5 associadas com as zonas de incidência UV, no entanto, não apresentaram uma correlação clara entre a intensidade da radiação UV e a diversidade haplotípica / Among the various existing genetic markers, some SNPs may be associated with the determination of a series of phenotypic characteristics (skin, eyes and hair color, height, face shape, hair thickness) and its prediction could have a great value in forensic investigations. Among the major genes known to control human pigmentation through melanin production, SLC24A5 (solute carrier family 24, member 5) and MC1R (melanocortin-1 receptor) have a major role in melanogenesis. The advent of next-generation sequencing has enabled processing of several individuals and genomic regions simultaneously while increasing the availability and accuracy of whole genomes data, such as 1000 Genomes Project has achieved. Our preliminary analysis showed that Phase 3 data from the 1000 Genomes are far more reliable than previous versions. Using this data, obtained for 26 global populations, several analyzes were performed (haplotype diversity, linkage disequilibrium, haplotype networks and molecular variance) for MC1R and SLC24A5 in order to better understand their diversity patterns, correlating them to natural selection events which may have shaped their evolutionary history, such as UV radiation intensity in different geographical regions. Some patterns were found between African, European and Asian groups, according to the evolutionary history already described for these groups. Nevertheless, a strong pattern of selective sweep was observed for SLC24A5 in our data, such as high linkage disequilibrium and low haplotype diversity. Analysis of SLC24A5 diversity associated to UV, however, did not show a clear correlation between the UV radiation intensity and haplotype diversity