Quantitative trait loci mapping in five new large recombinant inbred line populations of Arabidopsis thaliana genotyped with consensus single-nucleotide polymorphism markers

Quantitative approaches conducted in a single mapping population are limited by the extent of genetic variation distinguishing the parental genotypes. To overcome this limitation and allow a more complete dissection of the genetic architecture of complex traits, we built an integrated set of 15 new large Arabidopsis thaliana recombinant inbred line (RIL) populations optimized for quantitative trait loci (QTL) mapping, having Columbia as a common parent crossed to distant accessions. Here we present 5 of these populations that were validated by investigating three traits: flowering time, rosette size, and seed production as an estimate of fitness. The large number of RILs in each population (between 319 and 377 lines) and the high density of evenly spaced genetic markers scored ensure high power and precision in QTL mapping even under a minimal phenotyping framework. Moreover, the use of common markers across the different maps allows a direct comparison of the QTL detected within the different RIL sets. In addition, we show that following a selective phenotyping strategy by performing QTL analyses on genotypically chosen subsets of 164 RILs (core populations) does not impair the power of detection of QTL with phenotypic contributions >7%.     

Quantitative trait locus mapping using sets of near-isogenic lines: relative power comparisons and technical considerations

 The power to detect QTL using near-isogenic line (NIL) comparisons versus recombinant inbred (RI) populations was assessed. The power to detect QTL was found to be generally greater when using RI populations than when using NIL contrasts. Power to detect QTL with NIL contrasts never exceeded that of RI populations when the number of RI lines is maximized relative to replication of lines for a given number of experimental units. The relative power of NIL contrasts is highest for traits with high heritability and when a gain in precision is realized due to increased replication of entries. Although NIL populations are generally less powerful than RI populations of similar size, some practical considerations may enhance the value of these materials. Availability of NILs allows the genetic effect of a specific chromosome region to be determined by comparing two lines; all RI lines in a population need to be rescored for each new trait even if the effect of a specific chromosome region is suspected. NIL comparisons may allow genetic differences to be detected by visual inspection; genetic effects can only be expressed as means and variances with recombinant inbred populations. In summary, RI populations generally, and in some cases, substantially, provide better power for QTL detection than NIL comparisons. Practical considerations, however, indicate that many factors need to be considered when choosing a population structure to meet an experimental objective. Received: 12 December 1996 / Accepted: 21 March 1997     

Bay-0 x Shahdara recombinant inbred line population: a powerful tool for the genetic dissection of complex traits in Arabidopsis

Natural genetic variation in Arabidopsis is considerable, but has not yet been used extensively as a source of variants to identify new genes of interest. From the cross between two genetically distant ecotypes, Bay-0 and Shahdara, we generated a Recombinant Inbred Line (RIL) population dedicated to Quantitative Trait Locus (QTL) mapping. A set of 38 physically anchored microsatellite markers was created to construct a robust genetic map from the 420 F6 lines. These markers, evenly distributed throughout the five chromosomes, revealed a remarkable equilibrium in the segregation of parental alleles in the genome. As a model character, we have analysed the genetic basis of variation in flowering time in two different environments. The simultaneous mapping of both large- and small-effect QTLs responsible for this variation explained 90% of the total genotypic variance. Two of the detected QTLs colocalize very precisely with FRIGIDA and FLOWERING LOCUS C genes; we provide information on the polymorphism of genes confirming this hypothesis. Another QTL maps in a region where no QTL had been found previously for this trait. This confirms the accuracy of QTL detection using the Bay-0 x Shahdara RIL population, which constitutes the largest in size available so far in Arabidopsis. As an alternative to mutant analysis, this population represents a powerful tool which is currently being used to undertake the genetic dissection of complex metabolic pathways.     

Variation in seed dormancy quantitative trait loci in Arabidopsis thaliana originating from one site

