Analysis of Key Genes and miRNA-mRNA Networks Associated with Glucocorticoids Treatment in Chronic Obstructive Pulmonary Disease

Background

Some patients with chronic obstructive pulmonary disease (COPD) benefit from glucocorticoid (GC) treatment, but its mechanism is unclear.

Objective

With the help of the Gene Expression Omnibus (GEO) database, the key genes and miRNA-mRNA related to the treatment of COPD by GCs were discussed, and the potential mechanism was explained.

Methods

The miRNA microarray dataset (GSE76774) and mRNA microarray dataset (GSE36221) were downloaded, and differential expression analysis were performed. Gene Ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on the differentially expressed genes (DEGs). The protein interaction network of the DEGs in the regulatory network was constructed with the STRING database, and the key genes were screened through Cytoscape. Potential downstream target genes regulated by differentially expressed miRNAs (DEMs) were predicted by the miRWalk3.0 database, and miRNA-mRNA regulatory networks were constructed. Finally, some research results were validated.

Results

① Four DEMs and 83 DEGs were screened; ② GO and KEGG enrichment analysis mainly focused on the PI3K/Akt signalling pathway, ECM receptor interaction, etc.; ③ CD2, SLAMF7, etc. may be the key targets of GC in the treatment of COPD; ④ 18 intersection genes were predicted by the mirwalk 3.0 database, and 9 pairs of miRNA-mRNA regulatory networks were identified; ⑤ The expression of miR-320d-2 and TFCP2L1 were upregulated by dexamethasone in the COPD cell model, while the expression of miR-181a-2-3p and SLAMF7 were downregulated.

Conclusion

In COPD, GC may mediate the expression of the PI3K/Akt signalling pathway through miR-181a-2-3p, miR-320d-2, miR-650, and miR-155-5p, targeting its downstream signal factors. The research results provide new ideas for RNA therapy strategies of COPD, and also lay a foundation for further research.

Integrated analysis of miRNA-mRNA expression of newly emerging swine H3N2 influenza virus cross-species infection with tree shrews

BACKGROUND

Cross-species transmission of zoonotic IAVs to humans is potentially widespread and lethal, posing a great threat to human health, and their cross-species transmission mechanism has attracted much attention. miRNAs have been shown to be involved in the regulation of IAVs infection and immunity, however, few studies have focused on the molecular mechanisms underlying miRNAs and mRNAs expression after IAVs cross-species infection.

METHODS

We used tree shrews, a close relative of primates, as a model and used RNA-Seq and bioinformatics tools to analyze the expression profiles of DEMs and DEGs in the nasal turbinate tissue at different time points after the newly emerged swine influenza A virus SW2783 cross-species infection with tree shrews, and miRNA-mRNA interaction maps were constructed and verified by RT-qPCR, miRNA transfection and luciferase reporter assay.

RESULTS

14 DEMs were screened based on functional analysis and interaction map, miR-760-3p, miR-449b-2, miR-30e-3p, and miR-429 were involved in the signal transduction process of replication and proliferation after infection, miR-324-3p, miR-1301-1, miR-103-1, miR-134-5p, miR-29a, miR-31, miR-16b, miR-34a, and miR-125b participate in negative feedback regulation of genes related to the immune function of the body to activate the antiviral immune response, and miR-106b-3p may be related to the cross-species infection potential of SW2783, and the expression level of these miRNAs varies in different days after infection.

CONCLUSIONS

The miRNA regulatory networks were constructed and 14 DEMs were identified, some of them can affect the replication and proliferation of viruses by regulating signal transduction, while others can play an antiviral role by regulating the immune response. It indicates that abnormal expression of miRNAs plays a crucial role in the regulation of cross-species IAVs infection, which lays a solid foundation for further exploration of the molecular regulatory mechanism of miRNAs in IAVs cross-species infection and anti-influenza virus targets.

The Functional Role of microRNAs and mRNAs in Diabetic Kidney Disease: A Review

