An integrated analysis of mRNAs, lncRNAs, and miRNAs based on weighted gene co-expression network analysis involved in bovine endometritis

Enhanced expression of miR-204 attenuates LPS stimulated inflammatory injury through inhibiting the Wnt/β-catenin pathway via targeting CCND2

One of the most common bacterial diseases of the reproductive system in dairy cows is endometritis, which will cause huge economic loss. Here, we investigate the mechanisms of miR-204 on LPS-stimulated endometritis in vitro and in vivo. Experiments displayed that the expression of miR-204 was lower in bovine uterine tissue samples or bovine endometrial epithelial cell line (BEND) that stimulated by LPS. Compared with the negative group, miR-204 treatment significantly suppressed the production of proinflammatory factors and the Wnt/β-catenin pathway activation. Additionally, the result of the dual luciferase assay showed that miR-204 targeted cyclin D2. More importantly, up-regulation of miR-204 alleviated LPS induced uterine injury was confirmed in vivo studies. Molecular experiments indicated that the expression level of tight junctional proteins Claudin3 and cadherin1 were both enchanced by miR-204 treatment. Accordingly, miR-204 may serve as a new measure to prevent and treat endometritis caused by LPS.

Comprehensive transcriptomic analysis of long non-coding RNAs in bovine ovarian follicles and early embryos

Long non-coding RNAs (lncRNAs) have been the subject of numerous studies over the past decade. First thought to come from aberrant transcriptional events, lncRNAs are now considered a crucial component of the genome with roles in multiple cellular functions. However, the functional annotation and characterization of bovine lncRNAs during early development remain limited. In this comprehensive analysis, we review lncRNAs expression in bovine ovarian follicles and early embryos, based on a unique database comprising 468 microarray hybridizations from a single platform designed to target 7,724 lncRNA transcripts, of which 5,272 are intergenic (lincRNA), 958 are intronic, and 1,524 are antisense (lncNAT). Compared to translated mRNA, lncRNAs have been shown to be more tissue-specific and expressed in low copy numbers. This analysis revealed that protein-coding genes and lncRNAs are both expressed more in oocytes. Differences between the oocyte and the 2-cell embryo are also more apparent in terms of lncRNAs than mRNAs. Co-expression network analysis using WGCNA generated 25 modules with differing proportions of lncRNAs. The modules exhibiting a higher proportion of lncRNAs were found to be associated with fewer annotated mRNAs and housekeeping functions. Functional annotation of co-expressed mRNAs allowed attribution of lncRNAs to a wide array of key cellular events such as meiosis, translation initiation, immune response, and mitochondrial related functions. We thus provide evidence that lncRNAs play diverse physiological roles that are tissue-specific and associated with key cellular functions alongside mRNAs in bovine ovarian follicles and early embryos. This contributes to add lncRNAs as active molecules in the complex regulatory networks driving folliculogenesis, oogenesis and early embryogenesis all of which are necessary for reproductive success.

Integrated co-expression analysis of regulatory elements (miRNA, lncRNA, and TFs) in bovine monocytes induced by Str. uberis

Non-coding RNAs, including long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), together with transcription factors, are critical pre-, co-, and post-transcriptional regulators. In addition to their criteria as ideal biomarkers, they have great potential in disease prognosis, diagnosis, and treatment of complex diseases. Investigation of regulatory mechanisms in the context of bovine mastitis, as most common and economic disease in the dairy industry, to identify elements influencing the expression of candidate genes as key regulators of the mammary immune response is not yet fully understood. Transcriptome profiles (50 RNA-Seq and 50 miRNA-Seq samples) of bovine monocytes induced by Str. uberis were used for co-expression module detection and preservation analysis using the weighted gene co-expression network analysis (WGCNA) approach. Assigned mi-, lnc-, and m-modules used to construct the integrated regulatory networks and miRNA-lncRNA-mRNA regulatory sub-networks. Remarkably, we have identified 18 miRNAs, five lncRNAs, and seven TFs as key regulators of str. uberis-induced mastitis. Most of the genes introduced here, mainly involved in immune response, inflammation, and apoptosis, were new to mastitis. These findings may help to further elucidate the underlying mechanisms of bovine mastitis, and the discovered genes may serve as signatures for early diagnosis and treatment of the disease.