A Quantitative Trait Locus (QTL) analysis was performed using two novel Recombinant Inbred Line (RIL) populations, derived from the progeny between two Arabidopsis thaliana genotypes collected at the same site in Kyoto (Japan) crossed with the reference laboratory strain Landsberg erecta (Ler). We used these two RIL populations to determine the genetic basis of seed dormancy and flowering time, which are assumed to be the main traits controlling life history variation in Arabidopsis. The analysis revealed quantitative variation for seed dormancy that is associated with allelic variation at the seed dormancy QTL DOG1 (for Delay Of Germination 1) in one population and at DOG6 in both. These DOG QTL have been previously identified using mapping populations derived from accessions collected at different sites around the world. Genetic variation within a population may enhance its ability to respond accurately to variation within and between seasons. In contrast, variation for flowering time, which also segregated within each mapping population, is mainly governed by the same QTL.     

Precise mapping Fhb5, a major QTL conditioning resistance to Fusarium infection in bread wheat (Triticum aestivum L.)

Qfhi.nau-5A is a major quantitative trait locus (QTL) against Fusarium graminearum infection in the resistant wheat germplasm Wangshuibai. Genetic analysis using BC(3)F(2) and BC(4)F(2) populations, derived from selfing two near-isogenic lines (NIL) heterozygous at Qfhi.nau-5A that were developed, respectively, with Mianyang 99-323 and PH691 as the recurrent parent, showed that Qfhi.nau-5A inherited like a single dominant gene. This QTL was thus designated as Fhb5. To fine map it, these two backcross populations and a recombinant inbred line (RIL) population derived from Nanda2419?×?Wangshuibai were screened for recombinants occurring between its two flanking markers Xbarc56 and Xbarc100. Nineteen NIL recombinants were identified from the two backcross populations and nine from the RIL population. In the RIL recombinant selection process, selection against Fhb4 present in the RIL population was incorporated. Genotyping these recombinant lines with ten markers mapping to the Xbarc56-Xbarc100 interval revealed four types of Mianyang 99-323-derived NIL recombinants, three types of PH691-derived NIL recombinants, and four types of RIL recombinants. In different field trials, the percentage of infected spikes of these lines displayed a distinct two-peak distribution. The more resistant class had over 55% less infection than the susceptible class. Common to these resistant genotypes, the 0.3-cM interval flanked by Xgwm304 and Xgwm415 or one of these two loci was derived from Wangshuibai, while none of the susceptible recombinants had Wangshuibai chromatin in this interval. This interval harboring Fhb5 was mapped to the pericentromeric C-5AS3-0.75 bin through deletion bin mapping. The precise localization of Fhb5 will facilitate its utilization in marker-assisted wheat breeding programs.     

Genetic analysis of the physiological responses to low boron stress in Arabidopsis thaliana

Boron (B) is an essential micronutrient for higher plants. There is wide genetic variation in the response to B deficiency among plant species and cultivars. The objective of this study was to identify quantitative trait loci (QTL) that control B efficiency in natural Arabidopsis accessions. The B efficiency coefficient (BEC) and seed yield under low B conditions (SYLB) were investigated by solution culture in two separate experiments in an Arabidopsis recombinant inbred line (RIL) population. Both of the traits studied exhibited high transgressive variation in the RIL population, and, in total, five and three QTL were identified for BEC and SYLB, respectively. Three of the five QTL, including the QTL, AtBE1-2, that has a large effect on the BEC, were found at the interval of the corresponding QTL for SYLB in both experiments. The close genetic relationship between BEC and SYLB was further confirmed by conditional QTL mapping in the RIL population and unconditional QTL mapping in an AtBE1-2-segregated F(2) population. Epistatic interactions for the tested traits were analysed, and were found to be widespread in the detected QTL of Arabidopsis in the RIL population. Comparison of the QTL interval for B efficiency with reported B-related genes showed that 10 B-related genes, together with one BOR1 homolog (BOR5, At1g74810) were located in the QTL region of AtBE1-2. These results suggest that natural variation in B efficiency in Arabidopsis has a complex molecular basis. They also provide a basis for fine mapping and cloning of the B-efficiency genes, with the ultimate aim of discovering the physiological mechanism of action of the genes.     

Epistasis and genotype-environment interaction for quantitative trait loci affecting flowering time in Arabidopsis thaliana

