CUX1 attenuates the apoptosis of renal tubular epithelial cells induced by contrast media through activating the PI3K/AKT signaling pathway.

OBJECTIVE: Contrast media (CM) is a commonly applied drug in medical examination and surgery. However, contrast-induced acute kidney injury (CIAKI) poses a severe threat to human life and health. Notably, the CUT-like homeobox 1 (CUX1) gene shows protective effects in a variety of cells. Therefore, the objective of this study was to provide a new target for the treatment of CIAKI through exploring the role and possible molecular mechanism of CUX1 in CIAKI. METHOD: Blood samples were collected from 20 patients with CIAKI and healthy volunteers. Human kidney 2 (HK-2) cells were incubated with 200 mg/mL iohexol for 6 h to establish a contrast-induced injury model of HK-2 cells. Subsequently, qRT-PCR was used to detect the relative mRNA expression of CUX1; CCK-8 and flow cytometry to assess the proliferation and apoptosis of HK-2 cells; the levels of IL(interleukin)-1ß, tumor necrosis factor alpha (TNF-α) and malondialdehyde (MDA) in cells and lactate dehydrogenase (LDH) activity in cell culture supernatant were detect; and western blot to observe the expression levels of CUX1 and the PI3K/AKT signaling pathway related proteins [phosphorylated phosphoinositide 3-kinase (p-PI3K), PI3K, phosphorylated Akt (p-AKT), AKT]. RESULTS: CUX1 expression was significantly downregulated in blood samples of patients with CIAKI and contrast-induced HK-2 cells. Contrast media (CM; iohexol) treatment significantly reduced the proliferation of HK-2 cells, promoted apoptosis, stimulated inflammation and oxidative stress that caused cell damage. CUX1 overexpression alleviated cell damage by significantly improving the proliferation level of HK-2 cells induced by CM, inhibiting cell apoptosis, and reducing the level of LDH in culture supernatant and the expression of IL-1ß, TNF-α and MDA in cells. CM treatment significantly inhibited the activity of PI3K/AKT signaling pathway activity. Nevertheless, up-regulating CUX1 could activate the PI3K/AKT signaling pathway activity in HK-2 cells induced by CM. CONCLUSION: CUX1 promotes cell proliferation, inhibits apoptosis, and reduces inflammation and oxidative stress in CM-induced HK-2 cells to alleviate CM-induced damage. The mechanism of CUX1 may be correlated with activation of the PI3K/AKT signaling pathway.

Comprehensive analyses of 723 transcriptomes enhance genetic and biological interpretations for complex traits in cattle.

By uniformly analyzing 723 RNA-seq data from 91 tissues and cell types, we built a comprehensive gene atlas and studied tissue specificity of genes in cattle. We demonstrated that tissue-specific genes significantly reflected the tissue-relevant biology, showing distinct promoter methylation and evolution patterns (e.g., brain-specific genes evolve slowest, whereas testis-specific genes evolve fastest). Through integrative analyses of those tissue-specific genes with large-scale genome-wide association studies, we detected relevant tissues/cell types and candidate genes for 45 economically important traits in cattle, including blood/immune system (e.g., CCDC88C) for male fertility, brain (e.g., TRIM46 and RAB6A) for milk production, and multiple growth-related tissues (e.g., FGF6 and CCND2) for body conformation. We validated these findings by using epigenomic data across major somatic tissues and sperm. Collectively, our findings provided novel insights into the genetic and biological mechanisms underlying complex traits in cattle, and our transcriptome atlas can serve as a primary source for biological interpretation, functional validation, studies of adaptive evolution, and genomic improvement in livestock.

Dietary intake of microplastics impairs digestive performance, induces hepatic dysfunction, and shortens lifespan in the annual fish Nothobranchius guentheri.

