α-ketoglutarate orchestrates macrophage activation through metabolic and epigenetic reprogramming.

Glutamine metabolism provides synergistic support for macrophage activation and elicitation of desirable immune responses; however, the underlying mechanisms regulated by glutamine metabolism to orchestrate macrophage activation remain unclear. Here we show that the production of α-ketoglutarate (αKG) via glutaminolysis is important for alternative (M2) activation of macrophages, including engagement of fatty acid oxidation (FAO) and Jmjd3-dependent epigenetic reprogramming of M2 genes. This M2-promoting mechanism is further modulated by a high αKG/succinate ratio, whereas a low ratio strengthens the proinflammatory phenotype in classically activated (M1) macrophages. As such, αKG contributes to endotoxin tolerance after M1 activation. This study reveals new mechanistic regulations by which glutamine metabolism tailors the immune responses of macrophages through metabolic and epigenetic reprogramming.

A genetic basis of mitochondrial DNAJA3 in nonalcoholic steatohepatitis-related hepatocellular carcinoma.

BACKGROUND AND AIMS: NAFLD is the most common form of liver disease worldwide, but only a subset of individuals with NAFLD may progress to NASH. While NASH is an important etiology of HCC, the underlying mechanisms responsible for the conversion of NAFLD to NASH and then to HCC are poorly understood. We aimed to identify genetic risk genes that drive NASH and NASH-related HCC. APPROACH AND RESULTS: We searched genetic alleles among the 24 most significant alleles associated with body fat distribution from a genome-wide association study of 344,369 individuals and validated the top allele in 3 independent cohorts of American and European patients (N=1380) with NAFLD/NASH/HCC. We identified an rs3747579-TT variant significantly associated with NASH-related HCC and demonstrated that rs3747579 is expression quantitative trait loci of a mitochondrial DnaJ Heat Shock Protein Family (Hsp40) Member A3 ( DNAJA3 ). We also found that rs3747579-TT and a previously identified PNPLA3 as a functional variant of NAFLD to have significant additional interactions with NASH/HCC risk. Patients with HCC with rs3747579-TT had a reduced expression of DNAJA3 and had an unfavorable prognosis. Furthermore, mice with hepatocyte-specific Dnaja3 depletion developed NASH-dependent HCC either spontaneously under a normal diet or enhanced by diethylnitrosamine. Dnaja3 -deficient mice developed NASH/HCC characterized by significant mitochondrial dysfunction, which was accompanied by excessive lipid accumulation and inflammatory responses. The molecular features of NASH/HCC in the Dnaja3 -deficient mice were closely associated with human NASH/HCC. CONCLUSIONS: We uncovered a genetic basis of DNAJA3 as a key player of NASH-related HCC.

MED1 Regulates BMP/TGF-ß in Endothelium: Implication for Pulmonary Hypertension.

BACKGROUND: Dysregulated BMP (bone morphogenetic protein) or TGF-ß (transforming growth factor beta) signaling pathways are imperative in idiopathic and familial pulmonary arterial hypertension (PAH) as well as experimental pulmonary hypertension (PH) in rodent models. MED1 (mediator complex subunit 1) is a key transcriptional co-activator and KLF4 (Krüppel-like factor 4) is a master transcription factor in endothelium. However, MED1 and KLF4 epigenetic and transcriptional regulations of the BMP/TGF-ß axes in pulmonary endothelium and their dysregulations leading to PAH remain elusive. We investigate the MED1/KLF4 co-regulation of the BMP/TGF-ß axes in endothelium by studying the epigenetic regulation of BMPR2 (BMP receptor type II), ETS-related gene (ERG), and TGFBR2 (TGF-ß receptor 2) and their involvement in the PH. METHODS: High-throughput screening involving data from RNA-seq, MED1 ChIP-seq, H3K27ac ChIP-seq, ATAC-seq, and high-throughput chromosome conformation capture together with in silico computations were used to explore the epigenetic and transcriptional regulation of BMPR2, ERG, and TGFBR2 by MED1 and KLF4. In vitro experiments with cultured pulmonary arterial endothelial cells (ECs) and bulk assays were used to validate results from these in silico analyses. Lung tissue from patients with idiopathic PAH, animals with experimental PH, and mice with endothelial ablation of MED1 (EC-MED1-/-) were used to study the PH-protective effect of MED1. RESULTS: Levels of MED1 were decreased in lung tissue or pulmonary arterial endothelial cells from idiopathic PAH patients and rodent PH models. Mechanistically, MED1 acted synergistically with KLF4 to transactivate BMPR2, ERG, and TGFBR2 via chromatin remodeling and enhancer-promoter interactions. EC-MED1-/- mice showed PH susceptibility. In contrast, MED1 overexpression mitigated the PH phenotype in rodents. CONCLUSIONS: A homeostatic regulation of BMPR2, ERG, and TGFBR2 in ECs by MED1 synergistic with KLF4 is essential for the normal function of the pulmonary endothelium. Dysregulation of MED1 and the resulting impairment of the BMP/TGF-ß signaling is implicated in the disease progression of PAH in humans and PH in rodent models.

