Genetic mapping and analysis of quantitative trait loci affecting fiber and lignin content in maize

Plant cell walls of forage provide a major source of energy for ruminant animals. Digestion of cell walls is limited by the presence of lignin, therefore the improving the digestibility of forages by reducing lignin content is a major goal in forage crop breeding programs. A recombinant inbred line maize population was used to map quantitative trait loci (QTL) for neutral detergent fiber (NDF), acid detergent fiber (ADF), and acid detergent lignin (ADL) of leaf-sheath and stalk tissues. All traits were positively genetically correlated. The larger genetic correlations were between NDF and ADF in sheaths (r = 0.84), NDF and ADF (r = 0.96), ADF and ADL (r = 0.83), and NDF and ADL (r = 0.76) in stalks. Twelve QTL were detected for NDF and 11 QTL for ADF in leaf-sheaths. Eight QTL detected for both traits were defined by the same or linked marker loci. Eight QTL were associated with leaf-sheath ADL. Eleven QTL were detected for NDF and ADF, and 12 QTL for ADL in stalks. Nine of eleven QTL detected for both NDF and ADF in stalks coincided in their genomic position. A high proportion of QTL detected for these traits had the same parental effects and genomic locations, suggesting that it is only necessary to select on one fiber component (NDF or ADF) to improve digestibility. Favorable correlated responses of unselected fiber components are expected due to coincident genomic locations of QTL and the high genetic correlation between fiber components. Several QTL detected in this study coincided in their positions with putative cellulose synthase genes from maize.     

In silico mapping of quantitative trait loci in maize

Quantitative trait loci (QTL) are most often detected through designed mapping experiments. An alternative approach is in silico mapping, whereby genes are detected using existing phenotypic and genomic databases. We explored the usefulness of in silico mapping via a mixed-model approach in maize (Zea mays L.). Specifically, our objective was to determine if the procedure gave results that were repeatable across populations. Multilocation data were obtained from the 1995–2002 hybrid testing program of Limagrain Genetics in Europe. Nine heterotic patterns comprised 22,774 single crosses. These single crosses were made from 1,266 inbreds that had data for 96 simple sequence repeat (SSR) markers. By a mixed-model approach, we estimated the general combining ability effects associated with marker alleles in each heterotic pattern. The numbers of marker loci with significant effects—37 for plant height, 24 for smut [Ustilago maydis (DC.) Cda.] resistance, and 44 for grain moisture—were consistent with previous results from designed mapping experiments. Each trait had many loci with small effects and few loci with large effects. For smut resistance, a marker in bin 8.05 on chromosome 8 had a significant effect in seven (out of a maximum of 18) instances. For this major QTL, the maximum effect of an allele substitution ranged from 5.4% to 41.9%, with an average of 22.0%. We conclude that in silico mapping via a mixed-model approach can detect associations that are repeatable across different populations. We speculate that in silico mapping will be more useful for gene discovery than for selection in plant breeding programs.     

High-resolution association mapping of quantitative trait loci: a population-based approach

In this article, population-based regression models are proposed for high-resolution linkage disequilibrium mapping of quantitative trait loci (QTL). Two regression models, the "genotype effect model" and the "additive effect model," are proposed to model the association between the markers and the trait locus. The marker can be either diallelic or multiallelic. If only one marker is used, the method is similar to a classical setting by Nielsen and Weir, and the additive effect model is equivalent to the haplotype trend regression (HTR) method by Zaykin et al. If two/multiple marker data with phase ambiguity are used in the analysis, the proposed models can be used to analyze the data directly. By analytical formulas, we show that the genotype effect model can be used to model the additive and dominance effects simultaneously; the additive effect model takes care of the additive effect only. On the basis of the two models, F-test statistics are proposed to test association between the QTL and markers. By a simulation study, we show that the two models have reasonable type I error rates for a data set of moderate sample size. The noncentrality parameter approximations of F-test statistics are derived to make power calculation and comparison. By a simulation study, it is found that the noncentrality parameter approximations of F-test statistics work very well. Using the noncentrality parameter approximations, we compare the power of the two models with that of the HTR. In addition, a simulation study is performed to make a comparison on the basis of the haplotype frequencies of 10 SNPs of angiotensin-1 converting enzyme (ACE) genes.     

