LED and HS have been funded by startup funds from Stanford and a pilot grant to LED from the Stanford Center for Clinical and Translation Research and Education (UL1 TR001085, PI Greenberg). LED has also been funded by Cohen Veterans Bioscience (CVB), and she is part of the CVB Working Group for PTSD Adaptive Platform Trial. BG has been funded by the NIH (R01-HD-059835, PI Williams) and CVB. HH has been funded by the NIH (NIH K01DK114379 and NIH R21AI139012), the Zhengxu and Ying He Foundation, and the Stanley Center for Psychiatric Research. MBR received funds from WPA Congress Mexico City 2018, Guayaquil CEPAM 2019, Asunción X CONGRESO LATINOAMERICANO DE LA FLAPB 2018, Guayaquil 2019 (Bago), and Lancet Psychiatry, London (commission on Violence against women) 2019. SS declares no potential conflict of interest. / Genome-wide approaches including polygenic risk scores (PRSs) are now widely used in medical research; however, few studies have been conducted in low- and middle-income countries (LMICs), especially in South America. This study was designed to test the transferability of psychiatric PRSs to individuals with different ancestral and cultural backgrounds and to provide genome-wide association study (GWAS) results for psychiatric outcomes in this sample. The PrOMIS cohort (N = 3308) was recruited from prenatal care clinics at the Instituto Nacional Materno Perinatal (INMP) in Lima, Peru. Three major psychiatric outcomes (depression, PTSD, and suicidal ideation and/or self-harm) were scored by interviewers using valid Spanish questionnaires. Illumina Multi-Ethnic Global chip was used for genotyping. Standard procedures for PRSs and GWAS were used along with extra steps to rule out confounding due to ancestry. Depression PRSs significantly predicted depression, PTSD, and suicidal ideation/self-harm and explained up to 0.6% of phenotypic variation (minimum p = 3.9 × 10−6). The associations were robust to sensitivity analyses using more homogeneous subgroups of participants and alternative choices of principal components. Successful polygenic prediction of three psychiatric phenotypes in this Peruvian cohort suggests that genetic influences on depression, PTSD, and suicidal ideation/self-harm are at least partially shared across global populations. These PRS and GWAS results from this large Peruvian cohort advance genetic research (and the potential for improved treatments) for diverse global populations. / National Institutes of Health / Revisión por pares

No description available.

Yes / In many species, the offspring of related parents suffer reduced reproductive success, a phenomenon known as inbreeding depression. In humans, the importance of this effect has remained unclear, partly because reproduction between close relatives is both rare and frequently associated with confounding social factors. Here, using genomic inbreeding coefficients (FROH) for >1.4 million individuals, we show that FROH is significantly associated (p < 0.0005) with apparently deleterious changes in 32 out of 100 traits analysed. These changes are associated with runs of homozygosity (ROH), but not with common variant homozygosity, suggesting that genetic variants associated with inbreeding depression are predominantly rare. The effect on fertility is striking: FROH equivalent to the offspring of first cousins is associated with a 55% decrease [95% CI 44–66%] in the odds of having children. Finally, the effects of FROH are confirmed within full-sibling pairs, where the variation in FROH is independent of all environmental confounding.

No description available.

The purpose of this dissertation is to develop new statistical tools for mining big data in plant sciences. In particular, the dissertation consists of four inter-related projects to address various methodological and computational challenges in phylogenetic methods. Project 1 aims to systematically test different optimization tools and provide useful strategies to improve optimization in practice. Project 2 develops a new R package rPlant, which provides a friendly and convenient toolbox for users of iPlant. Project 3 presents a fast and effective group-screening method to identify important genetic factors in GWAS, with theoretical justifications and nice asymptotic properties. Project 4 develops a new statistical tool to identify gene-gene interactions, with the ability of handling the interactions between groups of covariates.

