R/qtl: QTL mapping in experimental crosses

SUMMARY: R/qtl is an extensible, interactive environment for mapping quantitative trait loci (QTLs) in experimental populations derived from inbred lines. It is implemented as an add-on package for the freely-available statistical software, R, and includes functions for estimating genetic maps, identifying genotyping errors, and performing single-QTL and two-dimensional, two-QTL genome scans by multiple methods, with the possible inclusion of covariates. AVAILABILITY: The package is freely available at http://www.biostat.jhsph.edu/~kbroman/qtl.     

Mapping multiple QTL of different effects: comparison of a simple sequential testing strategy and multiple QTL mapping

The aim of this study was to explore, by computer simulation, the mapping of QTLs in a realistic but complex situation of many (linked) QTLs with different effects, and to compare two QTL mapping methods. A novel method to dissect genetic variation on multiple chromosomes using molecular markers in backcross and F₂ populations derived from inbred lines was suggested, and its properties tested using simulations. The rationale for this sequential testing method was to explicitly test for alternative genetic models. The method consists of a series of four basic statistical tests to decide whether variance was due to a single QTL, two QTLs, multiple QTLs, or polygenes, starting with a test to detect genetic variance associated with a particular chromosome. The method was able to distinguish between different QTL configurations, in that the probability to `detect' the correct model was high, varying from 0.75 to 1. For example, for a backcross population of 200 and an overall heritability of 50%, in 78% of replicates a polygenic model was detected when that was the underlying true model. To test the method for multiple chromosomes, QTLs were simulated on 10 chromosomes, following a geometric series of allele effects, assuming positive alleles were in coupling in the founder lines For these simulations, the sequential testing method was compared to the established Multiple QTL Mapping (MQM) method. For a backcross population of 400 individuals, power to detect genetic variance was low with both methods when the heritability was 0.40. For example, the power to detect genetic variation on a chromosome on which 6 QTLs explained 12.6% of the genetic variance, was less than 60% for both methods. For a large heritability (0.90), the power of MQM to detect genetic variance and to dissect QTL configurations was generally better, due to the simultaneous fitting of markers on all chromosomes. It is concluded that when testing different QTL configurations on a single chromosome using the sequential testing procedure, regions of other chromosomes which explain a significant amount of variation should be fitted in the model of analysis. This study reinforces the need for large experiments in plants and other species if the aim of a genome scan is to dissect quantitative genetic variation.     

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.     

A comparison of regression interval mapping and multiple interval mapping for linked QTL

Regression interval mapping and multiple interval mapping are compared with regard to mapping linked quantitative trait loci (QTL) in inbred-line cross experiments. For that purpose, a simulation study was performed using genetic models with two linked QTL. Data were simulated for F(2) populations of different sizes and with all QTL and marker alleles fixed for alternative alleles in the parental lines. The criteria for comparison are power of QTL identification and the accuracy of the QTL position and effect estimates. Further, the estimates of the relative QTL variance are assessed. There are distinct differences in the QTL position estimates between the two methods. Multiple interval mapping tends to be more powerful as compared to regression interval mapping. Multiple interval mapping further leads to more accurate QTL position and QTL effect estimates. The superiority increased with wider marker intervals and larger population sizes. If QTL are in repulsion, the differences between the two methods are very pronounced. For both methods, the reduction of the marker interval size from 10 to 5 cM increases power and greatly improves QTL parameter estimates. This contrasts with findings in the literature for single QTL scenarios, where a marker density of 10 cM is generally considered as sufficient. The use of standard (asymptotic) statistical theory for the computation of the standard errors of the QTL position and effect estimates proves to give much too optimistic standard errors for regression interval mapping as well as for multiple interval mapping.     

Limitations of a two-step moment method for mapping linked multiple QTL

As an alternative to multiple-interval mapping a two-step moment method was recently proposed to map linked multiple quantitative trait loci (QTLs). The advantage of this moment method was supposed to be its simplicity and computational efficiency, especially in detecting closely linked QTLs within a marker bracket, but also in mapping QTLs in different marker intervals. Using simulations it is shown that the two-step moment method may give poor results compared with multiple-interval mapping, irrespective of whether the QTLs are in the same or in different marker intervals, especially if linked QTLs are in repulsion. The criteria of comparison are number of identified QTLs, likelihood ratio test statistics, means and empirical standard errors of the QTL position and QTL effects estimates, and the accuracy of the residual variance estimates. Further, the joint conditional probabilities of QTL genotypes for two putative QTLs within a marker interval were derived and compared with the unmodified approach ignoring the non-independence of the conditional probabilities.     

QTL mapping of multiple foliar disease and root-lesion nematode resistances in wheat

A genetic linkage map, based on a cross between the synthetic hexaploid CPI133872 and the bread wheat cultivar Janz, was established using 111 F₁-derived doubled haploid lines. The population was phenotyped in multiple years and/or locations for seven disease resistance traits, namely, Septoria tritici blotch (Mycosphaeralla graminicola), yellow leaf spot also known as tan spot (Pyrenophora tritici-repentis), stripe rust (Puccinia striiformis f. sp. tritici), leaf rust (Puccinia triticina), stem rust (Puccinia graminis f. sp. tritici) and two species of root-lesion nematode (Pratylenchyus thornei and P. neglectus). The DH population was also scored for coleoptile colour and the presence of the seedling leaf rust resistance gene Lr24. Implementation of a multiple-QTL model identified a tightly linked cluster of foliar disease resistance QTL in chromosome 3DL. Major QTL each for resistance to Septoria tritici blotch and yellow leaf spot were contributed by the synthetic hexaploid parent CPI133872 and linked in repulsion with the coincident Lr24/Sr24 locus carried by parent Janz. This is the first report of linked QTL for Septoria tritici blotch and yellow leaf spot contributed by the same parent. Additional QTL for yellow leaf spot were detected in 5AS and 5BL. Consistent QTL for stripe rust resistance were identified in chromosomes 1BL, 4BL and 7DS, with the QTL in 7DS corresponding to the Yr18/Lr34 region. Three major QTL for P. thornei resistance (2BS, 6DS, 6DL) and two for P. neglectus resistance (2BS, 6DS) were detected. The recombinants combining resistance to Septoria tritici blotch, yellow leaf spot, rust diseases and root-lesion nematodes from parents CPI133872 and Janz constitute valuable germplasm for the transfer of multiple disease resistance into new wheat cultivars.     

QTL mapping in rice

In the past 10 years, interest in applying the tools of molecular genetics to the problem of increasing world rice production has resulted in the generation of two highly saturated, molecular linkage maps of rice, and the localization of numerous genes and quantitative trait loci (QTLs). Primary studies have identified QTLs associated with disease resistance, abiotic stress tolerance and yield potential of rice in a range of ecosystems. The ability to identify, manipulate and potentially clone individual genes involved in quantitatively inherited characters, combined with the demonstrated conservation of numerous linkage blocks among members of the grass family, emphasizes the contribution of map-based genetic analyses both to applied and to basic crop research.     

Human QTL linkage mapping

Human quantitative trait locus (QTL) linkage mapping, although based on classical statistical genetic methods that have been around for many years, has been employed for genome-wide screening for only the last 10–15 years. In this time, there have been many success stories, ranging from QTLs that have been replicated in independent studies to those for which one or more genes underlying the linkage peak have been identified to a few with specific functional variants that have been confirmed in in vitro laboratory assays. Despite these successes, there is a general perception that linkage approaches do not work for complex traits, possibly because many human QTL linkage studies have been limited in sample size and have not employed the family configurations that maximize the power to detect linkage. We predict that human QTL linkage studies will continue to be productive for the next several years, particularly in combination with RNA expression level traits that are showing evidence of regulatory QTLs of large effect sizes and in combination with high-density genome-wide SNP panels. These SNP panels are being used to identify QTLs previously localized by linkage and linkage results are being used to place informative priors on genome-wide association studies.     

MCQTL: multi-allelic QTL mapping in multi-cross design

The aim of the MCQTL software package is to perform QTL mapping in multi-cross designs. It allows the analysis of the usual populations derived from inbred lines and can link the families by assuming that the QTL locations are the same in all of them. Moreover, a diallel modelling of the QTL genotypic effects is allowed in multiple related families. The implemented model is a linear regression model. A composite interval mapping and an iterative QTL mapping are implemented to deal with multiple QTL models. Marker cofactor selections by forward or backward stepwise methods are implemented as well as computation of threshold test value by permutation. AVAILABILITY: The program is available on request after signing a licence agreement; free of charge for academic and non-profit organizations at http://www.genoplante.org (Bioinformatics products).     

R/qtlbim: QTL with Bayesian Interval Mapping in experimental crosses

R/qtlbim is an extensible, interactive environment for the Bayesian Interval Mapping of QTL, built on top of R/qtl (Broman et al., 2003), providing Bayesian analysis of multiple interacting quantitative trait loci (QTL) models for continuous, binary and ordinal traits in experimental crosses. It includes several efficient Markov chain Monte Carlo (MCMC) algorithms for evaluating the posterior of genetic architectures, i.e. the number and locations of QTL, their main and epistatic effects and gene-environment interactions. R/qtlbim provides extensive informative graphical and numerical summaries, and model selection and convergence diagnostics of the MCMC output, illustrated through the vignette, example and demo capabilities of R (R Development Core Team 2006). Availability: The package is freely available from cran.r-project.org.     

QTL mapping for fiber quality traits across multiple generations and environments in upland cotton

Identification of quantitative trait loci (QTL) for fiber quality traits that are stable across multiple generations and environments could facilitate marker-assisted selection for improving cotton strains. In the present study, F₂, F_2:3, and recombinant inbred lines (RILs, F _6:8 ) populations derived from an upland cotton (Gossypium hirsutum L.) cross between strain 0-153, which has excellent fiber quality, and strain sGK9708, a commercial transgenic cultivar, were constructed for QTL tagging of fiber quality. We used 5,742 simple sequence repeat primer pairs to screen for polymorphisms between the two parent strains. Linkage maps of F₂ and RILs were constructed, containing 155 and 190 loci and with a total map distance of 959.4 centimorgans (cM) and 700.9?cM, respectively. We screened fiber quality QTL across multiple generations and environments through composite interval mapping of fiber quality data. Specifically, we studied F₂ and F_2:3 family lines from Anyang (Henan Province) in 2003 and 2004 and RILs in Anyang in 2007 and Anyang, Quzhou (Hebei Province), and Linqing (Shandong Province) in 2008. We identified 50 QTL for fiber quality: 10 for fiber strength, 10 for fiber length, 10 for micronaire, eight for fiber uniformity, and 12 for fiber elongation. Nine of these fiber quality QTL were identified in F₂, F_2:3 and RILs simultaneously. Two QTL for fiber strength on chromosomes C7 and C25 were detected in all three generations and all four environments and explained 16.67?C27.86% and 9.43?C21.36% of the phenotypic variation, respectively. These stable QTL for fiber quality traits could be used for marker assisted selection.     

Gene-environment interactions in complex diseases: genetic models and methods for QTL mapping in multiple half-sib populations

An interval quantitative trait locus (QTL) mapping method for complex polygenic diseases (as binary traits) showing QTL by environment interactions (QEI) was developed for outbred populations on a within-family basis. The main objectives, within the above context, were to investigate selection of genetic models and to compare liability or generalized interval mapping (GIM) and linear regression interval mapping (RIM) methods. Two different genetic models were used: one with main QTL and QEI effects (QEI model) and the other with only a main QTL effect (QTL model). Over 30 types of binary disease data as well as six types of continuous data were simulated and analysed by RIM and GIM. Using table values for significance testing, results show that RIM had an increased false detection rate (FDR) for testing interactions which was attributable to scale effects on the binary scale. GIM did not suffer from a high FDR for testing interactions. The use of empirical thresholds, which effectively means higher thresholds for RIM for testing interactions, could repair this increased FDR for RIM, but such empirical thresholds would have to be derived for each case because the amount of FDR depends on the incidence on the binary scale. RIM still suffered from higher biases (15-100% over- or under-estimation of true values) and high standard errors in QTL variance and location estimates than GIM for QEI models. Hence GIM is recommended for disease QTL mapping with QEI. In the presence of QEI, the model including QEI has more power (20-80% increase) to detect the QTL when the average QTL effect is small (in a situation where the model with a main QTL only is not too powerful). Top-down model selection is proposed in which a full test for QEI is conducted first and then the model is subsequently simplified. Methods and results will be applicable to human, plant and animal QTL mapping experiments.     

Advances in QTL mapping in pigs

Over the past 15 years advances in the porcine genetic linkage map and discovery of useful candidate genes have led to valuable gene and trait information being discovered. Early use of exotic breed crosses and now commercial breed crosses for quantitative trait loci (QTL) scans and candidate gene analyses have led to 110 publications which have identified 1,675 QTL. Additionally, these studies continue to identify genes associated with economically important traits such as growth rate, leanness, feed intake, meat quality, litter size, and disease resistance. A well developed QTL database called PigQTLdb is now as a valuable tool for summarizing and pinpointing in silico regions of interest to researchers. The commercial pig industry is actively incorporating these markers in marker-assisted selection along with traditional performance information to improve traits of economic performance. The long awaited sequencing efforts are also now beginning to provide sequence available for both comparative genomics and large scale single nucleotide polymorphism (SNP) association studies. While these advances are all positive, development of useful new trait families and measurement of new or underlying traits still limits future discoveries. A review of these developments is presented.     

Genetic mapping of a 7R Al tolerance QTL in triticale (x Triticosecale Wittmack)

Triticale (x Triticosecale Wittmack) is a relatively new cereal crop. In Poland, triticale is grown on 12 % of arable land (http://www.stat.gov.pl). There is an increasing interest in its cultivation due to lowered production costs and increased adaptation to adverse environmental conditions. However, it has an insufficient tolerance to the presence of aluminum ions (Al³⁺) in the soil. The number of genes controlling aluminum tolerance in triticale and their chromosomal location is not known. Two F2 mapping biparental populations (MP1 and MP15) segregating for aluminum (Al) tolerance were tested with AFLP, SSR, DArT, and specific PCR markers. Genetic mapping enabled the construction of linkage groups representing chromosomes 7R, 5R and 2B. Obtained linkage groups were common for both mapping populations and mostly included the same markers. Composite interval mapping (CIM) allowed identification of a single QTL that mapped to the 7R chromosome and explained 25 % (MP1) and 36 % (MP15) of phenotypic variation. The B1, B26 and Xscm150 markers were 0.04 cM and 0.02 cM from the maximum of the LOD function in the MP1 and MP15, respectively and were highly associated with aluminum tolerance as indicated by Kruskal–Wallis nonparametric test. Moreover, the molecular markers B1, B26, Xrems1162 and Xscm92, previously associated with the Alt4 locus that encoded an aluminum-activated malate transporter (ScALMT1) that was involved in Al tolerance in rye (Secale cereale) also mapped within QTL. Biochemical analysis of plants represented MP1 and MP15 mapping populations confirmed that the QTL located on 7R chromosome in both mapping populations is responsible for Al tolerance.     

QTL mapping of grain weight in rice and the validation of the QTL qTGW3.2

A recombinant inbred line (RIL) population bred from a cross between a javanica type (cv. D50) and an indica type (cv. HB277) rice was used to map seven quantitative trait loci (QTLs) for thousand grain weight (TGW). The loci were distributed on chromosomes 2, 3, 5, 6, 8 and 10. The chromosome 3 QTL qTGW3.2 was stably expressed over two years, and contributed 9–10% of the phenotypic variance. A residual heterozygous line (RHL) was selected from the RIL population and its selfed progeny was used to fine map qTGW3.2. In this “F₂” population, the QTL explained about 23% of the variance, rising to nearly 33% in the subsequent “F_2:3” generation. The physical location of qTGW3.2 was confined to a ~ 556 kb region flanked by the microsatellite loci RM16162 and RM16194. The region also contains other factors influencing certain yield-related traits, although it is also possible that qTGW3.2 affects these in a pleiotropic fashion.     

Detection of seed dormancy QTL in multiple mapping populations derived from crosses involving novel barley germplasm

Seed dormancy is one of the most important traits in germination process to control malting and pre-harvest sprouting in barley (Hordeum vulgare L.). EST based linkage maps were constructed on seven recombinant inbred (RI) and one doubled haploid (DH) populations derived from crosses including eleven cultivated and one wild barley strains showing the wide range of seed dormancy levels. Seed dormancy of each RI and DH line was estimated from the germination percentage at 5 and 10 weeks post-harvest after-ripening periods in 2003 and 2005. Quantitative trait loci (QTLs) controlling seed dormancy were detected by the composite interval mapping procedure on the RI and DH populations. A total of 38 QTLs clustered around 11 regions were identified on the barley chromosomes except 2H among the eight populations. Several QTL regions detected in the present study were reported on similar positions in the previous QTL studies. The QTL on at the centromeric region of long arm of chromosome 5H was identified in all the RI and DH populations with the different degrees of dormancy depth and period. The responsible gene of the QTL might possess a large allelic variation among the cross combinations, or can be multiple genes located on the same region. The various loci and their different effects in dormancy found in the barley germplasm in the present study enable us to control the practical level of seed dormancy in barley breeding programs. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.