Comparing transcription rate and mRNA abundance as parameters for biochemical pathway and network analysis

Next Generation Sequencing (NGS) Application in Multiparameter Gene Expression Analysis

Next generation sequencing (NGS) is routinely used in gene expression analyses. In particular, RNA-seq has been the method of choice for highly sensitive genome-wide quantification of RNA expression. The method can be used in a wide variety of model systems, including studies to gain insight into underlying mechanisms of toxicologic processes and disease development induced by environmental toxicants. RNA-seq has also been coupled to many other molecular biology protocols to monitor specific aspects of the gene expression process. Here, we describe two such coupling-(a) global run-on sequencing (GRO-seq) that coupled it to the nuclear run-on (NRO), and (b) polysome profiling that coupled it to sucrose-gradient-based polysome isolation. Simultaneous RNA-seq, GRO-seq, and polysome profiling analyses enabled genome-wide analysis of the mode of stability control of individual RNA molecules.

A Multi-Parameter Analysis of Cellular Coordination of Major Transcriptome Regulation Mechanisms

To understand cellular coordination of multiple transcriptome regulation mechanisms, we simultaneously measured transcription rate (TR), mRNA abundance (RA) and translation activity (TA). This revealed multiple insights. First, the three parameters displayed systematic statistical differences. Sequentially more genes exhibited extreme (low or high) expression values from TR to RA, and then to TA; that is, cellular coordination of multiple transcriptome regulatory mechanisms leads to sequentially enhanced gene expression selectivity as the genetic information flow from the genome to the proteome. Second, contribution of the stabilization-by-translation regulatory mechanism to the cellular coordination process was assessed. The data enabled an estimation of mRNA stability, revealing a moderate but significant positive correlation between mRNA stability and translation activity. Third, the proportion of mRNA occupied by un-translated regions (UTR) exhibited a negative relationship with the level of this correlation, and was thus a major determinant of the mode of regulation of the mRNA. High-UTR-proportion mRNAs tend to defy the stabilization-by-translation regulatory mechanism, staying out of the polysome but remaining stable; mRNAs with little UTRs largely followed this regulation. In summary, we quantitatively delineated the relationship among multiple transcriptome regulation parameters, i.e., cellular coordination of corresponding regulatory mechanisms.

Transcriptomic Analysis of Young and Old Erythrocytes of Fish

Understanding gene expression changes over the lifespan of cells is of fundamental interest and gives important insights into processes related to maturation and aging. This study was undertaken to understand the global transcriptome changes associated with aging in fish erythrocytes. Fish erythrocytes retain their nuclei throughout their lifetime and they are transcriptionally and translationally active. However, they lose important functions during their lifespan in the circulation. We separated rainbow trout (Oncorhynchus mykiss) erythrocytes into young and old fractions using fixed angle-centrifugation and analyzed transcriptome changes using RNA sequencing (RNA-seq) technology and quantitative real-time PCR. We found 930 differentially expressed between young and old erythrocyte fractions; 889 of these showed higher transcript levels in young, while only 34 protein-coding genes had higher transcript levels in old erythrocytes. In particular genes involved in ion binding, signal transduction, membrane transport, and those encoding various enzyme classes are affected in old erythrocytes. The transcripts with higher levels in old erythrocytes were associated with seven different GO terms within biological processes and nine within molecular functions and cellular components, respectively. Our study furthermore found several highly abundant transcripts as well as a number of differentially expressed genes (DEGs) for which the protein products are currently not known revealing the gaps of knowledge in most non-mammalian vertebrates. Our data provide the first insight into changes involved in aging on the transcriptional level and thus opens new perspectives for the study of maturation processes in fish erythrocytes.

The Pattern of microRNA Binding Site Distribution

