Effects of missing marker and segregation distortion on QTL mapping in F2 populations

Missing marker and segregation distortion are commonly encountered in actual quantitative trait locus (QTL) mapping populations. Our objective in this study was to investigate the impact of the two factors on QTL mapping through computer simulations. Results indicate that detection power decreases with increasing levels of missing markers, and the false discovery rate increases. Missing markers have greater effects on smaller effect QTL and smaller size populations. The effect of missing markers can be quantified by a population with a reduced size similar to the marker missing rate. As for segregation distortion, if the distorted marker is not closely linked with any QTL, it will not have significant impact on QTL mapping; otherwise, the impact of the distortion will depend on the degree of dominance of QTL, frequencies of the three marker types, the linkage distance between the distorted marker and QTL, and the mapping population size. Sometimes, the distortion can result in a higher genetic variance than that of non-distortion, and therefore benefits the detection of linked QTL. A formula of the ratio of genetic variance explained by QTL under distortion and non-distortion was given in this study, so as to easily determine whether the segregation distortion marker (SDM) increases or decreases the QTL detection power. The effect of SDM decreases rapidly as its linkage relationship with QTL becomes looser. In general, distorted markers will not have a great effect on the position and effect estimations of QTL, and their effects can be ignored in large-size mapping populations.     

Prevalence of segregation distortion in diploid alfalfa and its implications for genetics and breeding applications

Segregation distortion (SD) is often observed in plant populations; its presence can affect mapping and breeding applications. To investigate the prevalence of SD in diploid alfalfa (Medicago sativa L.), we developed two unrelated segregating F₁ populations and one F₂ population. We genotyped all populations with SSR markers and assessed SD at each locus in each population. The three maps were syntenic and largely colinear with the Medicago truncatula genome sequence. We found genotypic SD for 24 and 34% of markers in the F₁ populations and 68% of markers in the F₂ population; distorted markers were identified on every linkage group. The smaller percentage of genotypic SD in the F₁ populations could be because they were non-inbred and/or due to non-fully informative markers. For the F₂ population, 60 of 90 mapped markers were distorted, and they clustered into eight segregation distortion regions (SDR). Most SDR identified in the F₁ populations were also identified in the F₂ population. Genotypic SD was primarily due to zygotic rather than allelic distortion, suggesting zygotic not gametic selection is the main cause of SD. On the F₂ linkage map, distorted markers in all SDR except two showed heterozygote excess. The severe SD in the F₂ population likely biased genetic distances among markers and possibly also marker ordering and could affect QTL mapping of agronomic traits. To reduce the effects of SD and non-fully informative markers, we suggest constructing linkage maps and conducting QTL mapping in advanced generation populations.     

Multi-QTL mapping for quantitative traits using distorted markers

Marker segregation distortion is a common natural phenomenon. However, relatively little is known about utilizing distorted markers for detecting quantitative trait loci (QTL). Therefore, in this study we proposed a multi-QTL mapping approach that uses distorted markers. First, the information from all markers, including distorted markers, was used to detect segregation distortion loci (SDL). Second, the information from the detected SDL was used to correct the conditional probabilities of the QTL genotypes conditional on marker information, and these corrected probabilities were then incorporated into a multi-QTL mapping methodology. Finally, the proposed approach was validated by both Monte Carlo simulation studies and real data analysis. The results from the simulation studies show that as long as one or two SDL are placed around the simulated QTL, there are no differences between the new method and the ordinary interval mapping method in terms of the power of QTL detection or the estimates of the position and dominant effects of the QTL. However, the power of QTL detection is higher under the dominant genetic model of SDL than under the additive genetic model, and the estimate for the additive effect of QTL using the new method is significantly different from the estimate obtained using ordinary interval mapping. The above results were further confirmed by the detection of QTL for dried soymilk in 222 F_2:4 families in soybean.     

Model selection for quantitative trait loci mapping in a full-sib family