Despite many years of effort, linkage and candidate gene association studies have yielded disappointingly few risk loci for common human diseases such as diabetes, auto-immune disorders and cancers. Large sample sizes, increased understanding of the patterns of correlation in genetic variation, and plunging genotyping costs have enabled genome-wide association studies, which have good power to detect common risk alleles of modest effect. I present an evaluation of SNP choice in study design and show that overall, despite substantial differences in genotyping technologies, marker selection strategies and number of markers assayed, the first generation platforms all offer good levels of genome coverage (∼ 70%). I next describe the largest such project undertaken to date, the Wellcome Trust Case Control Consortium, which consisted of 2000 cases from each of seven common diseases and 3000 shared controls. It identified nearly two dozen new associations. I demonstrate the importance of careful data quality control, including both standard and unorthodox analyses. I next focus on the association results therein for Crohn’s disease. I present a replication experiment in over 1000 additional Crohn’s patients which unambiguously confirmed six previously published loci and four new loci. Next I describe, in a general context, several issues impeding the combination of genome-wide scans, including data annotation, population structure and differences in genotyping platform. Each of these problems is shown to be tractable with available methods, provided that these methods are applied prudently. I present the results of a meta-analysis of three genome-wide scans for Crohn’s disease. The data showed a striking excess of significant associations, and a replication experiment involving over 4000 independent Crohn’s patients verified twenty new risk loci. Finally, I discuss the early success of genome-wide association and its consequences for further understanding the biology of human disease.

Of central importance in the dissection of the components that govern complex traits is understanding the architecture of natural genetic variation. Genetic interaction, or epistasis, constitutes one aspect of this, but epistatic analysis has been largely avoided in genome wide association studies because of statistical and computational difficulties. This thesis explores both issues in the context of two-locus interactions. Initially, through simulation and deterministic calculations it was demonstrated that not only can epistasis maintain deleterious mutations at intermediate frequencies when under selection, but that it may also have a role in the maintenance of additive variance. Based on the epistatic patterns that are evolutionarily persistent, and the frequencies at which they are maintained, it was shown that exhaustive two dimensional search strategies are the most powerful approaches for uncovering both additive variance and the other genetic variance components that are co-precipitated. However, while these simulations demonstrate encouraging statistical benefits, two dimensional searches are often computationally prohibitive, particularly with the marker densities and sample sizes that are typical of genome wide association studies. To address this issue different software implementations were developed to parallelise the two dimensional triangular search grid across various types of high performance computing hardware. Of these, particularly effective was using the massively-multi-core architecture of consumer level graphics cards. While the performance will continue to improve as hardware improves, at the time of testing the speed was 2-3 orders of magnitude faster than CPU based software solutions that are in current use. Not only does this software enable epistatic scans to be performed routinely at minimal cost, but it is now feasible to empirically explore the false discovery rates introduced by the high dimensionality of multiple testing. Through permutation analysis it was shown that the significance threshold for epistatic searches is a function of both marker density and population sample size, and that because of the correlation structure that exists between tests the threshold estimates currently used are overly stringent. Although the relaxed threshold estimates constitute an improvement in the power of two dimensional searches, detection is still most likely limited to relatively large genetic effects. Through direct calculation it was shown that, in contrast to the additive case where the decay of estimated genetic variance was proportional to falling linkage disequilibrium between causal variants and observed markers, for epistasis this decay was exponential. One way to rescue poorly captured causal variants is to parameterise association tests using haplotypes rather than single markers. A novel statistical method that uses a regularised parameter selection procedure on two locus haplotypes was developed, and through extensive simulations it can be shown that it delivers a substantial gain in power over single marker based tests. Ultimately, this thesis seeks to demonstrate that many of the obstacles in epistatic analysis can be ameliorated, and with the current abundance of genomic data gathered by the scientific community direct search may be a viable method to qualify the importance of epistasis.

