Mapping quantitative trait loci with dominant markers in four-way crosses

 A common problem in mapping quantitative trait loci (QTLs) is that marker data are often incomplete. This includes missing data, dominant markers, and partially informative markers, arising in outbred populations. Here we briefly present an iteratively re-weighted least square method (IRWLS) to incorporate dominant and missing markers for mapping QTLs in four-way crosses under a heterogeneous variance model. The algorithm uses information from all markers in a linkage group to infer the QTL genotype. Monte Carlo simulations indicate that with half dominant markers, QTL detection is almost as efficient as with all co-dominant markers. However, the precision of the estimated QTL parameters generally decreases as more markers become missing or dominant. Notable differences are observed on the standard deviation of the estimated QTL position for varying levels of marker information content. The method is relatively simple so that more complex models including multiple QTLs or fixed effects can be fitted. Finally, the method can be readily extended to QTL mapping in full-sib families. Received: 16 June 1998 / Accepted: 29 September 1998     

Fine mapping of complex trait genes combining pedigree and linkage disequilibrium information: a Bayesian unified framework

We present a Bayesian method that combines linkage and linkage disequilibrium (LDL) information for quantitative trait locus (QTL) mapping. This method uses jointly all marker information (haplotypes) and all available pedigree information; i.e., it is not restricted to any specific experimental design and it is not required that phases are known. Infinitesimal genetic effects or environmental noise ("fixed") effects can equally be fitted. A diallelic QTL is assumed and both additive and dominant effects can be estimated. We have implemented a combined Gibbs/Metropolis-Hastings sampling to obtain the marginal posterior distributions of the parameters of interest. We have also implemented a Bayesian variant of usual disequilibrium measures like D' and r(2) between QTL and markers. We illustrate the method with simulated data in "simple" (two-generation full-sib families) and "complex" (four-generation) pedigrees. We compared the estimates with and without using linkage disequilibrium information. In general, using LDL resulted in estimates of QTL position that were much better than linkage-only estimates when there was complete disequilibrium between the mutant QTL allele and the marker. This advantage, however, decreased when the association was only partial. In all cases, additive and dominant effects were estimated accurately either with or without disequilibrium information.     

A linkage map of Atlantic salmon (Salmo salar) reveals an uncommonly large difference in recombination rate between the sexes

A genetic linkage map of the Atlantic salmon (Salmo salar) was constructed, using 54 microsatellites and 473 amplified fragment length polymorphism (AFLP) markers. The mapping population consisted of two full-sib families within one paternal half-sib family from the Norwegian breeding population. A mapping strategy was developed that facilitated the construction of separate male and female maps, while retaining all the information contributed by the dominant AFLP markers. By using this strategy, we were able to map a significant number of the AFLP markers for which all informative offspring had two heterozygous parents; these markers then served as bridges between the male and female maps. The female map spanned 901 cM and had 33 linkage groups, while the male spanned 103 cM and had 31 linkage groups. Twenty-five linkage groups were common between the two maps. The construction of the genetic map revealed a large difference in recombination rate between females and males. The ratio of female recombination rate vs. male recombination rate was 8.26, the highest ratio reported for any vertebrate. This map constitutes the first linkage map of Atlantic salmon, one of the most important aquaculture species worldwide.     

Construction of a genetic linkage map in tetraploid species using molecular markers

This article presents methodology for the construction of a linkage map in an autotetraploid species, using either codominant or dominant molecular markers scored on two parents and their full-sib progeny. The steps of the analysis are as follows: identification of parental genotypes from the parental and offspring phenotypes; testing for independent segregation of markers; partition of markers into linkage groups using cluster analysis; maximum-likelihood estimation of the phase, recombination frequency, and LOD score for all pairs of markers in the same linkage group using the EM algorithm; ordering the markers and estimating distances between them; and reconstructing their linkage phases. The information from different marker configurations about the recombination frequency is examined and found to vary considerably, depending on the number of different alleles, the number of alleles shared by the parents, and the phase of the markers. The methods are applied to a simulated data set and to a small set of SSR and AFLP markers scored in a full-sib population of tetraploid potato.     

