Transcriptome characterization by RNA-seq unravels the mechanisms of butyrate-induced epigenomic regulation in bovine cells

Gut microbiota mediates ambient PM2.5 exposure-induced abnormal glucose metabolism via short-chain fatty acids

PM2.5 exposure has been found to cause gut dysbiosis and impair glucose homeostasis in human and animals, yet their underlying biological connection remain unclear. In the present study, we aim to investigate the biological significance of gut microbiota in PM2.5-induced glucose metabolic abnormalities. Our results showed that microbiota depletion by antibiotics treatment significantly alleviated PM2.5-induced glucose intolerance and insulin resistance, as indicated by the intraperitoneal glucose tolerance test, glucose-induced insulin secretion, insulin tolerance test, insulin-induced phosphorylation levels of Akt and GSK-3β in insulin sensitive tissues. In addition, faecal microbiota transplantation (FMT) from PM2.5-exposed donor mice successfully remodeled the glucose metabolism abnormalities in recipient mice, while the transplantation of autoclaved faecal materials did not. Faecal microbiota analysis demonstrated that the composition and alpha diversity of the gut bacterial community were altered by PM2.5 exposure and in FMT recipient mice. Furthermore, short-chain fatty acids levels analysis showed that the circulating acetate was significantly decreased in PM2.5-exposed donor and FMT recipient mice, and supplementation of sodium acetate for 3 months successfully improved the glucose metabolism abnormalities induced by PM2.5 exposure. These results indicate that manipulating gut microbiota or its metabolites could be a potential strategy for preventing the adverse health effects of ambient PM2.5.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08758-4.

Microbial Metabolites Orchestrate a Distinct Multi-Tiered Regulatory Network in the Intestinal Epithelium That Directs P-Glycoprotein Expression

P-glycoprotein (P-gp) is a key component of the intestinal epithelium playing a pivotal role in removal of toxins and efflux of endocannabinoids to prevent excessive inflammation and sustain homeostasis. Recent studies revealed butyrate and secondary bile acids, produced by the intestinal microbiome, potentiate the induction of functional P-gp expression. We now aim to determine the molecular mechanism by which this functional microbiome output regulates P-gp. RNA sequencing of intestinal epithelial cells responding to butyrate and secondary bile acids in combination discovered a unique transcriptional program involving multiple pathways that converge on P-gp induction. Using shRNA knockdown and CRISPR/Cas9 knockout cell lines, as well as mouse models, we confirmed the RNA sequencing findings and discovered a role for intestinal HNF4α in P-gp regulation. These findings shed light on a sophisticated signaling network directed by intestinal microbial metabolites that orchestrate P-gp expression and highlight unappreciated connections between multiple pathways linked to colonic health. IMPORTANCE Preventing aberrant inflammation is essential to maintaining homeostasis in the mammalian intestine. Although P-glycoprotein (P-gp) expression in the intestine is critical for protecting the intestinal epithelium from toxins and damage due to neutrophil infiltration, its regulation in the intestine is poorly understood. Findings presented in our current study have now uncovered a sophisticated and heretofore unappreciated intracellular signaling network or "reactome" directed by intestinal microbial metabolites that orchestrate regulation of P-gp. Not only do we confirm the role of histone deacetylases (HDAC) inhibition and nuclear receptor activation in P-gp induction by butyrate and bile acids, but we also discovered new signaling pathways and transcription factors that are uniquely activated in response to the combination of microbial metabolites. Such findings shed new light into a multi-tiered network that maintains P-gp expression in the intestine in the context of the fluctuating commensal microbiome, to sustain a homeostatic tone in the absence of infection or insult.

Camptothecin Encapsulated in β-Cyclodextrin-EDTA-Fe3O4 Nanoparticles Induce Metabolic Reprogramming Repair in HT29 Cancer Cells through Epigenetic Modulation: A Bioinformatics Approach

Cancer progresses through a distinctive reprogramming of metabolic pathways directed by genetic and epigenetic modifications. The hardwired changes induced by genetic mutations are resilient, while epigenetic modifications are softwired and more vulnerable to therapeutic intervention. Colon cancer is no different. This gives us the need to explore the mechanism as an attractive therapeutic target to combat colon cancer cells. We have previously established the enhanced therapeutic efficacy of a newly formulated camptothecin encapsulated in β-cyclodextrin-EDTA-Fe₃O₄ nanoparticles (CPT-CEF) in colon cancer cells. We furthered this study by carrying out RNA sequencing (RNA-seq) to underscore specific regulatory signatures in the CPT-CEF treated versus untreated HT29 cells. In the study, we identified 95 upregulated and 146 downregulated genes spanning cellular components and molecular and metabolic functions. We carried out extensive bioinformatics analysis to harness genes potentially involved in epigenetic modulation as either the cause or effect of metabolic rewiring exerted by CPT-CEF. Significant downregulation of 13 genes involved in the epigenetic modulation and 40 genes from core metabolism was identified. Three genes, namely, DNMT-1, POLE3, and PKM-2, were identified as the regulatory overlap between epigenetic drivers and metabolic reprogramming in HT29 cells. Based on our results, we propose a possible mechanism that intercepts the two functional axes, namely epigenetic control, and metabolic modulation via CPT-CEF in colon cancer cells, which could skew cancer-induced metabolic deregulation towards metabolic repair. Thus, the study provides avenues for further validation of transcriptomic changes affected by these deregulated genes at epigenetic level, and ultimately may be harnessed as targets for regenerating normal metabolism in colon cancer with better treatment potential, thereby providing new avenues for colon cancer therapy.