Breeding has been the mainly responsible for the increase of poultry efficiency in the last decades. The breeding programs are geared towards higher meat yield and feed efficiency. Among the used genomic approaches, genome wide association studies (GWAS) identified quantitative trait loci (QTLs) associated with carcass traits in a meat-type population (TT Reference Population). GWAS analysis identifies variants in linkage disequilibrium with the possible causal mutation and with the aim of refining these results, association study with missense single nucleotide polymorphisms can be useful. A missense SNP can be predicted as deleterious via Sorting Intolerant From Tolerant (SIFT) score when the amino acid change has the potential to impact the protein function and consequently may affects the phenotype. Therefore, in this study, predicted deleterious SNPs within QTLs regions were identified and associated with thigh, drumstick, abdominal fat and breast weight and their yields. Mixed model was used with sex, incubation and SNPs genotypes as fixed effects and family as random effect. From the 20 SNPs analyzed, six were significantly associated (p <0.05) with weight and yield of thigh, breast and drumstick. Three of them s736010549, rs739508259 and rs313532967 are located in the genes WDR77, VWA8 and BARL, respectively. These genes are involved in biological process as steroid hormone signaling pathway, estrogen binding, and regulation of cell proliferation. We determined these genes as candidates for muscle growth. Our strategy allowed the identification of potential causal mutations associated with muscle growth and development. / O melhoramento genético é o principal responsável pelo aumento da eficiência da produção avícola nas últimas décadas e os programas de melhoramento de aves estão direcionados para um maior rendimento de carne e eficiência alimentar. Dentre as abordagens genômicas, estudos de associação genômica ampla (GWAS) identificaram loci associados com características quantitativas (QTLs) de carcaça em uma população de frangos de corte. Análise de GWAS identifica regiões em desequilíbrio de ligação com possíveis mutações causais e com o objetivo de refinar esses resultados, estudos de associações usando polimorfismos de base única (SNPs) não sinônimos podem ser úteis. O SNP não sinônimo pode ser predito como deletério por meio do Sorting Intolerant From Tolerant (SIFT) score quando a alteração de aminoácidos tem o potencial de impactar a função da proteína e consequentemente pode afetar o fenótipo. Portanto, neste estudo, SNPs preditos como deletérios localizados em regiões de QTLs foram identificados e associados com peso e rendimento de coxa, sobrecoxa, gordura abdominal e peito de frangos de corte. Modelo misto foi utilizado, com sexo, incubação e genótipos dos SNPs como efeitos fixos e família como efeito aleatório. De 20 SNPs analisados, seis foram associados significativamente (p<0,05) com peso e rendimento de coxa, sobrecoxa e peito, e três deles rs736010549, rs739508259 e rs313532967 estão presentes nos genes WDR77, VWA8 e BARL, respectivamente. Estes genes estão relacionados com processos biológicos como via de sinalização de esteroide, receptores de estrogênio e de ácidos biliares. Nossa estratégia permitiu a identificação de potenciais mutações causais associadas com crescimento e desenvolvimento muscular.

