CyTRANSFINDER: a Cytoscape 3.3 plugin for three-component (TF, gene, miRNA) signal transduction pathway construction

A comprehensive identification of potential molecular targets and small drugs candidate for melanoma cancer using bioinformatics and network-based screening approach

Melanoma is the third most common malignant skin tumor and has increased in morbidity and mortality over the previous decade due to its rapid spread into the bloodstream or lymphatic system. This study used integrated bioinformatics and network-based methodologies to reliably identify molecular targets and small molecular medicines that may be more successful for Melanoma diagnosis, prognosis and treatment. The statistical LIMMA approach utilized for bioinformatics analysis in this study found 246 common differentially expressed genes (cDEGs) between case and control samples from two microarray gene-expression datasets (GSE130244 and GSE15605). Protein-protein interaction network study revealed 15 cDEGs (PTK2, STAT1, PNO1, CXCR4, WASL, FN1, RUNX2, SOCS3, ITGA4, GNG2, CDK6, BRAF, AGO2, GTF2H1 and AR) to be critical in the development of melanoma (KGs). According to regulatory network analysis, the most important transcriptional and post-transcriptional regulators of DEGs and hub-DEGs are ten transcription factors and three miRNAs. We discovered the pathogenetic mechanisms of MC by studying DEGs' biological processes, molecular function, cellular components and KEGG pathways. We used molecular docking and dynamics modeling to select the four most expressed genes responsible for melanoma malignancy to identify therapeutic candidates. Then, utilizing the Connectivity Map (CMap) database, we analyzed the top 4-hub-DEGs-guided repurposable drugs. We validated four melanoma cancer drugs (Fisetin, Epicatechin Gallate, 1237586-97-8 and PF 431396) using molecular dynamics simulation with their target proteins. As a result, the results of this study may provide resources to researchers and medical professionals for the wet-lab validation of MC diagnosis, prognosis and treatments.Communicated by Ramaswamy H. Sarma.

Integrating miRNA and full-length transcriptome profiling to elucidate the mechanism of muscle growth in Muscovy ducks reveals key roles for miR-301a-3p/ANKRD1

BACKGROUND

The popularity of Muscovy ducks is attributed not only to their conformation traits but also to their slightly higher content of breast and leg meat, as well as their stronger-tasting meat compared to that of typical domestic ducks. However, there is a lack of comprehensive systematic research on the development of breast muscle in Muscovy ducks. In addition, since the number of skeletal muscle myofibers is established during the embryonic period, this study conducted a full-length transcriptome sequencing and microRNA sequencing of the breast muscle. Muscovy ducks at four developmental stages, namely Embryonic Day 21 (E21), Embryonic Day 27 (E27), Hatching Day (D0), and Post-hatching Day 7 (D7), were used to isolate total RNA for analysis.

RESULTS

A total of 68,161 genes and 472 mature microRNAs were identified. In order to uncover deeper insights into the regulation of mRNA by miRNAs, we conducted an integration of the differentially expressed miRNAs (known as DEMs) with the differentially expressed genes (referred to as DEGs) across various developmental stages. This integration allowed us to make predictions regarding the interactions between miRNAs and mRNA. Through this analysis, we identified a total of 274 DEGs that may serve as potential targets for the 68 DEMs. In the predicted miRNA‒mRNA interaction networks, let-7b, miR-133a-3p, miR-301a-3p, and miR-338-3p were the hub miRNAs. In addition, multiple DEMs also showed predicted target relationships with the DEGs associated with skeletal system development. These identified DEGs and DEMs as well as their predicted interaction networks involved in the regulation of energy homeostasis and muscle development were most likely to play critical roles in facilitating the embryo-to-hatchling transition. A candidate miRNA, miR-301a-3p, exhibited increased expression during the differentiation of satellite cells and was downregulated in the breast muscle tissues of Muscovy ducks at E21 compared to E27. A dual-luciferase reporter assay suggested that the ANKRD1 gene, which encodes a transcription factor, is a direct target of miR-301a-3p.

CONCLUSIONS

miR-301a-3p suppressed the posttranscriptional activity of ANKRD1, which is an activator of satellite cell proliferation, as determined with gain- and loss-of-function experiments. miR-301a-3p functions as an inducer of myogenesis by targeting the ANKRD1 gene in Muscovy ducks. These results provide novel insights into the early developmental process of black Muscovy breast muscles and will improve understanding of the underlying molecular mechanisms.

Hsa_circRNA_0040462: a sensor of cells' response to CAP treatment with double-edged roles on breast cancer malignancy

