Detection of quantitative trait loci affecting milk production traits on 10 chromosomes in Holstein cattle

Defining Candidate Imprinted loci in Bos taurus

Using a whole-genome assembly of Bos taurus, I applied my bioinformatics strategy to locate candidate imprinting control regions (ICRs) genome-wide. In mammals, genomic imprinting plays essential roles in embryogenesis. In my strategy, peaks in plots mark the locations of known, inferred, and candidate ICRs. Genes in the vicinity of candidate ICRs correspond to potential imprinted genes. By displaying my datasets on the UCSC genome browser, one could view peak positions with respect to genomic landmarks. I give two examples of candidate ICRs in loci that influence spermatogenesis in bulls: CNNM1 and CNR1. I also give examples of candidate ICRs in loci that influence muscle development: SIX1 and BCL6. By examining the ENCODE data reported for mice, I deduced regulatory clues about cattle. I focused on DNase I hypersensitive sites (DHSs). Such sites reveal accessibility of chromatin to regulators of gene expression. For inspection, I chose DHSs in chromatin from mouse embryonic stem cells (ESCs) ES-E14, mesoderm, brain, heart, and skeletal muscle. The ENCODE data revealed that the SIX1 promoter was accessible to the transcription initiation apparatus in mouse ESCs, mesoderm, and skeletal muscles. The data also revealed accessibility of BCL6 locus to regulatory proteins in mouse ESCs and examined tissues.

Direct Phenotyping and Principal Component Analysis of Type Traits Implicate Novel QTL in Bovine Mastitis through Genome-Wide Association

Simple Summary

It is well established that the physical conformation of a cow’s udder and teats may influence her susceptibility to mastitis, an inflammatory condition of the udder, which has 25% prevalence in the United States. Our aim was to improve the biological understanding of the genetics underlying mastitis by intensively characterizing cows for udder and teat conformation, including the novel traits of teat width and end shape, and directly associating those phenotypes with high-density genotypes for those exact same cows. We also generated a composite measure that accounts for multiple high-mastitis-risk udder and teat conformations in a single index for risk phenotypes. Using this approach, we identified novel genetic markers associated with udder and teat conformation, which may be good candidates for inclusion in national genetic evaluations for selection of mastitis-resistant cows. Mastitis is the costliest disease facing US dairy producers, and integrating genetic information regarding disease susceptibility into breeding programs may be an efficient way to mitigate economic loss, support the judicious use of antimicrobials, and improve animal welfare.

Abstract

Our objectives were to robustly characterize a cohort of Holstein cows for udder and teat type traits and perform high-density genome-wide association studies for those traits within the same group of animals, thereby improving the accuracy of the phenotypic measurements and genomic association study. Additionally, we sought to identify a novel udder and teat trait composite risk index to determine loci with potential pleiotropic effects related to mastitis. This approach was aimed at improving the biological understanding of the genetic factors influencing mastitis. Cows (N = 471) were genotyped on the Illumina BovineHD777k beadchip and scored for front and rear teat length, width, end shape, and placement; fore udder attachment; udder cleft; udder depth; rear udder height; and rear udder width. We used principal component analysis to create a single composite measure describing type traits previously linked to high odds of developing mastitis within our cohort of cows. Genome-wide associations were performed, and 28 genomic regions were significantly associated (Bonferroni-corrected p < 0.05). Interrogation of these genomic regions revealed a number of biologically plausible genes whicht may contribute to the development of mastitis and whose functions range from regulating cell proliferation to immune system signaling, including ZNF683, DHX9, CUX1, TNNT1, and SPRY1. Genetic investigation of the risk composite trait implicated a novel locus and candidate genes that have potentially pleiotropic effects related to mastitis.

Integrative analysis of genome-wide DNA methylation and gene expression profiles reveals important epigenetic genes related to milk production traits in dairy cattle

Epigenetic modification plays a critical role in establishing and maintaining cell differentiation, embryo development, tumorigenesis and many complex diseases. However, little is known about the epigenetic regulatory mechanisms for milk production in dairy cattle. Here, we conducted an epigenome-wide study, together with gene expression profiles to identify important epigenetic candidate genes related to the milk production traits in dairy cattle. Whole-genome bisulphite sequencing and RNA sequencing were employed to detect differentially methylated genes (DMG) and differentially expressed genes (DEG) in blood samples in dry period and lactation period between two groups of cows with extremely high and low milk production performance. A total of 10,877 and 6,617 differentially methylated regions were identified between the two groups in the two periods, which corresponded to 3,601 and 2,802 DMGs, respectively. Furthermore, 156 DEGs overlap with DMGs in comparison of the two groups, and 131 DEGs overlap with DMGs in comparison of the two periods. By integrating methylome, transcriptome and GWAS data, some potential candidate genes for milk production traits in dairy cattle were suggested, such as DOCK1, PTK2 and PIK3R1. Our studies may contribute to a better understanding of epigenetic modification on milk production traits of dairy cattle.

Identification of single nucleotide polymorphisms of PIK3R1 and DUSP1 genes and their genetic associations with milk production traits in dairy cows

