Bayesian model choice and search strategies for mapping interacting quantitative trait Loci

Most complex traits of animals, plants, and humans are influenced by multiple genetic and environmental factors. Interactions among multiple genes play fundamental roles in the genetic control and evolution of complex traits. Statistical modeling of interaction effects in quantitative trait loci (QTL) analysis must accommodate a very large number of potential genetic effects, which presents a major challenge to determining the genetic model with respect to the number of QTL, their positions, and their genetic effects. In this study, we use the methodology of Bayesian model and variable selection to develop strategies for identifying multiple QTL with complex epistatic patterns in experimental designs with two segregating genotypes. Specifically, we develop a reversible jump Markov chain Monte Carlo algorithm to determine the number of QTL and to select main and epistatic effects. With the proposed method, we can jointly infer the genetic model of a complex trait and the associated genetic parameters, including the number, positions, and main and epistatic effects of the identified QTL. Our method can map a large number of QTL with any combination of main and epistatic effects. Utility and flexibility of the method are demonstrated using both simulated data and a real data set. Sensitivity of posterior inference to prior specifications of the number and genetic effects of QTL is investigated.     

Bayesian analysis for genetic architecture of dynamic traits

The dissection of the genetic architecture of quantitative traits, including the number and locations of quantitative trait loci (QTL) and their main and epistatic effects, has been an important topic in current QTL mapping. We extend the Bayesian model selection framework for mapping multiple epistatic QTL affecting continuous traits to dynamic traits in experimental crosses. The extension inherits the efficiency of Bayesian model selection and the flexibility of the Legendre polynomial model fitting to the change in genetic and environmental effects with time. We illustrate the proposed method by simultaneously detecting the main and epistatic QTLs for the growth of leaf age in a doubled-haploid population of rice. The behavior and performance of the method are also shown by computer simulation experiments. The results show that our method can more quickly identify interacting QTLs for dynamic traits in the models with many numbers of genetic effects, enhancing our understanding of genetic architecture for dynamic traits. Our proposed method can be treated as a general form of mapping QTL for continuous quantitative traits, being easier to extend to multiple traits and to a single trait with repeat records.     

Bayesian mapping of genomewide interacting quantitative trait loci for ordinal traits

Development of statistical methods and software for mapping interacting QTL has been the focus of much recent research. We previously developed a Bayesian model selection framework, based on the composite model space approach, for mapping multiple epistatic QTL affecting continuous traits. In this study we extend the composite model space approach to complex ordinal traits in experimental crosses. We jointly model main and epistatic effects of QTL and environmental factors on the basis of the ordinal probit model (also called threshold model) that assumes a latent continuous trait underlies the generation of the ordinal phenotypes through a set of unknown thresholds. A data augmentation approach is developed to jointly generate the latent data and the thresholds. The proposed ordinal probit model, combined with the composite model space framework for continuous traits, offers a convenient way for genomewide interacting QTL analysis of ordinal traits. We illustrate the proposed method by detecting new QTL and epistatic effects for an ordinal trait, dead fetuses, in a F(2) intercross of mice. Utility and flexibility of the method are also demonstrated using a simulated data set. Our method has been implemented in the freely available package R/qtlbim, which greatly facilitates the general usage of the Bayesian methodology for genomewide interacting QTL analysis for continuous, binary, and ordinal traits in experimental crosses.     

Epistatic pleiotropy and the genetic architecture of covariation within early and late-developing skull trait complexes in mice

The role of epistasis as a source of trait variation is well established, but its role as a source of covariation among traits (i.e., as a source of "epistatic pleiotropy") is rarely considered. In this study we examine the relative importance of epistatic pleiotropy in producing covariation within early and late-developing skull trait complexes in a population of mice derived from an intercross of the Large and Small inbred strains. Significant epistasis was found for several pairwise combinations of the 21 quantitative trait loci (QTL) affecting early developing traits and among the 20 QTL affecting late-developing traits. The majority of the epistatic effects were restricted to single traits but epistatic pleiotropy still contributed significantly to covariances. Because of their proportionally larger effects on variances than on covariances, epistatic effects tended to reduce within-group correlations of traits and reduce their overall degree of integration. The expected contributions of single-locus and two-locus epistatic pleiotropic QTL effects to the genetic covariance between traits were analyzed using a two-locus population genetic model. The model demonstrates that, for single-locus or epistatic pleiotropy to contribute to trait covariances in the study population, both traits must show the same pattern of single-locus or epistatic effects. As a result, a large number of the cases where loci show pleiotropic effects do not contribute to the covariance between traits in this population because the loci show a different pattern of effect on the different traits. In general, covariance patterns produced by single-locus and epistatic pleiotropy predicted by the model agreed well with actual values calculated from the QTL analysis. Nearly all single-locus and epistatic pleiotropic effects contributed positive components to covariances between traits, suggesting that genetic integration in the skull is achieved by a complex combination of pleiotropic effects.     

