Exploring mediterranean durum wheat wild relatives and elite diversity panel for future breeding

Durum wheat (Triticum durum) is an important crop that has been used for millennia for human consumption, and modern breeding can take advantage of the wide variability useful for the adaptation to new challenges. Novel beneficial alleles can be found in wild relatives and landraces thus enhancing crop adaptation to many biotic and abiotic stresses. This dissertation considers the source of variability from both before and after wheat domestication, that caused a loss of potentially useful alleles. Chapter 1. is the thesis introduction, which outlines the importance of wheat in the world, providing an historical overview of the domestication, the evolution mechanisms that led to the current forms of durum wheat and the use of wild relatives as a source of germplasm for future breeding programs is crucial. Moreover, the emergence of Z. tritici has been considered as the main pathogen of wheat since it contains extremely high levels of genetic variability and is thus difficult to control. Chapter 2. Considers the contribution of the phenotypic diversity of 242 accessions of Aegilops tauschii from the Open Wild Wheat Consortium, involved in wheat domestication, provided with whole-genome resequencing. The accessions were phenotyped both in the field and in controlled conditions and A k-mer-based GWAS was performed to identify genomic regions involved in useful traits. Chapter 3. Describes the genetic basis of resistance to Z. tritici in a durum wheat elite diversity panel representative of the germplasm bred in Mediterranean. Quantitative trait loci (QTL) analysis results revealed several loci involved in the STB response that were found in several chromosome regions with a high infection rate. The genomic regions associated with STB resistance identified in this study could be of interest for marker assisted selection (MAS) in durum wheat breeding programs.

Genome-wide association analysis identifies 20 loci that influence adult height.

Adult height is a model polygenic trait, but there has been limited success in identifying the genes underlying its normal variation. To identify genetic variants influencing adult human height, we used genome-wide association data from 13,665 individuals and genotyped 39 variants in an additional 16,482 samples. We identified 20 variants associated with adult height (P < 5 x 10(-7), with 10 reaching P < 1 x 10(-10)). Combined, the 20 SNPs explain approximately 3% of height variation, with a approximately 5 cm difference between the 6.2% of people with 17 or fewer 'tall' alleles compared to the 5.5% with 27 or more 'tall' alleles. The loci we identified implicate genes in Hedgehog signaling (IHH, HHIP, PTCH1), extracellular matrix (EFEMP1, ADAMTSL3, ACAN) and cancer (CDK6, HMGA2, DLEU7) pathways, and provide new insights into human growth and developmental processes. Finally, our results provide insights into the genetic architecture of a classic quantitative trait.

Genome-Wide Mapping of Loci Explaining Variance in Scrotal Circumference in Nellore Cattle

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Power and Robustness of Linkage Tests for Quantitative Traits in General Pedigrees

There are numerous statistical methods for quantitative trait linkage analysis in human studies. An ideal such method would have high power to detect genetic loci contributing to the trait, would be robust to non-normality in the phenotype distribution, would be appropriate for general pedigrees, would allow the incorporation of environmental covariates, and would be appropriate in the presence of selective sampling. We recently described a general framework for quantitative trait linkage analysis, based on generalized estimating equations, for which many current methods are special cases. This procedure is appropriate for general pedigrees and easily accommodates environmental covariates. In this paper, we use computer simulations to investigate the power robustness of a variety of linkage test statistics built upon our general framework. We also propose two novel test statistics that take account of higher moments of the phenotype distribution, in order to accommodate non-normality. These new linkage tests are shown to have high power and to be robust to non-normality. While we have not yet examined the performance of our procedures in the context of selective sampling via computer simulations, the proposed tests satisfy all of the other qualities of an ideal quantitative trait linkage analysis method.

Cost and power analysis of linkage disequilibrium methods for locating genes influencing quantitative traits in humans

