DNA methylation and gene expression profiles to identify childhood atopic asthma associated genes

Background Asthma is a chronic inflammatory disorder of the airways involving many different factors. This study aimed to screen for the critical genes using DNA methylation/CpGs and miRNAs involved in childhood atopic asthma. Methods DNA methylation and gene expression data (Access Numbers GSE40732 and GSE40576) were downloaded from the Gene Expression Omnibus database. Each set contains 194 peripheral blood mononuclear cell (PBMC) samples of 97 children with atopic asthma and 97 control children. Differentially expressed genes (DEGs) with DNA methylation changes were identified. Pearson correlation analysis was used to select genes with an opposite direction of expression and differences in methylation levels, and then Gene Ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed. Protein–protein interaction network and miRNA–target gene regulatory networks were then constructed. Finally, important genes related to asthma were screened. Results A total of 130 critical DEGs with DNA methylation changes were screened from children with atopic asthma and compared with control samples from healthy children. GO and KEGG pathway enrichment analysis found that critical genes were primarily related to 24 GO terms and 10 KEGG pathways. In the miRNA–target gene regulatory networks, 9 KEGG pathways were identified. Analysis of the miRNA–target gene network noted an overlapping KEGG signaling pathway, hsa04060: cytokine-cytokine receptor interaction, in which the gene CCL2, directly related to asthma, was involved. This gene is targeted by eight asthma related miRNAs (hsa-miR-206, hsa-miR-19a, hsa-miR-9,hsa-miR-22, hsa-miR-33b, hsa-miR-122, hsa-miR-1, and hsa-miR-23b). The genes IL2RG and CCl4 were also involved in this pathway. Conclusions The present study provides a novel insight into the underlying molecular mechanism of childhood atopic asthma.

Identification of genes, pathways and transcription factor-miRNA-target gene networks and experimental verification in venous thromboembolism

Venous thromboembolism (VTE) is a complex, multifactorial life-threatening disease that involves vascular endothelial cell (VEC) dysfunction. However, the exact pathogenesis and underlying mechanisms of VTE are not completely clear. The aim of this study was to identify the core genes and pathways in VECs that are involved in the development and progression of unprovoked VTE (uVTE). The microarray dataset GSE118259 was downloaded from the Gene Expression Omnibus database, and 341 up-regulated and 8 down-regulated genes were identified in the VTE patients relative to the healthy controls, including CREB1, HIF1α, CBL, ILK, ESM1 and the ribosomal protein family genes. The protein–protein interaction (PPI) network and the transcription factor (TF)-miRNA-target gene network were constructed with these differentially expressed genes (DEGs), and visualized using Cytoscape software 3.6.1. Eighty-nine miRNAs were predicted as the targeting miRNAs of the DEGs, and 197 TFs were predicted as regulators of these miRNAs. In addition, 237 node genes and 4 modules were identified in the PPI network. The significantly enriched pathways included metabolic, cell adhesion, cell proliferation and cellular response to growth factor stimulus pathways. CREB1 was a differentially expressed TF in the TF-miRNA-target gene network, which regulated six miRNA-target gene pairs. The up-regulation of ESM1, HIF1α and CREB1 was confirmed at the mRNA and protein level in the plasma of uVTE patients. Taken together, ESM1, HIF1α and the CREB1-miRNA-target genes axis play potential mechanistic roles in uVTE development.

Genome-Wide Analysis of Coding and Non-coding RNA Reveals a Conserved miR164–NAC–mRNA Regulatory Pathway for Disease Defense in Populus

MicroRNAs (miRNAs) contribute to plant defense responses by increasing the overall genetic diversity; however, their origins and functional importance in plant defense remain unclear. Here, we employed Illumina sequencing technology to assess how miRNA and messenger RNA (mRNA) populations vary in the Chinese white poplar (Populus tomentosa) during a leaf black spot fungus (Marssonina brunnea) infection. We sampled RNAs from infective leaves at conidia germinated stage [12 h post-inoculation (hpi)], infective vesicles stage (24 hpi), and intercellular infective hyphae stage (48 hpi), three essential stages associated with plant colonization and biotrophic growth in M. brunnea fungi. In total, 8,938 conserved miRNA-target gene pairs and 3,901 Populus-specific miRNA-target gene pairs were detected. The result showed that Populus-specific miRNAs (66%) were more involved in the regulation of the disease resistance genes. By contrast, conserved miRNAs (>80%) target more whole-genome duplication (WGD)-derived transcription factors (TFs). Among the 1,023 WGD-derived TF pairs, 44.9% TF pairs had only one paralog being targeted by a miRNA that could be due to either gain or loss of a miRNA binding site after the WGD. A conserved hierarchical regulatory network combining promoter analyses and hierarchical clustering approach uncovered a miR164–NAM, ATAF, and CUC (NAC) transcription factor–mRNA regulatory module that has potential in Marssonina defense responses. Furthermore, analyses of the locations of miRNA precursor sequences reveal that pseudogenes and transposon contributed a certain proportion (∼30%) of the miRNA origin. Together, these observations provide evolutionary insights into the origin and potential roles of miRNAs in plant defense and functional innovation.

