Fine mapping and physical characterization of two linked quantitative trait loci affecting milk fat yield in dairy cattle on BTA26

Previously, a highly significant QTL affecting fat yield and protein yield and mapped to the bovine BTA26 chromosome has been reported to segregate in the French Holstein cattle population. To confirm and refine the location of this QTL, the original detection experiment was extended by adding 12 new families and genotyping 25 additional microsatellite markers (including 11 newly developed markers). Data were then analyzed by an approach combining both linkage and linkage disequilibrium information, making it possible to identify two linked QTL separated by 20 cM corresponding to approximately 29 Mb. The presence of a QTL affecting protein yield was confirmed but its position was found to be more telomeric than the two QTLunderlying fat yield. Each identified QTL affecting milk fat yield was physically mapped within a segment estimated to be <500 kb. Two strong functional candidate genes involved, respectively, in fatty acid metabolism and membrane permeability were found to be localized within this segment while other functional candidate genes were discarded. A haplotype comprising the favorable allele at each QTL position appears to be overrepresented in the artificial insemination bull population.     

Fine mapping QTL for female fertility on BTA04 and BTA13 in dairy cattle using HD SNP and sequence data

Background

Female fertility is important for the maintenance of the production in a dairy cattle herd. Two QTL regions on BTA04 and on BTA13 previously detected in Nordic Holstein (NH) and validated in the Danish Jersey (DJ) and Nordic Red (NR) were investigated further in the present study to further refine the QTL locations. Refined QTL regions were imputed to the full sequence data. The genes in the regions were then studied to ascertain their possible effect on fertility traits.

Results

BTA04 was screened for number of inseminations (AIS), 56-day non-return rate (NRR), days from first to last insemination (IFL), and the interval from calving to first insemination (ICF) in the range of 38,257,758 to 40,890,784 bp, whereas BTA13 was screened for ICF only in the range from 21,236,959 to 46,150,079 with the HD bovine SNP array for NH, DJ and NR. No markers in the DJ and NR breeds reached significance. By analyzing imputed sequence data the QTL position on BTA04 was narrowed down to two regions in the NH. In these two regions a total of 9 genes were identified. BTA13 was analyzed using sequence data for the NH breed. The highest –log₁₀(P-value) was 19.41 at 33,903,159 bp. Two regions were identified: Region 1: 33,900,143-33,908,994 bp and Region 2: 34,051,815-34,056,728 bp. SNPs within and between these two regions were annotated as intergenic.

Conclusion

Screening BTA04 and BTA13 for female fertility traits in NH, NR and DJ suggested that the QTL for female fertility were specific for NH. A missense mutation in CD36 showed the strongest association with fertility traits on BTA04. The annotated SNPs on BTA13 were all intergenic variants. It is possible that BTA13 at this stage is poorly annotated such that the associated polymorphisms are located in as-yet undiscovered genes. Fertility traits are complex traits as many different biological and physiological factors determine whether a cow is fertile. Therefore it is not expected that there is a simple explanation with an obvious candidate gene but it is more likely a network of genes and intragenic variants that explain the variation of these traits.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-790) contains supplementary material, which is available to authorized users.     

Genome-wide association mapping in Norwegian Red cattle identifies quantitative trait loci for fertility and milk production on BTA12

Reproductive performance is a critical trait in dairy cattle. Poor reproductive performance leads to prolonged calving intervals, higher culling rates and extra expenses related to multiple inseminations, veterinary treatments and replacements. Genetic gain for improved reproduction through traditional selection is often slow because of low heritability and negative correlations with production traits. Detection of DNA markers associated with improved reproductive performance through genome-wide association studies could lead to genetic gain that is more balanced between fertility and production. Norwegian Red cattle are well suited for such studies, as very large numbers of detailed reproduction records are available. We conducted a genome-wide association study for non-return rate, fertility treatments and retained placenta using almost 1 million records on these traits and 17 343 genome-wide single-nucleotide polymorphisms. Genotyping costs were minimized by genotyping the sires of the cows recorded and by using daughter averages as phenotypes. The genotyped sires were assigned to either a discovery or a validation population. Associations were only considered to be validated if they were significant in both groups. Strong associations were found and validated on chromosomes 1, 5, 8, 9, 11 and 12. Several of these were highly supported by findings in other studies. The most important result was an association for non-return rate in heifers in a region of BTA12 where several associations for milk production traits have previously been found. Subsequent fine-mapping verified the presence of a quantitative trait loci (QTL) having opposing effects on non-return rate and milk production at 18 Mb. The other reproduction QTL did not have pleiotropic effects on milk production, and these are therefore of considerable interest for use in marker-assisted selection.     