Microplastics (MPs) are ubiquitous in aquatic and terrestrial ecosystem, increasingly becoming a serious concern of human health. Many studies have explored the biological effects of MPs on animal and plant life in recent years. However, information regarding the effects of MPs on aging and lifespan is completely lacking in vertebrate species to date. Here we first confirm the bioavailability of MPs by oral delivery in the annual fish N. guentheri. We then show for the first time that administration of MPs not only shortens the lifespan but also accelerates the development of age-related biomarkers in N. guentheri. We also demonstrate that administration of MPs induces oxidative stress, suppresses antioxidant enzymes, reduces digestive enzymes, and causes hepatic dysfunction. Therefore, we propose that administration of MPs reduces lifespan of N. guentheri via induction of both suppressed antioxidant system and digestive disturbance as well as hepatic damage. Our results also suggest that smaller MPs appear more toxic to digestion, metabolism and growth of animals.

Atorvastatin reduces contrast media-induced pyroptosis of renal tubular epithelial cells by inhibiting the TLR4/MyD88/NF-κB signaling pathway.

BACKGROUND: Contrast-induced acute kidney injury (CI-AKI) is the third most common cause of hospital-acquired renal failure. However, there is no effective treatment of CI-AKI, and its mechanism is unknown. Interestingly, atorvastatin has been reported to be effective in renal injury. Therefore, the aim of this study was to explore the effect and possible molecular mechanism of atorvastatin in CI-AKI. METHODS: On the CI-AKI in vitro model, rat tubular epithelial cells (NRK-52E) were treated with 18 mg I/ml meglumine diatrizoate (MEG) and then pretreated with atorvastatin. pcDNA3.1-TLR4 treatment was performed to overexpress toll-like receptor 4 (TLR4) in NRK-52E cells. Cell Counting Kit-8 (CCK-8) and lactate dehydrogenase (LDH) kits were used to detect NRK-52E cell viability as well as LDH release in each group, respectively; qRT-PCR to determine mRNA expression of TLR4 in cells; western blot to detect protein expression levels of pyroptosis-related proteins (NLRP3, caspase-1, ASC, and GSDMD) and TLR4/MyD88/NF-κB signaling pathway-related proteins (TLR4, MyD88, NF-κBp65, and p-NF-κB p65) in cells. RESULTS: MEG treatment significantly inhibited the viability of NRK-52E cells, increased pro-inflammatory factor levels and promoted pyroptosis, representing successful establishment of a rat tubular epithelial cell (NRK-52E) CI-AKI in vitro model. Notably, atorvastatin increased the activity of MEG-treated NRK-52E cells and alleviated cell injury in a concentration-dependent manner. In addition, atorvastatin significantly down-regulated the expression of TLR4 in MEG-treated NRK-52E cells. However, overexpression of TLR4 inhibited the effects of atorvastatin on increasing cell viability, alleviating cell injury, reducing pro-inflammatory factors (IL-1ß, IL-6, and TNF-α) levels, and inhibiting apoptosis (by down-regulating the expression of NLRP3, caspase-1, ASC, and GSDMD). Furthermore, atorvastatin also inhibited the expression of TLR4/MyD88/NF-κB pathway-related proteins (TLR4, MyD88, and p-NF-κB p65). CONCLUSION: Atorvastatin can attenuate CI-AKI through increasing the activity of MEG-treated renal tubular epithelial cells, relieving cell injury, as well as inhibiting pyroptosis and inflammation. More importantly, the mechanism was achieved by inhibiting the TLR4//MyD88/NF-κB signaling pathway.

Functional annotation of regulatory elements in cattle genome reveals the roles of extracellular interaction and dynamic change of chromatin states in rumen development during weaning.

We profiled landscapes of bovine regulatory elements and explored dynamic changes of chromatin states in rumen development during weaning. The regulatory elements (15 chromatin states) and their coordinated activities in cattle were defined through genome-wide profiling of four histone modifications, CTCF-binding, DNA accessibility, DNA methylation, and transcriptome in rumen epithelial tissues. Each chromatin state presented specific enrichment for sequence ontology, methylation, trait-associated variants, transcription, gene expression-associated variants, selection signatures, and evolutionarily conserved elements. During weaning, weak enhancers and flanking active transcriptional start sites (TSS) were the most dynamic chromatin states and occurred in tandem with significant variations in gene expression and DNA methylation, significantly associated with stature, production, and reproduction economic traits. By comparing with in vitro cultured epithelial cells and in vivo rumen tissues, we showed the commonness and uniqueness of these results, especially the roles of cell interactions and mitochondrial activities in tissue development.