Statistical methods for mapping quantitative trait loci (QTLs) in full-sib forest trees, in which the number of alleles and linkage phase can vary from locus to locus, are still not well established. Previous studies assumed that the QTL segregation pattern was fixed throughout the genome in a full-sib family, despite the fact that this pattern can vary among regions of the genome. In this paper, we propose a method for selecting the appropriate model for QTL mapping based on the segregation of different types of markers and QTLs in a full-sib family. The QTL segregation patterns were classified into three types: test cross (1:1 segregation), F₂ cross (1:2:1 segregation) and full cross (1:1:1:1 segregation). Akaike’s information criterion (AIC), the Bayesian information criterion (BIC) and the Laplace-empirical criterion (LEC) were used to select the most likely QTL segregation pattern. Simulations were used to evaluate the power of these criteria and the precision of parameter estimates. A Windows-based software was developed to run the selected QTL mapping method. A real example is presented to illustrate QTL mapping in forest trees based on an integrated linkage map with various segregation markers. The implications of this method for accurate QTL mapping in outbred species are discussed.     

Multi-QTL Mapping for Quantitative Traits Using Epistatic Distorted Markers

The interaction between segregation distortion loci (SDL) has been often observed in all kinds of mapping populations. However, little has been known about the effect of epistatic SDL on quantitative trait locus (QTL) mapping. Here we proposed a multi-QTL mapping approach using epistatic distorted markers. Using the corrected linkage groups, epistatic SDL was identified. Then, these SDL parameters were used to correct the conditional probabilities of QTL genotypes, and these corrections were further incorporated into the new QTL mapping approach. Finally, a set of simulated datasets and a real data in 304 mouse F₂ individuals were used to validate the new method. As compared with the old method, the new one corrects genetic distance between distorted markers, and considers epistasis between two linked SDL. As a result, the power in the detection of QTL is higher for the new method than for the old one, and significant differences for estimates of QTL parameters between the two methods were observed, except for QTL position. Among two QTL for mouse weight, one significant difference for QTL additive effect between the above two methods was observed, because epistatic SDL between markers C66 and T93 exists (P = 2.94e-4).     

Genetic linkage maps of rose constructed with new microsatellite markers and locating QTL controlling flowering traits

New microsatellites markers [simple sequence repeat (SSR)] have been isolated from rose and integrated into an existing amplified fragment-length polymorphism genetic map. This new map was used to identify quantitative trait locus (QTL) controlling date of flowering and number of petals. From a rose bud expressed sequence tag (EST) database of 2,556 unigenes and a rose genomic library, 44 EST-SSRs and 20 genomic-SSR markers were developed, respectively. These new rose SSRs were used to expand genetic maps of the rose interspecific F₁ progeny. In addition, SSRs from other Rosaceae genera were also tested in the mapping progeny. Genetic maps for the two parents of the progeny were constructed using pseudo-testcross mapping strategy. The maps consist of seven linkage groups of 105 markers covering 432 cM for the maternal map and 136 markers covering 438 cM for the paternal map. Homologous relationships among linkage groups between the maternal and paternal maps were established using SSR markers. Loci controlling flowering traits were localised on genetic maps as a major gene and QTL for the number of petals and a QTL for the blooming date. New SSR markers developed in this study will provide tools for the establishment of a consensus linkage map for roses that combine traits and markers in various rose genetic maps.     

Linkage analysis of anther-derived monoploids showing distorted segregation of molecular markers

Monoploids can be obtained from several diploid plant species by anther culture. Mapping of molecular markers using monoploids is greatly facilitated by the simple 1:1 segregation ratio expected from all heterozygous loci in the genome. Distorted segregation of molecular markers, however, appears to be a common phenomenon in many crop species and hinders the use of monoploids for mapping purposes. This report examines the segregation pattern of two marker genes linked together with one locus or separately with two independent loci which are responsible for the observed distortion. Each of the loci exhibiting distorted segregation has one of the two alleles which inhibits regeneration of the gametic cells in vitro and disrupts the expected segregation ratio of the linked markers. All possible situations in which linkage occurs between markers and distortion-causing genes are considered. Theoretical results outlining the segregation pattern among these linkage types indicate that the distinguishable distorted ratios can be used for mapping purposes. A protocol is given for the mapping of distorted gene markers based on existing gene mapping software. An example is presented of the mapping of distorted RAPD markers of monoploids obtained from a diploid potato genotype. Received: 18 October 1999 / Accepted: 24 November 1999<@head-com-p1a.lf>Communicated by G. Wenzel     