Cold atmospheric plasma (CAP) represents a novel onco-therapeutic approach that has demonstrated its efficacy in many types of tumors. The efficacy of CAP is dose-dependent that determines the panel of tumors feasible for receiving CAP treatment under a certain parameter configuration. Identifying markers for easy and fast prognosis of tumors' sensitivity in response to CAP exposure is of critical value towards optimized therapeutic outcome, the lack of which has largely limited the translation of CAP into clinics. Circular RNAs represent a novel type of biomarkers for disease diagnosis that is featured by easy detection and stability. Through whole transcriptome sequencing, followed by in vitro validations, computational predictions and preliminary functional studies, we identified hsa_circRNA_0040462 as a sensor of breast cancer cells' response to CAP treatment. Yet we warrant the use of hsa_circRNA_0040462 as an onco-therapeutic target given its double-edged roles on breast cancer progression, i.e., suppressive on the growth and promotive on the migrative ability of triple negative breast cancer cells. Our study for the first time focused on markers prognostic of CAP's efficacy and tumors' sensitivity to CAP treatment under a certain parameter configuration, and reported hsa_circRNA_0040462 as a sensor of cells' response to CAP treatment. Also, the uncovered dual roles of hsa_circRNA_0040462 further advanced our knowledge on the complex yet critical regulatory functionalities of circular RNAs in cancer progression.

miRNAs Potentially Involved in Post Lung Transplant-Obliterative Bronchiolitis: The Role of miR-21-5p

Epigenetic changes, including miRNAs deregulation, have been suggested to play a significant role in development of obliterative bronchiolitis (OB) in transplanted lungs. Many studies have tried to identify ideal candidate miRNAs and the downstream pathways implicated in the bronchiolar fibro-obliterative process. Several candidate miRNAs, previously indicated as possibly being associated with OB, were analyzed by combining the quantitative real time-polymerase chain reaction (qRT-PCR) and in situ hybridization (ISH) of lung tissues of OB affected patients. Disease and OB-lesion-specific expression of miR-21-5p was confirmed and by computational analysis we were able to identify the network of genes most probably associated miR-21-5p in the context of OB fibrogenesis. Among all potentially associated genes, STAT3 had a very high probability score. Immunohistochemistry showed that STAT3/miR-21-5p were co-over expressed in OB lesions, thus, suggesting miR-21-5p could regulate STAT3 expression. However, miR-21-5p inhibition in cultures of bronchiolitis obliterans syndrome (BOS) derived myofibroblasts did not significantly affect STAT3 mRNA and protein expression levels. This study demonstrates the specificity of miR-21-5p over-expression in OB lesions and contributes to existing knowledge on the miR-21-5p downstream pathway. Activation of STAT3 is associated with miR-21-5p upregulation, however, STAT-3 network activation is most likely complex and miR-21-5p is not the sole regulator of STAT3.

'One DB to rule them all'-the RING: a Regulatory INteraction Graph combining TFs, genes/proteins, SNPs, diseases and drugs

In the last decade, genomics data have been largely adopted to sketch, study and better understand the complex mechanisms that underlie biological processes. The amount of publicly available data sources has grown accordingly, and several types of regulatory interactions have been collected and documented in literature. Unfortunately, often these efforts do not follow any data naming/interoperability/formatting standards, resulting in high-quality but often uninteroperable heterogeneous data repositories. To efficiently take advantage of the large amount of available data and integrate these heterogeneous sources of information, we built the RING (Regulatory Interaction Graph), an integrative standardized multilevel database of biological interactions able to provide a comprehensive and unmatched high-level perspective on several phenomena that take place in the regulatory cascade and that researchers can use to easily build regulatory networks around entities of interest.

Developmental Chromatin Restriction of Pro-Growth Gene Networks Acts as an Epigenetic Barrier to Axon Regeneration in Cortical Neurons

Axon regeneration in the central nervous system is prevented in part by a developmental decline in the intrinsic regenerative ability of maturing neurons. This loss of axon growth ability likely reflects widespread changes in gene expression, but the mechanisms that drive this shift remain unclear. Chromatin accessibility has emerged as a key regulatory mechanism in other cellular contexts, raising the possibility that chromatin structure may contribute to the age-dependent loss of regenerative potential. Here we establish an integrated bioinformatic pipeline that combines analysis of developmentally dynamic gene networks with transcription factor regulation and genome-wide maps of chromatin accessibility. When applied to the developing cortex, this pipeline detected overall closure of chromatin in sub-networks of genes associated with axon growth. We next analyzed mature CNS neurons that were supplied with various pro-regenerative transcription factors. Unlike prior results with SOX11 and KLF7, here we found that neither JUN nor an activated form of STAT3 promoted substantial corticospinal tract regeneration. Correspondingly, chromatin accessibility in JUN or STAT3 target genes was substantially lower than in predicted targets of SOX11 and KLF7. Finally, we used the pipeline to predict pioneer factors that could potentially relieve chromatin constraints at growth-associated loci. Overall this integrated analysis substantiates the hypothesis that dynamic chromatin accessibility contributes to the developmental decline in axon growth ability and influences the efficacy of pro-regenerative interventions in the adult, while also pointing toward selected pioneer factors as high-priority candidates for future combinatorial experiments. © 2018 Wiley Periodicals, Inc. Develop Neurobiol 00: 000-000, 2018.