Integrated analysis of inflammatory mRNAs, miRNAs, and lncRNAs elucidates the molecular interactome behind bovine mastitis

Mastitis is known as intramammary inflammation, which has a multifactorial complex phenotype. However, the underlying molecular pathogenesis of mastitis remains poorly understood. In this study, we utilized a combination of RNA-seq and miRNA-seq techniques, along with computational systems biology approaches, to gain a deeper understanding of the molecular interactome involved in mastitis. We retrieved and processed one hundred transcriptomic libraries, consisting of 50 RNA-seq and 50 matched miRNA-seq data, obtained from milk-isolated monocytes of Holstein-Friesian cows, both infected with Streptococcus uberis and non-infected controls. Using the weighted gene co-expression network analysis (WGCNA) approach, we constructed co-expressed RNA-seq-based and miRNA-seq-based modules separately. Module-trait relationship analysis was then performed on the RNA-seq-based modules to identify highly-correlated modules associated with clinical traits of mastitis. Functional enrichment analysis was conducted to understand the functional behavior of these modules. Additionally, we assigned the RNA-seq-based modules to the miRNA-seq-based modules and constructed an integrated regulatory network based on the modules of interest. To enhance the reliability of our findings, we conducted further analyses, including hub RNA detection, protein-protein interaction (PPI) network construction, screening of hub-hub RNAs, and target prediction analysis on the detected modules. We identified a total of 17 RNA-seq-based modules and 3 miRNA-seq-based modules. Among the significant highly-correlated RNA-seq-based modules, six modules showed strong associations with clinical characteristics of mastitis. Functional enrichment analysis revealed that the turquoise module was directly related to inflammation persistence and mastitis development. Furthermore, module assignment analysis demonstrated that the blue miRNA-seq-based module post-transcriptionally regulates the turquoise RNA-seq-based module. We also identified a set of different RNAs, including hub-hub genes, hub-hub TFs (transcription factors), hub-hub lncRNAs (long non-coding RNAs), and hub miRNAs within the modules of interest, indicating their central role in the molecular interactome underlying the pathogenic mechanisms of S. uberis infection. This study provides a comprehensive insight into the molecular crosstalk between immunoregulatory mRNAs, miRNAs, and lncRNAs during S. uberis infection. These findings offer valuable directions for the development of molecular diagnosis and biological therapies for mastitis.

Uterine Tissue Metabonomics Combined with 16S rRNA Gene Sequencing To Analyze the Changes of Gut Microbiota in Mice with Endometritis and the Intervention Effect of Tau Interferon

Endometritis is a common cow disease characterized by inflammation of endometrium, which leads to infertility or low fertility of cows and brings huge economic losses to the dairy industry. Tau interferon (IFN-τ) has many important biological functions, including an anti-inflammatory effect. The present study aimed to survey the effects of IFN-τ administration on gut microflora and body metabolism in mice with endometritis and to explore the potential relationship. The results indicated that IFN-τ obviously alleviated the damage and ultrastructural changes of mouse endometrium induced by Escherichia coli and enhanced tight junction protein's expression level. Through analysis by 16S rRNA gene sequencing, we found that IFN-τ, especially at 12 h, could regulate the composition of gut microbiota associated with Pediococcus, Staphylococcus, and Enterorhabdus in E. coli-induced mouse endometritis. Through histometabonomics, it was found that endometritis was related to 11 different metabolites and 4 potential metabolic pathways. These metabolites and metabolic pathways were major participants in metabolic pathways, cysteine and methionine metabolism, arachidonic acid metabolism, and pyrimidine metabolism. Correlation analysis of gut microbiota with uterine tissue metabolomics showed that changes in metabolic pathways might be affected by gut microbiota, such as Enterorhabdus in mouse endometritis. The above results indicated that the anti-inflammatory mechanism of IFN-τ might be reduction of the abundance of Enterorhabdus in the gut microbiota, affecting the expression level of important metabolites in uterine tissue and thus playing an anti-inflammatory role. IMPORTANCE The change in intestinal flora has been the focus of many disease studies in recent years, but the pathogenetic effect of interferon on endometritis is still unclear. The results of this study showed that IFN-τ alleviated the damage in mouse endometritis induced by E. coli and improved the endometrial tissue barrier. Its functional mechanism may be reduction of the abundance of Enterorhabdus in the intestinal microbiota, affecting the expression level of important metabolites in uterine tissue and thus playing an anti-inflammatory role.