Quantitative trait locus mapping with background control in genetic populations of clonal F1 and double cross

In this study, we considered five categories of molecular markers in clonal F₁ and double cross populations, based on the number of distinguishable alleles and the number of distinguishable genotypes at the marker locus. Using the completed linkage maps, incomplete and missing markers were imputed as fully informative markers in order to simplify the linkage mapping approaches of quantitative trait genes. Under the condition of fully informative markers, we demonstrated that dominance effect between the female and male parents in clonal F₁ and double cross populations can cause the interactions between markers. We then developed an inclusive linear model that includes marker variables and marker interactions so as to completely control additive effects of the female and male parents, as well as the dominance effect between the female and male parents. The linear model was finally used for background control in inclusive composite interval mapping (ICIM) of quantitative trait locus (QTL). The efficiency of ICIM was demonstrated by extensive simulations and by comparisons with simple interval mapping, multiple‐QTL models and composite interval mapping. Finally, ICIM was applied in one actual double cross population to identify QTL on days to silking in maize.     

Construction of integrated genetic linkage maps of the tiger shrimp (Penaeus monodon) using microsatellite and AFLP markers

The linkage maps of male and female tiger shrimp (P. monodon) were constructed based on 256 microsatellite and 85 amplified fragment length polymorphism (AFLP) markers. Microsatellite markers obtained from clone sequences of partial genomic libraries, tandem repeat sequences from databases and previous publications and fosmid end sequences were employed. Of 670 microsatellite and 158 AFLP markers tested for polymorphism, 341 (256 microsatellite and 85 AFLP markers) were used for genotyping with three F₁ mapping panels, each comprising two parents and more than 100 progeny. Chi‐square goodness‐of‐fit test (χ²) revealed that only 19 microsatellite and 28 AFLP markers showed a highly significant segregation distortion (P < 0.005). Linkage analysis with a LOD score of 4.5 revealed 43 and 46 linkage groups in male and female linkage maps respectively. The male map consisted of 176 microsatellite and 49 AFLP markers spaced every ～11.2 cM, with an observed genome length of 2033.4 cM. The female map consisted of 171 microsatellite and 36 AFLP markers spaced every ～13.8 cM, with an observed genome length of 2182 cM. Both maps shared 136 microsatellite markers, and the alignment between them indicated 38 homologous pairs of linkage groups including the linkage group representing the sex chromosome. The karyotype of P. monodon is also presented. The tentative assignment of the 44 pairs of P. monodon haploid chromosomes showed the composition of forty metacentric, one submetacentric and three acrocentric chromosomes. Our maps provided a solid foundation for gene and QTL mapping in the tiger shrimp.     

白桦AFLP遗传连锁图谱的构建

以80个中国白桦(Betula platyphylla Suk)×欧洲白桦(Betula pendula Roth)的F1个体为作图群体, 利用扩增片段长度多态性(Amplified fragment length polymorphism, AFLP)标记, 按照拟测交作图策略, 分别构建了中国白桦和欧洲白桦的分子标记遗传连锁图谱。从64对AFLP引物组合中筛选出34对多态性丰富的引物组合, 这些入选的引物组合在分离群体中共检测到451个多态性位点。χ2检验结果表明, 有362个位点符合1∶1分离(拟测交分离位点), 41个位点符合3∶1分离, 20个位点符合1∶3分离, 28个位点属偏分离位点。在符合拟测交分离的位点中, 201个位点来自中国白桦, 161个位点来自欧洲白桦。利用2点连锁分析, 来自中国白桦的201个标记构成了14个连锁群(4个以上标记), 10个三连体和14个连锁对, 45个为非连锁位点, 连锁标记覆盖的总图距为1 296.1 cM, 平均图距15.5 cM。而来自欧洲白桦的161个标记构成了17个不同的连锁群(4个以上标记), 8个三连体和4个连锁对, 15个为非连锁位点, 连锁标记覆盖的总图距为1 035.8 cM, 平均图距12 cM。     

The effect of an imprecise map on interval mapping QTLs