A comparison of single‐sample estimators of effective population sizes from genetic marker data

In molecular ecology and conservation genetics studies, the important parameter of effective population size (N_e) is increasingly estimated from a single sample of individuals taken at random from a population and genotyped at a number of marker loci. Several estimators are developed, based on the information of linkage disequilibrium (LD), heterozygote excess (HE), molecular coancestry (MC) and sibship frequency (SF) in marker data. The most popular is the LD estimator, because it is more accurate than HE and MC estimators and is simpler to calculate than SF estimator. However, little is known about the accuracy of LD estimator relative to that of SF and about the robustness of all single‐sample estimators when some simplifying assumptions (e.g. random mating, no linkage, no genotyping errors) are violated. This study fills the gaps and uses extensive simulations to compare the biases and accuracies of the four estimators for different population properties (e.g. bottlenecks, nonrandom mating, haplodiploid), marker properties (e.g. linkage, polymorphisms) and sample properties (e.g. numbers of individuals and markers) and to compare the robustness of the four estimators when marker data are imperfect (with allelic dropouts). Extensive simulations show that SF estimator is more accurate, has a much wider application scope (e.g. suitable to nonrandom mating such as selfing, haplodiploid species, dominant markers) and is more robust (e.g. to the presence of linkage and genotyping errors of markers) than the other estimators. An empirical data set from a Yellowstone grizzly bear population was analysed to demonstrate the use of the SF estimator in practice.     

A multistage testing strategy for detection of quantitative trait Loci affecting disease resistance in Atlantic salmon

A multistage testing strategy to detect QTL for resistance to infectious salmon anemia (ISA) in Atlantic salmon is proposed. First, genotyping of amplified fragment length polymorphisms (AFLP) and a transmission disequilibrium test (TDT) were carried out using dead offspring from a disease resistance challenge test. Second, AFLP genotyping among survivors followed by a Mendelian segregation test was performed. Third, within-family survival analyses using all offspring were developed and applied to significant TDT markers with Mendelian inheritance. Maximum-likelihood methodology was developed for TDT with dominant markers to exploit linkage disequilibrium within families. The strategy was tested with two full-sib families of Atlantic salmon sired by the same male and consisting of 79 offspring in total. All dead offspring from the two families were typed for 64 primer combinations, resulting in 340 scored markers. There were 26 significant results out of 401 TDTs using dead offspring. In the second stage, only 17 marker families showed Mendelian segregation and were tested in survival analysis. A permutation test was performed for all survival analyses to compute experimentwise P-values. Two markers, aaccac356 and agccta150, were significant at P < 0.05 when accounting for multiple testing in the survival analyses. The proposed strategy might be more powerful than current mapping strategies because it reduces the number of tests to be performed in the last testing stage.     

A graph-theoretic approach to the partition of individuals into full-sib families

We present an algorithm to partition a single generation of individuals into full-sib families using single-locus co-dominant marker data. Pairwise likelihood ratios are used to create a graph that represents the full-sib relationships within the data set. Connected-component and minimum-cut algorithms from the graph theory are then employed to find the full-sib families within the graph. The results of a large-scale simulation study show that the algorithm is able to produce accurate partitions when applied to data sets with eight or more loci. Although the algorithm performs best when the distribution of allele frequencies and family sizes in a data set is uniform, the inclusion of more loci or alleles per locus allows accurate partitions to be created from data sets in which these distributions are highly skewed.     

The development of two flanking SCAR markers linked to a sex determination locus in Salix viminalis L