Diabetes is a group of diseases marked by poor control of blood glucose levels. Diabetes mellitus (DM) occurs when pancreatic cells fail to make insulin, which is required to keep blood glucose levels stable, disorders, and so on. High glucose levels in the blood induce diabetic effects, which can cause catastrophic damage to bodily organs such as the eyes and lower extremities. Diabetes is classified into many forms, one of which is controlled by hyperglycemia or Diabetic Kidney Disease (DKD), and another that is not controlled by hyperglycemia (nondiabetic kidney disease or NDKD) and is caused by other factors such as hypertension, hereditary. DKD is associated with diabetic nephropathy (DN), a leading cause of chronic kidney disease (CKD) and end-stage renal failure. The disease is characterized by glomerular basement membrane thickening, glomerular sclerosis, and mesangial expansion, resulting in a progressive decrease in glomerular filtration rate, glomerular hypertension, and renal failure or nephrotic syndrome. It is also represented by some microvascular complications such as nerve ischemia produced by intracellular metabolic changes, microvascular illness, and the direct impact of excessive blood glucose on neuronal activity. Therefore, DKD-induced nephrotic failure is worse than NDKD. MicroRNAs (miRNAs) are important in the development and progression of several diseases, including diabetic kidney disease (DKD). These dysregulated miRNAs can impact various cellular processes, including inflammation, fibrosis, oxidative stress, and apoptosis, all of which are implicated during DKD. MiRNAs can alter the course of DKD by targeting several essential mechanisms. Understanding the miRNAs implicated in DKD and their involvement in disease development might lead to identifying possible therapeutic targets for DKD prevention and therapy. Therefore, this review focuses specifically on DKD-associated DN, as well as how in-silico approaches may aid in improving the management of the disease.

Integrated analysis of a miRNA-mRNA network related to immunity and autophagy in Macrobrachium rosenbergii infected with Aeromonas hydrophila

MicroRNAs (miRNAs) are a group of RNAs that regulate gene expression in the post-transcriptionally. miRNAs can regulate numerous processes, such as the immune response, due to their dynamic expression patterns. The giant freshwater prawn Macrobrachium rosenbergii is a major freshwater aquaculture prawn that is attacked by various bacteria, including Aeromonas hydrophila. For this study, we performed an analysis of the miRNA and mRNA transcriptome analysis of M. rosenbergii which was infected with A. hydrophila. We identified 56 differentially expressed miRNAs (DEMs) and 1542 differentially expressed mRNAs. Furthermore, an integrated analysis of miRNA-mRNA expression led to the identification of 729 differentially predicted target genes (DETGs) of the DEMs. Multiple functional categories related to immunity, apoptosis, and autophagy were found to be enriched in the DETGs. During the infection of M. rosenbergii by A. hydrophila, an elaborate regulatory network involving Toll and immune deficiency (IMD) signaling, mitogen-activated protein kinase (MAPK) signaling, lysosome, and cell apoptosis was formed by a complex interplay of 40 crucial DEMs and 22 DETGs, all associated with the immune and autophagy pathway. The findings suggest that infection with A. hydrophila triggers intricate responses in both miRNA and mRNA, significantly impacting immune and autophagy processes in M. rosenbergii.

Integration of miRNA and mRNA expression profiles in Asian spongy moth Lymantria dispar in response to cyantraniliprole

The Asian spongy moth, Lymantria dispar, is a worldwide forest pest that damages >500 plant species. Nowadays, chemical control is the most widely used method because of its rapidity and effectiveness, but the insecticide resistance is a growing concern for spongy moth. As important post-transcriptional regulators of gene expression, whether microRNAs (miRNAs) are involved in insecticide tolerance is little known in spongy moth. Therefore, an integrated analysis of miRNA and mRNA was performed on L. dispar larvae treated with cyantraniliprole. Compared to the control group, a total of 432 differentially expressed genes (DEGs) and 23 differentially expressed miRNAs (DEMs) were identified in L. dispar larvae under cyantraniliprole exposure. Among them, twelve DEGs encoding detoxification enzymes/proteins were further analyzed. Twenty-one genes related to insecticide tolerance were predicted by 11 DEMs, of which 25 miRNA-mRNA interactions were identified. In the miRNA-mRNA network, novel-miR-4 and mmu-miR-3475-3p were involved in the response of L. dispar to cyantraniliprole stress by regulating five genes associated with detoxification, respectively. The P450 gene CYP4C1 (c34384.graph_c0) was the only DEG related to detoxification in the network, which was regulated by novel-miR-4. The expression levels of ten DEMs were confirmed by quantitative reverse transcription PCR (RT-qPCR) and the trends were consistent with miRNA-seq. This study identified some candidate miRNAs and mRNAs related to cyantraniliprole tolerance in L. dispar, which provides valuable transcriptomic information for revealing the molecular mechanisms of insect tolerance and developing novel insecticides.

