Multiple Trait Analysis of Genetic Mapping for Quantitative Trait Loci

We present in this paper models and statistical methods for performing multiple trait analysis on mapping quantitative trait loci (QTL) based on the composite interval mapping method. By taking into account the correlated structure of multiple traits, this joint analysis has several advantages, compared with separate analyses, for mapping QTL, including the expected improvement on the statistical power of the test for QTL and on the precision of parameter estimation. Also this joint analysis provides formal procedures to test a number of biologically interesting hypotheses concerning the nature of genetic correlations between different traits. Among the testing procedures considered are those for joint mapping, pleiotropy, QTL by environment interaction, and pleiotropy vs. close linkage. The test of pleiotropy (one pleiotropic QTL at a genome position) vs. close linkage (multiple nearby nonpleiotropic QTL) can have important implications for our understanding of the nature of genetic correlations between different traits in certain regions of a genome and also for practical applications in animal and plant breeding because one of the major goals in breeding is to break unfavorable linkage. Results of extensive simulation studies are presented to illustrate various properties of the analyses.     

Quantitative Trait Loci for Metabolic Syndrome in the Hypertension Genetic Epidemiology Network Study

As part of the Hypertension Genetic Epidemiology Network study, genome scans were performed in two ethnicities on the categorical metabolic syndrome (MetS). Genome scans were performed also on the factor scores produced by factor analysis (quantitative MetS). Heritabilities were highest for the obesity‐insulin (INS) factor and lowest for blood pressure (BP) and central obesity. Seventeen unique putative quantitative trait loci (QTLs) yielded logarithm of the odds ratio (LOD) scores in excess of 1.7, 8 for blacks and 9 for whites. Important QTL findings in whites included an LOD score of 3.19 on chromosome 15q15 for the BP factor, 3.08 on chromosome 8p23 for the lipids‐INS factor, and 3.07 on chromosome 3p26 for the obesity‐INS factor. In blacks, after excluding type 2 diabetics, important QTLs were identified, including an LOD score of 2.77 on 13p12 for the obesity‐INS factor and 2.63 on chromosome 11q24 for the lipids‐INS factor. Categorical MetS had lower results than quantitative MetS. Notably, several loci identified overlap with those identified in other studies for a single or group of traits. The most promising candidate loci on 11q24 for lipids‐INS and 13p12 for obesity‐INS in blacks, 8p23 for lipids‐INS, 14q24 for obesity‐INS, and 15q15 for BP in whites warrant further investigation.     

水稻粒长QTL定位与主效基因的遗传分析

该研究利用短粒普通野生稻矮杆突变体和长粒栽培稻品种KJ01组配杂交组合F_1,构建分离群体F_2;并对该群体粒长进行性状遗传分析,利用平均分布于水稻的12条染色体上的132对多态分子标记对该群体进行QTL定位及主效QTLs遗传分析,为进一步克隆新的主效粒长基因奠定基础,并为水稻粒形育种提供理论依据。结果表明:(1)所构建的水稻杂交组合分离群体F_2的粒长性状为多基因控制的数量性状。(2)对543株F_2分离群体进行QTL连锁分析,构建了控制水稻粒长的连锁遗传图谱,总长为1 713.94 cM,共检测出24个QTLs,只有3个表现为加性遗传效应,其余位点均表现为遗传负效应。(3)检测到的3个主效QTLs分别位于3号染色体的分子标记PSM379～RID24455、RID24455～RM15689和RM571～RM16238之间,且三者对表型的贡献率分别为54.85%、31.02%和7.62%。(4)在标记PSM379～RID24455之间已克隆到的粒长基因为该研究新发现的主效QTL位点。     

Mapping Environment-Specific Quantitative Trait Loci

