Identification of quantitative trait loci affecting gastrointestinal parasite resistance in an experimental Angus population

DNA markers associated with quantitative trait loci (QTL) affecting host tolerance to gastrointestinal (GI) parasite infection are ideal targets for marker‐assisted selection. However, few studies in cattle have attempted to identify this type of QTL due to the difficulty of generating accurate phenotypic data from a resource population with adequate statistical power for detection. For this effort, we amassed fecal egg count (FEC) measures from annual natural field challenges with GI nematodes that spanned 12 different contemporary groups of Angus calves (1992–2000) derived from a closed breeding population. FEC and blood pepsinogen measures were taken weekly over a 26‐week period post‐weaning, and the FEC data were Box‐Cox transformed to normalize the distribution of phenotypes. These 305 test animals and more than 100 founding animals from the extended pedigree were genotyped across 190 microsatellites markers. The genome‐wide analyses identified a suggestive genome‐wide QTL on bovine chromosome (Chr) 8 (P < 0.002) and nominal QTL on Chr 4, 12 and 17 (P < 0.05). These findings were unique for cattle, and some corresponded to previously identified QTL locations for parasite‐related traits in sheep to provide genome locations for further fine mapping of parasite resistance/susceptibility in Angus cattle.     

A sampling algorithm for segregation analysis

Methods for detecting Quantitative Trait Loci (QTL) without markers have generally used iterative peeling algorithms for determining genotype probabilities. These algorithms have considerable shortcomings in complex pedigrees. A Monte Carlo Markov chain (MCMC) method which samples the pedigree of the whole population jointly is described. Simultaneous sampling of the pedigree was achieved by sampling descent graphs using the Metropolis-Hastings algorithm. A descent graph describes the inheritance state of each allele and provides pedigrees guaranteed to be consistent with Mendelian sampling. Sampling descent graphs overcomes most, if not all, of the limitations incurred by iterative peeling algorithms. The algorithm was able to find the QTL in most of the simulated populations. However, when the QTL was not modeled or found then its effect was ascribed to the polygenic component. No QTL were detected when they were not simulated.     

Detection of QTL for immune response to sheep red blood cells in laying hens

The aim of this study is to detect quantitative trait loci (QTL) involved in the regulation of the primary and the secondary immune response to sheep red blood cells (SRBC) in a resource population using microsatellite DNA markers. The F2 resource population originates from a cross of two divergently selected lines for either high (H line) or low (L line) primary antibody response to SRBC. The F2 population consisted of six half-sib families, three families per each of reciprocal crosses. Total antibody titres to SRBC were determined by agglutination in serum from all birds. F2, F1 and F0 generations were genotyped for 170 microsatellite markers, using a whole-genome scan approach. The half-sib and the line-cross analyses were performed to determine QTL regions associated with regulation of the immune response. In the half-sib analysis, four QTL for SRBC primary response have been identified: on GGA3, GGA5, GGA16 and GGA23. No QTL was identified for SRBC secondary response under the half-sib model. In the line-cross analysis, three QTL were identified on GGA10, GGA16 and GGA27 for SRBC primary response and five QTL were identified on GGA6, GGA9, GGA15, GGA16 and GGA27 for SRBC secondary response. Subsequently, the family contribution of individual families to the QTL was analysed. The family with the largest contribution was genotyped with additional microsatellite markers in the QTL region on GGA5. The extended half-sib analysis with additional genotype information results in narrowing down the QTL region on GGA5.     

The covariance between relatives conditional on genetic markers