Linkage disequilibrium methods can be used to find genes influencing quantitative trait variation in humans. Linkage disequilibrium methods can require smaller sample sizes than linkage equilibrium methods, such as the variance component approach to find loci with a specific effect size. The increase in power is at the expense of requiring more markers to be typed to scan the entire genome. This thesis compares different linkage disequilibrium methods to determine which factors influence the power to detect disequilibrium. The costs of disequilibrium and equilibrium tests were compared to determine whether the savings in phenotyping costs when using disequilibrium methods outweigh the additional genotyping costs.^ Nine linkage disequilibrium tests were examined by simulation. Five tests involve selecting isolated unrelated individuals while four involved the selection of parent child trios (TDT). All nine tests were found to be able to identify disequilibrium with the correct significance level in Hardy-Weinberg populations. Increasing linked genetic variance and trait allele frequency were found to increase the power to detect disequilibrium, while increasing the number of generations and distance between marker and trait loci decreased the power to detect disequilibrium. Discordant sampling was used for several of the tests. It was found that the more stringent the sampling, the greater the power to detect disequilibrium in a sample of given size. The power to detect disequilibrium was not affected by the presence of polygenic effects.^ When the trait locus had more than two trait alleles, the power of the tests maximized to less than one. For the simulation methods used here, when there were more than two-trait alleles there was a probability equal to 1-heterozygosity of the marker locus that both trait alleles were in disequilibrium with the same marker allele, resulting in the marker being uninformative for disequilibrium.^ The five tests using isolated unrelated individuals were found to have excess error rates when there was disequilibrium due to population admixture. Increased error rates also resulted from increased unlinked major gene effects, discordant trait allele frequency, and increased disequilibrium. Polygenic effects did not affect the error rates. The TDT, Transmission Disequilibrium Test, based tests were not liable to any increase in error rates.^ For all sample ascertainment costs, for recent mutations ($<$100 generations) linkage disequilibrium tests were less expensive than the variance component test to carry out. Candidate gene scans saved even more money. The use of recently admixed populations also decreased the cost of performing a linkage disequilibrium test. ^

The genetics of reproductive organ morphology in two Petunia species with contrasting pollination syndromes

Main conclusion Switches between pollination syndromes have happened frequently during angiosperm evolution. Using QTL mapping and reciprocal introgressions, we show that changes in reproductive organ morphology have a simple genetic basis. In animal-pollinated plants, flowers have evolved to optimize pollination efficiency by different pollinator guilds and hence reproductive success. The two Petunia species, P. axillaris and P. exserta, display pollination syndromes adapted to moth or hummingbird pollination. For the floral traits color and scent, genetic loci of large phenotypic effect have been well documented. However, such large-effect loci may be typical for shifts in simple biochemical traits, whereas the evolution of morphological traits may involve multiple mutations of small phenotypic effect. Here, we performed a quantitative trait locus (QTL) analysis of floral morphology, followed by an in-depth study of pistil and stamen morphology and the introgression of individual QTL into reciprocal parental backgrounds. Two QTLs, on chromosomes II and V, are sufficient to explain the interspecific difference in pistil and stamen length. Since most of the difference in organ length is caused by differences in cell number, genes underlying these QTLs are likely to be involved in cell cycle regulation. Interestingly, conservation of the locus on chromosome II in a different P. axillaris subspecies suggests that the evolution of organ elongation was initiated on chromosome II in adaptation to different pollinators. We recently showed that QTLs for pistil and stamen length on chromosome II are tightly linked to QTLs for petal color and volatile emission. Linkage of multiple traits will enable major phenotypic change within a few generations in hybridizing populations. Thus, the genomic architecture of pollination syndromes in Petunia allows for rapid responses to changing pollinator availability.

Localization of non-Mhc collagen-induced arthritis susceptibility loci in DBA/1j mice

One approach to understanding common human diseases is to determine the genetic defects responsible for similar diseases in animal models and place those defective genes in their corresponding biochemical pathways. Our laboratory is working with an animal model for human rheumatoid arthritis called collagen-induced arthritis (CIA). We are particularly interested in determining the location of disease-predisposing loci. To that end, we performed experiments to localize susceptibility loci for CIA in an F2 cross between the highly susceptible mouse strain DBA/1j and the highly resistant mouse strain SWR/j. Specifically, a quantitative trait locus analysis was performed to localize regions of the mouse genome responsible for susceptibility/severity to CIA. One susceptibility locus, Cia1 in the major histocompatibility locus, had been identified previously. Two additional loci were detected in our analysis that contribute to CIA severity (Cia2, Cia3) on chromosomes 2 and 6. A third locus was detected that contributes to the age of onset of the disease. This locus (Cia4) was located on chromosome 2 and was linked to the same region as Cia2. Determining the identity of these loci may provide insights into the etiology of human rheumatoid arthritis.

New strategies to detect and understand genotype-by-environment interactions and QTL-by-environment interactions

Dissertação para obtenção do Grau de Doutor em Estatística e Gestão do Risco, especialidade em Estatística

Rare genomic structural variants in complex disease: lessons from the replication of associations with obesity.

