MicroRNAs preferentially target the genes with high transcriptional regulation complexity

Biomarkers in Peripartum Cardiomyopathy-What We Know and What Is Still to Be Found

Peripartum cardiomyopathy (PPCM) is a form of heart failure, often severe, that occurs in previously healthy women at the end of their pregnancy or in the first few months after delivery. In PPCM, the recovery of heart function reaches 45-50%. However, the all-cause mortality in long-term observation remains high, reaching 20% irrespective of recovery status. The incidence of PPCM is increasing globally; therefore, effort is required to clarify the pathophysiological background of the disease, as well as to discover specific diagnostic and prognostic biomarkers. The etiology of the disease remains unclear, including oxidative stress; inflammation; hormonal disturbances; endothelial, microcirculatory, cardiomyocyte and extracellular matrix dysfunction; fibrosis; and genetic mutations. Currently, antiangiogenic 16-kDa prolactin (PRL), cleaved from standard 23-kDa PRL in the case of unbalanced oxidative stress, is recognized as the main trigger of the disease. In addition, 16-kDa PRL causes damage to cardiomyocytes, acting via microRNA-146a secreted from endothelial cells as a cause of the NF-κβ pathway. Bromocriptine, which inhibits the secretion of PRL from the pituitary gland, is now the only specific treatment for PPCM. Many different phenotypes of the disease, as well as cases of non-responders to bromocriptine treatment, indicate other pathophysiological pathways that need further investigation. Biomarkers in PPCM are not well established. There is a deficiency in specific diagnostic biomarkers. Pro-brain-type natriuretic peptide (BNP) and N-terminal BNP are the best, however unspecific, diagnostic biomarkers of heart failure at the moment. Therefore, more efforts should be engaged in investigating more specific biomolecules of a diagnostic and prognostic manner such as 16-kDa PRL, galectin-3, myeloperoxidase, or soluble Fms-like tyrosine kinase-1/placental growth factor ratio. In this review, we present the current state of knowledge and future directions of exploring PPCM pathophysiology, including microRNA and heat shock proteins, which may improve diagnosis, treatment monitoring, and the development of specific treatment strategies, and consequently improve patients' prognosis and outcome.

Panoramic view of key cross-talks underpinning the oral squamous cell carcinoma stemness - unearthing the future opportunities

The clinical management of oral cancer is often frequented with challenges that arise from relapse, recurrence, invasion and resistance towards the cornerstone chemo and radiation therapies. The recent conceptual advancement in oncology has substantiated the role of cancer stem cells (CSC) as a predominant player of these intricacies. CSC are a sub-group of tumor population with inherent adroitness to self-renew with high plasticity. During tumor evolution, the structural and functional reprogramming persuades the cancer cells to acquire stem-cell like properties, thus presenting them with higher survival abilities and treatment resistance. An appraisal on key features that govern the stemness is of prime importance to confront the current challenges encountered in oral cancer. The nurturing niche of CSC for maintaining its stemness characteristics is thought to be modulated by complex multi-layered components encompassing neoplastic cells, extracellular matrix, acellular components, circulatory vessels, various cascading signaling molecules and stromal cells. This review focuses on recapitulating both intrinsic and extrinsic mechanisms that impart the stemness. There are contemplating evidences that demonstrate the role of transcription factors (TF) in sustaining the neoplastic stem cell's pluripotency and plasticity alongside the miRNA in regulation of crucial genes involved in the transformation of normal oral mucosa to malignancy. This review illustrates the interplay between miRNA and various known TF of oral cancer such as c-Myc, SOX, STAT, NANOG and OCT in orchestrating the stemness and resistance features. Further, the cross-talks involved in tumor micro-environment inclusive of cytokines, macrophages, extra cellular matrix, angiogenesis leading pathways and influential factors of hypoxia on tumorigenesis and CSC survival have been elucidated. Finally, external factorial influence of oral microbiome gained due to the dysbiosis is also emphasized. There are growing confirmations of the possible roles of microbiomes in the progression of oral cancer. Given this, an attempt has been made to explore the potential links including EMT and signaling pathways towards resistance and stemness. This review provides a spectrum of understanding on stemness and progression of oral cancers at various regulatory levels along with their current therapeutic knowledge. These mechanisms could be exploited for future research to expand potential treatment strategies.

Whole transcriptome analysis resulted in the identification of Chinese sprangletop (Leptochloa chinensis) genes involved in cyhalofop-butyl tolerance

Background

Chinese sprangletop [Leptochloa chinensis (L.) Nees] is an annual malignant weed, which can often be found in paddy fields. Cyhalofop-butyl is a specialized herbicide which is utilized to control L. chinensis. However, in many areas, L. chinensis has become tolerant to this key herbicide due to its continuous long-term use.

Results

In this study, we utilized a tolerant (LC18002) and a sensitive (LC17041) L. chinensis populations previously identified in our laboratory, which were divided into four different groups. We then employed whole transcriptome analysis to identify candidate genes which may be involved in cyhalofop-butyl tolerance. This analysis resulted in the identification of six possible candidate genes, including three cytochrome P450 genes and three ATP-binding cassette transporter genes. We then carried out a phylogenetic analysis to identify homologs of the differentially expressed cytochrome P450 genes. This phylogenetic analysis indicated that all genes have close homologs in other species, some of which have been implicated in non-target site resistance (NTSR).

Conclusions

This study is the first to use whole transcriptome analysis to identify herbicide non-target resistance genes in L. chinensis. The differentially expressed genes represent promising targets for better understanding herbicide tolerance in L. chinensis. The six genes belonging to classes already associated in herbicide tolerance may play important roles in the metabolic resistance of L. chinensis to cyhalofop-butyl, although the exact mechanisms require further study.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07856-z.

Telomerase reverse transcriptase activates transcription of miR500A to inhibit Hedgehog signalling and promote cell invasiveness

Telomerase reverse transcriptase (TERT) maintains telomere homeostasis, thus ensuring chromosome stability and cell proliferation. In addition, several telomere-independent functions of human TERT have been described. In this study, we report that TERT binds directly to the TCF binding elements located upstream of the oncomiR miR500A, and induces its transcription. This function was independent of the telomerase activity, as shown with experiments using catalytically inactive TERT and inhibitors of TERT and the TERT RNA component. miR500A was in turn found to target three key components of the Hedgehog signalling pathway: Patched 1; Gli family zinc finger 3; and Cullin 3, thereby promoting cell invasion. Our results point to the crucial role of the TERT-miR500A-Hedgehog axis in tumour aggressiveness and highlight the therapeutic potential of targeting noncanonical TERT functions in cancer.

Analysis of Cells Proliferation and MicroRNAs Expression Profile in Human Chondrosarcoma SW1353 Cells Exposed to Iodine-125 Seeds Irradiation