Genetic analysis of two new quantitative trait loci for ear weight in maize inbred line Huangzao4

Ear weight is one of the most important agronomic traits considered necessary in maize (Zea mays L.) breeding projects. To determine its genetic basis, a population consisting of 239 recombinant inbred lines, derived from the cross Mo17 x Huangzao4, was used to detect quantitative trait loci (QTLs) for ear weight under two nitrogen regimes. Under a high nitrogen fertilization regime, one QTL was identified in chromosome bin 2.08-2.09, which explained 7.46% of phenotypic variance and an increase in ear weight of about 5.79 g, owing to an additive effect. Under a low nitrogen regime, another QTL was identified in chromosome bin 1.10-1.11; it accounted for 7.11% of phenotypic variance and a decrease of 5.24 g in ear weight, due to an additive effect. Based on comparisons with previous studies, these two QTLs are new loci associated with ear weight in maize. These findings contribute to our knowledge about the genetic basis of ear weight in maize.     

Computing power of quantitative trait locus association mapping for haploid loci

Background  

Statistical power calculations are a critical part of any study design for gene mapping. Most calculations assume that the locus of interest is biallelic. However, there are common situations in human genetics such as X-linked loci in males where the locus is haploid. The purpose of this work is to mathematically derive the biometric model for haploid loci, and to compute power for QTL mapping when the loci are haploid.     

Molecular mapping of quantitative trait loci for drought tolerance in maize plants

Drought tolerance is one of the most important but complex traits of crops. We looked for quantitative trait loci (QTLs) that affect drought tolerance in maize. Two maize inbreds and their advanced lines were evaluated for drought-related traits. A genetic linkage map developed using RFLP markers was used to identify QTLs associated with drought-related traits. Twenty-two QTLs were detected, with a minimum of one and a maximum of nine for drought-related traits. A single-QTL was detected for sugar concentration accounting for about 52.2% of the phenotypic variation on chromosome 6. A single-QTL was also identified for each of the traits root density, root dry weight, total biomass, relative water content, and leaf abscisic acid content, on chromosomes 1 and 7, contributing to 24, 0.2, 0.4, 7, and 19% of the phenotypic variance, respectively. Three QTLs were identified for grain yield on chromosomes 1, 5, and 9, explaining 75% of the observed phenotypic variability, whereas four QTLs were detected for osmotic potential on chromosomes 1, 3, and 9, together accounting for 50% of the phenotypic variance. Nine QTLs were detected for leaf surface area on chromosomes 3 and 9, with various degrees of phenotypic variance, ranging from 25.8 to 42.2%. Four major clusters of QTLs were identified on chromosomes 1, 3, 7, and 9. A QTL for yield on chromosome 1 was found co-locating with the QTLs for root traits, total biomass, and osmotic potential in a region of about 15 cM. A cluster of QTLs for leaf surface area were coincident with a QTL for osmotic potential on chromosome 3. The QTLs for leaf area also clustered on chromosome 9, whereas QTLs for leaf abscisic acid content and relative water content coincided on chromosome 7, 10 cM apart. Co-location of QTLs for different traits indicates potential pleiotropism or tight linkage, which may be useful for indirect selection in maize improvement for drought tolerance.     

Genetic mapping of quantitative trait loci for aseasonal reproduction in sheep