Environment-specific quantitative trait loci (QTL) refer to QTL that express differently in different environments, a phenomenon called QTL-by-environment (Q × E) interaction. Q × E interaction is a difficult problem extended from traditional QTL mapping. The mixture model maximum-likelihood method is commonly adopted for interval mapping of QTL, but the method is not optimal in handling QTL interacting with environments. We partitioned QTL effects into main and interaction effects. The main effects are represented by the means of QTL effects in all environments and the interaction effects are represented by the variances of the QTL effects across environments. We used the Markov chain Monte Carlo (MCMC) implemented Bayesian method to estimate both the main and the interaction effects. The residual error covariance matrix was modeled using the factor analytic covariance structure. A simulation study showed that the factor analytic structure is robust and can handle other structures as special cases. The method was also applied to Q × E interaction mapping for the yield trait of barley. Eight markers showed significant main effects and 18 markers showed significant Q × E interaction. The 18 interacting markers were distributed across all seven chromosomes of the entire genome. Only 1 marker had both the main and the Q × E interaction effects. Each of the other markers had either a main effect or a Q × E interaction effect but not both.GENOTYPE-BY-ENVIRONMENT (G × E) interaction is a very important phenomenon in quantitative genetics. With the advanced molecular technology and statistical methods for quantitative trait loci (QTL) mapping (Lander and Botstein 1989; Jansen 1993; Zeng 1994), G × E interaction analysis has shifted to QTL-by-environment (Q × E) interaction. In the early stage of QTL mapping, almost all statistical methods were developed in a single environment (Paterson et al. 1991; Stuber et al. 1992). Data from different environments were analyzed separately and the conclusions were drawn from the separate analyses of QTL across environments. These methods do not consider the correlation of data under different environments and thus may not extract maximum information from the data. Composite interval mapping for multiple traits can be used for Q × E interaction if different traits are treated as the same trait measured in different environments (Jiang and Zeng 1995). This multivariate composite interval mapping approach makes good use of all data simultaneously and increases statistical power of QTL detection and accuracy of the estimated QTL positions. However, the number of parameters of this method increases dramatically as the number of environments increases. Therefore, the method may not be applied when the number of environments is large. Several other models have been proposed to solve the problem of a large number of environments (Jansen et al. 1995; Beavis and Keim 1996; Romagosa et al. 1996). These methods were based on some special situations and assumptions. One typical assumption was independent errors or constant variances across environments. These assumptions are often violated in real QTL mapping experiments.Earlier investigators realized the problem and adopted the mixed-model methodology to solve the problem (Piepho 2000; Boer et al. 2007). Under the mixed-model framework, people can choose which model effects are random and which are fixed. The mixed-model methodology is very flexible, leading to an easy way to model genetic and environmental correlation between environments using a suitable error structure. Piepho (2000) proposed a mixed model to detect QTL main effect across environments. Similar to the composite interval mapping analysis, his model incorporated one putative QTL and a few cofactors. The Q × E effects in the model were assumed to be random, which greatly reduced the number of estimated parameters. However, the fact that only one QTL is included in the model means that Piepho''s (2000) model remains a single-QTL model rather than a multivariate model. Boer et al. (2007) proposed a step-by-step mixed-model approach to detecting QTL main effects, Q × E interaction effects, and QTL responses to specific environmental covariates. In the final step, Boer et al. (2007) rewrote the model to include all QTL in a multiple-QTL model and reestimated their effects.In this study, we extended the Bayesian shrinkage method (Xu 2003) to map Q × E interaction effects of QTL. In the original study (Xu 2003), we treated each marker as a putative QTL and used the shrinkage method to simultaneously estimate marker effects of the entire genome. In the multiple-environment case, we can still use this approach to simultaneously evaluate marker effects under multiple environments but we can further partition the marker effects into main and Q × E interaction effects. For any particular marker, the mean of the marker effects represents the main effect and the variance of the marker effects represents the Q × E interaction effect for that marker. Under the Bayesian framework, we assigned a normal prior with zero mean and an unknown variance to each marker main effect and a multivariate normal prior with zero vector mean and homogeneous diagonal variance–covariance matrix to the Q × E interaction effects of each maker. In multiple environments, the structure of the error terms might be very complicated since we need to consider the correlation of the same genotype under different environments. In our analysis, we used different variance–covariance structures to model the error terms. The simplest case was the homogeneous diagonal matrix, and the most complex choice was an unstructured matrix. We also used a heterogeneous diagonal matrix whose parameters are somewhere between the two models. Finally, we considered several factor analytic models. The reason to use the factor analytic structure is that it can separate genetic effects into common effects and environment-specific effects. In addition, the factor analytic structure is parsimonious and thus can substantially reduce the computational burden of the mixed-model analyses.     