Chondrosarcoma is the second most common bone malignancy in adults, and it is often resistant to traditional chemotherapy and radiation therapy. Permanent implantation of iodine-125 (¹²⁵I) seeds has been explored for the treatment of many types of cancer. In this study, the aim was to investigate the proliferative and microRNA (miRNA) effects of ¹²⁵I seeds irradiation on human chondrosarcoma SW1353 cells. First, a new in vitro ¹²⁵I seed irradiation model was established, and cell viability and miRNA microarray assays were performed before and after exposure to the ¹²⁵I seeds. Cell proliferation was inhibited, and miRNA expression was substantially altered by irradiation exposure. The inhibition of cell proliferation was positively correlated with increased radiation doses, with cells showing the highest total radiation dose 7 days after irradiation. A total of 2549 miRNAs were detected in the SW1353 cells after exposure to 6 Gy of radiation, which included 189 differentially expressed miRNAs (98 upregulated and 91 downregulated). Four miRNAs were found to play important roles in the inhibition of cell proliferation after irradiation exposure, including miR-1224-5p, miR-492, miR-135b-5p, and miR-6839-5p. The target genes of the associated miRNAs mentioned were vascular endothelial growth factor A (VEGFA), C-X-C motif chemokine 12 (CXCL12), mitogen-activated protein kinase kinase kinase kinase 3 (MAP4K3), and apoptosis facilitator Bcl-2-like protein 14 (BCL2L14). Hence, the mitogen-activated protein kinase signaling pathway may be involved in how chondrosarcoma cells respond to ¹²⁵I seed irradiation.

MicroRNAs as Emerging Regulators of Signaling in the Tumor Microenvironment

A myriad of signaling molecules in a heuristic network of the tumor microenvironment (TME) pose a challenge and an opportunity for novel therapeutic target identification in human cancers. MicroRNAs (miRs), due to their ability to affect signaling pathways at various levels, take a prominent space in the quest of novel cancer therapeutics. The role of miRs in cancer initiation, progression, as well as in chemoresistance, is being increasingly investigated. The canonical function of miRs is to target mRNAs for post-transcriptional gene silencing, which has a great implication in first-order regulation of signaling pathways. However, several reports suggest that miRs also perform non-canonical functions, partly due to their characteristic non-coding small RNA nature. Examples emerge when they act as ligands for toll-like receptors or perform second-order functions, e.g., to regulate protein translation and interactions. This review is a compendium of recent advancements in understanding the role of miRs in cancer signaling and focuses on the role of miRs as novel regulators of the signaling pathway in the TME.

Co-regulatory Network of Oncosuppressor miRNAs and Transcription Factors for Pathology of Human Hepatic Cancer Stem Cells (HCSC)

Hepatic cancer stem cells (HCSCs) are considered as main players for the hepatocellular carcinoma (HCC) initiation, metastasis, drug resistance and recurrence. There is a growing evidence supporting the down-regulated miRNAs in HCSCs as key suppressors for the stemness traits, but still more details are vague about how these miRNAs modulate the HCC development. To uncover some of these miRNA regulatory aspects in HCSC, we compiled 15 down-regulated miRNA and their validated and predicted up-regulated targets in HCSC. The targets were enriched for several cancer cell stemness hallmarks and CSC pre-metastatic niche, which support these miRNAs role in suppression of HCSCs neoplastic transformation. Further, we constructed miRNA-Transcription factor (TF) regulatory networks, which provided new insights on the role of the proposed miRNA-TF co-regulation in the cancer stemness axis and its cross talk with the surrounding microenvironment. Our analysis revealed HCSC important hubs as candidate regulators for targeting hepatic cancer stemness such as, miR-148a, miR-214, E2F family, MYC and SLC7A5. Finally, we proposed a possible model for miRNA and TF co-regulation of HCSC signaling pathways. Our study identified an HCSC signature and set bridges between the reported results to give guide for future validation of HCC therapeutic strategies avoiding drug resistance.

A Small RNA-Mediated Regulatory Network in Arabidopsis thaliana Demonstrates Connectivity Between phasiRNA Regulatory Modules and Extensive Co-Regulation of Transcription by miRNAs and phasiRNAs

Gene regulation involves the orchestrated action of multiple regulators to fine-tune the expression of genes. Hierarchical interactions and co-regulation among regulators are commonly observed in biological systems, leading to complex regulatory networks. Small RNA (sRNAs) have been shown to be important regulators of gene expression due to their involvement in multiple cellular processes. In plants, microRNA (miRNAs) and phased small interfering RNAs (phasiRNAs) correspond to two well-characterized types of sRNAs involved in the regulation of posttranscriptional gene expression, although information about their targets and interactions with other gene expression regulators is limited. We describe an extended sRNA-mediated regulatory network in Arabidopsis thaliana that provides a reference frame to understand sRNA biogenesis and activity at the genome-wide level. This regulatory network combines a comprehensive evaluation of phasiRNA production and sRNA targets supported by degradome data. The network includes ~17% of genes in the A. thaliana genome, representing ~50% annotated gene ontology (GO) functional categories. Approximately 14% of genes with GO annotations corresponding to regulation of gene expression were found to be under sRNA control. The unbiased bioinformatic approach used to produce the network was able to detect 107 PHAS loci (regions of phasiRNA production), 5,047 active phasiRNAs (~70% of which were non-canonical), and reconstruct 17 regulatory modules resulting from complex regulatory interactions between different sRNA-regulatory pathways. Known regulatory modules like miR173-TAS-PPR/TPR and miR390-TAS3-ARF/F-box were faithfully reconstructed and expanded, illustrating the accuracy and sensitivity of the methods and providing confidence for the validity of findings of previously unrecognized modules. The network presented here includes a 2X increase in the number of identified PHAS loci, a large complement (~70%) of non-canonical phasiRNAs, and the most comprehensive evaluation of sRNA cleavage activity in A. thaliana to date. Structural analysis showed similarities to networks of other biological systems and demonstrated connectivity between phasiRNA regulatory modules with extensive co-regulation of transcripts by miRNAs and phasiRNAs. The described regulatory network provides a reference that will facilitate global analyses of individual plant regulatory programs such as those that control homeostasis, development, and responses to biotic and abiotic environmental changes.

Establishment of microRNA, transcript and protein regulatory networks in Alport syndrome induced pluripotent stem cells

Alport syndrome (AS) is an inherited progressive disease caused by mutations in genes encoding for the α3, α4 and α5 chains, which are an essential component of type IV collagen and are required for formation of the glomerular basement membrane. However, the underlying etiology of AS remains largely unknown, and the aim of the present study was to examine the genetic mechanisms in AS. Induced pluripotent stem cells (iPSCs) were generated from renal tubular cells. The Illumina HiSeq™ 2000 system and iTRAQ‑coupled 2D liquid chromatography‑tandem mass spectrometry were used to generate the sequences of microRNAs (miRNAs), transcripts and proteins from AS‑iPSCs. Integration of miRNA, transcript and protein expression data was used to construct regulatory networks and identify specific miRNA targets amongst the transcripts and proteins. Relative quantitative proteomics using iTRAQ technology revealed 383 differentially abundant proteins, and high‑throughput sequencing identified 155 differentially expressed miRNAs and 1,168 differentially expressed transcripts. Potential miRNA targets were predicted using miRanda, TargetScan and Pictar. All target proteins and transcripts were subjected to network analysis with miRNAs. Gene ontology analysis of the miRNAs and their targets revealed functional information on the iPSCs, including biological process and cell signaling. Kyoto Encyclopedia of Genes and Genomes pathways analysis revealed that the transcripts and proteins were primarily enriched in metabolic and cell adhesion molecule pathways. In addition, the network maps identified hsa‑miRNA (miR)‑4775 as a prominent miRNA that was associated with a number of targets. Similarly, the prominent ELV‑like protein 1‑A and epidermal growth factor receptor (EGFR)‑associated transcripts were identified. Reverse transcription‑quantitative polymerase chain reaction analysis was used to confirm the upregulation of hsa‑miR‑4775 and EGFR. The integrated approach used in the present study provided a comprehensive molecular characterization of AS. The results may also further understanding of the genetic pathogenesis of AS and facilitate the identification of candidate biomarkers for AS.