Identification of quantitative trait loci affecting cattle temperament

In addition to its potential contribution to improving animal welfare, the study of the genetics of cattle behavior may provide more general insights into the genetic control of such complex traits. We carried out a genome scan in a Holstein x Charolais cross cattle population to identify quantitative trait loci (QTL) influencing temperament-related traits. Individuals belonging to the second-generation of this population (F(2) and backcross individuals) were subjected to 2 behavioral tests. The flight from feeder (FF) test measured the distance at which the animal moved away from an approaching human observer, whereas the social separation (SS) test categorized different activities which the animal engaged in when removed from its penmates. The entire population was genotyped with 165 microsatellite markers. A regression interval mapping analysis identified 29 regions exceeding the 5% chromosome-wide significance level, which individually explained a relatively small fraction of the phenotypic variance of the traits (from 3.8% to 8.4%). One of the significant associations influencing an FF test trait on chromosome 29 reached the 5% genome-wide significance level. Eight other QTL, all associated with an SS test trait, reached the 1% chromosome-wide significance level. The location of some QTL coincided with other previously reported temperament QTL in cattle, whereas those that are reported for the first time here may represent general loci controlling temperament differences between cattle breeds. No overlapping QTL were identified for the traits measured by the 2 different tests, supporting the hypothesis that different genetic factors influence behavioral responses to different situations.     

Fine mapping of quantitative trait loci using linkage disequilibria with closely linked marker loci

A multimarker linkage disequilibrium mapping method was developed for the fine mapping of quantitative trait loci (QTL) using a dense marker map. The method compares the expected covariances between haplotype effects given a postulated QTL position to the covariances that are found in the data. The expected covariances between the haplotype effects are proportional to the probability that the QTL position is identical by descent (IBD) given the marker haplotype information, which is calculated using the genedropping method. Simulation results showed that a QTL was correctly positioned within a region of 3, 1.5, or 0.75 cM in 70, 62, and 68%, respectively, of the replicates using markers spaced at intervals of 1, 0.5, and 0.25 cM, respectively. These results were rather insensitive to the number of generations since the QTL occurred and to the effective population size, except that 10 generations yielded rather poor estimates of the QTL position. The position estimates of this multimarker disequilibrium mapping method were more accurate than those from a single marker transmission disequilibrium test. A general approach for identifying QTL is suggested, where several stages of disequilibrium mapping are used with increasingly dense marker spacing.     

肉牛的数量性状基因座定位研究进展

遗传图谱的发展为寻找和定位影响重要数量性状变异的基因提供了便利。迄今为止,育种学家已经在肉牛的1、2、5、6、14、15、17、18、19、21、23、27、和29号常染色体上发现了QTL的踪迹。候选基因的研分显示肌肉生长抑制素基因等可能就是生长和屠宰重性状的QTL,基困组统计定位则揭示最有可能的QTL区域在2、5、6、15、19、27、29号染色体上。进一步的定位仍需遗传学家、分子生物学家及育种学家的共同努力。     

Fine mapping quantitative trait loci affecting milk production traits on bovine chromosome 6 in a Chinese Holstein population

To fine map the previously detected quantitative trait loci （QTLs） affecting milk production traits on bovine chromosome 6 （BTA6）, 15 microsatellite markers situated within an interval of 14.3 cM spanning from BMS690 to BM4528 were selected and 918 daughters of 8 sires were genotyped. Two mapping approaches, haplotype sharing based LD mapping and single marker regression mapping, were used to analyze the data. Both approaches revealed a quantitative trait locus （QTL） with significant effects on milk yield, fat yield and protein yield located in the segment flanked by markers BMS483 and MNB209, which spans a genetic distance of 0.6 cM and a physical distance of 1.5 Mb. In addition, the single marker regression mapping also revealed a QTL affecting fat percentage and protein percentage at marker DIK2291. Our fine mapping work will facilitate the cloning of candidate genes underlying the QTLs for milk production traits.     