The development of molecular genotyping techniques makes it possible to analyze quantitative traits on the basis of individual loci. With marker information, the classical theory of estimating the genetic covariance between relatives can be reformulated to improve the accuracy of estimation. In this study, an algorithm was derived for computing the conditional covariance between relatives given genetic markers. Procedures for calculating the conditional relationship coefficients for additive, dominance, additive by additive, additive by dominance, dominance by additive and dominance by dominance effects were developed. The relationship coefficients were computed based on conditional QTL allelic transmission probabilities, which were inferred from the marker allelic transmission probabilities. An example data set with pedigree and linked markers was used to demonstrate the methods developed. Although this study dealt with two QTLs coupled with linked markers, the same principle can be readily extended to the situation of multiple QTL. The treatment of missing marker information and unknown linkage phase between markers for calculating the covariance between relatives was discussed.     

Mapping of quantitative trait loci controlling timing of bud flush in Salix

Dormancy release is an important phenological stage, which determines plant growth and survival in northern temperate regions. Spring bud flushing was studied in a Salix pedigree (n=82) derived from a cross between the male hybrid clone "Bj?rn" (Salix viminalis x Salix schwerinii) and the female clone "78183" (Salix viminalis). The timing of bud flush was recorded outdoors in two consecutive years (1998, 1999) and indoor in the spring of 1998. Timing of bud flush was found to be under moderately strong genetic control (clonal mean heritabilities ranging from 0.43 to 0.72). Phenotypic correlations between height growth and bud flushing were negative but non-significant (r=0.1-0.3). Using a Salix linkage map composed of 325 AFLP and 38 RFLP markers, six quantitative trait loci (QTLs) and three unmapped marker loci associated with timing of bud flush were detected. Four QTLs were detected in the field experiment while two QTLs and three unmapped marker loci were identified in the indoor experiment. One QTL associated with indoor bud flushing coincided with one of the QTL detected from the field data. Individual QTL explained 6-16% of the phenotypic variance [corrected]. None of the bud flush QTLs coincided with QTLs controlling height growth identified previously in the same pedigree.     

A Triad of Highly Divergent Polymeric Immunoglobulin Receptor (PIGR) Haplotypes with Major Effect on IgA Concentration in Bovine Milk

The aim of this study was to determine a genetic basis for IgA concentration in milk of Bos taurus. We used a Holstein-Friesian x Jersey F2 crossbred pedigree to undertake a genome-wide search for QTL influencing IgA concentration and yield in colostrum and milk. We identified a single genome-wide significant QTL on chromosome 16, maximising at 4.8 Mbp. The polymeric immunoglobulin receptor gene (PIGR) was within the confidence interval of the QTL. In addition, mRNA expression analysis revealed a liver PIGR expression QTL mapping to the same locus as the IgA quantitative trait locus. Sequencing and subsequent genotyping of the PIGR gene revealed three divergent haplotypes that explained the variance of both the IgA QTL and the PIGR expression QTL. Genetic selection based on these markers will facilitate the production of bovine herds producing milk with higher concentrations of IgA.     

Genetic structure of a QTL hotspot on chromosome 2 in sweet cherry indicates positive selection for favorable haplotypes

A genomic region of particular interest for sweet cherry (Prunus avium L.) breeding is a quantitative trait locus (QTL) “hotspot” on chromosome 2. QTLs for fruit size, firmness, sweetness, and flowering time are reported to map to this region. An understanding of genetic diversity, allele sources, linkage relationships, and historical recombinations is critical to enable the combining of favorable alleles at multiple loci. The objectives of this study were to characterize, visualize, and interpret the genetic structure of this previously identified QTL hotspot within North American sweet cherry breeding germplasm, using a pedigree-based haploblocking approach. Across the 29.4 cM (6.3 Mbp) region defined by single nucleotide polymorphism (SNP) information from the RosBREED cherry 6K SNP array v1, a total of 12 recombination events falling into six inter-marker regions were traced within the pedigree of elite and wild germplasm (n = 55). These recombinations defined five haploblocks containing 5–15 markers and exhibiting 7–11 haplotypes each. Over the entire QTL hotspot, 30 extended haplotypes were identified for which parental gametes could be determined. When the haploblocks and their haplotypes were used to explore genetic diversity, ancestry, and recombination patterns, and then integrated with previous QTL results for fruit size, the results indicated that favorable alleles at this QTL hotspot are under positive selection in breeding. The genetic framework provided by a haploblock approach and knowledge of haplotype-level diversity sets the stage for assigning breeding utility to these haplotypes.     

