Confirmation and refinement of a QTL on BTA5 affecting milk production traits in the Fleckvieh dual purpose cattle breed

We analysed a QTL affecting milk yield (MY), milk protein yield (PY) and milk fat yield (FY) in the dual purpose cattle breed Fleckvieh on BTA5. Twenty-six microsatellite markers covering 135 cM were selected to analyse nine half-sib families containing 605 sons in a granddaughter design. We thereby assigned two new markers to the public linkage map using the CRI-MAP program. Phenotypic records were daughter yield deviations (DYD) originating from the routinely performed genetic evaluations of breeding animals. To determine the position of the QTL, three different approaches were applied: interval mapping (IM), linkage analysis by variance component analysis (LAVC), and combined linkage disequilibrium (LD) and linkage (LDL) analysis. All three methods mapped the QTL in the same marker interval ( BM2830-ETH152 ) with the greatest test-statistic value at 118, 119.33 and 119.33 cM respectively. The positive QTL allele simultaneously increases DYD in the first lactation by 272 kg milk, 7.1 kg milk protein and 7.0 kg milk fat. Although the mapping accuracy and the significance of a QTL effect increased from IM over LAVC to LDL, the confidence interval was large (13, 20 and 24 cM for FY, MY and PY respectively) for the positional cloning of the causal gene. The estimated averages of pair wise marker LD with a distance <5 cM were low (0.107) and reflect the large effective population size of the Fleckvieh subpopulation analysed. This low level of LD suggests a need for increase in marker density in following fine mapping steps.     

Identification of SNPs in Interferon Gamma,Interleukin-22, and Their Receptors and Associations with Health and Production-Related Traits in Canadian Holstein Bulls

Genetic variants in a number of immunoregulatory genes have been previously associated with health and production traits in dairy cattle. Therefore, in the following study, the genes coding interferon gamma (IFNG), IFNG receptor 1 and 2 domains, interleukin-22 (IL22), and IL22 receptor alpha 1, were investigated for single nucleotide polymorphisms (SNPs) in Holstein bulls. These SNPs, along with SNPs previously identified in IL10, IL10 receptor, and transforming growth factor beta 1 (TGFB1) genes, were evaluated for statistical associations to estimated breeding values for milk somatic cell score (SCS), a trait highly correlated to mastitis incidence, and various production-related traits, including milk yield, protein yield, fat yield, and lactation persistency. While no significant associations were found between these SNPs and SCS, SNPs in IL10 receptor beta subunit showed a significant effect on protein yield and lactation persistency. While there is evidence that IL10 plays an important role during lactation, it is also likely that the effects of SNPs in IL10 receptor beta subunit on protein yield and lactation persistency are due to linkage disequilibrium with a neighboring QTL.     

Identification of SNPs in interferon gamma, interleukin-22, and their receptors and associations with health and production-related traits in Canadian Holstein bulls

Genetic variants in a number of immunoregulatory genes have been previously associated with health and production traits in dairy cattle. Therefore, in the following study, the genes coding interferon gamma (IFNG), IFNG receptor 1 and 2 domains, interleukin-22 (IL22), and IL22 receptor alpha 1, were investigated for single nucleotide polymorphisms (SNPs) in Holstein bulls. These SNPs, along with SNPs previously identified in IL10, IL10 receptor, and transforming growth factor beta 1 (TGFB1) genes, were evaluated for statistical associations to estimated breeding values for milk somatic cell score (SCS), a trait highly correlated to mastitis incidence, and various production-related traits, including milk yield, protein yield, fat yield, and lactation persistency. While no significant associations were found between these SNPs and SCS, SNPs in IL10 receptor beta subunit showed a significant effect on protein yield and lactation persistency. While there is evidence that IL10 plays an important role during lactation, it is also likely that the effects of SNPs in IL10 receptor beta subunit on protein yield and lactation persistency are due to linkage disequilibrium with a neighboring QTL.     

GWA Analysis for Milk Production Traits in Dairy Sheep and Genetic Support for a QTN Influencing Milk Protein Percentage in the LALBA Gene