A major goal of evolutionary biology is to understand the genetic architecture of the complex quantitative traits that may lead to adaptations in natural populations. Of particular relevance is the evaluation of the frequency and magnitude of epistasis (gene–gene and gene–environment interaction) as it plays a controversial role in models of adaptation within and among populations. Here, we explore the genetic basis of flowering time in Arabidopsis thaliana using a series of quantitative trait loci (QTL) mapping experiments with two recombinant inbred line (RIL) mapping populations [Columbia (Col) x Landsberg erecta (Ler), Ler x Cape Verde Islands (Cvi)]. We focus on the response of RILs to a series of environmental conditions including drought stress, leaf damage, and apical damage. These data were explicitly evaluated for the presence of epistasis using Bayesian based multiple-QTL genome scans. Overall, we mapped fourteen QTL affecting flowering time. We detected two significant QTL–QTL interactions and several QTL–environment interactions for flowering time in the Ler x Cvi population. QTL–environment interactions were due to environmentally induced changes in the magnitude of QTL effects and their interactions across environments – we did not detect antagonistic pleiotropy. We found no evidence for QTL interactions in the Ler x Col population. We evaluate these results in the context of several other studies of flowering time in Arabidopsis thaliana and adaptive evolution in natural populations.     

Genetics of drought adaptation in Arabidopsis thaliana II. QTL analysis of a new mapping population, KAS-1 x TSU-1

Despite compelling evidence that adaptation to local climate is common in plant populations, little is known about the evolutionary genetics of traits that contribute to climatic adaptation. A screen of natural accessions of Arabidopsis thaliana revealed Tsu-1 and Kas-1 to be opposite extremes for water-use efficiency and climate at collection sites for these accessions differs greatly. To provide a tool to understand the genetic basis of this putative adaptation, Kas-1 and Tsu-1 were reciprocally crossed to create a new mapping population. Analysis of F(3) families showed segregating variation in both delta(13)C and transpiration rate, and as expected these traits had a negative genetic correlation (r(g)=- 0.3). 346 RILs, 148 with Kas-1 cytoplasm and 198 with Tsu-1 cytoplasm, were advanced to the F(9) and genotyped using 48 microsatellites and 55 SNPs for a total of 103 markers. This mapping population was used for QTL analysis of delta(13)C using F(9) RIL means. Analysis of this reciprocal cross showed a large effect of cytoplasmic background, as well as two QTL for delta(13)C. The Kas-1 x Tsu-1 mapping population provides a powerful new resource for mapping QTL underlying natural variation and for dissecting the genetic basis of water-use efficiency differences.     

Dissection and fine mapping of a major QTL for preharvest sprouting resistance in white wheat Rio Blanco

Preharvest sprouting (PHS) is a major constraint to white wheat production. Previously, we mapped quantitative trait loci (QTL) for PHS resistance in white wheat by using a recombinant inbred line (RIL) population derived from the cross Rio Blanco/NW97S186. One QTL, QPhs.pseru-3A, showed a major effect on PHS resistance, and three simple sequence repeat (SSR) markers were mapped in the QTL region. To determine the flanking markers for the QTL and narrow down the QTL to a smaller chromosome region, we developed a new fine mapping population of 1,874 secondary segregating F₂ plants by selfing an F6 RIL (RIL25) that was heterozygous in the three SSR marker loci. Segregation of PHS resistance in the population fitted monogenic inheritance. An additive effect of the QTL played a major role on PHS resistance, but a dominant effect was also observed. Fifty-six recombinants among the three SSR markers were identified in the population and selfed to produce homozygous recombinants or QTL near-isogenic lines (NIL). PHS evaluation of the recombinants delineated the QTL in the region close to Xbarc57 flanked by Xbarc321 and Xbarc12. To saturate the QTL region, 11 amplified fragment length polymorphism (AFLP) markers were mapped in the QTL region with 7 AFLP co-segregated with Xbarc57 by using the NIL population. Dissection of the QTL as a Mendelian factor and saturation of the QTL region with additional markers created a solid foundation for positional cloning of the major QTL.     

Quantitative trait loci and candidate genes underlying genotype by environment interaction in the response of Arabidopsis thaliana to drought

