Combined genetic and physiological analysis of a locus contributing to quantitative variation

The natural variation of many traits is controlled by multiple genes, individually referred to as quantitative trait loci (QTL), that interact with the environment to determine the ultimate phenotype of any individual. A QTL has yet to be described molecularly, in part because strategies to systematically identify them are underdeveloped and because the subtle nature of QTLs prevents the application of standard methods of gene identification. Therefore, it will be necessary to develop a systematic approach(es) for the identification of QTLs based upon the numerous positional data now being accumulated through molecular marker analyses. We have characterized a QTL by the following three-step approach: (1) identification of a QTL in complex populations, (2) isolation and genetic mapping of this QTL in near-isogenic lines, and (3) identification of a candidate gene using map position and physiological criteria. Using this approach we have characterized a plant height QTL in maize that maps to chromosome 9 near the centromere. Both map position and physiological criteria suggest the gibberillin biosynthesis gene dwarf3 as a candidate gene for this QTL.     

The analysis of QTL by simultaneous use of the full linkage map

An extension of interval mapping is presented that incorporates all intervals on the linkage map simultaneously. The approach uses a working model in which the sizes of putative QTL for all intervals across the genome are random effects. An outlier detection method is used to screen for possible QTL. Selected QTL are subsequently fitted as fixed effects. This screening and selection approach is repeated until the variance component for QTL sizes is not statistically significant. A comprehensive simulation study is conducted in which map uncertainty is included. The proposed method is shown to be superior to composite interval mapping in terms of power of detection of QTL. There is an increase in the rate of false positive QTL detected when using the new approach, but this rate decreases as the population size increases. The new approach is much simpler computationally. The analysis of flour milling yield in a doubled haploid population illustrates the improved power of detection of QTL using the approach, and also shows how vital it is to allow for sources of non-genetic variation in the analysis.     

Fine-Scale Mapping of Quantitative Trait Loci Using Historical Recombinations

With increasing popularity of QTL mapping in economically important animals and experimental species, the need for statistical methodology for fine-scale QTL mapping becomes increasingly urgent. The ability to disentangle several linked QTL depends on the number of recombination events. An obvious approach to increase the recombination events is to increase sample size, but this approach is often constrained by resources. Moreover, increasing the sample size beyond a certain point will not further reduce the length of confidence interval for QTL map locations. The alternative approach is to use historical recombinations. We use analytical methods to examine the properties of fine QTL mapping using historical recombinations that are accumulated through repeated intercrossing from an F(2) population. We demonstrate that, using the historical recombinations, both simple and multiple regression models can reduce significantly the lengths of support intervals for estimated QTL map locations and the variances of estimated QTL map locations. We also demonstrate that, while the simple regression model using historical recombinations does not reduce the variances of the estimated additive and dominant effects, the multiple regression model does. We further determine the power and threshold values for both the simple and multiple regression models. In addition, we calculate the Kullback-Leibler distance and Fisher information for the simple regression model, in the hope to further understand the advantages and disadvantages of using historical recombinations relative to F(2) data.     

Bayesian shrinkage analysis of quantitative trait Loci for dynamic traits

Many quantitative traits are measured repeatedly during the life of an organism. Such traits are called dynamic traits. The pattern of the changes of a dynamic trait is called the growth trajectory. Studying the growth trajectory may enhance our understanding of the genetic architecture of the growth trajectory. Recently, we developed an interval-mapping procedure to map QTL for dynamic traits under the maximum-likelihood framework. We fit the growth trajectory by Legendre polynomials. The method intended to map one QTL at a time and the entire QTL analysis involved scanning the entire genome by fitting multiple single-QTL models. In this study, we propose a Bayesian shrinkage analysis for estimating and mapping multiple QTL in a single model. The method is a combination between the shrinkage mapping for individual quantitative traits and the Legendre polynomial analysis for dynamic traits. The multiple-QTL model is implemented in two ways: (1) a fixed-interval approach where a QTL is placed in each marker interval and (2) a moving-interval approach where the position of a QTL can be searched in a range that covers many marker intervals. Simulation study shows that the Bayesian shrinkage method generates much better signals for QTL than the interval-mapping approach. We propose several alternative methods to present the results of the Bayesian shrinkage analysis. In particular, we found that the Wald test-statistic profile can serve as a mechanism to test the significance of a putative QTL.     

基于元分析的抗玉米丝黑穗病QTL比较定位

