Marker-assisted introgression using non-unique marker alleles I: selection on the presence of linked marker alleles

This paper investigates marker-assisted introgression of a major gene into an outbred line, where identification of the introgressed gene is incomplete because marker alleles are not unique to the base populations (the same marker allele can occur in both donor and recipient population). Those markers are used to identify the introgressed allele as well as the background genotype. The effect of using those markers, as if they were completely informative on the retention of the introgressed allele, was examined over five generations of backcrossing by using a single marker or a marker bracket for different starting frequencies of the marker alleles. Results were calculated by using both a deterministic approach, where selection is only for the desired allele, and by a stochastic approach, where selection is also on background genotype. When marker allele frequencies in donor and recipient population diverged from 1 and 0 (using a diallelic marker), the ability to retain the desired allele rapidly declined. Marker brackets performed notably better than single markers. If selection on background marker genotype was applied, the desired allele could be lost even more quickly than expected at random because the chance that the allele, which is common in the donor line, is present on the locus identifying the introgressed allele and is surrounded by alleles common in the recipient line on the background marker loci, will descend from the donor line (double recombination has taken place), is a lot smaller than the chance that this allele will stem from the recipient line (in which the allele occurs in low frequency). Marker brackets again performed better. Preselection against marker homozygotes (producing uninformative gametes) gave a slightly better retention of the introgressed allele.     

Effectiveness of bulking procedures in measuring population-pairwise similarity with dominant and co-dominant genetic markers

Two bulking procedures (bulking individuals before and after genotyping) are commonly applied in similarity based studies of genetic distance at the population or higher level, but their effectiveness is largely unknown. In this study, expected population-pairwise similarity for both bulking procedures is derived with dominant and co-dominant diallelic markers. Numerical examples for the derived formulae are given with up to ten individuals randomly selected from each population. The procedure of bulking individuals after genotyping with either marker system is generally more informative than the procedure of bulking individuals before genotyping, because the former incorporates the information from marker alleles of intermediate frequency. Both procedures are effective with 5–10 individuals selected randomly from either population, but the procedure of bulking before genotyping requires a genotyping effort several-fold less than the procedure of bulking after genotyping. For either bulking procedure, a co-dominant marker system is generally more informative than a dominant marker system. Received: 20 October 1999 / Accepted: 11 November 1999     

Use of near-isogenic lines derived by backcrossing or selfing to map qualitative traits

Near-isogenic lines (NILs) are a valuable resource for detecting linkages between qualitative trait loci and molecular markers. Molecular marker studies are expensive and methods that require genotyping fewer individuals, such as the NIL-analysis method, are desirable. We present a theory for using sets of NILs to detect linkages between molecular markers and introgressed loci. The probability that a marker a specific distance from the introgressed gene will have a donor parent allele in a near-isogenic line is a function of the distance between the marker and the gene, and the number of back-crosses and/or selfs used in deriving the NIL. The binomial probability formula is used to calculate the probability of having a donor parent allele at a given marker when sets of NILs are used. The formulae given allow calculation of the probability that a marker is linked to the introgressed gene, as well as the probability that a gene will be successfully detected when using given numbers of NILs, backcrosses, and molecular markers.     

Unequivocal determination of sire allele origin for multiallelic microsatellites when only the sire and progeny are genotyped

The effect of a segregating economic trait locus (ETL) can be detected with the aid of a linked genetic marker, if specific alleles of each locus are in association among the individuals genotyped for the genetic marker. For dairy cattle this can be achieved by application of the ‘granddaughter design’. If only the sires and their sons are genotyped for the genetic markers, then the allele origin of sons having the same genotypes as their sires cannot be determined. Seven sires and 101 sons were genotyped for five microsatellites. The mean frequency of heterozygous sires was 77%. The mean number of alleles per locus was 8.2. Frequency of informative sons per locus ranged from 60% to 80% with a mean of 72%. With highly polymorphic microsatellites, at least 60% more grandsire families can be included in the analysis, and the number of sons assayed can be reduced by 40%, as compared to diallelic markers.     