In this study, we used the Illumina OvineSNP50 BeadChip to conduct a genome-wide association (GWA) analysis for milk production traits in dairy sheep by analyzing a commercial population of Spanish Churra sheep. The studied population consisted of a total of 1,681 Churra ewes belonging to 16 half-sib families with available records for milk yield (MY), milk protein and fat yields (PY and FY) and milk protein and fat contents (PP and FP). The most significant association identified reached experiment-wise significance for PP and FP and was located on chromosome 3 (OAR3). These results confirm the population-level segregation of a previously reported QTL affecting PP and suggest that this QTL has a significant pleiotropic effect on FP. Further associations were detected at the chromosome-wise significance level on 14 other chromosomal regions. The marker on OAR3 showing the highest significant association was located at the third intron of the alpha-lactalbumin (LALBA) gene, which is a functional and positional candidate underlying this association. Sequencing this gene in the 16 Churra rams of the studied resource population identified additional polymorphisms. One out of the 31 polymorphisms identified was located within the coding gene sequence (LALBA_g.242T>C) and was predicted to cause an amino acid change in the protein (Val27Ala). Different approaches, including GWA analysis, a combined linkage and linkage disequilibrium study and a concordance test with the QTL segregating status of the sires, were utilized to assess the role of this mutation as a putative QTN for the genetic effects detected on OAR3. Our results strongly support the polymorphism LALBA_g.242T>C as the most likely causal mutation of the studied OAR3 QTL affecting PP and FP, although we cannot rule out the possibility that this SNP is in perfect linkage disequilibrium with the true causal polymorphism.     

Association of the heart fatty acid-binding protein (FABP3) gene with milk traits in Manchega breed sheep

The ovine fatty acid-binding protein type 3 gene has been chosen as a functional candidate gene for milk traits. Two different single nucleotide polymorphisms (SNPs) of ovine FABP3 gene have been tested in a daughter design comprising 13 families. No association was found between estimated breeding values for milk yield, protein and fat contents (FC) and genotypes across families using anova and transmission disequilibrium test (TDT). In within-family analysis, one family showed a significant effect for FC. These results could indicate linkage disequilibrium between the FABP3 gene and a quantitative trait loci (QTL) for FC, with the heterozygous genotype associated with a positive effect in this trait.     

Mapping of a milk production quantitative trait locus to a 1.056 Mb region on bovine chromosome 5 in the Fleckvieh dual purpose cattle breed

Background

In a previous study in the Fleckvieh dual purpose cattle breed, we mapped a quantitative trait locus (QTL) affecting milk yield (MY1), milk protein yield (PY1) and milk fat yield (FY1) during first lactation to the distal part of bovine chromosome 5 (BTA5), but the confidence interval was too large for positional cloning of the causal gene. Our objective here was to refine the position of this QTL and to define the candidate region for high-throughput sequencing.

Methods

In addition to those previously studied, new Fleckvieh families were genotyped, in order to increase the number of recombination events. Twelve new microsatellites and 240 SNP markers covering the most likely QTL region on BTA5 were analysed. Based on haplotype analysis performed in this complex pedigree, families segregating for the low frequency allele of this QTL (minor allele) were selected. Single- and multiple-QTL analyses using combined linkage and linkage disequilibrium methods were performed.

Results

Single nucleotide polymorphism haplotype analyses on representative family sires and their ancestors revealed that the haplotype carrying the minor QTL allele is rare and most probably originates from a unique ancestor in the mapping population. Analyses of different subsets of families, created according to the results of haplotype analysis and availability of SNP and microsatellite data, refined the previously detected QTL affecting MY1 and PY1 to a region ranging from 117.962 Mb to 119.018 Mb (1.056 Mb) on BTA5. However, the possibility of a second QTL affecting only PY1 at 122.115 Mb was not ruled out.

Conclusion

This study demonstrates that targeting families segregating for a less frequent QTL allele is a useful method. It improves the mapping resolution of the QTL, which is due to the division of the mapping population based on the results of the haplotype analysis and to the increased frequency of the minor allele in the families. Consequently, we succeeded in refining the region containing the previously detected QTL to 1 Mb on BTA5. This candidate region contains 27 genes with unknown or partially known function(s) and is small enough for high-throughput sequencing, which will allow future detailed analyses of candidate genes.     