Association between Expression Quantitative Trait Loci and Metabolic Traits in Two Korean Populations

Most genome-wide association studies consider genes that are located closest to single nucleotide polymorphisms (SNPs) that are highly significant for those studies. However, the significance of the associations between SNPs and candidate genes has not been fully determined. An alternative approach that used SNPs in expression quantitative trait loci (eQTL) was reported previously for Crohn’s disease; it was shown that eQTL-based preselection for follow-up studies was a useful approach for identifying risk loci from the results of moderately sized GWAS. In this study, we propose an approach that uses eQTL SNPs to support the functional relationships between an SNP and a candidate gene in a genome-wide association study. The genome-wide SNP genotypes and 10 biochemical measures (fasting glucose levels, BUN, serum albumin levels, AST, ALT, gamma GTP, total cholesterol, HDL cholesterol, triglycerides, and LDL cholesterol) were obtained from the Korean Association Resource (KARE) consortium. The eQTL SNPs were isolated from the SNP dataset based on the RegulomeDB eQTL-SNP data from the ENCODE projects and two recent eQTL reports. A total of 25,658 eQTL SNPs were tested for their association with the 10 metabolic traits in 2 Korean populations (Ansung and Ansan). The proportion of phenotypic variance explained by eQTL and non-eQTL SNPs showed that eQTL SNPs were more likely to be associated with the metabolic traits genetically compared with non-eQTL SNPs. Finally, via a meta-analysis of the two Korean populations, we identified 14 eQTL SNPs that were significantly associated with metabolic traits. These results suggest that our approach can be expanded to other genome-wide association studies.     

Mapping multiple Quantitative Trait Loci by Bayesian classification

We developed a classification approach to multiple quantitative trait loci (QTL) mapping built upon a Bayesian framework that incorporates the important prior information that most genotypic markers are not cotransmitted with a QTL or their QTL effects are negligible. The genetic effect of each marker is modeled using a three-component mixture prior with a class for markers having negligible effects and separate classes for markers having positive or negative effects on the trait. The posterior probability of a marker's classification provides a natural statistic for evaluating credibility of identified QTL. This approach performs well, especially with a large number of markers but a relatively small sample size. A heat map to visualize the results is proposed so as to allow investigators to be more or less conservative when identifying QTL. We validated the method using a well-characterized data set for barley heading values from the North American Barley Genome Mapping Project. Application of the method to a new data set revealed sex-specific QTL underlying differences in glucose-6-phosphate dehydrogenase enzyme activity between two Drosophila species. A simulation study demonstrated the power of this approach across levels of trait heritability and when marker data were sparse.     

水稻QTL定位研究进展

水稻的许多重要农艺性状均属于数量性状,研究水稻数量性状遗传对水稻育种具有十分重要的意义.近年来大量的研究揭示了水稻QTL的基本特征,剖析了重要农艺性状的遗传基础,给水稻遗传改良带来了新策略,不断深入的研究已经完成了水稻特定数量基因的精细定位和克隆,到目前为止已经有一万多个水稻QTL进行了定位,其中有19个进行了克隆,这对水稻育种具有十分重要的意义.本文主要综述了QTL定位的理论基础,水稻QTL定位的研究进展,并对水稻QTL研究的趋势进行了展望.     