Efficient marker-based recurrent selection for multiple quantitative trait loci

We studied the efficiency of recurrent selection based solely on marker genotypes (marker-based selection), in order to increase favourable allele frequency at 50 previously detected quantitative trait loci (QTLs). Two selection procedures were investigated, using computer simulations: (1) Truncation Selection (MTS), in which individuals are ranked based on marker score, and best individuals are selected for recombination; and (2) QTL Complementation Selection (QCS), in which individuals are selected such that their QTL composition complements those individuals already selected. Provided QTL locations are accurate, marker-based selection with a population size of 200 was very effective in rapidly increasing frequencies of favourable QTL alleles. QCS methods were more effective than MTS for improving the mean frequency and fixation of favourable QTL alleles. Marker-based selection was not very sensitive to a reduction in population size, and appears valuable to optimize the use of molecular markers in recurrent selection programmes.     

Selection in backcross programmes

Backcrossing is a well-known and long established breeding scheme where a characteristic is introgressed from a donor parent into the genomic background of a recurrent parent. The various uses of backcrossing in modern genetics, particularly with the help of molecular markers, are reviewed here. Selection in backcross programmes is used to either improve the genetic value of plant and animal populations or fine map quantitative trait loci. Both cases are helpful in our understanding of the genetic bases of quantitative traits variation.     

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.     

Mapping genes that underlie ethnic differences in disease risk: methods for detecting linkage in admixed populations, by conditioning on parental admixture.

Genes that underlie ethnic differences in disease risk can be mapped in affected individuals of mixed descent if the ancestry of the alleles at each marker locus can be assigned to one of the two founding populations. Linkage can be detected by testing for association of the disease with the ancestry of alleles at the marker locus, by conditioning on the admixture (defined as the proportion of genes that have ancestry from the high-risk population) of both parents. With regard to exploiting the effects of admixture, this test is more flexible and powerful than the transmission-disequilibrium test. Under the assumption of a multiplicative model, the statistical power for a given sample size depends only on parental admixture and the risk ratio r between populations that is generated by the locus. The most informative families are those in which mean parental admixture is .2-.7 and in which admixture is similar in both parents. The number of markers required for a genome search depends on the number of generations since admixture and on the information content for ancestry (f) of the markers, defined as a function of allele frequencies in the two founding populations. Simulations using a hidden Markov model suggest that, when admixture has occurred 2-10 generations earlier, a multipoint analysis using 2,000 biallelic markers, with f values of 30%, can extract 70%-90% of the ancestry information for each locus. Sets of such markers could be selected from libraries of single-nucleotide polymorphisms, when these become available.     

The length of the intact donor chromosome segment around a target gene in marker-assisted backcrossing

Recurrent backcrossing is an established procedure to transfer target genes from a donor into the genetic background of a recipient genotype. By assessing the parental origin of alleles at markers flanking the target locus one can select individuals with a short intact donor chromosome segment around the target gene and thus reduce the linkage drag. We investigated the probability distribution of the length of the intact donor chromosome segment around the target gene in recurrent backcrossing with selection for heterozygosity at the target locus and homozygosity for the recurrent parent allele at flanking markers for a diploid species. Assuming no interference in crossover formation, we derived the cumulative density function, probability density function, expected value, and variance of the length of the intact chromosome segment for the following cases: (1) backcross generations prior to detection of a recombinant individual between the target gene and the flanking marker; (2) the backcross generation in which for the first time a recombinant individual is detected, which is selected for further backcrossing; and (3) subsequent backcross generations after selection of a recombinant. Examples are given of how these results can be applied to investigate the efficiency of marker-assisted backcrossing for reducing the length of the intact donor chromosome segment around the target gene under various situations relevant in breeding and genetic research.     

Genetic structure of five susceptibility gene regions for coronary artery disease: disequilibria within and among regions

