Endurance exercise training-responsive miR-19b-3p improves skeletal muscle glucose metabolism.

Skeletal muscle is a highly adaptable tissue and remodels in response to exercise training. Using short RNA sequencing, we determine the miRNA profile of skeletal muscle from healthy male volunteers before and after a 14-day aerobic exercise training regime. Among the exercise training-responsive miRNAs identified, miR-19b-3p was selected for further validation. Overexpression of miR-19b-3p in human skeletal muscle cells increases insulin signaling, glucose uptake, and maximal oxygen consumption, recapitulating the adaptive response to aerobic exercise training. Overexpression of miR-19b-3p in mouse flexor digitorum brevis muscle enhances contraction-induced glucose uptake, indicating that miR-19b-3p exerts control on exercise training-induced adaptations in skeletal muscle. Potential targets of miR-19b-3p that are reduced after aerobic exercise training include KIF13A, MAPK6, RNF11, and VPS37A. Amongst these, RNF11 silencing potentiates glucose uptake in human skeletal muscle cells. Collectively, we identify miR-19b-3p as an aerobic exercise training-induced miRNA that regulates skeletal muscle glucose metabolism.

Circulating Exosomal miR-20b-5p Is Elevated in Type 2 Diabetes and Could Impair Insulin Action in Human Skeletal Muscle.

miRNAs are noncoding RNAs representing an important class of gene expression modulators. Extracellular circulating miRNAs are both candidate biomarkers for disease pathogenesis and mediators of cell-to-cell communication. We examined the miRNA expression profile of total serum and serum-derived exosome-enriched extracellular vesicles in people with normal glucose tolerance or type 2 diabetes. In contrast to total serum miRNA, which did not reveal any differences in miRNA expression, we identified differentially abundant miRNAs in patients with type 2 diabetes using miRNA expression profiles of exosome RNA (exoRNA). To validate the role of these differentially abundant miRNAs on glucose metabolism, we transfected miR-20b-5p, a highly abundant exoRNA in patients with type 2 diabetes, into primary human skeletal muscle cells. miR-20b-5p overexpression increased basal glycogen synthesis in human skeletal muscle cells. We identified AKTIP and STAT3 as miR-20b-5p targets. miR-20b-5p overexpression reduced AKTIP abundance and insulin-stimulated glycogen accumulation. In conclusion, exosome-derived extracellular miR-20b-5p is a circulating biomarker associated with type 2 diabetes that plays an intracellular role in modulating insulin-stimulated glucose metabolism via AKT signaling.

FANTOM5 CAGE profiles of human and mouse samples.

In the FANTOM5 project, transcription initiation events across the human and mouse genomes were mapped at a single base-pair resolution and their frequencies were monitored by CAGE (Cap Analysis of Gene Expression) coupled with single-molecule sequencing. Approximately three thousands of samples, consisting of a variety of primary cells, tissues, cell lines, and time series samples during cell activation and development, were subjected to a uniform pipeline of CAGE data production. The analysis pipeline started by measuring RNA extracts to assess their quality, and continued to CAGE library production by using a robotic or a manual workflow, single molecule sequencing, and computational processing to generate frequencies of transcription initiation. Resulting data represents the consequence of transcriptional regulation in each analyzed state of mammalian cells. Non-overlapping peaks over the CAGE profiles, approximately 200,000 and 150,000 peaks for the human and mouse genomes, were identified and annotated to provide precise location of known promoters as well as novel ones, and to quantify their activities.

microManaging glucose and lipid metabolism in skeletal muscle: Role of microRNAs.

MicroRNAs have been described as important regulators of skeletal muscle differentiation and development, but the role of microRNAs in glucose and lipid metabolism is less well established. Here we review the microRNAs involved in insulin resistance and glucose metabolism, as well as microRNAs regulating lipid metabolism and mitochondrial functions in skeletal muscle, with an emphasis on metabolic disorders such as type 2 diabetes and the adaptive response to exercise training. Finally, we raise some methodological considerations for studying microRNAs, as well as challenges investigators may face when elucidating the direct role of microRNAs in the regulation of glucose and lipid metabolism in skeletal muscle. This article is part of a Special Issue entitled: MicroRNAs and lipid/energy metabolism and related diseases edited by Carlos Fernández-Hernando and Yajaira Suárez.

miRNA let-7 expression is regulated by glucose and TNF-α by a remote upstream promoter.