In silico identification and characterization of a diverse subset of conserved microRNAs in bioenergy crop Arundo donax L

MicroRNAs (miRNAs) are small non-coding RNA molecules involved in the post-transcriptional regulation of gene expression in plants. Arundo donax L. is a perennial C₃ grass considered one of the most promising bioenergy crops. Despite its relevance, many fundamental aspects of its biology still remain to be elucidated. In the present study we carried out the first in silico mining and tissue-specific characterization of microRNAs and their putative targets in A. donax. We identified a total of 141 miRNAs belonging to 14 families along with the corresponding primary miRNAs, precursor miRNAs and a total of 462 high-confidence predicted targets and novel target sites were validated by 5′-race. Gene Ontology functional annotation showed that miRNA targets are constituted mainly by transcription factors, but three of the newly validated targets are enzymes involved in novel functions like RNA editing, acyl lipid metabolism and post-Golgi trafficking. Folding variability of pre-miRNA loops and phylogenetic analyses indicate variable selective pressure acting on the different miRNA families. The set of miRNAs identified in this study will pave the road to further miRNA research in Arundo donax and contribute towards a better understanding of miRNA-mediated gene regulatory processes in other bioenergy crops.

MicroRNAs in pancreatic cancer diagnosis and therapy

Pancreatic cancer remains a disease with very poor prognosis (only 5-6% of patients are still alive after five years). Attempts to improve the results of treatment of pancreatic cancer focus on a better understanding of the pathogenesis, and non-invasive diagnostic methods (genetic testing from peripheral blood), which would create the possibility of early diagnosis and early surgical treatment before the onset of metastasis. New hopes for the improvement of early diagnosis and treatment of pancreatic ductal adenocarcinoma (PDAC) are associated with genetic testing of microRNA expression changes. A large body of evidence has revealed that microRNAs are aberrantly expressed in the serum and in cancer tissues and elicit oncogenic or tumour-suppressive functions. Selected microRNAs can distinguish pancreatic ductal adenocarcinoma from non-cancerous lesions of the pancreas. This review focuses on the involvement of microRNAs in the early diagnosis of pancreatic cancer. Research results related to the development of a novel therapeutic strategy based on the modulation of microRNA expressions for a better outcome in patients with pancreatic cancer are also presented.

Transcriptomic signature of drought response in pearl millet (Pennisetum glaucum (L.) and development of web-genomic resources

Pearl millet, (Pennisetum glaucum L.), an efficient (C4) crop of arid/semi-arid regions is known for hardiness. Crop is valuable for bio-fortification combating malnutrition and diabetes, higher caloric value and wider climatic resilience. Limited studies are done in pot-based experiments for drought response at gene-expression level, but field-based experiment mimicking drought by withdrawal of irrigation is still warranted. We report de novo assembly-based transcriptomic signature of drought response induced by irrigation withdrawal in pearl millet. We found 19983 differentially expressed genes, 7595 transcription factors, gene regulatory network having 45 hub genes controlling drought response. We report 34652 putative markers (4192 simple sequence repeats, 12111 SNPs and 6249 InDels). Study reveals role of purine and tryptophan metabolism in ABA accumulation mediating abiotic response in which MAPK acts as major intracellular signal sensing drought. Results were validated by qPCR of 13 randomly selected genes. We report the first web-based genomic resource ( http://webtom.cabgrid.res.in/pmdtdb/ ) which can be used for candidate genes-based SNP discovery programs and trait-based association studies. Looking at climatic change, nutritional and pharmaceutical importance of this crop, present investigation has immense value in understanding drought response in field condition. This is important in germplasm management and improvement in endeavour of pearl millet productivity.

Integration analysis of microRNA and mRNA paired expression profiling identifies deregulated microRNA-transcription factor-gene regulatory networks in ovarian endometriosis

BACKGROUND

The etiology and pathophysiology of endometriosis remain unclear. Accumulating evidence suggests that aberrant microRNA (miRNA) and transcription factor (TF) expression may be involved in the pathogenesis and development of endometriosis. This study therefore aims to survey the key miRNAs, TFs and genes and further understand the mechanism of endometriosis.

METHODS

Paired expression profiling of miRNA and mRNA in ectopic endometria compared with eutopic endometria were determined by high-throughput sequencing techniques in eight patients with ovarian endometriosis. Binary interactions and circuits among the miRNAs, TFs, and corresponding genes were identified by the Pearson correlation coefficients. miRNA-TF-gene regulatory networks were constructed using bioinformatic methods. Eleven selected miRNAs and TFs were validated by quantitative reverse transcription-polymerase chain reaction in 22 patients.

RESULTS

Overall, 107 differentially expressed miRNAs and 6112 differentially expressed mRNAs were identified by comparing the sequencing of the ectopic endometrium group and the eutopic endometrium group. The miRNA-TF-gene regulatory network consists of 22 miRNAs, 12 TFs and 430 corresponding genes. Specifically, some key regulators from the miR-449 and miR-34b/c cluster, miR-200 family, miR-106a-363 cluster, miR-182/183, FOX family, GATA family, and E2F family as well as CEBPA, SOX9 and HNF4A were suggested to play vital regulatory roles in the pathogenesis of endometriosis.

CONCLUSION

Integration analysis of the miRNA and mRNA expression profiles presents a unique insight into the regulatory network of this enigmatic disorder and possibly provides clues regarding replacement therapy for endometriosis.

Integrated analysis of the potential roles of miRNA‑mRNA networks in triple negative breast cancer

Triple negative breast cancer (TNBC) is a type of breast cancer where the tumor cells are negative for the estrogen, progesterone and human epidermal growth factor 2 receptors. To date, expression profiling of microRNA (miRNA/miR) and mRNA sequences have been widely applied for the diagnosis of TNBC. In the present study, an integrated analysis of miRNA‑mRNA profiling arrays was performed. A total of five dysregulated miRNAs in patients with TNBC were identified, including upregulated miR‑558 expression and downregulated miR‑320d‑1, miR‑548v, miR‑99a and miR‑21 expression. In addition, 49 potential target mRNA sequences were identified. Bioinformatics analyses were performed on the identified miRNAs and mRNAs, including gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes pathway and miRNA‑mRNA network analyses. A total of 31 GO terms and three signaling pathways were identified. The results indicated that the differentially expressed miRNAs and their potential target mRNAs may affect the pathogenesis of TNBC, and may therefore be considered as promising biomarkers for the early diagnosis and targeted therapy of patients with TNBC.

Discovery of Putative Herbicide Resistance Genes and Its Regulatory Network in Chickpea Using Transcriptome Sequencing

Background: Chickpea (Cicer arietinum L.) contributes 75% of total pulse production. Being cheaper than animal protein, makes it important in dietary requirement of developing countries. Weed not only competes with chickpea resulting into drastic yield reduction but also creates problem of harboring fungi, bacterial diseases and insect pests. Chemical approach having new herbicide discovery has constraint of limited lead molecule options, statutory regulations and environmental clearance. Through genetic approach, transgenic herbicide tolerant crop has given successful result but led to serious concern over ecological safety thus non-transgenic approach like marker assisted selection is desirable. Since large variability in tolerance limit of herbicide already exists in chickpea varieties, thus the genes offering herbicide tolerance can be introgressed in variety improvement programme. Transcriptome studies can discover such associated key genes with herbicide tolerance in chickpea. Results: This is first transcriptomic studies of chickpea or even any legume crop using two herbicide susceptible and tolerant genotypes exposed to imidazoline (Imazethapyr). Approximately 90 million paired-end reads generated from four samples were processed and assembled into 30,803 contigs using reference based assembly. We report 6,310 differentially expressed genes (DEGs), of which 3,037 were regulated by 980 miRNAs, 1,528 transcription factors associated with 897 DEGs, 47 Hub proteins, 3,540 putative Simple Sequence Repeat-Functional Domain Marker (SSR-FDM), 13,778 genic Single Nucleotide Polymorphism (SNP) putative markers and 1,174 Indels. Randomly selected 20 DEGs were validated using qPCR. Pathway analysis suggested that xenobiotic degradation related gene, glutathione S-transferase (GST) were only up-regulated in presence of herbicide. Down-regulation of DNA replication genes and up-regulation of abscisic acid pathway genes were observed. Study further reveals the role of cytochrome P450, xyloglucan endotransglucosylase/hydrolase, glutamate dehydrogenase, methyl crotonoyl carboxylase and of thaumatin-like genes in herbicide resistance. Conclusion: Reported DEGs can be used as genomic resource for future discovery of candidate genes associated with herbicide tolerance. Reported markers can be used for future association studies in order to develop marker assisted selection (MAS) for refinement. In endeavor of chickpea variety development programme, these findings can be of immense use in improving productivity of chickpea germplasm.

The Binding Sites of miR-619-5p in the mRNAs of Human and Orthologous Genes

Background

Normally, one miRNA interacts with the mRNA of one gene. However, there are miRNAs that can bind to many mRNAs, and one mRNA can be the target of many miRNAs. This significantly complicates the study of the properties of miRNAs and their diagnostic and medical applications.

Results

The search of 2,750 human microRNAs (miRNAs) binding sites in 12,175 mRNAs of human genes using the MirTarget program has been completed. For the binding sites of the miR-619-5p the hybridization free energy of the bonds was equal to 100% of the maximum potential free energy. The mRNAs of 201 human genes have complete complementary binding sites of miR-619-5p in the 3’UTR (214 sites), CDS (3 sites), and 5’UTR (4 sites). The mRNAs of CATAD1, ICA1L, GK5, POLH, and PRR11 genes have six miR-619-5p binding sites, and the mRNAs of OPA3 and CYP20A1 genes have eight and ten binding sites, respectively. All of these miR-619-5p binding sites are located in the 3’UTRs. The miR-619-5p binding site in the 5’UTR of mRNA of human USP29 gene is found in the mRNAs of orthologous genes of primates. Binding sites of miR-619-5p in the coding regions of mRNAs of C8H8orf44, C8orf44, and ISY1 genes encode the WLMPVIP oligopeptide, which is present in the orthologous proteins. Binding sites of miR-619-5p in the mRNAs of transcription factor genes ZNF429 and ZNF429 encode the AHACNP oligopeptide in another reading frame. Binding sites of miR-619-5p in the 3’UTRs of all human target genes are also present in the 3’UTRs of orthologous genes of mammals. The completely complementary binding sites for miR-619-5p are conservative in the orthologous mammalian genes.

Conclusions

The majority of miR-619-5p binding sites are located in the 3’UTRs but some genes have miRNA binding sites in the 5’UTRs of mRNAs. Several genes have binding sites for miRNAs in the CDSs that are read in different open reading frames. Identical nucleotide sequences of binding sites encode different amino acids in different proteins. The binding sites of miR-619-5p in 3’UTRs, 5’UTRs and CDSs are conservative in the orthologous mammalian genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3811-6) contains supplementary material, which is available to authorized users.

High throughput deep degradome sequencing reveals microRNAs and their targets in response to drought stress in mulberry (Morus alba)

MicroRNAs (miRNAs) play important regulatory roles by targeting mRNAs for cleavage or translational repression. Identification of miRNA targets is essential to better understanding the roles of miRNAs. miRNA targets have not been well characterized in mulberry (Morus alba). To anatomize miRNA guided gene regulation under drought stress, transcriptome-wide high throughput degradome sequencing was used in this study to directly detect drought stress responsive miRNA targets in mulberry. A drought library (DL) and a contrast library (CL) were constructed to capture the cleaved mRNAs for sequencing. In CL, 409 target genes of 30 conserved miRNA families and 990 target genes of 199 novel miRNAs were identified. In DL, 373 target genes of 30 conserved miRNA families and 950 target genes of 195 novel miRNAs were identified. Of the conserved miRNA families in DL, mno-miR156, mno-miR172, and mno-miR396 had the highest number of targets with 54, 52 and 41 transcripts, respectively, indicating that these three miRNA families and their target genes might play important functions in response to drought stress in mulberry. Additionally, we found that many of the target genes were transcription factors. By analyzing the miRNA-target molecular network, we found that the DL independent networks consisted of 838 miRNA-mRNA pairs (63.34%). The expression patterns of 11 target genes and 12 correspondent miRNAs were detected using qRT-PCR. Six miRNA targets were further verified by RNA ligase-mediated 5' rapid amplification of cDNA ends (RLM-5' RACE). Gene Ontology (GO) annotations and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that these target transcripts were implicated in a broad range of biological processes and various metabolic pathways. This is the first study to comprehensively characterize target genes and their associated miRNAs in response to drought stress by degradome sequencing in mulberry. This study provides a framework for understanding the molecular mechanisms of drought resistance in mulberry.