Drought stress was imposed on two sets of Arabidopsis thaliana genotypes grown in sand under short‐day conditions and analysed for several shoot and root growth traits. The response to drought was assessed for quantitative trait locus (QTL) mapping in a genetically diverse set of Arabidopsis accessions using genome‐wide association (GWA) mapping, and conventional linkage analysis of a recombinant inbred line (RIL) population. Results showed significant genotype by environment interaction (G×E) for all traits in response to different watering regimes. For the RIL population, the observed G×E was reflected in 17 QTL by environment interactions (Q×E), while 17 additional QTLs were mapped not showing Q×E. GWA mapping identified 58 single nucleotide polymorphism (SNPs) associated with loci displaying Q×E and an additional 16 SNPs associated with loci not showing Q×E. Many candidate genes potentially underlying these loci were suggested. The genes for RPS3C and YLS7 were found to contain conserved amino acid differences when comparing Arabidopsis accessions with strongly contrasting drought response phenotypes, further supporting their candidacy. One of these candidate genes co‐located with a QTL mapped in the RIL population.     

Identification of QTL for flag leaf length in common wheat and their pleiotropic effects

Leaf size is an important factor contributing to the photosynthetic capability of wheat plants. It also significantly affects various agronomic traits. In particular, the flag leaves contribute significantly to grain yield in wheat. A recombinant inbred line (RIL) population developed between varieties with significant differences in flag leaf traits was used to map quantitative trait loci (QTL) of flag leaf length (FLL) and to evaluate its pleiotropic effects on five yield-related traits, including spike length (SL), spikelet number per spike (SPN), kernel number per spike (KN), kernel length (KL), and thousand-kernel weight (TKW). Two additional RIL populations were used to validate the detected QTL and reveal the relationships in different genetic backgrounds. Using the diversity arrays technology (DArT) genetic linkage map, three major QTL for FLL were detected, with single QTL in different environments explaining 8.6–23.3% of the phenotypic variation. All the QTL were detected in at least four environments, and validated in two related populations based on the designed primers. These QTL and the newly developed primers are expected to be valuable for fine mapping and marker-assisted selection in wheat breeding programs.     

Quantitative trait locus analysis of a recombinant inbred line population derived from a Lycopersicon esculentum x Lycopersicon cheesmanii cross

Quantitative trait loci influencing fruit traits were identified by restriction fragment length polymorphism (RFLP) analysis in a population of recombinant inbred lines (RIL) derived from a cross of the cultivated tomato, Lycopersicon esculentum with a related wild species Lycopersicon cheesmanii. One hundred thirty-two polymorphic RFLP loci spaced throughout the tomato genome were scored for 97 F₈ RIL families. Fruit weight and soluble solids were measured in replicated trials during 1991 and 1992. Seed weight was measured in 1992. Significant (P<0.01 level) quantitative trait locus (QTL) associations of marker loci were identified for each trait. A total of 73 significant marker locus-trait associations were detected for the three traits measured. Fifty-three of these associations were for fruit weight and soluble solids, many of which involved marker loci signficantly associated with both traits. QTL with large effects on all three traits were detected on chromosome 6. Greater homozygosity at many loci in the RIL population as compared to F₂ populations and greater genomic coverage resulted in increased precision in the estimation of QTL effects, and large proportions of the total phenotypic variance were explained by marker class variation at significant marker loci for many traits. The RIL population was effective in detecting and discriminating among QTL for these traits previously identified in other investigations despite skewed segregation ratios at many marker loci. Large additive effects were measured at significant marker loci. Lower fruit weight, higher soluble solids, and lower seed weight were generally associated with RFLP alleles from theL. cheesmanii parent.     

Genetic analysis and QTL mapping of the seed hardness trait in a black common bean (<Emphasis Type="Italic">Phaseolus vulgaris</Emphasis>) recombinant inbred line (RIL) population

Seed hardness trait has a profound impact on cooking time and canning quality in dry beans. This study aims to identify the unknown genetic factors and associated molecular markers to better understand and tag this trait. An F_2:7 recombinant inbred line (RIL) population was derived from a cross between the hard and soft seeded black bean parents (H68-4 and BK04-001). Eighty-five RILs and the parental lines were grown at two locations in southern Manitoba during years 2014–2016. Seed samples were harvested manually at maturity to test for seed hardness traits. The hydration capacity and stone seed count were estimated by soaking the seeds overnight at room temperature following AACC method 56-35.01. Seed samples from 2016 tests were also cooked to determine effect of seed hardness on cooking quality. For mapping of genomic regions contributing to the traits, the RIL population was genotyped using the genotype by sequencing (GBS) approach. The QTL mapping revealed that in addition to the major QTL on chromosome 7 at a genomic location previously reported to affect seed-hydration, two novel QTL with significant effects were also detected on chromosomes 1 and 2. In addition, a major QTL affecting the visual appeal of cooked bean was mapped on chromosome 4. This multi-year-site study shows that despite large environmental effects, seed hardness is an oligo-genic and highly heritable trait, which is inherited independently of the cooking quality scored as visual appeal of cooked beans. The identification of the QTLs and development of SNP markers associated with seed hardness can be applied for common bean variety improvement and genetic exploitation of these traits.     

