 Methodology article
 Open access
 Published:
Targeted Recombinant Progeny: a design for ultrahigh resolution mapping of Quantitative Trait Loci in crosses between inbred or pure lines
BMC Genetics volume 16, Article number: 76 (2015)
Abstract
Background
Highresolution mapping of the loci (QTN) responsible for genetic variation in quantitative traits is essential for positional cloning of candidate genes, and for effective marker assisted selection. The confidence interval (QTL) flanking the point estimate of QTNlocation is proportional to the number of individuals in the mapping population carrying chromosomes recombinant in the given interval. Consequently, many designs for high resolution QTN mapping are based on increasing the proportion of recombinants in the mapping population.
The "Targeted Recombinant Progeny" (TRP) design is a new design for high resolution mapping of a target QTN in crosses between pure, or inbred lines. It is a threegeneration procedure generating a large number of recombinant individuals within a QTL previously shown to contain a QTN. This is achieved by having individuals that carry chromosomes recombinant across the target QTL interval as parents of a large mapping population; most of whom will therefore carry recombinant chromosomes targeted to the given QTL. The TRP design is particularly useful for high resolution mapping of QTN that differentiate inbred or pure lines, and hence are not amenable to high resolution mapping by genomewide association tests.
Results
In the absence of residual polygenic variation, population sizes required for achieving given mapping resolution by the TRPF_{2} design relative to a standard F_{2} design ranged from 0.289 for a QTN with standardized allele substitution effect = 0.2, mapped to an initial QTL of 0.2 Morgan to 0.041 for equivalent QTN mapped to an initial QTL of 0.02 M. In the presence of residual polygenic variation, the relative effectiveness of the TRP design ranges from 1.068 to 0.151 for the same initial QTL intervals and QTN effect. Thus even in the presence of polygenic variation, the TRP can still provide major savings. Simulation showed that mapping by TRP should be based on 3050 markers spanning the initial interval; and on at least 50 or more G_{2} families representing this number of recombination points,.
Conclusions
The TRP design can be an effective procedure for achieving high and ultrahigh mapping resolution of a target QTN previously mapped to a known confidence interval (QTL).
Background
Genetic analysis of complex quantitative traits involves mapping the polymorphic sites (QTN – Quantitative Trait Nucleotide) responsible for genetic variation in these traits to their individual chromosomal regions (QTL – Quantitative Trait Loci). The QTL is defined here as a confidence interval (CI) within which the QTN is located with confidence (1α), where α is the probability that the QTN is found outside the CI. High resolution QTL mapping (i.e., localizing the QTL to a narrow confidence interval) is essential for positional cloning of candidate genes, and for effective marker assisted selection or marker assisted introgression. With the advent of high density SNP microarrays, high resolution mapping of segregating QTN within populations of outcrossing species (all animals, and many plants), is now based on markerQTN association due to populationwide linkage disequilibrium (LD), socalled “Whole Genome Association Studies” (WGAS) [1]. WGAS exploits for mapping the very large number of recombination events accumulated across an outcrossing population over the untold generations of sexual reproduction since the origin of the QTN. These recombination events limit populationwide LD between markers and QTN in the population to a very narrow region about the QTN, so that statistically significant association between marker and QTN indicates that the QTN is very close to the marker. Within pure lines of selfing species, or inbred lines of outcrossing species, however, markers and QTN are at fixation. Hence, WGAS is clearly not possible and other designs have been developed. Many of these were already reviewed in the classic Darvasi paper [2]. Basically they fall into two types: Group I, designs that increase mapping resolution across the entire genome; Group II, designs that are targeted to a specific QTL. Group I designs include the Advanced Intercross Line (AIL) design [3,4] which can be applied ad hoc to any F_{2} population. These designs also include specialized mouse stocks based on the AIL principle, developed specifically for high resolution mapping, the Heterogeneous Stock (HS) [5] and the Diversity Outcross [6]. In a happy stroke of good fortune, the commercially available MF1 outbred stock, was adventitiously found to be highly suitable for Group I mapping [7,8]. All three stocks have the advantage that mapping is on the basis of founder haplotypes instead of single markers, in this way avoiding dilution of effects when the same marker allele is associated with both alternative QTN alleles. On the other hand they all have the limitation that they can only access QTN that are segregating among the founder lines of the crosses. Resources comparable to the HS, Diversity outcross and MF1 are not available for plant species. In the Near Isogenic Line (NIL) design, analyzed in detail in [9], the entire genome of a target donor line is dissected into smaller segments, each isolated in a different NIL through successive backcrossing to a recipient line followed by a number of selfing generations. Each segment can then be tested against the corresponding donor segment for QTL detection. Construction of NIL libraries was pioneered by the tomato geneticists [10,11] and has since been applied with success in a wide assortment of plant species (listed in [9]). Also in Group I are new major resources consisting of sets of many Recombinant Inbred Lines (RIL) derived from intercrosses among a large set of founder parent lines chosen to maximize genetic diversity within the resource. These include the Collaborative Cross (CC) Mouse Resource population [12,13]; and the conceptually identical Multiparent Advanced Generation InterCross (MAGIC) resources in plants [14]. The MAGIC resource was first proposed by Cavanough [14], based on the discussions leading to the CC resource, and has already been applied to develop MAGIC populations in Arabidopsis [15], wheat [16] and rice [17]. The advantages and strengths of the CC and MAGIC resources are too many to detail here. When a large mapping population can be generated and genotyped but the limiting factor is phenotyping (as in microarray experiments), selective phenotyping (reviewed in [18]) can be employed to select a subset of individuals to maximize their mapping power, e.g., by maximizing their genotypic dissimilarity [19] or the complementariness of crossover sites within the sample [20].
Group II designs are aimed at high resolution mapping targeted to a specific QTL previously mapped with high power to a relatively large CI by standard F_{2}, BC, or AIL QTLmapping designs. Group II designs were reviewed by Darvasi [2], and new designs do not seem to have been proposed since then. All of these designs are based on the principle of “Chromosome Dissection” pioneered by the Drosophila geneticists over half a century ago (e.g., [21]). In Recombinant Progeny Testing (RPT) [2, 22] a series of individuals in the original F_{2} or BC mapping population that carry recombinant chromosomes with the points of recombination laddered across the target region are queried individually as to whether the QTN location is upstream or downstream of the recombination point. In the IntervalSpecific Congenic Strains (ISCS) design [23] the same class of individuals as for RPT are individually backcrossed repeatedly to one of the parental strains (chosen to be recessive for the QTL effect, if applicable) to retain only a segment containing the recombination point and small flanking donor intervals. This establishes a series of congenic strains that cover the target region. These are then queried individually to determine whether the QTN is located upstream or downstream of the recombination point. A great advantage of the ICSC design is the reduction in residual genetic variance by the repeated backcrossing to the inbred parent. This enables strain status with respect to the QTL to be determined with relatively few individuals. This design has been widely applied in mice (e.g., [24]). When large numbers of individuals are available for the mapping population, selective recombinant genotyping [25] can be employed to reduce genotyping costs.
In the present study, we propose a new, "Targeted Recombinant Progeny" (TRP) design for high and ultrahigh resolution mapping of specific target QTN that are at fixation within populations, but differ across populations, e.g., in crosses between pure lines of selfer species, or between highly inbred lines. The TRP design is a threegeneration procedure for generating a large number of recombinant individuals within a QTL shown by previous mapping to contain a QTN. Similar to RPT and ISCS this is achieved by identifying “founder” individuals in an F_{2}, BC, or AIL population that carry chromosomes that are recombinant in the target QTL. In contrast to RPT, there is no attempt to determine QTN status of individual founders. Rather, by having these founder individuals serve as parents of a large progeny population, a mapping population is generated most of whose members carry recombinant chromosomes targeted to the given QTL. Since the TRP mapping population consists primarily of recombinant individuals, this provides for high resolution mapping of the QTN within its QTL, with appreciably smaller total mapping populations than required for equal precision by classical F_{2}, BC or even AIL designs. The only condition is that the number of founder individuals is sufficient to provide enough points of recombination to refine the QTN location to the desired degree. It should be stressed, that the TRP design is aimed at high resolution mapping of a specific QTN previously mapped to a QTL. This contrasts to the classical F_{2}, BC and AIL designs for which a single set of progeny provides the required recombinants for mapping all QTN that are segregating in the mapping population.
Methods
(i) Notation and assumptions
Notation
To simplify derivations, map distance are given in units of Morgans (M), rather than centiMorgans (cM). Thus, 20 cM = 0.2 M.
Assumptions
We assume a QTN that has been previously mapped to a point location within a QTL of width C M, with confidence level (1α). The purpose of the TRP is to generate a new mapping population which is densely populated by recombinant chromosomes targeted to this QTL, so that a new point location with (1α) CI of width c = C/s is obtained, where s is the factor by which the original CI is reduced. The TRP achieves this by generating an F_{2}, BC or AIL population and identifying within it individuals that carry a recombinant chromosome in the target interval. These are then selfed or backcrossed to one of the parental lines to generate a TRP mapping population highly enriched with recombinant chromosomes in the interval of interest.
It is convenient to describe the TRP by starting with an F_{1} generation created by crossing two pure lines homozygous for alternative alleles at a large number of marker loci and at the target QTL. However, the TRP can equally be initiated from any population (denoted the G_{1} generation), with a similar F_{1}type genetic structure. On the basis of previous mapping studies, these G_{1} individuals will have genotypes of known haplotype composition, including phase of the QTN allele. By appropriate notation of marker and QTN alleles, these genotypes can be given the form: M_{U}QM_{D}/m_{U}qm_{D,} where M_{U} and M_{D} are markers defining the upstream (M_{U}) and downstream (M_{D}) (1α) CI boundaries of the original QTL of width C M (see Figure 1). Q and q are alternative alleles at the QTN, and a large number of markers with alternative alleles (not shown) are found spanning the region between M_{U} and M_{D}.
We now describe the construction of the TRP mapping population in detail. Based on Weller and Soller [26], we then derive the required size of the TRP mapping population to achieve target (1α) CI about the QTN point location and compare to that required by standard F_{2}, BC or AIL mapping populations for equivalent precision.
(ii) The G_{2} population and the G_{ 2R } individuals
Depending on the reproductive biology of the species, the G_{1} parent individuals are selfed (TRPF_{2} design) or backcrossed to one of the parental lines (TRPBC design) to produce a G_{2} progeny population some of whom carry a recombinant chromosome in the target QTL. In what follows, we continue with the analysis of the TRPF_{2} design. The TRPBC design is a simpler application of the same principles and will be briefly described in Appendix I.
The G_{2} progeny are genotyped for the markers M_{U} and M_{D}, and G_{2R} individuals, each carrying a single recombinant chromosome in the interval M_{U} to M_{D}, are identified (see Figure 1). In principle, this would also uncover triple recombinants but these will be in vanishingly small frequency and can safely be ignored. Since there are two types of recombinant haplotypes (M_{U}m_{D}, and m_{U}M_{D}) and two types of parental haplotypes (M_{U}M_{D}, and m_{U}m_{D}) there will be four types of G_{2R} individuals, as shown in Table 1. In addition, a small proportion (= C^{2}) of G_{2R} individuals carrying two independent recombinant chromosomes will also be produced. In practice, these individuals would not be included in the G_{2R} group, as the families they produce do not contain nonrecombinant progeny needed to correct for polygenic family effects (see section xv for further details). To make up for this, it will be sufficient to increase the size of the G_{2} population by the proportion C^{2}, e.g., by 4% if the original CI is 0.20 M. This will be considered as negligible and will not be taken into account in what follows.
The G_{2R} individuals are genotyped for the markers spanning the region between M_{U} and M_{D}, identifying the haplotype of the recombinant chromosome and the point of recombination. This information will allow some selection among the G_{2R} individuals chosen to produce the G_{3} generation to obtain a more evenly spaced distribution of recombination points across the original QTL. It also enables the complete genotypes of the G_{3} generation to be inferred from their G_{2R} parents, without requiring further genotyping of the G_{3} individuals (this is explained in the next section).
(iii) Producing the G_{ 3R } mapping population and identifying a new point location for the QTN
Continuing with the TRPF_{2} design, each parent G_{2R} individual will produce on selfing three classes of G_{3} progeny: Class I, Double recombinants in proportion 0.25, carrying two recombinant haplotypes; Class II, Single recombinants in proportion 0.50, carrying one recombinant and one nonrecombinant (parental) haplotype; and Class III, Double nonrecombinants in proportion 0.25, carrying two nonrecombinant parental haplotypes (Table 1). The three G_{3} genotype classes within each family are identified by genotyping the G_{3} population for the pair of flanking markers (M_{U} and M_{D}) that define the target region. The combined Class I and Class II recombinant progeny of all G_{2R} individuals form the G_{3R} mapping population. Thus, in contrast to the usual F_{2}, BC or even AIL designs, in which the mapping population consists of only a small proportion of informative recombinant progeny for any target QTL, the TRP mapping population consists primarily of informative recombinant progeny, targeted to a specific QTL. The nonrecombinant Class III G_{3} progeny group is not included in the G_{3R} mapping population. However, as will be shown later (section xv), it will contribute to evaluation and correction of family polygenic values.
Once G_{3R} recombinant status is determined, further genotyping is not required to obtain full G_{3R} genotypes, since (neglecting rare new recombination), the full genotype of any G_{3R} individual across the small target c interval, is determined by the genotype of the G_{2R} parent of that individual.
At some point along the way, either before or after genotyping, the G_{3} or G_{3R} population is phenotyped. A suitable ttest (or other appropriate test) is implemented in the G_{3R} population for each marker in the interval M_{U} to M_{D}, with data pooled across all G_{3R} families. The marker with the most significant ttest is identified as the new point location of the QTN, denoted M_{1}. Our goal is to determine the total G_{3} population size, N_{TG3}, required to have the (1 − α) CI about M_{1} equal to c = C/s. We do this in two steps. basing our approach on Weller and Soller [26]. As shown in that paper, for any given mapping population, the (1α) confidence level of a symmetrical interval of width c M, about the point location of a QTN is solely a function of: (i) α, (ii) the number, K, of recombinant chromosomes across that interval in the mapping population, and (iii) the standardized allele substitution effect, d, at the QTN. It follows, that for given α and d, K is a constant. Consequently, the required size, N_{G3}, of a mapping population to deliver a (1α) CI of width c will be the size of the population that will deliver K recombinant chromosomes across an interval of size c. In the first step, then, we calculate K for the given d and target α. N_{G3} will then depend solely on the proportion, R, of recombinant chromosomes in the mapping population such that, N_{G3} = K/R. In the second step, we note that since the new point location of the QTN in the original QTL is not known, the entire original (1α)CI = C M, is considered as composed of s consecutive smaller subCI each of width c M. Each of these s subCI will require generating a G_{3} population of size N_{G3}. Thus, N_{TG3}, the total number of G_{3} progeny required to map the QTL to the subCI of width c, within its original QTL will be:
We now derive expressions for N_{G3}, as a function of α and d for the TRPF_{2}design.
(iv) Required size of the G_{3} mapping population for given α and d
At this point it is convenient to shift our attention from Table 1, which gives the overall composition of the G_{3} population; to Table 2 which focuses our attention on the new point location of the QTN at M_{1}, and its boundary marker, M_{2}. With respect to these markers, the G_{3} population derived from the G_{2R} parents has the composition shown in Table 2. All told there are four G_{2R} family types and eight G_{3R} recombinant progeny groups, denoted A to H. For each G_{3R} recombinant group, Table 2 also shows the genotypic value, and the expected frequency in that fraction of the G_{3} population defining the new 95% CI = c about M_{1}, calculated on the assumption that each of the four G_{2R} parental types contributes equally to the G_{3} population.
Given the G_{3R} mapping population, with single markermapping the new point location of the QTN will be at the marker, M_{1}, with the greatest difference between alternative marker genotypes, with marker M_{2} defining the oneside limit of the target (1− \( \boldsymbol{\upalpha} \) )CI. As shown by Weller and Soller [26] the probability that the CI of QTL location includes the marker M_{2} located at a remove of L M from M_{1}, is equal to the probability of obtaining the value
where,
Z_{α/2} is the standard normal variable corresponding to a probability of α/2, D = E(M_{1})  E(M_{2}), where E(M_{1}) is the expected effect at M_{1} (QTN located at the marker), E(M_{2}) is the expected QTN effect at M_{2} (located at a remove from the QTN) considering recombinant individuals only; and SE(D) is the standard error of D. Some thought will show that only recombinant individuals in the region M_{1} to M_{2} can contribute to a difference in expectation for M_{1} and M_{2}, as nonrecombinants have the same value at both markers.
Letting italics denote the mean genotypic value of the corresponding marker genotype group (including recombinant genotypes only), we have
From Table 2 it is apparent that each of the four marker genotype groups is composed of three recombinant marker groups. For example the marker genotype group M_{2}M_{2} is composed of recombinant genotype groups C, D, and G of Table 2 with respective genotypes: m_{1}M_{2}/m_{1}M_{2}, m_{1}M_{2}/m_{1}M_{2}, and m_{1}M_{2}/M_{1}M_{2}; having genotypic values d, d, and h; frequencies 1/16, 1/16, and 2/16 of the G_{3} population; and relative frequencies 1/4, 1/4, and 1/2 within the M_{2}M_{2} genotype. The mean genotypic value of the M_{2}M_{2} genotype group, including recombinants only (M _{ 2 } M _{ 2 } ) is the mean of the genotype groups C, D, and G, weighted by their relative frequencies in the M_{2}M_{2} recombinant group, i.e. 1/4, 1/4 and 1/2, respectively.
On this basis, letting A, B, C, D, E, F, G, H, represent the genotypic values of the corresponding genotype groups, we have
Substituting in (1) and combining like terms, we obtain
Substituting genotypic values of the recombinant groups we have
To calculate SE(D), we assume that variance within the G_{3} families is the same as variance within the F_{2} generation (set to 1.0 for standardization). This is conservative, as the variance within G_{3} families will generally be less than in an F_{2} population, depending on degree of homzoygosity in their G_{2R} parent. On this basis, we can write:
σ ^{2} _{A} = σ ^{2} _{B} = σ ^{2} _{C} = σ ^{2} _{D} = 1/(N _{ G3}/16) = 16/N _{ G3} (as each of these genotype groups comprises 1/16 of the G_{3} mapping population of size N_{G3}
σ ^{2} _{E} = σ ^{2} _{F} = σ ^{2} _{G} = σ ^{2} _{H} = 1/(N _{ G3}/8) = 8/N _{ G3} (as each of these genotype groups comprises 1/8 of the entire G_{3} mapping population).
Then,
Substituting, we have
Solving for N_{G3}, we have
For example, if we set α = 0.05 and d = 0.2, we have Z_{α/2} = 1.96, and N_{G3} = 256.4
Note, however, that the interval between M_{1} and M_{2} defines only half of the CI of QTL location. Hence, to cover the entire CI would require twice this
On these assumed values for α and d, the total G_{3R} population required for CI(0.95) = N_{G3} = 512.85 This number is constant for given d, and α, and does not depend on the target size of the confidence interval e.g., for α = 0.05 and d = 0.25, N_{G3} = 328.22. If N_{G3} is achieved across a small interval, then the (1α) CI will be narrow; and if N_{G3} is achieved across a large interval, than the (1α) CI will be wide.
Note, that in contrast to the usual F_{2} or similar mapping populations where the difference between alternative marker genotype groups is a function of proportion of recombination between marker and QTL, this is not the case for the TRP design, for which, since it is based on recombinant progeny only, the difference between alternative marker genotype groups is independent of proportion of recombination.
(v) Required total G_{3} population size, N_{TG3}
Since the location of the QTN in the original QTL is not known, the entire original CI_{(1α)} = C M, is considered as composed of s consecutive smaller subCI_{(1α)}, each of width c M (the target CI_{(1α)}). Each of these s subintervals will require generating a G_{3} population of size N_{G3}. Thus, for the example above (α = 0.05, d = 0.2) the total number of G_{3} progeny required to map the QTL to the subCI of size c = C/s within its original QTLR will be :
As noted, this number is constant for given s, d, and α, and does not depend on the target size c. Continuing with the same example (α = 0.05, d = 0.20), if we are reducing an original CI of C = 0.20 M to a target CI of 0.10 M (s = 2), then N_{G3} = 512.85 , and N_{TG3} = 2*512.85 = 1025.70 Of this total, 75% are included in the actual G_{3R} mapping population, and 25% are produced but are nonrecombinant and not included in the mapping population.
(vi) Required number of G _{ 2R } parent individuals
The number of G_{2R} individuals needed to produce the required number of G_{3} progeny for each subCI of size c, is:
where
n is the number of progeny that are produced by selfing each G_{2R} parent.
Continuing our example, and assuming n = 50, N_{G2R} = 512.85/50 = 10.26. It is important that N_{G2R} is large enough to provide a sufficient density of recombination points across the target QTL to define the subCI boundaries with a degree of precision (see simulation for further details).
(vii) Required G _{ 2 } population size
By definition, when c is small the proportion of G_{1} gametes carrying a recombinant chromosome across the interval c M, will be c. In the G_{2} generation produced from these gametes, a proportion C^{2} will carry two recombinant chromosomes and are excluded from the G_{3R} mapping population as explained above. As C^{2} is generally small it is conveniently neglected in the calculations. In the remaining G_{3R} population, all recombinants across the region c, will be present in proportion 2c(1c), as heterozygotes carrying a recombinant chromosome together with a parental haplotype. Thus, as a close approximation, a bit on the high side, we have
\( {\mathsf{N}}_{\mathsf{G2}\mathsf{R}}=\mathsf{2}\mathsf{c}\left(\mathsf{1}\hbox{} \mathsf{c}\right){\mathsf{N}}_{\mathsf{G2}} \) and solving for N_{G2} we have
Noting that N_{G2R} = N_{G3}/n and N_{G2} = N_{G2R}/2c(1c), we can solve for N_{G2} as a function of N_{G3}
This is convenient as N_{G3} is the basic number determining all the remaining numbers.
The same total of G_{2} individuals needed to yield N_{G2R} recombinant individuals for one subCI of size c M, simultaneously produces N_{G2R} recombinant individuals for each subCI of size c M across the entire genome. Hence, the same G_{2} population of size N_{G2} can provide N_{G2R} individuals for all subCI of the designated QTL (and indeed, as noted previously, for any other QTL mapped in the original F_{2} population).
Thus, for our example, C = 0.2, c = 0.1, n = 50, N_{G3} = 512.75 , so that
That is, a G_{2} population of size 56.98 individuals (113.97 chromosomes), will contain 10.26 G_{2R} recombinants for each target interval of 0.1 M in an initial CI of size 0.2M.
(viii) Required G _{ 1 } population size
The number of G_{1} individuals required to produce the G_{2} population will obviously equal
For our example, N_{G1} = 512.85/450 = 1.14, so that 2 or 3 G_{1} individuals should suffice.
(ix) Total progeny required across the G_{1} , G_{ 2 } and G_{ 3 } generations (N_{TRP})
Combining expressions for N_{G1}, N_{G2} and N_{TG3}, as function of N_{G3}, we have
For our example, N_{G3} = 512.85, n = 50, c = 0.1, s = 2
The contribution of N_{G1} is always negligible, and can safely be ignored. For the present example with relatively large c, the contribution of N_{G2} is small, but when c is small, this is not the case, and N_{G2} makes an appreciable contribution to N_{TRP}.
To reduce the same 20 cM original QTL to 2 cM (s = 10), would require 5,196.71 individuals. All the increase will come from the increased G_{3} population; G_{1} and G_{2} remain the same.
(x) Total required genotyping and phenotyping
All individuals of the G_{2} generation are genotyped for the pair of markers flanking the original QTL to identify the G_{2R} individuals. These, making up 2c(1c) of the G_{2} generation, are genotyped for k internal markers to locate the QTN more precisely within the original QTL. Based on the simulation results, about k = 30 to 50 internal markers should be genotyped in the G_{2R} individuals to extract maximum mapping information from a given G_{3} population. Finally, the entire G_{3} population is genotyped for a single marker, to identify their genotype (Class I, Class II or Class III in Tables 1 and 2) with respect to the haplotype transmitted from their G_{2R} parent. Thus, total genotyping data points (g.d.p.) will be
For example, from the combination C = 0.02, s = 4, c = 0.005 we have N_{TG3} = 7,591, N_{G2} = 3,814. Then, taking k = 50, we have Total g.d.p. = 15,131: a remarkably low figure considering today’s genotyping capacities and costs.
With respect to phenotyping, the entire G_{3} population is phenotyped, but there is no need to phenotype the G_{2} generation. If phenotyping costs are a major consideration, some form of selective phenotyping [18,19] may be useful to reduce phenotyping data points at the expense of an increase in total population size and a negligible increase in genotyping costs
(xi) Total progeny required for an F_{2} mapping population to achieve equivalent map resolution:
Following Weller and Soller [26], population size required by an F_{2} population to achieve 95% CI = c M is:
For our example, d = 0.2, c = 0.1 we have
Noting that in our example N_{G3} = 5.34(Z_{α/2} ^{2})/d^{2}, the ratio of N_{TRP} to N_{F2} is given by
Examination of the final expression for N_{TRP}/N_{F2} shows that it is a function of n, c, and s only, and is not affected by α, nor, even more remarkably, by d.
(xii) Genotyping and phenotyping the F_{2}
All N_{F2} individuals of the F_{2} mapping population are genotyped for the two markers flanking the original QTL of width C. This will identify 2C(1C) F_{2R} individuals carrying a recombinant chromosome in the target region. These are then genotyped for all k internal markers. Thus, Total g.d.p. = 2N_{F2} + 2C(1C)kN_{F2}. For the example we used for the TRP design (C = 0.02, s = 4, c = 0.005), we have N_{F2} = 75,000, giving Total g.d.p. = 297,000; 20fold more than required by the TRP design.
Ordinarily, all N_{F2} individuals would be phenotyped. However, if genotyping precedes phenotyping there is opportunity to greatly reduce phenotyping numbers by phenotyping the recombinant F_{2} individuals only. Similarly, if phenotyping precedes genotyping it may be possible to reduce genotyping numbers by use of selective genotyping [27] or even selective DNA pooling [28,29].
(xiii) Mapping requirements in other TRP designs
Mapping requirements in an TRPAIL are the same as in an TRPF_{2}, except that the mapping population sizes are less by a factor of 0.5t, where t is the generation number of the AIL. Mapping population size requirements in a TRPBC depend strongly on degree of dominance at the QTL (h), and range from twice to half that of an TRPF_{2}, depending on whether h = 0 or h = 1, respectively. For the case h = 0, we expect the TRPBC to require about twice as many data points as a TRPF2 for equivalent power, and this is indeed what is found (Appendix I).
(xiv) Polygenic variance component of SE(D)
To this point we have calculated SE(D) on the assumption that the expected value of the G_{3} marker genotype group within families is μ + d , μd or μ + h, depending on the genotype at the QTL. The underlying assumption being that μ is the same for all families and marker groups within families. This is true for a standard F_{2} population, where all F_{1} parent individuals have the same polygenic value.
However, in the TRPF_{2} design, each G_{3} family is generated by selfing from a different G_{2R} individual. These will differ in polygenic value with standardized polygenic variance between the G_{3} families equal to the heritability (h^{2}) of the trait in the G_{2} generation plus a small dominance component [30]. This has a strong effect on the SE(D). In the absence of a polygenic effect, and using the previous notation, we have
Four families types are represented (Table 2). Groups A and E are derived from Family type 1; Groups B and F from Family type 2; groups C and G from Family type 3 and groups D and H from Family type 4. Then, in the presence of a polygenic effect, to each of the groups, in addition to the expected genotypic value, we need to add a polygenic value taken from a normal distribution with mean 0 and variance h^{2}.
Letting P _{ 1 } , P _{ 2 } , P _{ 3 } , P _{ 4 } be the polygenic effects of Families 1, 2, 3, 4 respectively, we have:
As before
σ ^{2}(2P _{1}) = σ ^{2}(2P _{2}) = σ ^{2}(2P _{3}) = σ ^{2}(2P _{4}) = 4h ^{2}/k, where k is the average number of replicate families of each type.
Then
The polygenic component has a powerful effect. For example, if N_{G3} = 1000, n = 50, then k = (1000/50)/4 = 5, Then if h^{2} = 0.25, we have SE^{2}(D) = 24/1000 = 0.024 in absence of a polygenic effect; while the polygenic effect will add 4*0.25/5 = 0.20, to give SE^{2}(D’) = 0.224. For SE^{2}(D’) to equal 0.024 in the presence of a polygenic effect, we would need to increase N_{G3} tenfold, from 1000 to 10,000, giving k = 50, and SE^{2}(D’) = 24/10,000 + 4h^{2}/50 = 0.0224.
(xv) Using nonrecombinant progeny to correct for polygenic effects
To deal with polygenic effects in a more effective manner, we propose to use the nonrecombinant group of each family as an estimate of the family polygenic effect, and express the recombinant groups as deviations from the nonrecombinant group. Examination of Table 2 shows that there are two classes of nonrecombinant groups. For Families 1 and 3, the nonrecombinant groups have genotype: M_{1}QM_{2}/M_{1}QM_{2}, with genotype values + d; for Families 2 and 4, the nonrecombinant groups have genotype m_{1}qm_{2}/ m_{1}qm_{2} , with genotype values –d. Thus, to bring all nonrecombinant families to the same expectation, we need to estimate d from the previous mapping experiment that defined the original QTL, and correct Families 2 and 4 for the effect of the QTN by adding 2d to the nonrecombinant groups of these families. When this is done, the variation among nonrecombinant groups of the different families, will be due to polygenic variation alone. Thus, the mean of the nonrecombinant group will represent the polygenic effect of the family and will be common to all genotype groups within the family. Consequently, the expected deviation of the recombinant groups from the nonrecombinant group mean will be due to genotype at the target QTL only, and not due to polygenic effects.
Although the expected polygenic effect of the nonrecombinant group and of the recombinant groups is the same, the mean of a nonrecombinant group will have a sampling variance, depending on the number of individuals in the group. Since proportion of total family in the nonrecombinant groups is 1/16, this will equal 16/N_{G3}. Note, in this case, k is not relevant as it does not make any difference how many subfamilies are within a given family, sampling variance of the nonrecombinant groups will depend on the total population size only, and not on how it is divided among replicate families.
If we let Q _{ 1 } , Q _{ 2 } , Q _{ 3 } , Q _{ 4 }, be the sampling deviation of the nonrecombinant group from expectation, we have
Where, as before
Then SE^{2}(D) = (1/4)(4*16/N_{G3} + 4*8/N_{G3} + 4*4*16/N_{G3}) = (16/N_{G3} + 8/N_{G3} + 64/N_{G3}) = 88/N_{G3}
Substituting and solving for N we have
And solving for N_{G3}, we have
N_{G3} = 88Z_{α/2} ^{2}/9d^{2} = 9.88Z_{α/2} ^{2}/d^{2} for one side of the CI, and twice this, 19.76Z_{α/2} ^{2}/d^{2} for twosided CI.
If we set α = 0.05, d = 0.2 as in our example, we have Z_{α/2} = 1.96, and total G_{3} population required for a twoside CI, N_{G3} = 1897.75. This can be compared to the corresponding value N_{G3} = 512.85 in the absence of polygenic effects. Since all other numbers in the analysis are functions of N_{G3}, the presence of polygenic effects of magnitude h^{2} = 0.25, increases required populations sizes 3.7fold all down the line. Relative to a standard F_{2} design, TRPF_{2} population sizes for equivalent power as a proportion of the required F_{2} population would increase from 0.288 for c = 0.1, s = 2 in the absence of polygenic effects, to 1.068 in their presence, i.e., no savings. For c = 0.01, s = 2, however, relative savings would be 0.152, which is still appreciable.
Simulation
The above calculations assume a saturated marker map, and saturation of the QTL by points of recombination. A simulation study was carried out to study the effect of marker spacing, and number of G_{2R} families (i.e., number of randomly spaced points of recombination per target CI) on the Standard Error of QTN map location (SEQTN). For simplicity we simulated a TRPBC (backcross) design, and assumed absence of polygenic effects.
As before, we assume that conventional QTL mapping in a BC, F_{2} or AIL population has detected a QTL of interest with a 95% confidence interval (CI) defined by a pair of flanking upstream and downstream markers, denoted M_{U} and M_{D}, respectively. A very large set of evenly spaced ordered makers, denoted M_{1} to M_{k} spanning the interval M_{U} to M_{D} is available, and haplotypes of the parental lines (denoted the G_{1} generation) with respect to these markers are known. In the proposed scheme, the original or a new mapping population (the G_{2} generation) is genotyped for the M_{U}  M_{D} marker pair, and G_{2R} individuals carrying recombinant chromosomes across this region (i.e., m_{U}M_{D} and M_{U}–m_{D}) are identified. These are backcrossed to one of the original founder lines, and a large simulated BC progeny population (the G_{3} generation) is formed, consisting of a number of G_{2R} families, and genotyped for the marker set spanning the QTL. A QTN is simulated at a specific location in the QTL, and the marker showing the largest difference between alternative marker genotypes is identified as the estimated QTN location. The distance between the estimated QTN location and the simulated location was determined in units of the initial QTL, and the accuracy of map location was evaluated as the SEQTN. For each combination of parameters, 1000 Monte Carlo simulations were run.
To investigate the effect of population size (N), number of internal markers (k) and QTN location (L) on SEQTN we set: N = 1000, 3000, 5000, 7000, 9000, 11000; k = 2, 5, 10, 23, 30, 40, 50 (k does not include the flanking M_{U} and M_{D} markers, e. g., k = 2 represents a total of 4 markers, the two flanking markers and two internal markers); allele substitution effect in standardized units, d = 0.2; number of G_{2R} families, F = 25; QTN location, L = 0.57, 0.77, again taking the width of the CI = 1.0 as the unit of measure. Thus, L = 0.57 is a bit distal to the center of the QTL, and L = 77 is slightly distal to the threequarter mark. QTN positions and the k = 23 marker number were chosen to ensure that in no instance did marker position and simulated QTN position coincide.
To investigate effect of number of families within a QTL on SEQTN we set N = 3000, 11000; F = 500, 250, 100, 50, 25, 10; L = 0.57; k and d as before.
The simulation results were obtained as SEQTN, while the results of the deterministic analysis are presented as the population size required to achieve a given factor of QTL reduction, s. To compare the two approaches we converted both of their outputs to the achieved reduction factor, s. For SEQTN, we assumed that the new confidence interval of QTN location would equal 4*SEQTN (i.e., two standard deviations to each side). Thus, taking the original QTL = 1.0, and the new QTL = 4(SEQTN), the reduction factor would be s = 1/4(SEQTN). For the deterministic analysis, we have N_{TG3} = sN_{G3}, so that s = N_{TG3}/N_{G3}.
For example, in the simulation, at N = 7000 (L = 0.57), we obtain SEQTN = 0.0575. Then s = 1/4(SEQTN) = 4.35. In the deterministic analysis of the TRPBC design, N_{G3} = 1536.64 (see Results). Then, s = 7000/1536.64 = 4.56.
Results
As shown in Methods, the basic parameter determining the required numbers of G_{3} progeny and G_{2} parents under the TRP design, is N_{G3}, the number of G_{3} individuals needed to define a single (1α) CI for given allele substitution effect, d. For a 95% CI in the absence of polygenic effects, we derived the expression N_{G3} = 5.34z_{α/2} ^{2}/d^{2}. Taking 95% CI as the standard, z_{α/2} = 1.96, and folding this into the constant, we have N_{G3} = 20.51/d^{2}. This shows starkly that N_{G3} will be very sensitive to the allele substitution effect, d, e.g., in absence of polygenic effects, N_{G3} = 2051.0, 512.75, and 227.9 for d = 0.1, 0.2 and 0.3, respectively. We stress again that these numbers are independent of the size of the original QTL (C) or of the target subQTL (c), or the reduction factor (s). It makes no difference if the reduction is from 20 cM to 2 cM, or from 2 cM to 0.2 cM, N_{G3} will be the same. The same holds true in the presence of polygenic effects, except that in this case, the basic expression for N_{G3} = 75.91/d^{2}
The total G_{3} population, N_{TG3} required for mapping under TRP, depends solely on N_{G3} and s, the desired reductionfactor in QTL size. A fivefold reduction means that the original QTL will be divided into 5 subQTL. This will require the same number, N_{G3} of individuals for each of the subQTL. Consequently, N_{TG3} = sN_{G3}, where s is the number of subQTL into which we divide our original QTL.
The total G_{2} population required to produce the G_{3} population in absence of polygenic effects, depends on N_{G3}, and on the reproductive potential (n) and target interval (c) as follows, N_{G2} = N_{G3}/2nc(1c). Thus, it stands in proportion to 1/2nc(1c). Assuming, n = 50, 1/2nc(1c) will be quite small when c is large, e.g., for C = 0.2 M, c = 0.1 M, 1/2nc(1c) = 0.111 and N_{G2} = 56.98 indeed very small (5%) relative to N_{TG3} = 1026 for this case (Table 3). However, 1/2nc(1c) will assume larger values when c is small; e.g., for C = 0.1, c = 0.01, 1/2nc(1c) = 0.99, N_{G2} = 507.72 and N_{TG3} = 5129 (10%); and for C = 0.2. c = 0.002, 1/2nc(1c) = 5.01, and N_{G2} = 2569 and N_{TG3} = 5129, (50%).
Table 3 shows population size required for the G_{2} and G_{3} stages of the TRP design and total numbers across both generations as a function of the width of the original QTL (C), the target subQTL (c), and the reduction factor, s, on the assumptions of no polygenic effect, d = 0.2, n = 50 and α = 0.05. Population size for TRP designs relative to F_{2} designs range from 0.289 for initial CI of 0.2 M to 0.041 for initial CI of 0.02 M. Thus, the clear conclusion is that the TRP design will be most useful when the initial CI of the QTL has been brought to fairly high resolution already, and the desired step is ultrahigh resolution. Within a given initial CI, the relative effectiveness of the TRP designs compared to the F_{2} designs does not depend on the target CI. Of course, the actual required numbers vary considerably, as seen in Table 3. But the relative numbers required remain the same.
The table shows clearly that N_{TG3} varies directly with and depends solely on the reduction factor, s; while N_{G2} varies inversely with the original QTL and target subQTL size, and directly with the reduction factor.
In contrast to the required number of G_{3} and G_{2R} individuals, which is a function only of the reduction factor, the total number of G_{2} individuals depends strongly on the target CI. This is due to the fact that it requires, for example, ten times as many total G_{2} individuals to uncover a given number of G_{2R} individuals in a region of 0.02 M, as compared to a region of 0.2 M. Thus, the number of G_{2} individuals ranges from 5% of the number of G_{2} individuals when initial CI is 0.2 M; to ten times this, or 50% when initial CI is 0.02 M.
In the absence of polygenic effects, the TRP delivers major reductions in mapping population size relative to a standard F_{2} design, depending on the initial QTL size. This is particularly evident for very high resolution, c = 0.01M or 0.002M, where the TRP ostensibly delivers ultrahigh resolution at very acceptable population sizes, while the F_{2} design requires very high numbers for mapping at this resolution. Sadly, this fine performance is markedly reduced when polygenic effects are taken into account (Table 4). In this case, the relative effectiveness of the TRP design depends strongly on the initial QTL start width. Indeed, when start point is at C = 0.2, the TRP design required a bit larger population size than an F_{2}. When the initial interval is small (e.g., 0.05 or 0.02 M), TRP requires only about 15% the population size of an F2 for equivalent mapping precision. In this case, TRP can close the gap to 0.01 M with manageable numbers (11,406 and 5,713, for C = 0.05 and C = 0.02, respectively).
For species with short generation interval (annual plants, mice) a two stage TRP can start with a large initial CI, say C = 0.2M, and yet reach c = 0.01 in two steps, with acceptable total population numbers. For example, starting with C = 0.2 M, first stage might reduce fourfold to 0.05 M requiring N_{TRP} = 7991; second stage would reduce fourfold again to 0.0125 M requiring N_{TRP} = 7991, total 15,982, spread fairly equally over two years; while F_{2} would require 30,000.
The number of recombination points per subQTL
The number of G_{2R} individuals needed to produce the required number of G_{3} progeny for each subCI of size c, is:
where
n is the number of progeny that are produced by selfing each G_{2R} parent.
Continuing our example, and assuming n = 50, then N_{G2R} = 512.75/50 = 10.26. The G_{3} progeny of each G_{2R} parent, present a single point of recombination per subQTL Thus, in our example, there will be on average 10.26 points of recombination within each subQTL. The corresponding number when polygenic effects are present, is 37.93. These numbers are constants that depend only on N_{G3}, and are not affected by C, c, or s. The number of recombination points per subQTL when taking polygenic effects into account is quite large, and in view of the simulation results, would seem more than sufficient so as not to be a limiting factor in the precision of setting QTL confidence intervals by the TRP design.
Simulation
Table 5 shows simulation results with respect to effect of mapping population size (N), number of markers spanning the QTL (k) and marker position within the QTL (L). Results are shown as the Standard error of estimated QTN location relative to the simulated location (SEQTN). Also shown are comparisons of reduction factor (s) as obtained from the simulation (sSIM) and as obtained from the deterministic analysis (sDET). Considering first the effect of marker number, except for N = 1000, at both QTN positions SEQTN is reduced by each step from k = 2 to k = 23, but there is no further decrease in going from k = 23 to k = 50. (For N = 1000 there is a further reduction going from k = 23 to k = 30, but not beyond this). This is somewhat surprising, as we would expect an interaction between population size and marker density. The lack of such interaction is probably due to the fact that number of families was set at 25, so that the number of recombination points across the QTL was a limiting factor in reducing SEQTN by increase in marker density.
Each step increase in N, resulted in a decrease in SEQTN. As could be expected, the marginal decrease in SEQTL was less with each additional step, but still appreciable until the last; the gain in going from N = 9000 to N = 11000 might appear not to be worth the cost. However, this may be a further consequence of the limit of G_{2R} families to 25. In practice, increase in N would be accompanied by an increase in N_{G2R}, so that the combined effect of both might be appreciable (see further discussion of Table 6).
Placing the QTN closer to the boundary of the QTL (L = 0.77 compared to L = 0.57) gave slightly smaller SEQTL at all population sizes. This was higher for smaller N and smaller k, and decreased as N and k increased. This is an artifact, due to the fact that in the simulation, the estimated QTN position could not fall outside the boundaries of the QTL. Hence, the boundary set an artificial upper limit to the simulated errors, reducing the SEQTN accordingly. However, the effect was small, suggesting that the results for the more central QTN location (L = 0.57) are not affected in a major way by boundary effects.
Table 5 also shows reduction factor (s) achieved at k = 50, F = 25 by the simulation (sSIM) and by the deterministic analysis (sDET). At marker saturation (k = 50), correspondence is rather close for intermediate population sizes (N = 5000 and N = 7000) but fall off to either side. On the low side (N = 1000 and 3000) this is probably due to the aforementioned boundary effect. This affects the simulation results, placing an upper boundary on the error values; but apparently does not affect the deterministic results. On the high side (N = 9000 and 11000) SEQTL for the simulation is probably limited by the small number of G_{2R} parents, and the resulting paucity of recombination points for mapping, while this is not a limitation for the deterministic analyses.
Table 6 shows simulation results with respect to number of G_{2R} families in interaction with population size and number of markers. Generally, there is a clear and major decrease in SEQTN with increase in number of families for given population size. The largest effect on SEQTL is given by the first step, from F = 10 to F = 25, but even at N = 11000, each additional step results in a decrease in SEQTN. SEQTN at F = 500, is just half the SEQTN at F = 25. At N = 3000 there appears to be a slight optimum at F = 250. This may be due to inverse relation of number of families and number of individuals within families. The expected interaction between population size and marker density can be observed: at N = 3000, increase of k from 30 to 50 does not reduce SEQTN; but at N = 11000, from F = 50 on, there is a clear reduction in SEQTN with increased marker density. Thus, as suggested above, by increasing F along with N, major reductions in SEQTN may be achieved, even by going beyond N = 11000. For N = 11000 and F = 50 or more, sSIM is greater than sDET, to an increasing degree with greater F. For example, with F = 500, sSIM = 10.5, while sDET = 7.16. We do not have an intuitive explanation for this effect.
Discussion
The TRP is aimed at mapping a specific target QTN, previously assigned to a rather broad confidence interval; to a smaller CI, reduced by some desired factor relative to the original CI. As such, it is a Group II design as defined by Darvasi [2]. A TRP mapping population of size 10,000 individuals can reduce CI by a factor of 4 or 5, even for a QTN of relatively small standardized effect (d = 0.2); for a strong QTN (d = 0.5), the same reduction would be provided by a TRP population size under 2000. A twostage experiment with fourfold reduction at each stage, would reduce a starting CI of 20 cM to 5 cM in the first stage, and to 1.25 cM in the second stage. For d = 0.2, this would require about 16,000 individuals in two cohorts of 8,000 each. For d = 0.5 this would be achieved with two cohorts of about 2000 each. Thus, the TRP design provides a useful solution to the challenge of achieving high resolution mapping for a known specific target QTN. A unique property of the TRP design, is the fact that mapping population size required for a given reduction factor is independent of the size of the starting CI. It will be the same whether the starting CI is 20 cM or 2 cM. Consequently, at a given mapping population size, the smaller the initial CI, the smaller the final CI in proportion. For maximum effectiveness, as shown by the simulation, the TRP design should be implemented with at least 3040 evenly spaced markers and at least 50 to 100 recombination points within the original CI.
When a QTL is mapped to a large CI, it is not possible to tell if the effect is due to a single QTN in the region, or to the summed effect of a number of closely linked QTN. By dividing the initial CI into subintervals, each with considerable mapping power, the TRP design can distinguish between an effect localized to a single point in the original CI, indicating that the total effect was produced by a single QTN, and an effect that is spread all through the QTL, indicating that the total effect represents the summed effects of a number of QTN.
For outcrossing species, WGAS appears to provide a satisfactory solution to the challenge of high resolution mapping of QTN with appreciable effects. This option is obviously not available for high resolution mapping in pure lines of selfing species, or inbred lines of outcrossing species. Although there are a plethora of Group I designs that can provide high resolution mapping across the entire genome (reviewed in Background), with the exception of the AIL design [3] these are all based on specialized resources specifically constructed for high resolution mapping and are perforce limited to the QTL segregating among the founder lines of the resource. For Group II designs that deal with a previously mapped QTL on an ad hoc basis, the choice of designs is limited to the ICSC and RPT designs. The TRP design is an addition to this group. Strictly speaking, the NIL design [9] is not a procedure for high or ultrahigh resolution mapping, since the isolated donor segments will likely be in the range 1020 cM. However, any of the NIL sublines, carrying a QTL of interest would be a superb start point for subsequent application of the ICSC, RPT or TRP designs. There is no doubt that when applicable the ICSC design is a very effective means for high resolution mapping, with minimal requirements for population size. The limitation of the ICSC design is the large number of backcross generations required to generate the series of congenic strains spanning the original CI. This is not a major limitation for species such as the mouse with rapid reproduction cycles; but it is a major limitation for the many plant species that have a single growing season a year. Selection within each backcross step for individuals carrying the least amount of donor genome, could probably reduce the number of backcross steps by one or two generations. The RPT design was presented only briefly in [2] and does not appear to have been subjected to detailed analysis, or to have been widely applied. It appears to be very similar to the TRP design proposed and analyzed here. In both designs, individuals that are recombinant within the original CI are chosen as parents of the mapping population. The major difference is that the RPT design works with a limited number of families, each representing one of a laddered series of recombination points spanning the original QTL at 1 cM intervals. Each of the families must be sufficiently powerful to give a decisive decision as to its QTL status and in this way determine whether the QTN is upstream or downstream of the recombination point in that family. In the TRP design, these requirements are relaxed, and the only requirement is that the total number of families is sufficient to give adequate coverage of the original QTL by recombination points, and the total number of recombinant progeny across all families is sufficient for high resolution mapping across the original CI. It should also be noted that the RPT design may have to base on a backcross design in order to obtain the very large numbers required for high power of the individual family, while the TRP design works well with an F_{2} design and fewer individuals per family. Since the F2 design is generally more powerful than the BC design with equivalent numbers (as seen in the present study), this might also favor the TRP. More exact comparison of the two designs remains for deeper analysis of the RPT. However, all in all RPT and TRP can be considered as variants of the same basic design.
Conclusions
TRP design allows reducing confidence interval of a known target QTL by some desired factor, independent of the original QTL size. The population size required to achieve this depends greatly on the allele substitution effect and also on the polygenic effect and the factor of reduction that was chosen. It is most effective for reducing CI from high resolution (CI = 0.02  0.05 Morgan) to ultrahigh resolution (0.002  0.005 Morgan) since, as opposed to the standard designs, the number of individuals required does not depend strongly on the size of the target CI. Even in the presence of polygenic effects, TRP provides opportunities to achieve CI reduction with a manageable population size where F_{2} and other designs fail. The TRP design also saves greatly in reducing the amount of genotyping required. Thus, the TRP design provides a useful solution to the problem of achieving high and ultrahighresolution mapping in crosses of inbred or pure lines, where genome wide associations tests are not applicable.
Appendix
Calculating N_{G3} for the TRPBC design
As shown by Weller and Soller [26] if M_{1} is a marker located at an estimated QTN location, the probability that the CI of QTN location includes the marker M_{2} located at a remove of L Morgan from M_{1}, is equal to the probability of obtaining the value
where ,
Z_{α/2} is the standard normal variable corresponding to a probability of α/2,
D = E(M_{1})  E(M_{2}), where E(M_{1}) is the expected effect at M_{1} (QTN located at the marker), E(M_{2}) is the expected QTN effect at M_{2} (QTL located at a remove from the QTN) considering recombinant individuals only; and SE(D) is the standard error of D.
From Table 7 it is apparent that there are two marker genotype groups, each of which is composed of a single recombinant marker group. Marker genotypegroup M_{1}m_{1} is composed of recombinant genotypegroup A of Table 7 with genotype: M_{1}m_{2}/m_{1}m_{2}, having genotypic value h; frequency 1/4 of the entire G3 population. Marker genotypegroup M_{2}m_{2} is composed of recombinant genotypegroup B of Table 7 with genotype: m_{1}M_{2}/m_{1}m_{2}, having genotypic value d; frequency 1/4 of the entire G3 population.
Letting italics denote the mean genotypic value of the corresponding marker genotype group (including recombinant genotypes only), we have
On this basis, letting A, B, represent the genotypic values of the corresponding genotype groups, we have
Substituting in (1) and combining like terms, we obtain
Substituting genotypic values of the recombinant groups we have
To calculate SE(D), we assume that variance within the G_{3} families is the same as variance within the F_{2} generation (set to 1.0 for standardization). This is conservative, as the genetic variance within G_{3} BC families will be less than in an F_{2} population, depending on degree of homzoygosity in their G_{2R} parent. On this basis, we can write:
(as each of these genotype groups comprises 1/4 of the G_{3} mapping population) of size N_{G3}
Then,
Substituting, we have
If we set α=0.05 and d=0.2, as we did for the TRPF2 design, we have Z_{α/2}=1.96, and
Since the interval between M_{1} and M_{2} defines only half of the CI of QTN location, to cover the entire CI would require twice this
On these assumed values for α and d, the basic G_{3R} population required for CI(0.95) = N_{G3}= 1536.64 This number is constant for given d, and α, and does not depend on the target size of the confidence interval. The corresponding value for TRPF_{2} design is 512.34. On general principles we would have expected NG3 for TRPBC with h=0 to be twice that for TRPF_{2} (that is, 1024,68). The difference is due to the fact that the mapping population for the TRPF_{2} is 75% of the total population; while the mapping population for the TRPBC is 50% of the total population. If we increase the TRPBC N_{G3} value by 50% to make up for this we obtain NG3= 1529.37, as found.
Abbreviations
 QTN:

Quantitative Trait Nucleotide
 QTL:

Quantitative Trait Locus
 TRP:

Targeted recombinant Progeny
 CI:

Confidence Interval
 MAS:

Marker assisted selection
 MAI:

Marker assisted introgression
 SNP:

Single Nucleotide Polymorphism
 WGAS:

Whole genome association studies
 LD:

Linkage disequilibrium
 AIL:

Advanced intercross line
 HS:

Heterogeneous Stock
 RIL:

Recombinant Inbred Lines
 MAGIC:

Multiparent Advanced Generation InterCross
 RPT:

Recombinant Progeny Testing
 ISCS:

IntervalSpecific Congenic Strains
 g.d.p.:

genotyping data points
 SEQTN:

Standard error of QTN
References
Carlson CS, Eberle MA, Kruglyak L, Nickerson DA. Mapping complex disease loci in wholegenome association studies. Nature. 2004;429:446–52.
Darvasi A. Experimental strategies for the genetic dissection of complex traits in animal models. Nat Genet. 1998;18:19–24.
Darvasi A, Soller M. Advanced intercross lines, an experimental population for fine genetic mapping. Genetics. 1995;141:1199–207.
Iraqi F, Clapcot S, Kuman P, Heley C, Kemp S, Teale A. Fine mapping of trypanosomiasis resistance QTLs in mice using advanced intercross lines. Mamm Genome. 2000;11:645–8.
Mott R, Talbot CJ, Turri MG, Collins AC, Flint J. A method for fine mapping quantitative trait loci in outbred animal stocks. PNAS. 2000;97:12649–54.
Svenson KL, Gatti DM, Valdar W, Welsh CE, Cheng R, Chesler EJ, et al. Highresolution genetic mapping using the mouse diversity outbred population. Genetics. 2012;190:437–47.
Darvasi A. Dissecting complex traits: the geneticists 'Around the world in 80 days'. Trends Genet. 2005;21:373–6.
Yalcin B, WillisOwen SAG, Fullerton J, Meesaq A, Deacon RM, Rawlins JNP, et al. Genetic dissection of a behavioral quantitative trait locus shows that Rgs2 modulates anxiety in mice. Nat Genet. 2004;36:1197–202.
Falke KC, Frisch M. Power and false positive rate in QTL detection with nearisogenic line libraries. Heredity. 2011;106:576–84.
Paterson AH, DeVerna JW, Lanini B, Tanksley SD. Fine mapping of quantitative trait loci using selected overlapping recombinant chromosomes in an interspecific cross of tomato. Genetics. 1990;124:735–42.
Eshed Y, Zamir D. An introgression line population of Lycopersicon penelli in the cultivated tomato enables the identification and fine mapping of yield associated QTL. Genetics. 1995;141:1147–62.
Churchill GA, Airey DC, Allayee H, Angel JM, Attie AD, Beatty J, et al. Complex trait consortium (2004). the collaborative cross, a community resource for the genetic analysis of complex traits. Nat Genet. 2004;36:1133–7.
Iraqi FA, Athamni H, Dorman A, Salymah Y, Tomlinson I, Nashif A, et al. Heritability and coefficient of genetic variation analyses of phenotypic traits provide strong basis for highresolution QTL mapping in the Collaborative Cross mouse genetic reference population. Mamm Genome. 2014;25:109–19.
Cavanagh C, Morell M, Mackay I, Powell W. From mutations to MAGIC: resources for gene discovery, validation and delivery in crop plants. Curr Opin Plant Biol. 2008;11:215–21.
Kover PX, Valdar W, Trakalo J, Scarcelli N, Ehrenreich IM, Purugganan MD, et al. A multiparent advanced GenerationIntercross to finemap quantitative traits in Arabidopsis thaliana. PLoS Genet. 2009;5:e1000551.
Huang BE, George AW, Forrest KL, Kilian A, Hayden MJ, Morell MK, et al. A multiparent advanced generation intercross population for genetic analysis in wheat. Plant Biotechnol J. 2012;10:826–39.
Bandillo N, Raghavan C, Muyco PA, Sevilla MAL, Lobina IT, DillaErmita CJ, et al. Multiparent advanced generation intercross (MAGIC) populations in rice: progress and potential for genetics research and breeding. Rice. 2013;6:11.
Rosa GJM, Leon ND, Rosa AJM. Review of microarray experimental design strategies for genetical genomics studies. Physiol Genomics. 2006;28:15–23.
Jin C, Lan H, Attie AD, Churchill GA, Bulutuglo D, Yandell S. Selective phenotyping for increased efficiency in genetic mapping studies. Genetics. 2004;168:2285–93.
Xu Z, Zou F, Vision TJ. Improving quantitative trait loci mapping resolution in experimental crosses by the Use of genotypically selected samples. Genetics. 2005;170:401–8.
Spickett SP, Thoday JM. Regular responses to selection; 3. Interaction between located polygenes. Genet Res. 1966;7:96–121.
Liu X, Oliver F, Brown SDM, Denny P, Keightley PD. Highresolution quantitative trait locus mapping for body weight in mice by recombinant progeny testing. Genet Res. 2001;77:191–7.
Darvasi A. Intervalspecific congenic strains (ISCS): an experimental design for mapping a QTL into a 1centimorgan interval. Mamm Genome. 1997;8:163–7.
Wang L, Jiao Y, Huang Y, Bennett B, Williams RW, Li D, et al. Ttn as a likely causal gene for QTL of alcohol preference on mouse chromosome 2. BMC Bioinformatics. 2014;15 Suppl 10:12.
Ronin Y, Korol A, Shtemberg M, Nevo E, Soller M. Highresolution mapping of quantitative trait loci by selective recombinant genotyping. Genetics. 2003;164:1657–66.
Weller JI, Soller M. An analytical formula to estimate confidence interval of QTL location with a saturated genetic map as a function of experimental design. Theor Appl Genet. 2004;109:1224–9.
Darvasi A, Soller M. Selective genotyping for determination of linkage between a marker locus and a quantitative trait locus. Theor Appl Genet. 1992;85:353–9.
Darvasi A, Soller M. Selective DNA pooling for determination of linkage between a molecular marker and a quantitative trait locus. Genetics. 1994;138:1365–73.
Lipkin E, Mosig MO, Darvasi A, Ezra E, Shalom A, Friedmann A, et al. Mapping loci controlling milk protein percentage in dairy cattle by means of selective milk DNA pooling using dinucleotide microsatellite markers. Genetics. 1998;149:1557–67.
Mather K, Jinks JL. Biometrical Genetics. London: Chapmand and Hall; 1970.
Acknowledgements
Except for EH and MS, no one else made a substantial contribution to the conception, design, acquisition of data, or analysis or was involved in drafting the manuscript or revising it critically for important intellectual content. Preparation of the manuscript was selffunded and did not include a language editor or scientific writer. The present study itself was selffunded and was not supported by a grant from a funding agency. Publication page costs were provided by the Research Authority of JCT  Lev Academic Center.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
MS conceived of the study, and participated in its design and implementation and helped to draft the manuscript. EMH wrote and ran the simulation, verified the analyses and participated actively in all stages of manuscript preparation. Both authors read and approved the final manuscript.
Rights and permissions
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Heifetz, E.M., Soller, M. Targeted Recombinant Progeny: a design for ultrahigh resolution mapping of Quantitative Trait Loci in crosses between inbred or pure lines. BMC Genet 16, 76 (2015). https://doi.org/10.1186/s128630150206z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s128630150206z