The productivity and economic prosperity of sheep farming could benefit greatly from more effective methods of selection for year-round lambing. Identification of QTL for aseasonal reproduction in sheep could lead to more accurate selection and faster genetic improvement. One hundred and twenty microsatellite markers were genotyped on 159 backcross ewes from a Dorset × East Friesian crossbred pedigree. Interval mapping was undertaken to map the QTL underlying several traits describing aseasonal reproduction including the number of oestrous cycles, maximum level of progesterone prior to breeding, pregnancy status determined by progesterone level, pregnancy status determined by ultrasound, lambing status and number of lambs born. Seven chromosomes (1, 3, 12, 17, 19, 20 and 24) were identified to harbour putative QTL for one or more component traits used to describe aseasonal reproduction. Ovine chromosomes 12, 17, 19 and 24 harbour QTL significant at the 5% chromosome-wide level, chromosomes 3 and 20 harbour QTL that exceeded the threshold at the 1% chromosome-wide level, while the QTL identified on chromosome 1 exceeded the 1% experiment-wide significance level. These results are a first step towards understanding the genetic mechanism of this complex trait and show that variation in aseasonal reproduction is associated with multiple chromosomal regions.     

Efficient association mapping of quantitative trait loci with selective genotyping

Selective genotyping (i.e., genotyping only those individuals with extreme phenotypes) can greatly improve the power to detect and map quantitative trait loci in genetic association studies. Because selection depends on the phenotype, the resulting data cannot be properly analyzed by standard statistical methods. We provide appropriate likelihoods for assessing the effects of genotypes and haplotypes on quantitative traits under selective-genotyping designs. We demonstrate that the likelihood-based methods are highly effective in identifying causal variants and are substantially more powerful than existing methods.     

Genetic dissection of tocopherol content and composition in maize grain using quantitative trait loci analysis and the candidate gene approach

Tocopherols are essential micronutrients for humans and animals, with several beneficial effects in plants. Among cereals, only maize grains contain high concentrations of tocopherols. In this investigation we analyzed, during 2004 and 2005, by high-performance liquid chromatography (HPLC), a population of 233 recombinant inbred lines (RIL) which were derived from two diverse parents and had extremely variable tocopherol content and composition. A genetic map was constructed using 208 polymorphic molecular markers including gene-targeted markers based on six candidate genes of the tocopherol biosynthesis pathway (HPPD, VTE1, VTE3, VTE4, P3VTE5, and P4VTE5). Thirty-one quantitative trait loci (QTL) associated with quantitative variation of tocopherol content and composition were identified by composite interval mapping (CIM); these were located on sixteen genomic regions covering all the chromosomes except chromosome 4. Most (65%) QTL were co-located, suggesting that in some cases the same QTL predominantly affected the amounts of more than one tocopherol. Two candidate genes, HPPD and VTE4 showed co-localization with major QTL for tocopherol content and composition whereas only one interval (umc1075–umc1304) on chromosome eight exhibited a QTL for α, δ, γ, and total tocopherols with high LOD and PVE values. The candidate genes associated with tocopherol content and with composition, especially VTE4 and HPPD, could be precisely used for alteration of the tocopherol content and composition of maize grains by development of functional markers. Other identified major QTL especially those on chromosomes 8, 1, and 2 (near candidate gene VTE5) can also be used for improvement of maize grain quality by marker-assisted selection.     

Genetic complexity and quantitative trait loci mapping of yeast morphological traits

Functional genomics relies on two essential parameters: the sensitivity of phenotypic measures and the power to detect genomic perturbations that cause phenotypic variations. In model organisms, two types of perturbations are widely used. Artificial mutations can be introduced in virtually any gene and allow the systematic analysis of gene function via mutants fitness. Alternatively, natural genetic variations can be associated to particular phenotypes via genetic mapping. However, the access to genome manipulation and breeding provided by model organisms is sometimes counterbalanced by phenotyping limitations. Here we investigated the natural genetic diversity of Saccharomyces cerevisiae cellular morphology using a very sensitive high-throughput imaging platform. We quantified 501 morphological parameters in over 50,000 yeast cells from a cross between two wild-type divergent backgrounds. Extensive morphological differences were found between these backgrounds. The genetic architecture of the traits was complex, with evidence of both epistasis and transgressive segregation. We mapped quantitative trait loci (QTL) for 67 traits and discovered 364 correlations between traits segregation and inheritance of gene expression levels. We validated one QTL by the replacement of a single base in the genome. This study illustrates the natural diversity and complexity of cellular traits among natural yeast strains and provides an ideal framework for a genetical genomics dissection of multiple traits. Our results did not overlap with results previously obtained from systematic deletion strains, showing that both approaches are necessary for the functional exploration of genomes.     