Integrated analysis of miRNA-mRNA regulatory network and functional verification of miR-338-3p in coronary heart disease

Coronary heart disease is a cardiovascular disease with high morbidity and mortality. Although great progress has been made in treatment, the prognosis is still very poor. Therefore, this project aims to screen potential diagnostic markers and therapeutic targets related to the progression of coronary heart disease. A total of 94 overlapping differentially expressed mRNAs and 70 differentially expressed miRNAs were identified from GSE20681, GSE12288, GSE49823, and GSE105449. Through a series of bioinformatics methods and experiment, we obtained 5 core miRNA-mRNA regulatory pairs, and selected miR-338-3p/RPS23 for functional analysis. Moreover, we found that RPS23 directly targets miR-338-3p by dual luciferase assay, western, and qPCR. And the expression of miR-338-3p and RPS23 is negatively correlated. The AUC value of miR-338-3p is 0.847. Downregulation of miR-338-3p can significantly inhibit the proliferation and migration of HUVEC. On the contrary, overexpression of miR-338-3p promoted the proliferation and migration of HUVEC. In addition, the interference of RPS23 expression can reverse the regulation of miR-338-3p on HUVEC proliferation. In conclusion, miR-338-3p/RPS23 may be involved in the progression of coronary heart disease, and miR-338-3p may be a diagnostic biomarker and therapeutic target for coronary heart disease.

Identification of immune-related genes and susceptible population of pulmonary tuberculosis by constructing TF-miRNA-mRNA regulatory network

We aimed to explore the potential biomarkers and susceptible population for early diagnosis and treatment of tuberculosis (TB). Ten hub differentially expressed TB-related genes (DETRGs) from GSE83456 dataset were screened with the "limma" package and the GeneCards database. Unsupervised clustering was utilized to identify susceptible population among TB patients based on 10 hub DETRGs. TRANSFAC, MirTarbase, miRanda and TargetScan was used to predict microRNAs and transcription factors (TFs) and construct TF-miRNA-mRNA regulatory network. The results showed that a total of 266 DEGs were identified. Functional analysis mainly enriched in interferon pathway, cytokine and receptor interaction and host defense response to virus, while the four-module genes screened were closely related to interferon-γ signal transduction pathway as well. Based on 10 DETRGs, TB patients were divided into two clusters with significant differences in neutrophil function and 16 hub miRNAs and 10 hub TFs were predicted. Finally, NFATc1- (miR145) - STAT1 regulatory pathway was identified as the critical regulatory pathway, which mediates cytokine receptor binding, interleukin-1 receptor binding and TNF signaling pathway. Hence, we concluded that immunoheterogeneity exists among TB patients and NFATC1-(miR145)-STAT1 regulatory pathway might be associated with tuberculosis infection, which may be valuable targets for prevention and treatment of tuberculosis.

Integrative analysis of miRNA-mRNA network in idiopathic membranous nephropathy by bioinformatics analysis

Background

Currently, several specific antigens, M-type receptor for secretory phospholipase A2(PLA2R1), thrombospondin type-1 domain-containing 7A(THSD7A), and neural epidermal growth factor-like 1 protein (NELL-1), are discovered associated with the onset of idiopathic membranous nephropathy (IMN). But the pathomechanisms of IMN still need to be further claried. Understanding the mechanisms of IMN is required to improve its diagnosis and treatment.

Methods

In this study, we constructed miRNA regulatory networks to investigate IMN development. Moreover, miRNAs and mRNAs that were differentially expressed between Idiopathic Membranous Nephropathy (IMN) patients and normal controls were examined using the GSE115857 dataset and our previous sequence study. DE miRNA target genes were determined based on the FUNRICH software, starBase, miRDB, and miRWalk, and an miRNA-mRNA network was designed using DE-mRNAs that were negatively correlated with DE-miRNAs. The miRNA-mRNA network contained 228 miRNA-mRNA pairs. Thereafter, we conducted KEGG pathway, GO functional annotation, immune-related gene screening, protein interaction networks, and potential hub gene analyses. Furthermore, 10 miRNAs and 10 genes were determined and preliminarily validated using the validation dataset from GEO. Finally, we identified which pair may offer more accurate diagnosis and therapeutic targets for IMN.