Background

Previously, phosphoinositide-3-kinase regulatory subunit 1 (PIK3R1) and dual specificity phosphatase 1 (DUSP1) were identified as promising candidate genes for milk production traits due to their being differentially expressed between the dry period and the peak of lactation in livers of dairy cows. Hence, in this study, the single nucleotide polymorphisms (SNPs) of PIK3R1 and DUSP1 genes were identified and their genetic associations with milk yield, fat yield, fat percentage, protein yield, and protein percentage, were investigated using 1067 Chinese Holstein cows from 40 sire families.

Results

By re-sequencing the entire coding region and 2000 bp of the 5′ and 3′ flanking regions of the two genes, one SNP in the 5′ untranslated region (UTR), three in the 3′ UTR, and two in the 3′ flanking region of PIK3R1 were identified, and one in the 5′ flanking region, one in the 3′ UTR, and two in the 3′ flanking region of DUSP1 were found. Subsequent single-locus association analyses showed that five SNPs in PIK3R1, rs42590258, rs210389799, rs208819656, rs41255622, rs133655926, and rs211408208, and four SNPs in DUSP1, rs207593520, rs208460068, rs209154772, and rs210000760, were significantly associated with milk, fat and protein yields in the first or second lactation (P values ≤ 0.0001 and 0.0461). In addition, by the Haploview 4.2 software, the six and four SNPs in PIK3R1 and DUSP1 respectively formed one haplotype block, and the haplotype-based association analyses showed significant associations between their haplotype combinations and the milk traits in both two lactations (P values ≤ 0.0001 and 0.0364). One SNP, rs207593520(T/G), was predicted to alter the transcription factor binding sites (TFBSs) in the 5′ flanking region of DUSP1. Further, the dual-luciferase assay showed that the transcription activity of allele T in rs207593520 was significantly higher than that of allele G, suggesting the activation of transcriptional activity of DUSP1 gene by allele T of rs207593520. Thus, the rs207593520 SNP was highlighted as a potential causal mutation that should be further verified.

Conclusions

We demonstrated novel and significant genetic effects of the PIK3R1 and DUSP1 genes on milk production traits in dairy cows, and our findings provide information for use in dairy cattle breeding.

SNV discovery and functional candidate gene identification for milk composition based on whole genome resequencing of Holstein bulls with extremely high and low breeding values

We have sequenced the whole genomes of eight proven Holstein bulls from the four half-sib or full-sib families with extremely high and low estimated breeding values (EBV) for milk protein percentage (PP) and fat percentage (FP) using Illumina re-sequencing technology. Consequently, 2.3 billion raw reads were obtained with an average effective depth of 8.1×. After single nucleotide variant (SNV) calling, total 10,961,243 SNVs were identified, and 57,451 of them showed opposite fixed sites between the bulls with high and low EBVs within each family (called as common differential SNVs). Next, we annotated the common differential SNVs based on the bovine reference genome, and observed that 45,188 SNVs (78.70%) were located in the intergenic region of genes and merely 11,871 SNVs (20.67%) located within the protein-coding genes. Of them, 13,099 common differential SNVs that were within or close to protein-coding genes with less than 5 kb were chosen for identification of candidate genes for milk compositions in dairy cattle. By integrated analysis of the 2,657 genes with the GO terms and pathways related to protein and fat metabolism, and the known quantitative trait loci (QTLs) for milk protein and fat traits, we identified 17 promising candidate genes: ALG14, ATP2C1, PLD1, C3H1orf85, SNX7, MTHFD2L, CDKN2D, COL5A3, FDX1L, PIN1, FIG4, EXOC7, LASP1, PGS1, SAO, GPLD1 and MGEA5. Our findings provided an important foundation for further study and a prompt for molecular breeding of dairy cattle.

Identification of quantitative trait loci for milk protein percentage in Murrah buffaloes

Milk protein is an important constituent of milk in buffaloes and is moderately heritable. The milk protein percentage varies significantly between breeds/herds/species. Buffaloes can be selected for higher milk protein percentage and this paper provides QTLs for marker assisted selection in buffaloes. The milk protein percentage records on 2,028 daughters belonging to 12 half sib families were analyzed for the identification of QTLs on 8 chromosomes in buffaloes using chromosome scans. The single marker analysis revealed 74 markers to be associated with milk protein percentage in 12 sire families. When common markers were removed from the analysis, 51 markers remained. The Interval mapping using R/qtl identified 69 QTLs in 12 half sib families on 8 chromosomes of buffalo. The meta QTL analysis defined 25 consensus QTL regions in buffaloes for milk protein percentage. Most of the QTLs identified have been reported for cattle however few new chromosomal locations were also identified to be associated with milk protein percentage in buffaloes. Comparative genomics revealed 1117 genes underlying the QTL regions associated with milk protein percentage. Among these, 109 genes were directly associated with protein metabolism. The protein-protein interaction among the genes and gene ontology analysis and pathways have been identified. These 109 genes have potential to be candidate genes for milk protein percentage in buffaloes.

Identification of quantitative trait loci for milk yield in Murrah buffaloes