Circulating microRNAs Associated With Reversible Cerebral Vasoconstriction Syndrome.

OBJECTIVE: The purpose of this study was to investigate the significance of circulating micro RNAs (miRNAs) in the pathogenesis of reversible cerebral vasoconstriction syndrome (RCVS). METHODS: We prospectively recruited 3 independent cohorts of patients with RCVS and age-matched and sex-matched controls in a single medical center. Next-generation small RNA sequencing followed by quantitative polymerase chain reaction (PCR) was used to identify and validate differentially expressed miRNAs, which was cross-validated in migraine patients in ictal stage or interictal stage. Computational analysis was used to predict the target genes of miRNAs, followed by in vitro functional analysis. RESULTS: We identified a panel of miRNAs including miR-130a-3p, miR-130b-3p, let-7a-5p, let-7b-5p, and let-7f-5p that well differentiated patients with RCVS from controls (area under the receiver operating characteristics curve [AUC] was 0.906, 0.890, and 0.867 in the 3 cohorts, respectively). The abundance of let-7a-5p, let-7b-5p, and let-7f-5p, but not miR-130a-3p nor miR-130b-3p, was significantly higher in patients with ictal migraine compared with that of controls and patients with interictal migraine. Target prediction and pathway enrichment analysis suggested that the transforming growth factor-ß signaling pathway and endothelin-1 responsible for vasomotor control might link these miRNAs to RCVS pathogenesis, which was confirmed in vitro by transfecting miRNAs mimics or incubating the patients' cerebrospinal fluid (CSF) in 3 different vascular endothelial cells. Moreover, miR-130a-3p was associated with imaging-proven disruption of the blood-brain barrier (BBB) in patients with RCVS and its overexpression led to reduced transendothelial electrical resistance (ie, increased permeability) in in vitro human BBB model. INTERPRETATION: We identified the circulating miRNA signatures associated with RCVS, which may be functionally linked to its headache, BBB integrity, and vasomotor function. ANN NEUROL 2021;89:459-473.

miRTarBase 2020: updates to the experimentally validated microRNA-target interaction database.

MicroRNAs (miRNAs) are small non-coding RNAs (typically consisting of 18-25 nucleotides) that negatively control expression of target genes at the post-transcriptional level. Owing to the biological significance of miRNAs, miRTarBase was developed to provide comprehensive information on experimentally validated miRNA-target interactions (MTIs). To date, the database has accumulated >13,404 validated MTIs from 11,021 articles from manual curations. In this update, a text-mining system was incorporated to enhance the recognition of MTI-related articles by adopting a scoring system. In addition, a variety of biological databases were integrated to provide information on the regulatory network of miRNAs and its expression in blood. Not only targets of miRNAs but also regulators of miRNAs are provided to users for investigating the up- and downstream regulations of miRNAs. Moreover, the number of MTIs with high-throughput experimental evidence increased remarkably (validated by CLIP-seq technology). In conclusion, these improvements promote the miRTarBase as one of the most comprehensively annotated and experimentally validated miRNA-target interaction databases. The updated version of miRTarBase is now available at http://miRTarBase.cuhk.edu.cn/.

Multi-omics profiling reveals microRNA-mediated insulin signaling networks.