Results

Two miRNA-mRNA pairs, miR-155-5p-FOS and miR-146a-5p-BTG2, were differentially expressed in IMN, indicating that these genes may affect IMN through immune processes. These findings may offer more accurate diagnoses and therapeutic targets for IMN.

miR-509-5p anti-infection response for mycoplasma pneumonia in sheep by targeting NF-κB pathway

MicroRNAs play a key role in Mannan-binding lectin-mediated resistance to Mycoplasma ovipneumoniae pneumonia, by regulating the translation of mRNAs of target genes, thereby regulating the immune response. Additionally, TRAF6 is a key molecule in Toll-like receptor signal transduction, which mediates inflammation and apoptosis signaling pathways and is widely involved in inflammation and immune response. While the molecular regulation mechanism has not been reported. In this study, we screened differentially expressed miRNAs and genes of Anti-infection for M. pneumonia on Sheep, through relevant bioinformatics analysis. Further, the effect of differential expression of NF-κB signaling pathway related genes on the molecular mechanism of M. pneumonia was detected. We used miRNA-mRNA integrated analysed, the target gene TRAF6 of miR-509-5p was selected. TRAF6 dual luciferase reporter vector was co-transfected into HEK 293T cells and primary sheep respiratory mucosal epithelial cells to detect changes in luciferase activity. qRT-PCR was used to analyze the effect of miR-509-5p on the expression and regulation of TRAF6 and other genes related to the NF-κB signaling pathway. The result confirmed that TRAF6 was a target gene of miR-509-5p. Compared with miR-509-5p-NC group, the luciferase activity of miR-509-5p group was significantly down-regulated (P < 0.01). Further, in sheep respiratory mucosal epithelial cells, miR-509-5p mimic could significantly down-regulate the fold change value of TRAF6 (P < 0.01). On the contrary, miR-509-5p-inhibitor up-regulated the fold change value of TRAF6 (P < 0.05). Interestingly, the expression levels of other genes were different. Among them, miR-509-5p mimic significantly up-regulated TLR4 and IRAK4 (P < 0.05), significantly down-regulated TAK1 (P < 0.05) and NF-κB (P < 0.01). miR-509-5p-inhibitor significantly up-regulated NF-κB (P < 0.05) and TAK1 (P < 0.01). miR-509-5p targets TRAF6 to affect the expression of downstream genes, which negatively regulates the NF-κB pathway, thereby affecting the inflammatory response.

Five Circular RNAs in Metabolism Pathways Related to Prostate Cancer

Prostate cancer (PCa) is the most common malignant tumor in men, and its incidence increases with age. Serum prostate-specific antigen and tissue biopsy remain the standard for diagnosis of suspected PCa. However, these clinical indicators may lead to aggressive overtreatment in patients who have been treated sufficiently with active surveillance. Circular RNAs (circRNAs) have been recently recognized as a new type of regulatory RNA that is not easily degraded by RNases and other exonucleases because of their covalent closed cyclic structure. Thus, we utilized high-throughput sequencing data and bioinformatics analysis to identify specifically expressed circRNAs in PCa and filtered out five specific circRNAs for further analysis—hsa_circ_0006410, hsa_circ_0003970, hsa_circ_0006754, hsa_circ_0005848, and a novel circRNA, hsa_circ_AKAP7. We constructed a circRNA-miRNA regulatory network and used miRNA and differentially expressed mRNA interactions to predict the function of the selected circRNAs. Furthermore, survival analysis of their cognate genes and PCR verification of these five circRNAs revealed that they are closely related to well-known PCa pathways such as the MAPK signaling pathway, P53 pathway, androgen receptor signaling pathway, cell cycle, hormone-mediated signaling pathway, and cellular lipid metabolic process. By understanding the related metabolism of circRNAs, these circRNAs could act as metabolic biomarkers, and monitoring their levels could help diagnose PCa. Meanwhile, the exact regulatory mechanism for AR-related regulation in PCa is still unclear. The circRNAs we found can provide new solutions for research in this field.

