Genome-wide association studies (GWAS) have identified many quantitative trait loci (QTL) associated with complex traits, predominantly in non-coding regions, posing challenges in pinpointing the causal variants and their target genes. Three types of evidence can help identify the gene through which QTL acts: (1) proximity to the most significant GWAS variant, (2) correlation of gene expression with the trait, and (3) the gene's physiological role in the trait. However, there is still uncertainty about the success of these methods in identifying the correct genes. Here, we test the ability of these methods in a comparatively simple series of traits associated with the concentration of polar lipids in milk. We conducted single-trait GWAS for ~14 million imputed variants and 56 individual milk polar lipid (PL) phenotypes in 336 cows. A multi-trait meta-analysis of GWAS identified 10,063 significant SNPs at FDR ≤ 10% (P ≤ 7.15E-5). Transcriptome data from blood (~12.5K genes, 143 cows) and mammary tissue (~12.2K genes, 169 cows) were analyzed using the genetic score omics regression (GSOR) method. This method links observed gene expression to genetically predicted phenotypes and was used to find associations between gene expression and 56 PL phenotypes. GSOR identified 2,186 genes in blood and 1,404 in mammary tissue associated with at least one PL phenotype (FDR ≤ 1%). We partitioned the genome into non-overlapping windows of 100 Kb to test for overlap between GSOR-identified genes and GWAS signals. We found a significant overlap between these two datasets, indicating that GSOR-significant genes were more likely to be located within 100 Kb windows that include GWAS signals than those that do not (P = 0.01; odds ratio = 1.47). These windows included 70 significant genes expressed in mammary tissue and 95 in blood. Compared to all expressed genes in each tissue, these genes were enriched for lipid metabolism gene ontology (GO). That is, seven of the 70 significant mammary transcriptome genes (P < 0.01; odds ratio = 3.98) and five of the 95 significant blood genes (P < 0.10; odds ratio = 2.24) were involved in lipid metabolism GO. The candidate causal genes include DGAT1, ACSM5, SERINC5, ABHD3, CYP2U1, PIGL, ARV1, SMPD5, and NPC2, with some overlap between the two tissues. The overlap between GWAS, GSOR, and GO analyses suggests that together, these methods are more likely to identify genes mediating QTL, though their power remains limited, as reflected by modest odds ratios. Larger sample sizes would enhance the power of these analyses, but issues like linkage disequilibrium would remain.

Segmental axonopathy is a recessively inherited neurodegenerative disorder that has affected Merino sheep since the early 1930s. Despite its long-standing recognition, the genetic basis of the condition remained unknown. This study aimed to identify the genetic cause of segmental axonopathy and confirm its pathological features to improve diagnostic accuracy and inform breeding strategies. Whole genome sequencing and genotyping of affected and unaffected Merino sheep identified a novel homozygous frameshift variant in the ALS2 gene that segregated with disease. RNA sequencing of cerebellar peduncle tissue confirmed the nonsense consequence on the ALS2 transcript. Histological analysis highlighted the hallmarks of the disease as large, foamy eosinophilic axonal swellings predominantly in the trigeminal ganglia, with additional degenerative changes in both the brain and spinal cord. These findings support the value of targeted sampling of sensory roots of the trigeminal nerve, spinal cord tracts, and dorsal nerve rootlets to enhance diagnostic accuracy. The same ALS2 variant was found across multiple unrelated flocks in both Australia and New Zealand, indicating a broader presence within the fine-wool Merino sheep population. This study identifies a novel ALS2 frameshift variant associated with segmental axonopathy in Merino sheep and provides both genetic and histological evidence supporting its role in disease pathology. The development of a DNA diagnostic test will enable more informed breeding decisions, reduce the prevalence of this condition, and improve animal welfare and productivity in the Merino industry. Moreover, the findings offer a potential large-animal model for exploring early-onset forms of human motor neuron diseases, including amyotrophic lateral sclerosis, in which ALS2 variants are implicated.