Association of Extracellular Signal-Regulated Kinase Genes With Myopia: A Longitudinal Study of Chinese Children

Graphical AbstractThe genes in the miRNA-target gene network represent the intersection of the target genes and the genes from String that had direct or indirect interaction relationships with significant genes.

MicroRNA-Mediated Responses to Cadmium Stress in Arabidopsis thaliana

In recent decades, the presence of cadmium (Cd) in the environment has increased significantly due to anthropogenic activities. Cd is taken up from the soil by plant roots for its subsequent translocation to shoots. However, Cd is a non-essential heavy metal and is therefore toxic to plants when it over-accumulates. MicroRNA (miRNA)-directed gene expression regulation is central to the response of a plant to Cd stress. Here, we document the miRNA-directed response of wild-type Arabidopsis thaliana (Arabidopsis) plants and the drb1, drb2 and drb4 mutant lines to Cd stress. Phenotypic and physiological analyses revealed the drb1 mutant to display the highest degree of tolerance to the imposed stress while the drb2 mutant was the most sensitive. RT-qPCR-based molecular profiling of miRNA abundance and miRNA target gene expression revealed DRB1 to be the primary double-stranded RNA binding (DRB) protein required for the production of six of the seven Cd-responsive miRNAs analyzed. However, DRB2, and not DRB1, was determined to be required for miR396 production. RT-qPCR further inferred that transcript cleavage was the RNA silencing mechanism directed by each assessed miRNA to control miRNA target gene expression. Taken together, the results presented here reveal the complexity of the miRNA-directed molecular response of Arabidopsis to Cd stress.

Pinpointing miRNA and genes enrichment over trait-relevant tissue network in Genome-Wide Association Studies

Background Understanding gene regulation is important but difficult. Elucidating tissue-specific gene regulation mechanism is even more challenging and requires gene co-expression network assembled from protein–protein interaction, transcription factor and gene binding, and post-transcriptional regulation (e.g., miRNA targeting) information. The miRNA binding affinity could therefore be changed by SNP(s) located at the 3′ untranslated regions (3′UTR) of the target messenger RNA (mRNA) which miRNA(s) interacts with. Genome-wide association study (GWAS) has reported significant numbers of loci hosting SNPs associated with many traits. The goal of this study is to pinpoint GWAS functional variants located in 3′UTRs and elucidate if the genes harboring these variants along with their targeting miRNAs are associated with genetic traits relevant to certain tissues. Methods By applying MIGWAS, CoCoNet, ANNOVAR, and DAVID bioinformatics software and utilizing the gene expression database (e.g. GTEx data) to study GWAS summary statistics for 43 traits from 28 GWAS studies, we have identified a list of miRNAs and targeted genes harboring 3′UTR variants, which could contribute to trait-relevant tissue over miRNA-target gene network. Results Our result demonstrated that strong association between traits and tissues exists, and in particular, the Primary Biliary Cirrhosis (PBC) trait has the most significant p-value for all 180 tissues among all 43 traits used for this study. We reported SNPs located in 3′UTR regions of genes (SFMBT2, ZC3HAV1, and UGT3A1) targeted by miRNAs for PBC trait and its tissue association network. After employing Gene Ontology (GO) analysis for PBC trait, we have also identified a very important miRNA targeted gene over miRNA-target gene network, PFKL, which encodes the liver subunit of an enzyme. Conclusions The non-coding variants identified from GWAS studies are casually assumed to be not critical to translated protein product. However, 3′ untranslated regions (3′UTRs) of genes harbor variants can often change the binding affinity of targeting miRNAs playing important roles in protein translation degree. Our study has shown that GWAS variants could play important roles on miRNA-target gene networks by contributing the association between traits and tissues. Our analysis expands our knowledge on trait-relevant tissue network and paves way for future human disease studies.

Identification of Novel Transcription Factor-microRNA-mRNA Co-regulatory Networks in Pulmonary Large-cell Neuroendocrine Carcinomas