We analyzed two-locus disequilibria for 16 polymorphic loci of seven susceptibility genes for coronary artery disease located in five chromosomal regions distributed across four chromosomes. Included were the genes coding for apolipoprotein B (ApoB, chromosome 2, four marker loci), lipoprotein lipase (LPL, chromosome 8, three marker loci), apolipoproteins AI, CIII, AIV (ApoAI–CIII–AIV, chromosome 11, three marker loci), apolipoprotein E (ApoE, chromosome 19, two marker loci), and the low density lipoprotein receptor (LDLR, chromosome 19, four marker loci). Our sample included 540 unrelated individuals from the Rochester, Minn. population. There were no statistically significant deviations of single-locus genotypes from Hardy-Weinberg equilibrium. The strongest associations within genes were for composite diallelic disequilibria; 17/19 were significant (13 at Pr <0.001, 1 at Pr <0.01, 3 at Pr <0.05). These observations suggest marker alleles within genes have a shared evolutionary history reflected by disequilibria that have not been dissipated by recombination. Disequilibrium was not generally concordant with the physical orderings of markers. Only two significant higher-order disequilibria were observed although 12 triallelic disequilibria were at maximum possible values. We observed 19 statistically significant disequilibria (Pr <0.05; 4 composite diallelic, 13 triallelic, and 2 quadriallelic) between 101 pairs of marker loci, where each locus in a pair was from a different unlinked region. These unexpected results are most likely explained by recent historical factors, including worldwide population expansion and amalgamation with continuous admixture, that influence the genetic structure (organization of alleles and non-alleles into genotypes) of a population. We conclude that disequilibria between loci from unlinked regions may be more extensive than is commonly assumed. Our findings also suggest that it is, on average, at least 15 times more likely to not detect significant disequilibrium among unlinked loci when it is really present than to make a false positive inference. Disequilibria between functional loci within or between regions will impact estimates of genetic variance associated with particular functional mutations within a susceptibility gene region. Received: 15 January 1998 / Accepted: 24 June 1998     

Segregation and conditional probability association of molecular markers with traits in autotetraploid alfalfa

The complications introduced by the autotetraploid, outcrossing nature of alfalfa (Medicago sativa L.) as related to detecting associations of marker loci and traits of interest are discussed, and a new method of detecting marker-trait associations is suggested. This method utilizes plant populations that are likely to have been produced through the plant breeding process: populations selected for one trait, and the base, unselected population. Marker allele frequency shifts between the populations are indicative of genomic regions involved in trait expression, and may indicate alleles that have reached the triplex or homozygous state and do not segregate in S₁ or F₁ populations. However, because many, perhaps hundreds, of sequential frequency comparisons are needed to detect fragments in significantly different proportions in the two populations, the type I error rate is very high. A resampling-based analysis method is proposed to address the concern of the type I error rate, and identify marker alleles associated with this trait of interest. The utility of marker-trait associations thus defined for identifying individual plants from heterogeneous populations was investigated through model-building and conditional probability studies. Factors investigated that influenced the utility of the marker associations and (in the base population) the frequencies of the trait and marker, and the frequencies of the markers in plants exhibiting the trait and in the plants not exhibiting the trait. The frequency of occurrence of a marker in undesirable plants profoundly influenced the efficiency with which the marker could be used to select desirable plants, however, under some circumstances, markers or combinations of markers can be highly efficient for selecting rare, desirable plants from a heterogeneous base population.     

Genome-enabled hitchhiking mapping identifies QTLs for stress resistance in natural Drosophila

