Marker-Assisted Introgression in Backcross Breeding Programs

The efficiency of marker-assisted introgression in backcross populations derived from inbred lines was investigated by simulation. Background genotypes were simulated assuming that a genetic model of many genes of small effects in coupling phase explains the observed breed difference and variance in backcross populations. Markers were efficient in introgression backcross programs for simultaneously introgressing an allele and selecting for the desired genomic background. Using a marker spacing of 10-20 cM gave an advantage of one to two backcross generations selection relative to random or phenotypic selection. When the position of the gene to be introgressed is uncertain, for example because its position was estimated from a trait gene mapping experiment, a chromosome segment should be introgressed that is likely to include the allele of interest. Even for relatively precisely mapped quantitative trait loci, flanking markers or marker haplotypes should cover ~10-20 cM around the estimated position of the gene, to ensure that the allele frequency does not decline in later backcross generations.     

Optimal positioning of markers to control genetic background in marker-assisted backcrossing

Molecular markers are commonly used in backcross breeding programs in plants. As genetic maps contain more and more markers, it is of interest to determine which markers are to be used for selection. Here we describe how one can compute an optimal positioning of markers resulting in a maximization of the expected proportion of recipient genome. This criterion allows us to take selection into account and to produce relevant results regarding the final efficiency of background selection in backcross programs.     

Marker-assisted introgression using non-unique marker alleles II: selection on probability of presence of the introgressed allele

If marker alleles that identify a gene for introgression are not completely unique to the different base populations, the trait allele can be lost quickly during the process of backcrossing. This study considers ways to deal with incompletely informative markers in order to retain the desired allele. Selection was based on the probability of the presence of the desired (introgressed) trait allele, which was calculated for each marker genotype, using a single marker or a diallelic or triallelic marker bracket. The percentage of individuals retaining the introgressed allele was calculated over five generations of backcrossing, for selected fractions between 0 and 1, for marker alleles that could occur in both base populations. The best results were obtained with a rather large selected fraction, when all individuals, heterozygous and homozygous for the most desirable allele at the marker loci, were selected. Additional selection against marker homozygotes (which might have the highest probability of carrying the desired-trait allele, but produce uninformative gametes) altered the optimum selected fraction, making the selected fraction more consistently inversely related to a better retention of the desired-trait allele. A marker bracket was found to give a better retention of the desired-trait allele than a single marker and triallelic markers were better than diallelic markers, giving a retention of almost 50%. The earlier that preselection of parents (on informativeness) took place the better the overall result; preselection should occur preferably in the base populations. Preselection could make marker alleles unique to alternative base populations and markers would effectively become fully informative. Selection in the base populations might not be possible or not desirable, for example, because of the available number of individuals. This is unlikely to be a problem when parents are paired up to exclude any common marker alleles.     

Selection theory for marker-assisted backcrossing

Marker-assisted backcrossing is routinely applied in breeding programs for gene introgression. While selection theory is the most important tool for the design of breeding programs for improvement of quantitative characters, no general selection theory is available for marker-assisted backcrossing. In this treatise, we develop a theory for marker-assisted selection for the proportion of the genome originating from the recurrent parent in a backcross program, carried out after preselection for the target gene(s). Our objectives were to (i) predict response to selection and (ii) give criteria for selecting the most promising backcross individuals for further backcrossing or selfing. Prediction of response to selection is based on the marker linkage map and the marker genotype of the parent(s) of the backcross population. In comparison to standard normal distribution selection theory, the main advantage of our approach is that it considers the reduction of the variance in the donor genome proportion due to selection. The developed selection criteria take into account the marker genotype of the candidates and consider whether these will be used for selfing or backcrossing. Prediction of response to selection is illustrated for model genomes of maize and sugar beet. Selection of promising individuals is illustrated with experimental data from sugar beet. The presented approach can assist geneticists and breeders in the efficient design of gene introgression programs.     

Validation of growth-related quantitative trait loci markers in different Exopalaemon carinicauda families for marker-assisted selection

We detected growth-related QTL and associated markers from the backcross population of Exopalaemon carinicauda in the previous study. Based on our previous study, the 47 SNP markers associated with candidate growth trait QTL were selected to analyze the association between these markers and body weight (BW), body length and abdominal segment length traits in four different populations including wild population, a full-sib family, a half-sib family and a backcross population for evaluating their potential application of marker-assisted selection in E. carinicauda. The general linear model (GLM) and mixed linear model were applied and the associations between SNP loci and three growth-related traits verified. The results showed that the Marker79268 and Marker100644 were significantly associated with the BW trait in more than three populations by the GLM method. The Marker100644 was significantly associated with BW in the full-sib family, half-sib family and backcross populations by the GLM and mixed linear model methods. Our findings will provide useful SNP markers to go forward to improve growth performance through more refined marker-assisted selection in E. carinicauda.     