In-depth systems biological evaluation of bovine alveolar macrophages suggests novel insights into molecular mechanisms underlying Mycobacterium bovis infection

OBJECTIVE

Bovine tuberculosis (bTB) is a chronic respiratory infectious disease of domestic livestock caused by intracellular Mycobacterium bovis infection, which causes ~$3 billion in annual losses to global agriculture. Providing novel tools for bTB managements requires a comprehensive understanding of the molecular regulatory mechanisms underlying the M. bovis infection. Nevertheless, a combination of different bioinformatics and systems biology methods was used in this study in order to clearly understand the molecular regulatory mechanisms of bTB, especially the immunomodulatory mechanisms of M. bovis infection.

METHODS

RNA-seq data were retrieved and processed from 78 (39 non-infected control vs. 39 M. bovis-infected samples) bovine alveolar macrophages (bAMs). Next, weighted gene co-expression network analysis (WGCNA) was performed to identify the co-expression modules in non-infected control bAMs as reference set. The WGCNA module preservation approach was then used to identify non-preserved modules between non-infected controls and M. bovis-infected samples (test set). Additionally, functional enrichment analysis was used to investigate the biological behavior of the non-preserved modules and to identify bTB-specific non-preserved modules. Co-expressed hub genes were identified based on module membership (MM) criteria of WGCNA in the non-preserved modules and then integrated with protein-protein interaction (PPI) networks to identify co-expressed hub genes/transcription factors (TFs) with the highest maximal clique centrality (MCC) score (hub-central genes).

RESULTS

As result, WGCNA analysis led to the identification of 21 modules in the non-infected control bAMs (reference set), among which the topological properties of 14 modules were altered in the M. bovis-infected bAMs (test set). Interestingly, 7 of the 14 non-preserved modules were directly related to the molecular mechanisms underlying the host immune response, immunosuppressive mechanisms of M. bovis, and bTB development. Moreover, among the co-expressed hub genes and TFs of the bTB-specific non-preserved modules, 260 genes/TFs had double centrality in both co-expression and PPI networks and played a crucial role in bAMs-M. bovis interactions. Some of these hub-central genes/TFs, including PSMC4, SRC, BCL2L1, VPS11, MDM2, IRF1, CDKN1A, NLRP3, TLR2, MMP9, ZAP70, LCK, TNF, CCL4, MMP1, CTLA4, ITK, IL6, IL1A, IL1B, CCL20, CD3E, NFKB1, EDN1, STAT1, TIMP1, PTGS2, TNFAIP3, BIRC3, MAPK8, VEGFA, VPS18, ICAM1, TBK1, CTSS, IL10, ACAA1, VPS33B, and HIF1A, had potential targets for inducing immunomodulatory mechanisms by M. bovis to evade the host defense response.

CONCLUSION

The present study provides an in-depth insight into the molecular regulatory mechanisms behind M. bovis infection through biological investigation of the candidate non-preserved modules directly related to bTB development. Furthermore, several hub-central genes/TFs were identified that were significant in determining the fate of M. bovis infection and could be promising targets for developing novel anti-bTB therapies and diagnosis strategies.

Identification of the hub genes related to adipose tissue metabolism of bovine