以玉米遗传连锁图谱IBM2 2005 Neighbors为参考图谱,通过映射整合不同试验中的抗玉米丝黑穗病QTL,构建QTL综合图谱。在国内外种质中,共发现22个抗病QTL,分布在除第7染色体外的9条玉米染色体上。采用元分析技术,获得2个“一致性”抗病QTL,图距分别为8.79 cM和18.92cM。从MaizeGDB网站下载“一致性”QTL区间内基因和标记的原始序列;采用NCBI网站在线软件BLASTx通过同源比对在2个“一致性”QTL区间内初步获得4个抗病位置候选基因。借助比较基因电子定位策略,将69个水稻和玉米抗性基因定位于玉米IBM2图谱上,在2个“一致性”QTL区间内分别发现1个水稻抗性基因,初步推断为抗病位置候选基因。本文结果为抗玉米丝黑穗病QTL精细定位和分子育种提供了基础。     

Recalculation of 23 mouse HDL QTL datasets improves accuracy and allows for better candidate gene analysis

In the past 15 years, the quantitative trait locus (QTL) mapping approach has been applied to crosses between different inbred mouse strains to identify genetic loci associated with plasma HDL cholesterol levels. Although successful, a disadvantage of this method is low mapping resolution, as often several hundred candidate genes fall within the confidence interval for each locus. Methods have been developed to narrow these loci by combining the data from the different crosses, but they rely on the accurate mapping of the QTL and the treatment of the data in a consistent manner. We collected 23 raw datasets used for the mapping of previously published HDL QTL and reanalyzed the data from each cross using a consistent method and the latest mouse genetic map. By utilizing this approach, we identified novel QTL and QTL that were mapped to the wrong part of chromosomes. Our new HDL QTL map allows for reliable combining of QTL data and candidate gene analysis, which we demonstrate by identifying Grin3a and Etv6, as candidate genes for QTL on chromosomes 4 and 6, respectively. In addition, we were able to narrow a QTL on Chr 19 to five candidates.     

Multivariate whole genome average interval mapping: QTL analysis for multiple traits and/or environments

A major aim in some plant-based studies is the determination of quantitative trait loci (QTL) for multiple traits or across multiple environments. Understanding these QTL by trait or QTL by environment interactions can be of great value to the plant breeder. A whole genome approach for the analysis of QTL is presented for such multivariate applications. The approach is an extension of whole genome average interval mapping in which all intervals on a linkage map are included in the analysis simultaneously. A random effects working model is proposed for the multivariate (trait or environment) QTL effects for each interval, with a variance-covariance matrix linking the variates in a particular interval. The significance of the variance-covariance matrix for the QTL effects is tested and if significant, an outlier detection technique is used to select a putative QTL. This QTL by variate interaction is transferred to the fixed effects. The process is repeated until the variance-covariance matrix for QTL random effects is not significant; at this point all putative QTL have been selected. Unlinked markers can also be included in the analysis. A simulation study was conducted to examine the performance of the approach and demonstrated the multivariate approach results in increased power for detecting QTL in comparison to univariate methods. The approach is illustrated for data arising from experiments involving two doubled haploid populations. The first involves analysis of two wheat traits, α-amylase activity and height, while the second is concerned with a multi-environment trial for extensibility of flour dough. The method provides an approach for multi-trait and multi-environment QTL analysis in the presence of non-genetic sources of variation.     

Quantifying evidence for candidate gene polymorphisms: Bayesian analysis combining sequence-specific and quantitative trait loci colocation information

We calculate posterior probabilities for candidate genes as a function of genomic location. Posterior probabilities for quantitative trait loci (QTL) presence in a small interval are calculated using a Bayesian model-selection approach based on the Bayesian information criterion (BIC) and used to combine QTL colocation information with sequence-specific evidence, e.g., from differential expression and/or association studies. Our method takes into account uncertainty in estimation of number and locations of QTL and estimated map position. Posterior probabilities for QTL presence were calculated for simulated data with n = 100, 300, and 1200 QTL progeny and compared with interval mapping and composite-interval mapping. Candidate genes that mapped to QTL regions had substantially larger posterior probabilities. Among candidates with a given Bayes factor, those that map near a QTL are more promising for further investigation with association studies and functional testing or for use in marker-aided selection. The BIC is shown to correspond very closely to Bayes factors for linear models with a nearly noninformative Zellner prior for the simulated QTL data with n > or = 100. It is shown how to modify the BIC to use a subjective prior for the QTL effects.     

QTL mapping in outbred half-sib families using Bayesian model selection

In this article, we propose a model selection method, the Bayesian composite model space approach, to map quantitative trait loci (QTL) in a half-sib population for continuous and binary traits. In our method, the identity-by-descent-based variance component model is used. To demonstrate the performance of this model, the method was applied to map QTL underlying production traits on BTA6 in a Chinese half-sib dairy cattle population. A total of four QTLs were detected, whereas only one QTL was identified using the traditional least square (LS) method. We also conducted two simulation experiments to validate the efficiency of our method. The results suggest that the proposed method based on a multiple-QTL model is efficient in mapping multiple QTL for an outbred half-sib population and is more powerful than the LS method based on a single-QTL model.     

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.     

Detection and modelling of time-dependent QTL in animal populations

