Background: Insights into the genetic basis of complex traits and disease in both human and livestock species have been achieved over the past decade through detection of genetic variants in genome-wide association studies (GWAS). A majority of such variants were found located in noncoding genomic regions, and though the involvement of numerous regulatory elements (REs) has been predicted across multiple tissues in domesticated animals, their evolutionary conservation and effects on complex traits have not been fully elucidated, particularly in ruminants. Here, we systematically analyzed 137 epigenomic and transcriptomic datasets of six mammals, including cattle, sheep, goats, pigs, mice, and humans, and then integrated them with large-scale GWAS of complex traits.Results: Using 40 ChIP-seq datasets of H3K4me3 and H3K27ac, we detected 68,479, 58,562, 63,273, 97,244, 111,881, and 87,049 REs in the liver of cattle, sheep, goats, pigs, humans and mice, respectively. We then systematically characterized the dynamic functional landscapes of these REs by integrating multi-omics datasets, including gene expression, chromatin accessibility, and DNA methylation. We identified a core set (n= 6359) of ruminant-specific REs that are involved in liver development, metabolism, and immune processes. Genes with more complex cis-REs exhibited higher gene expression levels and stronger conservation across species. Furthermore, we integrated expression quantitative trait loci (eQTLs) and GWAS from 44 and 52 complex traits/diseases in cattle and humans, respectively. These results demonstrated that REs with different degrees of evolutionary conservation across species exhibited distinct enrichments for GWAS signals of complex traits.Conclusions: We systematically annotated genome-wide functional REs in liver across six mammals and demonstrated the evolution of REs and their associations with transcriptional output and conservation. Detecting lineage-specific REs allows us to decipher the evolutionary and genetic basis of complex phenotypes in livestock and humans, which may benefit the discovery of potential biomedical models for functional variants and genes of specific human diseases.

Objective: RNA epigenetic modifications play an important role in regulating immune response of mammals. Bovine mastitis induced by Staphylococcus aureus (S. aureus) is a threat to the health of dairy cattle. There are numerous RNA modifications, and how these modification-associated enzymes systematically coordinate their immunomodulatory effects during bovine mastitis is not well reported. Therefore, the role of common RNA modificationrelated genes (RMRGs) in bovine S. aureus mastitis was investigated in this study.Methods: In total, 80 RMRGs were selected for this study. Four public RNA-seq data sets about bovine S. aureus mastitis were collected and one additional RNA-seq data set was generated by this study. Firstly, quantitative trait locus (QTL) database, transcriptome-wide association studies (TWAS) database and differential expression analyses were employed to characterize the potential functions of selected enzyme genes in bovine S. aureus mastitis. Correlation analysis and weighted gene co-expression network analysis (WGCNA) were used to further investigate the relationships of RMRGs from different types at the mRNA expression level. Interference experiments targeting the m6A demethylase FTO and utilizing public MeRIP-seq dataset from bovine Mac-T cells were used to investigate the potential interaction mechanisms among various RNA modifications.Results: Bovine QTL and TWAS database in cattle revealed associations between RMRGs and immune-related complex traits. S. aureus challenged and control groups were effectively distinguished by principal component analysis based on the expression of selected RMRGs. WGCNA and correlation analysis identified modules grouping different RMRGs, with highly correlated mRNA expression. The m6A modification gene FTO showed significant effects on the expression of m6A and other RMRGs (such as NSUN2, CPSF2, and METTLE), indicating complex co-expression relationships among different RNA modifications in the regulation of bovine S. aureus mastitis.Conclusion: RNA epigenetic modification genes play important immunoregulatory roles in bovine S. aureus mastitis, and there are extensive interactions of mRNA expression among different RMRGs. It is necessary to investigate the interactions between RNA modification genes regulating complex traits in the future.

Somatic cell count (SCC) is an important indicator of the health state of bovine udders. However, the exact cut-off value used for differentiating the cows with healthy quarters from the cows with subclinical mastitis remains controversial. Here, we collected composite milk (milk from four udder quarters) and peripheral blood samples from individual cows in two different dairy farms and used 16S rRNA gene sequencing combined with RNA-seq to explore the differences in the milk microbial composition and transcriptome of cows with three different SCC levels (LSCC: 200,000 cells/mL). Results showed that the milk microbial profiles and gene expression profiles of samples derived from cows in the MSCC group were indeed relatively easily discriminated from those from cows in the LSCC group. Discriminative analysis also uncovered some differentially abundant microbiota at the genus level, such as Bifidobacterium and Lachnospiraceae_AC2044_group, which were more abundant in milk samples from cows with SCC below 100,000 cells/mL. As for the transcriptome profiling, 79 differentially expressed genes (DEGs) were found to have the same direction of regulation in two sites, and functional analyses also showed that biological processes involved in inflammatory responses were more active in MSCC and HSCC cows. Overall, these results showed a similarity between the milk microbiota and gene expression profiles of MSCC and HSCC cows, which presented further evidence that 100,000 cells/ml is a more optimal cut-off value than 200,000 cells/mL for intramammary infection detection at the cow level.