Identification of genes underlying complex traits is an important problem. Quantitative trait loci (QTL) are mapped using marker-trait co-segregation in large panels of recombinant genotypes. Most frequently, recombinant inbred lines derived from two isogenic parents are used. Segregation patterns are also studied in pedigrees from multiple families. Great advances have been made through creative use of these techniques, but narrow sampling and inadequate power represent strong limitations. Here, we propose an approach combining the strengths of both techniques. We established a mapping population from a sample of natural genotypes, and applied artificial selection for a complex character. Selection changed the frequencies of alleles in QTLs contributing to the selection response. We infer QTLs with dense genotyping microarrays by identifying blocks of linked markers undergoing selective changes in allele frequency. We demonstrated this approach with an experimental population composed from 20 isogenic strains. Selection for starvation survival was executed in three replicated populations with three control non-selected populations. Three individuals per population were genotyped using Affymetrix GeneChips. Two regions of the genome, one each on the left arms of the second and third chromosomes, showed significant divergence between control and selected populations. For the former region, we inferred allele frequencies in selected and control populations by pyrosequencing. We conclude that the allele frequency difference, averaging approximately 40% between selected and control lines, contributed to selection response. Our approach can contribute to the fine scale decomposition of the genetics of direct and indirect selection responses, and genotype by environment interactions.     

Assessment of genetic variation and differentiation of hop genotypes by microsatellite and AFLP markers.

Microsatellites have many desirable marker properties and have been increasingly used in crop plants in genetic diversity studies. Here we report on the characterisation of microsatellite markers and on their use for the determination of genetic identities and the assessment of genetic variability among accessions from a germplasm collection of hop. Thirty-two polymorphic alleles were found in the 55 diploid genotypes, with an average number of eight alleles (3.4 effective alleles) for four microsatellite loci. Calculated polymorphic information content values classified three loci as informative markers and two loci as suitable for mapping. The average observed heterozygosity was 0.7 and the common probability of identical genotypes was 3.271 x 10(-4). An additional locus, amplified by one primer pair, was confirmed by segregation analysis of two crosses. The locus discovered was heterozygous, with a null allele in the segregating population. The same range of alleles was detected in nine triploid and five tetraploid hop genotypes. Cultivar heterozygosity varied among all 69 accessions, with only one cultivar being homozygous at four loci. Microsatellite allele polymorphisms distinguished 81% of all genotypes; the same allelic profile was found mainly in clonally selected cultivars. Cultivar-specific alleles were found in some genotypes, as well as a specific distribution of alleles in geographically distinct hop germplasms. The genetic relationship among 41 hop accessions was compared on the basis of microsatellite and AFLP polymorphisms. Genetic similarity dendrograms showed low correlation between the two marker systems. The microsatellite dendrogram grouped genetically related accessions reasonably well, while the AFLP dendrogram showed good clustering of closely related accessions and, additionally, separated two geographically distinct hop germplasms. The results of microsatellite and AFLP analysis are discussed from the point of view of the applicability of the two marker systems for different aspects of germplasm evaluation.     

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.     

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.     

A powerful likelihood method for the analysis of linkage disequilibrium between trait loci and one or more polymorphic marker loci.

Historically, most methods for detecting linkage disequilibrium were designed for use with diallelic marker loci, for which the analysis is straightforward. With the advent of polymorphic markers with many alleles, the normal approach to their analysis has been either to extend the methodology for two-allele systems (leading to an increase in df and to a corresponding loss of power) or to select the allele believed to be associated and then collapse the other alleles, reducing, in a biased way, the locus to a diallelic system. I propose a likelihood-based approach to testing for linkage disequilibrium, an approach that becomes more conservative as the number of alleles increases, and as the number of markers considered jointly increases in a multipoint test for linkage disequilibrium, while maintaining high power. Properties of this method for detecting associations and fine mapping the location of disease traits are investigated. It is found to be, in general, more powerful than conventional methods, and it provides a tractable framework for the fine mapping of new disease loci. Application to the cystic fibrosis data of Kerem et al, is included to illustrate the method.     

Multipoint linkage analysis with many multiallelic or dense diallelic markers: Markov chain-Monte Carlo provides practical approaches for genome scans on general pedigrees

