Molecular-Marker-Facilitated Investigations of Quantitative-Trait Loci in Maize. I. Numbers, Genomic Distribution and Types of Gene Action

Individual genetic factors which underlie variation in quantitative traits of maize were investigated in each of two F2 populations by examining the mean trait expressions of genotypic classes at each of 17-20 segregating marker loci. It was demonstrated that the trait expression of marker locus classes could be interpreted in terms of genetic behavior at linked quantitative trait loci (QTLs). For each of 82 traits evaluated, QTLs were detected and located to genomic sites. The numbers of detected factors varied according to trait, with the average trait significantly influenced by almost two-thirds of the marked genomic sites. Most of the detected associations between marker loci and quantitative traits were highly significant, and could have been detected with fewer than the 1800-1900 plants evaluated in each population. The cumulative, simple effects of marker-linked regions of the genome explained between 8 and 40% of the phenotypic variation for a subset of 25 traits evaluated. Single marker loci accounted for between 0.3% and 16% of the phenotypic variation of traits. Individual plant heterozygosity, as measured by marker loci, was significantly associated with variation in many traits. The apparent types of gene action at the QTLs varied both among traits and between loci for given traits, although overdominance appeared frequently, especially for yield-related traits. The prevalence of apparent overdominance may reflect the effects of multiple QTLs within individual marker-linked regions, a situation which would tend to result in overestimation of dominance. Digenic epistasis did not appear to be important in determining the expression of the quantitative traits evaluated. Examination of the effects of marked regions on the expression of pairs of traits suggests that genomic regions vary in the direction and magnitudes of their effects on trait correlations, perhaps providing a means of selecting to dissociate some correlated traits. Marker-facilitated investigations appear to provide a powerful means of examining aspects of the genetic control of quantitative traits. Modifications of the methods employed herein will allow examination of the stability of individual gene effects in varying genetic backgrounds and environments.     

Controlling the Type I and Type II Errors in Mapping Quantitative Trait Loci

Although the interval mapping method is widely used for mapping quantitative trait loci (QTLs), it is not very well suited for mapping multiple QTLs. Here, we present the results of a computer simulation to study the application of exact and approximate models for multiple QTLs. In particular, we focus on an automatic two-stage procedure in which in the first stage ``important' markers are selected in multiple regression on markers. In the second stage a QTL is moved along the chromosomes by using the pre-selected markers as cofactors, except for the markers flanking the interval under study. A refined procedure for cases with large numbers of marker cofactors is described. Our approach will be called MQM mapping, where MQM is an acronym for ``multiple-QTL models' as well as for ``marker-QTL-marker.' Our simulation work demonstrates the great advantage of MQM mapping compared to interval mapping in reducing the chance of a type I error (i.e., a QTL is indicated at a location where actually no QTL is present) and in reducing the chance of a type II error (i.e., a QTL is not detected).     

Quantitative trait loci for body growth and sex determination in the hermaphrodite teleost fish Sparus aurata L.

Gilthead sea bream (Sparus aurata L.) is an important marine fish in Mediterranean aquaculture. Sex determination by age and/or body weight is a critical life‐history trait, the genetic basis for which is largely unknown in this sequential hermaphrodite species. Herein, we performed a partial genome scan to map quantitative trait loci (QTL) affecting body weight and sex using 74 informative microsatellite markers from 10 paternal half‐sib families to construct nine linkage groups (LG). In total, four growth‐related QTL (two chromosome‐wide and two genome‐wide) and six QTL related to sex determination (three pairs in three different LGs) were detected (two chromosome‐wide and one genome‐wide). The proportion of phenotypic variation explained by the body‐weight QTL ranged from 9.3% to 17.2%, showing their potential for use in marker‐assisted selection. The results obtained offer solid ground to investigate the structure and function of the genomic regions involved in the mechanisms of sex reversal.     

基于F3种子的胚乳性状QTL区间定位

文章提出了包括胚乳效应和母体效应的胚乳性状QTL定位的统计方法,该方法的实验设计是分子标记基因型信息来自F2母体植株和F3种子胚(或植株),胚乳性状表型值来自F3单粒种子胚乳,称之为两步等级设计。同时,用计算机全面模拟以验证该模型的可行性,模拟结果表明,只要群体足够大,该模型能较有效地进行胚乳性状QTL定位并精确地估计出胚乳QTL的各种遗传效应和母体效应。     

