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

Pangenome genotyped structural variation improves molecular phenotype mapping in cattle

RNA sequencing variants are enriched for eQTL in cattle tissues

Large-effect pleiotropic or closely linked QTL segregate within and across ten US cattle breeds

High-resolution structural variants catalogue in a large-scale whole genome sequenced bovine family cohort data

Assembly of a pangenome for global cattle reveals missing sequences and novel structural variations, providing new insights into their diversity and evolutionary history

Epigenomics and Genotype-Phenotype Association Analyses Reveal Conserved Genetic Architecture of Complex Traits in Cattle and Human

Large-Scale in Silico Mapping of Complex Quantitative Traits in Inbred Mice

Identifying Pleiotropic Variants and Candidate Genes for Fertility and Reproduction Traits in Holstein Cattle Via Association Studies Based on Imputed Whole-Genome Sequence Genotypes

Expression quantitative trait loci in sheep liver and muscle contribute to variations in meat traits.

Allele-specific binding variants causing ChIP-seq peak height of histone modification are not enriched in expression QTL annotations

Expanded methylome and quantitative trait loci detection by long-read profiling of personal DNA

The eQTL colocalization and transcriptome-wide association study identify potentially causal genes responsible for economic traits in Simmental beef cattle

A whole-genome scan to map quantitative trait loci for conformation and functional traits in Canadian Holstein bulls.

SVLearn: a Dual-Reference Machine Learning Approach Enables Accurate Cross-Species Genotyping of Structural Variants

A large structural variant collection in Holstein cattle and associated database for variant discovery, characterization, and application

Comprehensive Analyses of 723 Transcriptomes Enhance Genetic and Biological Interpretations for Complex Traits in Cattle.

Genome-wide Detection of Copy Number Variations Using High-Density SNP Genotyping Platforms in Holsteins

Integrating CNVs into meta-QTL identified GBP4 as positional candidate for adult cattle stature

A Multi-Tissue Atlas of Regulatory Variants in Cattle

The first sheep graph-based pan-genome reveals the spectrum of structural variations and their effects on tail phenotypes

Molecular quantitative trait loci in reproductive tissues impact male fertility in cattle