Use of randomization testing to detect multiple epistatic QTLs

Here, we describe a randomization testing strategy for mapping interacting quantitative trait loci (QTLs). In a forward selection strategy, non-interacting QTLs and simultaneously mapped interacting QTL pairs are added to a total genetic model. Simultaneous mapping of epistatic QTLs increases the power of the mapping strategy by allowing detection of interacting QTL pairs where none of the QTL can be detected by their marginal additive and dominance effects. Randomization testing is used to derive empirical significance thresholds for every model selection step in the procedure. A simulation study was used to evaluate the statistical properties of the proposed randomization tests and for which types of epistasis simultaneous mapping of epistatic QTLs adds power. Least squares regression was used for QTL parameter estimation but any other QTL mapping method can be used. A genetic algorithm was used to search for interacting QTL pairs, which makes the proposed strategy feasible for single processor computers. We believe that this method will facilitate the evaluation of the importance at epistatic interaction among QTLs controlling multifactorial traits and disorders.     

Dissecting Genetic Architecture Underlying Seed Traits in Multiple Environments

The seeds of flowering plants develop from double fertilization and play a vital role in reproduction and supplying human and animal food. The genetic variation of seed traits is influenced by multiple genetic systems, e.g., maternal, embryo, and/or endosperm genomes. Understanding the genetic architecture of seed traits is a major challenge because of this complex mechanism of multiple genetic systems, especially the epistasis within or between different genomes and their interactions with the environment. In this study, a statistical model was proposed for mapping QTL with epistasis and QTL-by-environment (QE) interactions underlying endosperm and embryo traits. Our model integrates the maternal and the offspring genomes into one mapping framework and can accurately analyze maternal additive and dominant effects, endosperm/embryo additive and dominant effects, and epistatic effects of two loci in the same or two different genomes, as well as interaction effects of each genetic component of QTL with environment. Intensive simulations under different sampling strategies, heritabilities, and model parameters were performed to investigate the statistical properties of the model. A set of real cottonseed data was analyzed to demonstrate our methods. A software package, QTLNetwork-Seed-1.0.exe, was developed for QTL analysis of seed traits.     

Genome-wide association mapping of agronomic traits in sugar beet

Recent results indicate that association mapping in populations from applied plant breeding is a powerful tool to detect QTL which are of direct relevance for breeding. The focus of this study was to unravel the genetic architecture of six agronomic traits in sugar beet. To this end, we employed an association mapping approach, based on a very large population of 924 elite sugar beet lines from applied plant breeding, fingerprinted with 677 single nucleotide polymorphism (SNP) markers covering the entire genome. We show that in this population linkage disequilibrium decays within a short genetic distance and is sufficient for the detection of QTL with a large effect size. To increase the QTL detection power and the mapping resolution a much higher number of SNPs is required. We found that for QTL detection, the mixed model including only the kinship matrix performed best, even in the presence of a considerable population structure. In genome-wide scans, main effect QTL and epistatic QTL were detected for all six traits. Our full two-dimensional epistasis scan revealed that for complex traits there appear to be epistatic master regulators, loci which are involved in a large number of epistatic interactions throughout the genome.     

Inclusive composite interval mapping (ICIM) for digenic epistasis of quantitative traits in biparental populations

It has long been recognized that epistasis or interactions between non-allelic genes plays an important role in the genetic control and evolution of quantitative traits. However, the detection of epistasis and estimation of epistatic effects are difficult due to the complexity of epistatic patterns, insufficient sample size of mapping populations and lack of efficient statistical methods. Under the assumption of additivity of QTL effects on the phenotype of a trait in interest, the additive effect of a QTL can be completely absorbed by the flanking marker variables, and the epistatic effect between two QTL can be completely absorbed by the four marker-pair multiplication variables between the two pairs of flanking markers. Based on this property, we proposed an inclusive composite interval mapping (ICIM) by simultaneously considering marker variables and marker-pair multiplications in a linear model. Stepwise regression was applied to identify the most significant markers and marker-pair multiplications. Then a two-dimensional scanning (or interval mapping) was conducted to identify QTL with significant digenic epistasis using adjusted phenotypic values based on the best multiple regression model. The adjusted values retain the information of QTL on the two current mapping intervals but exclude the influence of QTL on other intervals and chromosomes. Epistatic QTL can be identified by ICIM, no matter whether the two interacting QTL have any additive effects. Simulated populations and one barley doubled haploids (DH) population were used to demonstrate the efficiency of ICIM in mapping both additive QTL and digenic interactions.     