Mapping quantitative trait loci using binned genotypes

Precise mapping of quantitative trait loci(QTLs)is critical for assessing genetic effects and identifying candidate genes for quantitative traits.Interval and composite interval mappings have been the methods of choice for several decades,which have provided tools for identifying genomic regions harboring causal genes for quantitative traits.Historically,the concept was developed on the basis of sparse marker maps where genotypes of loci within intervals could not be observed.Currently,genomes of many organisms have been saturated with markers due to the new sequencing technologies.Genotyping by sequencing usually generates hundreds of thousands of single nucleotide polymorphisms(SNPs),which often include the causal polymorphisms.The concept of interval no longer exists,prompting the necessity of a norm change in QTL mapping technology to make use of the high-volume genomic data.Here we developed a statistical method and a software package to map QTLs by binning markers into haplotype blocks,called bins.The new method detects associations of bins with quantitative traits.It borrows the mixed model methodology with a polygenic control from genome-wide association studies(GWAS)and can handle all kinds of experimental populations under the linear mixed model(LMM)framework.We tested the method using both simulated data and data from populations of rice.The results showed that this method has higher power than the current methods.An R package named binQTL is available from GitHub.     

Quantitative trait loci for spawning date and body weight in rainbow trout: testing for conserved effects across ancestrally duplicated chromosomes

We incorporated 69 microsatellite loci into an existing data set of 132 markers to test for quantitative trait loci (QTLs) affecting spawning date and body weight in a backcross between two outbred strains of rainbow trout (Oncorhynchus mykiss). Twenty-six linkage groups were identified and synteny of duplicated microsatellite markers was used to confirm 13 homeologous chromosome pairs. Gene-centromere data were used to localize the centromeres for 13 linkage groups whose orientations were previously unknown. We applied a combination of interval mapping and single marker analysis to the segregating maternal and paternal alleles at 201 microsatellite loci. Four spawning date QTLs with suggestive evidence for an additional two QTLs were detected in female trout spawning at 3 and 4 years of age. Similarly we detected three QTLs for body weight in females at 2 years of age plus four suggestive QTLs for this trait. We found marginal evidence that three pairs of ancestral homeologues contained detectable QTLs for the same trait. In one of the three pairs of homeologues, the duplicated QTL regions mapped to the same relative chromosomal location, while the exact localization of the QTL position in one of the other pairs was difficult to infer since it was based on data from a male-derived map. The existing data were unable to refute a hypothesis that duplicated functional genes will be maintained within the telomeric regions of salmonids due to preferential male-mediated crossing over in this region. Two of the four spawning date QTLs were detected on linkage groups with unknown homeologous relationships. QTLs with possible pleiotropic effects on both spawning date and body size were localized to two linkage groups.     

标记基因型中QTL基因型条件概率分布

随着分子数量遗传学的发展,人们提出了很多统计模型用于QTL定位分析。在这些模型中,首先得确定QTL在标记基因型中的条件概率分布,然后利用适当的统计方法对QTL在基因组中所处的位置进行估计。本文讨论了常见作图群体（如F2和回交群体）中在给定标记基因型下QTL的条件概率分布,提出了用Mathematics软件推导QTL基因型条件概率分布的方法。用该方法能够快速地、准确地推导出比较复杂的标记基因型中QTL基因型的条件概率分布。     

Relationship estimation by Markov-process models in a sib-pair linkage study

The results of sib-pair linkage studies may be compromised if a substantial number of putative sib pairs are not actually sib pairs. For classification of pairs in a sib-pair genome scan, I propose multipoint methods that are based on a Markov-process model of allele sharing along the chromosome. These methods can be implemented by standard algorithms that compute multipoint marker allele-sharing probabilities for sib pairs. When marker data from at least half the genome are used, misclassification rates are small. The methods will be implemented in an upcoming version of the computer software package S.A.G.E.     

Genetic variation of growth dynamics in maize (Zea mays L.) revealed through automated non‐invasive phenotyping