Background Genome sequence variants affecting complex traits (quantitative trait loci, QTL) are enriched in functional regions of the genome, such as those marked by certain histone modifications. These variants are believed to influence gene expression. However, due to the linkage disequilibrium among nearby variants, pinpointing the precise location of QTL is challenging. We aimed to identify allele-specific binding (ASB) QTL (asbQTL) that cause variation in the level of histone modification, as measured by the height of peaks assayed by ChIP-seq (chromatin immunoprecipitation sequencing). We identified DNA sequences that predict the difference between alleles in ChIP-seq peak height in H3K4me3 and H3K27ac histone modifications in the mammary glands of cows. Results We used a gapped k-mer support vector machine, a novel best linear unbiased prediction model, and a multiple linear regression model that combines the other two approaches to predict variant impacts on peak height. For each method, a subset of 1000 sites with the highest magnitude of predicted ASB was considered as candidate asbQTL. The accuracy of this prediction was measured by the proportion where the predicted direction matched the observed direction. Prediction accuracy ranged between 0.59 and 0.74, suggesting that these 1000 sites are enriched for asbQTL. Using independent data, we investigated functional enrichment in the candidate asbQTL set and three control groups, including non-causal ASB sites, non-ASB variants under a peak, and SNPs (single nucleotide polymorphisms) not under a peak. For H3K4me3, a higher proportion of the candidate asbQTL were confirmed as ASB when compared to the non-causal ASB sites (P < 0.01). However, these candidate asbQTL did not enrich for the other annotations, including expression QTL (eQTL), allele-specific expression QTL (aseQTL) and sites conserved across mammals (P > 0.05). Conclusions We identified putatively causal sites for asbQTL using the DNA sequence surrounding these sites. Our results suggest that many sites influencing histone modifications may not directly affect gene expression. However, it is important to acknowledge that distinguishing between putative causal ASB sites and other non-causal ASB sites in high linkage disequilibrium with the causal sites regarding their impact on gene expression may be challenging due to limitations in statistical power.

With decreasing wool values, interest is increasing regarding shedding sheep. To investigate this, two long-term studies introducing Wiltshire genes into Romney flocks were initiated. Data from these two studies provide phenotypic relationships between a range of shedding scores at different ages. The data included shedding scores (on a 0-5 scale) repeated on lambs (similar to 5 months), hoggets (similar to 14-18 months) and two-tooths (similar to 27 months), and lamb fleece weights. Positive relationships between shedding scores on the same animals were observed. Lamb fleece weight was negatively correlated with all shedding scores. Lamb shedding score in February had a correlation of 0.54 (P < 0.001) with the February score as a hogget at Riverside farm. Scoring wool shedding is a laborious activity requiring individual animals to be scored in the shearing position. Therefore, lamb fleece weight was investigated for its relationship with shedding scores, as a potentially easier alternative. Lamb fleece weight had a greater correlation with February hogget shedding score than with the lamb shedding score (-0.76 vs -0.52, P < 0.001). This study indicated that February scores are an accurate predictor of future shedding phenotypes, and when used in conjunction with fleece weight, are a good predictor of phenotypes expressed at later ages.

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.

In vitro investigation of bovine lactation processes is limited by a lack of physiologically representative cell models. This deficiency is most evident through the minimal or absent expression of lactation-specific genes in cultured bovine mammary tissues. Primary bovine mammary epithelial cells (pbMECs) extracted from lactating mammary tissue and grown in culture initially express milk protein transcripts at relatively representative levels. However, expression drops dramatically after only three or four passages, which greatly reduces the utility of primary cells to model and further examine lactogenesis. To investigate the effects of alternate alleles in pbMECs including effects on transcription, we have developed methods to deliver CRISPR-Cas9 gene editing reagents to primary mammary cells, resulting in very high editing efficiencies. We have also found that culturing the cells on an imitation basement membrane composed of Matrigel, results in the restoration of a more representative lactogenic gene expression profile and the cells forming three-dimensional structures in vitro. Here, we present data from four pbMEC lines recovered from pregnant cows and detail the expression profile of five key milk synthesis genes in these MECs grown on Matrigel. Additionally, we describe an optimised method for preferentially selecting CRISPR-Cas9-edited cells conferring a knock-out of DGAT1, using fluorescence-activated cell sorting (FACS). The combination of these techniques facilitates the use of pbMECs as a model to investigate the effects of gene introgressions and genetic variation in lactating mammary tissue.

Cerebellar hypoplasia is a heterogeneous neurological condition in which the cerebellum is smaller than usual or not completely developed. The condition can have genetic origins, with Mendelian-effect mutations described in several mammalian species. Here, we describe a genetic investigation of cerebellar hypoplasia in White Swiss Shepherd dogs, where two affected puppies were identified from a litter with a recent common ancestor on both sides of their pedigree. Whole genome sequencing was conducted for 10 dogs in this family, and filtering of these data based on a recessive transmission hypothesis highlighted five protein-altering candidate variants - including a frameshift-deletion of the Reelin (RELN) gene (p.Val947*). Given the status of RELN as a gene responsible for cerebellar hypoplasia in humans, sheep and mice, these data strongly suggest the loss-of-function variant as underlying these effects. This variant has not been found in other dog breeds nor in a cohort of European White Swiss Shepherds, suggesting a recent mutation event. This finding will support the genotyping of a more diverse sample of dogs, and should aid future management of the harmful allele through optimised mating schemes.