A reference family consisting of 12 half sib sire families were created for the identification of QTLs for milk yield in buffaloes. Daughters were recorded for monthly test day milk yield. The number of daughters per sire varied from 50 to 335 daughters per sire. Seventy nine polymorphic microsatellite markers located on 8 chromosomes were genotyped for 2281 daughters of the 12 sires. Whole chromosome scanning was done using single marker analysis and interval mapping using three different algorithms. The analysis was carried out sire family wise. QTLs (63) were identified in single marker analysis and 32 QTLs were identified using interval mapping. The significance of LOD score was tested using permutation tests. The metaQTL analysis was carried out to find out the consensus chromosomal regions associated with milk yield in buffaloes. Five models were utilised and the best was selected on the basis of Akaike Information content. Total 23 chromosomal regions were identified for milk yield in buffaloes. 2 metaQTL chromosomal regions were identified on buffalo chromosome BBU2q; 3 metaQTLs each on buffalo chromosomes BBU8, BBU10 and BBU15 and 4 metaQTL regions each on BBU1q, BBU6, BBU9.

Genetic association of the ACACB gene with milk yield and composition traits in dairy cattle

Previously, we re-sequenced the whole genomes of eight Holstein bulls with high or low milk protein and fat percentage, and we detected two indels in the ACACB (acetyl-CoA carboxylase beta) gene that were polymorphic between the two groups. Thus, we considered ACACB as a promising candidate gene potentially affecting milk composition traits. Herein, we verified the genetic effects of ACACB on five milk traits in a Chinese Holstein population. We identified six SNPs in the 5'-promoter region, five in the 5'- untranslated region (UTR), 11 in exons, four in the 3'-UTR and three in the 3'-flanking region by re-sequencing the entire coding and regulatory regions of ACACB. One of these SNPs (ss1987461005) is reported here for the first time, and three of the SNPs (rs109482081, rs110819816 and rs109281947) were predicted to result in amino acid replacements. Genotype-phenotype association analyses showed that all the identified SNPs, except for ss1987461005, rs208919019 and rs134447911, were significantly associated with milk yield, fat yield, fat percentage, protein yield or protein percentage (P < 0.0001 to 0.0484). Linkage disequilibrium analyses were conducted among the identified SNPs to confirm the genetic associations. Two SNPs-rs135874354 (g.66218726T>C) and rs210928430 (g.66218117G>A)-were predicted to alter transcription factor binding sites in the 5'-promoter region of ACACB. A luciferase activity assay showed that the promoter activity of haplotype TG was significantly higher than that of CG (P = 0.0002) and that the promoter activity of haplotype TA was remarkably higher than that of CA (P = 7.4285E-09), showing that the T allele of rs135874354 increased promoter activity. Thus, rs135874354 was considered to be a potentially functional mutation. Our findings have, for the first time, profiled the genetic effect of ACACB on milk production traits in dairy cattle and revealed a potentially causal mutation that requires further the in-depth validation.

Genome-wide association analyses based on a multiple-trait approach for modeling feed efficiency

Genome-wide association (GWA) of feed efficiency (FE) could help target important genomic regions influencing FE. Data provided by an international dairy FE research consortium consisted of phenotypic records on dry matter intakes (DMI), milk energy (MILKE), and metabolic body weight (MBW) on 6,937 cows from 16 stations in 4 counties. Of these cows, 4,916 had genotypes on 57,347 single nucleotide polymorphism (SNP) markers. We compared a GWA analysis based on the more classical residual feed intake (RFI) model with one based on a previously proposed multiple trait (MT) approach for modeling FE using an alternative measure (DMI|MILKE,MBW). Both models were based on a single-step genomic BLUP procedure that allowed the use of phenotypes from both genotyped and nongenotyped cows. Estimated effects for single SNP markers were small and not statistically important but virtually identical for either FE measure (RFI vs. DMI|MILKE,MBW). However, upon further refining this analysis to develop joint tests within nonoverlapping 1-Mb windows, significant associations were detected between either measure of FE with a window on each of Bos taurus autosomes BTA12 and BTA26. There was, as expected, no overlap between detected genomic regions for DMI|MILKE,MBW and genomic regions influencing the energy sink traits (i.e., MILKE and MBW) because of orthogonal relationships clearly defined between the various traits. Conversely, GWA inferences on DMI can be demonstrated to be partly driven by genetic associations between DMI with these same energy sink traits, thereby having clear implications when comparing GWA studies on DMI to GWA studies on FE-like measures such as RFI.

Annotation of the goat genome using next generation sequencing of microRNA expressed by the lactating mammary gland: comparison of three approaches

Background

MicroRNAs (miRNA) are small endogenous non-coding RNA involved in the post-transcriptional regulation of specific mRNA targets. The first whole goat genome sequence became available in 2013, with few annotations. Our goal was to establish a list of the miRNA expressed in the mammary gland of lactating goats, thus enabling implementation of the goat miRNA repertoire and considerably enriching annotation of the goat genome.

Results