Transcriptional Reprogramming in Rumen Epithelium during the Developmental Transition of Pre-Ruminant to the Ruminant in Cattle

Simple Summary

The rumen is the critical organ mediating nutrient uptake and use in cattle. Health development is essential to ensure animal feed efficiency. In this report, we present an analysis of gene expression dynamic in rumen epithelium during the transition from pre-ruminant to ruminant in cattle fed with hay or concentrated diets at weaning. The global shifts in gene expression and correlated transcription factors activities indicate transcriptional reprogramming during weaning. Transcriptional reprogramming in rumen epithelial tissue reflects critical nutrient-gene interactions occurring during the developmental progression. The results unveiled that nutrient-gene interactions compel transcriptional reprogramming. Our findings also suggest that this transcriptional reprogramming is the molecular basis of the transitional development of pre-ruminant to the ruminant in cattle.

Abstract

We present an analysis of transcriptomic dynamics in rumen epithelium of 18 Holstein calves during the transition from pre-rumination to rumination in cattle-fed hay or concentrated diets at weaning. Three calves each were euthanized at 14 and 42 d of age to exemplify preweaning, and six calves each were provided diets of either milk replacer and grass hay or calf starter to introduce weaning. The two distinct phases of rumen development and function in cattle are tightly regulated by a series of signaling events and clusters of effectors on critical pathways. The dietary shift from liquid to solid feeds prompted the shifting of gene activity. The number of differentially expressed genes increased significantly after weaning. Bioinformatic analysis revealed gene activity shifts underline the functional transitions in the ruminal epithelium and signify the transcriptomic reprogramming. Gene ontogeny (GO) term enrichment shows extensively activated biological functions of differentially expressed genes in the ruminal epithelium after weaning were predominant metabolic functions. The transcriptomic reprogramming signifies a correlation between gene activity and changes in metabolism and energy production in the rumen epithelium, which occur at weaning when transitioning from glucose use to VFA use by epithelium during the weaning.

Transcriptomic analyses show that 24-epibrassinolide (EBR) promotes cold tolerance in cotton seedlings

In plants, brassinosteroids (BRs) are a class of steroidal hormones that are involved in numerous physiological responses. However, the function of BRs in cold tolerance in cotton has not been explored. In this study, cotton seedlings were treated with five concentrations (0, 0.05, 0.1, 0.2, 0.5 and 1.0 mg/L) of 24-Epibrassinolide (EBR) at 4°C. We measured the electrolyte leakage, malondialdehyde (MDA) content, proline content, and net photosynthesis rate (Pn) of the seedlings, which showed that EBR treatment increased cold tolerance in cotton in a dose-dependent manner, and that 0.2 mg/L is an optimum concentration for enhancing cold tolerance. The function of EBR in cotton cotyledons was investigated in the control 0 mg/L (Cold+water) and 0.2 mg/L (Cold+EBR) treatments using RNA-Seq. A total of 4,001 differentially expressed genes (DEGs), including 2,591 up-regulated genes and 1,409 down-regulated genes were identified. Gene Ontology (GO) and biochemical pathway enrichment analyses showed that EBR is involved in the genetic information process, secondary metabolism, and also inhibits abscisic acid (ABA) and ethylene (ETH) signal transduction. In this study, physiological experiments showed that EBR can increase cold tolerance in cotton seedlings, and the comprehensive RNA-seq data shed light on the mechanisms through which EBR increases cold tolerance in cotton seedlings.

Establishment and transcriptomic analyses of a cattle rumen epithelial primary cells (REPC) culture by bulk and single-cell RNA sequencing to elucidate interactions of butyrate and rumen development

As a critical and high-value tool to study the development of rumen, we established a stable rumen epithelial primary cell (REPC) culture from a two-week-old Holstein bull calf rumen epithelial tissue. The transcriptomic profiling of the REPC and the direct effects of butyrate on gene expression were assessed. Correlated gene networks elucidated the putative roles and mechanisms of butyrate action in rumen epithelial development. The top networks perturbed by butyrate were associated with epithelial tissue development. Additionally, two critical upstream regulators, E2F1 and TGFB1, were identified to play critical roles in the differentiation, development, and growth of epithelial cells. Significant expression changes of upstream regulators and transcription factors provided further evidence in support that butyrate plays a specific and central role in regulating genomic and epigenomic activities influencing rumen development. This work is the essential component to obtain a complete global landscape of regulatory elements in cattle and to explore the dynamics of chromatin states in rumen epithelial cells induced by butyrate at early developmental stages.

Monobutyrin and Monovalerin Affect Brain Short-Chain Fatty Acid Profiles and Tight-Junction Protein Expression in ApoE-Knockout Rats Fed High-Fat Diets