Transmission of important chromosomal regions under selection revealed in rice pedigree breeding programs

Genetic analysis across a whole plant genome based on pedigree information offers considerable potential for enhancing genetic gain from plant breeding programs through quantitative trait loci (QTL) mapping and marker-assisted selection. Here, we report its application for graphically genotyping varieties used in Chinese japonica rice (Oryza sativa L.) pedigree breeding programs. We identified 34 important chromosomal regions from the founder parent that are under selection in the breeding programs, and by comparing donor genomic regions that are under selection with QTL locations of agronomic traits, we found that QTL clustered in important genomic regions, in accordance with association analyses of natural populations and other previous studies. The convergence of genomic regions under selection with QTL locations suggests that donor genomic regions harboring key genes/QTL for important agronomic traits have been selected by plant breeders since the 1950s from the founder rice plants. The results provide better understanding of the effects of selection in breeding programs on the traits of rice cultivars. They also provide potentially valuable information for enhancing rice breeding programs through screening candidate parents for targeted molecular markers, improving crop yield potential and identifying suitable genetic material for use in future breeding programs.     

Genetic mapping of quantitative trait loci for aseasonal reproduction in sheep

The productivity and economic prosperity of sheep farming could benefit greatly from more effective methods of selection for year-round lambing. Identification of QTL for aseasonal reproduction in sheep could lead to more accurate selection and faster genetic improvement. One hundred and twenty microsatellite markers were genotyped on 159 backcross ewes from a Dorset × East Friesian crossbred pedigree. Interval mapping was undertaken to map the QTL underlying several traits describing aseasonal reproduction including the number of oestrous cycles, maximum level of progesterone prior to breeding, pregnancy status determined by progesterone level, pregnancy status determined by ultrasound, lambing status and number of lambs born. Seven chromosomes (1, 3, 12, 17, 19, 20 and 24) were identified to harbour putative QTL for one or more component traits used to describe aseasonal reproduction. Ovine chromosomes 12, 17, 19 and 24 harbour QTL significant at the 5% chromosome-wide level, chromosomes 3 and 20 harbour QTL that exceeded the threshold at the 1% chromosome-wide level, while the QTL identified on chromosome 1 exceeded the 1% experiment-wide significance level. These results are a first step towards understanding the genetic mechanism of this complex trait and show that variation in aseasonal reproduction is associated with multiple chromosomal regions.     

Mapping of multiple quantitative trait loci affecting bovine spongiform encephalopathy

A whole-genome scan was conducted to map quantitative trait loci (QTL) for BSE resistance or susceptibility. Cows from four half-sib families were included and 173 microsatellite markers were used to construct a 2835-cM (Kosambi) linkage map covering 29 autosomes and the pseudoautosomal region of the sex chromosome. Interval mapping by linear regression was applied and extended to a multiple-QTL analysis approach that used identified QTL on other chromosomes as cofactors to increase mapping power. In the multiple-QTL analysis, two genome-wide significant QTL (BTA17 and X/Y(ps)) and four genome-wide suggestive QTL (BTA1, 6, 13, and 19) were revealed. The QTL identified here using linkage analysis do not overlap with regions previously identified using TDT analysis. One factor that may explain the disparity between the results is that a more extensive data set was used in the present study. Furthermore, methodological differences between TDT and linkage analyses may affect the power of these approaches.     

Validating a QTL region characterized by multiple haplotypes