A longitudinal approach is proposed to map QTL affecting function-valued traits and to estimate their effect over time. The method is based on fitting mixed random regression models. The QTL allelic effects are modelled with random coefficient parametric curves and using a gametic relationship matrix. A simulation study was conducted in order to assess the ability of the approach to fit different patterns of QTL over time. It was found that this longitudinal approach was able to adequately fit the simulated variance functions and considerably improved the power of detection of time-varying QTL effects compared to the traditional univariate model. This was confirmed by an analysis of protein yield data in dairy cattle, where the model was able to detect QTL with high effect either at the beginning or the end of the lactation, that were not detected with a simple 305 day model.     

High-resolution quantitative trait locus mapping for body weight in mice by recombinant progeny testing

A major obstacle to the positional cloning of quantitative trait loci (QTLs) lies in resolving genetic factors whose allelic effects are blurred by environmental and background genetic variation. We investigate a fine-mapping approach that combines the use of an interval-specific congenic strain with progeny testing of recombinants for markers flanking a QTL. We apply the approach to map a murine QTL with an approximately 20% effect on growth rate by progeny testing 39 recombinants in a 12 cM region of the X chromosome. We use a likelihood analysis in an attempt to maximize the information on QTL map location and effect. The major X-linked effect is mapped to an approximately 2 cM region flanked by markers about 5 cM apart, outside which LOD support for the QTL drops extremely steeply by about 80. Nearly unambiguous assignment of the QTL genotypic state is obtained for each recombinant. The resolution of individual recombinants in the region is therefore sufficiently high to facilitate the positional cloning of the locus, although progress has been hampered because the genomic region containing the QTL shows an exceptionally low level of polymorphism in comparison with recent studies.     

Linkage analysis of molecular markers and quantitative trait loci in populations of inbred backcross lines of Brassica napus L.

Backcross populations are often used to study quantitative trait loci (QTL) after they are initially discovered in balanced populations, such as F(2), BC(1), or recombinant inbreds. While the latter are more powerful for mapping marker loci, the former have the reduced background genetic variation necessary for more precise estimation of QTL effects. Many populations of inbred backcross lines (IBLs) have been developed in plant and animal systems to permit simultaneous study and dissection of quantitative genetic variation introgressed from one source to another. Such populations have a genetic structure that can be used for linkage estimation and discovery of QTL. In this study, four populations of IBLs of oilseed Brassica napus were developed and analyzed to map genomic regions from the donor parent (a winter-type cultivar) that affect agronomic traits in spring-type inbreds and hybrids. Restriction fragment length polymorphisms (RFLPs) identified among the IBLs were used to calculate two-point recombination fractions and LOD scores through grid searches. This information allowed the enrichment of a composite genetic map of B. napus with 72 new RFLP loci. The selfed and hybrid progenies of the IBLs were evaluated during two growing seasons for several agronomic traits. Both pedigree structure and map information were incorporated into the QTL analysis by using a regression approach. The number of QTL detected for each trait and the number of effective factors calculated by using biometrical methods were of similar magnitude. Populations of IBLs were shown to be valuable for both marker mapping and QTL analysis.     

Quantitative trait locus (QTL) isogenic recombinant analysis: a method for high-resolution mapping of QTL within a single population

In the quest for fine mapping quantitative trait loci (QTL) at a subcentimorgan scale, several methods that involve the construction of inbred lines and the generation of large progenies of such inbred lines have been developed (Complex Trait Consortium 2003). Here we present an alternative method that significantly speeds up QTL fine mapping by using one segregating population. As a first step, a rough mapping analysis is performed on a small part of the population. Once the QTL have been mapped to a chromosomal interval by standard procedures, a large population of 1000 plants or more is analyzed with markers flanking the defined QTL to select QTL isogenic recombinants (QIRs). QIRs bear a recombination event in the QTL interval of interest, while other QTL have the same homozygous genotype. Only these QIRs are subsequently phenotyped to fine map the QTL. By focusing at an early stage on the informative individuals in the population only, the efforts in population genotyping and phenotyping are significantly reduced as compared to prior methods. The principles of this approach are demonstrated by fine mapping an erucic acid QTL of rapeseed at a subcentimorgan scale.     

Family-based mapping of quantitative trait loci in plant breeding populations with resistance to Fusarium head blight in wheat as an illustration