BACKGROUND: MicroRNAs (miRNAs) play a key role in mediating the action of insulin on cell growth and the development of diabetes. However, few studies have been conducted to provide a comprehensive overview of the miRNA-mediated signaling network in response to glucose in pancreatic beta cells. In our study, we established a computational framework integrating multi-omics profiles analyses, including RNA sequencing (RNA-seq) and small RNA sequencing (sRNA-seq) data analysis, inverse expression pattern analysis, public data integration, and miRNA targets prediction to illustrate the miRNA-mediated regulatory network at different glucose concentrations in INS-1 pancreatic beta cells (INS-1), which display important characteristics of the pancreatic beta cells. RESULTS: We applied our computational framework to the expression profiles of miRNA/mRNA of INS-1, at different glucose concentrations. A total of 1437 differentially expressed genes (DEGs) and 153 differentially expressed miRNAs (DEmiRs) were identified from multi-omics profiles. In particular, 121 DEmiRs putatively regulated a total of 237 DEGs involved in glucose metabolism, fatty acid oxidation, ion channels, exocytosis, homeostasis, and insulin gene regulation. Moreover, Argonaute 2 immunoprecipitation sequencing, qRT-PCR, and luciferase assay identified Crem, Fn1, and Stc1 are direct targets of miR-146b and elucidated that miR-146b acted as a potential regulator and promising target to understand the insulin signaling network. CONCLUSIONS: In this study, the integration of experimentally verified data with system biology framework extracts the miRNA network for exploring potential insulin-associated miRNA and their target genes. The findings offer a potentially significant effect on the understanding of miRNA-mediated insulin signaling network in the development and progression of pancreatic diabetes.

miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions.

MicroRNAs (miRNAs) are small non-coding RNAs of ∼ 22 nucleotides that are involved in negative regulation of mRNA at the post-transcriptional level. Previously, we developed miRTarBase which provides information about experimentally validated miRNA-target interactions (MTIs). Here, we describe an updated database containing 422 517 curated MTIs from 4076 miRNAs and 23 054 target genes collected from over 8500 articles. The number of MTIs curated by strong evidence has increased ∼1.4-fold since the last update in 2016. In this updated version, target sites validated by reporter assay that are available in the literature can be downloaded. The target site sequence can extract new features for analysis via a machine learning approach which can help to evaluate the performance of miRNA-target prediction tools. Furthermore, different ways of browsing enhance user browsing specific MTIs. With these improvements, miRTarBase serves as more comprehensively annotated, experimentally validated miRNA-target interactions databases in the field of miRNA related research. miRTarBase is available at http://miRTarBase.mbc.nctu.edu.tw/.

The aquatic animals' transcriptome resource for comparative functional analysis.

BACKGROUND: Aquatic animals have great economic and ecological importance. Among them, non-model organisms have been studied regarding eco-toxicity, stress biology, and environmental adaptation. Due to recent advances in next-generation sequencing techniques, large amounts of RNA-seq data for aquatic animals are publicly available. However, currently there is no comprehensive resource exist for the analysis, unification, and integration of these datasets. This study utilizes computational approaches to build a new resource of transcriptomic maps for aquatic animals. This aquatic animal transcriptome map database dbATM provides de novo assembly of transcriptome, gene annotation and comparative analysis of more than twenty aquatic organisms without draft genome. RESULTS: To improve the assembly quality, three computational tools (Trinity, Oases and SOAPdenovo-Trans) were employed to enhance individual transcriptome assembly, and CAP3 and CD-HIT-EST software were then used to merge these three assembled transcriptomes. In addition, functional annotation analysis provides valuable clues to gene characteristics, including full-length transcript coding regions, conserved domains, gene ontology and KEGG pathways. Furthermore, all aquatic animal genes are essential for comparative genomics tasks such as constructing homologous gene groups and blast databases and phylogenetic analysis. CONCLUSION: In conclusion, we establish a resource for non model organism aquatic animals, which is great economic and ecological importance and provide transcriptomic information including functional annotation and comparative transcriptome analysis. The database is now publically accessible through the URL http://dbATM.mbc.nctu.edu.tw/ .

Extension of C. elegans lifespan using the ·NO-delivery dinitrosyl iron complexes.