Sodium Butyrate Induces Mitophagy and Apoptosis of Bovine Skeletal Muscle Satellite Cells through the Mammalian Target of Rapamycin Signaling Pathway.

Sodium butyrate (NaB) is one of the short-chain fatty acids and is notably produced in large amounts from dietary fiber in the gut. Recent evidence suggests that NaB induces cell proliferation and apoptosis. Skeletal muscle is rich in plenty of mitochondrial. However, it is unclear how NaB acts on host muscle cells and whether it is involved in mitochondria-related functions in myocytes. The present study aimed to investigate the role of NaB treatment on the proliferation, apoptosis, and mitophagy of bovine skeletal muscle satellite cells (BSCs). The results showed that NaB inhibited proliferation, promoted apoptosis of BSCs, and promoted mitophagy in a time- and dose-dependent manner in BSCs. In addition, 1 mM NaB increased the mitochondrial ROS level, decreased the mitochondrial membrane potential (MMP), increased the number of autophagic vesicles in mitochondria, and increased the mitochondrial DNA (mtDNA) and ATP level. The effects of the mTOR pathway on BSCs were investigated. The results showed that 1 mM NaB inhibited the mRNA and protein expression of mTOR and genes AKT1, FOXO1, and EIF4EBP1 in the mTOR signaling pathway. In contrast, the addition of PP242, an inhibitor of the mTOR signaling pathway also inhibited mRNA and protein expression levels of mTOR, AKT1, FOXO1, and EIF4EBP1 and promoted mitophagy and apoptosis, which were consistent with the effect of NaB treatment. NaB might promote mitophagy and apoptosis in BSCs by inhibiting the mTOR signaling pathway. Our results would expand the knowledge of sodium butyrate on bovine skeletal muscle cell state and mitochondrial function.

Investigation of rumen long noncoding RNA before and after weaning in cattle.

BACKGROUND: This study aimed to identify long non-coding RNA (lncRNA) from the rumen tissue in dairy cattle, explore their features including expression and conservation levels, and reveal potential links between lncRNA and complex traits that may indicate important functional impacts of rumen lncRNA during the transition to the weaning period. RESULTS: A total of six cattle rumen samples were taken with three replicates from before and after weaning periods, respectively. Total RNAs were extracted and sequenced with lncRNA discovered based on size, coding potential, sequence homology, and known protein domains. As a result, 404 and 234 rumen lncRNAs were identified before and after weaning, respectively. However, only nine of them were shared under two conditions, with 395 lncRNAs found only in pre-weaning tissues and 225 only in post-weaning samples. Interestingly, none of the nine common lncRNAs were differentially expressed between the two weaning conditions. LncRNA averaged shorter length, lower expression, and lower conservation scores than the genome overall, which is consistent with general lncRNA characteristics. By integrating rumen lncRNA before and after weaning with large-scale GWAS results in cattle, we reported significant enrichment of both pre- and after-weaning lncRNA with traits of economic importance including production, reproduction, health, and body conformation phenotypes. CONCLUSIONS: The majority of rumen lncRNAs are uniquely expressed in one of the two weaning conditions, indicating a functional role of lncRNA in rumen development and transition of weaning. Notably, both pre- and post-weaning lncRNA showed significant enrichment with a variety of complex traits in dairy cattle, suggesting the importance of rumen lncRNA for cattle performance in the adult stage. These relationships should be further investigated to better understand the specific roles lncRNAs are playing in rumen development and cow performance.

Single-cell transcriptomic and chromatin accessibility analyses of dairy cattle peripheral blood mononuclear cells and their responses to lipopolysaccharide.