Validation of a quantitative trait locus (QTL) for outcrossing perennial plants is rarely reported due to complexity of plausible genetic models and reliance on field designs already available. Here, a particular marker-QTL haplotype exerted a large, positive effect on height for Pinus taeda and its origin could be traced to a founder, GP₃, in a three-generation QTL pedigree. To validate this QTL effect, we used an extended GP₃-based pedigree. In the validation cross, each of the 46 offspring was clonally propagated from developing seeds using somatic embryogenesis technology. Subsequent analyses were conducted separately for seedlings and for other somatic emblings. For seedlings, the original QTL effect could not be fully validated. For somatic emblings, a strong negative QTL effect was detected in the validation cross; some evidence from another cross-supported the original positive QTL effect. From this part of the analysis, three distinct marker-QTL haplotypes at a single locus could be inferred. Validating QTL haplotypes in readily available field tests was feasible despite the genetic model complexity inherent to outcrossing long-lived perennials.     

Localization of a quantitative trait locus via a Bayesian approach

A Bayesian approach to the direct mapping of a quantitative trait locus (QTL), fully utilizing information from multiple linked gene markers, is presented in this paper. The joint posterior distribution (a mixture distribution modeling the linkage between a biallelic QTL and N gene markers) is computationally challenging and invites exploration via Markov chain Monte Carlo methods. The parameter's complete marginal posterior densities are obtained, allowing a diverse range of inferences. Parameters estimated include the QTL genotype probabilities for the sires and the offspring, the allele frequencies for the QTL, and the position and additive and dominance effects of the QTL. The methodology is applied through simulation to a half-sib design to form an outbred pedigree structure where there is an entire class of missing information. The capacity of the technique to accurately estimate parameters is examined for a range of scenarios.     

The role of pedigree information in combined linkage disequilibrium and linkage mapping of quantitative trait loci in a general complex pedigree

Combined linkage disequilibrium and linkage (LDL) mapping can exploit historical as well as recent and observed recombinations in a recorded pedigree. We investigated the role of pedigree information in LDL mapping and the performance of LDL mapping in general complex pedigrees. We compared using complete and incomplete genotypic data, spanning 5 or 10 generations of known pedigree, and we used bi- or multiallelic markers that were positioned at 1- or 5-cM intervals. Analyses carried out with or without pedigree information were compared. Results were compared with linkage mapping in some of the data sets. Linkage mapping or LDL mapping with sparse marker spacing ( approximately 5 cM) gave a poorer mapping resolution without considering pedigree information compared to that with considering pedigree information. The difference was bigger in a pedigree of more generations. However, LDL mapping with closely linked markers ( approximately 1 cM) gave a much higher mapping resolution regardless of using pedigree information. This study shows that when marker spacing is dense and there is considerable linkage disequilibrium generated from historical recombinations between flanking markers and QTL, the loss of power due to ignoring pedigree information is negligible and mapping resolution is very high.     

Haplotypic QTL mapping in an outbred pedigree

An offspring genome can be viewed as a mosaic of chromosomal segments or haplotypes contributed by multiple founders in any quantitative trait locus (QTL) detection study but tracing these is especially complex to achieve for outbred pedigrees. QTL haplotypes can be traced from offspring back to individual founders in outbred pedigrees by combining founder-origin probabilities with fully informative flanking markers. This haplotypic method was illustrated for QTL detection using a three-generation pedigree for a woody perennial plant, Pinus taeda L. Growth rate was estimated using height measurements from ages 2 to 10 years. Using simulated and actual datasets, power of the experimental design was shown to be efficient for detecting QTLs of large effect. Using interval mapping and fully informative markers, a large QTL accounting for 11.3% of the phenotypic variance in the growth rate was detected. This same QTL was expressed at all ages for height, accounting for 7.9-12.2% of the phenotypic variance. A mixed-model inheritance was more appropriate for describing genetic architecture of growth curves in P. taeda than a strictly polygenic model. The positive QTL haplotype was traced from the offspring to its contributing founder, GP3, then the haplotypic phase for GP3 was determined by assaying haploid megagametophytes. The positive QTL haplotype was a recombinant haplotype contributed by GP3. This study illustrates the combined power of fully informative flanking markers and founder origin probabilities for (1) estimating QTL haplotype magnitude, (2) tracing founder origin and (3) determining haplotypic transmission frequency.     