Quantitative Trait Loci Identification by Estimating the Genetic Model based on the Extremal Samples

Background In genetic association studies with quantitative trait loci (QTL), the association between a candidate genetic marker and the trait of interest is commonly examined by the omnibus F test or by the t-test corresponding to a given genetic model or mode of inheritance. It is known that the t-test with a correct model specification is more powerful than the F test. However, since the underlying genetic model is rarely known in practice, the use of a model-specific t-test may incur substantial power loss. Robust-efficient tests, such as the Maximin Efficiency Robust Test (MERT) and MAX3 have been proposed in the literature.Methods In this paper, we propose a novel two-step robust-efficient approach, namely, the genetic model selection (GMS) method for quantitative trait analysis. GMS selects a genetic model by testing Hardy-Weinberg disequilibrium (HWD) with extremal samples of the population in the first step and then applies the corresponding genetic model-specific t-test in the second step.Results Simulations show that GMS is not only more efficient than MERT and MAX3, but also has comparable power to the optimal t-test when the genetic model is known.Conclusion Application to the data from Alzheimer’s Disease Neuroimaging Initiative (ADNI) cohort demonstrates that the proposed approach can identify meaningful biological SNPs on chromosome 19.     

Marker-Assisted Introgression of Quantitative Trait Loci

The use of molecular markers for the introgression of one or several superior QTL alleles into a recipient line is investigated using analytic and simulation results. The positions of the markers devoted to the control of the genotype at the QTLs in a ``foreground selection' step are optimized given the confidence interval of the QTL position. Results demonstrate that using at least three markers per QTL allows a good control over several generations. Population sizes that should be recommended for various numbers of QTLs are calculated and are used to determine the limit in the number of QTLs that can be monitored simultaneously. If ``background selection' devoted to accelerate the return to the recipient parent genotype outside the QTL regions is applied, the positions of the markers devoted to the control of the QTLs have to be reconsidered. When several QTLs are monitored simultaneously, background selection among the limited number of individuals resulting from the foreground selection step accelerates the increase in genomic similarity with the recipient parent, with only limited costs. Background selection is even more efficient in a pyramidal backcross program where QTLs are first monitored one by one.     

作物数量性状基因研究进展

分子生物技术的发展对作物数量性状基因(QTL)研究提供了条件,不同的定位群体各有其特点,相继出现的QTL定位方法也逐步完善.大量的研究揭示了QTL的基本特征,剖析了重要农艺性状的遗传基础,给作物遗传改良带来了新的策略,不断深入的研究已经完成了特定QTL的精细定位和克隆.本文从QTL的定位群体,定位方法,研究现状,精细定位与克隆,以及QTL利用等方面对作物数量性状基因的研究进行了综述。 Abstract:With the rapid development of molecular biotechnology,QTL analyses were executed for a lot of important agronomic traits in many crops.Different experimental populations and mapping methods had their own advantages in QTL analysis.Amounts of studies paid attention to locate the QTLs for important traits,and others tried to disect the genetic bases using molecular markers.Near isogenic lines were the best populations for QTL fine mapping and positional cloning,A few studies had been reported their results on materials with improvement traits using marker-assisted selection.This paper summarizes the recent progress on QTL mapping populations and methods,the status of QTL locating,QTL fine mapping and positional cloning,and QTL.application in breeding.     

作物数量性状基因研究进展

分子生物技术的发展对作物数量性状基因（QTL）研究提供了条件,不同的定位群体各有其特点,相继出现的QTL定位也逐步完善。大量的研究揭示了QTL的基本特征,剖析了重要农艺4性状的遗传基础,给作物遗传改良带来了新的策略,不断深入的研究已经完成了特定的QTL的精细定位和克隆。本从QTL的定位群体,定位方法,研究现状,精细定位与克隆,以及QTL利用等方面对作物数量性状基因的研究进行了综述。     

