Development of a next-generation NIL library in Arabidopsis thaliana for dissecting complex traits
© Fletcher et al.; licensee BioMed Central Ltd. 2013
Received: 11 March 2013
Accepted: 20 September 2013
Published: 25 September 2013
The identification of the loci and specific alleles underlying variation in quantitative traits is an important goal for evolutionary biologists and breeders. Despite major advancements in genomics technology, moving from QTL to causal alleles remains a major challenge in genetics research. Near-isogenic lines are the ideal raw material for QTL validation, refinement of QTL location and, ultimately, gene discovery.
In this study, a population of 75 Arabidopsis thaliana near-isogenic lines was developed from an existing recombinant inbred line (RIL) population derived from a cross between physiologically divergent accessions Kas-1 and Tsu-1. First, a novel algorithm was developed to utilize genome-wide marker data in selecting RILs fully isogenic to Kas-1 for a single chromosome. Seven such RILs were used in 2 generations of crossing to Tsu-1 to create BC1 seed. BC1 plants were genotyped with SSR markers so that lines could be selected that carried Kas-1 introgressions, resulting in a population carrying chromosomal introgressions spanning the genome. BC1 lines were genotyped with 48 genome-wide SSRs to identify lines with a targeted Kas-1 introgression and the fewest genomic introgressions elsewhere. 75 such lines were selected and genotyped at an additional 41 SNP loci and another 930 tags using 2b-RAD genotyping by sequencing. The final population carried an average of 1.35 homozygous and 2.49 heterozygous introgressions per line with average introgression sizes of 5.32 and 5.16 Mb, respectively. In a simple case study, we demonstrate the advantage of maintaining heterozygotes in our library whereby fine-mapping efforts are conducted simply by self-pollination. Crossovers in the heterozygous interval during this single selfing generation break the introgression into smaller, homozygous fragments (sub-NILs). Additionally, we utilize a homozygous NIL for validation of a QTL underlying stomatal conductance, a low heritability trait.
The present results introduce a new and valuable resource to the Brassicaceae research community that enables rapid fine-mapping of candidate loci in parallel with QTL validation. These attributes along with dense marker coverage and genome-wide chromosomal introgressions make this population an ideal starting point for discovery of genes underlying important complex traits of agricultural and ecological significance.
Keywords2b-RAD Fine-mapping Quantitative trait loci Stomatal conductance
Linkage mapping of QTL is a common statistical approach in plant genetics where recombinant populations generated from crosses between inbred parent lines are used, in combination with molecular markers, to identify loci associated with variation in continuously distributed traits [1–8]. Mapping populations common to QTL analyses are many and include doubled haploids (DH), F2, backcross, advanced intercross, nested association mapping and RILs. Mapping QTL for complex traits is now routine, with the typical output being QTL spanning large confidence intervals encompassing many (hundreds or more) possible causal genes .
The steps following QTL identification frequently involve functional validation of the QTL, and refinement of location (fine-mapping) towards the goal of identification of a causal gene – the major challenge in quantitative genetics today . One of the most common approaches for accomplishing these objectives is through the development and phenotypic characterization of NILs . The generation and phenotyping of NILs is considered a laborious and time consuming process, but the robust design leads to a minimal false positive rate.
NILs are lines containing a single or small number of genomic introgressions from a donor parent in a different and otherwise homogeneous genomic background. By homogenizing all genetic factors outside of the focal genomic region, the true effect of the QTL on the phenotype can be estimated relative to the line into which the introgression was introduced (i.e. void of the chromosomal introgression) . In addition to the simplification of genetic analyses, NILs are considered genetically 'immortal’  which allows for replicated experiments across multiple environments resulting in more accurate estimates of effect size for complex traits. NILs have proven to be an effective resource for QTL validation and a logical starting point for the creation of fine-mapping populations [14–21].
Creation of a single near-isogenic line generally starts by crossing a line carrying the targeted QTL region to one of the parental lines of the population, thus creating a backcross population. Genome-wide genotyping of the backcross progeny is performed to identify recombination events allowing for selection of progeny which carry the target chromosomal introgression derived from the donor and recurrent parent genome elsewhere. Subsequent generations of self-pollination (selfing) are normally required to achieve homozygosity of the introgressed region and the process can take several backcrossing cycles to produce a NIL carrying an introgression of acceptable size and genomic location. An alternative approach has been the use of heterogeneous inbred families (HIFs) where NILs are selected from incompletely inbred lines which still harbour a small amount of heterozygosity at random intervals across the genome [22, 23]. Analysis of a HIF population with molecular markers allows for the selection of lines heterozygous at a candidate genomic location, which in combination with further selfing and genotyping, enables selection of NILs derived from several heterogeneous genetic backgrounds. Producing NILs with smaller introgressions requires greater effort. Large populations are needed to break up small chromosomal segments, and high-density genotyping is required to discover them.
A NIL library is a family of near-isogenic lines where each line carries a different donor parent fragment and the population carries introgressions spanning the entire genome . A NIL library is an ideal starting point for QTL validation, especially in cases where the library is derived from parent lines for which an immortal recombinant population (i.e. RILs, DH, etc.) already exists. In this case, QTL identified via traditional linkage mapping experiments performed on the mapping population can be immediately tested by selecting NIL(s) representing the QTL introgression and testing them for a phenotypic effect relative to the wild type recurrent parent. NIL libraries are also valuable starting material for fine-mapping QTL through the creation of sub-NILs , recombinant lines in which the original NIL introgression is broken into smaller genomic fragments. In this case, a candidate NIL is backcrossed to the recurrent parent and the progeny are genotyped using markers specific to the introgression region so that individuals carrying genomic fragments spanning the length of the original introgression can be identified. Subsequent phenotyping of the sub-NILs provides finer resolution of the region controlling the trait of interest, effectively narrowing the list of possible causal genes.
Several NIL populations are currently available to the Arabidopsis research community. Koumproglou et al. , using 31 simple sequence repeat (SSR) markers, created a population of Chromosome Substitution Strains by replacing chromosomes from the accession Columbia (Col-0) with homologous chromosomes from the accessions Landsberg erecta (Ler) and Niederzenz (Nd). Additionally, a population of more traditional NILs were created in a systematic approach where increasing lengths of chromosomal introgressions were introduced from Ler into the Col-0 background. Keurentjes et al.  generated a population of 92 NILs carrying genome-wide chromosomal introgression from the accession Cape Verde Islands (Cvi) into the Ler background. Selections were made from the genotyped RIL mapping population described by Alonso-Blanco et al.  and used in backcrosses to create the NIL library. The RIL population has been mapped for QTL underlying flowering time and carbon isotope ratio (δ13C) , recombination frequency , seed germination , seed mineral concentration  and fructose sensitivity . The same 321 AFLP (Amplified Fragment Length Polymorphism) markers used to build the RIL map were used in the NIL breeding scheme. Finally, Torjek et al.  created a population of 140 reciprocal NILs from the accessions Col-0 and C24 (78 NILs in the Col-0 background and 62 lines in the C24 background) utilizing a total of 125 markers . This NIL library has been used in subsequent studies of epistasis  and heterosis .
Here we report the development of a new population of 75 NILs constituting genome-wide chromosomal introgressions. The NIL population exploited inbred lines selected from the RIL population described in McKay et al.  as the starting material for backcrossing. Briefly, the RIL population is derived from a cross between the A. thaliana ecotypes Tsu-1 (CS1640), an accession originating from Tsushima, Japan and Kas-1 (CS903), an accession originating from Kashmir. These sites of collection are among the wettest and driest habitats, respectively, in the A. thaliana species range and the accessions differ in several aspects of drought physiology [39, 40]. Recombinant populations derived from these diverse accessions will therefore segregate alleles underlying variation in these physiological traits, providing a powerful resource for identifying functional genes.
We developed a population of 75 Arabidopsis thaliana NILs containing both homozygous and heterozygous introgressions, enabling simultaneous pursuit of QTL validation and fine-mapping. Genotyping the population with over 1,000 molecular markers has provided us with excellent resolution on the total number of introgressions existing in each NIL as well as their location and length. It is the most densely genotyped NIL population developed thus far by more than 3-fold. The utility of the NIL library is demonstrated in a simple case study where, in a single generation, we utilize a homozygous NIL to validate and localize a QTL for a low heritability physiological trait (g0; night-time stomatal conductance) while concurrently selfing heterozygous selections to create sub-NILs for further fine-mapping.
Marker-assisted NIL breeding program
Polymorphisms detected between Tsu-1 and Kas-1 by 2b-RAD genotyping
Restriction site–associated DNA (RAD) tag sequencing reduces genome complexity by focusing only on DNA flanking the recognition sites of the selected restriction endonuclease . The RAD method used in this study, described in  is a simple and effective means of discovering a large number of SNPs unique to the study population, avoiding the ascertainment bias associated with SNPs discovered via population surveys . The 2b-RAD method utilizes the type IIB restriction enzyme, AlfI, which operates by cleaving DNA both upstream and downstream of the recognition site. The resulting tags are uniform in length, making them ideal for amplification and sequencing on next-generation platforms. Following digestion, tags were labelled with sample-specific oligonucleotide barcodes for multiplexed sequencing. Finally, reads were quality filtered and aligned to a collection of AlfI sites in the Col-0 Arabidopsis reference genome (TAIR9) in order to assign a physical location to each SNP.
Genomics of chromosomal introgressions in the NIL population and the added value of increased marker resolution
Case study: utilizing selections from the NIL library for QTL validation and sub-NIL development
To demonstrate the value of this new resource, we analyzed the RIL population  for QTL for night-time leaf conductance (g0). g0 is a low-heritability, quantitative trait that is important for plant-water relations and mineral nutrition. While the adaptive value of g0 has yet to be fully understood, incomplete stomatal closure during the night can lead to substantial transpirational water loss . Variation in this trait has been found among and within species, and it correlates with some daytime gas-exchange traits such as water-use efficiency (the ratio of CO2 assimilation to transpiration) . Estimates of transpiration have been found to be particularly sensitive to g0, making it an interesting candidate for studies on the physiology and genetics of plant drought adaptation. In view of that, intraspecific variation in observed g0 has been found to have the largest effect on transpiration across a species’ native habitat (Bauerle, unpublished observations).
Results of QTL validation experiment comparing NIL g 0 values with Tsu-1
(mmol m-2 s-1)
NIL - (Tsu-1)
Nearly 1,500 genes are predicted to lie within the region spanning physical positions 505,086 to 5,273,972 of chromosome one. We have assembled a list of candidate genes based upon hits to gene ontology (GO) terms relevant to stomatal conductance: abscisic acid (ABA), stomata and water [see Additional file 7].
Maintenance of homozygous and heterozygous NILs facilitates simultaneous QTL validation and fine-mapping efforts
Near-isogenic lines remain the ideal starting material for validation of QTL as well as breeding schemes designed for fine-mapping with the end goal being the identification of candidate genes [48–51]. QTL validation is relatively straightforward and consists simply of phenotyping NILs with introgressions at the region of interest for the trait of interest. Creation of a suitable population for fine mapping is not as straightforward and is normally a three-generation process that starts with a cross between an inbred NIL and the recurrent parent. This is typically followed by a generation of self-pollination to allow for recombination in the introgression region. The seed harvested from these self-pollinated plants can then be genotyped with markers specific to the region so that homozygous sub-NILs can be identified. The process is fairly straightforward and inexpensive in the context of physical resource, but there is a time cost of at least 3 generations (equivalent to a minimum of 18 weeks).
Our case study illustrates the advantages of maintaining both homozygotes and heterozygotes in the NIL population, combining the benefits of traditional homozygous NILs with the advantages of HIFs [22, 23]. For example, measuring g0 on the homozygous NILs provided strong evidence for the presence of the QTL in a single generation, thus avoiding the process of generating homozygous lines that would be necessary in HIF populations. These results provided a better estimate of the QTL effect size relative to the results derived from our QTL mapping approach and have justified further investments in fine-mapping using heterozygous NILs. This emphasizes the power NILs create by isolating the genetic factors controlling a phenotype to a single locus as there were other loci worthy of consideration as contributors to variation in g0 in the RIL population. Analysis of the genes predicted to lie within this interval revealed a majority of them had GO annotations related to ABA, the major signalling molecule in stomatal regulation [52–54], but examination of the entire region with the AmiGO enrichment analysis tool  found it was not significantly enriched for ABA genes. Inspection of the physical location of these ABA-associated candidates reveals that they are clustered in a 1.2 Mb interval (At 1 physical interval: 712,473-1,894,148) which represents a relatively small portion of the 4.8 Mb introgression tested, thus providing an interesting focal region during fine-mapping of the g0 phenotype.
With regards to fine-mapping, selfing a heterozygous NIL selection from the population yielded several sub-NILs suitable for phenotyping or additional genotyping in future generations, an attribute common with HIF populations and advantageous over traditional NILs. This was accomplished using a modest population size of BC1S3 plants (n = 286) and the interval could be narrowed down further through genotyping at a higher number of loci and increasing the population size . Regardless, in a six-week period we have identified a population encompassing recombinants in the 4.8 Mb region identified as causal during the QTL validation experiment, translating to a 3-fold change in total time versus a breeding scheme utilizing inbred NILs.
2b-RAD is an efficient method for dense genotyping of recombinant populations
Arabidopsis thaliana recently celebrated its 25th anniversary as a model organism and now stands alone as the most thoroughly studied plant species on record http://www.arabidopsis.org/, . Recent efforts are producing comprehensive polymorphism databases (http://www.arabidopsis.org/, http://signal.salk.edu/cgi-bin/AtSFP). To interpret the significance and functional consequences of this natural variation, we need to understand the multivariate phenotypic consequences of these variants. NIL libraries, mutants and complementation studies are the tools required for this mechanistic understanding.
The 2b-RAD method added an additional 930 high confidence genotypes to our map providing a level of resolution not yet achieved in any of the Arabidopsis NIL populations described to date. The value of these additional markers is obvious as we compare the coarse and dense maps. The discovery of an additional 129 introgressions is clearly important when making selections for QTL validation. For instance, three additional homozygous introgressions were discovered in KT154_2_3, changing the estimate from one to four. This is a clear illustration of the risks associated with utilizing NILs genotyped at low density in experiments aimed at QTL validation. These offsite introgressions may have effects on the phenotype of interest, potentially resulting in erroneous or uncertain conclusions regarding the QTL effect size and location.
The Kas-1 × Tsu-1 RIL and NIL populations are a valuable resource for research on the genetics of drought adaptation in the Brassicaceae
Substantial variation for several traits relevant to drought adaptation have been observed in the Kas-1 × Tsu-1 RIL population including δ13C, leaf water content, instantaneous transpiration rate, flowering time, abscisic acid content and root mass , unpublished results]. Accordingly, the NIL population described herein is expected to vary for the same traits, providing a powerful resource for moving from QTL, encompassing thousands of genes, discovered in the RIL population towards a smaller list of putative functional candidates.
No other plant species has been more studied or characterized than Arabidopsis thaliana. A high degree of sequence collinearity between it and members of the agriculturally significant Brassica genus was discovered over a decade ago . Similar levels of synteny have been found in comparisons with other taxa in the Brassicaceae [59–62]. These results suggest that translational genomics, that is utilizing basic research findings in model organisms to answer practical research questions in species of higher economic value or importance [63, 64], could be a viable avenue in understanding complex traits. In this regard, we suggest the Kas-1 × Tsu-1 populations as the ideal starting point for basic research on the genetics and genomics of drought adaptation.
We have developed a population of 75 NILs that provides genetic resources for fine-mapping QTL as well as QTL corroboration. The high marker density used to construct the population provides a level of resolution not yet seen in a NIL population, thus minimizing ambiguity in fine-mapping and QTL validation studies caused by unidentified chromosomal introgressions elsewhere in the genome. The unique variation that exists between the parents used to construct this resource provides a valuable asset for research focused on identifying the genes responsible for drought adaptation in Arabidopsis and beyond.
Plant material & growth conditions
The A. thaliana accessions Kas-1 (CS903) and Tsu-1 (CS1640) were used as the original parent lines in developing the RIL population of 346 lines. Kas-1 and Tsu-1 were chosen as parents for developing this population as a result of their extreme differences in water use efficiency as measured by δ13C [39, 40]. RILs from this population served as the starting point for the NIL breeding program described below.
For the QTL experiment, seed of the RILs along with the parents were sown on soil (Fafard 4P mix, Conrad Fafard Inc., Agawam, MA) in 3-inch pots. Seeds were planted in a randomized complete block design consisting of 2 blocks, and then the pots were refrigerated at 4°C in darkness for 5 d to cold-stratify the seeds prior to commencement of a 8:16 h (light: dark) photoperiod in Conviron ATC60 growth chambers (Controlled Environments, Winnipeg, MB), at 23°C and 40% humidity during the day and 20°C and 50% humidity during the dark period. Light intensity was approximately 330 μmol m-2 s-1. Plants were grown for approximately 6 weeks prior to measurement. Stomatal conductance was measured in darkness on non-senescing leaves that were large enough to fully accommodate the leaf chamber (1 cm × 2 cm), using an infrared gas analyser (model Li-Cor 6400, LiCor Inc., Lincoln, NE). Prior to measurement the plants were dark adapted for 20 – 28 h. A humidifier was used to reduce variation in humidity over the course of the measurements. For each leaf 10 measurements were taken, with an interval of 10 s between measurements.
For the QTL validation experiment, plants were grown in a randomized complete block design consisting of 3 blocks where each genotype was replicated 6 times within each block. Plants were grown under exactly the same conditions as those described above except that the photoperiod was increased to 12;12 h (light:dark) to accommodate other experiments conducted in the same chamber. One major difference between the two experiments was the use of leaf porometers (model SC-1, Decagon Devices, Inc., Pullman, WA) rather than an infrared gas analyser for stomatal conductance estimates. Two non-senescing leaves were measured on each plant following the manufacturer’s recommended protocol.
Broad-sense heritability was estimated by calculating the ratio VG:VP, where VG is the among-RIL component of variance and VP is the total phenotypic variance. QTL mapping was performed in the R/qtl program of the R statistical package [65, 66], using Haley-Knott regression. Significance thresholds were determined using 1000 permutations. A penalized stepwise approach  was used for selection of a multiple-QTL model.
For the QTL validation experiment, data were analyzed with a linear mixed model using PROC MIXED in the SAS software package (SAS Institute Inc. 2003, Cary, NC) where block, row and column effects were treated as random.
Marker assisted NIL breeding program
To start, 7 RILs were selected from the original population of 346 using the code supplied in an additional file [see Additional file 1]. These 7 represented lines homozygous for Kas-1 alleles across one of each of the 5 chromosomes and all were crossed to Tsu-1 at least 10 times. Some attempted crosses may result in self-pollination due to technical error, thus we genotyped progeny to confirm they were F1s. In general, the real F1s were several times larger than the midparent value, so genotyping was almost unnecessary. Confirmed F1s were crossed back to Tsu-1 and each fruit was collected separately and considered a BC1 family, ultimately creating 25 families. 24 plants from each family were genotyped at the chromosome of interest and selected for selfing to generate BC1S1 seed. In addition to culling the occasional plant generated due to self-pollination, it was also necessary to remove individuals sired by (haploid) pollen from the F1 carrying Tsu-1 alleles for the chromosome of interest. In the next generation, 690 BC1S1 plants were genotyped with the 48 genome-wide SSRs described in . These were then ranked using an algorithm [see Additional file 2] to find lines that were largely Tsu-1, but carrying Kas-1 introgressions spanning the genome. In the end, 75 lines were selected which we screened at an additional 149 loci using the Sequenom MassARRAY® platform, of which 41 were polymorphic. 930 polymorphic loci were added to this marker data set via 2b-RAD  where class IIB restriction enzymes are used minimize genome complexity for a final total of 1011 genotyped.
DNA extraction and genotyping
Genomic DNA was isolated from lyophilized tissue collected from approximately 4-week-old, chamber grown plants using the DNeasy Plant Mini Kit (Qiagen, Valencia, CA) according to the manufacturer’s instructions.
The 48 polymorphic microsatellites used in this study were selected from the large number of those available in A. thaliana[68, 69], arabidopsis.org] due to easily distinguishable allele calls. Descriptions of the primers, PCR conditions and allele scoring are explained in .
DNA samples were used to prepare 2b-RAD libraries as previously described . A detailed protocol is available at the Meyer laboratory website (http://people.oregonstate.edu/~meyere/). Briefly, library preparation for 2b-RAD genotyping began with digestion of gDNA samples with AlfI (Fermentas) for 37°C for 3 h followed by ligation of adaptors at 4°C for 16 h. Ligation products were amplified by PCR and barcodes introduced to gel-extracted products in a second PCR reaction. Finally, libraries were pooled for multiplex sequencing on the SOLiD sequencing platform (Applied Biosystems). Raw sequences were processed to exclude low-quality reads, and the HQ reads that remained aligned in color-space using the SHRiMP software package  to AlfI sites extracted from the Arabidopsis genome (TAIR9). A custom Perl script was applied to eliminate short, statistically weak and ambiguous alignments (reads matching multiple sites equally well). Finally, genotypes were determined from nucleotide frequencies using custom Perl scripts to classify each locus as homozygous (minor allele frequencies [MAF] <1%), heterozygous (MAF > 25%), or undetermined (1% > MAF >25%). 20× coverage was required in the parental genomes to identify these alleles with high confidence, and a relaxed threshold of 10× in all other samples to maximize marker densities. Each polymorphic locus identified in these genotypes was compared with the parental genotypes (Tsu-1 and Kas-1) to assign it to one of these backgrounds, a comparison that would obviously not be possible for any loci genotyped in one parent but not the other as a result of variation in sequencing coverage. To reduce the effects of such missing data, we imported genotypes for Tsu-1 and Kas-1 from resequencing data (McKay, unpublished results) for any loci genotyped in one parent but not the other.
KASP SNP genotyping assays (LGC Genomics, Teddington, Middlesex, UK) were used for sub-NIL development. Primer sequences [see Additional file 8] were designed using sequence data from TAIR10  for amplification of SNPs identified and validated on the SNPlex genotyping system (Applied Biosystems) as described in . KASP is a novel allele-specific PCR assay that utilizes a FRET (Fluorescence Resonance Energy Transfer) system. In short, along with a common primer, allele-specific primers are designed to include a unique 18 bp sequence at the 5’ end. The unique sequences are identical to a pair of oligonucleotides with 3’ bound quenchers for a complement pair of 5’ fluorescently labelled oligos inside the reaction mix. During PCR, allele specific amplification leads to the generated product(s) outcompeting the quencher containing oligos for binding to the fluorescently labelled oligos, allowing for an observable signal to be measured using a light reader. The intensity of the signal(s) allows for a quantitative measure of SNP copy number.
Estimating chromosomal introgression length and number
The physical length of introgressions in the final NIL library was estimated using graphical genotypes . Physical length estimates of introgressions flanked by SSR markers were made using the location of the forward primers, SNP locations were determined by their location in the Col-0 reference genome. To avoid false-positives, an introgression was scored based on the presence of at least 3 consecutive markers with the Kas-1 genotype. Introgression boundaries were then defined by three consecutive markers with an alternative genotype. This helped avoid over-estimating introgression numbers due to occasional incorrect allele calls or differences in the location of loci in this population relative to the Col-0 genome used as a reference for mapping sequence reads. For the analysis of introgression discovery at varying marker densities an Excel Macro was written to sum the number of heterozygous and homozygous introgressions discovered. The loci included in replicated sampling were selected randomly using Excel’s RAND function.
Candidate gene identification
The full list of genes expected to lie within the QTL interval spanning physical positions 505,086 to 5,273,972 was assembled using TAIR10 . GO annotations for the full gene list were downloaded using the Bulk Data Retrieval and Analysis tool on TAIR10 and searched using the terms abscisic acid (ABA), stomata and water. Gene enrichment analysis was performed using the GO enrichment analysis tool in AmiGO .
This research was supported by NSF grants DEB-1022196 and DEB-0618302 to JKM and DEB-0618347 and IOS-09221457 to TEJ.
- Sax K: The association of size differences with seed-coat pattern and pigmentation in Phaseolus vulgaris. Genetics. 1923, 8: 552-560.PubMed CentralPubMedGoogle Scholar
- Soller M, Brody T: On the power of experimental designs for the detection of linkage between marker loci and quantitative loci in crosses between inbred lines. Theor Appl Genet. 1976, 47: 35-39.View ArticlePubMedGoogle Scholar
- Lander ES, Botstein D: Mapping Mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics. 1989, 121: 185-199.PubMed CentralPubMedGoogle Scholar
- Haley CS, Knott SA: A simple regression method for mapping quantitative trait loci in line crosses using flanking markers. Heredity. 1992, 69: 315-324. 10.1038/hdy.1992.131.View ArticlePubMedGoogle Scholar
- Jansen RC: Interval mapping of multiple quantitative trait loci. Genetics. 1993, 1: 205-211.Google Scholar
- Zeng ZB: Theoretical basis for separation of multiple linked gene effects in mapping quantitative trait loci. Proc Natl Acad Sci USA. 1993, 23: 10972-10976.View ArticleGoogle Scholar
- Zeng ZB: Precision mapping of quantitative trait loci. Genetics. 1994, 4: 1457-1468.Google Scholar
- Sen S, Churchill GA: A statistical framework for quantitative trait mapping. Genetics. 2001, 159: 371-387.PubMed CentralPubMedGoogle Scholar
- Price AH: Believe it or not, QTLs are accurate!. Trends Plant Sci. 2006, 11: 213-216. 10.1016/j.tplants.2006.03.006.View ArticlePubMedGoogle Scholar
- Rockman MV: Reverse engineering the genotype–phenotype map with natural genetic variation. Nature. 2008, 456: 738-744. 10.1038/nature07633.View ArticlePubMedGoogle Scholar
- Monforte AJ, Tanksley SD: Development of a set of near isogenic and backcross recombinant inbred lines containing most of the Lycopersicon hirsutum genome in a L. esculentum genetic background: a tool for gene mapping and gene discovery. Genome. 2000, 43: 803-813.View ArticlePubMedGoogle Scholar
- Landi P, Sanguineti MC, Salvi S, Giuliani S, Bellotti M, Maccaferri M, Conti S, Tuberosa R: Validation and characterization of a major QTL affecting leaf ABA concentration in maize. Mol Breeding. 2005, 15: 291-303. 10.1007/s11032-004-7604-7.View ArticleGoogle Scholar
- Iniguez-Luy FL, Lukens L, Farnham MW, Amasino RM, Osborn TC: Development of public immortal mapping populations, molecular markers and linkage maps for rapid cycling Brassica rapa and B. Oleracea. Theor Appl Genet. 2009, 120: 31-43. 10.1007/s00122-009-1157-4.View ArticlePubMedGoogle Scholar
- Eshed Y, Zamir D: A genomic library of Lycopersicon pennellii in L. esculentum: a tool for fine mapping of genes. Euphytica. 1994, 79: 175-179. 10.1007/BF00022516.View ArticleGoogle Scholar
- Frary A, Nesbitt TC, Grandillo S, van der Knaap E, Cong B, Liu JP, Meller J, Elber R, Alpert KB, Tanksley SD: fw2.2: a quantitative trait locus key to the evolution of tomato fruit size. Science. 2000, 289: 85-88. 10.1126/science.289.5476.85.View ArticlePubMedGoogle Scholar
- Takahashi Y, Shomura A, Sasaki T, Yano M: Hd6, a rice quantitative trait locus involved in photoperiod sensitivity, encodes the alpha subunit of protein kinase CK2. Proc Natl Acad Sci USA. 2001, 98: 7922-7927. 10.1073/pnas.111136798.PubMed CentralView ArticlePubMedGoogle Scholar
- Lin H, Liang ZW, Sasaki T, Yano M: Fine mapping and characterization of quantitative trait loci Hd4 and Hd5 controlling heading date in rice. Breeding Sci. 2003, 53: 51-59. 10.1270/jsbbs.53.51.View ArticleGoogle Scholar
- Juenger TE, Wayne T, Boles S, Symonds VV, McKay J, Coughlan SJ: Natural genetic variation in whole-genome expression in Arabidopsis thaliana: the impact of physiological QTL introgression. Mol Ecol. 2005, 15: 1351-1365.View ArticleGoogle Scholar
- Gerald JNF, Lehti-Shiu MD, Ingram PA, Deak KI, Biesiada T, Malamy JE: Identification of quantitative trait loci that regulate Arabidopsis root system size and plasticity. Genetics. 2006, 172: 485-498.View ArticleGoogle Scholar
- Xing YZ, Tang WJ, Xue WY, Xu CG, Zhang Q: Fine mapping of a major quantitative trait loci, qSSP7, controlling the number of spikelets per panicle as a single Mendelian factor in rice. Theor Appl Genet. 2008, 116: 789-796. 10.1007/s00122-008-0711-9.View ArticlePubMedGoogle Scholar
- Zhou L, Zeng Y, Zheng W, Tang B, Yang S, Zhang H, Li J, Li Z: Fine mapping a QTL qCTB7 for cold tolerance at the booting stage on rice chromosome 7 using a near-isogenic line. Theor Appl Genet. 2010, 121: 895-905. 10.1007/s00122-010-1358-x.View ArticlePubMedGoogle Scholar
- Tuinstra MR, Ejeta G, Goldsbrough PB: Heterogeneous inbred family (HIF) analysis: a method for developing near-isogenic lines that differ at quantitative trait loci. Theor Appl Genet. 1997, 95: 1005-1011. 10.1007/s001220050654.View ArticleGoogle Scholar
- Loudet O, Gaudon V, Trubuil A, Daniel-Vedele F: Quantitative trait loci controlling root growth and architecture in Arabidopsis thaliana confirmed by heterogeneous inbred family. Theor Appl Genet. 2005, 110: 742-753. 10.1007/s00122-004-1900-9.View ArticlePubMedGoogle Scholar
- Eshed Y, Zamir D: An introgression line population of Lycopersicon pennellii in the cultivated tomato enables the identification and fine mapping of yield-associated QTL. Genetics. 1995, 141: 1147-1162.PubMed CentralPubMedGoogle Scholar
- Monforte AJ, Tanksley SD: Fine mapping of a quantitative trait locus (QTL) from Lycopersicon hirsutum chromosome 1 affecting fruit characteristics and agronomic traits: breaking linkage among QTLs affecting different traits and dissection of heterosis for yield. Theor Appl Genet. 2000, 100: 471-479. 10.1007/s001220050061.View ArticleGoogle Scholar
- Koumproglou R, Wilkes TM, Townson P, Wang XY, Beynon J, Pooni HS, Newbury HJ, Kearsey MJ: STAIRS: a new genetic resource for functional genomic studies of Arabidopsis. Plant J. 2002, 31: 355-364. 10.1046/j.1365-313X.2002.01353.x.View ArticlePubMedGoogle Scholar
- Keurentjes JJB, Bentsink L, Alonso-Blanco C, Hanhart CJ, De Blankestijn Vries H, Effgen S, Vreugdenhil D, Koornneef M: Development of a near-isogenic line population of Arabidopsis thaliana and comparison of mapping power with a recombinant inbred line population. Genetics. 2007, 175: 891-905. 10.1534/genetics.106.066423.PubMed CentralView ArticlePubMedGoogle Scholar
- Alonso-Blanco C, Peeters AJM, Koornneef M, Lister C, Dean C, van den Bosch N, Pot J, Kuiper MTR: Development of an AFLP based linkage map of Ler, Col and Cvi Arabidopsis thaliana ecotypes and construction of a Ler/Cvi recombinant inbred line population. Plant J. 1998, 14: 259-271. 10.1046/j.1365-313X.1998.00115.x.View ArticlePubMedGoogle Scholar
- Juenger TE, Mckay JK, Hausmann N, Keurentjes JJB, Sen S, Stowe KA, Dawson TE, Simms EL, Richards JH: Identification and characterization of QTL underlying whole-plant physiology in Arabidopsis thaliana: δ13C, stomatal conductance and transpiration efficiency. Plant Cell Environ. 2005, 28: 697-708. 10.1111/j.1365-3040.2004.01313.x.View ArticleGoogle Scholar
- Esch E, Szymaniak JM, Yates H, Pawlowski WP, Buckler ES: Using crossover breakpoints in recombinant inbred lines to identify quantitative trait loci controlling the global recombination frequency. Genetics. 2007, 177: 1851-1858. 10.1534/genetics.107.080622.PubMed CentralView ArticlePubMedGoogle Scholar
- Laserna MP, Sanchez RA, Botto JF: Light-related Loci Controlling Seed Germination in Ler × Cvi and Bay-0 × Sha Recombinant Inbred-line Populations of Arabidopsis thaliana. Ann Bot-London. 2008, 102: 631-642. 10.1093/aob/mcn138.View ArticleGoogle Scholar
- Waters BM, Grusak MA: Quantitative trait locus mapping for seed mineral concentrations in two Arabidopsis thaliana recombinant inbred populations. New Phytol. 2008, 179: 1033-1047. 10.1111/j.1469-8137.2008.02544.x.View ArticlePubMedGoogle Scholar
- Li P, Wind JJ, Shi X, Zhang H, Hanson J, Smeekens SC, Teng S: Fructose sensitivity is suppressed in Arabidopsis by the transcription factor ANAC089 lacking the membrane-bound domain. Proc Natl Acad Sci USA. 2011, 108: 3436-3441. 10.1073/pnas.1018665108.PubMed CentralView ArticlePubMedGoogle Scholar
- Torjek O, Meyer RC, Zehnsdorf M, Teltow M, Strompen G, Witucka-Wall H, Blacha A, Altmann T: Construction and analysis of 2 reciprocal introgression line populations. J Hered. 2008, 99: 396-406. 10.1093/jhered/esn014.View ArticlePubMedGoogle Scholar
- Torjek O, Berger D, Meyer RC, Mussig C, Schmid KJ, Sorensen TR, Weisshaar B, Mitchell-Olds T, Altmann T: Establishment of a high-efficiency SNP-based framework marker set for Arabidopsis. Plant J. 2003, 36: 122-140. 10.1046/j.1365-313X.2003.01861.x.View ArticlePubMedGoogle Scholar
- Reif JC, Kusterer B, Piepho HP, Meyer RC, Altmann T, Schon CC, Melchinger AE: Unraveling epistasis with triple testcross progenies of near-isogenic lines. Genetics. 2009, 181: 247-257.PubMed CentralView ArticlePubMedGoogle Scholar
- Jan Lisec J, Steinfath M, Meyer RC, Selbig J, Melchinger AE, Willmitzer L, Altmann T: Identification of heterotic metabolite QTL in Arabidopsis thaliana RIL and IL populations. Plant J. 2009, 59: 777-788. 10.1111/j.1365-313X.2009.03910.x.View ArticlePubMedGoogle Scholar
- McKay JK, Richards JH, Nemali KS, Sen S, Mitchell-Olds T, Boles S, Stahl EA, Wayne T, Juenger TE: Genetics of drought adaptation in Arabidopsis thaliana II. QTL analysis of a new mapping population, Kas-1 × Tsu-1. Evolution. 2008, 62–12: 3014-3026.View ArticleGoogle Scholar
- McKay JK, Richards JH, Mitchell-Olds T: Genetics of drought adaptation in Arabidopsis thaliana: I. Pleiotropy contributes to genetic correlations among ecological traits. Mol Ecol. 2003, 12: 1137-1151. 10.1046/j.1365-294X.2003.01833.x.View ArticlePubMedGoogle Scholar
- Juenger TE, Sen S, Bray E, Stahl E, Wayne T, McKay J, Richards JH: Exploring genetic and expression differences between physiologically extreme ecotypes: comparative genomic hybridization and gene expression studies of Kas-1 and Tsu-1 accessions of Arabidopsis thaliana. Plant Cell Environ. 2010, 33: 1268-1284. 10.1111/j.1365-3040.2010.02146.x.View ArticlePubMedGoogle Scholar
- Wang S, Meyer E, McKay JK, Matz MV: 2b-RAD: a simple and flexible method for genome-wide genotyping. Nat Methods. 2012, 9: 808-810. 10.1038/nmeth.2023.View ArticlePubMedGoogle Scholar
- Davey JW, Hohenlohe PA, Etter PD, Boone JQ, Catchen JM, Blaxter ML: Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nat Rev Genet. 2011, 12: 499-510. 10.1038/nrg3012.View ArticlePubMedGoogle Scholar
- Clark AG, Hubisz MJ, Bustamante CD, Williamson SH, Nielsen R: Ascertainment bias in studies of human genome-wide polymorphism. Genome Res. 2005, 15: 1496-1502. 10.1101/gr.4107905.PubMed CentralView ArticlePubMedGoogle Scholar
- Young ND, Tanksley SD: Restriction fragment length polymorphism maps and the concept of graphical genotypes. Theor Appl Genet. 1989, 77: 95-101. 10.1007/BF00292322.View ArticlePubMedGoogle Scholar
- Caird MA, Richards JH, Donovan LA: Nighttime stomatal conductance and transpiration in C3 and C4 plants. Plant Physiol. 2007, 143: 4-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Christman MA, Richards JH, McKay JK, Stahl EA, Juenger TE, Donovan LA: Genetic variation in Arabidopsis thaliana for night-time leaf conductance. Plant Cell Environ. 2008, 31: 1170-1178. 10.1111/j.1365-3040.2008.01833.x.View ArticlePubMedGoogle Scholar
- Bauerle WL, Bowden JD: Separating foliar physiology from morphology reveals the relative roles of vertically structured transpiration factors within red maple crowns and limitations of larger scale models. J Exp Bot. 2011, 62: 4295-4307. 10.1093/jxb/err156.PubMed CentralView ArticlePubMedGoogle Scholar
- Yan WH, Wang P, Chen HX, Zhou HJ, Li QP, Wang CR, Ding ZH, Zhang YS, Yu SB, Xing YZ, Zhang QF: A major QTL, Ghd8, plays pleiotropic roles in regulating grain productivity, plant height, and heading date in rice. Mol Plant. 2011, 4: 319-330. 10.1093/mp/ssq070.View ArticlePubMedGoogle Scholar
- Wang S, Wu K, Yuan Q, Liu X, Liu Z, Lin X, Zeng R, Zhu H, Dong G, Qian Q, Zhang G, Fu X: Control of grain size, shape and quality by OsSPL16 in rice. Nat Genet. 2012, 44: 950-955. 10.1038/ng.2327.View ArticlePubMedGoogle Scholar
- Paterson AH, DeVerna JW, Lanini B, Tanksley SD: Fine mapping of quantitative trait loci using selected overlapping recombinant chromosomes, in an interspecies cross of tomato. Genetics. 1990, 124: 735-742.PubMed CentralPubMedGoogle Scholar
- Ducrocq S, Giauffret C, Madur D, Combes V, Dumas F, Jouanne S, Coubriche D, Jamin P, Moreau L, Charcosset A: Fine mapping and haplotype structure analysis of a major flowering time quantitative trait locus on maize chromosome 10. Genetics. 2009, 183: 1555-1563. 10.1534/genetics.109.106922.PubMed CentralView ArticlePubMedGoogle Scholar
- Acharya BR, Assmann SM: Hormone interactions in stomatal function. Plant Mol Biol. 2009, 69: 451-462. 10.1007/s11103-008-9427-0.View ArticlePubMedGoogle Scholar
- Nilson SE, Assmann SM: The control of transpiration. Insights from Arabidopsis. Plant Physiol. 2007, 143: 19-27.PubMed CentralView ArticlePubMedGoogle Scholar
- Schroeder JI, Allen GJ, Hugouvieux V, Kwak JM, Waner D: Guard cell signal transduction. Annu Rev Plant Physiol Plant Mol Biol. 2001, 52: 627-658. 10.1146/annurev.arplant.52.1.627.View ArticlePubMedGoogle Scholar
- Carbon S, Ireland A, Mungall CJ, Shu S, Marshall B, Lewis S: the AmiGO Hub, the Web Presence Working Group: AmiGO: online access to ontology and annotation data.http://bioinformatics.oxfordjournals.org/cgi/content/abstract/25/2/288,
- Dinka SJ, Campbell MA, Demers T, Raizada MN: Predicting the size of the progeny mapping population required to positionally clone a gene. Genetics. 2007, 176: 2035-2054. 10.1534/genetics.107.074377.PubMed CentralView ArticlePubMedGoogle Scholar
- Koornneef M, Meinke D: The development of Arabidopsis as a model plant. Plant J. 2010, 61: 909-921. 10.1111/j.1365-313X.2009.04086.x.View ArticlePubMedGoogle Scholar
- Cavell AC, Lydiate DJ, Parkin IAP, Dean C, Trick M: Collinearity between a 30-centimorgan segment of Arabidopsis thaliana chromosome 4 and duplicated regions within the Brassica napus genome. Genome. 1998, 41: 62-69.View ArticlePubMedGoogle Scholar
- Boivin K, Acarkan A, Mbulu R, Clarenz O, Schmidt R: The Arabidopsis genome sequence as a tool for genome sequence analysis in Brassicaceae. A comparison of the Arabidopsis and Capsella rubella genomes. Plant Physiol. 2004, 135: 735-744. 10.1104/pp.104.040030.PubMed CentralView ArticlePubMedGoogle Scholar
- Yogeeswaran K, Frary A, York TL, Amenta A, Lesser AH, Nasrallah JB, Tanksley SD, Nasrallah ME: Comparative genome analyses of Arabidopsis spp.: Inferring chromosomal rearrangement events in the evolutionary history of A. thaliana. Genome Res. 2005, 15: 505-515. 10.1101/gr.3436305.PubMed CentralView ArticlePubMedGoogle Scholar
- Lysak MA, Berr A, Pecinka A, Schmidt R, McBreen K, Schubert I: Mechanisms of chromosome number reduction in Arabidopsis thaliana and related Brassicaceae species. Proc Natl Acad Sci USA. 2006, 103: 5224-5229. 10.1073/pnas.0510791103.PubMed CentralView ArticlePubMedGoogle Scholar
- Schranz ME, Windsor AJ, Song B, Lawton-Rauh A, Mitchell-Olds T: Comparative Genetic Mapping in Boechera stricta, a Close Relative of Arabidopsis. Plant Physiol. 2007, 144: 286-298. 10.1104/pp.107.096685.PubMed CentralView ArticlePubMedGoogle Scholar
- Stacey G, VandenBosch K: "Translational" legume biology. Models to crops. Plant Physiol. 2005, 137: 1173-10.1104/pp.104.900141.PubMed CentralView ArticlePubMedGoogle Scholar
- Koebener RMD, Varshney RK: Development and application of genomic models for large-crop plant genomes. Model Plants and Crop Improvement. Edited by: Varshney RK, Robert MD. 2006, Koebner. Boca Raton: FL CRC Press, 1-9.Google Scholar
- Broman KW, Wu H, Sen S, Churchill GA: R/qtl: QTL mapping in experimental crosses. Bioinformatics. 2003, 19: 889-890. 10.1093/bioinformatics/btg112.View ArticlePubMedGoogle Scholar
- Broman KW, Sen S: A Guide to QTL Mapping with R/qtl. 2009, New York: SpringerView ArticleGoogle Scholar
- Manichaikul A, Moon JY, Sen S, Yandell BS, Broman KW: A model selection approach for the identification of quantitative trait loci in experimental crosses, allowing epistasis. Genetics. 2009, 181: 1077-1086. 10.1534/genetics.108.094565.PubMed CentralView ArticlePubMedGoogle Scholar
- Bell CJ, Ecker J: Assignment of 30 microsatellite loci to the linkage map of Arabidopsis. Genomics. 1994, 19: 137-144. 10.1006/geno.1994.1023.View ArticlePubMedGoogle Scholar
- Symonds VV, Lloyd AM: An analysis of microsatellite loci in Arabidopsis thaliana: mutational dynamics and application. Genetics. 2003, 165: 1475-1488.PubMed CentralPubMedGoogle Scholar
- Rumble SM, Lacroute P, Dalca AV, Fiume M, Sidow A, Brudno M: SHRiMP: Accurate mapping of short color-space reads. PLoS Comput Biol. 2009, 5: e1000386-10.1371/journal.pcbi.1000386.PubMed CentralView ArticlePubMedGoogle Scholar
- Lamesch P, Berardini TZ, Li D, Swarbreck D, Wilks C, Sasidharan R, Muller R, Dreher K, Alexander DL, Garcia-Hernandez M, Karthikeyan AS, Lee CH, Nelson WD, Ploetz L, Singh S, Wensel A, Huala E: The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools. Nucleic Acids Res. 2012, 40: D1201-D1210.View ArticleGoogle Scholar