Extent and structure of linkage disequilibrium in canola quality winter rapeseed (Brassica napus L.)

Linkage disequilibrium was investigated in canola quality winter rapeseed to analyze (1) the prospects for whole-genome association analyses and (2) the impact of the recent breeding history of rapeseed on linkage disequilibrium. A total of 845 mapped AFLP markers with allele frequencies ≥0.1 were used for the analysis of linkage disequilibrium in a population of 85 canola quality winter rapeseed genotypes. A low overall level of linkage disequilibrium was found with a mean r ² of only 0.027 over all 356,590 possible marker pairs. At a significance threshold of P = 2.8 × 10⁻⁷, which was derived by a Bonferroni correction from a global α-level of 0.1, only 0.78% of the marker pairs were in significant linkage disequilibrium. Among physically linked marker pairs, the level of linkage disequilibrium was about five times higher with more than 10% of marker pairs in significant linkage disequilibrium. Linkage disequilibrium decayed rapidly with distance between linked markers with high levels of linkage disequilibrium extending only for about 2 cM. Owing to the rapid decay of linkage disequilibrium with distance association analyses in canola quality rapeseed will have a significantly higher resolution than QTL analyses in segregating populations by interval mapping, but much larger number of markers will be necessary to cover the whole genome. A major impact of the recent breeding history of rapeseed on linkage disequilibrium could not be observed.     

Simultaneous mining of linkage and linkage disequilibrium to fine map quantitative trait loci in outbred half-sib pedigrees: revisiting the location of a quantitative trait locus with major effect on milk production on bovine chromosome 14

A maximum-likelihood QTL mapping method that simultaneously exploits linkage and linkage disequilibrium and that is applicable in outbred half-sib pedigrees is described. The method is applied to fine map a QTL with major effect on milk fat content in a 3-cM marker interval on proximal BTA14. This proximal location is confirmed by applying a haplotype-based association method referred to as recombinant ancestral haplotype analysis. The origin of the discrepancy between the QTL position derived in this work and that of a previous analysis is examined and shown to be due to the existence of distinct marker haplotypes associated with QTL alleles having large substitution effects.     

Joint linkage and linkage disequilibrium mapping of quantitative trait loci in natural populations

Linkage analysis and allelic association (also referred to as linkage disequilibrium) studies are two major approaches for mapping genes that control simple or complex traits in plants, animals, and humans. But these two approaches have limited utility when used alone, because they use only part of the information that is available for a mapping population. More recently, a new mapping strategy has been designed to integrate the advantages of linkage analysis and linkage disequilibrium analysis for genome mapping in outcrossing populations. The new strategy makes use of a random sample from a panmictic population and the open-pollinated progeny of the sample. In this article, we extend the new strategy to map quantitative trait loci (QTL), using molecular markers within the EM-implemented maximum-likelihood framework. The most significant advantage of this extension is that both linkage and linkage disequilibrium between a marker and QTL can be estimated simultaneously, thus increasing the efficiency and effectiveness of genome mapping for recalcitrant outcrossing species. Simulation studies are performed to test the statistical properties of the MLEs of genetic and genomic parameters including QTL allele frequency, QTL effects, QTL position, and the linkage disequilibrium of the QTL and a marker. The potential utility of our mapping strategy is discussed.     

Replication and refinement of a quantitative trait locus influencing milk protein percentage on ovine chromosome 3