Prediction of identity by descent probabilities from marker-haplotypes

The prediction of identity by descent (IBD) probabilities is essential for all methods that map quantitative trait loci (QTL). The IBD probabilities may be predicted from marker genotypes and/or pedigree information. Here, a method is presented that predicts IBD probabilities at a given chromosomal location given data on a haplotype of markers spanning that position. The method is based on a simplification of the coalescence process, and assumes that the number of generations since the base population and effective population size is known, although effective size may be estimated from the data. The probability that two gametes are IBD at a particular locus increases as the number of markers surrounding the locus with identical alleles increases. This effect is more pronounced when effective population size is high. Hence as effective population size increases, the IBD probabilities become more sensitive to the marker data which should favour finer scale mapping of the QTL. The IBD probability prediction method was developed for the situation where the pedigree of the animals was unknown (i.e. all information came from the marker genotypes), and the situation where, say T, generations of unknown pedigree are followed by some generations where pedigree and marker genotypes are known.     

Fruit size QTL identification and the prediction of parental QTL genotypes and breeding values in multiple pedigreed populations of sweet cherry

Large fruit size is a critical trait for any new sweet cherry (Prunus avium L.) cultivar, as it is directly related to grower profitability. Therefore, determining the genetic control of fruit size in relevant breeding germplasm is a high priority. The objectives of this study were (1) to determine the number and positions of quantitative trait loci (QTL) for sweet cherry fruit size utilizing data simultaneously from multiple families and their pedigreed ancestors, and (2) to estimate fruit size QTL genotype probabilities and genomic breeding values for the plant materials. The sweet cherry material used was a five-generation pedigree consisting of 23 founders and parents and 424 progeny individuals from four full-sib families, which were phenotyped for fruit size and genotyped with 78 RosCOS single nucleotide polymorphism and 86 simple sequence repeat markers. These data were analyzed by a Bayesian approach implemented in FlexQTL? software. Six QTL were identified: three on linkage group (G) 2 with one each on groups 1, 3, and 6. Of these QTL, the second G2 QTL and the G6 QTL were previously discovered while other QTL were novel. The predicted QTL genotypes show that some QTL were segregating in all families while other QTL were segregating in a subset of the families. The progeny varied for breeding value, with some progeny having higher breeding values than their parents. The results illustrate the use of multiple pedigree-linked families for integrated QTL mapping in an outbred crop to discover novel QTL and predict QTL genotypes and breeding values.     

Using dominance relationship coefficients based on linkage disequilibrium and linkage with a general complex pedigree to increase mapping resolution

Dominance (intralocus allelic interactions) plays often an important role in quantitative trait variation. However, few studies about dominance in QTL mapping have been reported in outbred animal or human populations. This is because common dominance effects can be predicted mainly for many full sibs, which do not often occur in outbred or natural populations with a general pedigree. Moreover, incomplete genotypes for such a pedigree make it infeasible to estimate dominance relationship coefficients between individuals. In this study, identity-by-descent (IBD) coefficients are estimated on the basis of population-wide linkage disequilibrium (LD), which makes it possible to track dominance relationships between unrelated founders. Therefore, it is possible to use dominance effects in QTL mapping without full sibs. Incomplete genotypes with a complex pedigree and many markers can be efficiently dealt with by a Markov chain Monte Carlo method for estimating IBD and dominance relationship matrices (D(RM)). It is shown by simulation that the use of D(RM) increases the likelihood ratio at the true QTL position and the mapping accuracy and power with complete dominance, overdominance, and recessive inheritance modes when using 200 genotyped and phenotyped individuals.     

A new test for family-based association mapping with inbred lines from plant breeding programs