The ubiquitous and emerging physiology function of endogenous nitric oxide in vascular, myocardial, immune, and neuronal systems prompts chemists to develop a prodrug for the controlled delivery of ·NO in vivo and for the translational biomedical application. Inspired by the discovery of natural [Fe(NO)2] motif, herein, we develop the synthetic dinitrosyl iron complexes (DNICs) [Fe2(µ-SR)2(NO)4] (1) as a universal platform for the O2-triggered release of ·NO, for the regulation of ·NO-release kinetics (half-life = 0.6-27.4 h), and for the activation of physiological function of ·NO. Using C. elegans as a model organism, the ·NO-delivery DNIC 1 regulates IIS signaling pathway, AMPK signaling pathway, and mitochondrial function pathway to extend the lifespan and to delay the aging process based on the lifespan analysis, SA-ßgal activity assay, and next-generation RNA sequencing analysis. This study unveils the anti-aging effect of ·NO and develops DNICs as a chemical biology probe for the continued discovery of unprecedented NO physiology.

Synthesis and Photophysical Characterization of 2,3-Dihydroquinolin-4-imines: New Fluorophores with Color-Tailored Emission.

In this study, a series of variously substituted 2,3-dihydroquinolin-4-imines (DQIs) were synthesized from N-substituted propargylanilines by copper(I)-catalyzed annulation. The approach adopted in this study under mild, effective conditions exhibited broad substrate tolerance, particularly for functional groups substituted on anilines. Most of the DQI derivatives synthesized under optimal conditions were obtained in good isolated yields of 63-88 %. 2,3-Dihydroquinolinimine thus obtained was easily converted to important structures like 2,3-dihydroquinolone and tetrahydrobenzodiazepin-5-one, confirming the importance of this strategy in constructing various heterocycles. Surprisingly, 2,3-dihydroquinolinimines thus obtained exhibited bright fluorescence with quantum yields up to 66 %. The density functional theory (DFT) and time-dependent DFT (TD-DFT) calculations were performed for understanding the excited-state nature of DQI system. Accordingly, a tailored DQI derivative bearing methoxy group at C-6 position and acetoxy group at C-7 position was designed and synthesized to give emission at 559 nm with redshift compared to the 7-methoxy substituted DQI. A detailed study of DQI structures with their photophysical properties was performed with five control molecules and consequently demonstrated the uniqueness of the chemical structures of DQIs.

Glycan-Directed Grafting-from Polymerization of Immunoglobulin G: Site-Selectively Modified IgG-Polymer Conjugates with Preserved Biological Activity.

Antibody and related antibody drugs for the treatment of malignancies have led to progress in targeted cancer therapy. Preparation of diverse antibody conjugates is critical for preclinical and clinical applications. However, precise control in tagging molecules at specific locations on antibodies is essential to preserve their native function. In this study, a synthetic boronic acid (BA)-tosyl initiator was used to trigger a glycan-directed modification of IgGs, and the obtained IgG macroinitiators allowed a growth of the poly N-isopropylacrylamide (PNIPAAm) chains specifically at Fc-domains. Therefore, the PNIPAAm chains are located away from the critical antigen-binding domains (Fab), which could reasonably prevent the loss of biological activity after the attachment of polymer chains. According to the proposed strategy, a site-selectively modified anticoncanavalin A (Con A) antibody-PNIPAAm conjugate showed 6-times higher efficiency in the binding of targeted Con A antigen to a randomly conjugated anti-Con A antibody-PNIPAAm conjugate. In this study, we developed the first chemical strategy for the site-specific preparation of IgG-polymer conjugates with conserved biological activity as well as intact glycan structures.

miRTarBase 2016: updates to the experimentally validated miRNA-target interactions database.