Monobutyrin (MB) and monovalerin (MV), esters of short-chain fatty acids (SCFAs), have previously been shown to reduce liver cholesterol and inflammation in conventional rats fed high-fat diets. This study explored the potential effects of MB and MV in hypercholesterolemic apolipoprotein E-knockout (ApoE-/-) rats. ApoE-/- rats were fed three high-fat (HF) diets, pure or supplemented with MB or MV (1%), for 5 weeks. One group of conventional rats (C) was also fed the pure high-fat diet and another group of ApoE-/- rats a low-fat (LF) diet. Blood and liver lipids, urinary lactulose/mannitol, SCFAs (blood and brain), tight junction proteins (small intestine and brain), and inflammation-related markers (blood, brain, and liver) were analyzed. MV supplementation elevated serum high-density lipoprotein (HDL) cholesterol and valeric acid concentration (p < 0.05), while the amounts of isovaleric acid in the brain were reduced (p < 0.05). MB increased butyric acid amounts in the brain, while the plasma concentration of interleukin 10 (IL-10) was lowered (p < 0.05). Both MV and MB upregulated the expression of occludin and zonula occludens-1 (ZO-1) in the brain (p < 0.05). Supplementation of MB or MV affected HDL cholesterol, the expression of tight junction proteins, and SCFA profiles. MB and MV may therefore be promising supplements to attenuate lipid metabolic disorders caused by high-fat intake and genetic deficiency.

Immunomodulating Activity and Therapeutic Effects of Short Chain Fatty Acids and Tryptophan Post-biotics in Inflammatory Bowel Disease

Crohn's disease (CD) and Ulcerative colitis (UC) are grouped as Inflammatory Bowel Diseases (IBD). The IBD is associated to a multifaceted interplay between immunologic, microbial, genetic, and environmental factors. Nowadays, the gut microbiota (GM) dysbiosis has been indicated as a cause in the IBD development, affecting the impaired cross-talk between GM and immune cells. Moreover, recent studies have uncovered a crucial role for bacterial post-biotics (metabolites) in the orchestration of the host immune response, as they could be messengers between the GM and the immune system. In addition, transgenic mouse models showed that SCFAs (Short Chain Fatty Acids) and Tryptophan (Trp) post-biotics play important immunomodulatory effects, regulating both innate and adaptive immune cell generation, their function and trafficking. Here, we present an overview on the main microbial post-biotics and their effects on the gut mucosa with specific emphasis on their relevance for IBD. Finally, we discuss the therapeutic potential of SCFA and Trp post-biotics on IBD through approaches based on the “immunonutrition,” defined as a modulation of the immune system provided by specific interventions that modify dietary nutrients.

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.

Bacterial endotoxin decreased histone H3 acetylation of bovine mammary epithelial cells and the adverse effect was suppressed by sodium butyrate

BACKGROUND

In practical production, dairy cows are frequently exposed to bacterial endotoxin (lipopolysaccharide, LPS) when they are subjected to high-concentrate diets, poor hygienic environments, as well as mastitis and metritis. Histone acetylation is an important epigenetic control of DNA transcription and a higher histone acetylation is associated with facilitated transcription. LPS might reduce histone acetylation in the mammary epithelial cells, resulting in lower transcription and mRNA expression of lactation-related genes. This study was conducted to investigate the effect of LPS on histone acetylation in bovine mammary epithelial cells and the efficacy of sodium butyrate (SB) in suppressing the endotoxin-induced adverse effect. Firstly, the bovine mammary epithelial cell line MAC-T cells were treated for 48 h with LPS at different doses of 0, 1, 10, 100, and 1000 endotoxin units (EU)/mL (1 EU = 0.1 ng), and the acetylation levels of histones H3 and H4 as well as the histone deacetylase (HDAC) activity were measured. Secondly, the MAC-T cells were treated for 48 h as follows: control, LPS (100 EU/mL), and LPS (100 EU/mL) plus SB (10 mmol/L), and the acetylation levels of histones H3 and H4 as well as milk gene mRNA expressions were determined.

RESULTS

The results showed that HDAC activity increased linearly with increasing LPS doses (P < 0.01). The histone H3 acetylation levels were significantly reduced by LPS, while the histone H4 acetylation levels were not affected by LPS (P > 0.05). Sodium butyrate, an inhibitor of HDAC, effectively suppressed the endotoxin-induced decline of histone H3 acetylation (P < 0.05). As a result, SB significantly enhanced the mRNA expression of lactation-related genes (P < 0.05).

CONCLUSIONS

The results suggest one of the adverse effects of LPS on the lactation of bovine mammary gland epithelial cells was due to decreasing histone H3 acetylation through increasing HDAC activity, whereas the endotoxin-induced adverse effects were effectively suppressed by SB.

Sodium Butyrate as a Histone Deacetylase Inhibitor Affects Toll-Like Receptor 4 Expression in Colorectal Cancer Cell Lines

We assessed the effect of sodium butyrate (SB) as a histone deacetylase inhibitor (HDACi) on Toll-like receptor 4 (TLR4) gene expression levels, in low TLR4 expressing (HCT116) and high TLR4 expressing (SW480) colorectal cancer cells. The cytotoxic effect of SB was assessed by culturing SW480 and HCT116 cell lines using a broad spectrum of times and concentrations of SB. The MTT assay was done to check the cytotoxic properties of different SB concentrations. Gene expression levels of TLR4 was then evaluated for non-cytotoxic SB concentrations. Morphological analysis and MTT assay confirmed that SB concentrations equal to or less than 5mM were not cytotoxic for both cell lines. At 5mM concentration of SB in SW480 cell line and 1mM concentration of SB in HCT116 cell line, TLR4 gene expression level significantly increased from 24 to 48 hrs and decreased significantly from 48 to 72 hrs with an "early increased and late decreased pattern". At 1mM concentration of SB in SW480 cell line and 5mM concentration of SB in HCT116 cell line, TLR4 expression had a "gradually increased pattern". This study focuses on the dose-time-effect of SB in the pathogenesis of colorectal cancer. SB alters the expression level of TLR4 in colorectal cancer cells. This effect may depend on the cell type, treatment duration and SB concentration. The alterations in TLR4 expression may be due to the direct effect of SB on TLR4 and/or the expression changes of in other genes which may indirectly affect the TLR4 expression. Abbreviations: TLR4: Toll-like receptor 4; HDACi: histone deacetylase inhibitor; SB: sodium Butyrate; CRC: colorectal cancer; SCFA: short-chain fatty acid; hrs: hours.