Association mapping holds great promise for the detection of quantitative trait loci (QTL) in plant breeding populations. The main objectives of this study were to (1) adapt the quantitative pedigree disequilibrium test to typical pedigrees of inbred lines produced in plant breeding programs, (2) compare the newly developed quantitative inbred pedigree disequilibrium test (QIPDT) with the commonly employed logistic regression ratio test (LRRT), with respect to the power and type I error rate of QTL detection, and (3) demonstrate the use of the QIPDT by applying it to flowering data of European elite maize inbreds. QIPDT and LRRT were compared based on computer simulations modeling 55 years of hybrid maize breeding in Central Europe. Furthermore, we applied QIPDT to a cross-section of 49 European elite maize inbred lines genotyped with 722 amplified fragment length polymorphism markers and phenotyped in four environments for days to anthesis. Compared to LRRT, the power to detect QTL was higher with QIPDT when using data collected routinely in plant breeding programs. Application of QIPDT to the 49 European maize inbreds resulted in a significant (P < 0.05) association located at a position for which a consensus QTL was detected in a previous study. The results of our study suggested that QIPDT is a promising QTL detection method for data collected routinely in plant breeding programs.     

Mapping quantitative trait loci for bovine ovulation rate

An elite, three-generation family from the USDA Meat Animal Research Center twinning population was examined for evidence of ovulation rate quantitative trait loci (QTL). This work was both a continuation of previously reported results suggesting evidence for ovulation rate QTL on bovine Chromosome (Chr) 7 and an extension of a genome-wide search for QTL. Additional markers were typed on Chr 7 to facilitate interval mapping and testing of the hypothesis of one versus two QTL on that chromosome. In addition, 14 other informative markers were added to a selective genotyping genome screening of this family, and markers exhibiting nominal significance were used to identify chromosomal regions that were then subjected to more exhaustive analysis. For Chr 7, a total of 12 markers were typed over a region spanning the proximal two-thirds of the chromosome. Results from interval mapping analyses indicated evidence suggestive of the presence of QTL (nominal P < 0.00077) within this region. Subsequent analysis with a model postulating two QTL provided evidence (P < 0.05) for two rather than one QTL on this chromosome. Preliminary analysis with additional markers indicated nominal significance (P < 0.05) for regions of Chrs 5, 10, and 19. Each of these regions was then typed with additional markers for the entire three-generation pedigree. Significant evidence (P < 0.000026) of ovulation rate QTL was found for Chrs 5 and 19, while support on Chr 10 failed to exceed a suggestive linkage threshold (P > 0.00077). Received: 14 May 1999 / Accepted: 14 October 1999     

Combining two Meishan F2 crosses improves the detection of QTL on pig chromosomes 2, 4 and 6

Background

In pig, a number of experiments have been set up to identify QTL and a multitude of chromosomal regions harbouring genes influencing traits of interest have been identified. However, the mapping resolution remains limited in most cases and the detected QTL are rather inaccurately located. Mapping accuracy can be improved by increasing the number of phenotyped and genotyped individuals and/or the number of informative markers. An alternative approach to overcome the limited power of individual studies is to combine data from two or more independent designs.

Methods

In the present study we report a combined analysis of two independent design (a French and a Dutch F2 experimental designs), with 2000 F2 individuals. The purpose was to further map QTL for growth and fatness on pig chromosomes 2, 4 and 6. Using QTL-map software, uni- and multiple-QTL detection analyses were applied separately on the two pedigrees and then on the combination of the two pedigrees.

Results

Joint analyses of the combined pedigree provided (1) greater significance of shared QTL, (2) exclusion of false suggestive QTL and (3) greater mapping precision for shared QTL.

Conclusions

Combining two Meishan x European breeds F2 pedigrees improved the mapping of QTL compared to analysing pedigrees separately. Our work was facilitated by the access to raw phenotypic data and DNA of animals from both pedigrees and the combination of the two designs with the addition of new markers allowed us to fine map QTL without phenotyping additional animals.