Computations for genome scans need to adapt to the increasing use of dense diallelic markers as well as of full-chromosome multipoint linkage analysis with either diallelic or multiallelic markers. Whereas suitable exact-computation tools are available for use with small pedigrees, equivalent exact computation for larger pedigrees remains infeasible. Markov chain-Monte Carlo (MCMC)-based methods currently provide the only computationally practical option. To date, no systematic comparison of the performance of MCMC-based programs is available, nor have these programs been systematically evaluated for use with dense diallelic markers. Using simulated data, we evaluate the performance of two MCMC-based linkage-analysis programs--lm_markers from the MORGAN package and SimWalk2--under a variety of analysis conditions. Pedigrees consisted of 14, 52, or 98 individuals in 3, 5, or 6 generations, respectively, with increasing amounts of missing data in larger pedigrees. One hundred replicates of markers and trait data were simulated on a 100-cM chromosome, with up to 10 multiallelic and up to 200 diallelic markers used simultaneously for computation of multipoint LOD scores. Exact computation was available for comparison in most situations, and comparison with a perfectly informative marker or interprogram comparison was available in the remaining situations. Our results confirm the accuracy of both programs in multipoint analysis with multiallelic markers on pedigrees of varied sizes and missing-data patterns, but there are some computational differences. In contrast, for large numbers of dense diallelic markers, only the lm_markers program was able to provide accurate results within a computationally practical time. Thus, programs in the MORGAN package are the first available to provide a computationally practical option for accurate linkage analyses in genome scans with both large numbers of diallelic markers and large pedigrees.     

Mapping quantitative trait loci with DNA microsatellites in a commercial dairy cattle population

Individual loci affecting economically important traits can be located using genetic linkage between quantitative trait loci and genetic markers. In the ‘granddaughter’ experimental design, heterozygous grandsires and their sons are genotyped for the genetic marker, while the quantitative trait records of the granddaughters are used for statistical analysis. Ten DNA microsatellite markers were used to look for associations with quantitative trait loci affecting milk production traits in seven Israeli Holstein grandsire families. At least 60% more grandsires were heterozygous, and 40% fewer individuals were discarded because of unknown paternal allele origin, as compared with diallelic markers. The effects of paternal alleles for locus D21S4 on kg milk and protein were significant (P < 0.025). The allele substitution effects for sire 783 were 283 kg milk and 5.7 kg protein. For both traits, progeny of sire 783 that inherited allele ‘18’ had higher evaluations than progeny that inherited allele ‘21’. These results were verified by genotyping 151 of his daughters. Thus, the rate of genetic gain for protein production can be increased by selecting progeny of sire 783 carrying allele ‘18’ at this locus.     

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.     

Developing rice with high yield under phosphorus deficiency: Pup1 sequence to application

The major quantitative trait locus (QTL) Phosphorus uptake1 (Pup1) confers tolerance of phosphorus deficiency in soil and is currently one of the most promising QTLs for the development of tolerant rice (Oryza sativa) varieties. To facilitate targeted introgression of Pup1 into intolerant varieties, the gene models predicted in the Pup1 region in the donor variety Kasalath were used to develop gene-based molecular markers that are evenly distributed over the fine-mapped 278-kb QTL region. To validate the gene models and optimize the markers, gene expression analyses and partial allelic sequencing were conducted. The markers were tested in more than 80 diverse rice accessions revealing three main groups with different Pup1 allele constitution. Accessions with tolerant (group I) and intolerant (group III) Pup1 alleles were distinguished from genotypes with Kasalath alleles at some of the analyzed loci (partial Pup1; group II). A germplasm survey additionally confirmed earlier data showing that Pup1 is largely absent from irrigated rice varieties but conserved in varieties and breeding lines adapted to drought-prone environments. A core set of Pup1 markers has been defined, and sequence polymorphisms suitable for single-nucleotide polymorphism marker development for high-throughput genotyping were identified. Following a marker-assisted backcrossing approach, Pup1 was introgressed into two irrigated rice varieties and three Indonesian upland varieties. First phenotypic evaluations of the introgression lines suggest that Pup1 is effective in different genetic backgrounds and environments and that it has the potential to significantly enhance grain yield under field conditions.