The statistical analysis of quantitative trait locus (QTL) experiments relies on the use of a linkage map of the markers genotyped. Such a map is, at best, a good estimate of the true map. Resources might be diverted into developing better marker maps or improved maps become available after the analysis, raising concerns over the original analysis. It is therefore important to understand the sensitivity of QTL analysis to map inaccuracy. We have used simulation methods to investigate the consequences of an incorrect map on the results of a QTL analysis using interval mapping. Backcross data sets were generated with a particular map and then analysed with both the correct map and incorrect maps. If the incorrect maps maintained the true linkage groups (i.e. no markers were incorrectly assigned to another linkage group), the accuracy of the map had little or no impact on the ability to detect QTLs, the true significance levels of the tests or the relative placement of QTLs. When a marker was incorrectly placed on another linkage group, there was a small increase in the level of the test. After adjusting for this increase, there was a decrease in power to detect a QTL near the misplaced marker. This decrease was of a similar magnitude to that found when using a single-marker analysis compared with interval mapping. These results mean that QTL analyses can proceed without the need for very accurate marker maps, and that estimated QTL positions can be translated onto updated maps without the need for reanalysis.     

A model for quantitative trait loci mapping,linkage phase,and segregation pattern estimation for a full-sib progeny

Quantitative trait loci (QTL) mapping is an important approach for the study of the genetic architecture of quantitative traits. For perennial species, inbred lines cannot be obtained due to inbreed depression and a long juvenile period. Instead, linkage mapping can be performed by using a full-sib progeny. This creates a complex scenario because both markers and QTL alleles can have different segregation patterns as well as different linkage phases between them. We present a two-step method for QTL mapping using full-sib progeny based on composite interval mapping (i.e., interval mapping with cofactors), considering an integrated genetic map with markers with different segregation patterns and conditional probabilities obtained by a multipoint approach. The model is based on three orthogonal contrasts to estimate the additive effect (one in each parent) and dominance effect. These estimatives are obtained using the EM algorithm. In the first step, the genome is scanned to detect QTL. After, segregation pattern and linkage phases between QTL and markers are estimated. A simulated example is presented to validate the methodology. In general, the new model is more effective than existing approaches, because it can reveal QTL present in a full-sib progeny that segregates in any pattern present and can also identify dominance effects. Also, the inclusion of cofactors provided more statistical power for QTL mapping.     

Comparison between Poncirus and Citrus genetic linkage maps

Five genetic linkage maps were constructed for the parents of three progenies: Citrus aurantium (A) x Poncirus trifoliata var. Flying Dragon (Pa), C. volkameriana (V) x P. trifoliata var. Rubidoux (Pv) and a self-pollination of P. trifoliata var. Flying Dragon (Pp). The number of polymorphic markers assayed ranged from 48 for Pa to 120 for A according to the heterozygosity of each parental. As our focus was on genome comparison, most of the markers were newly generated simple sequence repeats. Inter-retrotransposon amplified polymorphisms based on four retrotransposon sequences isolated from Citrus spp were also used to saturate the maps. These polymorphisms were much more frequent in A (53) than in Pa (15) and randomly distributed throughout both genomes. Since comparative genomics and quantitative trait locus analysis applicability depends on the reliability of marker ordering, the causes of variation in marker order were investigated. Around 25% of the markers showed gametal segregation distortions. Segregation distortions were also observed at the zygotic level towards a reduction in the observed frequency of homozygotes from that expected in linkage groups 5 and 7. The presence of balanced lethal factors or gametal incompatibility genes in those genomic regions would explain a zygotic advantage of heterozygotes at these specific regions. Four differences in genomic organization were observed; three are putative translocations and affect homeologous linkage groups 3, 7 and 11, where highly distorted markers are found. Other causes of variation in marker order are also discussed: the introduction of new markers in the map, lowering the LOD score and the mapping software. These results represent the first comparative mapping analysis among Citrus and Poncirus species.     

Identification of QTLs Conferring Resistance to Deltamethrin in Culex pipiens pallens