Most studies of sex determination systems in plants involve dioecious annuals that have known sex chromosomes. Despite the absence of such structures in the majority of dioecious plants, gender seems to be under relatively strict genetic control in some species. Genetic markers linked to a female sex-determination locus in Salix viminalis L. have been discovered through bulked segregant analysis of three full-sib families using approximately 1,000 arbitrary primers. Two RAPD markers that were present in the common female parent as well as in predominantly female progeny of these families were subsequently sequenced and converted to sequence characterized amplified region (SCAR) markers. The two SCAR markers are correlated with gender in the three full-sib families and are present in 96.4% of the female progeny and 2.2% of the males, providing evidence of linkage to a putative female-specific locus associated with gender determination in S. viminalis. Estimates of recombination suggest that the two markers are flanking a putative sex determination locus, SDL-II, in certain families of S. viminalis.     

SSRs,SNPs and DArTs comparison on estimation of relatedness and genetic parameters’ precision from a small half-sib sample population of <Emphasis Type="Italic">Eucalyptus grandis</Emphasis>

Simple sequence repeats (SSR) are the most widely used molecular markers for relatedness inference due to their multi-allelic nature and high informativeness. However, there is a growing trend toward using high-throughput and inter-specific transferable single-nucleotide polymorphisms (SNP) and Diversity Arrays Technology (DArT) in forest genetics owing to their wide genome coverage. We compared the efficiency of 15 SSRs, 181 SNPs and 2816 DArTs to estimate the relatedness coefficients, and their effects on genetic parameters’ precision, in a relatively small data set of an open-pollinated progeny trial of Eucalyptus grandis (Hill ex Maiden) with limited relationship from the pedigree. Both simulations and real data of Eucalyptus grandis were used to study the statistical performance of three relatedness estimators based on co-dominant markers. Relatedness estimates in pairs of individuals belonging to the same family (related) were higher for DArTs than for SNPs and SSRs. DArTs performed better compared to SSRs and SNPs in estimated relatedness coefficients in pairs of individuals belonging to different families (unrelated) and showed higher ability to discriminate unrelated from related individuals. The likelihood-based estimator exhibited the lowest root mean squared error (RMSE); however, the differences in RMSE among the three estimators studied were small. For the growth traits, heritability estimates based on SNPs yielded, on average, smaller standard errors compared to those based on SSRs and DArTs. Estimated relatedness in the realized relationship matrix and heritabilities can be accurately inferred from co-dominant or sufficiently dense dominant markers in a relatively small E. grandis data set with shallow pedigree.     

Results on quantitative trait loci for flushing date in oaks can be transferred to different segregating progenies

Flushing date (bud burst) is one of the most important traits for the adaptation to different environments and climates in the temperate zone. Because of their wide geographic distribution, Quercus robur L. and Q. petraea (Matt.) Liebl. are suitable as model plants to study the genetic basis of bud burst. QTLs (Quantitative Trait Loci) with comparatively large effects have been mapped in a former study in a Q. robur x Q. robur full-sib family (French cross). In the present study, we performed a Bulked Segregant Analysis (BSA) in the F (1) progeny comprising 144 seedlings derived from a cross between a single Q. robur tree as common seed parent and five different pollen donors both from Q. robur and Q. petraea (Q. robur x Q. spp., Diekholzen crosses). In addition, markers linked to two bud burst QTLs with comparably strong effect in the above-mentioned full-sib family (French cross) were tested for their association with bud burst in the Q. robur x Q. spp. (Diekholzen) progeny. Using three microsatellite markers as anchor points, we could map QTLs on linkage group 7 and on linkage group 2, together explaining 16.2 % of the total phenotypic variance (PVE) in 1999 and 38.1 % in 2003. Out of 10 markers that segregated in both mapping progenies, four markers including the two microsatellite markers, showed a significant effect on bud burst in both materials. At microsatellite loci ssrQpZAG1/5 (linkage group 7) and ssrQpZAG119 (linkage group 2) alleles associated with early (allele 166 bp in ssrQpZAG1/5) and late bud burst (allele 57 bp in ssrQpZAG119) in the Q. robur x Q. robur full-sib family (French cross) showed a highly significant association with the same polarity of the effect in the Q. robur x Q. spp. (Diekholzen) progeny. The usefulness of these markers for marker-assisted selection in full-sib and half-sib families is discussed.     