Marker pair selection for mapping quantitative trait loci

Mapping of quantitative trait loci (QTL) for backcross and F(2) populations may be set up as a multiple linear regression problem, where marker types are the regressor variables. It has been shown previously that flanking markers absorb all information on isolated QTL. Therefore, selection of pairs of markers flanking QTL is useful as a direct approach to QTL detection. Alternatively, selected pairs of flanking markers can be used as cofactors in composite interval mapping (CIM). Overfitting is a serious problem, especially if the number of regressor variables is large. We suggest a procedure denoted as marker pair selection (MPS) that uses model selection criteria for multiple linear regression. Markers enter the model in pairs, which reduces the number of models to be considered, thus alleviating the problem of overfitting and increasing the chances of detecting QTL. MPS entails an exhaustive search per chromosome to maximize the chance of finding the best-fitting models. A simulation study is conducted to study the merits of different model selection criteria for MPS. On the basis of our results, we recommend the Schwarz Bayesian criterion (SBC) for use in practice.     

Approximate Thresholds of Interval Mapping Tests for Qtl Detection

A general method is proposed for calculating approximate thresholds of interval mapping tests for quantitative trait loci (QTL) detection. Simulation results show that this method, when applied to backcross and F(2) populations, gives good approximations and is useful for any situation. Programs which calculate these thresholds for backcross, recombinant inbreds and F(2) for any given level and any chromosome with any given distribution of codominant markers were written in Fortran 77 and are available under request. The approach presented here could be used to obtain, after suitable calculations, thresholds for most segregating populations used in QTL mapping experiments.     

Size of donor chromosome segments around introgressed loci and reduction of linkage drag in marker-assisted backcross programs.

This article investigates the efficiency of marker-assisted selection in reducing the length of the donor chromosome segment retained around a locus held heterozygous by backcrossing. First, the efficiency of marker-assisted selection is evaluated from the length of the donor segment in backcrossed individuals that are (double) recombinants for two markers flanking the introgressed gene on each side. Analytical expressions for the probability density function, the mean, and the variance of this length are given for any number of backcross generations, as well as numerical applications. For a given marker distance, the number of backcross generations performed has little impact on the reduction of donor segment length, except for distant markers. In practical situations, the most important parameter is the distance between the introgressed gene and the flanking markers, which should be chosen to be as closely linked as possible to the introgressed gene. Second, the minimal population sizes required to obtain double recombinants for such closely linked markers are computed and optimized in the context of a multigeneration backcross program. The results indicate that it is generally more profitable to allow for three or more successive backcross generations rather than to favor recombinations in early generations.     

The construction of a substitution library of recombinant backcross lines in Brassica oleracea for the precision mapping of quantitative trait loci.

The currently available methods for locating quantitative trait loci (QTLs) and measuring their effects in segregating populations lack precision unless individual QTLs have very high heritabilities. The use of recombinant backcross lines containing short regions of donor chromosome introgressed into a constant recipient background permits QTLs to be located with greater precision. The present paper describes the use of molecular markers to introgress defined short regions of chromosome from a donor doubled haploid calabrese line of Brassica oleracea (var. italica) into a recipient short generation variety (Brassica oleracea var. alboglabra). We demonstrate that in just two or three generations of backcrossing, combined with selection for mapped molecular markers, the generation of a library of recombinant backcross lines is feasible. The possible use and refinement of these lines are discussed. Key words : backcrossing, Brassica oleracea, introgression, molecular markers, near-isogenic lines, QTL mapping, recombinant backcross lines, substitution lines.     

Genotype Selection to Rapidly Breed Congenic Strains

Congenic strains can now be constructed guided by the transmission of DNA markers. This allows not only selection for transmission of a desired, donor-derived differential region but also selection against the transmission of unwanted donor origin genomic material. The additional selection capacity should allow congenic strains to be produced in fewer generations than is possible with random backcrosses. Here, we consider modifications of a standard backcross breeding scheme to produce congenic mice by the inclusion of genotype-based selective breeding strategies. Simulation is used to evaluate the consequences of each strategy on the number of chromosomes that contain unwanted, donor-derived genetic material and the average length of this unwanted donor DNA for each backcross generation. Our prototypic strategy was to choose a single mouse to sire each generation using criteria designed to select against the transmission of chromosomes, other than the one containing the replacement genomic region, that contain any donor origin sequence at all. This chromosome elimination strategy resulted in an average of 16.4 chromosomes free of donor DNA in mice of the third backcross (N(3)) generation. A strategy based solely on positive selection for the replacement region required six backcross generations to achieve the same results.     