Here, we performed high throughput RNA sequencing on 10 lactating goat mammary glands. The bioinformatic detection of miRNA was carried out using miRDeep2 software. Three different methods were used to predict, quantify and annotate the sequenced reads. The first was a de novo approach based on the prediction of miRNA from the goat genome only. The second approach used bovine miRNA as an external reference whereas the last one used recently available goat miRNA. The three methods enabled the prediction and annotation of hundreds of miRNA, more than 95% were commonly identified. Using bovine miRNA, 1,178 distinct miRNA were detected, together with the annotation of 88 miRNA for which corresponding precursors could not be retrieved in the goat genome, and which were not detected using the de novo approach or with the use of goat miRNA. Each chromosomal coordinate of the precursors determined here were generated and depicted on a reference localisation map. Forty six goat miRNA clusters were also reported. The study revealed 263 precursors located in goat protein-coding genes, amongst which the location of 43 precursors was conserved between human, mouse and bovine, revealing potential new gene regulations in the goat mammary gland. Using the publicly available cattle QTL database, and cow precursors conserved in the goat and expressed in lactating mammary gland, 114 precursors were located within known QTL regions for milk production and composition.

Conclusions

The results reported here represent the first major identification study on miRNA expressed in the goat mammary gland at peak lactation. The elements generated by this study will now be used as references to decipher the regulation of miRNA expression in the goat mammary gland and to clarify their involvement in the lactation process.

Electronic supplementary material

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

Genetics of bovine abomasal displacement

Displacement of the abomasum (DA) is a common inherited condition in Holstein cows. This article reviews the genetics of DA including risk factors, genetic parameters and molecular genetic results. Breeds other than Holsteins affected by DA include Guernseys, Jerseys, Brown Swiss, Ayrshires and Simmental-Red Holsteins. In most DA cases, left displacements of the abomasum (LDA) are seen. Lactation incidence rates are higher for DA in first lactation Holsteins compared to later lactations. For Holstein cows, heritability estimates for DA are between 0.03 and 0.53. Genetic correlation estimates among DA and milk production traits range from positive to negative. Genome-wide significant genomic regions associated with LDA are located on bovine chromosomes (BTA) 1, 3, 11, 20 and 23. Motilin-associated single nucleotide polymorphisms on BTA23 exhibit a functional relationship with LDA. Pathways for deposition of calcium, insulin-dependent diabetes mellitus and synaptic transmission are significantly related to LDA in Holsteins. Deciphering the DA-associated genomic regions and genes may be an important step in the quest to understand the underlying disease-causing mechanisms and in unravelling mutations with a causal relationship to DA.

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.

In silico QTL mapping of maternal nurturing ability with the mouse diversity panel

Significant variation exists for maternal nurturing ability in inbred mice. Although classical mapping approaches have identified quantitative trait loci (QTL) that may account for this variation, the underlying genes are unknown. In this study, lactation performance data among the mouse diversity panel were used to map genomic regions associated with this variation. Females from each of 32 inbred strains (n = 8-19 dams/strain) were studied during the first 8 days of lactation by allowing them to raise weight- and size-normalized cross-foster litters (10 pups/litter). Average daily weight gain (ADG) of litters served as the primary indicator of milk production. The number of pups successfully reared to 8 days (PNUM8) also served as a related indicator of maternal performance. Initial haplotype association analysis using a Bonferroni-corrected, genome-wide threshold revealed 10 and 15 associations encompassing 11 and 13 genes for ADG and PNUM8, respectively. The most significant of these associated haplotype blocks were found on MMU 8, 11, and 19 and contained the genes Nr3c2, Egfr, Sec61g, and Gnaq. Lastly, two haplotype blocks on MMU9 were detected in association with PNUM8. These overlapped with the previously described maternal performance QTL, Neogq1. These results suggest that the application of in silico QTL mapping is a useful tool in discovering the presence of novel candidate genes involved in determining lactation capacity in mice.

Genetic effects of stearoyl-coenzyme A desaturase (SCD) polymorphism on milk production traits in the Chinese dairy population