Culex pipiens pallens is the most abundant Culex mosquito species in northern China and is an important vector of bancroftian filariasis and West Nile virus. Deltamethrin is an insecticide that is widely used for mosquito control, however resistance to this and other insecticides has become a major challenge in the control of vector-borne diseases that appear to be inherited quantitatively. Furthermore, the genetic basis of insecticide resistance remains poorly understood. In this study, quantitative trait loci (QTL) mapping of resistance to deltamethrin was conducted in F₂ intercross segregation populations using bulked segregation analysis (BSA) and amplified fragment length polymorphism markers (AFLP) in Culex pipiens pallens. A genetic linkage map covering 381 cM was constructed and a total of seven QTL responsible for resistance to deltamethrin were detected by composite interval mapping (CIM), which explained 95% of the phenotypic variance. The major QTL in linkage group 2 accounted for 62% of the variance and is worthy of further study. 12 AFLP markers in the map were cloned and the genomic locations of these marker sequences were determined by applying the Basic Local Alignment Search Tool (BLAST) tool to the genome sequence of the closely related Culex quinquefasciatus. Our results suggest that resistance to deltamethrin is a quantitative trait under the control of a major QTL in Culex pipiens pallens. Cloning of related AFLP markers confirm the potential utility for anchoring the genetic map to the physical map. The results provide insight into the genetic architecture of the trait.     

A genetic linkage map of water yam ( Dioscorea alata L.) based on AFLP markers and QTL analysis for anthracnose resistance

A genetic linkage map of the tetraploid water yam (Dioscorea alata L.) genome was constructed based on 469 co-dominantly scored amplified fragment length polymorphism (AFLP) markers segregating in an intraspecific F₁ cross. The F₁ was obtained by crossing two improved breeding lines, TDa 95/00328 as female parent and TDa 87/01091 as male parent. Since the mapping population was an F₁ cross between presumed heterozygous parents, marker segregation data from both parents were initially split into maternal and paternal data sets, and separate genetic linkage maps were constructed. Later, data analysis showed that this was not necessary and thus the combined markers from both parents were used to construct a genetic linkage map. The 469 markers were mapped on 20 linkage groups with a total map length of 1,233 cM and a mean marker spacing of 2.62 cM. The markers segregated like a diploid cross-pollinator population suggesting that the water yam genome is allo-tetraploid (2n = 4x = 40). QTL mapping revealed one AFLP marker E-14/M52-307 located on linkage group 2 that was associated with anthracnose resistance, explaining 10% of the total phenotypic variance. This map covers 65% of the yam genome and is the first linkage map reported for D. alata. The map provides a tool for further genetic analysis of traits of agronomic importance and for using marker-assisted selection in D. alata breeding programmes. QTL mapping opens new avenues for accumulating anthracnose resistance genes in preferred D. alata cultivars.     

Mapping of the QTL (quantitative trait locus) conferring partial resistance to leaf blast in rice cultivar Chubu 32

The rice cultivar Chubu 32 possesses a high level of partial resistance to leaf blast. The number and chromosomal location of genes conferring this resistance were detected by restriction fragment length polymorphism (RFLP) linkage mapping and quantitative trait locus (QTL) analysis. For the mapping, 149 F₃ lines derived from the cross between rice cultivar Norin 29, with a low level of partial resistance, and Chubu 32 were used, and their partial resistance to leaf blast was assessed in upland nurseries. A linkage map covering six chromosomes and consisting of 36 RFLP markers was constructed. In the map, only one significant QTL (LOD>2.0) for partial resistance was detected on chromosome 11. This QTL explained 45.6% of the phenotypic variation. The segregation ratio of the F3 lines was 3:1 for partial resistance to susceptibility. These results suggest that the partial resistance in Chubu 32 is controlled by a major gene. Received: 15 March 2001 / Accepted: 13 August 2001     

A genetic linkage map based on AFLP and SSR markers and mapping of QTL for dry-matter content in sweetpotato

We developed a genetic linkage map of sweetpotato using amplified fragment length polymorphism (AFLP) and simple sequence repeat (SSR) markers and a mapping population consisting of 202 individuals derived from a broad cross between Xushu 18 and Xu 781, and mapped quantitative trait loci (QTL) for the storage root dry-matter content. The linkage map for Xushu 18 included 90 linkage groups with 2077 markers (1936 AFLP and 141 SSR) and covered 8,184.5 cM with an average marker distance of 3.9 cM, and the map for Xu 781 contained 90 linkage groups with 1954 markers (1824 AFLP and 130 SSR) and covered 8,151.7 cM with an average marker distance of 4.2 cM. The maps described herein have the best coverage of the sweetpotato genome and the highest marker density reported to date. These are the first maps developed that have 90 complete linkage groups, which is in agreement with the actual number of chromosomes. Duplex and triplex markers were used to detect the homologous groups, and 13 and 14 homologous groups were identified in Xushu 18 and Xu 781 maps, respectively. Interval mapping was performed first and, subsequently, a multiple QTL model was used to refine the position and magnitude of the QTL. A total of 27 QTL for dry-matter content were mapped, explaining 9.0–45.1 % of the variation; 77.8 % of the QTL had a positive effect on the variation. This work represents an important step forward in genomics and marker-assisted breeding of sweetpotato.     