Varying coefficient models for mapping quantitative trait loci using recombinant inbred intercrosses

There has been a great deal of interest in the development of methodologies to map quantitative trait loci (QTL) using experimental crosses in the last 2 decades. Experimental crosses in animal and plant sciences provide important data sources for mapping QTL through linkage analysis. The Collaborative Cross (CC) is a renewable mouse resource that is generated from eight genetically diverse founder strains to mimic the genetic diversity in humans. The recombinant inbred intercrosses (RIX) generated from CC recombinant inbred (RI) lines share similar genetic structures of F(2) individuals but with up to eight alleles segregating at any one locus. In contrast to F(2) mice, genotypes of RIX can be inferred from the genotypes of their RI parents and can be produced repeatedly. Also, RIX mice typically do not share the same degree of relatedness. This unbalanced genetic relatedness requires careful statistical modeling to avoid false-positive findings. Many quantitative traits are inherently complex with genetic effects varying with other covariates, such as age. For such complex traits, if phenotype data can be collected over a wide range of ages across study subjects, their dynamic genetic patterns can be investigated. Parametric functions, such as sigmoidal or logistic functions, have been used for such purpose. In this article, we propose a flexible nonparametric time-varying coefficient QTL mapping method for RIX data. Our method allows the QTL effects to evolve with time and naturally extends classical parametric QTL mapping methods. We model the varying genetic effects nonparametrically with the B-spline bases. Our model investigates gene-by-time interactions for RIX data in a very flexible nonparametric fashion. Simulation results indicate that the varying coefficient QTL mapping has higher power and mapping precision compared to parametric models when the assumption of constant genetic effects fails. We also apply a modified permutation procedure to control overall significance level.     

The Statistical Power of Inclusive Composite Interval Mapping in Detecting Digenic Epistasis Showing Common F2 Segregation Ratios

Epistasis is a commonly observed genetic phenomenon and an important source of variation of complex traits,which could maintain additive variance and therefore assure the long-term genetic gain in breeding.Inclusive composite interval mapping（ICIM） is able to identify epistatic quantitative trait loci（QTLs） no matter whether the two interacting QTLs have any additive effects.In this article,we conducted a simulation study to evaluate detection power and false discovery rate（FDR） of ICIM epistatic mapping,by considering F₂ and doubled haploid（DH） populations,different F₂ segregation ratios and population sizes.Results indicated that estimations of QTL locations and effects were unbiased,and the detection power of epistatic mapping was largely affected by population size,heritability of epistasis,and the amount and distribution of genetic effects.When the same likelihood of odd（LOD） threshold was used,detection power of QTL was higher in F₂ population than power in DH population;meanwhile FDR in F₂ was also higher than that in DH.The increase of marker density from 10 cM to 5 cM led to similar detection power but higher FDR.In simulated populations,ICIM achieved better mapping results than multiple interval mapping（MIM） in estimation of QTL positions and effect.At the end,we gave epistatic mapping results of ICIM in one actual population in rice（Oryza sativa L.）.     

Statistical methods for dissecting triploid endosperm traits using molecular markers. An autogamous model

The endosperm, a result of double fertilization in flowering plants, is a triploid tissue whose genetic composition is more complex than diploid tissue. We present a new maximum-likelihood-based statistical method for mapping quantitative trait loci (QTL) underlying endosperm traits in an autogamous plant. Genetic mapping of quantitative endosperm traits is qualitatively different from traits for other plant organs because the endosperm displays complicated trisomic inheritance and represents a younger generation than its mother plant. Our endosperm mapping method is based on two different experimental designs: (1) a one-stage design in which marker information is derived from the maternal genome and (2) a two-stage hierarchical design in which marker information is derived from both the maternal and offspring genomes (embryos). Under the one-stage design, the position and additive effect of a putative QTL can be well estimated, but the estimates of the dominant and epistatic effects are upward biased and imprecise. The two-stage hierarchical design, which extracts more genetic information from the material, typically improves the accuracy and precision of the dominant and epistatic effects for an endosperm trait. We discuss the effects on the estimation of QTL parameters of different sampling strategies under the two-stage hierarchical design. Our method will be broadly useful in mapping endosperm traits for many agriculturally important crop plants and also make it possible to study the genetic significance of double fertilization in the evolution of higher plants.     

远交群体动态性状基因定位的似然分析Ⅰ.理论方法