Integrated transcriptome, small RNA and degradome sequencing approaches provide insights into Ascochyta blight resistance in chickpea

Ascochyta blight (AB) is one of the major biotic stresses known to limit the chickpea production worldwide. To dissect the complex mechanisms of AB resistance in chickpea, three approaches, namely, transcriptome, small RNA and degradome sequencing were used. The transcriptome sequencing of 20 samples including two resistant genotypes, two susceptible genotypes and one introgression line under control and stress conditions at two time points (3rd and 7th day post inoculation) identified a total of 6767 differentially expressed genes (DEGs). These DEGs were mainly related to pathogenesis-related proteins, disease resistance genes like NBS-LRR, cell wall biosynthesis and various secondary metabolite synthesis genes. The small RNA sequencing of the samples resulted in the identification of 651 miRNAs which included 478 known and 173 novel miRNAs. A total of 297 miRNAs were differentially expressed between different genotypes, conditions and time points. Using degradome sequencing and in silico approaches, 2131 targets were predicted for 629 miRNAs. The combined analysis of both small RNA and transcriptome datasets identified 12 miRNA-mRNA interaction pairs that exhibited contrasting expression in resistant and susceptible genotypes and also, a subset of genes that might be post-transcriptionally silenced during AB infection. The comprehensive integrated analysis in the study provides better insights into the transcriptome dynamics and regulatory network components associated with AB stress in chickpea and, also offers candidate genes for chickpea improvement.

Analysis of clock gene-miRNA correlation networks reveals candidate drivers in colorectal cancer

Altered functioning of the biological clock is involved in cancer onset and progression. MicroRNAs (miRNAs) interact with the clock genes modulating the function of genetically encoded molecular clockworks. Collaborative interactions may take place within the coding-noncoding RNA regulatory networks. We aimed to evaluate the cross-talk among miRNAs and clock genes in colorectal cancer (CRC). We performed an integrative analysis of miRNA-miRNA and miRNA-mRNA interactions on high-throughput molecular profiling of matched human CRC tissue and non-tumor mucosa, pinpointing core clock genes and their targeting miRNAs. Data obtained in silico were validated in CRC patients and human colon cancer cell lines. In silico we found severe alterations of clock gene–related coding-noncoding RNA regulatory networks in tumor tissues, which were later corroborated by the analysis of human CRC specimens and experiments performed in vitro. In conclusion, specific miRNAs target and regulate the transcription/translation of clock genes and clock gene-related miRNA-miRNA as well as mRNA-miRNA interactions are altered in colorectal cancer. Exploration of the interplay between specific miRNAs and genes, which are critically involved in the functioning of the biological clock, provides a better understanding of the importance of the miRNA-clock genes axis and its derangement in colorectal cancer.

Micronome revealed miR-19a/b as key regulator of SOCS3 during cancer related inflammation of oral squamous cell carcinoma

Although significant advances have been established in molecular biology of Oral squamous cell carcinoma (OSCC), innovative strategies are still required to further understand detailed molecular mechanisms. Using bioinformatic approach, we aim to explore the potential miRNA-mRNA pairs in cancer related inflammatory response and investigate their potential roles as signature miRNA and proteins in the signaling pathway. Firstly, the differentially expressed genes of OSCC were selected which then underwent gene ontology to identify genes engaged in inflammatory response and its regulation. Validated miRNAs were retrieved and miRNAs with complete complementarily with their targets were visualized for miRNA-mRNA regulatory network. Protein-protein interactions of inflammatory and its regulatory genes were analyzed for interacting genes involved in signaling pathway. Eight universal miRNAs were obtained for inflammation and its regulation. miRNA-19a/b showed significant influence in controlling inflammatory response in OSCC. Therefore, micronome on deregulated genes in inflammation identifies miRNA-mRNA pairs which have high potential to be targeted for diagnostic and treatment applications in OSCC.

Association Studies in Populus tomentosa Reveal the Genetic Interactions of Pto-MIR156c and Its Targets in Wood Formation