We report the generation of an organism-wide catalog of 976,813 cis-acting regulatory elements for the bovine detected by the assay for transposase accessible chromatin using sequencing (ATAC-seq). We regroup these regulatory elements in 16 components by nonnegative matrix factorization. Correlation between the genome-wide density of peaks and transcription start sites, correlation between peak accessibility and expression of neighboring genes, and enrichment in transcription factor binding motifs support their regulatory potential. Using a previously established catalog of 12,736,643 variants, we show that the proportion of single-nucleotide polymorphisms mapping to ATAC-seq peaks is higher than expected and that this is owing to an approximately 1.3-fold higher mutation rate within peaks. Their site frequency spectrum indicates that variants in ATAC-seq peaks are subject to purifying selection. We generate eQTL data sets for liver and blood and show that variants that drive eQTL fall into liver- and blood-specific ATAC-seq peaks more often than expected by chance. We combine ATAC-seq and eQTL data to estimate that the proportion of regulatory variants mapping to ATAC-seq peaks is approximately one in three and that the proportion of variants mapping to ATAC-seq peaks that are regulatory is approximately one in 25. We discuss the implication of these findings on the utility of ATAC-seq information to improve the accuracy of genomic selection.

Differentiating causative variants from linked variants is a major challenge of modern genetics. Here, we describe development of a massively parallel reporter assay (MPRA), an in vitro method that enables simultaneous functional testing of thousands of candidate causative regulatory variants. Candidates for the MPRA were identified through analysis of two large molecular quantitative trait loci (QTL) datasets derived from sequencing of bovine mammary biopsies, comprising an expression QTL resource (n=371 cows), and a histone QTL dataset (n=99 cows). Candidates that overlapped with lactation QTL discovered from genome-wide association studies were further prioritised to yield a final library of 242,340 oligos representing 24,234 regulatory sequences for expression testing. This paper gives a brief description of the MPRA method, our candidate selection criteria, cloning approach, and the composition of the final oligo pool for expression analysis.

Additional file 2: Table S1. GWAS results for the top 250 SNPs in the region of the BLG gene that were associated with total BLG concentration in milk (n = 4140). Two SNPs that differentiate the A and B protein isoforms are highlighted. Table S2. GWAS results for the top 250 SNPs in the region of the BLG gene that were associated with BLG B isoform concentrations in milks from heterozygous AB animals (n = 1966). The top SNP (highlighted) is B’ (chr11: 103256256G > A).

Charcot-Marie-Tooth disease (CMT) is a hereditary sensory and motor peripheral neuropathy that is one of the most common inherited neurological diseases of humans and may be caused by mutations in a number of different genes. The subtype Charcot-Marie-Tooth disease type 4H (CMT4H) is caused by homozygous mutations in the FGD4 (FYVE, RhoGEF, and PH domain-containing 4) gene. A previous genome-wide association study involving 130,783 dairy cows found 6 novel variants, one of which was a homozygous splice site mutation in the FGD4 gene. Descendants of carriers were genotyped to identify 9 homozygous Holstein Friesian calves that were raised to maturity, of which 5 were euthanized and sampled for histopathology and electron microscopy at 2 and 2.5 years of age. Three control Holstein Friesian animals were raised with the calves and euthanized at the same time points. No macroscopic lesions consistent with CMT4H were seen at necropsy. Microscopically, peripheral nerves were hypercellular due to hyperplasia of S100-positive Schwann cells, and there was onion bulb formation, axonal degeneration with demyelination, and increased thickness of the endoneurium. On electron microscopy, decreased axonal density, onion bulb formations, myelin outfoldings, and increased numbers of mitochondria were present. These changes are consistent with those described in mouse models and humans with CMT4H, making these cattle a potential large animal model for CMT.