受动物遗传育种中用来估计动态性状育种值的随机回归测定日模型思想的启发 ,将关于时间 (测定日期 )的Legendre多项式镶嵌在遗传模型的每个遗传效应中 ,以刻画QTL对动态性状变化过程的作用 ,从而建立起动态性状基因定位的数学模型。利用远交设计群体 ,阐述了动态性状基因定位的似然分析原理 ,推导了定位参数似然估计的EM法两步求解过程。结合动态性状遗传分析的特点和普通数量性状基因定位研究进展 ,还提出了有关动态性状基因定位进一步研究的设想     

QTL mapping of yield and fiber traits based on a four-way cross population in <Emphasis Type="Italic">Gossypium hirsutum</Emphasis> L.

Four-way cross (4WC) involving four different inbred lines frequently appears in the cotton breeding programs. However, linkage analysis and quantitative trait loci (QTL) mapping with molecular markers in cotton has largely been applied to populations derived from a cross between two inbred lines, and few results of QTL dissection were conducted in a 4WC population. In this study, an attempt was made to construct a linkage map and identify QTL for yield and fiber quality traits in 4WC derived from four different inbred lines in Gossypium hirsutum L. A linkage map was constructed with 285 SSR loci and one morphological locus, covering 2113.3 cM, approximately 42% of the total recombination length of the cotton genome. A total of 31 QTL with 5.1–25.8% of the total phenotypic variance explained were detected. Twenty-four common QTL across environments showed high stability, and six QTL were environment-specific. Several genomic segments affecting multiple traits were identified. The advantage of QTL mapping using a 4WC were discussed. This study presents the first example of QTL mapping using a 4WC population in upland cotton. The results presented here will enhance the understanding of the genetic basis of yield and fiber quality traits and enable further marker-assisted selection in cultivar populations in upland cotton.     

A random model approach to mapping quantitative trait loci for complex binary traits in outbred populations.

Mapping quantitative trait loci (QTL) for complex binary traits is more challenging than for normally distributed traits due to the nonlinear relationship between the observed phenotype and unobservable genetic effects, especially when the mapping population contains multiple outbred families. Because the number of alleles of a QTL depends on the number of founders in an outbred population, it is more appropriate to treat the effect of each allele as a random variable so that a single variance rather than individual allelic effects is estimated and tested. Such a method is called the random model approach. In this study, we develop the random model approach of QTL mapping for binary traits in outbred populations. An EM-algorithm with a Fisher-scoring algorithm embedded in each E-step is adopted here to estimate the genetic variances. A simple Monte Carlo integration technique is used here to calculate the likelihood-ratio test statistic. For the first time we show that QTL of complex binary traits in an outbred population can be scanned along a chromosome for their positions, estimated for their explained variances, and tested for their statistical significance. Application of the method is illustrated using a set of simulated data.     

Mapping quantitative trait loci with epistatic effects

Epistatic variance can be an important source of variation for complex traits. However, detecting epistatic effects is difficult primarily due to insufficient sample sizes and lack of robust statistical methods. In this paper, we develop a Bayesian method to map multiple quantitative trait loci (QTLs) with epistatic effects. The method can map QTLs in complicated mating designs derived from the cross of two inbred lines. In addition to mapping QTLs for quantitative traits, the proposed method can even map genes underlying binary traits such as disease susceptibility using the threshold model. The parameters of interest are various QTL effects, including additive, dominance and epistatic effects of QTLs, the locations of identified QTLs and even the number of QTLs. When the number of QTLs is treated as an unknown parameter, the dimension of the model becomes a variable. This requires the reversible jump Markov chain Monte Carlo algorithm. The utility of the proposed method is demonstrated through analysis of simulation data.     

Simultaneous mapping of epistatic QTL in DU6i × DBA/2 mice

We have mapped epistatic quantitative trait loci (QTL) in an F₂ cross between DU6i × DBA/2 mice. By including these epistatic QTL and their interaction parameters in the genetic model, we were able to increase the genetic variance explained substantially (8.8%–128.3%) for several growth and body composition traits. We used an analysis method based on a simultaneous search for epistatic QTL pairs without assuming that the QTL had any effect individually. We were able to detect several QTL that could not be detected in a search for marginal QTL effects because the epistasis cancelled out the individual effects of the QTL. In total, 23 genomic regions were found to contain QTL affecting one or several of the traits and eight of these QTL did not have significant individual effects. We identified 44 QTL pairs with significant effects on the traits, and, for 28 of the pairs, an epistatic QTL model fit the data significantly better than a model without interactions. The epistatic pairs were classified by the significance of the epistatic parameters in the genetic model, and visual inspection of the two-locus genotype means identified six types of related genotype–phenotype patterns among the pairs. Five of these patterns resembled previously published patterns of QTL interactions.     