KLF6 and STAT3 co-occupy regulatory DNA and functionally synergize to promote axon growth in CNS neurons

The failure of axon regeneration in the CNS limits recovery from damage and disease. Members of the KLF family of transcription factors can exert both positive and negative effects on axon regeneration, but the underlying mechanisms are unclear. Here we show that forced expression of KLF6 promotes axon regeneration by corticospinal tract neurons in the injured spinal cord. RNA sequencing identified 454 genes whose expression changed upon forced KLF6 expression in vitro, including sub-networks that were highly enriched for functions relevant to axon extension including cytoskeleton remodeling, lipid synthesis, and bioenergetics. In addition, promoter analysis predicted a functional interaction between KLF6 and a second transcription factor, STAT3, and genome-wide footprinting using ATAC-Seq data confirmed frequent co-occupancy. Co-expression of the two factors yielded a synergistic elevation of neurite growth in vitro. These data clarify the transcriptional control of axon growth and point the way toward novel interventions to promote CNS regeneration.

Reconstructing context-specific gene regulatory network and identifying modules and network rewiring through data integration

Reconstructing context-specific transcriptional regulatory network is crucial for deciphering principles of regulatory mechanisms underlying various conditions. Recently studies that reconstructed transcriptional networks have focused on individual organisms or cell types and relied on data repositories of context-free regulatory relationships. Here we present a comprehensive framework to systematically derive putative regulator-target pairs in any given context by integrating context-specific transcriptional profiling and public data repositories of gene regulatory networks. Moreover, our framework can identify core regulatory modules and signature genes underlying global regulatory circuitry, and detect network rewiring and core rewired modules in different contexts by considering gene modules and edge (gene interaction) modules collaboratively. We applied our methods to analyzing Autism RNA-seq experiment data and produced biologically meaningful results. In particular, all 11 hub genes in a predicted rewired autistic regulatory subnetwork have been linked to autism based on literature review. The predicted rewired autistic regulatory network may shed some new insight into disease mechanism.

Visualization and Analysis of MicroRNAs within KEGG Pathways using VANESA

MicroRNAs (miRNAs) are small RNA molecules which are known to take part in post-transcriptional regulation of gene expression. Here, VANESA, an existing platform for reconstructing, visualizing, and analysis of large biological networks, has been further expanded to include all experimentally validated human miRNAs available within miRBase, TarBase and miRTarBase. This is done by integrating a custom hybrid miRNA database to DAWIS-M.D., VANESA’s main data source, enabling the visualization and analysis of miRNAs within large biological pathways such as those found within the Kyoto Encyclopedia of Genes and Genomes (KEGG). Interestingly, 99.15 % of human KEGG pathways either contain genes which are targeted by miRNAs or harbor them. This is mainly due to the high number of interaction partners that each miRNA could have (e.g.: hsa-miR-335-5p targets 2544 genes and 71 miRNAs target NUFIP2). We demonstrate the usability of our system by analyzing the measles virus KEGG pathway as a proof-of-principle model and further highlight the importance of integrating miRNAs (both experimentally validated and predicted) into biological networks for the elucidation of novel miRNA-mRNA interactions of biological importance.

Identification of miRNAs Potentially Involved in Bronchiolitis Obliterans Syndrome: A Computational Study

The pathogenesis of Bronchiolitis Obliterans Syndrome (BOS), the main clinical phenotype of chronic lung allograft dysfunction, is poorly understood. Recent studies suggest that epigenetic regulation of microRNAs might play a role in its development. In this paper we present the application of a complex computational pipeline to perform enrichment analysis of miRNAs in pathways applied to the study of BOS. The analysis considered the full set of miRNAs annotated in miRBase (version 21), and applied a sequence of filtering approaches and statistical analyses to reduce this set and to score the candidate miRNAs according to their potential involvement in BOS development. Dysregulation of two of the selected candidate miRNAs–miR-34a and miR-21 –was clearly shown in in-situ hybridization (ISH) on five explanted human BOS lungs and on a rat model of acute and chronic lung rejection, thus definitely identifying miR-34a and miR-21 as pathogenic factors in BOS and confirming the effectiveness of the computational pipeline.