Hereford bulls are being used in dairy farming systems to sire calves that can be finished as beef cattle, bred to cows that would not otherwise be used for replacement cows. Their characteristically white faces offer simple markers for parentage, however we observed that some dairy × Hereford crossbred calves have splotchy faces. We conducted a genome-wide association analysis for face colour in 128 Hereford-cross calves using genotype data from the GMK50kv1 chip and found that the adulteration of the white face trait is caused by an epistatic interaction between mutations at the KIT and MITF loci. We found that the one copy of the G allele at rs209784468 was sufficient to rescue pigmentation in the face in Hereford-cross calves. These results present a novel gene × gene interaction and suggest that the white face trait is not technically dominant.

Causal variants are enriched in functional non-coding regions of the genome, where it is hypothesised that they influence regulation of gene expression. Histone modifications are thought to mark many noncoding functional regions. Additionally, the height of some histone modification peaks has been found to be correlated with nearby gene expression, demonstrating a potential causal relationship between the two. Therefore, SNP which associate with differences in gene expression and with differences in peak height could be strong candidate causal variants. We found that QTL for peak height for three histone modifications (H3K4Me1, H3K4Me3 and H3K27ac) in ~100 mammary samples in dairy cows overlapped substantially with QTL for gene expression in ~400 mammary samples, thus identifying putative regulatory variants in the bovine genome.

Additional file 2 Supplementary 2 Summary of the 340 whole genome sequence variants identified from 302 goats located within the chromosome 19 QTL interval with high linkage disequilibrium, R2 > 0.8 with the most significant SNP marker rs268292132

Abstract 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.

Background Causal variants for complex traits, such as eQTL are often found in non-coding regions of the genome, where they are hypothesised to influence phenotypes by regulating gene expression. Many regulatory regions are marked by histone modifications, which can be assayed by chromatin immunoprecipitation followed by sequencing (ChIP-seq). Sequence reads from ChIP-seq form peaks at putative regulatory regions, which may reflect the amount of regulatory activity at this region. Therefore, eQTL which are also associated with differences in histone modifications are excellent candidate causal variants. Results We assayed the histone modifications H3K4Me3, H3K4Me1 and H3K27ac and mRNA in the mammary gland of up to 400 animals. We identified QTL for peak height (histone QTL), exon expression (eeQTL), allele specific expression (aseQTL) and allele specific binding (asbQTL). By intersecting these results, we identify variants which may influence gene expression by altering regulatory regions of the genome, and may be causal variants for other traits. Lastly, we find that these variants are found in putative transcription factor binding sites, identifying a mechanism for the effect of many eQTL. Conclusions We find that allele specific and traditional QTL analysis often identify the same genetic variants and provide evidence that many eQTL are regulatory variants which alter activity at regulatory regions of the bovine genome. Our work provides methodological and biological updates on how regulatory mechanisms interplay at multi-omics levels.

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.

Background Economically important milk production traits including milk volume, milk fat and protein yield vary considerably across dairy goats in New Zealand. A significant portion of the variation is attributable to genetic variation. Discovery of genetic markers linked to milk production traits can be utilised to drive selection of high-performance animals. A previously reported genome wide association study across dairy goats in New Zealand identified a quantitative trait locus (QTL) located on chromosome 19. The most significantly associated single nucleotide polymorphism (SNP) marker for this locus is located at position 26,610,610 (SNP marker rs268292132). This locus is associated with multiple milk production traits including fat, protein and volume. The predicted effect of selection for the beneficial haplotype would result in an average production increase of 2.2 kg fat, 1.9 kg protein and 73.6 kg milk yield. An outstanding question was whether selection for the beneficial allele would co-select for any negative pleiotropic effects. An adverse relationship between milk production and udder health traits has been reported at this locus. Therefore, a genome wide association study was undertaken looking for loci associated with udder traits. Results The QTL and production associated marker rs268292132 was identified in this study to also be associated with several goat udder traits including udder depth (UD), fore udder attachment (FUA) and rear udder attachment (RUA). Our study replicates the negative relationship between production and udder traits with the high production allele at position 19:26,610,610 (SNP marker rs268292132) associated with an adverse change in UD, FUA and RUA. Conclusions Our study has confirmed the negative relationship between udder traits and production traits in the NZ goat population. We have found that the frequency of the high production allele is relatively high in the NZ goat population, indicating that its effect on udder conformation is not significantly detrimental on animal health. It will however be important to monitor udder conformation as the chromosome 19 locus is progressively implemented for marker assisted selection. It will also be of interest to determine if the gene underlying the production QTL has a direct effect on mammary gland morphology or whether the changes observed are a consequence of the increased milk volume.