Key message
The number of SNPs required for QTL discovery is justified by the distance at which linkage disequilibrium has decayed. Simulations and real potato SNP data showed how to estimate and ...interpret LD decay.
The magnitude of linkage disequilibrium (LD) and its decay with genetic distance determine the resolution of association mapping, and are useful for assessing the desired numbers of SNPs on arrays. To study LD and LD decay in tetraploid potato, we simulated autotetraploid genotypes and used it to explore the dependence on: (1) the number of haplotypes in the population (the amount of genetic variation) and (2) the percentage of haplotype specific SNPs (hs-SNPs). Several estimators for short-range LD were explored, such as the average
r
2
, median
r
2
, and other percentiles of
r
2
(80, 90, and 95 %). For LD decay, we looked at LD
½,90
, the distance at which the short-range LD is halved when using the 90 % percentile of
r
2
at short range, as estimator for LD. Simulations showed that the performance of various estimators for LD decay strongly depended on the number of haplotypes, although the real value of LD decay was not influenced very much by this number. The estimator LD
½,90
was chosen to evaluate LD decay in 537 tetraploid varieties. LD
½,90
values were 1.5 Mb for varieties released before 1945 and 0.6 Mb in varieties released after 2005. LD
½,90
values within three different subpopulations ranged from 0.7 to 0.9 Mb. LD
½,90
was 2.5 Mb for introgressed regions, indicating large haplotype blocks. In pericentromeric heterochromatin, LD decay was negligible. This study demonstrates that several related factors influencing LD decay could be disentangled, that no universal approach can be suggested, and that the estimation of LD decay has to be performed with great care and knowledge of the sampled material.
Genome-wide association studies (GWAS) based on linkage disequilibrium (LD) provide a promising tool for the detection and fine mapping of quantitative trait loci (QTL) underlying complex agronomic ...traits. In this study we explored the genetic basis of variation for the traits heading date, plant height, thousand grain weight, starch content and crude protein content in a diverse collection of 224 spring barleys of worldwide origin. The whole panel was genotyped with a customized oligonucleotide pool assay containing 1536 SNPs using Illumina's GoldenGate technology resulting in 957 successful SNPs covering all chromosomes. The morphological trait "row type" (two-rowed spike vs. six-rowed spike) was used to confirm the high level of selectivity and sensitivity of the approach. This study describes the detection of QTL for the above mentioned agronomic traits by GWAS.
Population structure in the panel was investigated by various methods and six subgroups that are mainly based on their spike morphology and region of origin. We explored the patterns of linkage disequilibrium (LD) among the whole panel for all seven barley chromosomes. Average LD was observed to decay below a critical level (r2-value 0.2) within a map distance of 5-10 cM. Phenotypic variation within the panel was reasonably large for all the traits. The heritabilities calculated for each trait over multi-environment experiments ranged between 0.90-0.95. Different statistical models were tested to control spurious LD caused by population structure and to calculate the P-value of marker-trait associations. Using a mixed linear model with kinship for controlling spurious LD effects, we found a total of 171 significant marker trait associations, which delineate into 107 QTL regions. Across all traits these can be grouped into 57 novel QTL and 50 QTL that are congruent with previously mapped QTL positions.
Our results demonstrate that the described diverse barley panel can be efficiently used for GWAS of various quantitative traits, provided that population structure is appropriately taken into account. The observed significant marker trait associations provide a refined insight into the genetic architecture of important agronomic traits in barley. However, individual QTL account only for a small portion of phenotypic variation, which may be due to insufficient marker coverage and/or the elimination of rare alleles prior to analysis. The fact that the combined SNP effects fall short of explaining the complete phenotypic variance may support the hypothesis that the expression of a quantitative trait is caused by a large number of very small effects that escape detection. Notwithstanding these limitations, the integration of GWAS with biparental linkage mapping and an ever increasing body of genomic sequence information will facilitate the systematic isolation of agronomically important genes and subsequent analysis of their allelic diversity.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Genotype-by-environment interaction (GEI) is an important phenomenon in plant breeding. This paper presents a series of models for describing, exploring, understanding, and predicting GEI. All models ...depart from a two-way table of genotype by environment means. First, a series of descriptive and explorative models/approaches are presented: Finlay-Wilkinson model, AMMI model, GGE biplot. All of these approaches have in common that they merely try to group genotypes and environments and do not use other information than the two-way table of means. Next, factorial regression is introduced as an approach to explicitly introduce genotypic and environmental covariates for describing and explaining GEI. Finally, QTL modeling is presented as a natural extension of factorial regression, where marker information is translated into genetic predictors. Tests for regression coefficients corresponding to these genetic predictors are tests for main effect QTL expression and QTL by environment interaction (QEI). QTL models for which QEI depends on environmental covariables form an interesting model class for predicting GEI for new genotypes and new environments. For realistic modeling of genotypic differences across multiple environments, sophisticated mixed models are necessary to allow for heterogeneity of genetic variances and correlations across environments. The use and interpretation of all models is illustrated by an example data set from the CIMMYT maize breeding program, containing environments differing in drought and nitrogen stress. To help readers to carry out the statistical analyses, GenStat® programs, 15th Edition and Discovery® version, are presented as "Appendix."
High throughput phenotyping (HTP) platforms and devices are increasingly used for the characterization of growth and developmental processes for large sets of plant genotypes. Such HTP data require ...challenging statistical analyses in which longitudinal genetic signals need to be estimated against a background of spatio-temporal noise processes. We propose a two-stage approach for the analysis of such longitudinal HTP data. In a first stage, we correct for design features and spatial trends per time point. In a second stage, we focus on the longitudinal modelling of the spatially corrected data, thereby taking advantage of shared longitudinal features between genotypes and plants within genotypes. We propose a flexible hierarchical three-level P-spline growth curve model, with plants/plots nested in genotypes, and genotypes nested in populations. For selection of genotypes in a plant breeding context, we show how to extract new phenotypes, like growth rates, from the estimated genotypic growth curves and their first-order derivatives. We illustrate our approach on HTP data from the PhenoArch greenhouse platform at INRAE Montpellier and the outdoor Field Phenotyping platform at ETH Zürich.
•We present various genotype to phenotype models to predict complex phenotypes with G×E as function of genotypic and environmental inputs.•We show how genotype-to-phenotype models can be generalized ...to incorporate additional phenotypic information to improve yield predictions.•We discuss how to evaluate the utility of information from new phenotyping techniques in the context of predictive genotype-to-phenotype models.
New types of phenotyping tools generate large amounts of data on many aspects of plant physiology and morphology with high spatial and temporal resolution. These new phenotyping data are potentially useful to improve understanding and prediction of complex traits, like yield, that are characterized by strong environmental context dependencies, i.e., genotype by environment interactions. For an evaluation of the utility of new phenotyping information, we will look at how this information can be incorporated in different classes of genotype-to-phenotype (G2P) models. G2P models predict phenotypic traits as functions of genotypic and environmental inputs. In the last decade, access to high-density single nucleotide polymorphism markers (SNPs) and sequence information has boosted the development of a class of G2P models called genomic prediction models that predict phenotypes from genome wide marker profiles. The challenge now is to build G2P models that incorporate simultaneously extensive genomic information alongside with new phenotypic information. Beyond the modification of existing G2P models, new G2P paradigms are required. We present candidate G2P models for the integration of genomic and new phenotyping information and illustrate their use in examples. Special attention will be given to the modelling of genotype by environment interactions. The G2P models provide a framework for model based phenotyping and the evaluation of the utility of phenotyping information in the context of breeding programs.
A good statistical analysis of genotype × environment interactions (G × E) is a key requirement for progress in any breeding program. Data for G × E analyses traditionally come from multi‐environment ...trials. In recent years, increasingly data are generated from managed stress trials, phenotyping platforms, and high throughput phenotyping techniques in the field. Simultaneously, and complementary to the phenotyping, more elaborate genotyping and envirotyping occur. All of these developments further increase the importance of a sound statistical framework for analyzing G × E. This paper presents considerations on such a framework from the point of view of the choices that need to be made with respect to the content of short academic courses on statistical methods for G × E. Based on our experiences in teaching statistical methods to plant breeders, for specialized G × E courses between three and 5 d are reserved. The audience in such courses includes MSc students, PhD students, postdocs, and researchers at breeding companies. For such specialized courses, we propose a collection of topics to be covered. Our outlook on G × E analyses is two‐fold. On the one hand, we see the G × E problem as the building of predictive models for genotype‐specific reaction norms. On the other hand, the G × E problem consists in the identification of suitable variance‐covariance models to describe heterogeneity of genetic variance and correlations across environments. Our preferred class of statistical models is the class of mixed linear‐bilinear models. These statistical models allow us to answer breeding questions on adaptation, adaptability, stability, and the identification and subdivision of the target population of environments. By a citation analysis of the literature on G × E, we show that our preference for mixed linear‐bilinear models for analyzing G × E is supported by recent trends in the types of methods for G × E analysis that are most frequently cited.
Prediction of the phenotypes for a set of genotypes across multiple environments is a fundamental task in any plant breeding program. Genomic prediction (GP) can assist selection decisions by ...combining incomplete phenotypic information over multiple environments (MEs) with dense sets of markers. We compared a range of ME‐GP models differing in the way environment‐specific genetic effects were modeled. Information among environments was shared either implicitly via the response variable, or by the introduction of explicit environmental covariables. We discuss the models not only in the light of their accuracy, but also in their ability to predict the different parts of the incomplete genotype × environment interaction (G × E) table: (Gt; Et), (Gu; Et), (Gt; Eu), and (Gu; Eu), where G is genotype, E is environment, both tested (t; in one or more instances) and untested (u). Using the ‘Steptoe’ × ‘Morex’ barley (Hordeum vulgare L.) population as an example, we show the advantage of ME‐GP models that account for G × E. In addition, for our example data set, we show that for prediction in the most challenging scenario of untested environments (Eu), the use of explicit environmental information is preferable over the simpler approach of predicting from a main effects model. Besides producing the most general ME‐GP model, the use of environmental covariables naturally links with ecophysiological and crop‐growth models (CGMs) for G × E. We conclude with a list of future research topics in ME‐GP, where we see CGMs playing a central role.
Root phenotyping Kuijken, René C.P.; van Eeuwijk, Fred. A.; Marcelis, Leo F.M. ...
Journal of experimental botany,
09/2015, Letnik:
66, Številka:
18
Journal Article
Recenzirano
In the last decade cheaper and faster sequencing methods have resulted in an enormous increase in genomic data. High throughput genotyping, genotyping by sequencing and genomic breeding are becoming ...a standard in plant breeding. As a result, the collection of phenotypic data is increasingly becoming a limiting factor in plant breeding. Genetic studies on root traits are being hampered by the complexity of these traits and the inaccessibility of the rhizosphere. With an increasing interest in phenotyping, breeders and scientists try to overcome these limitations, resulting in impressive developments in automated phenotyping platforms. Recently, many such platforms have been thoroughly described, yet their efficiency to increase genetic gain often remains undiscussed. This efficiency depends on the heritability of the phenotyped traits as well as the correlation of these traits with agronomically relevant breeding targets. This review provides an overview of the latest developments in root phenotyping and describes the environmental and genetic factors influencing root phenotype and heritability. It also intends to give direction to future phenotyping and breeding strategies for optimizing root system functioning. A quantitative framework to determine the efficiency of phenotyping platforms for genetic gain is described. By increasing heritability, managing effects caused by interactions between genotype and environment and by quantifying the genetic relation between traits phenotyped in platforms and ultimate breeding targets, phenotyping platforms can be utilized to their maximum potential.
Associations between markers and complex quantitative traits were investigated in a collection of 146 modern two-row spring barley cultivars, representing the current commercial germ plasm in Europe. ...Using 236 AFLP markers, associations between markers were found for markers as far apart as 10 cM. Subsequently, for the 146 cultivars the complex traits mean yield, adaptability (Finlay-Wilkinson slope), and stability (deviations from regression) were estimated from the analysis of variety trial data. Regression of those traits on individual marker data disclosed marker-trait associations for mean yield and yield stability. Support for identified associations was obtained from association profiles, i.e., from plots of P-values against chromosome positions. In addition, many of the associated markers were located in regions where earlier QTL were found for yield and yield components. To study the oligogenic genetic base of the traits in more detail, multiple linear regression of the traits on markers was carried out, using stepwise selection. By this procedure, 18-20 markers that accounted for 40-58% of the variation were selected. Our results indicate that association mapping approaches can be a viable alternative to classical QTL approaches based on crosses between inbred lines, especially for complex traits with costly measurements.
Use of whole-genome sequence data (WGS) is expected to improve identification of quantitative trait loci (QTL). However, this requires imputation to WGS, often with a limited number of sequenced ...animals for the target population. The objective of this study was to investigate imputation to WGS in two pig lines using a multi-line reference population and, subsequently, to investigate the effect of using these imputed WGS (iWGS) for GWAS.
Phenotypes and genotypes were available on 12,184 Large White pigs (LW-line) and 4943 Dutch Landrace pigs (DL-line). Imputed 660 K and 80 K genotypes for the LW-line and DL-line, respectively, were imputed to iWGS using Beagle v.4.1. Since only 32 LW-line and 12 DL-line boars were sequenced, 142 animals from eight commercial lines were added. GWAS were performed for each line using the 80 K and 660 K SNPs, the genotype scores of iWGS SNPs that had an imputation accuracy (Beagle R
) higher than 0.6, and the dosage scores of all iWGS SNPs.
For the DL-line (LW-line), imputation of 80 K genotypes to iWGS resulted in an average Beagle R
of 0.39 (0.49). After quality control, 2.5 × 10
(3.5 × 10
) SNPs had a Beagle R
higher than 0.6, resulting in an average Beagle R
of 0.83 (0.93). Compared to the 80 K and 660 K genotypes, using iWGS led to the identification of 48.9 and 64.4% more QTL regions, for the DL-line and LW-line, respectively, and the most significant SNPs in the QTL regions explained a higher proportion of phenotypic variance. Using dosage instead of genotype scores improved the identification of QTL, because the model accounted for uncertainty of imputation, and all SNPs were used in the analysis.
Imputation to WGS using the multi-line reference population resulted in relatively poor imputation, especially when imputing from 80 K (DL-line). In spite of the poor imputation accuracies, using iWGS instead of a lower density SNP chip increased the number of detected QTL and the estimated proportion of phenotypic variance explained by these QTL, especially when dosage scores were used instead of genotype scores. Thus, iWGS, even with poor imputation accuracy, can be used to identify possible interesting regions for fine mapping.