miRNAs regulate protein abundance and control diverse aspects of cellular processes and biological functions in metabolic diseases, such as obesity and type 2 diabetes (T2D). Let (lethal)-7 miRNAs specifically targets genes associated with T2D and have been implicated in the regulation of peripheral glucose metabolism, yet the direct regulators of let-7 miRNA expression are unknown. In the present study, we report on a putative promoter region for the let-7a-1, let-7f-1 and let-7d gene cluster on chromosome 9 and characterize the promoter activity of this novel area. We show that promoter activity and let-7 miRNA expression is dynamically regulated in response to different factors including serum, glucose, tumour necrosis factor (TNF)-α and caffeine. These findings will contribute to understanding the interaction between precise promoter elements to control the transcription and translation of let-7 miRNA genes.

Temporal analysis of reciprocal miRNA-mRNA expression patterns predicts regulatory networks during differentiation in human skeletal muscle cells.

microRNAs (miRNAs) are short noncoding RNAs that regulate gene expression through posttranscriptional repression of target genes. miRNAs exert a fundamental level of control over many developmental processes, but their role in the differentiation and development of skeletal muscle from myogenic progenitor cells in humans remains incompletely understood. Using primary cultures established from human skeletal muscle satellite cells, we performed microarray profiling of miRNA expression during differentiation of myoblasts (day 0) into myotubes at 48 h intervals (day 2, 4, 6, 8, and 10). Based on a time-course analysis, we identified 44 miRNAs with altered expression [false discovery rate (FDR) < 5%, fold change > ±1.2] during differentiation, including the marked upregulation of the canonical myogenic miRNAs miR-1, miR-133a, miR-133b, and miR-206. Microarray profiling of mRNA expression at day 0, 4, and 10 identified 842 and 949 genes differentially expressed (FDR < 10%) at day 4 and 10, respectively. At day 10, 42% of altered transcripts demonstrated reciprocal expression patterns in relation to the directional change of their in silico predicted regulatory miRNAs based on analysis using Ingenuity Pathway Analysis microRNA Target Filter. Bioinformatic analysis predicted networks of regulation during differentiation including myomiRs miR-1/206 and miR-133a/b, miRNAs previously established in differentiation including miR-26 and miR-30, and novel miRNAs regulated during differentiation of human skeletal muscle cells such as miR-138-5p and miR-20a. These reciprocal expression patterns may represent new regulatory nodes in human skeletal muscle cell differentiation. This analysis serves as a reference point for future studies of human skeletal muscle differentiation and development in healthy and disease states.

Proteasome inhibition in skeletal muscle cells unmasks metabolic derangements in type 2 diabetes.

Two-dimensional difference gel electrophoresis (2-D DIGE)-based proteome analysis has revealed intrinsic insulin resistance in myotubes derived from type 2 diabetic patients. Using 2-D DIGE-based proteome analysis, we identified a subset of insulin-resistant proteins involved in protein turnover in skeletal muscle of type 2 diabetic patients, suggesting aberrant regulation of the protein homeostasis maintenance system underlying metabolic disease. We then validated the role of the ubiquitin-proteasome system (UPS) in myotubes to investigate whether impaired proteasome function may lead to metabolic arrest or insulin resistance. Myotubes derived from muscle biopsies obtained from people with normal glucose tolerance (NGT) or type 2 diabetes were exposed to the proteasome inhibitor bortezomib (BZ; Velcade) without or with insulin. BZ exposure increased protein carbonylation and lactate production yet impaired protein synthesis and UPS function in myotubes from type 2 diabetic patients, marking the existence of an insulin-resistant signature that was retained in cultured myotubes. In conclusion, BZ treatment further exacerbates insulin resistance and unmasks intrinsic features of metabolic disease in myotubes derived from type 2 diabetic patients. Our results highlight the existence of a confounding inherent abnormality in cellular protein dynamics in metabolic disease, which is uncovered through concurrent inhibition of the proteasome system.

Comparison of CAGE and RNA-seq transcriptome profiling using clonally amplified and single-molecule next-generation sequencing.

CAGE (cap analysis gene expression) and RNA-seq are two major technologies used to identify transcript abundances as well as structures. They measure expression by sequencing from either the 5' end of capped molecules (CAGE) or tags randomly distributed along the length of a transcript (RNA-seq). Library protocols for clonally amplified (Illumina, SOLiD, 454 Life Sciences [Roche], Ion Torrent), second-generation sequencing platforms typically employ PCR preamplification prior to clonal amplification, while third-generation, single-molecule sequencers can sequence unamplified libraries. Although these transcriptome profiling platforms have been demonstrated to be individually reproducible, no systematic comparison has been carried out between them. Here we compare CAGE, using both second- and third-generation sequencers, and RNA-seq, using a second-generation sequencer based on a panel of RNA mixtures from two human cell lines to examine power in the discrimination of biological states, detection of differentially expressed genes, linearity of measurements, and quantification reproducibility. We found that the quantified levels of gene expression are largely comparable across platforms and conclude that CAGE and RNA-seq are complementary technologies that can be used to improve incomplete gene models. We also found systematic bias in the second- and third-generation platforms, which is likely due to steps such as linker ligation, cleavage by restriction enzymes, and PCR amplification. This study provides a perspective on the performance of these platforms, which will be a baseline in the design of further experiments to tackle complex transcriptomes uncovered in a wide range of cell types.