BACKGROUND: Gram-negative bacteria are important pathogens in cattle, causing severe infectious diseases, including mastitis. Lipopolysaccharides (LPS) are components of the outer membrane of Gram-negative bacteria and crucial mediators of chronic inflammation in cattle. LPS modulations of bovine immune responses have been studied before. However, the single-cell transcriptomic and chromatin accessibility analyses of bovine peripheral blood mononuclear cells (PBMCs) and their responses to LPS stimulation were never reported. RESULTS: We performed single-cell RNA sequencing (scRNA-seq) and single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) in bovine PBMCs before and after LPS treatment and demonstrated that seven major cell types, which included CD4 T cells, CD8 T cells, and B cells, monocytes, natural killer cells, innate lymphoid cells, and dendritic cells. Bioinformatic analyses indicated that LPS could increase PBMC cell cycle progression, cellular differentiation, and chromatin accessibility. Gene analyses further showed significant changes in differential expression, transcription factor binding site, gene ontology, and regulatory interactions during the PBMC responses to LPS. Consistent with the findings of previous studies, LPS induced activation of monocytes and dendritic cells, likely through their upregulated TLR4 receptor. NF-κB was observed to be activated by LPS and an increased transcription of an array of pro-inflammatory cytokines, in agreement that NF-κB is an LPS-responsive regulator of innate immune responses. In addition, by integrating LPS-induced differentially expressed genes (DEGs) with large-scale GWAS of 45 complex traits in Holstein, we detected trait-relevant cell types. We found that selected DEGs were significantly associated with immune-relevant health, milk production, and body conformation traits. CONCLUSION: This study provided the first scRNAseq and scATAC-seq data for cattle PBMCs and their responses to the LPS stimulation to the best of our knowledge. These results should also serve as valuable resources for the future study of the bovine immune system and open the door for discoveries about immune cell roles in complex traits like mastitis at single-cell resolution.

Towards the detection of copy number variation from single sperm sequencing in cattle.

BACKGROUND: Copy number variation (CNV) has been routinely studied using bulk-cell sequencing. However, CNV is not well studied on the single-cell level except for humans and a few model organisms. RESULTS: We sequenced 143 single sperms of two Holstein bulls, from which we predicted CNV events using 14 single sperms with deep sequencing. We then compared the CNV results derived from single sperms with the bulk-cell sequencing of one bull's family trio of diploid genomes. As a known CNV hotspot, segmental duplications were also predicted using the bovine ARS-UCD1.2 genome. Although the trio CNVs validated only some single sperm CNVs, they still showed a distal chromosomal distribution pattern and significant associations with segmental duplications and satellite repeats. CONCLUSION: Our preliminary results pointed out future research directions and highlighted the importance of uniform whole genome amplification, deep sequence coverage, and dedicated software pipelines for CNV detection using single cell sequencing data.

Genome-wide recombination map construction from single sperm sequencing in cattle.

BACKGROUND: Meiotic recombination is one of the important phenomena contributing to gamete genome diversity. However, except for human and a few model organisms, it is not well studied in livestock, including cattle. RESULTS: To investigate their distributions in the cattle sperm genome, we sequenced 143 single sperms from two Holstein bulls. We mapped meiotic recombination events at high resolution based on phased heterozygous single nucleotide polymorphism (SNP). In the absence of evolutionary selection pressure in fertilization and survival, recombination events in sperm are enriched near distal chromosomal ends, revealing that such a pattern is intrinsic to the molecular mechanism of meiosis. Furthermore, we further validated these findings in single sperms with results derived from sequencing its family trio of diploid genomes and our previous studies of recombination in cattle. CONCLUSIONS: To our knowledge, this is the first large-scale single sperm whole-genome sequencing effort in livestock, which provided useful information for future studies of recombination, genome instability, and male infertility.

Dietary intake of GDF11 delays the onset of several biomarkers of aging in male mice through anti-oxidant system via Smad2/3 pathway.