A previous genome scan that was conducted in Spanish Churra sheep identified a significant quantitative trait locus (QTL) for milk protein percentage (PP) on chromosome 3 (OAR3), between markers KD103 and OARVH34. The aim of this study was to replicate these results and to refine the mapped position of this QTL. To accomplish this goal, we analysed 14 new half‐sib families of Spanish Churra sheep including 1661 ewes from 29 different flocks. These animals were genotyped for 21 microsatellite markers mapping to OAR3. In addition to a classical linkage analysis (LA), a combined linkage disequilibrium and linkage analysis (LDLA) was performed with the aim of enhancing the resolution of the QTL mapping. The LA that was performed in this sheep population identified the presence of a highly significant QTL for PP near marker KD103 (P_c < 0.001; P_exp < 0.001). The phenotypic variance that was owing to the QTL was 2.74%. Two segregating families for the target QTL were identified in this population with QTL effect estimates of 0.47 and 0.95 SD. The LDLA identified the same QTL as the previous analyses with a high level of statistical significance (P = 9.184 E‐11) and narrowed the confidence interval (CI) to a 13 cM region. These results confirm the segregation of the previously identified OAR3 QTL that influences PP in Spanish Churra sheep. Future research will aim to increase the marker density across the refined CI and to analyse the corresponding candidate genes to identify the allelic variant or variants that underlie this genetic effect.     

Molecular dissection of a quantitative trait locus: a phenylalanine-to-tyrosine substitution in the transmembrane domain of the bovine growth hormone receptor is associated with a major effect on milk yield and composition

We herein report on our efforts to improve the mapping resolution of a QTL with major effect on milk yield and composition that was previously mapped to bovine chromosome 20. By using a denser chromosome 20 marker map and by exploiting linkage disequilibrium using two distinct approaches, we provide strong evidence that a chromosome segment including the gene coding for the growth hormone receptor accounts for at least part of the chromosome 20 QTL effect. By sequencing individuals with known QTL genotype, we identify an F to Y substitution in the transmembrane domain of the growth hormone receptor gene that is associated with a strong effect on milk yield and composition in the general population.     

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.     

Linkage disequilibrium fine mapping of quantitative trait loci: A simulation study

Recently, the use of linkage disequilibrium (LD) to locate genes which affect quantitative traits (QTL) has received an increasing interest, but the plausibility of fine mapping using linkage disequilibrium techniques for QTL has not been well studied. The main objectives of this work were to (1) measure the extent and pattern of LD between a putative QTL and nearby markers in finite populations and (2) investigate the usefulness of LD in fine mapping QTL in simulated populations using a dense map of multiallelic or biallelic marker loci. The test of association between a marker and QTL and the power of the test were calculated based on single-marker regression analysis. The results show the presence of substantial linkage disequilibrium with closely linked marker loci after 100 to 200 generations of random mating. Although the power to test the association with a frequent QTL of large effect was satisfactory, the power was low for the QTL with a small effect and/or low frequency. More powerful, multi-locus methods may be required to map low frequent QTL with small genetic effects, as well as combining both linkage and linkage disequilibrium information. The results also showed that multiallelic markers are more useful than biallelic markers to detect linkage disequilibrium and association at an equal distance.     

Mapping quantitative trait loci (QTL) in sheep. IV. Analysis of lactation persistency and extended lactation traits in sheep

Background

In sheep dairy production, total lactation performance, and length of lactation of lactation are of economic significance. A more persistent lactation has been associated with improved udder health. An extended lactation is defined by a longer period of milkability. This study is the first investigation to examine the presence of quantitative trait loci (QTL) for extended lactation and lactation persistency in sheep.

Methods

An (Awassi × Merino) × Merino single-sire backcross family with 172 ewes was used to map QTL for lactation persistency and extended lactation traits on a framework map of 189 loci across all autosomes. The Wood model was fitted to data from multiple lactations to estimate parameters of ovine lactation curves, and these estimates were used to derive measures of lactation persistency and extended lactation traits of milk, protein, fat, lactose, useful yield, and somatic cell score. These derived traits were subjected to QTL analyses using maximum likelihood estimation and regression analysis.

Results

Overall, one highly significant (LOD > 3.0), four significant (2.0 < LOD < 3.0) and five suggestive (1.7 < LOD < 2.0) QTL were detected across all traits in common by both mapping methods. One additional suggestive QTL was identified using maximum likelihood estimation, and four suggestive (0.01 < P < 0.05) and two significant (P < 0.01) QTL using the regression approach only. All detected QTL had effect sizes in the range of 0.48 to 0.64 SD, corresponding to QTL heritabilities of 3.1 to 8.9%. The comparison of the detected QTL with results in cattle showed conserved linkage regions. Most of the QTL identified for lactation persistency and extended lactation did not coincide. This suggests that persistency and extended lactation for the same as well as different milk yield and component traits are not controlled by the same genes.