MicroRNAs (miRNAs) regulate gene expression in many biological processes, but the significance of the interaction between a miRNA and its targets in perennial trees remains largely unknown. Here, we employed transcript profiling and association studies in Populus tomentosa (Pto) to decipher the effect of genetic variation and interactions between Pto-miR156c and its potential targets (Pto-SPL15, Pto-SPL20, and Pto-SPL25) in 435 unrelated individuals from a natural population of P. tomentosa. Single-SNP (single-nucleotide polymorphism) based association studies with analysis of the underlying additive and dominant effects identified 69 significant associations (P < 0.01), representing 51 common SNPs (minor allele frequency > 0.05) from Pto-MIR156c and its three potential targets, with six wood and growth traits, revealing their common roles in wood formation. Epistasis analysis uncovered 129 significant SNP-SNP associations with ten traits, indicating the potential genetic interactions of Pto-MIR156c and its three putative targets. Interestingly, expression analysis in stem (phloem, cambium, and xylem) revealed that Pto-miR156c expression showed strong negative correlations with Pto-SPL20 (r = -0.90, P < 0.01) and Pto-SPL25 (r = -0.65, P < 0.01), and a positive correlation with Pto-SPL15 (r = 0.40, P < 0.01), which also indicated the putative interactions of Pto-miR156c and its potential targets and their common roles in wood formation. Thus, our study provided an alternative approach to decipher the interaction between miRNAs and their targets and to dissect the genetic architecture of complex traits in trees.

Computational and in vitro Investigation of miRNA-Gene Regulations in Retinoblastoma Pathogenesis: miRNA Mimics Strategy

PURPOSE

Retinoblastoma (RB), a primary pediatric intraocular tumor, arises from primitive retinal layers. Several novel molecular strategies are being developed for the clinical management of RB. miRNAs are known to regulate cancer-relevant biological processes. Here, the role of selected miRNAs, namely, miR-532-5p and miR-486-3p, has been analyzed for potential therapeutic targeting in RB.

METHODS

A comprehensive bioinformatic analysis was performed to predict the posttranscriptional regulators (miRNAs) of the select panel of genes [Group 1: oncogenes (HMGA2, MYCN, SYK, FASN); Group 2: cancer stem cell markers (TACSTD, ABCG2, CD133, CD44, CD24) and Group 3: cell cycle regulatory proteins (p53, MDM2)] using Microcosm, DIANALAB, miRBase v 18, and REFSEQ database, and RNA hybrid. The expressions of five miRNAs, namely, miR-146b-5p, miR-532-5p, miR-142-5p, miR-328, and miR-486-3p, were analyzed by qRT–PCR on primary RB tumor samples (n = 30; including 17 invasive RB tumors and 13 noninvasive RB tumors). Detailed complementary alignment between 5’ seed sequence of differentially expressed miRNAs and the sequence of target genes was determined. Based on minimum energy level and piCTAR scores, the gene targets were selected. Functional roles of these miRNA clusters were studied by using mimics in cultured RB (Y79, Weri Rb-1) cells in vitro. The gene targets (SYK and FASN) of the studied miRNAs were confirmed by qRT-PCR and western blot analysis. Cell proliferation and apoptotic studies were performed.

RESULTS

Nearly 1948 miRNAs were identified in the in silico analysis, From this list, only 9 upregulated miRNAs (miR-146b-5p, miR-305, miR-663b, miR-299, miR-532-5p, miR-892b, miR-501, miR-142-5p, and miR-513b) and 10 downregulated miRNAs (miR-1254, miR-328, miR-133a, miR-1287, miR-1299, miR-375, miR-486-3p, miR-720, miR-98, and miR-122*) were found to be common with the RB serum miRNA profile. Downregulation of five miRNAs (miR-146b-5p, miR-532-5p, miR-142-5p, miR-328, and miR-486-3p) was confirmed experimentally. Predicted common oncogene targets (SYK and FASN) of miR-486-3p and miR-532-5p were evaluated for their mRNA and protein expression in these miRNA mimic-treated RB cells. Experimental overexpression of these miRNAs mediated apoptotic cell death without significantly altering the cell cycle in RB cells.

CONCLUSION

Key miRNAs in RB pathogenesis were identified by an in silico approach. Downregulation of miR-486-3p and miR-532-5p in primary retinoblastoma tissues implicates their role in tumorigenesis. Prognostic and therapeutic potential of these miRNA was established by the miRNA mimic strategy.