Current studies have generated controversy over the age-related change in concentration of growth differentiation factor 11 (GDF11) and its role in the genesis of rejuvenation conditions. In this study, we displayed rGDF11 on the surface of Yarrowic Lipolytica (Y. lipolytica), and proved the bioavailability of the yeast-displayed rGDF11 by oral delivery in aged male mice. On the basis of these findings, we started to explore the anti-aging activity and underlying mechanisms of displayed rGDF11. It was found that dietary intake of displayed rGDF11 had little influence on the body weight and biochemical parameters of aged male mice, but delayed the occurrence and development of age-related biomarkers such as lipofuscin (LF) and senescence-associated-ß-galactosidase, and to some extent, prolonged the lifespan of aged male mice. Moreover, we demonstrated once again that dietary intake of displayed rGDF11 enhanced the activity of anti-oxidant enzymes, including catalase (CAT), superoxide dismutase (SOD) and glutathione peroxidase (GPX), reduced the reactive oxygen species (ROS) level, and slowed down the protein oxidation and lipid peroxidation. Importantly, we showed for the first time that rGDF11 enhanced the activity of CAT, SOD and GPX through activation of the Smad2/3 signaling pathway. Our study also provided a simple and safe route for delivery of recombinant GDF11, facilitating its therapeutic application in the future.

Dietary intake of diosgenin delays aging of male fish Nothobranchius guentheri through modulation of multiple pathways that play prominent roles in ROS production.

Oxidative stress including DNA damage, increased lipid and protein oxidation, is an important feature of aging. Diosgenin (DG) has been shown to have diverse biological effects, including amelioration of aging-related cognition deficits, but the anti-aging activity of DG has not been tested before in animal models. In the present study, we clearly demonstrated that dietary intake of DG extended both mean and maximum lifespans of the male fish Nothobranchius guentheri by approximately 3.23 and 3.67 weeks, respectively, reduced the accumulation of lipofuscin (LF) in the gills and senescence-associated-ß-galactosidase (SA-ß-Gal) in the caudal fins, and lowered the levels of protein oxidation, lipid peroxidation and reactive oxygen species (ROS) in the muscles, indicating that DG possesses rejuvenation and anti-aging property. We also showed that DG enhanced the activity of antioxidant enzymes, including catalase, superoxide dismutase and glutathione peroxidase, promoted the proteolytic activity of the ubiquitin-proteasome pathway, and suppressed the phosphatidylinositol 3-kinase/protein kinase/molecular target of rapamycin (PI3K/AKT/mTOR) signaling pathway. Altogether, this study highlights for the first time the rejuvenation and anti-aging property of the naturally occurring steroidal sapogenin DG. It also suggests that DG exerts its rejuvenation and anti-aging activity through modulation of multiple signaling pathways that play prominent roles in ROS production.

Identification of Isthmin1 in the small annual fish, Nothobranchius guentheri, as a novel biomarker of aging and its potential rejuvenation activity.

Isthmin 1 (Ism1) has been shown to play roles in multiple biological processes including morphogenesis, hematopoiesis, antiviral immune response and suppression of tumor growth. However, it remains unknown if it plays any role in aging process. Here we showed for the first time that Ism1 was a new age-related biomarker, which decreased with age in fish, mice and humans. Interestingly, Ism1 was also useful to measure the "rejuvenated" age of fish Nothobranchius guentheri reversed by salidroside treatment and temperature reduction, providing additional evidence that Ism1 was an aging biomarker. In addition, we clearly showed that dietary intake of recombinant Ism1 had little effects on the body length and weight of aging N. guentheri, but it retarded the onset of age-related biomarkers and prolonged both the maximum and median lifespan of the fish. We also showed that Ism1 exerted its rejuvenation activity via the enhancement of antioxidant system. Collectively, our results indicate that Ism1 is not only is a novel biomarker of aging but also a potential rejuvenation factor capable of reversing aging of N. guentheri.

Comprehensive transcriptome and methylome analysis delineates the biological basis of hair follicle development and wool-related traits in Merino sheep.