Effect of consuming endophyte-infected fescue seed on transcript abundance in the mammary gland of lactating and dry cows, as assessed by RNA sequencing

Ergot alkaloids in endophyte-infected grasses inhibit prolactin secretion and reduce milk production in lactating cows. However, we previously showed that prepartum consumption of infected seed throughout the dry period did not inhibit subsequent milk production and prior exposure to bromocriptine (ergot peptide) actually increased production in the next lactation. To identify changes in the transcriptome and molecular pathways mediating the mammary gland's response to ergot alkaloids in the diet, RNA sequencing (RNA-seq) was performed on mammary tissues obtained from 22 multiparous Holstein cows exposed to 1 of 3 treatments. Starting at 90 ± 4 d prepartum, cows were fed endophyte-free fescue seed (control; CON), endophyte-free fescue seed plus 3×/wk subcutaneous injections of bromocriptine (BROMO; 0.1 mg/kg of BW), or endophyte-infected fescue seed (INF) as 10% of the diet. Cows were dried off 60 ± 2 d prepartum. Mammary biopsies from 4 (BROMO, INF) or 5 (CON) cows/treatment at each of the 3 phases were obtained: 7 d before dry off during the initial lactation (L1), mid-dry period (D), and 10 d postpartum (L2). Although tissue from the same cow was preferentially used at 3 phases (L1, D, L2), tissue from additional cows were used to as necessary to provide RNA of sufficient quality. Individual samples were used to generate individual RNA-seq libraries. Normalized reads of the RNA-seq data were organized into technical and biological replicates before processing with the RSEM software package. Each lactation phase was processed separately and genes that differed between any of 3 treatments were identified. A large proportion of genes differentially expressed in at least 1 treatment (n = 866) were found to be similarly expressed in BROMO and INF treatments, but differentially expressed from CON (n = 575, total for 3 phases). Of genes differentially expressed compared with CON, 104 genes were common to the L1 and L2 phases. Consistent with the production findings, networks most affected by treatments in L1 and L2 included lipid metabolism, small molecule biochemistry, and molecular transport, whereas networks related more to developmental and cellular functions and maintenance were evident during D phase. Similar patterns of expression in BROMO and INF during late and early lactation suggest involvement of similar cell signaling pathways or mechanisms of action for BROMO and INF and the importance of prolactin messaging pathways.

Transcriptome analysis revealed anti-obesity effects of the Sodium Alginate in high-fat diet -induced obese mice

Human obesity and overweight, caused by accumulated of fat, is the most commonly phenomenon from all over the world, especially in Western countries and Chinese mainland during the past three decades. Sodium Alginate, a polysaccharide extracted from brown seaweeds, has been proved its strong ability on body weight loss and anti-inflammatory response. However, no studies have been explored the effects of Sodium Alginate on colonic transcriptome, especially in obese individuals. Therefore, the current study was designed to detect whether Sodium Alginate could remit obesity and ease chronic metabolism disease through strengthening the bio-functionality of the lower intestine, particularly in colon. The data showed after Sodium Alginate gavaged for four weeks, the body weight, fat accumulation, triglyceride and total cholesterol were ameliorated in high fat diet induced obese mice. Sodium Alginate also improved the blood glucose level and lipopolysaccharides in serum. Furthermore, data from RNA sequence indicated that there were significantly changes in several genes, which involved in lipid metabolism and carbohydrate metabolism. In conclusion, these results suggested that Sodium Alginate could effectively suppress obesity and obesity related metabolic syndromes, due to the colonic transcriptome changes.

A Cross-Talk Between Microbiota-Derived Short-Chain Fatty Acids and the Host Mucosal Immune System Regulates Intestinal Homeostasis and Inflammatory Bowel Disease

Gut microbiota has a fundamental role in the energy homeostasis of the host and is essential for proper "education" of the immune system. Intestinal microbial communities are able to ferment dietary fiber releasing short-chain fatty acids (SCFAs). The SCFAs, particularly butyrate (BT), regulate innate and adaptive immune cell generation, trafficing, and function. For example, BT has an anti-inflammatory effect by inhibiting the recruitment and proinflammatory activity of neutrophils, macrophages, dendritic cells, and effector T cells and by increasing the number and activity of regulatory T cells. Gut microbial dysbiosis, ie, a microbial community imbalance, has been suggested to play a role in the development of inflammatory bowel disease (IBD). The relationship between dysbiosis and IBD has been difficult to prove, especially in humans, and is probably complex and dynamic, rather than one of a simple cause and effect relationship. However, IBD patients have dysbiosis with reduced numbers of SCFAs-producing bacteria and reduced BT concentration that is linked to a marked increase in the number of proinflammatory immune cells in the gut mucosa of these patients. Thus, microbial dysbiosis and reduced BT concentration may be a factor in the emergence and severity of IBD. Understanding the relationship between microbial dysbiosis and reduced BT concentration to IBD may lead to novel therapeutic interventions.