Understanding the combinatorial action of transcription factors and microRNA regulation from regions of open chromatin

Transcriptional regulatory cascades are always triggered through the combinatorial interplay between transcription factors (TFs) and microRNAs (miRNAs) in eukaryotes. However, it is still a very substantial undertaking to dynamically profile their coordinated actions. In this work, we compared the differences in TFBS numbers between miRNA targets and non-targets, and found that miRNA targets tend to have more TFBS numbers. With the attempt to comprehensively understand the combinatorial action of TF and miRNA regulation from regions of open chromatin, we retrieved recently published DNase I hypersensitive sites (DHSs) across different human cell lines. The result showed that the differences are more statistically significant in DHS regions than non-DHS regions. Next, we trained classifiers for miRNA targets and non-targets. The result showed that TFBSs located in DHS regions achieved a competitive performance when discriminating miRNA targets and non-targets, whereas the performance of classifiers using TFBSs located in non-DHS regions is close to that of a random classifier. After the DHSs were divided into intergenic, transcription start sites (TSSs) and gene body DHS regions based on their genomic locations, only TFBSs located in TSS DHS regions provided a competitive performance. Our results provide us a clue that the coordinated activity of miRNAs and TFs describing the mechanism of gene expression control should be examined in a dynamic perspective.

Describing a Transcription Factor Dependent Regulation of the MicroRNA Transcriptome

While the transcription regulation of protein coding genes was extensively studied, little is known on how transcription factors are involved in transcription of non-coding RNAs, specifically of microRNAs. Here, we propose a strategy to study the potential role of transcription factor in regulating transcription of microRNAs using publically available data, computational resources and high throughput data. We use the H3K4me3 epigenetic signature to identify microRNA promoters and chromatin immunoprecipitation (ChIP)-sequencing data from the ENCODE project to identify microRNA promoters that are enriched with transcription factor binding sites. By transfecting cells of interest with shRNA targeting a transcription factor of interest and subjecting the cells to microRNA array, we study the effect of this transcription factor on the microRNA transcriptome. As an illustrative example we use our study on the effect of STAT3 on the microRNA transcriptome of chronic lymphocytic leukemia (CLL) cells.

Inferring the determinants of protein evolutionary rates in mammals