BACKGROUND: Characterization of the molecular mechanisms underlying hair follicle development is of paramount importance in the genetic improvement of wool-related traits in sheep and skin-related traits in humans. The Merino is the most important breed of fine-wooled sheep in the world. In this study, we systematically investigated the complexity of sheep hair follicle development by integrating transcriptome and methylome datasets from Merino sheep skin. RESULTS: We analysed 72 sequence datasets, including DNA methylome and the whole transcriptome of four gene types, i.e. protein-coding genes (PCGs), lncRNAs, circRNAs, and miRNAs, across four embryonic days (E65, E85, E105, and E135) and two postnatal days (P7 and P30) from the skin tissue of 18 Merino sheep. We revealed distinct expression profiles of these four gene types across six hair follicle developmental stages, and demonstrated their complex interactions with DNA methylation. PCGs with stage-specific expression or regulated by stage-specific lncRNAs, circRNAs, and miRNAs were significantly enriched in epithelial differentiation and hair follicle morphogenesis. Regulatory network and gene co-expression analyses identified key transcripts controlling hair follicle development. We further predicted transcriptional factors (e.g. KLF4, LEF1, HOXC13, RBPJ, VDR, RARA, and STAT3) with stage-specific involvement in hair follicle morphogenesis. Through integrating these stage-specific genomic features with results from genome-wide association studies (GWAS) of five wool-related traits in 7135 Merino sheep, we detected developmental stages and genes that were relevant with wool-related traits in sheep. For instance, genes that were specifically upregulated at E105 were significantly associated with most of wool-related traits. A phenome-wide association study (PheWAS) demonstrated that candidate genes of wool-related traits (e.g. SPHK1, GHR, PPP1R27, CSRP2, EEF1A2, and PTPN1) in sheep were also significantly associated with dermatological, metabolic, and immune traits in humans. CONCLUSIONS: Our study provides novel insights into the molecular basis of hair follicle morphogenesis and will serve as a foundation to improve breeding for wool traits in sheep. It also indicates the importance of studying gene expression in the normal development of organs in understanding the genetic architecture of economically important traits in livestock. The datasets generated here are useful resources for functionally annotating the sheep genome, and for elucidating early skin development in mammals, including humans.

Single-cell transcriptomic analyses of dairy cattle ruminal epithelial cells during weaning.

Using the 10× Genomics Chromium Controller, we obtained scRNA-seq data of 5064 and 1372 individual cells from two Holstein calf ruminal epithelial tissues before and after weaning, respectively. We detected six distinct cell clusters, designated their cell types, and reported their marker genes. We then examined these clusters' underlining cell types and relationships by performing cell cycle, pseudotime trajectory, regulatory network, weighted gene co-expression network and gene ontology analyses. By integrating these cell marker genes with Holstein GWAS signals, we found they were enriched for animal production and body conformation traits. Finally, we confirmed their cell identities by comparing them with human and mouse stomach epithelial cells. This study presents an initial effort to implement single-cell transcriptomic analysis in cattle, and demonstrates ruminal tissue epithelial cell types and their developments during weaning, opening the door for new discoveries about tissue/cell type roles in complex traits at single-cell resolution.

Differentially CTCF-Binding Sites in Cattle Rumen Tissue during Weaning.