Genetic dissection of grain weight in bread wheat through quantitative trait locus interval and association mapping

Genetic dissection of grain weight in bread wheat was undertaken through both genome-wide quantitative trait locus (QTL) interval mapping and association mapping. QTL interval mapping involved preparation of a framework linkage map consisting of 294 loci {194 simple sequence repeats (SSRs), 86 amplified fragment length polymorphisms (AFLPs) and 14 selective amplifications of microsatellite polymorphic loci (SAMPL)} using a bi-parental recombinant inbred line (RIL) mapping population derived from Rye Selection111 × Chinese Spring. Using the genotypic data and phenotypic data on grain weight (GW) of RILs collected over six environments, genome-wide single locus QTL analysis was conducted to identify main effect QTL. This led to identification of as many as ten QTL including four major QTL (three QTL were stable), each contributing >20% phenotypic variation (PV) for GW. The above study was supplemented with association mapping, which allowed identification of 11 new markers in the genomic regions that were not reported earlier to harbour any QTL for GW. It also allowed identification of closely linked markers for six known QTL, and validation of eight QTL reported earlier. The QTL identified through QTL interval mapping and association mapping may prove useful in marker-assisted selection (MAS) for the development of cultivars with high GW in bread wheat.     

Bayesian LASSO for quantitative trait loci mapping

The mapping of quantitative trait loci (QTL) is to identify molecular markers or genomic loci that influence the variation of complex traits. The problem is complicated by the facts that QTL data usually contain a large number of markers across the entire genome and most of them have little or no effect on the phenotype. In this article, we propose several Bayesian hierarchical models for mapping multiple QTL that simultaneously fit and estimate all possible genetic effects associated with all markers. The proposed models use prior distributions for the genetic effects that are scale mixtures of normal distributions with mean zero and variances distributed to give each effect a high probability of being near zero. We consider two types of priors for the variances, exponential and scaled inverse-chi(2) distributions, which result in a Bayesian version of the popular least absolute shrinkage and selection operator (LASSO) model and the well-known Student's t model, respectively. Unlike most applications where fixed values are preset for hyperparameters in the priors, we treat all hyperparameters as unknowns and estimate them along with other parameters. Markov chain Monte Carlo (MCMC) algorithms are developed to simulate the parameters from the posteriors. The methods are illustrated using well-known barley data.     

Genetic mapping and confirmation of quantitative trait loci for grain chalkiness in rice

Grain chalkiness is a highly undesirable trait affecting rice grain quality and milled rice yield. In order to clarify the genetic basis of chalkiness, a recombinant inbred line population (RIL) derived from a cross between Beilu130 (a japonica cultivar with high chalkiness) and Jin23B (an indica cultivar with low chalkiness) was developed for quantitative trait locus (QTL) mapping. A total of 10 QTLs for white belly rate (WBR) and white core rate (WCR) were detected on eight different chromosomes over 2 years. Two QTLs for WBR (qWBR2 and qWBR5) showed similar chromosomal locations with GW2 and qSW5/GW5, which have been cloned previously to control the grain width and should be the right candidate genes. Three novel minor QTLs controlling WCR, qWCR1, qWCR3, and qWCR4 were further validated in near isogenic F₂ populations (NIL-F₂) and explained 26.1, 18.3, and 21.1% of the phenotypic variation, respectively. These QTLs could be targets for map-based cloning of the candidate genes to elucidate the molecular mechanism of chalkiness and for marker-assisted selection (MAS) in rice grain quality improvement.     