Integrative analysis of DNA methylation and gene expression in butyrate-treated CHO cells

The cellular mechanisms responsible for the versatile properties of CHO cells as the major production cell line for biopharmaceutical molecules are not entirely understood yet, although several 'omics' data facilitate the understanding of CHO cells and their reactions to environmental conditions. However, genome-wide studies of epigenetic processes such as DNA methylation are still limited. To prove the applicability and usefulness of integrating DNA methylation and gene expression data in a biotechnological context, we exemplarily analyzed the time course of cellular reactions upon butyrate addition in antibody-producing CHO cells by whole-genome bisulfite sequencing and CHO-specific cDNA microarrays. Gene expression and DNA methylation analyses showed that pathways known to be affected by butyrate, including cell cycle and apoptosis, as well as pathways potentially involved in butyrate-induced hyperproductivity such as central energy metabolism and protein biosynthesis were affected. Differentially methylated regions were furthermore found to contain binding-site motifs of specific transcription factors and were hypothesized to represent regulatory regions closely connected to the cellular response to butyrate. Generally, our experiment underlines the benefit of integrating DNA methylation and gene expression data, as it provided potential novel candidate genes for rational cell line development and allowed for new insights into the butyrate effect on CHO cells.

Low Concentration of Sodium Butyrate from Ultrabraid+NaBu suture, Promotes Angiogenesis and Tissue Remodelling in Tendon-bones Injury

Sodium butyrate (NaBu), a form of short-chain fatty acid (SCFA), acts classically as a potent anti-angiogenic agent in tumour angiogenesis models, some authors demonstrated that low concentrations of NaBu may contribute to healing of tendon-bone injury in part at least through promotion of tissue remodelling. Here, we investigated the effects of low-range concentrations of NaBu using in vitro and in vivo assays using angiogenesis as the primary outcome measure and the mechanisms through which it acts. We demonstrated that NaBu, alone or perfused from the UltraBraid+NaBu suture was pro-angiogenic at very low-range doses promoting migration, tube formation and cell invasion in bovine aortic endothelial cells (BAECs). Furthermore, cell exposure to low NaBu concentrations increased expression of proteins involved in angiogenic cell signalling, including p-PKCβ1, p-FAK, p-ERK1/2, p-NFκβ, p-PLCγ1 and p-VEGFR2. In addition, inhibitors of both VEGFR2 and PKCβ1 blocked the angiogenic response. In in vivo assays, low concentrations of NaBu induced neovascularization in sponge implants in mice, evidenced by increased numbers of vessels and haemoglobin content in these implants. The findings in this study indicate that low concentrations of NaBu could be an important compound to stimulate angiogenesis at a site where vasculature is deficient and healing is compromised.

TRIENNIAL LACTATION SYMPOSIUM: Nutrigenomics in livestock: Systems biology meets nutrition

The advent of high-throughput technologies to study an animal's genome, proteome, and metabolome (i.e., "omics" tools) constituted a setback to the use of reductionism in livestock research. More recent development of "next-generation sequencing" tools was instrumental in allowing in-depth studies of the microbiome in the rumen and other sections of the gastrointestinal tract. Omics, along with bioinformatics, constitutes the foundation of modern systems biology, a field of study widely used in model organisms (e.g., rodents, yeast, humans) to enhance understanding of the complex biological interactions occurring within cells and tissues at the gene, protein, and metabolite level. Application of systems biology concepts is ideal for the study of interactions between nutrition and physiological state with tissue and cell metabolism and function during key life stages of livestock species, including the transition from pregnancy to lactation, in utero development, or postnatal growth. Modern bioinformatic tools capable of discerning functional outcomes and biologically meaningful networks complement the ever-increasing ability to generate large molecular, microbial, and metabolite data sets. Simultaneous visualization of the complex intertissue adaptations to physiological state and nutrition can now be discerned. Studies to understand the linkages between the microbiome and the absorptive epithelium using the integrative approach are emerging. We present examples of new knowledge generated through the application of functional analyses of transcriptomic, proteomic, and metabolomic data sets encompassing nutritional management of dairy cows, pigs, and poultry. Published work to date underscores that the integrative approach across and within tissues may prove useful for fine-tuning nutritional management of livestock. An important goal during this process is to uncover key molecular players involved in the organismal adaptations to nutrition.

Transcriptomic Sequencing Reveals a Set of Unique Genes Activated by Butyrate-Induced Histone Modification

Butyrate is a nutritional element with strong epigenetic regulatory activity as a histone deacetylase inhibitor. Based on the analysis of differentially expressed genes in the bovine epithelial cells using RNA sequencing technology, a set of unique genes that are activated only after butyrate treatment were revealed. A complementary bioinformatics analysis of the functional category, pathway, and integrated network, using Ingenuity Pathways Analysis, indicated that these genes activated by butyrate treatment are related to major cellular functions, including cell morphological changes, cell cycle arrest, and apoptosis. Our results offered insight into the butyrate-induced transcriptomic changes and will accelerate our discerning of the molecular fundamentals of epigenomic regulation.

Epigenetics and transgenerational inheritance in domesticated farm animals