The weaning transition in calves is characterized by major structural changes such as an increase in the rumen capacity and surface area due to diet changes. Studies evaluating rumen development in calves are vital to identify genetic mechanisms affected by weaning. This study aimed to provide a genome-wide characterization of CTCF-binding sites and differentially CTCF-binding sites (DCBS) in rumen tissue during the weaning transition of four Holstein calves to uncover regulatory elements in rumen epithelial tissue using ChIP-seq. Our study generated 67,280 CTCF peaks for the before weaning (BW) and 39,891 for after weaning (AW). Then, 7401 DCBS were identified for the AW vs. BW comparison representing 0.15% of the cattle genome, comprising ~54% of induced DCBS and ~46% of repressed DCBS. Most of the induced and repressed DCBS were in distal intergenic regions, showing a potential role as insulators. Gene ontology enrichment revealed many shared GO terms for the induced and the repressed DCBS, mainly related to cellular migration, proliferation, growth, differentiation, cellular adhesion, digestive tract morphogenesis, and response to TGFß. In addition, shared KEGG pathways were obtained for adherens junction and focal adhesion. Interestingly, other relevant KEGG pathways were observed for the induced DCBS like gastric acid secretion, salivary secretion, bacterial invasion of epithelial cells, apelin signaling, and mucin-type O-glycan biosynthesis. IPA analysis further revealed pathways with potential roles in rumen development during weaning, including TGFß, Integrin-linked kinase, and Integrin signaling. When DCBS were further integrated with RNA-seq data, 36 putative target genes were identified for the repressed DCBS, including KRT84, COL9A2, MATN3, TSPAN1, and AJM1. This study successfully identified DCBS in cattle rumen tissue after weaning on a genome-wide scale and revealed several candidate target genes that may have a role in rumen development, such as TGFß, integrins, keratins, and SMADs. The information generated in this preliminary study provides new insights into bovine genome regulation and chromatin landscape.

Comparative whole genome DNA methylation profiling across cattle tissues reveals global and tissue-specific methylation patterns.

BACKGROUND: Efforts to improve animal health, and understand genetic bases for production, may benefit from a comprehensive analysis of animal genomes and epigenomes. Although DNA methylation has been well studied in humans and other model species, its distribution patterns and regulatory impacts in cattle are still largely unknown. Here, we present the largest collection of cattle DNA methylation epigenomic data to date. RESULTS: Using Holstein cattle, we generated 29 whole genome bisulfite sequencing (WGBS) datasets for 16 tissues, 47 corresponding RNA-seq datasets, and 2 whole genome sequencing datasets. We did read mapping and DNA methylation calling based on two different cattle assemblies, demonstrating the high quality of the long-read-based assembly markedly improved DNA methylation results. We observed large differences across cattle tissues in the methylation patterns of global CpG sites, partially methylated domains (PMDs), hypomethylated regions (HMRs), CG islands (CGIs), and common repeats. We detected that each tissue had a distinct set of PMDs, which showed tissue-specific patterns. Similar to human PMD, cattle PMDs were often linked to a general decrease of gene expression and a decrease in active histone marks and related to long-range chromatin organizations, like topologically associated domains (TADs). We tested a classification of the HMRs based on their distributions relative to transcription start sites (TSSs) and detected tissue-specific TSS-HMRs and genes that showed strong tissue effects. When performing cross-species comparisons of paired genes (two opposite strand genes with their TSS located in the same HMR), we found out they were more consistently co-expressed among human, mouse, sheep, goat, yak, pig, and chicken, but showed lower consistent ratios in more divergent species. We further used these WGBS data to detect 50,023 experimentally supported CGIs across bovine tissues and found that they might function as a guard against C-to-T mutations for TSS-HMRs. Although common repeats were often heavily methylated, some young Bov-A2 repeats were hypomethylated in sperm and could affect the promoter structures by exposing potential transcription factor binding sites. CONCLUSIONS: This study provides a comprehensive resource for bovine epigenomic research and enables new discoveries about DNA methylation and its role in complex traits.

Functional annotation of the cattle genome through systematic discovery and characterization of chromatin states and butyrate-induced variations.

