Michael Keehan

Variants modulating the expression of a chromosome domain encompassing PLAG1 influence bovine stature

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.

A high density linkage map of the bovine genome

BackgroundRecent technological advances have made it possible to efficiently genotype large numbers of single nucleotide polymorphisms (SNPs) in livestock species, allowing the production of high-density linkage maps. Such maps can be used for quality control of other SNPs and for fine mapping of quantitative trait loci (QTL) via linkage disequilibrium (LD).ResultsA high-density bovine linkage map was constructed using three types of markers. The genotypic information was obtained from 294 microsatellites, three milk protein haplotypes and 6769 SNPs. The map was constructed by combining genetic (linkage) and physical information in an iterative mapping process. Markers were mapped to 3,155 unique positions; the 6,924 autosomal markers were mapped to 3,078 unique positions and the 123 non-pseudoautosomal and 19 pseudoautosomal sex chromosome markers were mapped to 62 and 15 unique positions, respectively. The linkage map had a total length of 3,249 cM. For the autosomes the average genetic distance between adjacent markers was 0.449 cM, the genetic distance between unique map positions was 1.01 cM and the average genetic distance (cM) per Mb was 1.25.ConclusionThere is a high concordance between the order of the SNPs in our linkage map and their physical positions on the most recent bovine genome sequence assembly (Btau 4.0). The linkage maps provide support for fine mapping projects and LD studies in bovine populations. Additionally, the linkage map may help to resolve positions of unassigned portions of the bovine genome.

Sequence-based Association Analysis Reveals an MGST1 eQTL with Pleiotropic Effects on Bovine Milk Composition.

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.

Genome scan for the degree of white spotting in dairy cattle

White spotting is one of the most distinguishing visual characters in dairy cattle. There is considerable variation within and between breeds of cattle. The objective of this study was to map quantitative trait loci (QTL) affecting the degree of white spotting in dairy cattle based on an F(2) experimental design using Holstein-Friesian and Jersey crossbred cows. The genome scan was implemented using half-sib and line-of-descent approaches with high density markers. Significant QTL were found on chromosomes 6, 18 and 22. The mapped region on BTA6 confirmed the widely conserved KIT locus affecting mammalian pigmentation. Haplotype information linked the highly significant QTL on BTA22 to the Microphthalmia-associated transcription factor (MITF) gene, which has been reported to be associated with pigmentation traits in some other mammals.

Genome scan of pigmentation traits in Friesian-Jersey crossbred cattle.

Pigmentation traits expressed in animals are visual characteristics that allow us to distinguish between breeds and between strains within breed. The objective of this study was to map quantitative trait loci (QTLs) affecting the pigmentation traits in approximately 800 F(2) grand daughter dairy cattle from a Holstein-Friesian and Jersey cross breed cattle. Traits analyzed included pigmentation phenotypes on the body, teat and hoop. The phenoypes were collected from digital photos or visual inspection of live animals. QTL mapping was implemented using half-sib and line-of-descent inheritance models. Our analysis initially detected a number of significant QTLs on chromosomes: 2, 6, 13, 15, 18 and 22. The significant QTLs were divided into two groups: one group influencing the pigmentation color and the other group affecting the absence or level of pigmentation. The most significant QTL peaks were observed on Bovine taurus autosome 18 (BTA18) close to melanocortin 1 receptor (MC1R) for the color traits, on BTA6 close to the receptor tyrosine kinase (KIT) and BTA22 close to microphthalmia-associated transcription factor (MITF) gene for the spotting traits. Association studies were conducted for candidate regions or genes known to affect pigmentation in dairy cattle.

Detection and assessment of copy number variation using PacBio long read and Illumina sequencing in New Zealand dairy cattle