Epigenetics provides a molecular mechanism of inheritance that is not solely dependent on DNA sequence and that can account for non-Mendelian inheritance patterns. Epigenetic changes underlie many normal developmental processes, and can lead to disease development as well. While epigenetic effects have been studied in well-characterized rodent models, less research has been done using agriculturally important domestic animal species. This review will present the results of current epigenetic research using farm animal models (cattle, pigs, sheep and chickens). Much of the work has focused on the epigenetic effects that environmental exposures to toxicants, nutrients and infectious agents has on either the exposed animals themselves or on their direct offspring. Only one porcine study examined epigenetic transgenerational effects; namely the effect diet micronutrients fed to male pigs has on liver DNA methylation and muscle mass in grand-offspring (F2 generation). Healthy viable offspring are very important in the farm and husbandry industry and epigenetic differences can be associated with production traits. Therefore further epigenetic research into domestic animal health and how exposure to toxicants or nutritional changes affects future generations is imperative.

A high-resolution whole-genome map of the distinctive epigenomic landscape induced by butyrate in bovine cells

This report presents a study utilizing next-generation sequencing technology, combined with chromatin immunoprecipitation (ChIP-seq) technology to analyze histone modification induced by butyrate and to construct a high-definition map of the epigenomic landscape with normal histone H3 and H4 and their variants in bovine cells at the whole-genome scale. A total of 10 variants of histone H3 and H4 modifications were mapped at the whole-genome scale (acetyl-H3K18-ChIP-seq, trimethy-H3K9, histone H4 ChIP-seq, acetyl-H4K5 ChIP-seq, acetyl-H4K12 ChIP-seq, acetyl-H4K16 ChIP-seq, histone H3 ChIP-seq, acetyl H3H9 ChIP-seq, acetyl H3K27 ChIP-seq and tetra-acetyl H4 ChIP-seq). Integrated experiential data and an analysis of histone and histone modification at a single base resolution across the entire genome are presented. We analyzed the enriched binding regions in the proximal promoter (within 5 kb upstream or at the 5'-untranslated region from the transcriptional start site (TSS)), and the exon, intron and intergenic regions (defined by regions 25 kb upstream and 10 kb downstream from the TSS). A de novo search for the binding motif of the 10 ChIP-seq datasets discovered numerous motifs from each of the ChIP-seq datasets. These consensus sequences indicated that histone modification at different locations changes the histone H3 and H4 binding preferences. Nevertheless, a high degree of conservation in histone binding also was presented in these motifs. This first extensive epigenomic landscape mapping in bovine cells offers a new framework and a great resource for testing the role of epigenomes in cell function and transcriptomic regulation.

Comparative epigenetics: relevance to the regulation of production and health traits in cattle

With the development of genomic, transcriptomic and bioinformatic tools, recent advances in molecular technologies have significantly impacted bovine bioscience research and are revolutionising animal selection and breeding. Integration of epigenetic information represents yet another challenging molecular frontier. Epigenetics is the study of biochemical modifications to DNA and to histones, the proteins that provide stability to DNA. These epigenetic changes are induced by environmental stimuli; they alter gene expression and are potentially heritable. Epigenetics research holds the key to understanding how environmental factors contribute to phenotypic variation in traits of economic importance in cattle including development, nutrition, behaviour and health. In this review, we discuss the potential applications of epigenetics in bovine research, using breakthroughs in human and murine research to signpost the way.

Bioinformatic Dissecting of TP53 Regulation Pathway Underlying Butyrate-induced Histone Modification in Epigenetic Regulation

Butyrate affects cell proliferation, differentiation, and motility. Butyrate inhibits histone deacetylase (HDAC) activities and induces cell-cycle arrest and apoptosis. TP53 is one of the most active upstream regulators discovered by ingenuity pathways analysis (IPA) in our RNA-sequencing data set. TP53 signaling pathway plays key role in many cellular processes. TP53 pathway and their involvement in cellular functions modified by butyrate treatment were scrutinized in this report by data mining the RNA-sequencing data using IPA (Ingenuity System^®). The TP53 mechanistic pathway targets more than 600 genes. Downstream analysis predicted the activation of the TP53 pathway after butyrate treatment. The data mining also revealed that nine transcription factors are downstream regulators in TP53 signaling pathways. The analysis results also indicated that butyrate not only inhibits the HDAC activities, but also regulates genes encoding the HDAC enzymes through modification of histones and epigenomic landscape.

The microbiome and colorectal neoplasia: environmental modifiers of dysbiosis

The etiology of colon cancer is complex, yet it is undoubtedly impacted by intestinal microbiota. Whether the contribution to colon carcinogenesis is generated through the presence of an overall dysbiosis or by specific pathogens is still a matter for debate. However, it is apparent that interactions between microbiota and the host are mediated by a variety of processes, including signaling cascades, the immune system, host metabolism, and regulation of gene transcription. To fully appreciate the role of microbiota in colon carcinogenesis, it will be necessary to expand efforts to define populations in niche environments, such as colonic crypts, explore cross talk between the host and the microbiota, and more completely define the metabolomic profile of the microbiota. These efforts must be pursued with appreciation that dietary substrates and other environmental modifiers mediate changes in the microbiota, as well as their metabolism and functional characteristics.

Identification of novel transcripts and noncoding RNAs in bovine skin by deep next generation sequencing

BACKGROUND

Deep RNA sequencing (RNAseq) has opened a new horizon for understanding global gene expression. The functional annotation of non-model mammalian genomes including bovines is still poor compared to that of human and mouse. This particularly applies to tissues without direct significance for milk and meat production, like skin, in spite of its multifunctional relevance for the individual. Thus, applying an RNAseq approach, we performed a whole transcriptome analysis of pigmented and nonpigmented bovine skin to describe the comprehensive transcript catalogue of this tissue.

RESULTS

A total of 39,577 unique primary skin transcripts were mapped to the bovine reference genome assembly. The majority of the transcripts were mapped to known transcriptional units (65%). In addition to the reannotation of known genes, a substantial number (10,884) of unknown transcripts (UTs) were discovered, which had not previously been annotated. The classification of UTs was based on the prediction of their coding potential and comparative sequence analysis, subsequently followed by meticulous manual curation. The classification analysis and experimental validation of selected UTs confirmed that RNAseq data can be used to amend the annotation of known genes by providing evidence for additional exons, untranslated regions or splice variants, by approving genes predicted in silico and by identifying novel bovine loci. A large group of UTs (4,848) was predicted to potentially represent long noncoding RNA (lncRNA). Predominantly, potential lncRNAs mapped in intergenic chromosome regions (4,365) and therefore, were classified as potential intergenic lncRNA. Our analysis revealed that only about 6% of all UTs displayed interspecies conservation and discovered a variety of unknown transcripts without interspecies homology but specific expression in bovine skin.

CONCLUSIONS

The results of our study demonstrate a complex transcript pattern for bovine skin and suggest a possible functional relevance of novel transcripts, including lncRNA, in the modulation of pigmentation processes. The results also indicate that the comprehensive identification and annotation of unknown transcripts from whole transcriptome analysis using RNAseq data remains a tremendous future challenge.

Establishment of a CpG island microarray for analyses of genome-wide DNA methylation in Chinese hamster ovary cells

Optimizing productivity and growth rates of recombinant Chinese hamster ovary (CHO) cells requires insight into the regulation of cellular processes. In this regard, the elucidation of the epigenetic process of DNA methylation, known to influence transcription by a differential occurrence in CpG islands in promoter regions, is increasingly gaining importance. However, DNA methylation has not yet been investigated on a genomic scale in CHO cells and suitable tools have not existed until now. Based on the genomic and transcriptomic CHO data currently available, we developed a customized oligonucleotide microarray covering 19598 CpG islands (89 % of total bioinformatically identified CpG islands) in the CHO genome. We applied our CHO-specific CpG island microarray to investigate the effect of butyrate treatment on differential DNA methylation in CHO cultures in a time-dependent approach. Supplementation of butyrate is known to enhance cell specific productivities in CHO cells and leads to alterations of epigenetic silencing events. Gene ontology clusters regarding, e.g., chromatin modification or DNA repair, were significantly overrepresented 24 h after butyrate addition. Functional classifications furthermore indicated that several major signaling systems such as the Wnt/β-catenin pathway were affected by butyrate treatment. Our novel CHO-specific CpG island microarray will provide valuable information in future studies of cellular processes associated with productivity and product characteristics.

The anti-tumor effect of HDAC inhibition in a human pancreas cancer model is significantly improved by the simultaneous inhibition of cyclooxygenase 2

Pancreatic ductal adenocarcinoma is the fourth leading cause of cancer death worldwide, with no satisfactory treatment to date. In this study, we tested whether the combined inhibition of cyclooxygenase-2 (COX-2) and class I histone deacetylase (HDAC) may results in a better control of pancreatic ductal adenocarcinoma. The impact of the concomitant HDAC and COX-2 inhibition on cell growth, apoptosis and cell cycle was assessed first in vitro on human pancreas BxPC-3, PANC-1 or CFPAC-1 cells treated with chemical inhibitors (SAHA, MS-275 and celecoxib) or HDAC1/2/3/7 siRNA. To test the potential antitumoral activity of this combination in vivo, we have developed and characterized, a refined chick chorioallantoic membrane tumor model that histologically and proteomically mimics human pancreatic ductal adenocarcinoma. The combination of HDAC1/3 and COX-2 inhibition significantly impaired proliferation of BxPC-3 cells in vitro and stalled entirely the BxPC-3 cells tumor growth onto the chorioallantoic membrane in vivo. The combination was more effective than either drug used alone. Consistently, we showed that both HDAC1 and HDAC3 inhibition induced the expression of COX-2 via the NF-kB pathway. Our data demonstrate, for the first time in a Pancreatic Ductal Adenocarcinoma (PDAC) model, a significant action of HDAC and COX-2 inhibitors on cancer cell growth, which sets the basis for the development of potentially effective new combinatory therapies for pancreatic ductal adenocarcinoma patients.

Dietary inhibitors of histone deacetylases in intestinal immunity and homeostasis

Intestinal epithelial cells (IECs) are integral players in homeostasis of immunity and host defense in the gut and are under influence of the intestinal microbiome. Microbial metabolites and dietary components, including short chain fatty acids (acetate, propionate, and butyrate, SCFAs), have an impact on the physiology of IECs at multiple levels, including the inhibition of deacetylases affecting chromatin remodeling and global changes in transcriptional activity. The number and diversity of butyrate-producing bacteria is subject to factors related to age, disease, and to diet. At physiological levels, SCFAs are inhibitors of histone deacetylases (HDACs) which may explain the transcriptional effects of SCFAs on epithelial cells, although many effects of SCFAs on colonic mucosa can be ascribed to mechanisms beyond HDAC inhibition. Interference with this type of post-translational modification has great potential in cancer and different inflammatory diseases, because HDAC inhibition has anti-proliferative and anti-inflammatory effects in vitro, and in in vivo models of intestinal inflammation. Hence, the influence of dietary modulators on HDAC activity in epithelia is likely to be an important determinant of its responses to inflammatory and microbial challenges.

RNA-Seq analysis of Cocos nucifera: transcriptome sequencing and de novo assembly for subsequent functional genomics approaches

BACKGROUND

Cocos nucifera (coconut), a member of the Arecaceae family, is an economically important woody palm grown in tropical regions. Despite its agronomic importance, previous germplasm assessment studies have relied solely on morphological and agronomical traits. Molecular biology techniques have been scarcely used in assessment of genetic resources and for improvement of important agronomic and quality traits in Cocos nucifera, mostly due to the absence of available sequence information.

METHODOLOGY/PRINCIPAL FINDINGS

To provide basic information for molecular breeding and further molecular biological analysis in Cocos nucifera, we applied RNA-seq technology and de novo assembly to gain a global overview of the Cocos nucifera transcriptome from mixed tissue samples. Using Illumina sequencing, we obtained 54.9 million short reads and conducted de novo assembly to obtain 57,304 unigenes with an average length of 752 base pairs. Sequence comparison between assembled unigenes and released cDNA sequences of Cocos nucifera and Elaeis guineensis indicated that the assembled sequences were of high quality. Approximately 99.9% of unigenes were novel compared to the released coconut EST sequences. Using BLASTX, 68.2% of unigenes were successfully annotated based on the Genbank non-redundant (Nr) protein database. The annotated unigenes were then further classified using the Gene Ontology (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases.

CONCLUSIONS/SIGNIFICANCE

Our study provides a large quantity of novel genetic information for Cocos nucifera. This information will act as a valuable resource for further molecular genetic studies and breeding in coconut, as well as for isolation and characterization of functional genes involved in different biochemical pathways in this important tropical crop species.

Modulation of antimicrobial host defense peptide gene expression by free fatty acids

Routine use of antibiotics at subtherapeutic levels in animal feed drives the emergence of antimicrobial resistance. Development of antibiotic-alternative approaches to disease control and prevention for food animals is imperatively needed. Previously, we showed that butyrate, a major species of short-chain fatty acids (SCFAs) fermented from undigested fiber by intestinal microflora, is a potent inducer of endogenous antimicrobial host defense peptide (HDP) genes in the chicken (PLoS One 2011, 6: e27225). In the present study, we further revealed that, in chicken HD11 macrophages and primary monocytes, induction of HDPs is largely in an inverse correlation with the aliphatic hydrocarbon chain length of free fatty acids, with SCFAs being the most potent, medium-chain fatty acids moderate and long-chain fatty acids marginal. Additionally, three SCFAs, namely acetate, propionate, and butyrate, exerted a strong synergy in augmenting HDP gene expression in chicken cells. Consistently, supplementation of chickens with a combination of three SCFAs in water resulted in a further reduction of Salmonella enteritidis in the cecum as compared to feeding of individual SCFAs. More importantly, free fatty acids enhanced HDP gene expression without triggering proinflammatory interleukin-1β production. Taken together, oral supplementation of SCFAs is capable of boosting host immunity and disease resistance, with potential for infectious disease control and prevention in animal agriculture without relying on antibiotics.

Alternative splicing regulated by butyrate in bovine epithelial cells

As a signaling molecule and an inhibitor of histone deacetylases (HDACs), butyrate exerts its impact on a broad range of biological processes, such as apoptosis and cell proliferation, in addition to its critical role in energy metabolism in ruminants. This study examined the effect of butyrate on alternative splicing in bovine epithelial cells using RNA-seq technology. Junction reads account for 11.28 and 12.32% of total mapped reads between the butyrate-treated (BT) and control (CT) groups. 201,326 potential splicing junctions detected were supported by ≥3 junction reads. Approximately 94% of these junctions conformed to the consensus sequence (GT/AG) while ∼3% were GC/AG junctions. No AT/AC junctions were observed. A total of 2,834 exon skipping events, supported by a minimum of 3 junction reads, were detected. At least 7 genes, their mRNA expression significantly affected by butyrate, also had exon skipping events differentially regulated by butyrate. Furthermore, COL5A3, which was induced 310-fold by butyrate (FDR <0.001) at the gene level, had a significantly higher number of junction reads mapped to Exon#8 (Donor) and Exon#11 (Acceptor) in BT. This event had the potential to result in the formation of a COL5A3 mRNA isoform with 2 of the 69 exons missing. In addition, 216 differentially expressed transcript isoforms regulated by butyrate were detected. For example, Isoform 1 of ORC1 was strongly repressed by butyrate while Isoform 2 remained unchanged. Butyrate physically binds to and inhibits all zinc-dependent HDACs except HDAC6 and HDAC10. Our results provided evidence that butyrate also regulated deacetylase activities of classical HDACs via its transcriptional control. Moreover, thirteen gene fusion events differentially affected by butyrate were identified. Our results provided a snapshot into complex transcriptome dynamics regulated by butyrate, which will facilitate our understanding of the biological effects of butyrate and other HDAC inhibitors.