Hitherto, most quantitative trait loci of maize growth and biomass yield have been identified for a single time point, usually the final harvest stage. Through this approach cumulative effects are detected, without considering genetic factors causing phase‐specific differences in growth rates. To assess the genetics of growth dynamics, we employed automated non‐invasive phenotyping to monitor the plant sizes of 252 diverse maize inbred lines at 11 different developmental time points; 50 k SNP array genotype data were used for genome‐wide association mapping and genomic selection. The heritability of biomass was estimated to be over 71%, and the average prediction accuracy amounted to 0.39. Using the individual time point data, 12 main effect marker‐trait associations (MTAs) and six pairs of epistatic interactions were detected that displayed different patterns of expression at various developmental time points. A subset of them also showed significant effects on relative growth rates in different intervals. The detected MTAs jointly explained up to 12% of the total phenotypic variation, decreasing with developmental progression. Using non‐parametric functional mapping and multivariate mapping approaches, four additional marker loci affecting growth dynamics were detected. Our results demonstrate that plant biomass accumulation is a complex trait governed by many small effect loci, most of which act at certain restricted developmental phases. This highlights the need for investigation of stage‐specific growth affecting genes to elucidate important processes operating at different developmental phases.     

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.     

Genome-wide search for markers associated with osteochondrosis in Hanoverian warmblood horses

A genome-wide scan was performed to detect quantitative trait loci (QTLs) for osteochondrosis (OC) and osteochondrosis dissecans (OCD) in horses. The marker set comprised 260 microsatellites. We collected data from 211 Hanoverian warmblood horses consisting of 14 paternal half-sib families. Traits used were OC (fetlock and/or hock joints affected), OCD (fetlock and/or hock joints affected), fetlock OC, fetlock OCD, hock OC, and hock OCD. The first genome scan included 172 microsatellite markers. In a second step 88 additional markers were chosen to refine putative QTLs found in the first scan. Genome-wide significant QTLs were located on equine chromosomes 2, 4, 5, and 16. QTLs for fetlock OC and hock OC partly overlapped on the same chromosomes, indicating that these traits may be genetically related. QTLs reached the chromosome-wide significance level on eight different equine chromosomes: 2, 3, 4, 5, 15, 16, 19, and 21. This whole-genome scan was a first step toward the identification of candidate genome regions harboring genes responsible for equine OC. Further investigations are necessary to refine the map positions of the QTLs already identified for OC. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The power of interval mapping of quantitative trait loci, using selected sib pairs.

The interval-mapping procedure of Fulker and Cardon for analysis of a quantitative-trait loci (QTL) is extended for application to selected samples of sib pairs. Phenotypic selection of sib pairs, which is known to yield striking increases in power when a single marker is used, provides further increases in power when the interval-mapping approach is used. The greatest benefits of the combined approach are apparent with coarse maps, where QTLs of relatively modest (15%-20%) heritability can be detected with widely spaced markers (40-60 cM apart) in reasonably sized sibling samples. Useful information concerning QTL location is afforded by interval mapping in both selected and unselected samples.     

Genomewide scan of hoarding in sib pairs in which both sibs have Gilles de la Tourette syndrome

A genome scan of the hoarding phenotype (a component of obsessive-compulsive disorder) was conducted on 77 sib pairs collected by the Tourette Syndrome Association International Consortium for Genetics (TSAICG). All sib pairs were concordant for a diagnosis of Gilles de la Tourette syndrome (GTS). However, the analyses reported here were conducted for hoarding as both a dichotomous trait and a quantitative trait. Not all sib pairs in the sample were concordant for hoarding. Standard linkage analyses were performed using GENEHUNTER and Haseman-Elston methods. In addition, novel analyses with a recursive-partitioning technique were employed. Significant allele sharing was observed for both the dichotomous and the quantitative hoarding phenotypes for markers at 4q34-35 (P=.0007), by use of GENEHUNTER, and at 5q35.2-35.3 (P=.000002) and 17q25 (P=.00002), by use of the revisited Haseman-Elston method. The 4q site is in proximity to D4S1625, which was identified by the TSAICG as a region linked to the GTS phenotype. The recursive-partitioning technique examined multiple markers simultaneously. Results suggest joint effects of specific loci on 5q and 4q, with an overall P value of.000003. Although P values were not adjusted for multiple comparison, nearly all were much smaller than the customary significance level of.0001 for genomewide scans.     