Conclusion

This study identified ten novel QTL for lactation persistency and extended lactation in sheep, but results suggest that lactation persistency and extended lactation do not have a major gene in common. These results provide a basis for further validation in extended families and other breeds as well as targeting regions for genome-wide association mapping using high-density SNP arrays.     

Mapping and Analysis of Dairy Cattle Quantitative Trait Loci by Maximum Likelihood Methodology Using Milk Protein Genes as Genetic Markers

Maximum likelihood methodology was used to estimate effects of both a marker gene and a linked quantitative trait locus (QTL) on quantitative traits in a segregating population. Two alleles were assumed for the QTL. In addition to the effects of genotypes at both loci on the mean of the quantitative trait, recombination frequency between the loci, frequency of the QTL alleles and the residual standard deviation were also estimated. Thus six parameters were estimated in addition to the marker genotype means. The statistical model was tested on simulated data, and used to estimate direct and linked effects of the milk protein genes, β-lactoglobulin, κcasein, and β-casein, on milk, fat, and protein production and fat and protein percent in the Dutch dairy cattle population. β-Lactoglobulin had significant direct effects on milk yield and fat percent. κ-Casein had significant direct effects on milk yield, protein percent and fat yield. β-Casein had significant direct effects on milk yield, fat and protein percent and fat and protein yield. Linked QTL with significant effects on fat percent were found for κ-casein and β-casein. Since the β-casein and κ-casein genes are closely linked, it is likely that the same QTL was detected for those two markers. Further, a QTL with a significant effect on fat yield was found to be linked to κ-casein and a QTL with a significant effect on protein yield was linked to β-lactoglobulin.     

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.     

Detection of loci affecting milk production and health traits in an elite US Holstein population using microsatellite markers

Quantitative trait loci (QTL) affecting health and milk production traits were studied in seven large half-sib US Holstein families by using the granddaughter design. Genotyping for 16 markers was completed and marker allele differences within and pooled-across families were analysed. Potential QTL were identified for somatic cell score (SCS), fat yield, fat percentage, protein yield and protein percentage. Three markers (BM203, BM4505 and BM2078) were associated with significant effects for different traits and, after further analysis, may be useful in marker-assisted selection in specific families. Comparisons between these data and previously identified QTL support the location of a QTL for milk yield and protein yield on chromosome 21.     

Fine-mapping QTL for mastitis resistance on BTA9 in three Nordic red cattle breeds

A QTL affecting clinical mastitis and/or somatic cell score (SCS) has been reported previously on chromosome 9 from studies in 16 families from the Swedish Red and White (SRB), Finnish Ayrshire (FA) and Danish Red (DR) breeds. In order to refine the QTL location, 67 markers were genotyped over the whole chromosome in the 16 original families and 18 additional half-sib families. This enabled linkage disequilibrium information to be used in the analysis. Data were analysed by an approach that combines information from linkage and linkage disequilibrium, which allowed the QTL affecting clinical mastitis to be mapped to a small interval (<1 cM) between the markers BM4208 and INRA084 . This QTL showed a pleiotropic effect on SCS in the DR and SRB breeds. Haplotypes associated with variations in mastitis resistance were identified. The haplotypes were predictive in the general population and can be used in marker-assisted selection. Pleiotropic effects of the mastitis QTL were studied for three milk production traits and eight udder conformation traits. This QTL was also associated with yield traits in DR but not in FA or SRB. No QTL were found for udder conformation traits on chromosome 9.     

人类复杂遗传疾病QTL分析方法的理论探讨

利用连锁不平衡理论,人类遗传学家已能把影响人类疾病的质量基因定位在小至1cM区域内,有些基因已被克隆出来。罗泽伟等进一步发展统计分析方法检测及估算分子标记与QTL之间的连锁不平衡系数,从而提出了人类复杂遗传病高解析度基因定位的理论策略。以此为基础,进一步探讨了供试群体在双亲基因频率存在差异时检测QTL和检测QTL互作的方法,给出了有关的理论结果。     

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.