The limited ability of common variants to account for the genetic contribution to complex disease has prompted searches for rare variants of large effect, to partly explain the 'missing heritability'. Analyses of genome-wide genotyping data have identified genomic structural variants (GSVs) as a source of such rare causal variants. Recent studies have reported multiple GSV loci associated with risk of obesity. We attempted to replicate these associations by similar analysis of two familial-obesity case-control cohorts and a population cohort, and detected GSVs at 11 out of 18 loci, at frequencies similar to those previously reported. Based on their reported frequencies and effect sizes (OR≥25), we had sufficient statistical power to detect the large majority (80%) of genuine associations at these loci. However, only one obesity association was replicated. Deletion of a 220 kb region on chromosome 16p11.2 has a carrier population frequency of 2×10(-4) (95% confidence interval [9.6×10(-5)-3.1×10(-4)]); accounts overall for 0.5% [0.19%-0.82%] of severe childhood obesity cases (P = 3.8×10(-10); odds ratio = 25.0 [9.9-60.6]); and results in a mean body mass index (BMI) increase of 5.8 kg.m(-2) [1.8-10.3] in adults from the general population. We also attempted replication using BMI as a quantitative trait in our population cohort; associations with BMI at or near nominal significance were detected at two further loci near KIF2B and within FOXP2, but these did not survive correction for multiple testing. These findings emphasise several issues of importance when conducting rare GSV association, including the need for careful cohort selection and replication strategy, accurate GSV identification, and appropriate correction for multiple testing and/or control of false discovery rate. Moreover, they highlight the potential difficulty in replicating rare CNV associations across different populations. Nevertheless, we show that such studies are potentially valuable for the identification of variants making an appreciable contribution to complex disease.

Methods for testing association between uncertain genotypes and quantitative traits.

Interpretability and power of genome-wide association studies can be increased by imputing unobserved genotypes, using a reference panel of individuals genotyped at higher marker density. For many markers, genotypes cannot be imputed with complete certainty, and the uncertainty needs to be taken into account when testing for association with a given phenotype. In this paper, we compare currently available methods for testing association between uncertain genotypes and quantitative traits. We show that some previously described methods offer poor control of the false-positive rate (FPR), and that satisfactory performance of these methods is obtained only by using ad hoc filtering rules or by using a harsh transformation of the trait under study. We propose new methods that are based on exact maximum likelihood estimation and use a mixture model to accommodate nonnormal trait distributions when necessary. The new methods adequately control the FPR and also have equal or better power compared to all previously described methods. We provide a fast software implementation of all the methods studied here; our new method requires computation time of less than one computer-day for a typical genome-wide scan, with 2.5 M single nucleotide polymorphisms and 5000 individuals.

Discovery and fine mapping of serum protein loci through transethnic meta-analysis.

Many disorders are associated with altered serum protein concentrations, including malnutrition, cancer, and cardiovascular, kidney, and inflammatory diseases. Although these protein concentrations are highly heritable, relatively little is known about their underlying genetic determinants. Through transethnic meta-analysis of European-ancestry and Japanese genome-wide association studies, we identified six loci at genome-wide significance (p < 5 × 10(-8)) for serum albumin (HPN-SCN1B, GCKR-FNDC4, SERPINF2-WDR81, TNFRSF11A-ZCCHC2, FRMD5-WDR76, and RPS11-FCGRT, in up to 53,190 European-ancestry and 9,380 Japanese individuals) and three loci for total protein (TNFRS13B, 6q21.3, and ELL2, in up to 25,539 European-ancestry and 10,168 Japanese individuals). We observed little evidence of heterogeneity in allelic effects at these loci between groups of European and Japanese ancestry but obtained substantial improvements in the resolution of fine mapping of potential causal variants by leveraging transethnic differences in the distribution of linkage disequilibrium. We demonstrated a functional role for the most strongly associated serum albumin locus, HPN, for which Hpn knockout mice manifest low plasma albumin concentrations. Other loci associated with serum albumin harbor genes related to ribosome function, protein translation, and proteasomal degradation, whereas those associated with serum total protein include genes related to immune function. Our results highlight the advantages of transethnic meta-analysis for the discovery and fine mapping of complex trait loci and have provided initial insights into the underlying genetic architecture of serum protein concentrations and their association with human disease.

CUBN is a gene locus for albuminuria.

Identification of genetic risk factors for albuminuria may alter strategies for early prevention of CKD progression, particularly among patients with diabetes. Little is known about the influence of common genetic variants on albuminuria in both general and diabetic populations. We performed a meta-analysis of data from 63,153 individuals of European ancestry with genotype information from genome-wide association studies (CKDGen Consortium) and from a large candidate gene study (CARe Consortium) to identify susceptibility loci for the quantitative trait urinary albumin-to-creatinine ratio (UACR) and the clinical diagnosis microalbuminuria. We identified an association between a missense variant (I2984V) in the CUBN gene, which encodes cubilin, and both UACR (P = 1.1 × 10(-11)) and microalbuminuria (P = 0.001). We observed similar associations among 6981 African Americans in the CARe Consortium. The associations between this variant and both UACR and microalbuminuria were significant in individuals of European ancestry regardless of diabetes status. Finally, this variant associated with a 41% increased risk for the development of persistent microalbuminuria during 20 years of follow-up among 1304 participants with type 1 diabetes in the prospective DCCT/EDIC Study. In summary, we identified a missense CUBN variant that associates with levels of albuminuria in both the general population and in individuals with diabetes.