Micro-RNA (miRNA or miR) regulates at least 60% of the genes in the human genome through their target sites at mRNA 3'-untranslated regions (UTR), and defects in miRNA expression regulation and target sites are frequently observed in cancers. We report here a systematic analysis of the distribution of miRNA target sites. Using the evolutionarily conserved miRNA binding sites in the TargetScan database (release 7.1), we constructed a miRNA co-regulation network by connecting genes sharing common miRNA target sites. The network possesses characteristics of the ubiquitous small-world network. Non-hub genes in the network-those sharing miRNA target sites with small numbers of genes-tend to form small cliques with their neighboring genes, while hub genes exhibit high levels of promiscuousness in their neighboring genes. Additionally, miRNA target site distribution is extremely uneven. Among the miRNAs, the distribution concentrates on a small number of miRNAs, in that their target sites occur in an extraordinarily large number of genes, that is, they have large numbers of target genes. The distribution across the genes follows a similar pattern; the mRNAs of a small proportion of the genes contain extraordinarily large numbers of miRNA binding sites. Quantitatively, the patterns fit into the P_(K) ∝ K^-α relationship (P_(K): the number of miRNAs with K target genes or genes with K miRNA sites; α: a positive constant), the mathematical description of connection distribution among the nodes and a defining characteristic of the so-called scale-free networks-a subset of small-world networks. Notably, well-known tumor-suppressive miRNAs (Let-7, miR-15/16, 26, 29, 31, 34, 145, 200, 203-205, 223, and 375) collectively have more than expected target genes, and well-known cancer genes contain more than expected miRNA binding sites. In summary, miRNA target site distribution exhibits characteristics of the small-world network. The potential to use this pattern to better understand miRNA function and their oncological roles is discussed.

Integrated Genomics Reveals Convergent Transcriptomic Networks Underlying Chronic Obstructive Pulmonary Disease and Idiopathic Pulmonary Fibrosis

RATIONALE

Despite shared environmental exposures, idiopathic pulmonary fibrosis (IPF) and chronic obstructive pulmonary disease are usually studied in isolation, and the presence of shared molecular mechanisms is unknown.

OBJECTIVES

We applied an integrative genomic approach to identify convergent transcriptomic pathways in emphysema and IPF.

METHODS

We defined the transcriptional repertoire of chronic obstructive pulmonary disease, IPF, or normal histology lungs using RNA-seq (n = 87).

MEASUREMENTS AND MAIN RESULTS

Genes increased in both emphysema and IPF relative to control were enriched for the p53/hypoxia pathway, a finding confirmed in an independent cohort using both gene expression arrays and the nCounter Analysis System (n = 193). Immunohistochemistry confirmed overexpression of HIF1A, MDM2, and NFKBIB members of this pathway in tissues from patients with emphysema or IPF. Using reads aligned across splice junctions, we determined that alternative splicing of p53/hypoxia pathway-associated molecules NUMB and PDGFA occurred more frequently in IPF or emphysema compared with control and validated these findings by quantitative polymerase chain reaction and the nCounter Analysis System on an independent sample set (n = 193). Finally, by integrating parallel microRNA and mRNA-Seq data on the same samples, we identified MIR96 as a key novel regulatory hub in the p53/hypoxia gene-expression network and confirmed that modulation of MIR96 in vitro recapitulates the disease-associated gene-expression network.

CONCLUSIONS

Our results suggest convergent transcriptional regulatory hubs in diseases as varied phenotypically as chronic obstructive pulmonary disease and IPF and suggest that these hubs may represent shared key responses of the lung to environmental stresses.

A SNP in the HTT promoter alters NF-κB binding and is a bidirectional genetic modifier of Huntington disease

Cis-regulatory variants that alter gene expression can modify disease expressivity, but none have previously been identified in Huntington disease (HD). Here we provide in vivo evidence in HD patients that cis-regulatory variants in the HTT promoter are bidirectional modifiers of HD age of onset. HTT promoter analysis identified a NF-κB binding site that regulates HTT promoter transcriptional activity. A non-coding SNP, rs13102260:G > A, in this binding site impaired NF-κB binding and reduced HTT transcriptional activity and HTT protein expression. The presence of the rs13102260 minor (A) variant on the HD disease allele was associated with delayed age of onset in familial cases, whereas the presence of the rs13102260 (A) variant on the wild-type HTT allele was associated with earlier age of onset in HD patients in an extreme case-based cohort. Our findings suggest a previously unknown mechanism linking allele-specific effects of rs13102260 on HTT expression to HD age of onset and have implications for HTT silencing treatments that are currently in development.

Dynamic network of transcription and pathway crosstalk to reveal molecular mechanism of MGd-treated human lung cancer cells

Recent research has revealed various molecular markers in lung cancer. However, the organizational principles underlying their genetic regulatory networks still await investigation. Here we performed Network Component Analysis (NCA) and Pathway Crosstalk Analysis (PCA) to construct a regulatory network in human lung cancer (A549) cells which were treated with 50 uM motexafin gadolinium (MGd), a metal cation-containing chemotherapeutic drug for 4, 12, and 24 hours. We identified a set of key TFs, known target genes for these TFs, and signaling pathways involved in regulatory networks. Our work showed that putative interactions between these TFs (such as ESR1/Sp1, E2F1/Sp1, c-MYC-ESR, Smad3/c-Myc, and NFKB1/RELA), between TFs and their target genes (such as BMP41/Est1, TSC2/Myc, APE1/Sp1/p53, RARA/HOXA1, and SP1/USF2), and between signaling pathways (such as PPAR signaling pathway and Adipocytokines signaling pathway). These results will provide insights into the regulatory mechanism of MGd-treated human lung cancer cells.