Selection and subsequent analysis of sib pair data for QTL detection.

Haseman and Elston (1972) developed a robust regression method for the detection of linkage between a marker and a quantitative trait locus (QTL) using sib pair data. The principle underlying this method is that the difference in phenotypes between pairs of sibs becomes larger as they share a decreasing number of alleles at a particular QTL identical by descent (IBD) from their parents. In this case, phenotypically very different sibs will also on average share a proportion of alleles IBD at any marker linked to the QTL that is lower than the expected value of 0.5. Thus, the deviation of the proportion of marker alleles IBD from the expected value in pairs of sibs selected to be phenotypically different (i.e. discordant) can provide a test for the presence of a QTL. A simple regression method for QTL detection in sib pairs selected for high phenotypic differences is presented here. The power of the analytical method was found to be greater than the power obtained using the standard analysis when samples of sib pairs with high phenotypic differences were used. However, the use of discordant sib pairs was found to be less powerful for QTL detection than alternative selective genotyping schemes based on the phenotypic values of the sibs except with intense selection, when its advantage was only marginal. The most effective selection scheme overall was the use of sib pairs from entire families selected on the basis of high within-family variance for the trait in question. There is little effect of selection on QTL position estimates, which are in good agreement with the simulated values. However, QTL variance estimates are biased to a greater or lesser degree, depending on the selection method.     

Bayesian multilocus association mapping on ordinal and censored traits and its application to the analysis of genetic variation among Oryza sativa L. germplasms

Association mapping can be a powerful tool for detecting quantitative trait loci (QTLs) without requiring line-crossing experiments. We previously proposed a Bayesian approach for simultaneously mapping multiple QTLs by a regression method that directly incorporates estimates of the population structure. In the present study, we extended our method to analyze ordinal and censored traits, since both types of traits are common in the evaluation of germplasm collections. Ordinal-probit and tobit models were employed to analyze ordinal and censored traits, respectively. In both models, we postulated the existence of a latent continuous variable associated with the observable data, and we used a Markov-chain Monte Carlo algorithm to sample the latent variable and determine the model parameters. We evaluated the efficiency of our approach by using simulated- and real-trait analyses of a rice germplasm collection. Simulation analyses based on real marker data showed that our models could reduce both false-positive and false-negative rates in detecting QTLs to reasonable levels. Simulation analyses based on highly polymorphic marker data, which were generated by coalescent simulations, showed that our models could be applied to genotype data based on highly polymorphic marker systems, like simple sequence repeats. For the real traits, we analyzed heading date as a censored trait and amylose content and the shape of milled rice grains as ordinal traits. We found significant markers that may be linked to previously reported QTLs. Our approach will be useful for whole-genome association mapping of ordinal and censored traits in rice germplasm collections.     

Effectiveness of selective genotyping for detection of quantitative trait loci: an analysis of grain and malt quality traits in three barley populations.

Marker genotype data and grain and malt quality phenotype data from three barley (Hordeum vulgare L.) mapping populations were used to investigate the feasibility of selective genotyping for detection of quantitative trait loci (QTLs). With selective genotyping, only individuals with high and low phenotypic values for the trait of interest are genotyped. Here, genotyping of 10 to 70% of each population (i.e., 5 to 35% in each tail of the phenotypic distribution) was considered. Genomic positions detected by selective genotyping were compared to QTL position estimates from interval mapping analysis using marker genotype data from the entire population. Selective genotyping reliably detected almost all of the mapped QTLs, often with only 10% of the population genotyped. Selective genotyping also detected spurious QTLs in regions of the genome where no significant QTL had been mapped. Even with additional genotyping to verify putative QTLs, the total genotyping effort for detection of QTLs for a single trait by selective genotyping was usually less than 30% of that required for conventional interval mapping. Simultaneous investigation of two or more traits by selective genotyping would require additional genotyping effort, but could still be worthwhile.     

A genomewide linkage scan for quantitative-trait loci for obesity phenotypes