Precision Mapping of Quantitative Trait Loci

Adequate separation of effects of possible multiple linked quantitative trait loci (QTLs) on mapping QTLs is the key to increasing the precision of QTL mapping. A new method of QTL mapping is proposed and analyzed in this paper by combining interval mapping with multiple regression. The basis of the proposed method is an interval test in which the test statistic on a marker interval is made to be unaffected by QTLs located outside a defined interval. This is achieved by fitting other genetic markers in the statistical model as a control when performing interval mapping. Compared with the current QTL mapping method (i.e., the interval mapping method which uses a pair or two pairs of markers for mapping QTLs), this method has several advantages. (1) By confining the test to one region at a time, it reduces a multiple dimensional search problem (for multiple QTLs) to a one dimensional search problem. (2) By conditioning linked markers in the test, the sensitivity of the test statistic to the position of individual QTLs is increased, and the precision of QTL mapping can be improved. (3) By selectively and simultaneously using other markers in the analysis, the efficiency of QTL mapping can be also improved. The behavior of the test statistic under the null hypothesis and appropriate critical value of the test statistic for an overall test in a genome are discussed and analyzed. A simulation study of QTL mapping is also presented which illustrates the utility, properties, advantages and disadvantages of the method.     

Quantitative Trait Loci for Murine Growth

Body size is an archetypal quantitative trait with variation due to the segregation of many gene loci, each of relatively minor effect, and the environment. We examine the effects of quantitative trait loci (QTLs) on age-specific body weights and growth in the F(2) intercross of the LG/J and SM/J strains of inbred mice. Weekly weights (1-10 wk) and 75 microsatellite genotypes were obtained for 535 mice. Interval mapping was used to locate and measure the genotypic effects of QTLs on body weight and growth. QTL effects were detected on 16 of the 19 autosomes with several chromosomes carrying more than one QTL. The number of QTLs for age-specific weights varied from seven at 1 week to 17 at 10 wk. The QTLs were each of relatively minor, subequal effect. QTLs affecting early and late growth were generally distinct, mapping to different chromosomal locations indicating separate genetic and physiological systems for early and later murine growth.     

作物数量性状(QTL)基因研究进展

从作物数量性状基因座QTL(quantitative trait locus)作图群体类型及特点,QTL定位的原理和方法,作物QTL研究现状,以及QTL精细定位、克隆、利用等方面进行了综述。对作物QTL分子标记辅助选择育种进行了探讨,并对目前QTL定位中存在的问题和今后QTL的研究方向提出了一些思考。     

构建数量性状基因图谱的几种统计方法

以线性数学模型为线索,概述了用于构建数量性状基因图谱的几种主要统计方法,包括方差分析法、标记回归法、区间作图法、复合区间作图法、Jansen的复合区间作图法、双侧标记回归法以及新近发展的多区间作图法和多亲本作图法等.讨论了各种方法的优缺点. Abstract:Statistical methods for mapping QTLs were summarized, including one marker analysis, arker regression analysis,interval mapping (IM),composite interval mapping (CIM),Jansen's composite interval mapping, flanking marker regression analysis,multiple interval mapping (MIM) and multiple families mapping.Their advantages and disadvantages were discussed.     

苹果叶片相关性状的QTL定位及其遗传效应分析

利用苹果栽培品种‘红富士’和新疆野苹果优系‘红肉苹果’杂交的110个F1株系为作图群体,构建了苹果的分子遗传图谱,采用区间作图法对苹果9个叶片相关性状(叶片长度、叶片宽度、叶片厚度、叶柄长度、叶片面积、总叶绿素含量、叶绿素a含量、叶绿素b含量和类胡萝卜素含量)进行了QTL定位分析。结果显示:从110个F1株系中共检测到20个控制叶片相关性状的QTL位点,分布在第1、2、3、4、5、7、8、10、11、12、16、17连锁群上;各QTL位点的LOD值在2.58~3.55之间,其中主效QTL位点2个(LOD≥3.5),可解释11.63%~16.36%的表型变异。获得紧密连锁的特异标记(CH05d11-435m、CH04c06-201m)为进一步进行QTL精细定位提供了参考。     