Autocrine role of interleukin-13 on skeletal muscle glucose metabolism in type 2 diabetic patients involves microRNA let-7.

Low-grade inflammation associated with type 2 diabetes (T2DM) is postulated to exacerbate insulin resistance. We report that serum levels, as well as IL-13 secreted from cultured skeletal muscle, are reduced in T2DM vs. normal glucose-tolerant (NGT) subjects. IL-13 exposure increases skeletal muscle glucose uptake, oxidation, and glycogen synthesis via an Akt-dependent mechanism. Expression of microRNA let-7a and let-7d, which are direct translational repressors of the IL-13 gene, was increased in skeletal muscle from T2DM patients. Overexpression of let-7a and let-7d in cultured myotubes reduced IL-13 secretion. Furthermore, basal glycogen synthesis was reduced in cultured myotubes exposed to an IL-13-neutralizing antibody. Thus, IL-13 is synthesized and released by skeletal muscle through a mechanism involving let-7, and this effect is attenuated in skeletal muscle from insulin-resistant T2DM patients. In conclusion, IL-13 plays an autocrine role in skeletal muscle to increase glucose uptake and metabolism, suggesting a role in glucose homeostasis in metabolic disease.

Automated workflow for preparation of cDNA for cap analysis of gene expression on a single molecule sequencer.

BACKGROUND: Cap analysis of gene expression (CAGE) is a 5' sequence tag technology to globally determine transcriptional starting sites in the genome and their expression levels and has most recently been adapted to the HeliScope single molecule sequencer. Despite significant simplifications in the CAGE protocol, it has until now been a labour intensive protocol. METHODOLOGY: In this study we set out to adapt the protocol to a robotic workflow, which would increase throughput and reduce handling. The automated CAGE cDNA preparation system we present here can prepare 96 'HeliScope ready' CAGE cDNA libraries in 8 days, as opposed to 6 weeks by a manual operator.We compare the results obtained using the same RNA in manual libraries and across multiple automation batches to assess reproducibility. CONCLUSIONS: We show that the sequencing was highly reproducible and comparable to manual libraries with an 8 fold increase in productivity. The automated CAGE cDNA preparation system can prepare 96 CAGE sequencing samples simultaneously. Finally we discuss how the system could be used for CAGE on Illumina/SOLiD platforms, RNA-seq and full-length cDNA generation.

LRRN4 and UPK3B are markers of primary mesothelial cells.

BACKGROUND: Mesothelioma is a highly malignant tumor that is primarily caused by occupational or environmental exposure to asbestos fibers. Despite worldwide restrictions on asbestos usage, further cases are expected as diagnosis is typically 20-40 years after exposure. Once diagnosed there is a very poor prognosis with a median survival rate of 9 months. Considering this the development of early pre clinical diagnostic markers may help improve clinical outcomes. METHODOLOGY: Microarray expression arrays on mesothelium and other tissues dissected from mice were used to identify candidate mesothelial lineage markers. Candidates were further tested by qRTPCR and in-situ hybridization across a mouse tissue panel. Two candidate biomarkers with the potential for secretion, uroplakin 3B (UPK3B), and leucine rich repeat neuronal 4 (LRRN4) and one commercialized mesothelioma marker, mesothelin (MSLN) were then chosen for validation across a panel of normal human primary cells, 16 established mesothelioma cell lines, 10 lung cancer lines, and a further set of 8 unrelated cancer cell lines. CONCLUSIONS: Within the primary cell panel, LRRN4 was only detected in primary mesothelial cells, but MSLN and UPK3B were also detected in other cell types. MSLN was detected in bronchial epithelial cells and alveolar epithelial cells and UPK3B was detected in retinal pigment epithelial cells and urothelial cells. Testing the cell line panel, MSLN was detected in 15 of the 16 mesothelioma cells lines, whereas LRRN4 was only detected in 8 and UPK3B in 6. Interestingly MSLN levels appear to be upregulated in the mesothelioma lines compared to the primary mesothelial cells, while LRRN4 and UPK3B, are either lost or down-regulated. Despite the higher fraction of mesothelioma lines positive for MSLN, it was also detected at high levels in 2 lung cancer lines and 3 other unrelated cancer lines derived from papillotubular adenocarcinoma, signet ring carcinoma and transitional cell carcinoma.