Lineage-specific mapping of quantitative trait loci

We present an approach for quantitative trait locus (QTL) mapping, termed as ‘lineage-specific QTL mapping'', for inferring allelic changes of QTL evolution along with branches in a phylogeny. We describe and analyze the simplest case: by adding a third taxon into the normal procedure of QTL mapping between pairs of taxa, such inferences can be made along lineages to a presumed common ancestor. Although comparisons of QTL maps among species can identify homology of QTLs by apparent co-location, lineage-specific mapping of QTL can classify homology into (1) orthology (shared origin of QTL) versus (2) paralogy (independent origin of QTL within resolution of map distance). In this light, we present a graphical method that identifies six modes of QTL evolution in a three taxon comparison. We then apply our model to map lineage-specific QTLs for inbreeding among three taxa of yellow monkey-flower: Mimulus guttatus and two inbreeders M. platycalyx and M. micranthus, but critically assuming outcrossing was the ancestral state. The two most common modes of homology across traits were orthologous (shared ancestry of mutation for QTL alleles). The outbreeder M. guttatus had the fewest lineage-specific QTL, in accordance with the presumed ancestry of outbreeding. Extensions of lineage-specific QTL mapping to other types of data and crosses, and to inference of ancestral QTL state, are discussed.     

Combined linkage and association mapping of quantitative trait loci by multiple markers

Using multiple diallelic markers, variance component models are proposed for high-resolution combined linkage and association mapping of quantitative trait loci (QTL) based on nuclear families. The objective is to build a model that may fully use marker information for fine association mapping of QTL in the presence of prior linkage. The measures of linkage disequilibrium and the genetic effects are incorporated in the mean coefficients and are decomposed into orthogonal additive and dominance effects. The linkage information is modeled in variance-covariance matrices. Hence, the proposed methods model both association and linkage in a unified model. On the basis of marker information, a multipoint interval mapping method is provided to estimate the proportion of allele sharing identical by descent (IBD) and the probability of sharing two alleles IBD at a putative QTL for a sib-pair. To test the association between the trait locus and the markers, both likelihood-ratio tests and F-tests can be constructed on the basis of the proposed models. In addition, analytical formulas of noncentrality parameter approximations of the F-test statistics are provided. Type I error rates of the proposed test statistics are calculated to show their robustness. After comparing with the association between-family and association within-family (AbAw) approach by Abecasis and Fulker et al., it is found that the method proposed in this article is more powerful and advantageous based on simulation study and power calculation. By power and sample size comparison, it is shown that models that use more markers may have higher power than models that use fewer markers. The multiple-marker analysis can be more advantageous and has higher power in fine mapping QTL. As an application, the Genetic Analysis Workshop 12 German asthma data are analyzed using the proposed methods.     

Genetic mapping of quantitative trait loci for resistance to Haemonchus contortus in sheep

This paper presents results from a mapping experiment to detect quantitative trait loci (QTL) for resistance to Haemonchus contortus infestation in merino sheep. The primary trait analysed was faecal worm egg count in response to artificial challenge at 6 months of age. In the first stage of the experiment, whole genome linkage analysis was used for broad-scale mapping. The animal resource used was a designed flock comprising 571 individuals from four half-sib families. The average marker spacing was about 20 cM. For the primary trait, 11 QTL (as chromosomal/family combinations) were significant at the 5% chromosome-wide level, with allelic substitution effects of between 0.19 and 0.38 phenotypic standard deviation units. In general, these QTL did not have a significant effect on faecal worm egg count recorded at 13 months of age. In the second stage of the experiment, three promising regions (located on chromosomes 1, 3 and 4) were fine-mapped. This involved typing more closely spaced markers on individuals from the designed flock as well as an additional 495 individuals selected from a related population with a deeper pedigree. Analysis was performed using a linkage disequilibrium–linkage approach, under additive, dominant and multiple QTL models. Of these, the multiple QTL model resulted in the most refined QTL positions, with resolutions of <10 cM achieved for two regions. Because of the moderate size of effect of the QTL, and the apparent age and/or immune status specificity of the QTL, it is suggested that a panel of QTL will be required for significant genetic gains to be achieved within industry via marker-assisted selection.     