MicroRNAs (miRNAs) are small non-coding RNAs of approximately 22 nucleotides, which negatively regulate the gene expression at the post-transcriptional level. This study describes an update of the miRTarBase (http://miRTarBase.mbc.nctu.edu.tw/) that provides information about experimentally validated miRNA-target interactions (MTIs). The latest update of the miRTarBase expanded it to identify systematically Argonaute-miRNA-RNA interactions from 138 crosslinking and immunoprecipitation sequencing (CLIP-seq) data sets that were generated by 21 independent studies. The database contains 4966 articles, 7439 strongly validated MTIs (using reporter assays or western blots) and 348 007 MTIs from CLIP-seq. The number of MTIs in the miRTarBase has increased around 7-fold since the 2014 miRTarBase update. The miRNA and gene expression profiles from The Cancer Genome Atlas (TCGA) are integrated to provide an effective overview of this exponential growth in the miRNA experimental data. These improvements make the miRTarBase one of the more comprehensively annotated, experimentally validated miRNA-target interactions databases and motivate additional miRNA research efforts.

A Three-MicroRNA Signature as a Potential Biomarker for the Early Detection of Oral Cancer.

Oral squamous cell carcinoma (OSCC) is often diagnosed at a late stage and may be malignantly transformed from oral leukoplakia (OL). This study aimed to identify potential plasma microRNAs (miRNAs) for the early detection of oral cancer. Plasma from normal, OL, and OSCC patients were evaluated. Small RNA sequencing was used to screen the differently expressed miRNAs among the groups. Next, these miRNAs were validated with individual samples by quantitative real-time polymerase chain reaction (qRT-PCR) assays in the training phase (n = 72) and validation phase (n = 178). The possible physiological roles of the identified miRNAs were further investigated using bioinformatics analysis. Three miRNAs (miR-222-3p, miR-150-5p, and miR-423-5p) were identified as differentially expressed among groups; miR-222-3p and miR-423-5p negatively correlated with T stage, lymph node metastasis status, and clinical stage. A high diagnostic accuracy (Area under curve = 0.88) was demonstrated for discriminating OL from OSCC. Bioinformatics analysis reveals that miR-423-5p and miR-222-3p are significantly over-expressed in oral cancer tissues and involved in various cancer pathways. The three-plasma miRNA panel may be useful to monitor malignant progression from OL to OSCC and as potential biomarkers for early detection of oral cancer.

MethHC: a database of DNA methylation and gene expression in human cancer.

We present MethHC (http://MethHC.mbc.nctu.edu.tw), a database comprising a systematic integration of a large collection of DNA methylation data and mRNA/microRNA expression profiles in human cancer. DNA methylation is an important epigenetic regulator of gene transcription, and genes with high levels of DNA methylation in their promoter regions are transcriptionally silent. Increasing numbers of DNA methylation and mRNA/microRNA expression profiles are being published in different public repositories. These data can help researchers to identify epigenetic patterns that are important for carcinogenesis. MethHC integrates data such as DNA methylation, mRNA expression, DNA methylation of microRNA gene and microRNA expression to identify correlations between DNA methylation and mRNA/microRNA expression from TCGA (The Cancer Genome Atlas), which includes 18 human cancers in more than 6000 samples, 6548 microarrays and 12 567 RNA sequencing data.

Integrated MicroRNA-mRNA Analysis Reveals miR-204 Inhibits Cell Proliferation in Gastric Cancer by Targeting CKS1B, CXCL1 and GPRC5A.

Gastric cancer (GC) is the second most frequent cause of cancer-related deaths worldwide. MicroRNAs are single-stranded RNA molecules of 21-23 nucleotides that regulate target gene expression through specific base-pairing interactions between miRNA and untranslated regions of targeted mRNAs. In this study, we generated a multistep approach for the integrated analysis of miRNA and mRNA expression. First, both miRNA and mRNA expression profiling datasets in gastric cancer from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) identified 79 and 1042 differentially expressed miRNAs and mRNAs, respectively, in gastric cancer. Second, inverse correlations between miRNA and mRNA expression levels identified 3206 miRNA-mRNA pairs combined with 79 dysregulated miRNAs and their 774 target mRNAs predicted by three prediction tools, miRanda, PITA, and RNAhybrid. Additionally, miR-204, which was found to be down-regulated in gastric cancer, was ectopically over-expressed in the AGS gastric cancer cell line and all down-regulated targets were identified by RNA sequencing (RNA-seq) analysis. Over-expression of miR-204 reduced the gastric cancer cell proliferation and suppressed the expression of three targets which were validated by qRT-PCR and luciferase assays. For the first time, we identified that CKS1B, CXCL1, and GPRC5A are putative targets of miR-204 and elucidated that miR-204 acted as potential tumor suppressor and, therefore, are useful as a promising therapeutic target for gastric cancer.

Transcriptome sequencing based annotation and homologous evidence based scaffolding of Anguilla japonica draft genome.

BACKGROUND: Anguilla japonica (Japanese eel) is currently one of the most important research subjects in eastern Asia aquaculture. Enigmatic life cycle of the organism makes study of artificial reproduction extremely limited. Henceforth genomic and transcriptomic resources of eels are urgently needed to help solving the problems surrounding this organism across multiple fields. We hereby provide a reconstructed transcriptome from deep sequencing of juvenile (glass eels) whole body samples. The provided expressed sequence tags were used to annotate the currently available draft genome sequence. Homologous information derived from the annotation result was applied to improve the group of scaffolds into available linkage groups. RESULTS: With the transcriptome sequence data combined with publicly available expressed sequence tags evidences, 18,121 genes were structurally and functionally annotated on the draft genome. Among them, 3,921 genes were located in the 19 linkage groups. 137 scaffolds covering 13 million bases were grouped into the linkage groups in additional to the original partial linkage groups, increasing the linkage group coverage from 13 to 14%. CONCLUSIONS: This annotation provide information of the coding regions of the genes supported by transcriptome based evidence. The derived homologous evidences pave the way for phylogenetic analysis of important genetic traits and the improvement of the genome assembly.

miRTarBase update 2014: an information resource for experimentally validated miRNA-target interactions.

MicroRNAs (miRNAs) are small non-coding RNA molecules capable of negatively regulating gene expression to control many cellular mechanisms. The miRTarBase database (http://mirtarbase.mbc.nctu.edu.tw/) provides the most current and comprehensive information of experimentally validated miRNA-target interactions. The database was launched in 2010 with data sources for >100 published studies in the identification of miRNA targets, molecular networks of miRNA targets and systems biology, and the current release (2013, version 4) includes significant expansions and enhancements over the initial release (2010, version 1). This article reports the current status of and recent improvements to the database, including (i) a 14-fold increase to miRNA-target interaction entries, (ii) a miRNA-target network, (iii) expression profile of miRNA and its target gene, (iv) miRNA target-associated diseases and (v) additional utilities including an upgrade reminder and an error reporting/user feedback system.

The evolutionary dynamics of tRNA-gene copy number and codon-use in E. coli.

BACKGROUND: The introduction of foreign DNA by Lateral Gene Transfer (LGT) can quickly and drastically alter genome composition. Problems can arise if the genes introduced by LGT use codons that are not suited to the host's translational machinery. Here we investigate compensatory adaptation of E. coli in response to the introduction of large volumes of codons that are rarely used by the host genome. RESULTS: We analyze genome sequences from the E. coli/Shigella complex, and find that certain tRNA genes are present in multiple copies in two pathogenic Shigella and O157:H7 subgroups of E. coli. Furthermore, we show that the codons that correspond to these multi-copy number tRNA genes are enriched in the high copy number Selfish Genetic Elements (SGE's) in Shigella and laterally introduced genes in O157:H7. We analyze the duplicate copies and find evidence for the selective retention of tRNA genes introduced by LGT in response to the changed codon content of the genome. CONCLUSION: These data support a model where the relatively rapid influx of LGT genes and SGE's introduces a large number of genes maladapted to the host's translational machinery. Under these conditions, it becomes advantageous for the host to retain tRNA genes that are required for the incorporation of amino acids at these codons. Subsequently, the increased number of copies of these specific tRNA genes adjusts the cellular tRNA pool to the demands set by global shifts in codon usage.

Integrated analyses to reconstruct microRNA-mediated regulatory networks in mouse liver using high-throughput profiling.

BACKGROUND: MicroRNAs (miRNAs) simultaneously target many transcripts through partial complementarity binding, and have emerged as a key type of post-transcriptional regulator for gene expression. How miRNA accomplishes its pleiotropic effects largely depends on its expression and its target repertoire. Previous studies discovered thousands of miRNAs and numerous miRNA target genes mainly through computation and prediction methods which produced high rates of false positive prediction. The development of Argonaute cross-linked immunoprecipitation coupled with high-throughput sequencing (CLIP-Seq) provides a system to effectively determine miRNA target genes. Likewise, the accuracy of dissecting the transcriptional regulation of miRNA genes has been greatly improved by chromatin immunoprecipitation of the transcription factors coupled with sequencing (ChIP-Seq). Elucidation of the miRNA target repertoire will provide an in-depth understanding of the functional roles of microRNA pathways. To reliably reconstruct a miRNA-mediated regulatory network, we established a computational framework using publicly available, sequence-based transcription factor-miRNA databases, including ChIPBase and TransmiR for the TF-miRNA interactions, along with miRNA-target databases, including miRTarBase, TarBase and starBase, for the miRNA-target interactions. We applied the computational framework to elucidate the miRNA-mediated regulatory network in the Mir122a⁻/⁻ mouse model, which has an altered transcriptome and progressive liver disease. RESULTS: We applied our computational framework to the expression profiles of miRNA/mRNA of Mir122a⁻/⁻ mutant mice and wild-type mice. The miRNA-mediated network involves 40 curated TFs contributing to the aberrant expression of 65 miRNAs and 723 curated miRNA target genes, of which 56% was found in the differentially-expressed genes of Mir122a--mice. Hence, the regulatory network disclosed previously-known and also many previously-unidentified miRNA-mediated regulations in mutant mice. Moreover, we demonstrate that loss of imprinting at the chromosome 12qF1 region is associated with miRNA overexpression in human hepatocellular carcinoma and stem cells, suggesting initiation of precancerous changes in young mice deficient in miR-122. A group of 9 miRNAs was found to share miR-122 target genes, indicating synergy between miRNAs and target genes by way of multiplicity and cooperativity. CONCLUSIONS: The study provides significant insight into miRNA-mediated regulatory networks. Based on experimentally verified data, this network is highly reliable and effective in revealing previously-undetermined disease-associated molecular mechanisms. This computational framework can be applied to explore the significant TF-miRNA-miRNA target interactions in any complex biological systems with high degrees of confidence.

A computational approach for identifying microRNA-target interactions using high-throughput CLIP and PAR-CLIP sequencing.

BACKGROUND: MicroRNAs (miRNAs) play a critical role in down-regulating gene expression. By coupling with Argonaute family proteins, miRNAs bind to target sites on mRNAs and employ translational repression. A large amount of miRNA-target interactions (MTIs) have been identified by the crosslinking and immunoprecipitation (CLIP) and the photoactivatable-ribonucleoside-enhanced CLIP (PAR-CLIP) along with the next-generation sequencing (NGS). PAR-CLIP shows high efficiency of RNA co-immunoprecipitation, but it also lead to T to C conversion in miRNA-RNA-protein crosslinking regions. This artificial error obviously reduces the mappability of reads. However, a specific tool to analyze CLIP and PAR-CLIP data that takes T to C conversion into account is still in need. RESULTS: We herein propose the first CLIP and PAR-CLIP sequencing analysis platform specifically for miRNA target analysis, namely miRTarCLIP. From scratch, it automatically removes adaptor sequences from raw reads, filters low quality reads, reverts C to T, aligns reads to 3'UTRs, scans for read clusters, identifies high confidence miRNA target sites, and provides annotations from external databases. With multi-threading techniques and our novel C to T reversion procedure, miRTarCLIP greatly reduces the running time comparing to conventional approaches. In addition, miRTarCLIP serves with a web-based interface to provide better user experiences in browsing and searching targets of interested miRNAs. To demonstrate the superior functionality of miRTarCLIP, we applied miRTarCLIP to two public available CLIP and PAR-CLIP sequencing datasets. miRTarCLIP not only shows comparable results to that of other existing tools in a much faster speed, but also reveals interesting features among these putative target sites. Specifically, we used miRTarCLIP to disclose that T to C conversion within position 1-7 and that within position 8-14 of miRNA target sites are significantly different (p value = 0.02), and even more significant when focusing on sites targeted by top 102 highly expressed miRNAs only (p value = 0.01). These results comply with previous findings and further suggest that combining miRNA expression and PAR-CLIP data can improve accuracy of the miRNA target prediction. CONCLUSION: To sum up, we devised a systematic approach for mining miRNA-target sites from CLIP-seq and PAR-CLIP sequencing data, and integrated the workflow with a graphical web-based browser, which provides a user friendly interface and detailed annotations of MTIs. We also showed through real-life examples that miRTarCLIP is a powerful tool for understanding miRNAs. Our integrated tool can be accessed online freely at http://miRTarCLIP.mbc.nctu.edu.tw.