Unamplified cap analysis of gene expression on a single-molecule sequencer.

We report the development of a simplified cap analysis of gene expression (CAGE) protocol adapted for single-molecule sequencers that avoids second strand synthesis, ligation, digestion, and PCR. HeliScopeCAGE directly sequences the 3' end of cap trapped first-strand cDNAs. As with previous versions of CAGE, we better define transcription start sites (TSS) than known models, identify novel regions of transcription and alternative promoters, and find two major classes of TSS signal, sharp peaks and broad regions. However, using this protocol, we observe reproducible evidence of regulation at the much finer level of individual TSS positions. The libraries are quantitative over 5 orders of magnitude and highly reproducible (Pearson's correlation coefficient of 0.987). We have also scaled down the sample requirement to 5 µg of total RNA for a standard HeliScopeCAGE library and 100 ng for a low-quantity version. When the same RNA was run as 5-µg and 100-ng versions, the 100 ng was still able to detect expression for ∼60% of the 13,468 loci detected by a 5-µg library using the same threshold, allowing comparative analysis of even rare cell populations. Testing the protocol for differential gene expression measurements on triplicate HeLa and THP-1 samples, we find that the log fold change compared to Illumina microarray measurements is highly correlated (0.871). In addition, HeliScopeCAGE finds differential expression for thousands more loci including those with probes on the array. Finally, although the majority of tags are 5' associated, we also observe a low level of signal on exons that is useful for defining gene structures.

The combination of gene perturbation assay and ChIP-chip reveals functional direct target genes for IRF8 in THP-1 cells.

Gene regulatory networks in living cells are controlled by the interaction of multiple cell type-specific transcription regulators with DNA binding sites in target genes. Interferon regulatory factor 8 (IRF8), also known as interferon consensus sequence binding protein (ICSBP), is a transcription factor expressed predominantly in myeloid and lymphoid cell lineages. To find the functional direct target genes of IRF8, the gene expression profiles of siRNA knockdown samples and genome-wide binding locations by ChIP-chip were analyzed in THP-1 myelomonocytic leukemia cells. Consequently, 84 genes were identified as functional direct targets. The ETS family transcription factor PU.1, also known as SPI1, binds to IRF8 and regulates basal transcription in macrophages. Using the same approach, we identified 53 direct target genes of PU.1; these overlapped with 19 IRF8 targets. These 19 genes included key molecules of IFN signaling such as OAS1 and IRF9, but excluded other IFN-related genes amongst the IRF8 functional direct target genes. We suggest that IRF8 and PU.1 can have both combined, and independent actions on different promoters in myeloid cells.

An atlas of combinatorial transcriptional regulation in mouse and man.

Combinatorial interactions among transcription factors are critical to directing tissue-specific gene expression. To build a global atlas of these combinations, we have screened for physical interactions among the majority of human and mouse DNA-binding transcription factors (TFs). The complete networks contain 762 human and 877 mouse interactions. Analysis of the networks reveals that highly connected TFs are broadly expressed across tissues, and that roughly half of the measured interactions are conserved between mouse and human. The data highlight the importance of TF combinations for determining cell fate, and they lead to the identification of a SMAD3/FLI1 complex expressed during development of immunity. The availability of large TF combinatorial networks in both human and mouse will provide many opportunities to study gene regulation, tissue differentiation, and mammalian evolution.

Cross-mapping and the identification of editing sites in mature microRNAs in high-throughput sequencing libraries.

MicroRNAs (miRNAs) are short (20-23 nt) RNAs that are sequence-specific mediators of transcriptional and post-transcriptional regulation of gene expression. Modern high-throughput technologies enable deep sequencing of such RNA species on an unprecedented scale. We find that the analysis of small RNA deep-sequencing libraries can be affected by cross-mapping, in which RNA sequences originating from one locus are inadvertently mapped to another. Similar to cross-hybridization on microarrays, cross-mapping is prevalent among miRNAs, as they tend to occur in families, are similar or derived from repeat or structural RNAs, or are post-transcriptionally modified. Here, we develop a strategy to correct for cross-mapping, and apply it to the analysis of RNA editing in mature miRNAs. In contrast to previous reports, our analysis suggests that RNA editing in mature miRNAs is rare in animals.