多个连锁数量性状位点的多性状最优化选择

一种扩展的方法能够在多个世代对具有多个数量性状位点的多性状选择进行最优化。这种方法的基础是在目标雨数中用综合遗传值替代单个性状遗传值,并在整个规划期内最大化所有世代选择反应的加权和。利用多阶段系统优化控制理论,整个最优化问题通过一个向前和向后的迭代循环解决。用一个实际育种猪群的育种参数来评价该方法的选择效果,并和标准QTL选择和常规BLUP选择进行比较。结果表明,优化选择要优于标准QTL选择和常规BLUP选择。经济权重对优化选择的影响较明显,随着达100kg日龄赋予的经济权重的增加,优化选择的优势越明显。优化选择通过两种方式增加总选择反应：1）选择早期减少一部分QTL选择反应;2）对达100kgH龄给予更大的权重。选择后期优化累积贴现选择比优化终端选择给予达100kgH龄更大的权重。     

Mapping and Analysis of Dairy Cattle Quantitative Trait Loci by Maximum Likelihood Methodology Using Milk Protein Genes as Genetic Markers

Maximum likelihood methodology was used to estimate effects of both a marker gene and a linked quantitative trait locus (QTL) on quantitative traits in a segregating population. Two alleles were assumed for the QTL. In addition to the effects of genotypes at both loci on the mean of the quantitative trait, recombination frequency between the loci, frequency of the QTL alleles and the residual standard deviation were also estimated. Thus six parameters were estimated in addition to the marker genotype means. The statistical model was tested on simulated data, and used to estimate direct and linked effects of the milk protein genes, β-lactoglobulin, κcasein, and β-casein, on milk, fat, and protein production and fat and protein percent in the Dutch dairy cattle population. β-Lactoglobulin had significant direct effects on milk yield and fat percent. κ-Casein had significant direct effects on milk yield, protein percent and fat yield. β-Casein had significant direct effects on milk yield, fat and protein percent and fat and protein yield. Linked QTL with significant effects on fat percent were found for κ-casein and β-casein. Since the β-casein and κ-casein genes are closely linked, it is likely that the same QTL was detected for those two markers. Further, a QTL with a significant effect on fat yield was found to be linked to κ-casein and a QTL with a significant effect on protein yield was linked to β-lactoglobulin.     

Interval Mapping of Quantitative Trait Loci Employing Correlated Trait Complexes

An approach to increase the resolution power of interval mapping of quantitative trait (QT) loci is proposed, based on analysis of correlated trait complexes. For a given set of QTs, the broad sense heritability attributed to a QT locus (QTL) (say, A/ a) is an increasing function of the number of traits. Thus, for some traits x and y, H(xy)(2) (A/ a) >/= H(x)(2) (A/ a). The last inequality holds even if y does not depend on A/ a at all, but x and y are correlated within the groups AA, Aa and aa due to nongenetic factors and segregation of genes from other chromosomes. A simple relationship connects H(2) (both in single trait and two-trait analysis) with the expected LOD value, ELOD = -1/2N log(1 - H(2)). Thus, situations could exist that from the inequality H(xy)(2) (A/ a) >/= H(x)(2) (A/ a) a higher resolution is provided by the two-trait analysis as compared to the single-trait analysis, in spite of the increased number of parameters. Employing LOD-score procedure to simulated backcross data, we showed that the resolution power of the QTL mapping model can be elevated if correlation between QTs is taken into account. The method allows us to test numerous biologically important hypotheses concerning manifold effects of genomic segments on the defined trait complex (means, variances and correlations).