Further studies on using multiple-cross mapping (MCM) to map quantitative trait loci

We have completed whole-genome scans for quantitative trait loci (QTLs) associated with acute ethanol-induced activation in the six F₂ intercrosses that can be formed from the C57BL/6J (B6), DBA/2J (D2) , BALB/cJ (C), and LP/J (LP) inbred strains. The goal was to test the hypothesis that given the relatively simple structure of the laboratory mouse genome, the same QTLs will be detected in multiple crosses which in turn will provide support for the strategy of multiple-cross mapping (MCM). QTLs with LOD scores greater than 4 were detected on Chrs 1, 2, 3, 8, 9, 13, 14, and 16. Only for the QTL on distal Chr 1 was there convincing evidence that the same or at least a very similar QTL was detected in multiple crosses. We also mapped the Chr 2 QTL directly in heterogeneous stock (HS) animals derived from the four inbred strains. At G₁₉ the QTL was mapped to an approximately 3-Mbp interval and this interval was associated with a haplotype block with a largely biallelic structure: B6-L:C-D2. We conclude that mapping in HS animals not only provides significantly greater QTL resolution, at least in some cases it provides significantly more information about the QTL haplotype structure.     

Quantitative trait loci mapping in dairy cattle: review and meta-analysis

From an extensive review of public domain information on dairy cattle quantitative trait loci (QTL), we have prepared a draft online QTL map for dairy production traits. Most publications (45 out of 55 reviewed) reported QTL for the major milk production traits (milk, fat and protein yield, and fat and protein concentration (%)) and somatic cell score. Relatively few QTL studies have been reported for more complex traits such as mastitis, fertility and health. The collated QTL map shows some chromosomal regions with a high density of QTL, as well as a substantial number of QTL at single chromosomal locations. To extract the most information from these published records, a meta-analysis was conducted to obtain consensus on QTL location and allelic substitution effect of these QTL. This required modification and development of statistical methodologies. The meta-analysis indicated a number of consensus regions, the most striking being two distinct regions affecting milk yield on chromosome 6 at 49 cM and 87 cM explaining 4.2 and 3.6 percent of the genetic variance of milk yield, respectively. The first of these regions (near marker BM143) affects five separate milk production traits (protein yield, protein percent, fat yield, fat percent, as well as milk yield).     

Deficiency mapping of quantitative trait loci affecting longevity in Drosophila melanogaster

In a previous study, sex-specific quantitative trait loci (QTL) affecting adult longevity were mapped by linkage to polymorphic roo transposable element markers, in a population of recombinant inbred lines derived from the Oregon and 2b strains of Drosophila melanogaster. Two life span QTL were each located on chromosomes 2 and 3, within sections 33E-46C and 65D-85F on the cytological map, respectively. We used quantitative deficiency complementation mapping to further resolve the locations of life span QTL within these regions. The Oregon and 2b strains were each crossed to 47 deficiencies spanning cytological regions 32F-44E and 64C-76B, and quantitative failure of the QTL alleles to complement the deficiencies was assessed. We initially detected a minimum of five and four QTL in the chromosome 2 and 3 regions, respectively, illustrating that multiple linked factors contribute to each QTL detected by recombination mapping. The QTL locations inferred from deficiency mapping did not generally correspond to those of candidate genes affecting oxidative and thermal stress or glucose metabolism. The chromosome 2 QTL in the 35B-E region was further resolved to a minimum of three tightly linked QTL, containing six genetically defined loci, 24 genes, and predicted genes that are positional candidates corresponding to life span QTL. This region was also associated with quantitative variation in life span in a sample of 10 genotypes collected from nature. Quantitative deficiency complementation is an efficient method for fine-scale QTL mapping in Drosophila and can be further improved by controlling the background genotype of the strains to be tested.     

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.     