O transtorno de déficit de atenção e hiperatividade (TDAH) é um dos transtornos psiquiátricos mais freqüentes da infância e adolescência, sendo caracterizado por sintomas de desatenção, hiperatividade e impulsividade. A contribuição genética na etiologia do TDAH é uma das mais altas já verificadas para transtornos psiquiátricos, com herdabilidade média estimada de 76%. Dentre os fatores genéticos que contribuiriam para o desenvolvimento da doença, genes que codificam componentes do sistema dopaminérgico estão entre os principais candidatos. Entre estes, o gene que codifica o receptor D4 de dopamina (DRD4) é o loco mais intensamente investigado nos estudos moleculares com o TDAH. O polimorfismo mais estudado no DRD4 é um VNTR de 48 pb localizado no exon 3; porém outros polimorfismos, localizados na região promotora do gene – uma duplicação de 120pb e os SNPs -521C>T e - 616C>G – também vêm sendo propostos como polimorfismos de suscetibilidade ao TDAH. Além desses, novas variantes em regiões regulatórias do gene, os SNPs rs11246227 e rs11246228, foram observados recentemente em associação com sintomas de desatenção do TDAH. O objetivo geral do presente trabalho foi aumentar a compreensão acerca da participação do gene DRD4 na etiologia do TDAH na nossa população Para tanto, foi testada inicialmente a possibilidade de associação do SNP rs11246227, sendo em seguida investigado o significado funcional dos SNPs rs11246227 e rs11246228, e sua possível relação com a doença, através de ferramentas de bioinformática. O estudo de associação foi realizado em uma amostra composta por 478 pacientes com TDAH, diagnosticados segundo os critérios do DSM-IV, e seus pais biológicos. O rs112466227 foi investigado por abordagens baseada em família (FBAT) e dimensional (PBAT, ANOVA). A possibilidade de desequilíbrio de ligação (DL) com polimorfismos previamente investigados na presente amostra foi estimada pelo programa MLocus. A análise in silico foi realizada utilizando diferentes bases de dados genômicos e programas de predição de sítios alvo para miRNAs e de funcionalidade. A análise pelo FBAT mostrou um desvio significativo da transmissão do alelo C nos pacientes do subtipo desatento. Foram observadas evidências de DL com a duplicação de 120bp e com o VNTR do exon 3. As análises de bioinformática mostraram que os SNPs rs11246227 e rs11246228 estão localizadas na região 3’ do gene DRD4, e não na região 5’, como previamente descrito. Diferenças entre os alelos, com perda ou ganho de sítios de ligação para diferentes miRNAs, foram detectados em ambos os SNPs pelos programas MicroInspector, 5 smiRNAdb e miRecords, e apenas no rs11246227 pelos programas Human miRNA Target e Mirò. A grande variabilidade e a complexidade genética marcante do gene DRD4 aliada à heterogeneidade fenotípica do TDAH provavelmente contribuíram para nossos resultados de associação, divergentes dos descritos na literatura, os quais necessitam de replicação em estudos futuros. Nossos achados em bioinformática sugerem um possível envolvimento dos SNPs investigados com a ligação de miRNAs relacionados aos processos de neurogênese e neuroplasticidade. Genes envolvidos com estes processos vêm sendo identificados nos genome-wide association studies realizados com o TDAH, o que apóia nossos resultados in silico. Entretanto, mais estudos funcionais são necessários, tanto in silico como in vitro, para esclarecer o envolvimento dos polimorfismos analisados na regulação da expressão do gene DRD4 via miRNAs e, consequentemente, do possível efeito desses elementos na etiologia da doença. / Attention-deficit/hyperactivity disorder (ADHD) is one of the most common psychiatric disorders of childhood and adolescence, characterized by inattentive, hyperactive and impulsive symptoms. Genetic contribution to ADHD etiology is one of the highest ever recorded for psychiatric disorders, with a mean heritability of 76%. Among genetic factors that could contribute to disorder development, genes encoding components from dopaminergic system are the main candidate. Of these, the dopamine D4 receptor gene (DRD4) is the most extensively investigated locus in molecular studies of ADHD. The most studied polymorphism in DRD4 gene is a variable number of tandem repeats (VNTR) of 48bp, located at exon 3, although other polymorphisms, located in promoter region – a 120bp duplication and the SNPs -521C> T and-616C> G – have also been proposed as susceptibility polymorphisms for ADHD. In addition, new variants in regulatory regions, the SNPs rs11246227 and rs11246228, have recently been associated with inattentive symptoms of the disorder. The overall objective of this study was to increase the understanding on the involvement of DRD4 gene in ADHD etiology in our population For this purpose, the possibility of association with the SNP rs11246227 was initially tested, being afterwards investigated the functional effect of both rs11246227 and rs11246228 and their possible relation to ADHD through bioinformatics approach. The association study was performed in a sample composed by 478 ADHD patients, diagnosed according to DSM-IV criteria, and their biological parents. The rs112466227 was investigated by both family-based (FBAT) and dimensional (PBAT, ANOVA) approaches. The possibility of linkage disequilibrium (LD) with polymorphisms previously investigated in the present sample was estimated by MLocus software. In silico analysis was conducted using different genomic databases and programs to predict miRNA target sites and functionality. FBAT analysis showed a significant excess of C allele transmission in inattentive subtype patients. Evidences of LD with both 120bp tandem duplication and exon 3 VNTR were observed. Bioinformatics analyses showed that both SNPs rs11246227 and rs11246228 are located in the 3' region of DRD4 gene, and not at 5’ region, as previously described. Differences between alleles, with loss or gain of binding sites, were detected in both SNPs by MicroInspector, smiRNAdb and miRecords, and only in rs11246227 by Human miRNA Targets and miRò DRD4 huge variability and marked genetic complexity allied to ADHD phenotypic heterogeneity might have contributed to our 7 association results, distinct from the ones reported in literature, what needs to be replicated in future studies. Our bioinformatics findings suggest a possible involvement of investigated SNPs in binding properties of miRNAs related to processes of neurogenesis and neuronal plasticity. Genes involved in these processes have been identified in ADHD genome-wide association studies, reinforcing our in silico results. However, new functional studies, using both in silico and in vitro approaches, are needed to clarify the involvement of the investigated polymorphisms in DRD4 expression control mediated by miRNAs and, consequently, the possible effect of these elements in ADHD etiology.