Multiple-line cross quantitative trait locus mapping in sugar beet (Beta vulgaris L.)

Linkage mapping based on multiple-line crosses is a promising strategy for mapping quantitative trait loci (QTL) underlying important agronomic traits. The main goal of this survey was to study the advantages of QTL mapping across versus within biparental populations using experimental data from three connected sugar beet (Beta vulgaris L.) populations evaluated for beet yield and potassium and sodium content. For the combined analysis across populations, we used two approaches for cofactor selection. In Model A, we assumed identical cofactors for every segregating population. In contrast, in Model B we selected cofactors specific for every segregating population. Model A performed better than Model B with respect to the number of QTL detected and the total proportion of phenotypic variance explained. The QTL analyses across populations revealed a substantially higher number of QTL compared to the analyses of single biparental populations. This clearly emphasizes the potential to increase QTL detection power with a joint analysis across biparental populations.     

Quantitative trait loci for inflorescence development in Arabidopsis thaliana

Variation in inflorescence development patterns is a central factor in the evolutionary ecology of plants. The genetic architectures of 13 traits associated with inflorescence developmental timing, architecture, rosette morphology, and fitness were investigated in Arabidopsis thaliana, a model plant system. There is substantial naturally occurring genetic variation for inflorescence development traits, with broad sense heritabilities computed from 21 Arabidopsis ecotypes ranging from 0.134 to 0.772. Genetic correlations are significant for most (64/78) pairs of traits, suggesting either pleiotropy or tight linkage among loci. Quantitative trait locus (QTL) mapping indicates 47 and 63 QTL for inflorescence developmental traits in Ler x Col and Cvi x Ler recombinant inbred mapping populations, respectively. Several QTL associated with different developmental traits map to the same Arabidopsis chromosomal regions, in agreement with the strong genetic correlations observed. Epistasis among QTL was observed only in the Cvi x Ler population, and only between regions on chromosomes 1 and 5. Examination of the completed Arabidopsis genome sequence in three QTL regions revealed between 375 and 783 genes per region. Previously identified flowering time, inflorescence architecture, floral meristem identity, and hormone signaling genes represent some of the many candidate genes in these regions.     

Genetic Architecture of Natural Variation of Telomere Length in Arabidopsis thaliana

Telomeres represent the repetitive sequences that cap chromosome ends and are essential for their protection. Telomere length is known to be highly heritable and is derived from a homeostatic balance between telomeric lengthening and shortening activities. Specific loci that form the genetic framework underlying telomere length homeostasis, however, are not well understood. To investigate the extent of natural variation of telomere length in Arabidopsis thaliana, we examined 229 worldwide accessions by terminal restriction fragment analysis. The results showed a wide range of telomere lengths that are specific to individual accessions. To identify loci that are responsible for this variation, we adopted a quantitative trait loci (QTL) mapping approach with multiple recombinant inbred line (RIL) populations. A doubled haploid RIL population was first produced using centromere-mediated genome elimination between accessions with long (Pro-0) and intermediate (Col-0) telomere lengths. Composite interval mapping analysis of this population along with two established RIL populations (Ler-2/Cvi-0 and Est-1/Col-0) revealed a number of shared and unique QTL. QTL detected in the Ler-2/Cvi-0 population were examined using near isogenic lines that confirmed causative regions on chromosomes 1 and 2. In conclusion, this work describes the extent of natural variation of telomere length in A. thaliana, identifies a network of QTL that influence telomere length homeostasis, examines telomere length dynamics in plants with hybrid backgrounds, and shows the effects of two identified regions on telomere length regulation.     