Multipoint identity-by-descent prediction using dense markers to map quantitative trait loci and estimate effective population size

A novel multipoint method, based on an approximate coalescence approach, to analyze multiple linked markers is presented. Unlike other approximate coalescence methods, it considers all markers simultaneously but only two haplotypes at a time. We demonstrate the use of this method for linkage disequilibrium (LD) mapping of QTL and estimation of effective population size. The method estimates identity-by-descent (IBD) probabilities between pairs of marker haplotypes. Both LD and combined linkage and LD mapping rely on such IBD probabilities. The method is approximate in that it considers only the information on a pair of haplotypes, whereas a full modeling of the coalescence process would simultaneously consider all haplotypes. However, full coalescence modeling is computationally feasible only for few linked markers. Using simulations of the coalescence process, the method is shown to give almost unbiased estimates of the effective population size. Compared to direct marker and haplotype association analyses, IBD-based QTL mapping showed clearly a higher power to detect a QTL and a more realistic confidence interval for its position. The modeling of LD could be extended to estimate other LD-related parameters such as recombination rates.     

Marker-based mapping of quantitative trait loci using replicated progenies

Summary When heritability of the trait under investigation is low, replicated progenies can bring about a major reduction in the number of individuals that need to be scored for marker genotype in determining linkage between marker loci and quantitative trait loci (QTL). Savings are greatest when heritability of the trait is low, but are much reduced when heritability of the quantitative trait is moderate to high. Required numbers for recombinant inbred lines will be greater than those required for a simple F₂ population when heritabilities are moderate to high and the proportion of recombination between marker locus and quantitative trait locus is substantial.Contribution No. 2613-E of the Agricultural Research Organization, 1989 series     

Mapping quantitative trait loci affecting female reproductive traits on porcine chromosome 8

An understanding of the genetic control of porcine female reproductive performance would offer the opportunity to utilize natural variation and improve selective breeding programs through marker-assisted selection. The Chinese Meishan is one of the most prolific pig breeds known, farrowing three to five more viable piglets per litter than the European Large White breed. This difference in prolificacy is attributed to the Meishan's superior prenatal survival levels. The present study utilized a three-generation cross in which the founder grandparental animals were purebred Meishan and Large White pigs in a scan for quantitative trait loci (QTL) on porcine chromosome 8 (SSC8) associated with reproductive performance. Reproductive traits, including number of corpora lutea (ovulation rate), teat number, litter size, and prenatal survival, were recorded for as many as 220 F2 females. Putative QTL for the related traits of litter size and prenatal survival were identified at the distal end of the long arm of SSC8. A physiological candidate gene, SPP1, was found to lie within the 95% confidence interval of these QTL. A suggestive QTL for teat number was revealed on the short arm of SSC8. The present study demonstrates, to our knowledge, the first independent confirmation of QTL for fecundity on SSC8, and these QTL regions provide a crucial starting point in the search for the causal genetic variants.     

Fine mapping of quantitative trait loci for mastitis resistance on bovine chromosome 11

Quantitative trait loci (QTL) affecting clinical mastitis (CM) and somatic cell score (SCS) were mapped on bovine chromosome 11. The mapping population consisted of 14 grandsire families belonging to three Nordic red cattle breeds: Finnish Ayrshire (FA), Swedish Red and White (SRB) and Danish Red. The families had previously been shown to segregate for udder health QTL. A total of 524 progeny tested bulls were included in the analysis. A linkage map including 33 microsatellite and five SNP markers was constructed. We performed combined linkage disequilibrium and linkage analysis (LDLA) using the whole data set. Further analyses were performed for FA and SRB separately to study the origin of the identified QTL/haplotype and to examine if it was common in both populations. Finally, different two-trait models were fitted. These postulated either a pleiotropic QTL affecting both traits; two linked QTL, each affecting one trait; or one QTL affecting a single trait. A QTL affecting CM was fine-mapped. In FA, a haplotype having a strong association with a high negative effect on mastitis resistance was identified. The mapping precision of an earlier detected SCS-QTL was not improved by the LDLA analysis because of lack of linkage disequilibrium between the markers used and the QTL in the region.     