The transcriptional network that controls growth arrest and differentiation in a human myeloid leukemia cell line.

Using deep sequencing (deepCAGE), the FANTOM4 study measured the genome-wide dynamics of transcription-start-site usage in the human monocytic cell line THP-1 throughout a time course of growth arrest and differentiation. Modeling the expression dynamics in terms of predicted cis-regulatory sites, we identified the key transcription regulators, their time-dependent activities and target genes. Systematic siRNA knockdown of 52 transcription factors confirmed the roles of individual factors in the regulatory network. Our results indicate that cellular states are constrained by complex networks involving both positive and negative regulatory interactions among substantial numbers of transcription factors and that no single transcription factor is both necessary and sufficient to drive the differentiation process.

Gemin2 plays an important role in stabilizing the survival of motor neuron complex.

The survival of motor neuron (SMN) protein, responsible for the neurodegenerative disease spinal muscular atrophy (SMA), oligomerizes and forms a stable complex with seven other major components, the Gemin proteins. Besides the SMN protein, Gemin2 is a core protein that is essential for the formation of the SMN complex, although the mechanism by which it drives formation is unclear. We have found a novel interaction, a Gemin2 self-association, using the mammalian two-hybrid system and the in vitro pull-down assays. Using in vitro dissociation assays, we also found that the self-interaction of the amino-terminal SMN protein, which was confirmed in this study, became stable in the presence of Gemin2. In addition, Gemin2 knockdown using small interference RNA treatment revealed a drastic decrease in SMN oligomer formation and in the assembly activity of spliceosomal small nuclear ribonucleoprotein (snRNP). Taken together, these results indicate that Gemin2 plays an important role in snRNP assembly through the stabilization of the SMN oligomer/complex via novel self-interaction. Applying the results/techniques to amino-terminal SMN missense mutants that were recently identified from SMA patients, we successfully showed that amino-terminal self-association, Gemin2 binding, the stabilization effect of Gemin2, and snRNP assembly activity were all lowered in the mutant SMN(D44V), suggesting that instability of the amino-terminal SMN self-association may cause SMA in patients carrying this allele.

Transcript annotation in FANTOM3: mouse gene catalog based on physical cDNAs.

The international FANTOM consortium aims to produce a comprehensive picture of the mammalian transcriptome, based upon an extensive cDNA collection and functional annotation of full-length enriched cDNAs. The previous dataset, FANTOM2, comprised 60,770 full-length enriched cDNAs. Functional annotation revealed that this cDNA dataset contained only about half of the estimated number of mouse protein-coding genes, indicating that a number of cDNAs still remained to be collected and identified. To pursue the complete gene catalog that covers all predicted mouse genes, cloning and sequencing of full-length enriched cDNAs has been continued since FANTOM2. In FANTOM3, 42,031 newly isolated cDNAs were subjected to functional annotation, and the annotation of 4,347 FANTOM2 cDNAs was updated. To accomplish accurate functional annotation, we improved our automated annotation pipeline by introducing new coding sequence prediction programs and developed a Web-based annotation interface for simplifying the annotation procedures to reduce manual annotation errors. Automated coding sequence and function prediction was followed with manual curation and review by expert curators. A total of 102,801 full-length enriched mouse cDNAs were annotated. Out of 102,801 transcripts, 56,722 were functionally annotated as protein coding (including partial or truncated transcripts), providing to our knowledge the greatest current coverage of the mouse proteome by full-length cDNAs. The total number of distinct non-protein-coding transcripts increased to 34,030. The FANTOM3 annotation system, consisting of automated computational prediction, manual curation, and final expert curation, facilitated the comprehensive characterization of the mouse transcriptome, and could be applied to the transcriptomes of other species.

Genome-wide analysis of mammalian promoter architecture and evolution.

Mammalian promoters can be separated into two classes, conserved TATA box-enriched promoters, which initiate at a well-defined site, and more plastic, broad and evolvable CpG-rich promoters. We have sequenced tags corresponding to several hundred thousand transcription start sites (TSSs) in the mouse and human genomes, allowing precise analysis of the sequence architecture and evolution of distinct promoter classes. Different tissues and families of genes differentially use distinct types of promoters. Our tagging methods allow quantitative analysis of promoter usage in different tissues and show that differentially regulated alternative TSSs are a common feature in protein-coding genes and commonly generate alternative N termini. Among the TSSs, we identified new start sites associated with the majority of exons and with 3' UTRs. These data permit genome-scale identification of tissue-specific promoters and analysis of the cis-acting elements associated with them.