Understanding the determinants of protein evolutionary rates is one of the most fundamental evolutionary questions. Previous studies have revealed that many biological variables are tightly associated with protein evolutionary rates in mammals. However, the dominant role of these biological variables and their combinatorial effects to evolutionary rates of mammalian proteins are still less understood. In this work, we derived a quantitative model to correlate protein evolutionary rates with the levels of these variables. The result showed that only a small number of variables are necessary to accurately predict protein evolutionary rates, among which miRNA regulation plays the most important role. Our result suggested that biological variables are extensively interrelated and suffer from hidden redundancies in determining protein evolutionary rates. Various variables should be considered in a natural ensemble to comprehensively assess the determinants of protein evolutionary rate.

A novel framework for inferring condition-specific TF and miRNA co-regulation of protein-protein interactions

Recent studies have shown that transcription factors (TFs) and microRNAs (miRNAs), while independently regulate their downstream targets, collaborate with each other to regulate gene expression. However, their synergistic roles in protein-protein interactions (PPIs) remain mostly unknown. In this paper, we present a novel framework (called CoRePPI) for inferring TF and miRNA co-regulation of PPIs. Particularly, CoRePPI is aimed at discovering the co-regulation specific to a condition of interest, by using heterogeneous data, including miRNA and messenger RNA (mRNA) expression profiles, putative miRNA targets, TF targets and PPIs. CoRePPI firstly finds the network motifs indicating the co-regulation of PPIs by TFs and miRNAs in tumor and normal conditions separately. Then by identifying the differential motifs found in one condition but not in the other, it builds the networks consisting of TFs, miRNAs and their co-regulated PPIs specific to different conditions respectively. To validate CoRePPI, we apply it to the Pan-Cancer dataset which includes the expression profiles of 12 cancer types from TCGA. Through network topology analysis, we found that the tumor and normal CoRePPI networks are scale-free. Furthermore, the results of differential and intersected network analysis between the tumor and normal CoRePPI networks suggest that only a small fraction of the regulatory relationships between TFs and miRNAs are conserved in both conditions but they co-regulate different downstream PPIs in tumor and normal conditions; and in different conditions the majority of the regulatory relationships between TFs and miRNAs are different although they may regulate the same PPIs in their respective conditions. The CoRePPI sub-networks constructed for the three types of cancers (breast cancer, lung cancer and ovarian cancer) are all scale-free, and the intersection of these CoRePPI sub-networks can be utilized as the biomarker CoRePPI sub-network of the three types of cancers. The PPI enrichment analyses of the tumor and normal CoRePPI networks suggest that the co-regulating TFs and miRNAs are significantly associated with the specific biological processes, diseases and pathways. In addition, comparing with the two non-condition-specific approaches, the tumor CoRePPI network is found to have the most enriched cancer-related PPIs. Altogether, the results uncover the combined regulatory patterns of TFs and miRNAs on the PPIs, and may provide new insights for research in cancer-associated TFs and miRNAs.

Cell type-selective disease-association of genes under high regulatory load

We previously showed that disease-linked metabolic genes are often under combinatorial regulation. Using the genome-wide ChIP-Seq binding profiles for 93 transcription factors in nine different cell lines, we show that genes under high regulatory load are significantly enriched for disease-association across cell types. We find that transcription factor load correlates with the enhancer load of the genes and thereby allows the identification of genes under high regulatory load by epigenomic mapping of active enhancers. Identification of the high enhancer load genes across 139 samples from 96 different cell and tissue types reveals a consistent enrichment for disease-associated genes in a cell type-selective manner. The underlying genes are not limited to super-enhancer genes and show several types of disease-association evidence beyond genetic variation (such as biomarkers). Interestingly, the high regulatory load genes are involved in more KEGG pathways than expected by chance, exhibit increased betweenness centrality in the interaction network of liver disease genes, and carry longer 3' UTRs with more microRNA (miRNA) binding sites than genes on average, suggesting a role as hubs integrating signals within regulatory networks. In summary, epigenetic mapping of active enhancers presents a promising and unbiased approach for identification of novel disease genes in a cell type-selective manner.

Determining Associations between Human Diseases and non-coding RNAs with Critical Roles in Network Control

Deciphering the association between life molecules and human diseases is currently an important task in systems biology. Research over the past decade has unveiled that the human genome is almost entirely transcribed, producing a vast number of non-protein-coding RNAs (ncRNAs) with potential regulatory functions. More recent findings suggest that many diseases may not be exclusively linked to mutations in protein-coding genes. The combination of these arguments poses the question of whether ncRNAs that play a critical role in network control are also enriched with disease-associated ncRNAs. To address this question, we mapped the available annotated information of more than 350 human disorders to the largest collection of human ncRNA-protein interactions, which define a bipartite network of almost 93,000 interactions. Using a novel algorithmic-based controllability framework applied to the constructed bipartite network, we found that ncRNAs engaged in critical network control are also statistically linked to human disorders (P-value of P = 9.8 × 10⁻¹⁰⁹). Taken together, these findings suggest that the addition of those genes that encode optimized subsets of ncRNAs engaged in critical control within the pool of candidate genes could aid disease gene prioritization studies.

Use of miRNAs in biofluids as biomarkers in dietary and lifestyle intervention studies

The selection of biomarkers in nutrigenomics needs to reflect subtle changes in homoeostasis representing the relation between nutrition and health, or nutrition and disease. It is believed that noncoding RNAs, such as circulating microRNAs (miRNAs), may represent such a new class of integrative biomarkers. Until now, the most relevant body fluids for miRNA quantification in response to nutrition have not been clearly defined, but recent studies listed in this review indicate that miRNAs from plasma or serum, PBMC and faeces might be relevant biomarkers to quantify the physiological impacts of dietary or lifestyle intervention studies. In addition, a number of recent studies also indicate that miRNAs could permit to monitor the impact of diet on gut microbiota. We also discuss the main preanalytical considerations that are important to take into account before miRNA screening which can affect the reproducibility of the data.

Prunus transcription factors: breeding perspectives

Many plant processes depend on differential gene expression, which is generally controlled by complex proteins called transcription factors (TFs). In peach, 1533 TFs have been identified, accounting for about 5.5% of the 27,852 protein-coding genes. These TFs are the reference for the rest of the Prunus species. TF studies in Prunus have been performed on the gene expression analysis of different agronomic traits, including control of the flowering process, fruit quality, and biotic and abiotic stress resistance. These studies, using quantitative RT-PCR, have mainly been performed in peach, and to a lesser extent in other species, including almond, apricot, black cherry, Fuji cherry, Japanese apricot, plum, and sour and sweet cherry. Other tools have also been used in TF studies, including cDNA-AFLP, LC-ESI-MS, RNA, and DNA blotting or mapping. More recently, new tools assayed include microarray and high-throughput DNA sequencing (DNA-Seq) and RNA sequencing (RNA-Seq). New functional genomics opportunities include genome resequencing and the well-known synteny among Prunus genomes and transcriptomes. These new functional studies should be applied in breeding programs in the development of molecular markers. With the genome sequences available, some strategies that have been used in model systems (such as SNP genotyping assays and genotyping-by-sequencing) may be applicable in the functional analysis of Prunus TFs as well. In addition, the knowledge of the gene functions and position in the peach reference genome of the TFs represents an additional advantage. These facts could greatly facilitate the isolation of genes via QTL (quantitative trait loci) map-based cloning in the different Prunus species, following the association of these TFs with the identified QTLs using the peach reference genome.

TMREC: A Database of Transcription Factor and MiRNA Regulatory Cascades in Human Diseases