Objective: Abnormally increased somatic cell counts (SCCs) in milk is usually a sign of bovine subclinical mastitis. Mutual interaction between the host and its associated microbiota plays an important role in developing such diseases. The main objective of this study was to explore the difference between cows with elevated SCCs and healthy cattle from the perspective of host-microbe interplay.Methods: A total of 31 milk samples and 23 bovine peripheral blood samples were collected from Holstein dairy cattle to conduct an integrated analysis of transcriptomic and metagenomics.Results: The results showed that Ralstonia and Sphingomonas were enriched in cows with subclinical mastitis. The relative abundance of the two bacteria was positively correlated with the expression level of bovine transcobalamin 1 and uridine phosphorylase 1 encoding gene. Moreover, functional analysis revealed a distinct alternation in some important microbial biological processes.Conclusion: These results reveal the relative abundance of Ralstonia and Sphingomonas other than common mastitis-causing pathogens varied from healthy cows to those with subclinical mastitis and might be associated with elevated SCCs. Potential association was observed between bovine milk microbiota composition and the transcriptional pattern of some genes, thus providing new insights to understand homeostasis of bovine udder.

Mastitis caused by Staphylococcus aureus (S. aureus) infection is one of the most difficult diseases to treat in dairy cattle. Exploring the biological progression of S. aureus mastitis via the interaction between host, pathogen, and environment is the key to an effective and sustainable improvement of animal health. Here, two strains of S. aureus and a strain of MRSA (Methicillin-resistant Staphylococcus aureus) isolated from cows with different inflammation phenotypes were used to challenge Mac-T cells and to investigate their effects on the global transcriptome of the cells, then to explore the potential regulatory mechanisms of folic acid on S. aureus mastitis prevention. Differential gene expression or splicing analysis showed that different strains of S. aureus led to distinct transcriptional responses from the host immune system. Folic acid could protect host defense against the challenge of S. aureus and MRSA partially through activating cytoplasmic DNA sensing and tight junction pathway. ZBP1 at the upstream of cytoplasmic DNA sensing pathway was verified and related to anti-pathogen through RNA interference. Further enrichment analysis using these transcriptome data with cattle large-scale genome-wide association study (GWAS) data confirmed that ZBP1 gene is highly associated with bovine somatic cell score (SCS) trait. Our data shed light on the potential effect of FA through regulating key gene and then protect host cells’ defense against S. aureus and MRSA.

Chinese Simmental cattle serve as an important dual-purpose breed in sustainable livestock systems. Despite their economic value, the genetic architecture underlying milk-production traits in this breed under temperate conditions remains poorly characterized. In this study, we estimated genetic parameters and identified associated genomic loci for 9 milk-production traits in a Chinese Simmental population. Our dataset consisted of 17,556 test-day records from 1,788 cows (parities 1 to 3), including whole-genome sequencing data for 781 individuals. Using a random regression test-day model with Legendre polynomials, we estimated variance components, heritabilities, breeding values (EBVs), and genetic correlations between different DIM. Heritability estimates ranged from 0.09 (fat-to-protein ratio) to 0.52 (protein percentage), with intermediate values for fat percentage (0.28), lactose percentage (0.35), and total milk solids (0.35). Within-lactation genetic correlations varied by trait, ranging from -0.88 to 0.99, with the strongest correlations between adjacent DIM, which weakened as intervals widened. Genome-wide association analysis (using the mixed linear model in the SLEMM package) and Bayesian fine-mapping analysis identified significant SNPs near known candidate genes (BOP1, MROH1, NEGR1). These analyses also revealed putative novel associations with CACNB4, MTHFD2L, and SGMS1. Overall, the results enhance our understanding of the genetic regulation of milk production and provide practical targets for genomic selection to improve dairy performance in dual-purpose breeding programs.