Advanced backcross QTL analysis: a method for the simultaneous discovery and transfer of valuable QTLs from unadapted germplasm into elite breeding lines

Advanced backcross QTL analysis is proposed as a method of combining QTL analysis with variety development. It is tailored for the discovery and transfer of valuable QTL alleles from unadapted donor lines (e.g., land races, wild species) into established elite inbred lines. Following this strategy, QTL analysis is delayed until the BC2 or BC3 generation and, during the development of these populations, negative selection is exercised to reduce the frequency of deleterious donor alleles. Simulations suggest that advanced backcross QTL analysis will be effective in detecting additive, dominant, partially dominant, or overdominant QTLs. Epistatic QTLs or QTLs with gene actions ranging from recessive to additive will be detected with less power than in selfing generations. QTL-NILs can be derived from advanced backcross populations in one or two additional generations and utilized to verify QTL activity. These same QTL-NILs also represent commercial inbreds improved (over the original recurrent inbred line) for one or more quantitative traits. The time lapse from QTL discovery to construction and testing of improved QTL-NILs is minimal (1–2 years). If successfully employed, advanced backcross QTL analysis can open the door to exploiting unadapted and exotic germplasm for the quantitative trait improvement of a number of crop plants.     

Preselection statistics and Random Forest classification identify population informative single nucleotide polymorphisms in cosmopolitan and autochthonous cattle breeds

Commercial single nucleotide polymorphism (SNP) arrays have been recently developed for several species and can be used to identify informative markers to differentiate breeds or populations for several downstream applications. To identify the most discriminating genetic markers among thousands of genotyped SNPs, a few statistical approaches have been proposed. In this work, we compared several methods of SNPs preselection (Delta, F_st and principal component analyses (PCA)) in addition to Random Forest classifications to analyse SNP data from six dairy cattle breeds, including cosmopolitan (Holstein, Brown and Simmental) and autochthonous Italian breeds raised in two different regions and subjected to limited or no breeding programmes (Cinisara, Modicana, raised only in Sicily and Reggiana, raised only in Emilia Romagna). From these classifications, two panels of 96 and 48 SNPs that contain the most discriminant SNPs were created for each preselection method. These panels were evaluated in terms of the ability to discriminate as a whole and breed-by-breed, as well as linkage disequilibrium within each panel. The obtained results showed that for the 48-SNP panel, the error rate increased mainly for autochthonous breeds, probably as a consequence of their admixed origin lower selection pressure and by ascertaining bias in the construction of the SNP chip. The 96-SNP panels were generally more able to discriminate all breeds. The panel derived by PCA-chrom (obtained by a preselection chromosome by chromosome) could identify informative SNPs that were particularly useful for the assignment of minor breeds that reached the lowest value of Out Of Bag error even in the Cinisara, whose value was quite high in all other panels. Moreover, this panel contained also the lowest number of SNPs in linkage disequilibrium. Several selected SNPs are located nearby genes affecting breed-specific phenotypic traits (coat colour and stature) or associated with production traits. In general, our results demonstrated the usefulness of Random Forest in combination to other reduction techniques to identify population informative SNPs.     

Development of sequence characterized DNA markers linked to leaf rust (<Emphasis Type="Italic">Hemileia vastatrix</Emphasis>) resistance in coffee (<Emphasis Type="Italic">Coffea arabica</Emphasis> L.)

Coffee leaf rust due to Hemileia vastatrix is one of the most serious diseases in Arabica coffee (Coffea arabica). A resistance gene (S_H3) has been transferred from C. liberica into C. arabica. The present work aimed at developing sequence-characterized genetic markers for leaf rust resistance. Linkage between markers and leaf rust resistance was tested by analysing two segregating populations, one F₂ population of 101 individuals and one backcross (BC₂) population of 43 individuals, derived from a cross between a susceptible and a SH3-introgressed resistant genotype. A total of ten sequence-characterized genetic markers closely associated with the S_H3 leaf rust resistance gene were generated. These included simple sequence repeats (SSR) markers, sequence-characterised amplified regions (SCAR) markers resulting from the conversion of amplified fragment length polymorphism (AFLP) markers previously identified and SCAR markers derived from end-sequences of bacterial artificial chromosome (BAC) clones. Those BAC clones were identified by screening of C. arabica genomic BAC library using a cloned AFLP-marker as probe. The markers we developed are easy and inexpensive to run, requiring one PCR step followed by gel separation. While three markers were linked in repulsion with the S_H3 gene, seven markers were clustered in coupling around the S_H3 gene. Notably, two markers appeared to co-segregate perfectly with the S_H3 gene in the two plant populations analyzed. These markers are suitable for marker-assisted selection for leaf rust resistance and to facilitate pyramiding of the S_H3 gene with other leaf rust resistance genes.     