Over the past decades, studies have reported that the combinatorial regulation of transcription factors (TFs) and microRNAs (miRNAs) is essential for the appropriate execution of biological events and developmental processes. Dysregulations of these regulators often cause diseases. However, there are no available resources on the regulatory cascades of TFs and miRNAs in the context of human diseases. To fulfill this vacancy, we established the ^TMREC database in this study. First, we integrated curated transcriptional and post-transcriptional regulations to construct the TF and miRNA regulatory network. Next, we identified all linear paths using the Breadth First Search traversal method. Finally, we used known disease-related genes and miRNAs to measure the strength of association between cascades and diseases. Currently, ^TMREC consists of 74,248 cascades and 25,194 cascade clusters, involving in 412 TFs, 266 miRNAs and 545 diseases. With the expanding of experimental support regulation data, we will regularly update the database. ^TMREC aims to help experimental biologists to comprehensively analyse gene expression regulation, to understand the aetiology and to predict novel therapeutic targets.^TMREC is freely available at http://bioinfo.hrbmu.edu.cn/TMREC/.

MiRNA and TF co-regulatory network analysis for the pathology and recurrence of myocardial infarction

Myocardial infarction (MI) is a leading cause of death in the world and many genes are involved in it. Transcription factor (TFs) and microRNAs (miRNAs) are key regulators of gene expression. We hypothesized that miRNAs and TFs might play combinatory regulatory roles in MI. After collecting MI candidate genes and miRNAs from various resources, we constructed a comprehensive MI-specific miRNA-TF co-regulatory network by integrating predicted and experimentally validated TF and miRNA targets. We found some hub nodes (e.g. miR-16 and miR-26) in this network are important regulators, and the network can be severed as a bridge to interpret the associations of previous results, which is shown by the case of miR-29 in this study. We also constructed a regulatory network for MI recurrence and found several important genes (e.g. DAB2, BMP6, miR-320 and miR-103), the abnormal expressions of which may be potential regulatory mechanisms and markers of MI recurrence. At last we proposed a cellular model to discuss major TF and miRNA regulators with signaling pathways in MI. This study provides more details on gene expression regulation and regulators involved in MI progression and recurrence. It also linked up and interpreted many previous results.

Characterization of microRNA expression profiles in Leishmania-infected human phagocytes

Leishmania are intracellular protozoa that influence host immune responses eliciting parasite species-specific pathologies. MicroRNAs (miRNAs) are short single-stranded ribonucleic acids that complement gene transcripts to block protein translation and have been shown to regulate immune system molecular mechanisms. Human monocyte-derived dendritic cells (DC) and macrophages (MP) were infected in vitro with Leishmania major or Leishmania donovani parasites. Small RNAs were isolated from total RNA and sequenced to identify mature miRNAs associated with leishmanial infections. Normalized sequence read count profiles revealed a global downregulation in miRNA expression among host cells following infection. Most identified miRNAs were expressed at higher levels in L. donovani-infected cells relative to L. major-infected cells. Pathway enrichments using in silico-predicted gene targets of differentially expressed miRNAs showed evidence of potentially universal MAP kinase signalling pathway effects. Whereas JAK-STAT and TGF-β signalling pathways were more highly enriched using targets of miRNAs upregulated in L. donovani-infected cells, these data provide evidence in support of a selective influence on host cell miRNA expression and regulation in response to differential Leishmania infections.

Integrative analysis of differentially expressed microRNAs of pulmonary alveolar macrophages from piglets during H1N1 swine influenza A virus infection

H1N1 swine influenza A virus (H1N1 SwIV) is one key subtype of influenza viruses with pandemic potential. MicroRNAs (miRNAs) are endogenous small RNA molecules that regulate gene expression. MiRNAs relevant with H1N1 SwIV have rarely been reported. To understand the biological functions of miRNAs during H1N1 SwIV infection, this study profiled differentially expressed (DE) miRNAs in pulmonary alveolar macrophages from piglets during the H1N1 SwIV infection using a deep sequencing approach, which was validated by quantitative real-time PCR. Compared to control group, 70 and 16 DE miRNAs were respectively identified on post-infection day (PID) 4 and PID 7. 56 DE miRNAs were identified between PID 4 and PID 7. Our results suggest that most host miRNAs are down-regulated to defend the H1N1 SwIV infection during the acute phase of swine influenza whereas their expression levels gradually return to normal during the recovery phase to avoid the occurrence of too severe porcine lung damage. In addition, targets of DE miRNAs were also obtained, for which bioinformatics analyses were performed. Our results would be useful for investigating the functions and regulatory mechanisms of miRNAs in human influenza because pig serves as an excellent animal model to study the pathogenesis of human influenza.

Systems approach for the selection of micro-RNAs as therapeutic biomarkers of anti-EGFR monoclonal antibody treatment in colorectal cancer

miRNA plays an important role in tumourgenesis by regulating expression of oncogenes and tumour suppressors. Thus affects cell proliferation and differentiation, apoptosis, invasion and angiogenesis. miRNAs are potential biomarkers for diagnosis, prognosis and therapies of different forms of cancer. However, relationship between response of cancer patients towards targeted therapy and the resulting modifications of the miRNA transcriptome in the context of pathway regulation is poorly understood. With ever-increasing pathways and miRNA-mRNA interaction databases, freely available mRNA and miRNA expression data in multiple cancer therapy have produced an unprecedented opportunity to decipher the role of miRNAs in early prediction of therapeutic efficacy in diseases. Efficient translation of -omics data and accumulated knowledge to clinical decision-making are of paramount scientific and public health interest. Well-structured translational algorithms are needed to bridge the gap from databases to decisions. Herein, we present a novel SMARTmiR algorithm to prospectively predict the role of miRNA as therapeutic biomarker for an anti-EGFR monoclonal antibody i.e. cetuximab treatment in colorectal cancer.

Identifying direct miRNA-mRNA causal regulatory relationships in heterogeneous data

Discovering the regulatory relationships between microRNAs (miRNAs) and mRNAs is an important problem that interests many biologists and medical researchers. A number of computational methods have been proposed to infer miRNA-mRNA regulatory relationships, and are mostly based on the statistical associations between miRNAs and mRNAs discovered in observational data. The miRNA-mRNA regulatory relationships identified by these methods can be both direct and indirect regulations. However, differentiating direct regulatory relationships from indirect ones is important for biologists in experimental designs. In this paper, we present a causal discovery based framework (called DirectTarget) to infer direct miRNA-mRNA causal regulatory relationships in heterogeneous data, including expression profiles of miRNAs and mRNAs, and miRNA target information. DirectTarget is applied to the Epithelial to Mesenchymal Transition (EMT) datasets. The validation by experimentally confirmed target databases suggests that the proposed method can effectively identify direct miRNA-mRNA regulatory relationships. To explore the upstream regulators of miRNA regulation, we further identify the causal feedforward patterns (CFFPs) of TF-miRNA-mRNA to provide insights into the miRNA regulation in EMT. DirectTarget has the potential to be applied to other datasets to elucidate the direct miRNA-mRNA causal regulatory relationships and to explore the regulatory patterns.

Dynamic gene expression patterns in animal models of early and late heart failure reveal biphasic-bidirectional transcriptional activation of signaling pathways