High-density linkage maps and sex-linked markers for the black tiger shrimp (Penaeus monodon)

We report on the construction of sex-specific high-density linkage maps and identification of sex-linked markers for the black tiger shrimp (Penaeus monodon). Overall, we identified 44 male and 43 female linkage groups (2n = 88) from the analysis of 2,306 AFLP markers segregating in three full-sib families, covering 2,378 and 2,362 cM, respectively. Twenty-one putatively homologous linkage groups, including the sex-linkage groups, were identified between the female and male linkage maps. Six sex-linked AFLP marker alleles were inherited from female parents in the three families, suggesting that the P. monodon adopts a WZ-ZZ sex-determining system. Two sex-linked AFLP markers, one of which we converted into an allele-specific assay, confirmed their association with sex in a panel of 52 genetically unrelated animals.     

Inheritance and molecular variations of PCR-SSCP fragments in pedunculate oak (Quercus robur L.)

Single-strand conformaiton polymorphism (SSCP) profiles of six PCR-amplified fragments (250–800 bp) were analyzed in three full-sib families of pedunculate oak (Quercus robur L.) and their parents. Among the six fragments, four were polymorphic and one exhibited complex patterns that were not changed by varying the SSCP conditions. The number of bands for the analyzed fragments varied between two and four among individuals regardless of fragment size. As shown by segregation data, the variation in the number of bands between trees could only be attributed to the allelic composition (homozygotes vs heterozygotes): a genotype that exhibited two bands was presumptively homozygous, wheras a genotype exhibiting three or four bands was heterozygous. Mendelian proportions were observed in all crosses for each polymorphic fragment. In one cross, we could clearly identify a null allele due to a possible mutation at a primer site. Single-base mutations and short insertion-deletions were shown to be the molecular causes of the SSCP polymorphism observed between different alleles. The use of SSCP as a technique to identify co-dominant markers of PCR fragments (up to 800 bp) is recommended for gene diversity studies or for gene mapping.     

A consensus linkage map for sugi (Cryptomeria japonica) from two pedigrees, based on microsatellites and expressed sequence tags

A consensus map for sugi (Cryptomeria japonica) was constructed by integrating linkage data from two unrelated third-generation pedigrees, one derived from a full-sib cross and the other by self-pollination of F1 individuals. The progeny segregation data of the first pedigree were derived from cleaved amplified polymorphic sequences, microsatellites, restriction fragment length polymorphisms, and single nucleotide polymorphisms. The data of the second pedigree were derived from cleaved amplified polymorphic sequences, isozyme markers, morphological traits, random amplified polymorphic DNA markers, and restriction fragment length polymorphisms. Linkage analyses were done for the first pedigree with JoinMap 3.0, using its parameter set for progeny derived by cross-pollination, and for the second pedigree with the parameter set for progeny derived from selfing of F1 individuals. The 11 chromosomes of C. japonica are represented in the consensus map. A total of 438 markers were assigned to 11 large linkage groups, 1 small linkage group, and 1 nonintegrated linkage group from the second pedigree; their total length was 1372.2 cM. On average, the consensus map showed 1 marker every 3.0 cM. PCR-based codominant DNA markers such as cleaved amplified polymorphic sequences and microsatellite markers were distributed in all linkage groups and occupied about half of mapped loci. These markers are very useful for integration of different linkage maps, QTL mapping, and comparative mapping for evolutional study, especially for species with a large genome size such as conifers.     

Linkage information content and efficiency of full-sib and half-sib designs for gene mapping