Traditional quantitative trait loci (QTL) mapping approaches are typically based on early or advanced generation analysis of bi-parental populations. A limitation associated with this methodology is the fact that mapping populations rarely give rise to new cultivars. Additionally, markers linked to the QTL of interest are often not immediately available for use in breeding and they may not be useful within diverse genetic backgrounds. Use of breeding populations for simultaneous QTL mapping, marker validation, marker assisted selection (MAS), and cultivar release has recently caught the attention of plant breeders to circumvent the weaknesses of conventional QTL mapping. The first objective of this study was to test the feasibility of using family-pedigree based QTL mapping techniques generally used with humans and animals within plant breeding populations (PBPs). The second objective was to evaluate two methods (linkage and association) to detect marker-QTL associations. The techniques described in this study were applied to map the well characterized QTL, Fhb1 for Fusarium head blight resistance in wheat (Triticum aestivum L.). The experimental populations consisted of 82 families and 793 individuals. The QTL was mapped using both linkage (variance component and pedigree-wide regression) and association (using quantitative transmission disequilibrium test, QTDT) approaches developed for extended family-pedigrees. Each approach successfully identified the known QTL location with a high probability value. Markers linked to the QTL explained 40–50% of the phenotypic variation. These results show the usefulness of a human genetics approach to detect QTL in PBPs and subsequent use in MAS.     

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.     

QTL Express: mapping quantitative trait loci in simple and complex pedigrees

QTL Express is the first application for Quantitative Trait Locus (QTL) mapping in outbred populations with a web-based user interface. User input of three files containing a marker map, trait data and marker genotypes allows mapping of single or multiple QTL by the regression approach, with the option to perform permutation or bootstrap tests.     

Mouse inbred strain sequence information and yin-yang crosses for quantitative trait locus fine mapping

The shared ancestry of mouse inbred strains, together with the availability of sequence and phenotype information, is a resource that can be used to map quantitative trait loci (QTL). The difficulty in using only sequence information lies in the fact that in most instances the allelic state of the QTL cannot be unambiguously determined in a given strain. To overcome this difficulty, the performance of multiple crosses between various inbred strains has been proposed. Here we suggest and evaluate a general approach, which consists of crossing the two strains used initially to map the QTL and any new strain. We have termed these crosses "yin-yang," because they are complementary in nature as shown by the fact that the QTL will necessarily segregate in only one of the crosses. We used the publicly available SNP database of chromosome 16 to evaluate the mapping resolution achievable through this approach. Although on average the improvement of mapping resolution using only four inbred strains was relatively small (i.e., reduction of the QTL-containing interval by half at most), we found a great degree of variability among different regions of chromosome 16 with regard to mapping resolution. This suggests that with a large number of strains in hand, selecting a small number of strains may provide a significant contribution to the fine mapping of QTL.     

Pathotype-specific QTL for stem rust resistance in Lolium perenne

A genetic map populated with RAD and SSR markers was created from F1 progeny of a stem rust-susceptible and stem rust-resistant parent of perennial ryegrass (Lolium perenne). The map supplements a previous map of this population by having markers in common with several other Lolium spp. maps including EST-SSR anchor markers from a consensus map published by other researchers. A QTL analysis was conducted with disease severity and infection type data obtained by controlled inoculation of the population with each of two previously characterized pathotypes of Puccinia graminis subsp. graminicola that differ in virulence to different host plant genotypes in the F1 population. Each pathotype activated a specific QTL on one linkage group (LG): qLpPg1 on LG7 for pathotype 101, or qLpPg2 on LG1 for pathotype 106. Both pathotypes also activated a third QTL in common, qLpPg3 on LG6. Anchor markers, present on a consensus map, were located in proximity to each of the three QTL. These QTL had been detected also in previous experiments in which a genetically heterogeneous inoculum of the stem rust pathogen activated all three QTL together. The results of this and a previous study are consistent with the involvement of the pathotype-specific QTL in pathogen recognition and the pathotype-nonspecific QTL in a generalized resistance response. By aligning the markers common to other published reports, it appears that two and possibly all three of the stem rust QTL reported here are in the same general genomic regions containing some of the L. perenne QTL reported to be activated in response to the crown rust pathogen (P. coronata).     

The use of genotyping by sequencing in blackcurrant (Ribes nigrum): developing high-resolution linkage maps in species without reference genome sequences

A genotyping by sequencing (GbS) approach is reported in blackcurrant (Ribes nigrum L.) using a de novo read assembly method developed because of the current absence of a reference genome sequence for this species. A new approach to single nucleotide polymorphism (SNP) genotype calling is described, where individual genotypes for a large number of SNPs were characterised from the GbS counts using a novel method based on a functional regression of major and minor allele read counts. The high-quality GbS SNPs were combined with SNPs and simple sequence repeats generated from other technologies to develop a linkage map with increased marker density and improved genome coverage, containing up to 204 SNPs on each linkage group. SNPs of lower quality were then located on the map using quantitative trait locus (QTL) interval mapping of the proportion of the major allele. Two QTL each for 100-berry weight and Brix scores, measured over three years, were identified using the map. The use of this approach to identify and map a significant number of novel SNPs in a woody species with hitherto limited genomic resources may have generic application to other under-resourced and minor species in the development of cost-effective and efficient high-density genetic maps.