Altered cardiac gene expression in heart failure (HF) has mostly been identified by single-point analysis of end-stage disease. This may miss earlier changes in gene expression that are transient and/or directionally opposite to those observed later. Myocardial datasets from the largest microarray data repository (Gene Expression Omnibus) yielded six HF studies with time-course data. Differentially expressed transcripts between nonfailing controls, early HF (<3 days after cardiac insult) and late HF (usually >2 wk) were determined, and analysis of KEGG pathways and predicted regulatory control elements performed. We found that gene expression followed varying patterns: Downregulation of metabolic pathways occurred early and was sustained into late-stage HF. In contrast, most signaling pathways undergo a complex biphasic pattern: Calcium signaling, p53, apoptosis, and MAPK pathways displayed a bidirectional response, declining early but rising late. These profiles were compatible with specific microRNA (miRNA) and transcription regulators: Estrogen-related receptor-α and myocyte-enhancer factor-2 binding sites were overrepresented in the promoter regions of downregulated transcripts. Concurrently, there were overrepresented binding sites for E2f and ETS family members (E-Twenty Six, including Gabp, Elf1, and Ets2), serum response and interferon regulated factor in biphasic-bidirectional and late-upregulated transcripts. Binding sites for miRNAs downregulated by HF were more common in upregulated transcripts (e.g., miRNA-22,-133a/b, and -150 in early HF and miRNA-1,-9,-499 in late HF). During the development of HF, gene expression is characterized by dynamic overlapping sets of transcripts controlled by specific interrelated regulatory mechanisms. While metabolic gene classes show early and sustained downregulation in HF, signaling pathways undergo a complex biphasic pattern with early down- and more pronounced late upregulation.

Novel approach for coexpression analysis of E2F1-3 and MYC target genes in chronic myelogenous leukemia

BACKGROUND

Chronic myelogenous leukemia (CML) is characterized by tremendous amount of immature myeloid cells in the blood circulation. E2F1-3 and MYC are important transcription factors that form positive feedback loops by reciprocal regulation in their own transcription processes. Since genes regulated by E2F1-3 or MYC are related to cell proliferation and apoptosis, we wonder if there exists difference in the coexpression patterns of genes regulated concurrently by E2F1-3 and MYC between the normal and the CML states.

RESULTS

We proposed a method to explore the difference in the coexpression patterns of those candidate target genes between the normal and the CML groups. A disease-specific cutoff point for coexpression levels that classified the coexpressed gene pairs into strong and weak coexpression classes was identified. Our developed method effectively identified the coexpression pattern differences from the overall structure. Moreover, we found that genes related to the cell adhesion and angiogenesis properties were more likely to be coexpressed in the normal group when compared to the CML group.

CONCLUSION

Our findings may be helpful in exploring the underlying mechanisms of CML and provide useful information in cancer treatment.

A systems biology approach identified different regulatory networks targeted by KSHV miR-K12-11 in B cells and endothelial cells

Background

Kaposi’s sarcoma associated herpes virus (KSHV) is associated with tumors of endothelial and lymphoid origin. During latent infection, KSHV expresses miR-K12-11, an ortholog of the human tumor gene hsa-miR-155. Both gene products are microRNAs (miRNAs), which are important post-transcriptional regulators that contribute to tissue specific gene expression. Advances in target identification technologies and molecular interaction databases have allowed a systems biology approach to unravel the gene regulatory networks (GRNs) triggered by miR-K12-11 in endothelial and lymphoid cells. Understanding the tissue specific function of miR-K12-11 will help to elucidate underlying mechanisms of KSHV pathogenesis.

Results

Ectopic expression of miR-K12-11 differentially affected gene expression in BJAB cells of lymphoid origin and TIVE cells of endothelial origin. Direct miRNA targeting accounted for a small fraction of the observed transcriptome changes: only 29 genes were identified as putative direct targets of miR-K12-11 in both cell types. However, a number of commonly affected biological pathways, such as carbohydrate metabolism and interferon response related signaling, were revealed by gene ontology analysis. Integration of transcriptome profiling, bioinformatic algorithms, and databases of protein-protein interactome from the ENCODE project identified different nodes of GRNs utilized by miR-K12-11 in a tissue-specific fashion. These effector genes, including cancer associated transcription factors and signaling proteins, amplified the regulatory potential of a single miRNA, from a small set of putative direct targets to a larger set of genes.

Conclusions

This is the first comparative analysis of miRNA-K12-11’s effects in endothelial and B cells, from tissues infected with KSHV in vivo. MiR-K12-11 was able to broadly modulate gene expression in both cell types. Using a systems biology approach, we inferred that miR-K12-11 establishes its GRN by both repressing master TFs and influencing signaling pathways, to counter the host anti-viral response and to promote proliferation and survival of infected cells. The targeted GRNs are more reproducible and informative than target gene identification, and our approach can be applied to other regulatory factors of interest.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-668) contains supplementary material, which is available to authorized users.

The properties of binding sites of miR-619-5p, miR-5095, miR-5096, and miR-5585-3p in the mRNAs of human genes

The binding of 2,578 human miRNAs with the mRNAs of 12,175 human genes was studied. It was established that miR-619-5p, miR-5095, miR-5096, and miR-5585-3p bind with high affinity to the mRNAs of the 1215, 832, 725, and 655 genes, respectively. These unique miRNAs have binding sites in the coding sequences and untranslated regions of mRNAs. The mRNAs of many genes have multiple miR-619-5p, miR-5095, miR-5096, and miR-5585-3p binding sites. Groups of mRNAs in which the ordering of the miR-619-5p, miR-5095, miR-5096, and miR-5585-3p binding sites differ were established. The possible functional and evolutional properties of unique miRNAs are discussed.

Using graphical adaptive lasso approach to construct transcription factor and microRNA's combinatorial regulatory network in breast cancer

Discovering the regulation of cancer-related gene is of great importance in cancer biology. Transcription factors and microRNAs are two kinds of crucial regulators in gene expression, and they compose a combinatorial regulatory network with their target genes. Revealing the structure of this network could improve the authors' understanding of gene regulation, and further explore the molecular pathway in cancer. In this article, the authors propose a novel approach graphical adaptive lasso (GALASSO) to construct the regulatory network in breast cancer. GALASSO use a Gaussian graphical model with adaptive lasso penalties to integrate the sequence information as well as gene expression profiles. The simulation study and the experimental profiles verify the accuracy of the authors' approach. The authors further reveal the structure of the regulatory network, and explore the role of feedforward loops in gene regulation. In addition, the authors discuss the combinatorial regulatory effect between transcription factors and microRNAs, and select miR-155 for detailed analysis of microRNA's role in cancer. The proposed GALASSO approach is an efficient method to construct the combinatorial regulatory network. It also provides a new way to integrate different data sources and could find more applications in meta-analysis problem.

Not miR-ly micromanagers: the functions and regulatory networks of microRNAs in mammalian skin

The microRNA (miRNA) pathway is a widespread mechanism of post-transcriptional gene regulation in eukaryotic cells. In animals, each miRNA species can regulate hundreds of protein-coding genes, resulting in pervasive functions for miRNAs in numerous cellular processes. Since the identification of the first mammalian miRNA, the function of miRNAs in mammals has been a topic of great interest, both because of the versatile roles of miRNAs in biological systems, as well as the clinical potential of these regulatory RNAs. With well-defined cell lineages and the availability of versatile tools for both in vivo and in vitro studies, mammalian skin has emerged as an important system in which to examine miRNAs' functions in adult tissues. In this review, we discuss recent insights into the functions and regulatory networks of miRNAs in mammals, with a specific focus on murine skin development as a model system. We first introduce functional analyses of the miRNA biogenesis pathway in the skin, then highlight the functions of individual miRNAs in skin development, followed by an examination of miRNA roles in skin stress responses. We finish with a discussion of miRNA regulatory networks and emphasize future challenges and emerging technologies that permit the genome-wide study of miRNA functions and regulatory mechanisms in mammalian skin.