Background: Large-cell neuroendocrine carcinoma (LCNEC) of the lung is a rare neuroendocrine neoplasm. Previous studies have shown that microRNAs can be widely involved in tumor regulation through targeting critical genes. However, it is unclear about which miRNAs play vital roles in LCNEC’s pathogenesis and how they interact with transcription factors (TFs) to regulate the cancer-related genes. Methods: To clarify the new TFs-miRNA-target gene feed-forward loops (FFLs) model of LCNEC, we integrated multi-omics data from GEO, TRRUST, TRED, and miRTarBase database. First, expression profile datasets for mRNAs (GSE1037) and miRNAs (GSE19945) were downloaded from the GEO database. Overlapping differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs) were identified through integrative analysis. The target genes of the FFL being obtained from the mirTarBase database, the GO and KEGG functional enrichment analysis was performed on the target gene. Then, we screened for key miRNAs in the FFL and perform gene regulatory network analysis based on key miRNAs. Finally, the TF-miRNA-target gene FFLs were constructed by the hypergeometric test.Results: A total of 343 DEGs and 60 DEMs were identified in LCNECs compared to normal tissues, including 210 downregulated, 133 upregulated genes and 29 downregulated, 31 upregulated miRNAs. Additionally, 55 DEMs were predicted through miRanda, mirDB, and TargetScan to be associated with the LCNEC 77 key genes. Finally, ETS1-miR195-CD36 and E2F3-miR195-CD36 from the TF-miRNA-mRNA FFLs were constructed and are significantly enriched in focal adhesion, cytokine-receptor interaction, hematopoietic cell lineage, and transforming growth factor (TGF)-β signaling pathways. Survival analysis showed that the differential expression of the multi-effect factors ETS1, CD36, and hsa-miR-195 in the FFL significantly affected the survival time of LCNEC patients.Conclusion: In summary, we construct the TF-miRNA-target regulatory network, which is helpful to understand the complex LCNEC regulatory mechanisms. The expression difference of regulated multi-effect target factors ETS1, CD36, and hsa-miR-195 has a significant impact on the survival and development of LCNEC.

MiR-124 synergism with ELAVL3 enhances target gene expression to promote neuronal maturity

SummaryNeuron-enriched microRNAs (miRNAs), miR-9/9* and miR-124 (miR-9/9*-124), direct cell fate switching of human fibroblasts to neurons when ectopically expressed by repressing anti-neurogenic genes. How these miRNAs function after the onset of the transcriptome switch to a neuronal fate remains unclear. Here, we identified direct targets of miRNAs by Argonaute (AGO) HITS-CLIP as reprogramming cells activate the neuronal program and reveal the role of miR-124 that directly promotes the expression of its target genes associated with neuronal development and function. The mode of miR-124 as a positive regulator is determined by a neuron-enriched RNA-binding protein, ELAVL3, that interacts with AGO and binds target transcripts, whereas the non-neuronal ELAVL1 counterpart fails to elevate the miRNA-target gene expression. Although existing literature indicate that miRNA-ELAVL1 interaction can result in either target gene upregulation or downregulation in a context-dependent manner, we specifically identified neuronal ELAVL3 as the driver for miRNA target gene upregulation in neurons. In primary human neurons, repressing miR-124 and ELAVL3 led to the downregulation of genes involved in neuronal function and process outgrowth, and cellular phenotypes of reduced inward currents and neurite outgrowth. Results from our study support the role of miR-124 promoting neuronal function through positive regulation of its target genes.

The integrative knowledge base for miRNA-mRNA expression in colorectal cancer

Abstract“miRNA colorectal cancer” (https://mirna-coadread.omics.si/) is a freely available web application for studying microRNA and mRNA expression and their correlation in colorectal cancer. To the best of our knowledge, “miRNA colorectal cancer” has the largest knowledge base of miRNA-target gene expressions and correlations in colorectal cancer, based on the largest available sample size from the same source of data. Data from high-throughput molecular profiling of 295 colon and rectum adenocarcinoma samples from The Cancer Genome Atlas was analyzed and integrated into our knowledge base. The objective of developing this web application was to help researchers to discover the behavior and role of miRNA-target gene interactions in colorectal cancer. For this purpose, results of differential expression and correlation analyses of miRNA and mRNA data collected in our knowledge base are available through web forms. To validate our knowledge base experimentally, we selected genes FN1, TGFB2, RND3, ZEB1 and ZEB2 and miRNAs hsa-miR-200a/b/c-3p, hsa-miR-141-3p and hsa-miR-429. Both approaches revealed a negative correlation between miRNA hsa-miR-200b/c-3p and its target gene FN1 and between hsa-miR-200a-3p and its target TGFB2, thus supporting the usefulness of the developed knowledge base.