Application of AFLP, RAPD and ISSR markers to genetic mapping of European and Japanese larch

Genetic linkage maps have been increasingly developed for a wide variety of plants, using segregating populations such as F₂s or backcrosses between inbred lines. These pedigrees are rarely available in outbred species like forest trees which have long generation times. Thus genetic mapping studies have to use peculiar pedigrees and markers in appropriate configurations. We constructed single-tree genetic linkage maps of European larch (Larix decidua Mill.) and Japanese larch [Larix kaempferi (Lamb.) Carr.] using segregation data from 112 progeny individuals of an hybrid family. A total of 266 markers (114 AFLP, 149 RAPD and 3 ISSR loci) showing a testcross configuration, i.e.heterozygous in one parent and null in the other parent, were grouped at LOD 4.0, θ=0.3. The maternal parent map (L. decidua)consisted of 117 markers partitioned within 17 linkage groups (1152 cM) and the paternal parent map (L. kaempferi) had 125 markers assembled into 21 linkage groups (1206 cM). The map distance covered by markers was determined by adding a 34.7-cM independence distance at the end of each group and unlinked marker. It reached 2537 cM and 2997 cM respectively for European larch and Japanese larch, and represented respectively a 79.6% and 80.8% coverage of the overall genome. A few 3:1 segregating markers were used to identify homologous linkage groups between the European larch and the Japanese larch genetic maps. The PCR-based molecular markers allowed the construction of genetic maps, thus ensuring a good coverage of the larch genome for further QTL detection and mapping studies. Received: 15 March 1999 / Accepted: 29 March 1999     

Detection and validation of EST-SSR markers associated with sugar-related traits in sugarcane using linkage and association mapping

Sugar-related traits are of great importance in sugarcane breeding. In the present study, quantitative trait loci (QTL) mapping validated with association mapping was used to identify expressed sequence tag-simple sequence repeats (EST-SSRs) associated with sugar-related traits. For linkage mapping, 524 EST-SSRs, 241 Amplified Fragment Length Polymorphisms, and 10 genomic SSR markers were mapped using 283 F₁ progenies derived from an interspecific cross. Six regions were identified using Multiple QTL Mapping, and 14 unlinked markers using single marker analysis. Association analysis was performed on a set of 200 accessions, based on the mixed linear model. Validation of the EST-SSR markers using association mapping within the target QTL genomic regions identified two EST-SSR markers showing a putative relationship with uridine diphosphate (UDP) glycosyltransferase, and beta-amylase, which are associated with pol and sugar yield. These functional markers can be used for marker-assisted selection of sugarcane.     

基于SRAP分子标记的冰草遗传连锁图谱构建

构建高密度遗传连锁图谱是冰草抗性、品质、产量等重要性状QTL精细定位及标记辅助育种研究的基础。该试验以四倍体杂交冰草F2群体的202个分离单株及其亲本为材料,利用SRAP分子标记技术和Join Map 4.0作图软件对冰草的遗传连锁图谱进行了构建。结果表明:(1)共筛选出22对多态性好、标记位点清晰稳定的SRAP适宜引物,对冰草杂种F2分离单株的基因组DNA进行PCR扩增,共获得510个SRAP多态性标记位点,其比率占88.2%。(2)偏分离分析表明,偏分离标记比率仅为14.12%,符合遗传作图的要求。(3)成功构建了冰草的SRAP分子标记遗传连锁图谱,该图谱有14个连锁群、510个标记,连锁群间长度范围86.4～179.0cM,覆盖基因组总长度1 912.9cM,标记间平均间距3.75cM,为高密度遗传图谱。