Impacts of pretranscriptional DNA methylation, transcriptional transcription factor, and posttranscriptional microRNA regulations on protein evolutionary rate

Gene expression is largely regulated by DNA methylation, transcription factor (TF), and microRNA (miRNA) before, during, and after transcription, respectively. Although the evolutionary effects of TF/miRNA regulations have been widely studied, evolutionary analysis of simultaneously accounting for DNA methylation, TF, and miRNA regulations and whether promoter methylation and gene body (coding regions) methylation have different effects on the rate of gene evolution remain uninvestigated. Here, we compared human–macaque and human–mouse protein evolutionary rates against experimentally determined single base-resolution DNA methylation data, revealing that promoter methylation level is positively correlated with protein evolutionary rates but negatively correlated with TF/miRNA regulations, whereas the opposite was observed for gene body methylation level. Our results showed that the relative importance of these regulatory factors in determining the rate of mammalian protein evolution is as follows: Promoter methylation ≈ miRNA regulation > gene body methylation > TF regulation, and further indicated that promoter methylation and miRNA regulation have a significant dependent effect on protein evolutionary rates. Although the mechanisms underlying cooperation between DNA methylation and TFs/miRNAs in gene regulation remain unclear, our study helps to not only illuminate the impact of these regulatory factors on mammalian protein evolution but also their intricate interaction within gene regulatory networks.

MicroRNA buffering and altered variance of gene expression in response to Salmonella infection

One potential role of miRNAs is to buffer variation in gene expression, although conflicting results have been reported. To investigate the buffering role of miRNAs in response to Salmonella infection in pigs, we sequenced miRNA and mRNA in whole blood from 15 pig samples before and after Salmonella challenge. By analyzing inter-individual variation in gene expression patterns, we found that for moderately and lowly expressed genes, putative miRNA targets showed significantly lower expression variance compared with non-miRNA-targets. Expression variance between highly expressed miRNA targets and non-miRNA-targets was not significantly different. Further, miRNA targets demonstrated significantly reduced variance after challenge whereas non-miRNA-targets did not. RNA binding proteins (RBPs) are significantly enriched among the miRNA targets with dramatically reduced variance of expression after Salmonella challenge. Moreover, we found evidence that targets of young (less-conserved) miRNAs showed lower expression variance compared with targets of old (evolutionarily conserved) miRNAs. These findings point to the importance of a buffering effect of miRNAs for relatively lowly expressed genes, and suggest that the reduced expression variation of RBPs may play an important role in response to Salmonella infection.

Computational developments in microRNA-regulated protein-protein interactions

Protein-protein interaction (PPI) is one of the most important functional components of a living cell. Recently, researchers have been interested in investigating the correlation between PPI and microRNA, which has been found to be a regulator at the post-transcriptional level. Studies on miRNA-regulated PPI networks will not only facilitate an understanding of the fine tuning role that miRNAs play in PPI networks, but will also provide potential candidates for tumor diagnosis. This review describes basic studies on the miRNA-regulated PPI network in the way of bioinformatics which includes constructing a miRNA-target protein network, describing the features of miRNA-regulated PPI networks and overviewing previous findings based on analysing miRNA-regulated PPI network features.

Human-specific microRNA regulation of FOXO1: implications for microRNA recognition element evolution

MicroRNAs (miRNAs) have been established as important negative post-transcriptional regulators for gene expression. Within the past decade, miRNAs targeting transcription factors (TFs) has emerged as an important mechanism for gene expression regulation. Here, we tested the hypothesis that in TF 3'UTRs, human-specific single nucleotide change(s) that create novel miRNA recognition elements (MREs) contribute to species-specific differences in TF expression. From several potential human-specific TF MREs, one candidate, a member of the Forkhead Box O (FOXO) subclass in the Forkhead family known as Forkhead Box O1 (FOXO1; FKHR; NM_002015) was tested further. Human FOXO1 contains two sites predicted to confer miR-183-mediated post-transcriptional regulation: one specific to humans and the other conserved. Utilizing dual luciferase expression reporters, we show that only the human FOXO1 3'UTR contains a functional miR-183 site, not found in chimpanzee or mouse 3'untranslated regions (UTRs). Site-directed mutagenesis supports functionality of the human-specific miR-183 site, but not the conserved miR-183 site. Via overexpression and target site protection assays, we show that human FOXO1 is regulated by miR-183, but mouse FOXO1 is not. Finally, FOXO1-regulated cellular phenotypes, including cell invasion and proliferation, are impacted by miR-183 targeting only in human cells. These results provide strong evidence for human-specific gain of TF MREs, a process that may underlie evolutionary differences between phylogenic groups.

Genome-wide characterization of microRNA in foxtail millet (Setaria italica)

BACKGROUND

MicroRNAs (miRNAs) are a class of short non-coding, endogenous RNAs that play key roles in many biological processes in both animals and plants. Although many miRNAs have been identified in a large number of organisms, the miRNAs in foxtail millet (Setaria italica) have, until now, been poorly understood.

RESULTS

In this study, two replicate small RNA libraries from foxtail millet shoots were sequenced, and 40 million reads representing over 10 million unique sequences were generated. We identified 43 known miRNAs, 172 novel miRNAs and 2 mirtron precursor candidates in foxtail millet. Some miRNA*s of the known and novel miRNAs were detected as well. Further, eight novel miRNAs were validated by stem-loop RT-PCR. Potential targets of the foxtail millet miRNAs were predicted based on our strict criteria. Of the predicted target genes, 79% (351) had functional annotations in InterPro and GO analyses, indicating the targets of the miRNAs were involved in a wide range of regulatory functions and some specific biological processes. A total of 69 pairs of syntenic miRNA precursors that were conserved between foxtail millet and sorghum were found. Additionally, stem-loop RT-PCR was conducted to confirm the tissue-specific expression of some miRNAs in the four tissues identified by deep-sequencing.

CONCLUSIONS

We predicted, for the first time, 215 miRNAs and 447 miRNA targets in foxtail millet at a genome-wide level. The precursors, expression levels, miRNA* sequences, target functions, conservation, and evolution of miRNAs we identified were investigated. Some of the novel foxtail millet miRNAs and miRNA targets were validated experimentally.

miRNA regulation in the context of functional protein networks: principles and applications

MicroRNAs (miRNAs) are small non-coding endogenous regulatory RNAs that fine-tune gene expression in a wide range of biological processes and diseases. miRNAs exert their function by targeting mRNAs to trigger their degradation or inhibit protein translation. The proteins encoded by the genes targeted by miRNAs may act as key components of cellular networks, thus the use of biological molecular network information for the purposes of elucidating the role of miRNAs in molecular disease mechanism is a key objective in systems biomedicine. The crosstalk layer between miRNA-target networks and functional protein is rich sources of information to explore the function of miRNAs at the system level. Characterizing the influence of miRNAs in the context of the target (protein interactors of the target) is in the early stages with potential to help better understand how miRNAs function within the cellular networks. In this article, the latest research on the cross-talk between miRNAs and protein networks, particularly physical protein interactions and gene regulatory networks is summarized. This article also covers recent research on understanding the biology of miRNAs at the system level and defines principles of miRNA regulation of protein and gene regulatory networks. The second part of the article highlights the promise of considering the protein context of the miRNA target when searching for functional miRNA-target interactions. Some of the applications of integrating protein networks with miRNA-targets that have clinical and functional utility are described.