Marker-assisted introgression of 4 Phytophthora capsici resistance QTL alleles into a bell pepper line: validation of additive and epistatic effects

The aim of the present study was to transfer resistance to P. capsici alleles at four quantitative trait loci (QTLs) from a small fruited pepper into a bell pepper recipient line using markers. The marker-assisted selection program was initiated from a doubled-haploid line issued from the mapping population and involved three cycles of marker-assisted backcross (MAB). Two populations, derived by selfing the plants selected after the first selection cycle, were genotyped and evaluated phenotypically for their resistance level. The additive and epistatic effects of the four resistance factors were re-detected and validated in these populations, indicating that introgression of 4 QTLs in this MAB program was successful. A decrease of the effect for the moderate-effect QTLs and of the epistatic interaction was observed. Phenotypic evaluations of horticultural traits were performed on sample of each backcross generation. The results indicated an efficient return to the recipient phenotype using this MAB strategy.     

RFLP analysis of the size of chromosomal segments retained around the Tm-2 locus of tomato during backcross breeding

Summary Genes introduced into cultivated plants by backcross breeding programs are flanked by introgressed segments of DNA derived from the donor parent. This phenomenon is known as linkage drag and is frequently thought to affect traits other than the one originally targeted. The Tm-2 gene of Lycopersicon peruvianum, which confers resistance to tobacco mosaic virus, was introduced into several different tomato cultivars (L. esculentum) by repeated backcrossing. We have measured the sizes of the introgressed segments flanking the Tm-2 locus in several of these cultivars using a high density map of restriction fragment length polymorphic (RFLP) markers. The smallest introgressed segment is estimated to be 4 cM in length, while the longest is over 51 cM in length and contains the entire short arm of chromosome 9. Additionally, RFLP analysis was performed on remnant seed from different intermediate generations corresponding to two different backcross breeding programs for TMV resistance. The results reveal that plants containing desirable recombination near the resistance gene were rarely selected during backcrossing and, as a result, the backcross breeding method was largely ineffective in reducing the size of linked DNA around the resistance gene. We propose that, by monitoring recombination around genes of interest with linked RFLP markers, one can quickly and efficiently reduce the amount of linkage drag associated with introgression. Using such a procedure, it is estimated that an introgressed segment can be obtained in two generations that is as small as that which would otherwise require 100 backcross generations without RFLP selection.     

Unilateral incompatibility as a major cause of skewed segregation in the cross between Lycopersicon esculentum and L. pennellii

Skewed segregations are frequent events in segregating populations derived from different interspecific crosses in tomato. To determine a basis for skewed segregations in the progeny of the cross between Lycopersicon esculentum and L. pennellii, monogenic segregations of 16 isozyme loci were analyzed in an F₂ and two backcross populations of this cross. In the F₂, 9 loci mapping to chromosomes 1, 2, 4, 9, 10 and 12 exhibited skewed segregations and in all cases there was an excess of L. pennellii homozygotes. The genotypic frequencies at all but one locus were at Hardy-Weinberg equilibria. In the backcross populations, all except two loci exhibited normal Mendelian segregations. No post-zygotic selection model could statistically or biologically explain the observed segregation patterns in the F₂ and backcross populations. A pre-zygotic selection model, assuming selective elimination of the male gametophytes during pollen function (i.e., from pollination to karyogamy), could adequately explain the observed segregations in all three populations. The direction of the skewed segregations in the F₂ population was consistent with that expected based on the effects of unilateral incompatibility reactions between the two species. In addition, the chromosomal locations of 5 of the 9 markers that exhibited skewed segregations coincided with the locations of several known compatibility-related genes in tomato. Multigenic unilateral incompatibility reactions between L. esculentum pollen and the stigma or style of L. pennellii (or its hybrid derivatives) are suggested to be the major cause of the skewed segregations in the F₂ progeny of this cross.     

Marker-Assisted Introgression of Quantitative Trait Loci