Dissection of QTL effects for root traits using a chromosome arm-specific mapping population in bread wheat

A high-resolution chromosome arm-specific mapping population was used in an attempt to locate/detect gene(s)/QTL for different root traits on the short arm of rye chromosome 1 (1RS) in bread wheat. This population consisted of induced homoeologous recombinants of 1RS with 1BS, each originating from a different crossover event and distinct from all other recombinants in the proportions of rye and wheat chromatin present. It provides a simple and powerful approach to detect even small QTL effects using fewer progeny. A promising empirical Bayes method was applied to estimate additive and epistatic effects for all possible marker pairs simultaneously in a single model. This method has an advantage for QTL analysis in minimizing the error variance and detecting interaction effects between loci with no main effect. A total of 15 QTL effects, 6 additive and 9 epistatic, were detected for different traits of root length and root weight in 1RS wheat. Epistatic interactions were further partitioned into inter-genomic (wheat and rye alleles) and intra-genomic (rye–rye or wheat–wheat alleles) interactions affecting various root traits. Four common regions were identified involving all the QTL for root traits. Two regions carried QTL for almost all the root traits and were responsible for all the epistatic interactions. Evidence for inter-genomic interactions is provided. Comparison of mean values supported the QTL detection.     

植物QTL分析的理论研究进展

数量性状的表型是由数量性状基因座 ( Quantitative trait locus,QTL)和环境效应共同作用的结果。传统的数量遗传学采用统计学的方法由一级统计量和二级统计量描述处理 QTL的复合作用 ,估计各种遗传参数 (例如遗传力、遗传相关、遗传进度、有效因子数等 ) ,用于指导遗传育种实践。然而 ,在传统的数量遗传学分析中 ,往往假设数量性状受微效多基因控制 ,这些基因具有相同的并且是较微小的效应 ,所估计的遗传参数反映的是数量性状多基因系统的整体特征 ,其理论方法不能用于追踪研究和描述单个数量性状基因的作用。近年来 ,由于分子生物学技…     

Genetic expectations of quantitative trait loci main and interaction effects obtained with the triple testcross design and their relevance for the analysis of heterosis

Interpretation of experimental results from quantitative trait loci (QTL) mapping studies on the predominant type of gene action can be severely affected by the choice of statistical model, experimental design, and provision of epistasis. In this study, we derive quantitative genetic expectations of (i) QTL effects obtained from one-dimensional genome scans with the triple testcross (TTC) design and (ii) pairwise interactions between marker loci using two-way analyses of variance (ANOVA) under the F(2)- and the F(infinity)-metric model. The theoretical results show that genetic expectations of QTL effects estimated with the TTC design are complex, comprising both main and epistatic effects, and that genetic expectations of two-way marker interactions are not straightforward extensions of effects estimated in one-dimensional scans. We also demonstrate that the TTC design can partially overcome the limitations of the design III in separating QTL main effects and their epistatic interactions in the analysis of heterosis and that dominance x additive epistatic interactions of individual QTL with the genetic background can be estimated with a one-dimensional genome scan. Furthermore, we present genetic expectations of variance components for the analysis of TTC progeny tested in a split-plot design, assuming digenic epistasis and arbitrary linkage.     

An algorithm for molecular dissection of tumor progression

The volumetric growth of tumor cells as a function of time is most often likely to be a complex trait, controlled by the combined influences of multiple genes and environmental influences. Genetic mapping has proven to be a powerful tool for detecting and identifying specific genes affecting complex traits, i.e., quantitative trait loci (QTL), based on polymorphic markers. In this article, we present a novel statistical model for genetic mapping of QTL governing tumor growth trajectories in humans. In principle, this model is a combination of functional mapping proposed to map function-valued traits and linkage disequilibrium mapping designed to provide high resolution mapping of QTL by making use of recombination events created at a historic time. We implement an EM-simplex hybrid algorithm for parameter estimation, in which a closed-form solution for the EM algorithm is derived to estimate the population genetic parameters of QTL including the allele frequencies and the coefficient of linkage disequilibrium, and the simplex algorithm incorporated to estimate the curve parameters describing the dynamic changes of cancer cells for different QTL genotypes. Extensive simulations are performed to investigate the statistical properties of our model. Through a number of hypothesis tests, our model allows for cutting-edge studies aimed to decipher the genetic mechanisms underlying cancer growth, development and differentiation. The implications of our model in gene therapy for cancer research are discussed.