The accuracy of a genetic map depends on the amount of linkage information contained in the data set used for construction of the map. The amount of linkage information is related to the designs employed for linkage analysis. The purpose of this study was to provide general formulations for various genotyping schemes and family structures in order to evaluate the amount of linkage information in a data set. Linkage information content (LIC) was defined as the frequency of fully informative gametes, which are gametes from doubly heterozygous parents with known linkage phases. Depending on the design, LIC is based on two generations if the parental phases are determined statistically, or three generations if the parental phases are determined genetically. Different schemes were considered in deriving LIC: (1) genotyping of one parent or two parents, and (2) genotyping of two or three generation families. The LIC for a full-sib design was found to be generally greater than for a half-sib design but requires typing a large number of individuals when at least one locus has only two alleles. The efficiency of the full-sib design is reduced significantly if a sex-specific linkage map is sought.     

Comparing the distribution of RAPD and RFLP markers in a high density linkage map of sugar beet.

The distribution of RAPD markers was compared with that of RFLP markers in a high density linkage map of sugar beet. The same mapping population of 161 F2 individuals was used to generate all the marker data. The total map comprises 160 RAPD and 248 RFLP markers covering 508 cM. Both the RAPD and the RFLP markers show a high degree of clustering over the nine linkage groups. The pattern is compatible with a strong distal localization of recombination in the sugar beet. It leads generally to one major cluster of markers in the centre of each linkage group. In regions of high marker density, dominant RAPD markers present in either linkage phase and codominant RFLP markers are subclustered relative to each other. This phenomenon is shown to be attributable to: (i) effects of the mapping procedure when dominant and codominant data are combined, (ii) effects of the mapping procedure when dominant data in both linkage phases are combined, and (iii) genuine differences in the way RAPD and RFLP markers are recruited.     

Comparative mapping in loblolly and radiata pine using RFLP and microsatellite markers

Genetic linkage maps were constructed for loblolly pine (Pinus taeda L.) and radiata pine (P. radiata D. Don) using a common set of RFLP and microsatellite markers. The map for loblolly pine combined data from two full-sib families and consisted of 20 linkage groups covering 1281 cM. The map for radiata pine had 14 linkage groups and covered 1223 cM. All of the RFLP probes readily hybridise between loblolly and radiata pine often producing similar hybridisation patterns. There were in total 60 homologous RFLP loci mapped in both species which could be used for comparative purposes. A set of 20 microsatellite markers derived from radiata pine were also assayed; however, only 9 amplified and revealed polymorphic loci in both species. Single-locus RFLP and microsatellite markers were used to match up linkage groups and compare order between species. Twelve syntenic groups were obtained each consisting of from 3 to 9 homologous loci. The order of homologous loci was colinear in most cases, suggesting no major chromosomal rearrangements in the evolution of these species. Comparative mapping between loblolly and radiata pine should facilitate genetic research in both species and provide a framework for mapping in other pine species. Received: 25 November 1998 / /Accepted: 19 December 1998     

Genetic analysis using trans-dominant linked markers in an F2 family

Trans-dominant linked markers pairs (trans referring to the repulsion linkage phase) provide a model for inferring the F₂ progeny genotype based upon both the conditional probabilities of F₂ genotypes, given the F₂ phenotype, and prior information on marker arrangement. Prior information of marker arrangement can be readily obtained from a linkage analysis performed on marker segregation data in a family resulting by crossing the F₁ individual to a tester parent or else can be obtained directly from the gametes of the F₁, or from recombinant inbred lines. We showed that a trans-dominant linked marker (TDLM) pair can be recoded as a co-dominant megalocus when the recombination fraction, r₁, for apair of TDLMs is less than 0.05. We obtained a maximum-likelihood estimator (MLE) of the recombination frequency, r₂, between a TDLM pair and a co-dominant marker in an F₂ family using the EM algorithm. The MLE was biased. Mean bias increased as r₁ and r₂ increased, and decreased as sample size increased. The information content for r₂ was compared to the information content of dominant and co-dominant markers segregating in an F₂ family. It was almost identical with two co-dominant markers when r₁0.01 and r₂0.05. For larger values of r₁, (0.05r₁0.15) a TDLM pair provided 75%–66% of the information content of two co-dominant markers. Although dominant markers can be converted to co-dominant markers by a laborious process of cloning, sequencing, and PCR, TDLM pairs could easily substitute for co-dominant markers in order to detect quantitative trait loci (QTLs) and estimate gene action in an F₂ family.     