The use of molecular markers for the introgression of one or several superior QTL alleles into a recipient line is investigated using analytic and simulation results. The positions of the markers devoted to the control of the genotype at the QTLs in a ``foreground selection' step are optimized given the confidence interval of the QTL position. Results demonstrate that using at least three markers per QTL allows a good control over several generations. Population sizes that should be recommended for various numbers of QTLs are calculated and are used to determine the limit in the number of QTLs that can be monitored simultaneously. If ``background selection' devoted to accelerate the return to the recipient parent genotype outside the QTL regions is applied, the positions of the markers devoted to the control of the QTLs have to be reconsidered. When several QTLs are monitored simultaneously, background selection among the limited number of individuals resulting from the foreground selection step accelerates the increase in genomic similarity with the recipient parent, with only limited costs. Background selection is even more efficient in a pyramidal backcross program where QTLs are first monitored one by one.     

单位点孢子体-配子体互作模型中不育基因的位置和效应的估计方法

水稻籼粳亚种间杂交F1通常表现为高度不育,这种不育性的一种遗传学解释称为单位点孢子体-配子体互作模型.为了研究这种不育性,提出了一种统计方法,可以估计单位点孢子体-配子体互作模型中不育基因位点的位置和效应.该方法利用回交群体中呈现异常分离的标记位点,用最大似然法对不育基因与标记位点之间的重组率和雌配子存活率进行估计.由于所依据的是非连续变异的遗传标记的分离,而不是连续分布的配子育性指标,因此可以避免由育性直接估计所带来的重组率结果的不稳定.     

A nearest-neighboring-end algorithm for genetic mapping

MOTIVATION: High-throughput methods are beginning to make possible the genotyping of thousands of loci in thousands of individuals, which could be useful for tightly associating phenotypes to candidate loci. Current mapping algorithms cannot handle so many data without building hierarchies of framework maps. RESULTS: A version of Kruskal's minimum spanning tree algorithm can solve any genetic mapping problem that can be stated as marker deletion from a set of linkage groups. These include backcross, recombinant inbred, haploid and double-cross recombinational populations, in addition to conventional deletion and radiation hybrid populations. The algorithm progressively joins linkage groups at increasing recombination fractions between terminal markers, and attempts to recognize and correct erroneous joins at peaks in recombination fraction. The algorithm is O (mn3) for m individuals and n markers, but the mean run time scales close to mn2. It is amenable to parallel processing and has recovered true map order in simulations of large backcross, recombinant inbred and deletion populations with up to 37,005 markers. Simulations were used to investigate map accuracy in response to population size, allelic dominance, segregation distortion, missing data and random typing errors. It produced accurate maps when marker distribution was sufficiently uniform, although segregation distortion could induce translocated marker orders. The algorithm was also used to map 1003 loci in the F7 ITMI population of bread wheat, Triticum aestivum L. emend Thell., where it shortened an existing standard map by 16%, but it failed to associate blocks of markers properly across gaps within linkage groups. This was because it depends upon the rankings of recombination fractions at individual markers, and is susceptible to sampling error, typing error and joint selection involving the terminal markers of nearly finished linkage groups. Therefore, the current form of the algorithm is useful mainly to improve local marker ordering in linkage groups obtained in other ways. AVAILABILITY: The source code and supplemental data are http://www.iubio.bio.indiana.edu/soft/molbio/qtl/flipper/ CONTACT: ccrane@purdue.edu.     

Mapping 28 erythrocyte antigen, plasma protein and enzyme polymorphisms using an efficient genomic scan of the porcine genome

One hundred and fifty-four microsatellite markers were selected for genomic scanning of the porcine genome and were grouped into amplification sets to reduce the cost and labour required. Thirty amplification sets had two markers (duplex), 20 sets had three markers (triplex) and five sets had four markers (quadruplex) while 14 markers were analysed separately. The selection criteria for microsatellites were: ease of scoring, level of polymorphism, genetic location and ability to be genotyped in a multiplexed polymerase chain reaction (PCR). The selected microsatellites were chosen to span the entire genome flanked by the porcine linkage map with intervals between adjacent markers of 15–20 cM where possible. The utility of this set of markers was demonstrated by linkage analyses with loci controlling blood plasma protein and red cell enzyme polymorphisms ( n = 13), erythrocyte antigens ( n = 15), the S blood group, coat colour and ryanodine receptor from 174 backcross Meishan-White Composite pigs. These loci displayed various forms of inheritance and most (24 loci) have been placed in linkage groups. Significant two-point linkages (lod > 3·0) were detected for each polymorphic marker. These results provide the first linkage assignments for phosphoglucomutase (PGM2) and erythrocyte antigen F (EAF) to SSC8; and serum amylase (AMY) and erythrocyte antigen I (EAI) to SSC18. All of the remaining polymorphic loci ( n = 24) mapped to previously identified regions confirming earlier results. Most of the markers used in this study should be useful in resource populations of various breed crosses as the number of alleles detected in a multibreed reference population was one of the selection criteria.