Quantitative trait loci affecting clinical mastitis and somatic cell count in dairy cattle

Norway has a field recording system for dairy cattle that includes recording of all veterinary treatments on an individual animal basis from 1978 onwards. Application of these data in a genome search for quantitative trait loci (QTL) verified genome-wise significant QTL affecting clinical mastitis on Chromosome (Chr) 6. Additional putative QTL for clinical mastitis were localized to Chrs. 3, 4, 14, and 27. The comprehensive field recording system includes information on somatic cell count as well. This trait is often used in selection against mastitis when direct information on clinical mastitis is not available. The absence of common QTL positions for the two traits in our study indicates that the use of somatic cell count data in QTL studies aimed for reducing the incidence of mastitis should be carefully evaluated. Received: 9 May 2001 / Accepted: 6 July 2001     

Fine mapping of collagen-induced arthritis quantitative trait loci in an advanced intercross line

The generation of advanced intercross lines (AIL) is a powerful approach for high-resolution fine mapping of quantitative trait loci (QTLs), because they accumulate much more recombination events compared with conventional F2 intercross and N2 backcross. However, the application of this approach is severely hampered by the requirements of excessive resources to maintain such crosses, i.e., in terms of animal care, space, and time. Therefore, in this study, we produced an AIL to fine map collagen-induced arthritis (CIA) QTLs using comparatively limited resources. We used only 308 (DBA/1 x FVB/N)F11/12 AIL mice to refine QTLs controlling the severity and onset of arthritis as well as the Ab response and T cell subset in CIA, namely Cia2, Cia27, and Trmq3. These QTLs were originally identified in (DBA/1 x FVB/N)F2 progeny. The confidence intervals of the three QTLs were refined from 40, 43, and 48 Mb to 12, 4.1, and 12 Mb, respectively. The data were complemented by the use of another QTL fine-mapping approach, haplotype analysis, to further refine Cia2 into a 2-Mb genomic region. To aid in the search for candidate genes for the QTLs, genome-wide expression profiling was performed to identify strain-specific differentially expressed genes within the confidence intervals. Of the 1396 strain-specific differentially expressed genes, 3, 3, and 12 genes were within the support intervals of the Cia2, Cia27, and Trmq3, respectively. In addition, this study revealed that Cia27 and Trmq3 controlling anti-CII IgG2a Ab and CD4:CD8 T cell ratio, respectively, also regulated CIA clinical phenotypes.     

High-resolution mapping of quantitative trait loci affecting increased life span in Drosophila melanogaster

Limited life span and senescence are near-universal characteristics of eukaryotic organisms, controlled by many interacting quantitative trait loci (QTL) with individually small effects, whose expression is sensitive to the environment. Analyses of mutations in model organisms have shown that genes affecting stress resistance and metabolism affect life span across diverse taxa. However, there is considerable segregating variation for life span in nature, and relatively little is known about the genetic basis of this variation. Replicated lines of Drosophila that have evolved increased longevity as a correlated response to selection for postponed senescence are valuable resources for identifying QTL affecting naturally occurring variation in life span. Here, we used deficiency complementation mapping to identify at least 11 QTL on chromosome 3 that affect variation in life span between five old (O) lines selected for postponed senescence and their five base (B) population control lines. Most QTL were sex specific, and all but one affected multiple O lines. The latter observation is consistent with alleles at intermediate frequency in the base population contributing to the response to selection for postponed senescence. The QTL were mapped with high resolution and contained from 12 to 170 positional candidate genes.     

Joint mapping of quantitative trait loci using F2 populations

In this paper, the theory of joint mapping of quantitative trait loci is extended to F₂ populations. Two independent regression equations, related to the additive and dominance effects respectively, are derived. Therefore, there are three alternative strategies for mapping QTLs, called additive-based mapping (ABM), dominance-based mapping (DBM) and additive-dominance-based mapping (ADBM). Simulation results have shown that ADBM is the most appropriate in most situations.