Vincent Segura

An efficient multi-locus mixed-model approach for genome-wide association studies in structured populations

Population structure causes genome-wide linkage disequilibrium between unlinked loci, leading to statistical confounding in genome-wide association studies. Mixed models have been shown to handle the confounding effects of a diffuse background of large numbers of loci of small effect well, but they do not always account for loci of larger effect. Here we propose a multi-locus mixed model as a general method for mapping complex traits in structured populations. Simulations suggest that our method outperforms existing methods in terms of power as well as false discovery rate. We apply our method to human and Arabidopsis thaliana data, identifying new associations and evidence for allelic heterogeneity. We also show how a priori knowledge from an A. thaliana linkage mapping study can be integrated into our method using a Bayesian approach. Our implementation is computationally efficient, making the analysis of large data sets (n > 10,000) practicable.

A mixed-model approach for genome-wide association studies of correlated traits in structured populations

Genome-wide association studies (GWAS) are a standard approach for studying the genetics of natural variation. A major concern in GWAS is the need to account for the complicated dependence structure of the data, both between loci as well as between individuals. Mixed models have emerged as a general and flexible approach for correcting for population structure in GWAS. Here, we extend this linear mixed-model approach to carry out GWAS of correlated phenotypes, deriving a fully parameterized multi-trait mixed model (MTMM) that considers both the within-trait and between-trait variance components simultaneously for multiple traits. We apply this to data from a human cohort for correlated blood lipid traits from the Northern Finland Birth Cohort 1966 and show greatly increased power to detect pleiotropic loci that affect more than one blood lipid trait. We also apply this approach to an Arabidopsis thaliana data set for flowering measurements in two different locations, identifying loci whose effect depends on the environment.

Genome-Wide Association in Tomato Reveals 44 Candidate Loci for Fruit Metabolic Traits

Genome-wide association sheds light on the genetic architecture of metabolic traits underlying fruit quality in tomato. Genome-wide association studies have been successful in identifying genes involved in polygenic traits and are valuable for crop improvement. Tomato (Solanum lycopersicum) is a major crop and is highly appreciated worldwide for its health value. We used a core collection of 163 tomato accessions composed of S. lycopersicum, S. lycopersicum var cerasiforme, and Solanum pimpinellifolium to map loci controlling variation in fruit metabolites. Fruits were phenotyped for a broad range of metabolites, including amino acids, sugars, and ascorbate. In parallel, the accessions were genotyped with 5,995 single-nucleotide polymorphism markers spread over the whole genome. Genome-wide association analysis was conducted on a large set of metabolic traits that were stable over 2 years using a multilocus mixed model as a general method for mapping complex traits in structured populations and applied to tomato. We detected a total of 44 loci that were significantly associated with a total of 19 traits, including sucrose, ascorbate, malate, and citrate levels. These results not only provide a list of candidate loci to be functionally validated but also a powerful analytical approach for finding genetic variants that can be directly used for crop improvement and deciphering the genetic architecture of complex traits.

GWAPP: A Web Application for Genome-Wide Association Mapping in Arabidopsis

A user-friendly, interactive Web-based application is presented for conducting genome-wide association studies in Arabidopsis. Genome-wide scans for association between phenotype and ∼206,000 single nucleotide polymorphisms in 1386 public accessions can be completed in minutes. The application combines a state-of-the-art mixed model with interactive Manhattan and linkage disequilibrium plots, making it easy to carry out exploratory analyses without programming skills. Arabidopsis thaliana is an important model organism for understanding the genetics and molecular biology of plants. Its highly selfing nature, small size, short generation time, small genome size, and wide geographic distribution make it an ideal model organism for understanding natural variation. Genome-wide association studies (GWAS) have proven a useful technique for identifying genetic loci responsible for natural variation in A. thaliana. Previously genotyped accessions (natural inbred lines) can be grown in replicate under different conditions and phenotyped for different traits. These important features greatly simplify association mapping of traits and allow for systematic dissection of the genetics of natural variation by the entire A. thaliana community. To facilitate this, we present GWAPP, an interactive Web-based application for conducting GWAS in A. thaliana. Using an efficient implementation of a linear mixed model, traits measured for a subset of 1386 publicly available ecotypes can be uploaded and mapped with a mixed model and other methods in just a couple of minutes. GWAPP features an extensive, interactive, and user-friendly interface that includes interactive Manhattan plots and linkage disequilibrium plots. It also facilitates exploratory data analysis by implementing features such as the inclusion of candidate polymorphisms in the model as cofactors.

Phenotyping progenies for complex architectural traits: a strategy for 1-year-old apple trees ( Malus x domestica Borkh.)

The aim of this study was to define a methodology for describing architectural traits in a quantitative way on tree descendants. Our strategy was to collect traits related to both tree structural organization, resulting from growth and branching, and tree form and then to select among these traits relevant descriptors on the basis of their genetic parameters. Because the complexity of tree architecture increases with tree age, we chose to describe the trees in the early stages of development. The study was carried out on a 1-year-old apple progeny derived from two parent cultivars with contrasted architecture. A large number of variables were collected at different positions and scales within the trees. Broad-sense heritability and genetic correlations were estimated and the within tree variability was analyzed for variables measured on long sylleptic axillary shoots (LSAS). These results were combined to select heritable and not correlated variables. Finally, the selection of variables proposed combines topological with geometric traits measured on both trunks and LSAS: (1) on the trunk, mean internode length, and number of sylleptic axillary shoots; (2) on axillary shoots, conicity, bending, and number of sylleptic axillary shoots born at order 3. The trees of the progeny were partitioned on the basis of these variables. The putative agronomic interest of the selected variables with respect to the subsequent tree development is discussed.