Bayesian quantitative trait loci mapping for multiple traits

Most quantitative trait loci (QTL) mapping experiments typically collect phenotypic data on multiple correlated complex traits. However, there is a lack of a comprehensive genomewide mapping strategy for correlated traits in the literature. We develop Bayesian multiple-QTL mapping methods for correlated continuous traits using two multivariate models: one that assumes the same genetic model for all traits, the traditional multivariate model, and the other known as the seemingly unrelated regression (SUR) model that allows different genetic models for different traits. We develop computationally efficient Markov chain Monte Carlo (MCMC) algorithms for performing joint analysis. We conduct extensive simulation studies to assess the performance of the proposed methods and to compare with the conventional single-trait model. Our methods have been implemented in the freely available package R/qtlbim (http://www.qtlbim.org), which greatly facilitates the general usage of the Bayesian methodology for unraveling the genetic architecture of complex traits.     

Sib mating designs for mapping quantitative trait loci

The power to separate the variance of a quantitative trait locus (QTL) from the polygenic variance is determined by the variability of genes identical by descent (IBD) at the QTL. This variability may increase with inbreeding. Selfing, the most extreme form of inbreeding, increases the variability of the IBD value shared by siblings, and thus has a higher efficiency for QTL mapping than random mating. In self-incompatible organisms, sib mating is the closest form of inbreeding. Similar to selfing, sib mating may also increase the power of QTL detection relative to random mating. In this study, we develop an IBD-based method under sib mating designs for QTL mapping. The efficiency of sib mating is then compared with random mating. Monte Carlo simulations show that sib mating designs notably increase the power for QTL detection. When power is intermediate, the power to detect a QTL using full-sib mating is, on average, 7% higher than under random mating. In addition, the IBD-based method proposed in this paper can be used to combine data from multiple families. As a result, the estimated QTL parameters can be applied to a wide statistical inference space relating to the entire reference population. This revised version was published online in July 2006 with corrections to the Cover Date.     

Multiple interval mapping for quantitative trait loci.

A new statistical method for mapping quantitative trait loci (QTL), called multiple interval mapping (MIM), is presented. It uses multiple marker intervals simultaneously to fit multiple putative QTL directly in the model for mapping QTL. The MIM model is based on Cockerham's model for interpreting genetic parameters and the method of maximum likelihood for estimating genetic parameters. With the MIM approach, the precision and power of QTL mapping could be improved. Also, epistasis between QTL, genotypic values of individuals, and heritabilities of quantitative traits can be readily estimated and analyzed. Using the MIM model, a stepwise selection procedure with likelihood ratio test statistic as a criterion is proposed to identify QTL. This MIM method was applied to a mapping data set of radiata pine on three traits: brown cone number, tree diameter, and branch quality scores. Based on the MIM result, seven, six, and five QTL were detected for the three traits, respectively. The detected QTL individually contributed from approximately 1 to 27% of the total genetic variation. Significant epistasis between four pairs of QTL in two traits was detected, and the four pairs of QTL contributed approximately 10.38 and 14.14% of the total genetic variation. The asymptotic variances of QTL positions and effects were also provided to construct the confidence intervals. The estimated heritabilities were 0.5606, 0.5226, and 0. 3630 for the three traits, respectively. With the estimated QTL effects and positions, the best strategy of marker-assisted selection for trait improvement for a specific purpose and requirement can be explored. The MIM FORTRAN program is available on the worldwide web (http://www.stat.sinica.edu.tw/chkao/).