Obesity is an increasingly serious health problem in the world. Body mass index (BMI), percentage fat mass, and body fat mass are important indices of obesity. For a sample of pedigrees that contains >10,000 relative pairs (including 1,249 sib pairs) that are useful for linkage analyses, we performed a whole-genome linkage scan, using 380 microsatellite markers to identify genomic regions that may contain quantitative-trait loci (QTLs) for obesity. Each pedigree was ascertained through a proband who has extremely low bone mass, which translates into a low BMI. A major QTL for BMI was identified on 2q14 near the marker D2S347 with a LOD score of 4.04 in two-point analysis and a maximum LOD score (MLS) of 4.44 in multipoint analysis. The genomic region near 2q14 also achieved an MLS >2.0 for percentage of fat mass and body fat mass. For the putative QTL on 2q14, as much as 28.2% of BMI variation (after adjustment for age and sex) may be attributable to this locus. In addition, several other genomic regions that may contain obesity-related QTLs are suggested. For example, 1p36 near the marker D1S468 may contain a QTL for BMI variation, with a LOD score of 2.75 in two-point analysis and an MLS of 2.09 in multipoint analysis. The genomic regions identified in this and earlier reports are compared for further exploration in extension studies that use larger samples and/or denser markers for confirmation and fine-mapping studies, to eventually identify major functional genes involved in obesity.     

A Score-Statistic Approach for the Mapping of Quantitative-Trait Loci with Sibships of Arbitrary Size

The Haseman-Elston method is widely used for the mapping of quantitative-trait loci. However, this method does not use all the information in the data, because it only considers the sib-pair trait-value difference. In addition, the Haseman-Elston method was developed for independent sib pairs; its generalization to nonindependent sib pairs is not straightforward. Here we introduce a score test statistic derived from a normal likelihood based on multiplex sibship data, conditional on identical-by-descent sharing statuses. This score test is asymptotically equivalent to the corresponding likelihood-ratio test, but it is much easier to implement. Because the proposed test uses all of the trait values, it makes more efficient use of the data than does the Haseman-Elston method. The proposed test is naturally applicable to sibships of arbitrary size. The finite-sample properties of the proposed score statistic are evaluated via simulations.     

QTL analysis of morphological traits in a tomato recombinant inbred line population.

Quantitative trait loci (QTLs) influencing morphological traits were identified by restriction fragment length polymorphism (RFLP) analysis in a population of recombinant inbred lines (RILs) derived from a cross of the cultivated tomato Lycopersicon esculentum with a related wild species, Lycopersicon cheesmanii. One hundred and thirty-two RFLP loci spaced throughout the tomato genome were used as DNA probes on genomic DNA from 97 RIL families. Morphological traits, including plant height, plant fresh mass, number of branches, number of nodes, first flower-bearing node, and leaf length, were evaluated in two controlled environment trials in 1992 and 1993. QTLs were detected via regression analyses at multiple marker loci for each morphological trait. A total of 41 markers were significantly associated with the traits examined. Large additive effects were measured at many of these loci. QTLs for multiple traits were detected on chromosomes 3 (TG74) and 4 (CT188), suggesting the possible association of these chromosome segments with genes controlling growth and development in tomato. These chromosomal regions were also associated with multiple morphological traits in a L. esculentum x Lycopersicon pennellii cross. A total of 13% of the QTLs identified for traits common to both studies occupied similar map positions.     

Complete multipoint sib-pair analysis of qualitative and quantitative traits.

Sib-pair analysis is an increasingly important tool for genetic dissection of complex traits. Current methods for sib-pair analysis are primarily based on studying individual genetic markers one at a time and thus fail to use the full inheritance information provided by multipoint linkage analysis. In this paper, we describe how to extract the complete multipoint inheritance information for each sib pair. We then describe methods that use this information to map loci affecting traits, thereby providing a unified approach to both qualitative and quantitative traits. Specifically, complete multipoint approaches are presented for (1) exclusion mapping of qualitative traits; (2) maximum-likelihood mapping of qualitative traits; (3) information-content mapping, showing the extent to which all inheritance information has been extracted at each location in the genome; and (4) quantitative-trait mapping, by two parametric methods and one nonparametric method. In addition, we explore the effects of marker density, marker polymorphism, and availability of parents on the information content of a study. We have implemented the analysis methods in a new computer package, MAPMAKER/SIBS. With this computer package, complete multipoint analysis with dozens of markers in hundreds of sib pairs can be carried out in minutes.