Stearoyl-CoA desaturase (SCD) is a multifunctional complex enzyme important in the cellular biosynthesis of fatty acids. The present study was to investigate the association of the SCD gene with milk production traits in dairy cattle. Two single nucleotide polymorphisms (SNPs) (g.6926A>G and g.8646A>G) in introns 3 and 4, and three SNPs (g.10153A>G, g.10213T>C and g.10329C>T) in exon 5 were identified with pooled DNA sequencing and genotyped using matrix-assisted laser desorption/ionization time of flight mass spectrometry assay in 752 Chinese Holstein cows. Polymorphism g.10329C>T was predicted to result in an amino acid replacement from alanine to valine in the SCD protein. With a mixed animal model, the significant associations of the five SNPs with 305-day milk, fat and protein yields and protein percentage were determined. We further demonstrated cows with heterozygous genotypes (A/G or C/T) had highest 305 day milk yield, fat yield, protein yield and lowest protein percentage. Heterozygous cows with genotype AG at the g.6926A>G locus showed the greatest milk yield (P < 0.0001), fat yield (P < 0.0001) and protein yield (P < 0.0001) among other heterozygous genotypes at any of the loci. Dominance effects of all identified SNPs on milk, fat and protein yields and protein percentage were significant. Moreover, significant allele substitution effects at g.6926A>G locus on milk yield and at g.10213T>C on protein yield were observed. Five-locus haplotypes and strong linkage disequilibrium (D' > 0.9) between the five SNPs were also observed. The results suggest that identified polymorphisms could be potential genetic markers to improve the production performance of Chinese Holstein.

Genome-wide associations for milk production and somatic cell score in Holstein-Friesian cattle in Ireland

BACKGROUND

Contemporary dairy breeding goals have broadened to include, along with milk production traits, a number of non-production-related traits in an effort to improve the overall functionality of the dairy cow. Increased indirect selection for resistance to mastitis, one of the most important production-related diseases in the dairy sector, via selection for reduced somatic cell count has been part of these broadened goals. A number of genome-wide association studies have identified genetic variants associated with milk production traits and mastitis resistance, however the majority of these studies have been based on animals which were predominantly kept in confinement and fed a concentrate-based diet (i.e. high-input production systems). This genome-wide association study aims to detect associations using genotypic and phenotypic data from Irish Holstein-Friesian cattle fed predominantly grazed grass in a pasture-based production system (low-input).

RESULTS

Significant associations were detected for milk yield, fat yield, protein yield, fat percentage, protein percentage and somatic cell score using separate single-locus, frequentist and multi-locus, Bayesian approaches. These associations were detected using two separate populations of Holstein-Friesian sires and cows. In total, 1,529 and 37 associations were detected in the sires using a single SNP regression and a Bayesian method, respectively. There were 103 associations in common between the sires and cows across all the traits. As well as detecting associations within known QTL regions, a number of novel associations were detected; the most notable of these was a region of chromosome 13 associated with milk yield in the population of Holstein-Friesian sires.

CONCLUSIONS

A total of 276 of novel SNPs were detected in the sires using a single SNP regression approach. Although obvious candidate genes may not be initially forthcoming, this study provides a preliminary framework upon which to identify the causal mechanisms underlying the various milk production traits and somatic cell score. Consequently this will deepen our understanding of how these traits are expressed.

Multisite haplotype on cattle chromosome 3 is associated with quantitative trait locus effects on lactation traits

The goal of this study was to identify candidate genes and DNA polymorphisms for quantitative trait loci (QTL) affecting milk yield (MY), fat yield (FY), and protein yield (PY) previously mapped to bovine chromosome 3 (BTA3). To accomplish this, 373 half-siblings sired by three bulls previously shown to be segregating for lactation trait QTL, and 263 additional sires in the U.S. Dairy Bull DNA Repository (DBDR) were genotyped for 2,500 SNPs within a 16.3 Mbp QTL critical region on BTA3. Targeted resequencing of ∼1.8 Mbp within the QTL critical region of one of the QTL heterozygous sires identified additional polymorphisms useful for association studies. Twenty-three single nucleotide polymorphisms (SNPs) within a fine-mapped region were associated with effects on breeding values for MY, FY, or PY in DBDR sires, of which five SNPs were in strong linkage disequilibrium in the population. This multisite haplotype included SNPs located within exons or promoters of four tightly linked genes: RAP1A, ADORA3, OVGP1, and C3H1orf88. An SNP within RAP1A showed strong evidence of a recent selective sweep based on integrated haplotype score and was also associated with breeding value for PY. Because of its known function in alveolar lumen formation in the mammary gland, RAP1A is thus a strong candidate gene for QTL effects on lactation traits. Our results provide a detailed assessment of a QTL region that will be a useful guide for complex traits analysis in humans and other noninbred species.

Genetic diversity, structure, and breed relationships in Iberian cattle

In Iberia there are 51 officially recognized cattle breeds of which 15 are found in Portugal and 38 in Spain. We present here a comprehensive analysis of the genetic diversity and structure of Iberian cattle. Forty of these breeds were genotyped with 19 highly polymorphic microsatellite markers. Asturiana de los Valles displayed the greatest allelic diversity and Mallorquina the least. Unbiased heterozygosity values ranged from 0.596 to 0.787. The network based on Reynolds distances was star-shaped with few pairs of interrelated breeds and a clear cluster of 4 breeds (Alistana/Arouquesa/Marinhoa/Mirandesa). The analysis of the genetic structure of Iberian cattle indicated that the most probable number of population clusters included in the study would be 36. Distance results were supported by the STRUCTURE software indicating a relatively recent origin or possible crossbreeding or both between pairs or small groups of breeds. Five clusters included 2 different breeds (Betizu/Pirenaica, Morucha/Avileña, Parda de Montaña/Bruna de los Pirineos, Barrosã/Cachena, and Toro de Lidia/Brava de Lide), 3 breeds (Berrenda en Negro, Negra Andaluza, and Mertolenga) were divided in 2 independent clusters each, and 2 breeds were considered admixed (Asturiana de los Valles and Berrenda en Colorado). Individual assignation to breeds was not possible in the 2 admixed breeds and the pair Parda de Montaña/Bruna de los Pirineos. The relationship between Iberian cattle reflects their geographical origin rather than their morphotypes. Exceptions to this geographic clustering are most probably a consequence of crossbreeding with foreign breeds. The relative genetic isolation within their geographical origin, the consequent genetic drift, the adaptation to specific environment and production systems, and the influence of African and European cattle have contributed to the current genetic status of Iberian cattle, which are grouped according to their geographical origin. The greater degree of admixture observed in some breeds should be taken into account before using molecular markers for genetic assignment of individuals to breeds.

QTL mapping in outbred half-sib families using Bayesian model selection

In this article, we propose a model selection method, the Bayesian composite model space approach, to map quantitative trait loci (QTL) in a half-sib population for continuous and binary traits. In our method, the identity-by-descent-based variance component model is used. To demonstrate the performance of this model, the method was applied to map QTL underlying production traits on BTA6 in a Chinese half-sib dairy cattle population. A total of four QTLs were detected, whereas only one QTL was identified using the traditional least square (LS) method. We also conducted two simulation experiments to validate the efficiency of our method. The results suggest that the proposed method based on a multiple-QTL model is efficient in mapping multiple QTL for an outbred half-sib population and is more powerful than the LS method based on a single-QTL model.

A genome-wide scan for signatures of recent selection in Holstein cattle

The data from the newly available 50 K SNP chip was used for tagging the genome-wide footprints of positive selection in Holstein-Friesian cattle. For this purpose, we employed the recently described Extended Haplotype Homozygosity test, which detects selection by measuring the characteristics of haplotypes within a single population. To assess formally the significance of these results, we compared the combination of frequency and the Relative Extended Haplotype Homozygosity value of each core haplotype with equally frequent haplotypes across the genome. A subset of the putative regions showing the highest significance in the genome-wide EHH tests was mapped. We annotated genes to identify possible influence they have in beneficial traits by using the Gene Ontology database. A panel of genes, including FABP3, CLPN3, SPERT, HTR2A5, ABCE1, BMP4 and PTGER2, was detected, which overlapped with the most extreme P-values. This panel comprises some interesting candidate genes and QTL, representing a broad range of economically important traits such as milk yield and composition, as well as reproductive and behavioural traits. We also report high values of linkage disequilibrium and a slower decay of haplotype homozygosity for some candidate regions harbouring major genes related to dairy quality. The results of this study provide a genome-wide map of selection footprints in the Holstein genome, and can be used to better understand the mechanisms of selection in dairy cattle breeding.

The bovine annexin 9 gene (ANXA9) is significantly associated with milk-fat yield in a Spanish Holstein-Friesian population

On the basis of QTL studies for milk-fat yield trait on BTA3, annexin 9 protein (ANXA9), fatty acid transport protein type 3 (SLC27A3) and diacylglycerol O-acyltransferase 1 (DGAT1) were selected as candidate genes. Three different single nucleotide polymorphisms (SNPs) of bovine ANXA9, SLC27A3 and DGAT1 genes have been tested in a selective genotyping design for milk-fat yield. Significant allele frequency differences were found for ANXA9 (p=0.02), in Holstein-Friesian animals with high and low breeding values for milk-fat yield. Regression analysis also showed a significant effect (p=0.0207) between estimated breeding values (EBVs) for fat milk content and ANXA9 polymorphism. So ANXA9 gene falls into a significant quantitative trait loci interval for milk-fat yield that was previously reported on bovine chromosome 3 in other dairy populations. Our results suggest that the ANXA9 gene polymorphism or a linked segregating QTL contributes to variation in milk-fat yield.

Genome-wide association study of growth in crossbred beef cattle

Chromosomal regions harboring variation affecting cattle birth weight and BW gain to 1 yr of age were identified by marker association using the highly parallel BovineSNP50 BeadChip (50K) assay composed of 54,001 individual SNP. Genotypes were obtained from progeny (F(1); 590 steers) and 2-, 3-, and 4-breed cross grandprogeny (F(1)(2) = F(1) x F(1); 1,306 steers and 707 females) of 150 AI sires representing 7 breeds (22 sires per breed; Angus, Charolais, Gelbvieh, Hereford, Limousin, Red Angus, and Simmental). Genotypes and birth, weaning, and yearling BW records were used in whole-genome association analyses to estimate effects of individual SNP on growth. Traits analyzed included growth component traits: birth weight (BWT), 205-d adjusted birth to weaning BW gain (WG), 160-d adjusted postweaning BW gain (PWG); cumulative traits: 205-d adjusted weaning weight (WW = BWT + WG) and 365-d adjusted yearling weight (YW = BWT + WG + PWG); and indexes of relative differences between postnatal growth and birth weight. Modeled fixed effects included additive effects of calf and dam SNP genotype, year-sex-management contemporary groups, and covariates for calf and dam breed composition and heterosis. Direct and maternal additive polygenic effects and maternal permanent environment effects were random. Missing genotypes, including 50K genotypes of most dams, were approximated with a single-locus BLUP procedure from pedigree relationships and known 50K genotypes. Various association criteria were applied: stringent tests to account for multiple testing but with limited power to detect associations with small effects, and relaxed nominal P that may detect SNP associated with small effects but include excessive false positive associations. Genomic locations of the 231 SNP meeting stringent criteria generally coincided with described previously QTL affecting growth traits. The 12,425 SNP satisfying relaxed tests were located throughout the genome. Most SNP associated with BWT and postnatal growth affected components in the same direction, although detection of SNP associated with one component independent of others presents a possible opportunity for SNP-assisted selection to increase postnatal growth relative to BWT.

Bayesian shrinkage mapping for multiple QTL in half-sib families

Recently, an effective Bayesian shrinkage estimation method has been proposed for mapping QTL in inbred line crosses. However, with regard to outbred populations, such as half-sib populations with maternal information unavailable, it is not straightforward to utilize such a shrinkage estimation for QTL mapping. The reasons are: (1) the linkage phase of markers in the outbred population is usually unknown; and (2) only paternal genotypes can be used for inferring QTL genotypes of offspring. In this article, a novel Bayesian shrinkage method was proposed for mapping QTL under the half-sib design using a mixed model. A simulation study clearly demonstrated that the proposed method was powerful for detecting multiple QTL. In addition, we applied the proposed method to map QTL for economic traits in the Chinese dairy cattle population. Two or more novel QTL harbored in the chromosomal region were detected for each trait of interest, whereas only one QTL was found using traditional maximum likelihood analyses in our earlier studies. This further validated that our shrinkage estimation method could perform well in empirical data analyses and had practical significance in the field of linkage studies for outbred populations.

QTL detection for milk production traits in goats using a longitudinal model

Eight paternal half-sib families were used to identify chromosomal regions associated with variation in the lactation curves of dairy goats. DNA samples from 162 animals were amplified by PCR for 37 microsatellite markers, from Capra hircus autosomes CHI3, CHI6, CHI14 and CHI20. Milk samples were collected during 6 years, and there were 897 records for milk yield (MY) and 814 for fat (FP) and protein percentage (PP). The analysis was conducted in two stages. First, a random regression model with several fixed effects was fitted to describe the lactation function, using a scale (alpha) plus four shape parameters: beta and gamma, both associated with a decrease in the slope of the curve, and delta and phi that are related to the increase in slope. Predictions of alpha, beta, gamma, delta and phi were regressed using an interval mapping model, and F-tests were used to test for quantitative trait loci (QTL) effects. Significant (p < 0.05) QTLs were found for: (i) MY: CHI6 at 70-80 cM for all parameters; CHI14 at 14 cM for delta and phi; (ii) FP: CHI14, at 63 cM was associated with beta; CHI20, at 72 cM, showed association with alpha; (iii) PP: chromosomal regions associated with beta were found at 59 cM in CHI3 and at 55 cM in CHI20 with alpha and gamma. Analyses using more families and more animals will be useful to confirm or to reject these findings.

Isolation, mapping and identification of SNPs for four genes (ACP6, CGN, ANXA9, SLC27A3) from a bovine QTL region on BTA3

On the basis of fine mapping of a quantitative trait loci region of BTA3 for milk fat content, an examination of the comparative map between cattle and human indicates that the annexin 9 protein gene (ANXA9) and the fatty acid transport protein type 3 gene (SLC27A3) are two strong candidate genes. The objective of the present study is to isolate, map and characterize these genes and identify polymorphisms that could be further utilized in linkage or association studies. Furthermore, two new genes which are in the same region, cingulin protein gene (CGN) and lysophosphatidic acid phosphatase protein gene (ACP6) were studied. DNA fragments (869, 1778, 1933 and 2618 bp) corresponding to partial sequences of ACP6,CGN,ANXA9 and SLC27A3 genes were isolated. Direct sequencing of PCR products amplified from different cattle breeds revealed 1, 4, 4 and 2 SNPs for ACP6, CGN,ANXA9 and SLC27A3, respectively. For ANXA9 one SNP was located in exon 5 (A-->G 951) resulting in an amino acid change from histidine to arginine. Finally, ACP6,CGN,ANXA9 and SLC27A3 genes were located on chromosome 3 between ILSTS096 and BMS819 markers, in a region in which quantitative trait loci (QTL) for several milk traits have been described.

Fine mapping of genes on sheep chromosome 1 and their association with milk traits

On the basis of comparative mapping between cattle/sheep and human for milk trait quantitative trait loci (QTL) on BTA3/OAR1, annexin A9 (ANXA9) and solute carrier family 27 (fatty acid transporter), member 3 (SLC27A3) were selected as candidate genes for fat content (FC) in sheep milk. Two other genes in the same region, cingulin (CGN) and acid phosphatase 6, lysophosphatidic (ACP6), were also considered. DNA fragments of 1931 and 2790 bp corresponding to ANXA9 and SLC27A3 respectively were isolated, and 14 and 6 single nucleotide polymorphisms (SNPs) respectively were found in each gene. ANXA9, SLC27A3, CGN and ACP6 were localized to chromosome 1 between INRA006 and AE57 by linkage mapping using the International Mapping Flock. Across-family analyses of a daughter design comprising 13 sire families revealed significant sire and SLC27A3 genotype-nested-within-sire effects for FC. Within-family analyses indicated significant regression coefficients for FC in four of six heterozygous sires. These results could reflect the existence of a QTL for FC linked to SLC27A3 in sheep.

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.

An independent confirmation of a quantitative trait locus for milk yield and composition traits on bovine chromosome 26

Several reports have demonstrated that bovine chromosome 26 (BTA26) harbours significant or suggestive quantitative trait loci (QTL) for milk production and composition traits in dairy cattle. Our previous study showed that a C/T substitution in the bovine TCF7L2 gene on BTA26 was significantly linked to QTL for protein yield (PY) in a Canadian dairy cattle population. Actually, this polymorphism was one of the markers derived from a genome-wide screening of QTL for milk PY using an amplified fragment length polymorphism technique combined with a DNA pooling strategy. In the present study, 990 Holstein bulls with complete genotype and phenotype data from 14 sire families were analysed to confirm, if the QTL effects exist in other populations. Statistical analysis revealed that this marker was significantly associated with PY, protein percentage, milk yield and fat yield (FY) (p < 0.001) in the US Holstein population. These results indicate that this QTL region has a pleiotrophic effecton different milk traits and is portable in different populations.

Quantitative Trait Loci Affecting Health Traits in Swedish Dairy Cattle

The purpose of this study was to map quantitative trait loci (QTL) affecting health traits in Swedish dairy cattle. A genome scan covering 17 chromosomes was performed. Ten grandsire families were used in a granddaughter design. Nine of the families belonged to the Swedish Red and White breed, which is related to other Nordic Ayrshire breeds, and one family was of the Swedish Holstein breed. A total of 417 bulls were genotyped for 116 microsatellite markers distributed over 17 chromosomes. Daughter yield deviations for clinical mastitis, somatic cell count (SCC), and other diseases (OD) were included in the analysis. Least squares interval mapping using putative QTL as cofactors was applied both within and across grandsire families. Significance thresholds were set by permutation tests. In the across-family analysis, we detected 8 suggestive QTL and 3 QTL significant at the genome level. The QTL affecting clinical mastitis were found on 3 chromosomes (9, 11, and 25), 4 QTL for SCC were found (on chromosomes 5, 9, 11, and 23), and we detected 4 QTL for OD (on chromosomes 9, 11, 15, and 25). In addition, we found several QTL that segregated within single families but where the QTL effect was not significant in the across-family analysis. In conclusion, we were able to locate QTL for all 3 analyzed traits, and overlapping QTL for several traits were observed.

Quantitative trait loci mapping for dairy cattle production traits using a maximum likelihood method

A maximum likelihood method was developed for QTL mapping in half-sib designs and compared to the regression method in analyses of both field and simulated data. The field data consisted of milk production evaluations of 433 progeny tested sons of 6 sires and 64 microsatellite markers distributed over 12 chromosomes. Based on permutation tests, 5 significant QTL were detected in the field data by the regression method compared with 10 by the maximum likelihood method (P < 0.05). In field data analysis, the maximum likelihood method detected more significant QTL and had a smaller residual variance than the regression method. The simulation included 9 scenarios differing in number of families, family size, QTL variance, and marker density, each replicated 100 times. The simulation results suggested that, as for the regression method, the precision of estimating QTL from the maximum likelihood method improves with increasing number of sons per sire, increasing the ratio of QTL to phenotypic variance, and decreasing marker interval. The maximum likelihood method had a smaller dispersion of estimated QTL positions than the regression method in 6 of 9 scenarios simulated. Overall, the maximum likelihood method shows potential advantage in QTL detection over the regression method, especially in the situations with less favorable conditions for QTL detection.

A Complete Genome Scan of the Israeli Holstein Population for Quantitative Trait Loci by a Daughter Design

Eleven Israeli Holstein families including 5221 cows were analyzed by a daughter design for eight economic traits: milk, fat and protein production, fat and protein percentage, somatic cell score (SCS), herd-life, and female fertility. The cows were genotyped for 73 microsatellites with maximum spacing between markers of 53 cM. There were 86,304 informative genotypes. Preliminary analysis was by ANOVA of each trait, with the marker effect nested within sire. Significance was determined by controlling the false discovery rate at 0.4, after excluding markers with genome-wide significance for at least a single trait, and traits without any significant effects at this level. Thus, four markers on chromosomes 6 and 14 and female fertility were excluded. There remained 40 significant marker-trait combinations, and it is expected that 24 of these are true effects. To perform interval mapping for the families with significant contrasts, 21 additional markers were genotyped on chromosomes 2, 7, and 27. The bootstrap confidence intervals for gene effect did not include zero for protein percent on chromosome 2 and fat yield, protein yield, and SCS on chromosome 7. Quantitative trait locus heterozygosity was 33%, which is consistent with the hypothesis that only two alleles are segregating with unequal allele frequency.