An update on the strategies in multicomponent activity monitoring within the phytopharmaceutical field

Background

To-date modern drug research has focused on the discovery and synthesis of single active substances. However, multicomponent preparations are gaining increasing importance in the phytopharmaceutical field by demonstrating beneficial properties with respect to efficacy and toxicity.

Discussion

In contrast to single drug combinations, a botanical multicomponent therapeutic possesses a complex repertoire of chemicals that belong to a variety of substance classes. This may explain the frequently observed pleiotropic bioactivity spectra of these compounds, which may also suggest that they possess novel therapeutic opportunities. Interestingly, considerable bioactivity properties are exhibited not only by remedies that contain high doses of phytochemicals with prominent pharmaceutical efficacy, but also preparations that lack a sole active principle component. Despite that each individual substance within these multicomponents has a low molar fraction, the therapeutic activity of these substances is established via a potentialization of their effects through combined and simultaneous attacks on multiple molecular targets. Although beneficial properties may emerge from such a broad range of perturbations on cellular machinery, validation and/or prediction of their activity profiles is accompanied with a variety of difficulties in generic risk-benefit assessments. Thus, it is recommended that a comprehensive strategy is implemented to cover the entirety of multicomponent-multitarget effects, so as to address the limitations of conventional approaches.

Summary

An integration of standard toxicological methods with selected pathway-focused bioassays and unbiased data acquisition strategies (such as gene expression analysis) would be advantageous in building an interaction network model to consider all of the effects, whether they were intended or adverse reactions.

Genome-wide studies of mRNA synthesis and degradation in eukaryotes

In recent years, the use of genome-wide technologies has revolutionized the study of eukaryotic transcription producing results for thousands of genes at every step of mRNA life. The statistical analyses of the results for a single condition, different conditions, different transcription stages, or even between different techniques, is outlining a totally new landscape of the eukaryotic transcription process. Although most studies have been conducted in the yeast Saccharomyces cerevisiae as a model cell, others have also focused on higher eukaryotes, which can also be comparatively analyzed. The picture which emerges is that transcription is a more variable process than initially suspected, with large differences between genes at each stage of the process, from initiation to mRNA degradation, but with striking similarities for functionally related genes, indicating that all steps are coordinately regulated. This article is part of a Special Issue entitled: Nuclear Transport and RNA Processing.

Genomic insights into the different layers of gene regulation in yeast

The model organism Saccharomyces cerevisiae has allowed the development of new functional genomics techniques devoted to the study of transcription in all its stages. With these techniques, it has been possible to find interesting new mechanisms to control gene expression that act at different levels and for different gene sets apart from the known cis-trans regulation in the transcription initiation step. Here we discuss a method developed in our laboratory, Genomic Run-On, and other new methods that have recently appeared with distinct technical features. A comparison between the datasets generated by them provides interesting genomic insights into the different layers of gene regulation in eukaryotes.

Heat shock response in yeast involves changes in both transcription rates and mRNA stabilities

We have analyzed the heat stress response in the yeast Saccharomyces cerevisiae by determining mRNA levels and transcription rates for the whole transcriptome after a shift from 25 °C to 37 °C. Using an established mathematical algorithm, theoretical mRNA decay rates have also been calculated from the experimental data. We have verified the mathematical predictions for selected genes by determining their mRNA decay rates at different times during heat stress response using the regulatable tetO promoter. This study indicates that the yeast response to heat shock is not only due to changes in transcription rates, but also to changes in the mRNA stabilities. mRNA stability is affected in 62% of the yeast genes and it is particularly important in shaping the mRNA profile of the genes belonging to the environmental stress response. In most cases, changes in transcription rates and mRNA stabilities are homodirectional for both parameters, although some interesting cases of antagonist behavior are found. The statistical analysis of gene targets and sequence motifs within the clusters of genes with similar behaviors shows that both transcriptional and post-transcriptional regulons apparently contribute to the general heat stress response by means of transcriptional factors and RNA binding proteins.