Breeding has been the mainly responsible for the increase of poultry efficiency in the last decades. The breeding programs are geared towards higher meat yield and feed efficiency. Among the used genomic approaches, genome wide association studies (GWAS) identified quantitative trait loci (QTLs) associated with carcass traits in a meat-type population (TT Reference Population). GWAS analysis identifies variants in linkage disequilibrium with the possible causal mutation and with the aim of refining these results, association study with missense single nucleotide polymorphisms can be useful. A missense SNP can be predicted as deleterious via Sorting Intolerant From Tolerant (SIFT) score when the amino acid change has the potential to impact the protein function and consequently may affects the phenotype. Therefore, in this study, predicted deleterious SNPs within QTLs regions were identified and associated with thigh, drumstick, abdominal fat and breast weight and their yields. Mixed model was used with sex, incubation and SNPs genotypes as fixed effects and family as random effect. From the 20 SNPs analyzed, six were significantly associated (p <0.05) with weight and yield of thigh, breast and drumstick. Three of them s736010549, rs739508259 and rs313532967 are located in the genes WDR77, VWA8 and BARL, respectively. These genes are involved in biological process as steroid hormone signaling pathway, estrogen binding, and regulation of cell proliferation. We determined these genes as candidates for muscle growth. Our strategy allowed the identification of potential causal mutations associated with muscle growth and development. / O melhoramento genético é o principal responsável pelo aumento da eficiência da produção avícola nas últimas décadas e os programas de melhoramento de aves estão direcionados para um maior rendimento de carne e eficiência alimentar. Dentre as abordagens genômicas, estudos de associação genômica ampla (GWAS) identificaram loci associados com características quantitativas (QTLs) de carcaça em uma população de frangos de corte. Análise de GWAS identifica regiões em desequilíbrio de ligação com possíveis mutações causais e com o objetivo de refinar esses resultados, estudos de associações usando polimorfismos de base única (SNPs) não sinônimos podem ser úteis. O SNP não sinônimo pode ser predito como deletério por meio do Sorting Intolerant From Tolerant (SIFT) score quando a alteração de aminoácidos tem o potencial de impactar a função da proteína e consequentemente pode afetar o fenótipo. Portanto, neste estudo, SNPs preditos como deletérios localizados em regiões de QTLs foram identificados e associados com peso e rendimento de coxa, sobrecoxa, gordura abdominal e peito de frangos de corte. Modelo misto foi utilizado, com sexo, incubação e genótipos dos SNPs como efeitos fixos e família como efeito aleatório. De 20 SNPs analisados, seis foram associados significativamente (p<0,05) com peso e rendimento de coxa, sobrecoxa e peito, e três deles rs736010549, rs739508259 e rs313532967 estão presentes nos genes WDR77, VWA8 e BARL, respectivamente. Estes genes estão relacionados com processos biológicos como via de sinalização de esteroide, receptores de estrogênio e de ácidos biliares. Nossa estratégia permitiu a identificação de potenciais mutações causais associadas com crescimento e desenvolvimento muscular.