Genome-wide association analysis of blood-pressure traits in African-ancestry individuals reveals common associated genes in African and non-African populations.

High blood pressure (BP) is more prevalent and contributes to more severe manifestations of cardiovascular disease (CVD) in African Americans than in any other United States ethnic group. Several small African-ancestry (AA) BP genome-wide association studies (GWASs) have been published, but their findings have failed to replicate to date. We report on a large AA BP GWAS meta-analysis that includes 29,378 individuals from 19 discovery cohorts and subsequent replication in additional samples of AA (n = 10,386), European ancestry (EA) (n = 69,395), and East Asian ancestry (n = 19,601). Five loci (EVX1-HOXA, ULK4, RSPO3, PLEKHG1, and SOX6) reached genome-wide significance (p < 1.0 × 10(-8)) for either systolic or diastolic BP in a transethnic meta-analysis after correction for multiple testing. Three of these BP loci (EVX1-HOXA, RSPO3, and PLEKHG1) lack previous associations with BP. We also identified one independent signal in a known BP locus (SOX6) and provide evidence for fine mapping in four additional validated BP loci. We also demonstrate that validated EA BP GWAS loci, considered jointly, show significant effects in AA samples. Consequently, these findings suggest that BP loci might have universal effects across studied populations, demonstrating that multiethnic samples are an essential component in identifying, fine mapping, and understanding their trait variability.

Sex-stratified genome-wide association studies including 270,000 individuals show sexual dimorphism in genetic loci for anthropometric traits.

Given the anthropometric differences between men and women and previous evidence of sex-difference in genetic effects, we conducted a genome-wide search for sexually dimorphic associations with height, weight, body mass index, waist circumference, hip circumference, and waist-to-hip-ratio (133,723 individuals) and took forward 348 SNPs into follow-up (additional 137,052 individuals) in a total of 94 studies. Seven loci displayed significant sex-difference (FDR<5%), including four previously established (near GRB14/COBLL1, LYPLAL1/SLC30A10, VEGFA, ADAMTS9) and three novel anthropometric trait loci (near MAP3K1, HSD17B4, PPARG), all of which were genome-wide significant in women (P<5×10(-8)), but not in men. Sex-differences were apparent only for waist phenotypes, not for height, weight, BMI, or hip circumference. Moreover, we found no evidence for genetic effects with opposite directions in men versus women. The PPARG locus is of specific interest due to its role in diabetes genetics and therapy. Our results demonstrate the value of sex-specific GWAS to unravel the sexually dimorphic genetic underpinning of complex traits.

A genomewide linkage study on suicidality in major depressive disorder confirms evidence for linkage to 2p12.

Suicidal behavior is commonly associated with depression. Twin studies indicate that both suicidality and major depressive disorder (MDD) are heritable. However, epidemiological evidence suggests that the inheritance of suicidality is likely to be independent of the underlying psychiatric disorder, implying a distinct genetic contribution to suicidality. We conducted a genomewide linkage search aiming to detect genomic loci that may harbor susceptibility genes contributing to risk for suicidality in recurrent MDD. Affected sibling pair (ASP) variance components analysis was performed using the Depression Network cohort of 971 ASPs. The quantitative trait measuring suicidality as a broad phenotype, encompassing ideation and suicide attempts, was established from Schedules for Clinical Assessment in Neuropsychiatry interview items. We examined 1,060 genotyped microsatellite markers with an average spacing of 3.3 cM. Empirical thresholds for linkage evidence were set by whole-genome simulations (LOD = 2.71 for genomewide significance, 1.71 for suggestive linkage). No genomewide significant findings were found. Marker D3S1234 on 3p14 achieved suggestive linkage and yielded a maximum LOD of 1.853 (P = 0.0017), loci 9p24.3 and 18q22-q23 achieved LOD scores >1.5. We found some support for linkage to 2p12 (LOD = 1.2, P = 0.0087) which was previously implicated in linkage studies of suicidality. Our follow-up meta-analysis of five studies showed strong linkage to this region (P = 2 × 10(-6) ). In conclusion, this study analyzed suicidality as a continuous trait in MDD. We found modest evidence for linkage on 3p14. Our meta-analysis supports previous evidence of linkage to suicidality on 2p12. Some candidate genes in these regions may plausibly be implicated in suicidality.