Single nucleotide polymorphisms have been the DNA variant of choice for genomic prediction, largely because of the ease of single nucleotide polymorphism genotype collection. In contrast, structural variants (SV), which include copy number variants (CNV), translocations, insertions, and inversions, have eluded easy detection and characterization, particularly in nonhuman species. However, evidence increasingly shows that SV not only contribute a substantial proportion of genetic variation but also have significant influence on phenotypes. Here we present the discovery of CNV in a prominent New Zealand dairy bull using long-read PacBio (Pacific Biosciences, Menlo Park, CA) sequencing technology and the Sniffles SV discovery tool (version 0.0.1; https://github.com/fritzsedlazeck/Sniffles). The CNV identified from long reads were compared with CNV discovered in the same bull from Illumina sequencing using CNVnator (read depth-based tool; Illumina Inc., San Diego, CA) as a means of validation. Subsequently, further validation was undertaken using whole-genome Illumina sequencing of 556 cattle representing the wider New Zealand dairy cattle population. Very limited overlap was observed in CNV discovered from the 2 sequencing platforms, in part because of the differences in size of CNV detected. Only a few CNV were therefore able to be validated using this approach. However, the ability to use CNVnator to genotype the 557 cattle for copy number across all regions identified as putative CNV allowed a genome-wide assessment of transmission level of copy number based on pedigree. The more highly transmissible a putative CNV region was observed to be, the more likely the distribution of copy number was multimodal across the 557 sequenced animals. Furthermore, visual assessment of highly transmissible CNV regions provided evidence supporting the presence of CNV across the sequenced animals. This transmission-based approach was able to confirm a subset of CNV that segregates in the New Zealand dairy cattle population. Genome-wide identification and validation of CNV is an important step toward their inclusion in genomic selection strategies.

DNA and RNA-sequence based GWAS highlights membrane-transport genes as key modulators of milk lactose content

BackgroundLactose 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.ResultsUsing 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.ConclusionsOverall, 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.

Multiple QTL underlie milk phenotypes at the CSF2RB locus

Background Bovine milk provides an important source of nutrition in much of the Western world, forming components of many food products. Over many years, artificial selection has substantially improved milk production by cows. However, the genes underlying milk production quantitative trait loci (QTL) remain relatively poorly characterised. Here, we investigate a previously-reported QTL located at the CSF2RB locus, for several milk production phenotypes, to better understand its underlying genetic and molecular causes. Results Using a population of 29,350 taurine dairy cattle, 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 eleven 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 the cause of 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 implicating the CSF2RB gene. Conclusions This chromosome 5 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.

Mating strategies to maximize genetic merit in dairy cattle herds

The genetic merit of a herd is a key determinant in productivity for dairy farmers. However, making breeding decisions to maximize the rate of genetic gain can be complex because there is no certainty about which cows will become pregnant with a heifer calf. In this study, breeding worth (BrW) was used as a measure of genetic merit, and several mating strategies were evaluated. These strategies included randomly mating whole herds to the entire bull team, excluding low-ranked cows from producing replacement heifers, and nominating high-ranked cows to the most highly ranked bulls. Simulations were undertaken using 4 bull teams generated from bulls currently marketed in New Zealand and a selection of New Zealand dairy herds. Average replacement heifer BrW was calculated for 1,000 iterations of each combination of mating strategy, herd, and bull team (scenario). Variation in resulting average replacement heifer BrW within scenarios was due to random sampling of which cows became pregnant with a heifer calf. Relative to mating the whole herd to an entire bull team, excluding the lowest ranked cows from producing replacements resulted in the greatest increase in average replacement heifer BrW across all herds and bull teams, with a gain of approximately 0.4 BrW point for each 1% of cows excluded. Nominating top-ranking cows to the highest ranking bulls in the team had little effect (0.06-0.13 BrW increase for each 1% of top cows nominated) in improving BrW of replacement heifers. The number of top bulls nominated had a variable effect depending on the BrW spread of the entire bull team. Although excluding cows with the lowest BrW from producing replacement heifers is most effective for improving BrW, it is important to ensure that the number of heifers born is sufficient to replace cows leaving the herd. It is likely that optimal strategies for improving BrW will vary from farm to farm depending not only on the BrW structure of the herd, the bull team available, and the reproduction success on farm but also on farm management practices. This simulation study provides expected outcomes from a variety of mating strategies to allow informed decision making on farm.

Bovine mammary gland X chromosome inactivation

X chromosome inactivation (XCI) is a process by which 1 of the 2 copies of the X chromosomes present in female mammals is inactivated. The transcriptional silencing of one X chromosome achieves dosage compensation between XX females and XY males and ensures equal expression of X-linked genes in both sexes. Although all mammals use this form of dosage compensation, the complex mechanisms that regulate XCI vary between species, tissues, and development. These mechanisms include not only varying levels of inactivation, but also the nature of inactivation, which can range from being random in nature to driven by parent of origin. To date, no data describing XCI in calves or adult cattle have been reported and we are reliant on data from mice to infer potential mechanisms and timings for this process. In the context of dairy cattle breeding and genomic prediction, the implications of X chromosome inheritance and XCI in the mammary gland are particularly important where a relatively small number of bulls pass their single X chromosome on to all of their daughters. We describe here the use of RNA-seq, whole genome sequencing and Illumina BovineHD BeadChip (Illumina, San Diego, CA) genotypes to assess XCI in lactating mammary glands of dairy cattle. At a population level, maternally and paternally inherited copies of the X chromosome are expressed equally in the lactating mammary gland consistent with random inactivation of the X chromosome. However, average expression of the paternal chromosome ranged from 10 to 90% depending on the individual animal. These results suggest that either the mammary gland arises from 1 or 2 stem cells, or a nongenetic mechanism that skews XCI exists. Although a considerable amount of future work is required to fully understand XCI in cattle, the data reported here represent an initial step in ensuring that X chromosome variation is captured and used in an appropriate manner for future genomic selection.