ABSTRACT
BACKGROUND: Bovine lactoferrin (Lf) is an iron absorbing whey protein with antibacterial, antiviral, and antifungal activity. Lactoferrin is economically valuable and has an extremely variable concentration in milk, partly driven by environmental influences such as milking frequency, involution, or mastitis. A significant genetic influence has also been previously observed to regulate lactoferrin content in milk. Here, we conducted genetic mapping of lactoferrin protein concentration in conjunction with RNA-seq, ChIP-seq, and ATAC-seq data to pinpoint candidate causative variants that regulate lactoferrin concentrations in milk. RESULTS: We identified a highly-significant lactoferrin protein quantitative trait locus (pQTL), as well as a cis lactotransferrin (LTF) expression QTL (cis-eQTL) mapping to the LTF locus. Using ChIP-seq and ATAC-seq datasets representing lactating mammary tissue samples, we also report a number of regions where the openness of chromatin is under genetic influence. Several of these also show highly significant QTL with genetic signatures similar to those highlighted through pQTL and eQTL analysis. By performing correlation analysis between these QTL, we revealed an ATAC-seq peak in the putative promotor region of LTF, that highlights a set of 115 high-frequency variants that are potentially responsible for these effects. One of the 115 variants (rs110000337), which maps within the ATAC-seq peak, was predicted to alter binding sites of transcription factors known to be involved in lactation-related pathways. CONCLUSIONS: Here, we report a regulatory haplotype of 115 variants with conspicuously large impacts on milk lactoferrin concentration. These findings could enable the selection of animals for high-producing specialist herds.
Subject(s)
Lactation , Lactoferrin , Milk , Animals , Female , Haplotypes , Lactation/genetics , Lactoferrin/genetics , Lactoferrin/analysis , Lactoferrin/metabolism , Milk/chemistry , Milk/metabolism , CattleABSTRACT
BACKGROUND: Deleterious recessive conditions have been primarily studied in the context of Mendelian diseases. Recently, several deleterious recessive mutations with large effects were discovered via non-additive genome-wide association studies (GWAS) of quantitative growth and developmental traits in cattle, which showed that quantitative traits can be used as proxies of genetic disorders when such traits are indicative of whole-animal health status. We reasoned that lactation traits in cattle might also reflect genetic disorders, given the increased energy demands of lactation and the substantial stresses imposed on the animal. In this study, we screened more than 124,000 cows for recessive effects based on lactation traits. RESULTS: We discovered five novel quantitative trait loci (QTL) that are associated with large recessive impacts on three milk yield traits, with these loci presenting missense variants in the DOCK8, IL4R, KIAA0556, and SLC25A4 genes or premature stop variants in the ITGAL, LRCH4, and RBM34 genes, as candidate causal mutations. For two milk composition traits, we identified several previously reported additive QTL that display small dominance effects. By contrasting results from milk yield and milk composition phenotypes, we note differing genetic architectures. Compared to milk composition phenotypes, milk yield phenotypes had lower heritabilities and were associated with fewer additive QTL but had a higher non-additive genetic variance and were associated with a higher proportion of loci exhibiting dominance. CONCLUSIONS: We identified large-effect recessive QTL which are segregating at surprisingly high frequencies in cattle. We speculate that the differences in genetic architecture between milk yield and milk composition phenotypes derive from underlying dissimilarities in the cellular and molecular representation of these traits, with yield phenotypes acting as a better proxy of underlying biological disorders through presentation of a larger number of major recessive impacts.
Subject(s)
Genome-Wide Association Study , Quantitative Trait Loci , Animals , Cattle/genetics , Female , Lactation/genetics , Milk , PhenotypeABSTRACT
BACKGROUND: Milk samples from 10,641 dairy cattle were screened by a mass spectrometry method for extreme concentrations of the A or B isoforms of the whey protein, ß-lactoglobulin (BLG), to identify causative genetic variation driving changes in BLG concentration. RESULTS: A cohort of cows, from a single sire family, was identified that produced milk containing a low concentration of the BLG B protein isoform. A genome-wide association study (GWAS) of BLG B protein isoform concentration in milk from AB heterozygous cows, detected a group of highly significant single nucleotide polymorphisms (SNPs) within or close to the BLG gene. Among these was a synonymous G/A variation at position + 78 bp in exon 1 of the BLG gene (chr11:103256256G > A). The effect of the A allele of this SNP (which we named B') on BLG expression was evaluated in a luciferase reporter assay in transfected CHO-K1 and MCF-7 cells. In both cell types, the presence of the B' allele in a plasmid containing the bovine BLG gene from -922 to + 898 bp (relative to the transcription initiation site) resulted in a 60% relative reduction in mRNA expression, compared to the plasmid containing the wild-type B sequence allele. Examination of a mammary RNAseq dataset (n = 391) identified 14 heterozygous carriers of the B' allele which were homozygous for the BLG B protein isoform (BB'). The level of expression of the BLG B' allele was 41.9 ± 1.0% of that of the wild-type BLG B allele. Milk samples from three cows, homozygous for the A allele at chr11:103,256,256 (B'B'), were analysed (HPLC) and showed BLG concentrations of 1.04, 1.26 and 1.83 g/L relative to a mean of 4.84 g/L in milk from 16 herd contemporaries of mixed (A and B) BLG genotypes. The mechanism by which B' downregulates milk BLG concentration remains to be determined. CONCLUSIONS: High-throughput screening and identification of outliers, enabled the discovery of a synonymous G > A mutation in exon 1 of the B allele of the BLG gene (B'), which reduced the milk concentration of ß-lactoglobulin B protein isoform, by more than 50%. Milk from cows carrying the B' allele is expected to have improved processing characteristics, particularly for cheese-making.
Subject(s)
Lactoglobulins , Milk , Polymorphism, Single Nucleotide , Animals , Cattle/genetics , Female , Genome-Wide Association Study , Lactoglobulins/analysis , Milk/chemistry , Protein Isoforms/analysisABSTRACT
Fourier-transform mid-infrared (FT-MIR) spectroscopy is a high-throughput and inexpensive methodology used to evaluate concentrations of fat and protein in dairy cattle milk samples. The objective of this study was to compare the genetic characteristics of FT-MIR predicted fatty acids and individual milk proteins with those that had been measured directly using gas and liquid chromatography methods. The data used in this study was based on 2,005 milk samples collected from 706 Holstein-Friesian × Jersey animals that were managed in a seasonal, pasture-based dairy system, with milk samples collected across 2 consecutive seasons. Concentrations of fatty acids and protein fractions in milk samples were directly determined by gas chromatography and high-performance liquid chromatography, respectively. Models to predict each directly measured trait based on FT-MIR spectra were developed using partial least squares regression, with spectra from a random selection of half the cows used to train the models, and predictions for the remaining cows used as validation. Variance parameters for each trait and genetic correlations for each pair of measured/predicted traits were estimated from pedigree-based bivariate models using REML procedures. A genome-wide association study was undertaken using imputed whole-genome sequence, and quantitative trait loci (QTL) from directly measured traits were compared with QTL from the corresponding FT-MIR predicted traits. Cross-validation prediction accuracies based on partial least squares for individual and grouped fatty acids ranged from 0.18 to 0.65. Trait prediction accuracies in cross-validation for protein fractions were 0.53, 0.19, and 0.48 for α-casein, ß-casein, and κ-casein, 0.31 for α-lactalbumin, 0.68 for ß-lactoglobulin, and 0.36 for lactoferrin. Heritability estimates for directly measured traits ranged from 0.07 to 0.55 for fatty acids; and from 0.14 to 0.63 for individual milk proteins. For FT-MIR predicted traits, heritability estimates were mostly higher than for the corresponding measured traits, ranging from 0.14 to 0.46 for fatty acids, and from 0.30 to 0.70 for individual proteins. Genetic correlations between directly measured and FT-MIR predicted protein fractions were consistently above 0.75, with the exceptions of C18:0 and C18:3 cis-3, which had genetic correlations of 0.72 and 0.74, respectively. The GWAS identified trait QTL for fatty acids with likely candidates in the DGAT1, CCDC57, SCD, and GPAT4 genes. Notably, QTL for SCD were largely absent in the FT-MIR predicted traits, and QTL for GPAT4 were absent in directly measured traits. Similarly, for directly measured individual proteins, we identified QTL with likely candidates in the CSN1S1, CSN3, PAEP, and LTF genes, but the QTL for CSN3 and LTF were absent in the FT-MIR predicted traits. Our study indicates that genetic correlations between directly measured and FT-MIR predicted fatty acid and protein fractions are typically high, but that phenotypic variation in these traits may be underpinned by differing genetic architecture.
Subject(s)
Fatty Acids , Genome-Wide Association Study , Female , Cattle/genetics , Animals , Fatty Acids/metabolism , Genome-Wide Association Study/veterinary , Milk/chemistry , Milk Proteins/analysis , Caseins/analysisABSTRACT
BACKGROUND: Fourier-transform mid-infrared (FT-MIR) spectroscopy provides a high-throughput and inexpensive method for predicting milk composition and other novel traits from milk samples. While there have been many genome-wide association studies (GWAS) conducted on FT-MIR predicted traits, there have been few GWAS for individual FT-MIR wavenumbers. Using imputed whole-genome sequence for 38,085 mixed-breed New Zealand dairy cattle, we conducted GWAS on 895 individual FT-MIR wavenumber phenotypes, and assessed the value of these direct phenotypes for identifying candidate causal genes and variants, and improving our understanding of the physico-chemical properties of milk. RESULTS: Separate GWAS conducted for each of 895 individual FT-MIR wavenumber phenotypes, identified 450 1-Mbp genomic regions with significant FT-MIR wavenumber QTL, compared to 246 1-Mbp genomic regions with QTL identified for FT-MIR predicted milk composition traits. Use of mammary RNA-seq data and gene annotation information identified 38 co-localized and co-segregating expression QTL (eQTL), and 31 protein-sequence mutations for FT-MIR wavenumber phenotypes, the latter including a null mutation in the ABO gene that has a potential role in changing milk oligosaccharide profiles. For the candidate causative genes implicated in these analyses, we examined the strength of association between relevant loci and each wavenumber across the mid-infrared spectrum. This revealed shared association patterns for groups of genomically-distant loci, highlighting clusters of loci linked through their biological roles in lactation and their presumed impacts on the chemical composition of milk. CONCLUSIONS: This study demonstrates the utility of FT-MIR wavenumber phenotypes for improving our understanding of milk composition, presenting a larger number of QTL and putative causative genes and variants than found from FT-MIR predicted composition traits. Examining patterns of significance across the mid-infrared spectrum for loci of interest further highlighted commonalities of association, which likely reflects the physico-chemical properties of milk constituents.
Subject(s)
Cattle/genetics , Milk/chemistry , Quantitative Trait Loci , Animals , Genome-Wide Association Study , Hybridization, Genetic , Milk/standards , Oligosaccharides/analysis , Spectroscopy, Fourier Transform InfraredABSTRACT
BACKGROUND: The DGAT1 gene encodes an enzyme responsible for catalysing the terminal reaction in mammary triglyceride synthesis, and underpins a well-known pleiotropic quantitative trait locus (QTL) with a large influence on milk composition phenotypes. Since first described over 15 years ago, a protein-coding variant K232A has been assumed as the causative variant underlying these effects, following in-vitro studies that demonstrated differing levels of triglyceride synthesis between the two protein isoforms. RESULTS: We used a large RNAseq dataset to re-examine the underlying mechanisms of this large milk production QTL, and hereby report novel expression-based functions of the chr14 g.1802265AA > GC variant that encodes the DGAT1 K232A substitution. Using expression QTL (eQTL) mapping, we demonstrate a highly-significant mammary eQTL for DGAT1, where the K232A mutation appears as one of the top associated variants for this effect. By conducting in vitro expression and splicing experiments in bovine mammary cell culture, we further show modulation of splicing efficiency by this mutation, likely through disruption of an exon splice enhancer as a consequence of the allele encoding the 232A variant. CONCLUSIONS: The relative contributions of the enzymatic and transcription-based mechanisms now attributed to K232A remain unclear; however, these results suggest that transcriptional impacts contribute to the diversity of lactation effects observed at the DGAT1 locus.
Subject(s)
Diacylglycerol O-Acyltransferase , Lactation , Animals , Cattle , Diacylglycerol O-Acyltransferase/genetics , Exons , Female , Gene Expression , Milk , MutationABSTRACT
BACKGROUND: Over many years, artificial selection has substantially improved milk production by cows. However, the genes that underlie milk production quantitative trait loci (QTL) remain relatively poorly characterised. Here, we investigate a previously reported QTL located at the CSF2RB locus on chromosome 5, for several milk production phenotypes, to better understand its underlying genetic and molecular causes. RESULTS: Using a population of 29,350 taurine dairy cows, we conducted association analyses for milk yield and composition traits, and identified highly significant QTL for milk yield, milk fat concentration, and milk protein concentration. Strikingly, protein concentration and milk yield appear to show co-located yet genetically distinct QTL. To attempt to understand the molecular mechanisms that might be mediating these effects, gene expression data were used to investigate eQTL for 11 genes in the broader interval. This analysis highlighted genetic impacts on CSF2RB and NCF4 expression that share similar association signatures to those observed for lactation QTL, strongly implicating one or both of these genes as responsible for these effects. Using the same gene expression dataset representing 357 lactating cows, we also identified 38 novel RNA editing sites in the 3' UTR of CSF2RB transcripts. The extent to which two of these sites were edited also appears to be genetically co-regulated with lactation QTL, highlighting a further layer of regulatory complexity that involves the CSF2RB gene. CONCLUSIONS: This locus presents a diversity of molecular and lactation QTL, likely representing multiple overlapping effects that, at a minimum, highlight the CSF2RB gene as having a causal role in these processes.
Subject(s)
Cattle/genetics , Cytokine Receptor Common beta Subunit/genetics , Lactation/genetics , Phenotype , Quantitative Trait Loci , 3' Untranslated Regions , Animals , Cytokine Receptor Common beta Subunit/metabolism , Female , Male , Milk/metabolism , Phosphoproteins/geneticsABSTRACT
Inflammation of the mammary gland following bacterial infection, commonly known as mastitis, affects all mammalian species. Although the aetiology and epidemiology of mastitis in the dairy cow are well described, the genetic factors mediating resistance to mammary gland infection are not well known, due in part to the difficulty in obtaining robust phenotypic information from sufficiently large numbers of individuals. To address this problem, an experimental mammary gland infection experiment was undertaken, using a Friesian-Jersey cross breed F2 herd. A total of 604 animals received an intramammary infusion of Streptococcus uberis in one gland, and the clinical response over 13 milkings was used for linkage mapping and genome-wide association analysis. A quantitative trait locus (QTL) was detected on bovine chromosome 11 for clinical mastitis status using micro-satellite and Affymetrix 10 K SNP markers, and then exome and genome sequence data used from the six F1 sires of the experimental animals to examine this region in more detail. A total of 485 sequence variants were typed in the QTL interval, and association mapping using these and an additional 37 986 genome-wide markers from the Illumina SNP50 bovine SNP panel revealed association with markers encompassing the interleukin-1 gene cluster locus. This study highlights a region on bovine chromosome 11, consistent with earlier studies, as conferring resistance to experimentally induced mammary gland infection, and newly prioritises the IL1 gene cluster for further analysis in genetic resistance to mastitis.
Subject(s)
Mastitis, Bovine/genetics , Mastitis, Bovine/immunology , Streptococcal Infections/veterinary , Animals , Cattle , Chromosome Mapping/veterinary , Disease Models, Animal , Female , Genetic Predisposition to Disease , Genetic Variation , Genome-Wide Association Study/veterinary , Genotype , Hybridization, Genetic , Mastitis, Bovine/microbiology , Polymorphism, Single Nucleotide , Quantitative Trait Loci/geneticsABSTRACT
BACKGROUND: Lactose provides an easily-digested energy source for neonates, and is the primary carbohydrate in milk in most species. Bovine lactose is also a key component of many human food products. However, compared to analyses of other milk components, the genetic control of lactose has been little studied. Here we present the first GWAS focussed on analysis of milk lactose traits. RESULTS: Using a discovery population of 12,000 taurine dairy cattle, we detail 27 QTL for lactose concentration and yield, and subsequently validate the effects of 26 of these loci in a distinct population of 18,000 cows. We next present data implicating causative genes and variants for these QTL. Fine mapping of these regions using imputed, whole genome sequence-resolution genotypes reveals protein-coding candidate causative variants affecting the ABCG2, DGAT1, STAT5B, KCNH4, NPFFR2 and RNF214 genes. Eleven of the remaining QTL appear to be driven by regulatory effects, suggested by the presence of co-locating, co-segregating eQTL discovered using mammary RNA sequence data from a population of 357 lactating cows. Pathway analysis of genes representing all lactose-associated loci shows significant enrichment of genes located in the endoplasmic reticulum, with functions related to ion channel activity mediated through the LRRC8C, P2RX4, KCNJ2 and ANKH genes. A number of the validated QTL are also found to be associated with additional milk volume, fat and protein phenotypes. CONCLUSIONS: Overall, these findings highlight novel candidate genes and variants involved in milk lactose regulation, whose impacts on membrane transport mechanisms reinforce the key osmo-regulatory roles of lactose in milk.
Subject(s)
Lactose/metabolism , Membrane Transport Proteins/genetics , Milk/metabolism , Quantitative Trait Loci , Alleles , Animals , Cattle , Female , Gene Expression , Genetic Variation , Genome-Wide Association Study , Ion Transport/genetics , Lactation/genetics , Phenotype , Polymorphism, Single Nucleotide , Sequence Analysis, DNA , Sequence Analysis, RNAABSTRACT
The objective of this study was to estimate heritability and crossbreeding parameters (breed and heterosis effects) of various fatty acid (FA) concentrations in milk fat of New Zealand dairy cattle. For this purpose, calibration equations to predict concentration of each of the most common FAs were derived with partial least squares (PLS) using mid-infrared (MIR) spectral data from milk samples (n=850) collected in the 2003-04 season from 348 second-parity crossbred cows during peak, mid and late lactation. The milk samples produced both, MIR spectral data and concentration of the most common FAs determined using gas chromatography (GC). The concordance correlation coefficients (CCC) between the concentration of a FA determined by GC and the PLS equation ranged from 0.63 to 0.94, suggesting that some prediction equations can be considered to have substantial predictive ability. The PLS calibration equations were then used to predict the concentration of each of the fatty acids in 26,769 milk samples from 7385 cows that were herd-tested during the 2007-08 season. Data were analysed using a single-trait repeatability animal model. Shorter chain FA (16:0 and below) were significantly higher (P<0.05) in Jersey cows, while longer chain, including unsaturated longer chain FA were higher in Holstein-Friesian cows. The estimates of heritabilities ranged from 0.17 to 0.41 suggesting that selective breeding could be used to ensure milk fat composition stays aligned to consumer, market and manufacturing needs.
Subject(s)
Cattle/genetics , Fatty Acids/analysis , Milk/chemistry , Animals , Chromatography, Gas/veterinary , Female , Hybridization, Genetic/genetics , Quantitative Trait, Heritable , Spectrophotometry, Infrared/veterinaryABSTRACT
Mammalian species carry ~100 loss-of-function variants per individual1,2, where ~1-5 of these impact essential genes and cause embryonic lethality or severe disease when homozygous3. The functions of the remainder are more difficult to resolve, although the assumption is that these variants impact fitness in less manifest ways. Here we report one of the largest sequence-resolution screens of cattle to date, targeting discovery and validation of non-additive effects in 130,725 animals. We highlight six novel recessive loci with impacts generally exceeding the largest-effect variants identified from additive genome-wide association studies, presenting analogs of human diseases and hitherto-unrecognized disorders. These loci present compelling missense (PLCD4, MTRF1 and DPF2), premature stop (MUS81) and splice-disrupting (GALNT2 and FGD4) mutations, together explaining substantial proportions of inbreeding depression. These results demonstrate that the frequency distribution of deleterious alleles segregating in selected species can afford sufficient power to directly map novel disorders, presenting selection opportunities to minimize the incidence of genetic disease.
Subject(s)
Cattle Diseases/diagnosis , Cattle Diseases/etiology , Loss of Function Mutation , Phenotype , Alleles , Animals , Biomarkers , Cattle , Cattle Diseases/epidemiology , Genome-Wide Association Study , Genotype , Inbreeding , Incidence , SyndromeABSTRACT
A major pleiotropic quantitative trait locus (QTL) located at ~25 Mbp on bovine chromosome 14 affects a myriad of growth and developmental traits in Bos taurus and indicus breeds. These QTL have been attributed to two functional variants in the bidirectional promoter of PLAG1 and CHCHD7. Although PLAG1 is a good candidate for mediating these effects, its role remains uncertain given that these variants are also associated with expression of five additional genes at the broader locus. In the current study, we conducted expression QTL (eQTL) mapping of this region using a large, high depth mammary RNAseq dataset representing 375 lactating cows. Here we show that of the seven previously implicated genes, only PLAG1 and LYN are differentially expressed by QTL genotype, and only PLAG1 bears the same association signature of the growth and body weight QTLs. For the first time, we also report significant association of PLAG1 genotype with milk production traits, including milk fat, volume, and protein yield. Collectively, these data strongly suggest PLAG1 as the causative gene underlying this diverse range of traits, and demonstrate new effects for the locus on lactation phenotypes.
Subject(s)
Body Weight/genetics , Cattle/genetics , DNA-Binding Proteins/genetics , Genetic Pleiotropy , Milk/chemistry , Animals , Chromosomes, Mammalian/genetics , Dairying , Gene Expression Regulation , Genome-Wide Association Study , Polymorphism, Single Nucleotide/genetics , Quantitative Trait Loci/genetics , Quantitative Trait, Heritable , Reproducibility of ResultsABSTRACT
The mammary gland is a prolific lipogenic organ, synthesising copious amounts of triglycerides for secretion into milk. The fat content of milk varies widely both between and within species, and recent independent genome-wide association studies have highlighted a milk fat percentage quantitative trait locus (QTL) of large effect on bovine chromosome 5. Although both EPS8 and MGST1 have been proposed to underlie these signals, the causative status of these genes has not been functionally confirmed. To investigate this QTL in detail, we report genome sequence-based imputation and association mapping in a population of 64,244 taurine cattle. This analysis reveals a cluster of 17 non-coding variants spanning MGST1 that are highly associated with milk fat percentage, and a range of other milk composition traits. Further, we exploit a high-depth mammary RNA sequence dataset to conduct expression QTL (eQTL) mapping in 375 lactating cows, revealing a strong MGST1 eQTL underpinning these effects. These data demonstrate the utility of DNA and RNA sequence-based association mapping, and implicate MGST1, a gene with no obvious mechanistic relationship to milk composition regulation, as causally involved in these processes.
Subject(s)
Genome-Wide Association Study/methods , Glutathione Transferase/genetics , Milk/chemistry , Quantitative Trait Loci , Animals , Bayes Theorem , Cattle , Chromosome Mapping/veterinary , Chromosomes, Mammalian/genetics , Genome-Wide Association Study/veterinary , Polymorphism, Single Nucleotide , Sequence Analysis, DNA/veterinary , Sequence Analysis, RNA/veterinaryABSTRACT
A mechanobullous skin disorder was identified in the progeny of a 3-y-old Friesian-Jersey crossbred bull. The condition presented as loss of skin and mucosa from contact areas and inflammation. Examination of skin samples under light microscopy revealed separation of the epidermis from the dermis. Electron microscopic analysis refined the site of cleavage to above the basement membrane involving lysis of basal keratinocytes. These observations were consistent with the simplex form of epidermolysis bullosa (EB) in humans. Candidate genes based on human gene mutations were assessed, resulting in keratin 5 being identified as the most likely candidate gene. The sequence of bovine keratin 5 was established and sequencing led to identification of a G to A substitution in all affected animals. This mutation leads to an amino acid change of glutamic acid to lysine in the final E (478) of the KLLEGE motif of the protein. The sire carried a de novo mutation and was mosaic, explaining his asymptomatic status and the less than expected frequency of affected offspring. Remarkably, the same mutation has been previously described in EB simplex in humans.
Subject(s)
Cattle Diseases/genetics , Epidermolysis Bullosa Simplex/veterinary , Keratins/genetics , Mosaicism , Mutation , Amino Acid Motifs , Amino Acid Substitution , Animals , Animals, Newborn , Cattle , Cattle Diseases/pathology , Cattle Diseases/physiopathology , Chromosome Mapping , Female , Glutamic Acid , Keratin-5 , Lysine , Male , Microscopy, Electron , PedigreeABSTRACT
Selective breeding has strongly reduced the genetic diversity in livestock species, and contemporary breeding practices exclude potentially beneficial rare genetic variation from the future gene pool. Here we test whether important traits arising by new mutations can be identified and rescued in highly selected populations. We screened milks from 2.5 million cows to identify an exceptional individual which produced milk with reduced saturated fat content, and improved unsaturated and omega-3 fatty acid concentrations. The milk traits were transmitted dominantly to her offspring, and genetic mapping and genome sequencing revealed a new mutation in a previously unknown splice enhancer of the DGAT1 gene. Homozygous carriers show features of human diarrheal disorders, and may be useful for the development of therapeutic strategies. Our study demonstrates that high-throughput phenotypic screening can uncover rich genetic diversity even in inbred populations, and introduces a novel strategy to develop novel milks with improved nutritional properties.
Subject(s)
Diacylglycerol O-Acyltransferase/genetics , Milk/metabolism , Mutation, Missense , Animals , Base Sequence , Cattle/genetics , Fatty Acids/biosynthesis , Female , Genetic Association Studies , Lipid Metabolism/genetics , Male , Pedigree , Phenotype , Polymorphism, Single NucleotideABSTRACT
Milk is composed of a complex mixture of lipids, proteins, carbohydrates and various vitamins and minerals as a source of nutrition for young mammals. The composition of milk varies between individuals, with lipid composition in particular being highly heritable. Recent reports have highlighted a region of bovine chromosome 27 harbouring variants affecting milk fat percentage and fatty acid content. We aimed to further investigate this locus in two independent cattle populations, consisting of a Holstein-Friesian x Jersey crossbreed pedigree of 711 F2 cows, and a collection of 32,530 mixed ancestry Bos taurus cows. Bayesian genome-wide association mapping using markers imputed from the Illumina BovineHD chip revealed a large quantitative trait locus (QTL) for milk fat percentage on chromosome 27, present in both populations. We also investigated a range of other milk composition phenotypes, and report additional associations at this locus for fat yield, protein percentage and yield, lactose percentage and yield, milk volume, and the proportions of numerous milk fatty acids. We then used mammary RNA sequence data from 212 lactating cows to assess the transcript abundance of genes located in the milk fat percentage QTL interval. This analysis revealed a strong eQTL for AGPAT6, demonstrating that high milk fat percentage genotype is also additively associated with increased expression of the AGPAT6 gene. Finally, we used whole genome sequence data from six F1 sires to target a panel of novel AGPAT6 locus variants for genotyping in the F2 crossbreed population. Association analysis of 58 of these variants revealed highly significant association for polymorphisms mapping to the 5'UTR exons and intron 1 of AGPAT6. Taken together, these data suggest that variants affecting the expression of AGPAT6 are causally involved in differential milk fat synthesis, with pleiotropic consequences for a diverse range of other milk components.
Subject(s)
Cattle/genetics , Genetic Variation , Glycerol-3-Phosphate O-Acyltransferase/genetics , Lipogenesis/genetics , Milk/metabolism , Animals , Base Sequence , Chromosomes, Mammalian/genetics , Crosses, Genetic , Female , Genome-Wide Association Study , Genotype , Lactation/genetics , Lipids/genetics , Male , Mammary Glands, Animal/metabolism , Molecular Sequence Data , Phenotype , Quantitative Trait Loci/genetics , RNA, Messenger/genetics , RNA, Messenger/metabolismABSTRACT
Lactation, hair development and homeothermy are characteristic evolutionary features that define mammals from other vertebrate species. Here we describe the discovery of two autosomal dominant mutations with antagonistic, pleiotropic effects on all three of these biological processes, mediated through the prolactin signalling pathway. Most conspicuously, mutations in prolactin (PRL) and its receptor (PRLR) have an impact on thermoregulation and hair morphology phenotypes, giving prominence to this pathway outside of its classical roles in lactation.
Subject(s)
Body Temperature Regulation/genetics , Hair/metabolism , Lactation/genetics , Mutation , Prolactin/genetics , Signal Transduction , Amino Acid Sequence , Animals , Base Sequence , Body Temperature , Cattle , Female , Gene Expression , Genetic Pleiotropy , Genotype , Hair/ultrastructure , Homeostasis , Humans , Molecular Sequence Data , Phenotype , Prolactin/metabolism , Receptors, Prolactin/genetics , Receptors, Prolactin/metabolism , Skin/anatomy & histology , Skin/metabolismABSTRACT
We report mapping of a quantitative trait locus (QTL) with a major effect on bovine stature to a â¼780-kb interval using a Hidden Markov Model-based approach that simultaneously exploits linkage and linkage disequilibrium. We re-sequenced the interval in six sires with known QTL genotype and identified 13 clustered candidate quantitative trait nucleotides (QTNs) out of >9,572 discovered variants. We eliminated five candidate QTNs by studying the phenotypic effect of a recombinant haplotype identified in a breed diversity panel. We show that the QTL influences fetal expression of seven of the nine genes mapping to the â¼780-kb interval. We further show that two of the eight candidate QTNs, mapping to the PLAG1-CHCHD7 intergenic region, influence bidirectional promoter strength and affect binding of nuclear factors. By performing expression QTL analyses, we identified a splice site variant in CHCHD7 and exploited this naturally occurring null allele to exclude CHCHD7 as single causative gene.
Subject(s)
Cattle/anatomy & histology , Cattle/genetics , Genetic Variation , Alleles , Animals , Base Sequence , Chromosome Mapping , Crosses, Genetic , DNA Primers/genetics , Female , Genome-Wide Association Study , Haplotypes , Linkage Disequilibrium , Male , Open Reading Frames , Polymorphism, Single Nucleotide , Promoter Regions, Genetic , Quantitative Trait Loci , RNA Splice Sites , Regulon , Sequence Homology, Nucleic AcidABSTRACT
The imprints of domestication and breed development on the genomes of livestock likely differ from those of companion animals. A deep draft sequence assembly of shotgun reads from a single Hereford female and comparative sequences sampled from six additional breeds were used to develop probes to interrogate 37,470 single-nucleotide polymorphisms (SNPs) in 497 cattle from 19 geographically and biologically diverse breeds. These data show that cattle have undergone a rapid recent decrease in effective population size from a very large ancestral population, possibly due to bottlenecks associated with domestication, selection, and breed formation. Domestication and artificial selection appear to have left detectable signatures of selection within the cattle genome, yet the current levels of diversity within breeds are at least as great as exists within humans.