Linkage disequilibrium with linkage analysis of multiline crosses reveals different multiallelic QTL for hybrid performance in the flint and dent heterotic groups of maize.

Multiparental designs combined with dense genotyping of parents have been proposed as a way to increase the diversity and resolution of quantitative trait loci (QTL) mapping studies, using methods combining linkage disequilibrium information with linkage analysis (LDLA). Two new nested association mapping designs adapted to European conditions were derived from the complementary dent and flint heterotic groups of maize (Zea mays L.). Ten biparental dent families (N = 841) and 11 biparental flint families (N = 811) were genotyped with 56,110 single nucleotide polymorphism markers and evaluated as test crosses with the central line of the reciprocal design for biomass yield, plant height, and precocity. Alleles at candidate QTL were defined as (i) parental alleles, (ii) haplotypic identity by descent, and (iii) single-marker groupings. Between five and 16 QTL were detected depending on the model, trait, and genetic group considered. In the flint design, a major QTL (R2 = 27%) with pleiotropic effects was detected on chromosome 10, whereas other QTL displayed milder effects (R2 < 10%). On average, the LDLA models detected more QTL but generally explained lower percentages of variance, consistent with the fact that most QTL display complex allelic series. Only 15% of the QTL were common to the two designs. A joint analysis of the two designs detected between 15 and 21 QTL for the five traits. Of these, between 27 for silking date and 41% for tasseling date were significant in both groups. Favorable allelic effects detected in both groups open perspectives for improving biomass production.

Genome-Wide Prediction Methods in Highly Diverse and Heterozygous Species: Proof-of-Concept through Simulation in Grapevine

Nowadays, genome-wide association studies (GWAS) and genomic selection (GS) methods which use genome-wide marker data for phenotype prediction are of much potential interest in plant breeding. However, to our knowledge, no studies have been performed yet on the predictive ability of these methods for structured traits when using training populations with high levels of genetic diversity. Such an example of a highly heterozygous, perennial species is grapevine. The present study compares the accuracy of models based on GWAS or GS alone, or in combination, for predicting simple or complex traits, linked or not with population structure. In order to explore the relevance of these methods in this context, we performed simulations using approx 90,000 SNPs on a population of 3,000 individuals structured into three groups and corresponding to published diversity grapevine data. To estimate the parameters of the prediction models, we defined four training populations of 1,000 individuals, corresponding to these three groups and a core collection. Finally, to estimate the accuracy of the models, we also simulated four breeding populations of 200 individuals. Although prediction accuracy was low when breeding populations were too distant from the training populations, high accuracy levels were obtained using the sole core-collection as training population. The highest prediction accuracy was obtained (up to 0.9) using the combined GWAS-GS model. We thus recommend using the combined prediction model and a core-collection as training population for grapevine breeding or for other important economic crops with the same characteristics.

Comparison of tree architecture using tree edit distances: application to 2-year-old apple hybrids

In fruit trees, understanding genetic determinisms of architectural traits is considered as a promising manner to control vegetative development and yield regularity. In this context, our study aimed to classify 2-year-old apple hybrids on the basis of their architectural traits. From a fine phenotyping, trees were described as tree graphs, including topological and geometric information. To evaluate the similarity between trees, comparison methods based on edit operations (substitution, insertion and deletion) were carried out. Distance between two tree graphs was computed by minimising the sum of the costs of the edit operations applied to transform one tree into another. Two algorithms for the comparison of unordered and partially ordered tree graphs were applied to a sub-sample of the population, taking into account several geometric attributes. For each comparison, a dissimilarity matrix was computed, and subsequently trees were clustered. A local interpretation of the matched entities was proposed through schematic representations of the trees, and similarities between trees were analysed within and between clusters. The tree graphs, both unordered or partially ordered and whether the attributes were considered or not, were grouped, by clustering, according to the number of entities per tree. The robustness of the unordered comparison was demonstrated by its application to the whole population, since it provided results similar to those obtained on the sub-sample. Further developments towards a higher relative weight of geometric versus topological information are discussed in the perspective to define an architectural ideotype in apple.

Association mapping reveals the genetic architecture of tomato response to water deficit: focus on major fruit quality traits

Highlight Tomato quality could be improved under deficit irrigation while maintaining yield. The underlying genetic architecture is polygenic and varies with water availability. Candidate genes related to primary metabolism were identified.

Changes in the epigenome and transcriptome of the poplar shoot apical meristem in response to water availability affect preferentially hormone pathways

The adaptive capacity of long-lived organisms such as trees to the predicted climate changes, including severe and successive drought episodes, will depend on the presence of genetic diversity and phenotypic plasticity. Here, the involvement of epigenetic mechanisms in phenotypic plasticity toward soil water availability was examined in Populus×euramericana. This work aimed at characterizing (i) the transcriptome plasticity, (ii) the genome-wide plasticity of DNA methylation, and (iii) the function of genes affected by a drought-rewatering cycle in the shoot apical meristem. Using microarray chips, differentially expressed genes (DEGs) and differentially methylated regions (DMRs) were identified for each water regime. The rewatering condition was associated with the highest variations of both gene expression and DNA methylation. Changes in methylation were observed particularly in the body of expressed genes and to a lesser extent in transposable elements. Together, DEGs and DMRs were significantly enriched in genes related to phytohormone metabolism or signaling pathways. Altogether, shoot apical meristem responses to changes in water availability involved coordinated variations in DNA methylation, as well as in gene expression, with a specific targeting of genes involved in hormone pathways, a factor that may enable phenotypic plasticity.