BACKGROUND: The functional annotation of genomes, including chromatin accessibility and modifications, is important for understanding and effectively utilizing the increased amount of genome sequences reported. However, while such annotation has been well explored in a diverse set of tissues and cell types in human and model organisms, relatively little data are available for livestock genomes, hindering our understanding of complex trait variation, domestication, and adaptive evolution. Here, we present the first complete global landscape of regulatory elements in cattle and explore the dynamics of chromatin states in rumen epithelial cells induced by the rumen developmental regulator-butyrate. RESULTS: We established the first global map of regulatory elements (15 chromatin states) and defined their coordinated activities in cattle, through genome-wide profiling for six histone modifications, RNA polymerase II, CTCF-binding sites, DNA accessibility, DNA methylation, and transcriptome in rumen epithelial primary cells (REPC), rumen tissues, and Madin-Darby bovine kidney epithelial cells (MDBK). We demonstrated that each chromatin state exhibited specific enrichment for sequence ontology, transcription, methylation, trait-associated variants, gene expression-associated variants, selection signatures, and evolutionarily conserved elements, implying distinct biological functions. After butyrate treatments, we observed that the weak enhancers and flanking active transcriptional start sites (TSS) were the most dynamic chromatin states, occurred concomitantly with significant alterations in gene expression and DNA methylation, which was significantly associated with heifer conception rate and stature economic traits. CONCLUSION: Our results demonstrate the crucial role of functional genome annotation for understanding genome regulation, complex trait variation, and adaptive evolution in livestock. Using butyrate to induce the dynamics of the epigenomic landscape, we were able to establish the correlation among nutritional elements, chromatin states, gene activities, and phenotypic outcomes.

Analyses of inter-individual variations of sperm DNA methylation and their potential implications in cattle.

BACKGROUND: DNA methylation has been shown to be involved in many biological processes, including X chromosome inactivation in females, paternal genomic imprinting, and others. RESULTS: Based on the correlation patterns of methylation levels of neighboring CpG sites among 28 sperm whole genome bisulfite sequencing (WGBS) data (486 × coverage), we obtained 31,272 methylation haplotype blocks (MHBs). Among them, we defined conserved methylated regions (CMRs), variably methylated regions (VMRs) and highly variably methylated regions (HVMRs) among individuals, and showed that HVMRs might play roles in transcriptional regulation and function in complex traits variation and adaptive evolution by integrating evidence from traditional and molecular quantitative trait loci (QTL), and selection signatures. Using a weighted correlation network analysis (WGCNA), we also detected a co-regulated module of HVMRs that was significantly associated with reproduction traits, and enriched for glycosyltransferase genes, which play critical roles in spermatogenesis and fertilization. Additionally, we identified 46 VMRs significantly associated with reproduction traits, nine of which were regulated by cis-SNPs, implying the possible intrinsic relationships among genomic variations, DNA methylation, and phenotypes. These significant VMRs were co-localized (± 10 kb) with genes related to sperm motility and reproduction, including ZFP36L1, CRISP2 and HGF. We provided further evidence that rs109326022 within a predominant QTL on BTA18 might influence the reproduction traits through regulating the methylation level of nearby genes JOSD2 and ASPDH in sperm. CONCLUSION: In summary, our results demonstrated associations of sperm DNA methylation with reproduction traits, highlighting the potential of epigenomic information in genomic improvement programs for cattle.

Transcriptome analysis unraveled potential mechanisms of resistance to Haemonchus contortus infection in Merino sheep populations bred for parasite resistance.

Haemonchus contortus is one of the most pathogenic gastrointestinal nematodes in small ruminants. To understand molecular mechanisms underlying host resistance to this parasite, we used RNA-sequencing technology to compare the transcriptomic response of the abomasal tissue, the site of the host-parasite interaction, of Merino sheep bred to be either genetically resistant or susceptible to H. contortus infection. Two different selection flocks, the Haemonchus selection flock (HSF) and the Trichostrongylus selection flock (TSF), and each contains a resistant and susceptible line, were studied. The TSF flock was seemingly more responsive to both primary and repeated infections than HSF. A total of 127 and 726 genes displayed a significant difference in abundance between resistant and susceptible animals in response to a primary infection in HSF and TSF, respectively. Among them, 38 genes were significantly affected by infection in both flocks. Gene ontology (GO) enrichment of the differentially expressed genes identified in this study predicted the likely involvement of extracellular exosomes in the immune response to H. contortus infection. While the resistant lines in HSF and TSF relied on different mechanisms for the development of host resistance, adhesion and diapedesis of both agranulocytes and granulocytes, coagulation and complement cascades, and multiple pathways related to tissue repair likely played critical roles in the process. Our results offered a quantitative snapshot of changes in the host transcriptome induced by H. contortus infection and provided novel insights into molecular mechanisms of host resistance.