Power and false-positive rate in QTL detection with near-isogenic line libraries

Libraries of near-isogenic lines (NILs) were used for quantitative trait locus (QTL) detection in model species and economically important crops. The experimental design and genetic architecture of the considered traits determine the statistical properties of QTL detection. The objectives of our simulation study were to (i) investigate the population sizes required to develop NIL libraries in barley and maize, (ii) compare NIL libraries with nonoverlapping and overlapping donor segments and (iii) study the number of QTLs and the size of their effects with respect to the power and the false-positive rate of QTL detection. In barley, the development of NIL libraries with target segment lengths of 10 c and marker distances of 5 c was possible using a BC₃S₂ backcrossing scheme and population sizes of 140. In maize, population sizes larger than 200 were required. Selection for the recipient parent genome at markers flanking the target segments with distances between 5 and 10 c was required for an efficient control of the false-positive rate. NIL libraries with nonoverlapping donor chromosome segments had a greater power of QTL detection and a smaller false-positive rate than libraries with overlapping segments. Major genes explaining 30% of the genotypic difference between the donor and recipient were successfully detected even with low heritabilities of 0.5, whereas for minor genes explaining 5 !or 10%, high heritabilities of 0.8 or 0.9 were required. The presented results can assist geneticists and breeders in the efficient development of NIL libraries for QTL detection.     

A multiparental cross population for mapping QTL for agronomic traits in durum wheat (Triticum turgidum ssp. durum)

Multiparental cross designs for mapping quantitative trait loci (QTL) provide an efficient alternative to biparental populations because of their broader genetic basis and potentially higher mapping resolution. We describe the development and deployment of a recombinant inbred line (RIL) population in durum wheat (Triticum turgidum ssp. durum) obtained by crossing four elite cultivars. A linkage map spanning 2664 cM and including 7594 single nucleotide polymorphisms (SNPs) was produced by genotyping 338 RILs. QTL analysis was carried out by both interval mapping on founder haplotype probabilities and SNP bi‐allelic tests for heading date and maturity date, plant height and grain yield from four field experiments. Sixteen QTL were identified across environments and detection methods, including two yield QTL on chromosomes 2BL and 7AS, with the former mapped independently from the photoperiod response gene Ppd‐B1, while the latter overlapped with the vernalization locus VRN‐A3. Additionally, 21 QTL with environment‐specific effects were found. Our results indicated a prevalence of environment‐specific QTL with relatively small effect on the control of grain yield. For all traits, functionally different QTL alleles in terms of direction and size of genetic effect were distributed among parents. We showed that QTL results based on founder haplotypes closely matched functional alleles at known heading date loci. Despite the four founders, only 2.1 different functional haplotypes were estimated per QTL, on average. This durum wheat population provides a mapping resource for detailed genetic dissection of agronomic traits in an elite background typical of breeding programmes.     

Identification of quantitative trait loci controlling symptom development during viral infection in Arabidopsis thaliana

In compatible interactions between plants and viruses that result in systemic infection, symptom development is a major phenotypic trait. However, host determinants governing this trait are mostly unknown, and the mechanisms underlying it are still poorly understood. In a previous study on the Arabidopsis thaliana-Plum pox virus (PPV) pathosystem, we showed a large degree of variation in symptom development among susceptible accessions. In particular, Cvi-1 (Cape Verde islands) accumulates viral particules but remains symptomless, Col-0 (Columbia) sometimes shows weak symptoms compared with Ler (Landsberg erecta), which always shows severe symptoms. Genetic analyses of Col x Ler and Cvi x Ler F2 and recombinant inbred line (RIL) populations suggested that symptom development as well as viral accumulation traits are polygenic and quantitative. Three of the symptom quantitative trait loci (QTL) identified could be confirmed in near-isogenic lines, including PSI1 (PPV symptom induction 1), which was identified on the distal part of chromosome 1 in both RIL populations. With respect to viral accumulation, several factors have been detected and, interestingly, in the Col x Ler population, two out of three viral accumulation QTL colocalized with loci controlling symptom development, although correlation analysis showed weak linearity between symptom severity and virus accumulation. In addition, in the Cvi x Ler RIL population, a digenic recessive determinant controlling PPV infection was identified.