Comparative analysis of sequences of mitochondrial genomes of wild abortive male sterile (WA-CMS) and male fertile lines of rice,development of functional markers for WA-CMS trait and their use in assessment of genetic purity of seeds of WA-CMS lines

WA-CMS system based rice hybrids are widely adopted in many rice growing countries, including India. Even though it is well known that the trait is controlled by mitochondria, the genes underpinning the trait remain enigmatic. In the present study, a complete genome-wide comparative sequence analysis was performed using draft mitochondrial genomes of WA-CMS and male fertile lines in a step-wise manner, progressively covering 5–10 kb every time through BLASTN tool. The sequence polymorphisms identified in different mitochondrial regions were targeted to develop two different sets of dominant PCR-based markers, one consisting of six markers targeting WA-CMS mitochondria, the other set consisting of five markers targeting male fertile mitochondria in addition to development of a set of eight co-dominant PCR-based markers targeting both the genomes. When a set of candidate genes/ORFs reported earlier to be associated with WA-CMS trait in rice were analyzed through RT-PCR of RNA isolated from immature rice florets, it was observed that the chimeric ORF, WA352 is expressed only in WA-CMS line and hybrid (i.e. genotypes containing sterile mitochondria), indicating it’s candidacy for the WA-CMS trait. Targeting the functional nucleotide polymorphism between WA-CMS and maintainer mitochondria with respect to WA352, two dominant markers, one targeting sterile and another targeting fertile mitochondria were developed. In addition, a robust, co-dominant functional marker targeting the candidate gene was also developed and validated for its utility in identification of genetic impurities in seed lots of WA-CMS lines.     

Simultaneous maximum likelihood estimation of linkage and linkage phases in outcrossing species

With the advent of new molecular marker technologies, it is now feasible to initiate genome projects for outcrossing plant species, which have not received much attention in genetic research, despite their great agricultural and environmental value. Because outcrossing species typically have heterogeneous genomes, data structure for molecular markers representing an entire genome is complex: some markers may have more alleles than others, some markers are codominant whereas others are dominant, and some markers are heterozygous in one parent but fixed in the other parent whereas the opposite can be true for other markers. A major difficulty in analyzing these different types of marker at the same time arises from uncertainty about parental linkage phases over markers. In this paper, we present a general maximum-likelihood-based algorithm for simultaneously estimating linkage and linkage phases for a mixed set of different marker types containing fully informative markers (segregating 1:1:1:1) and partially informative markers (or missing markers, segregating 1:2:1, 3:1, and 1:1) in a full-sib family derived from two outbred parent plants. The characterization of linkage phases is based on the posterior probability distribution of the assignment of alternative alleles at given markers to two homologous chromosomes of each parent, conditional on the observed phenotypes of the markers. Two- and multi-point analyses are performed to estimate the recombination fraction and determine the most likely linkage phase between different types of markers. A numerical example is presented to demonstrate the statistical properties of the model for characterizing the linkage phase between markers.     

Genetic linkage studies with cleft lip and palate: report of two family studies

Genetic linkage studies are reported on two families with cleft lip +/- cleft palate. For the first family (LP01) the etiology of the clefting is unknown, and the linkage analyses were done assuming both autosomal dominant and autosomal recessive inheritance. Close linkage is rejected with the Duffy blood group under the dominant model and with four loci (Duffy, Kidd, and ABO blood groups and haptoglobin) under the recessive model. The second family (LP02) is a Mexican-American family segregating the van der Woude syndrome with lip pits. The linkage analyses for this autosomal dominant trait excluded close linkage with seven genetic markers, including three on chromosome one. The maximum lod scores were 0.6 with BF (chromosome 6) and 0.4 with the P blood group, which is not yet mapped.