Due to the demand for high-quality animal protein, there has been consistent interest in how to obtain more high-quality beef. As well-known, the adipose content of beef has a close connection with the taste and quality of beef, and cattle with different energy or protein diet have corresponding effects on the lipid metabolism of beef. Thus, we performed weighted gene co-expression network analysis (WGCNA) with subcutaneous adipose genes from Norwegian red heifers fed different diets to identify hub genes regulating bovine lipid metabolism. For this purpose, the RNA sequencing data of subcutaneous adipose tissue of 12-month-old Norwegian red heifers (n = 48) with different energy or protein levels were selected from the GEO database, and 7,630 genes with the largest variation were selected for WGCNA analysis. Then, three modules were selected as hub genes candidate modules according to the correlation between modules and phenotypes, including pink, magenta and grey60 modules. GO and KEGG enrichment analysis showed that genes were related to metabolism, and participated in Rap, MAPK, AMPK, VEGF signaling pathways, and so forth. Combined gene interaction network analysis using Cytoscape software, eight hub genes of lipid metabolism were identified, including TIA1, LOC516108, SNAPC4, CPSF2, ZNF574, CLASRP, MED15 and U2AF2. Further, the expression levels of hub genes in the cattle tissue were also measured to verify the results, and we found hub genes in higher expression in muscle and adipose tissue in adult cattle. In summary, we predicted the key genes of lipid metabolism in the subcutaneous adipose tissue that were affected by the intake of various energy diets to find the hub genes that coordinate lipid metabolism, which provide a theoretical basis for regulating beef quality.

In Silico Identification of lncRNAs Regulating Sperm Motility in the Turkey (Meleagris gallopavo L.)

Long non-coding RNAs (lncRNAs) are transcripts not translated into proteins with a length of more than 200 bp. LncRNAs are considered an important factor in the regulation of countless biological processes, mainly through the regulation of gene expression and interactions with proteins. However, the detailed mechanism of interaction as well as functions of lncRNAs are still unclear and therefore constitute a serious research challenge. In this study, for the first time, potential mechanisms of lncRNA regulation of processes related to sperm motility in turkey were investigated and described. Customized bioinformatics analysis was used to detect and identify lncRNAs, and their correlations with differentially expressed genes and proteins were also investigated. Results revealed the expression of 863 new/unknown lncRNAs in ductus deferens, testes and epididymis of turkeys. Moreover, potential relationships of the lncRNAs with the coding mRNAs and their products were identified in turkey reproductive tissues. The results obtained from the OMICS study may be useful in describing and characterizing the way that lncRNAs regulate genes and proteins as well as signaling pathways related to sperm motility.

Weighted Gene Co-expression Network Analysis Identifies Specific Modules and Hub Genes Related to Subacute Ruminal Acidosis

Weighted gene co-expression network analysis (WGCNA) was used to understand the pathogenesis of subacute ruminal acidosis (SARA) and identify potential genes related to the disease. Microarray data from dataset GSE143765, which included 22 cows with and nine cows without SARA, were downloaded from the NCBI Gene Expression Omnibus (GEO). Results of WGCNA identified highly correlated modules of sample genes, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses allowed further biological insights into SARA-related modules. The protein-protein interaction (PPI) network, modules from the PPI network, and cistron annotation enrichment of modules were also analyzed. A total of 14,590 DEGs were used for the WGCNA. Construction of a protein-protein network identified DCXR, MMP15, and MMP17 as hub genes. Functional annotation showed that DCXR mainly exhibited L-xylulose reductase (NADP+) activity, glucose metabolic process, xylulose metabolic process, and carbonyl reductase (NADPH) activity, which are involved in the pentose and glucuronate interconversion pathways. MMP15 and MMP17 mainly have a collagen catabolic process. Overall, the results of this study aid the clarification of the biological and metabolic processes associated with SARA at the molecular level. The data highlight potential mechanisms for the future development of intervention strategies to reduce or alleviate the risk of SARA.

miR-484: A Potential Biomarker in Health and Disease

Disorders of miR-484 expression are observed in cancer, different diseases or pathological states. There is accumulating evidence that miR-484 plays an essential role in the development as well as the regression of different diseases, and miR-484 has been reported as a key regulator of common cancer and non-cancer diseases. The miR-484 targets that have effects on inflammation, apoptosis and mitochondrial function include SMAD7, Fis1, YAP1 and BCL2L13. For cancer, identified targets include VEGFB, VEGFR2, MAP2, MMP14, HNF1A, TUSC5 and KLF12. The effects of miR-484 on these targets have been documented separately. Moreover, miR-484 is typically described as an oncosuppressor, but this claim is simplistic and one-sided. This review will combine relevant basic and clinical studies to find that miR-484 promotes tumorigenesis and metastasis in liver, prostate and lung tissues. It will provide a basis for the possible mechanisms of miR-484 in early tumor diagnosis, prognosis determination, disease assessment, and as a potential therapeutic target for tumors.