Coupling and coordination in gene expression processes: a systems biology view

Nanodynamo quantifies subcellular RNA dynamics revealing extensive coupling between steps of the RNA life cycle

The coordinated action of transcriptional and post-transcriptional machineries shapes gene expression programs at steady state and determines their concerted response to perturbations. We have developed Nanodynamo, an experimental and computational workflow for quantifying the kinetic rates of nuclear and cytoplasmic steps of the RNA life cycle. Nanodynamo is based on mathematical modelling following sequencing of native RNA from cellular fractions and polysomes. We have applied this workflow to triple-negative breast cancer cells, revealing widespread post-transcriptional RNA processing that is mutually exclusive with its co-transcriptional counterpart. We used Nanodynamo to unravel the coupling between transcription, processing, export, decay and translation machineries. We have identified a number of coupling interactions within and between the nucleus and cytoplasm that largely contribute to coordinating how cells respond to perturbations that affect gene expression programs. Nanodynamo will be instrumental in unravelling the determinants and regulatory processes involved in the coordination of gene expression responses.

Predicting splicing patterns from the transcription factor binding sites in the promoter with deep learning

BACKGROUND

Alternative splicing is a pivotal mechanism of post-transcriptional modification that contributes to the transcriptome plasticity and proteome diversity in metazoan cells. Although many splicing regulations around the exon/intron regions are known, the relationship between promoter-bound transcription factors and the downstream alternative splicing largely remains unexplored.

RESULTS

In this study, we present computational approaches to unravel the regulatory relationship between promoter-bound transcription factor binding sites (TFBSs) and the splicing patterns. We curated a fine dataset that includes DNase I hypersensitive site sequencing and transcriptomes across fifteen human tissues from ENCODE. Specifically, we proposed different representations of TF binding context and splicing patterns to examine the associations between the promoter and downstream splicing events. While machine learning models demonstrated potential in predicting splicing patterns based on TFBS occupancies, the limitations in the generalization of predicting the splicing forms of singleton genes across diverse tissues was observed with carefully examination using different cross-validation methods. We further investigated the association between alterations in individual TFBS at promoters and shifts in exon splicing efficiency. Our results demonstrate that the convolutional neural network (CNN) models, trained on TF binding changes in the promoters, can predict the changes in splicing patterns. Furthermore, a systemic in silico substitutions analysis on the CNN models highlighted several potential splicing regulators. Notably, using empirical validation using K562 CTCFL shRNA knock-down data, we showed the significant role of CTCFL in splicing regulation.

CONCLUSION

In conclusion, our finding highlights the potential role of promoter-bound TFBSs in influencing the regulation of downstream splicing patterns and provides insights for discovering alternative splicing regulations.

Gene regulatory networks in disease and ageing

The precise control of gene expression is required for the maintenance of cellular homeostasis and proper cellular function, and the declining control of gene expression with age is considered a major contributor to age-associated changes in cellular physiology and disease. The coordination of gene expression can be represented through models of the molecular interactions that govern gene expression levels, so-called gene regulatory networks. Gene regulatory networks can represent interactions that occur through signal transduction, those that involve regulatory transcription factors, or statistical models of gene-gene relationships based on the premise that certain sets of genes tend to be coexpressed across a range of conditions and cell types. Advances in experimental and computational technologies have enabled the inference of these networks on an unprecedented scale and at unprecedented precision. Here, we delineate different types of gene regulatory networks and their cell-biological interpretation. We describe methods for inferring such networks from large-scale, multi-omics datasets and present applications that have aided our understanding of cellular ageing and disease mechanisms.

LowTempGAL: a highly responsive low temperature-inducible GAL system in Saccharomyces cerevisiae

Temperature is an important control factor for biologics biomanufacturing in precision fermentation. Here, we explored a highly responsive low temperature-inducible genetic system (LowTempGAL) in the model yeast Saccharomyces cerevisiae. Two temperature biosensors, a heat-inducible degron and a heat-inducible protein aggregation domain, were used to regulate the GAL activator Gal4p, rendering the leaky LowTempGAL systems. Boolean-type induction was achieved by implementing a second-layer control through low-temperature-mediated repression on GAL repressor gene GAL80, but suffered delayed response to low-temperature triggers and a weak response at 30°C. Application potentials were validated for protein and small molecule production. Proteomics analysis suggested that residual Gal80p and Gal4p insufficiency caused suboptimal induction. 'Turbo' mechanisms were engineered through incorporating a basal Gal4p expression and a galactose-independent Gal80p-supressing Gal3p mutant (Gal3Cp). Varying Gal3Cp configurations, we deployed the LowTempGAL systems capable for a rapid stringent high-level induction upon the shift from a high temperature (37-33°C) to a low temperature (≤30°C). Overall, we present a synthetic biology procedure that leverages 'leaky' biosensors to deploy highly responsive Boolean-type genetic circuits. The key lies in optimisation of the intricate layout of the multi-factor system. The LowTempGAL systems may be applicable in non-conventional yeast platforms for precision biomanufacturing.

Directional integration and pathway enrichment analysis for multi-omics data

Omics techniques generate comprehensive profiles of biomolecules in cells and tissues. However, a holistic understanding of underlying systems requires joint analyses of multiple data modalities. We present DPM, a data fusion method for integrating omics datasets using directionality and significance estimates of genes, transcripts, or proteins. DPM allows users to define how the input datasets are expected to interact directionally given the experimental design or biological relationships between the datasets. DPM prioritises genes and pathways that change consistently across the datasets and penalises those with inconsistent directionality. To demonstrate our approach, we characterise gene and pathway regulation in IDH-mutant gliomas by jointly analysing transcriptomic, proteomic, and DNA methylation datasets. Directional integration of survival information in ovarian cancer reveals candidate biomarkers with consistent prognostic signals in transcript and protein expression. DPM is a general and adaptable framework for gene prioritisation and pathway analysis in multi-omics datasets.

Norovirus NS1/2 protein increases glutaminolysis for efficient viral replication

Viruses are obligate intracellular parasites that rely on host cell metabolism for successful replication. Thus, viruses rewire host cell pathways involved in central carbon metabolism to increase the availability of building blocks for successful propagation. However, the underlying mechanisms of virus-induced alterations to host metabolism are largely unknown. Noroviruses (NoVs) are highly prevalent pathogens that cause sporadic and epidemic viral gastroenteritis. In the present study, we uncovered several strain-specific and shared host cell metabolic requirements of three murine norovirus (MNV) strains, MNV-1, CR3, and CR6. While all three strains required glycolysis, glutaminolysis, and the pentose phosphate pathway for optimal infection of macrophages, only MNV-1 relied on host oxidative phosphorylation. Furthermore, the first metabolic flux analysis of NoV-infected cells revealed that both glycolysis and glutaminolysis are upregulated during MNV-1 infection of macrophages. Glutamine deprivation affected the viral lifecycle at the stage of genome replication, resulting in decreased non-structural and structural protein synthesis, viral assembly, and egress. Mechanistic studies further showed that MNV infection and overexpression of the non-structural protein NS1/2 increased the enzymatic activity of the rate-limiting enzyme glutaminase. In conclusion, the inaugural investigation of NoV-induced alterations to host glutaminolysis identified NS1/2 as the first viral molecule for RNA viruses that regulates glutaminolysis either directly or indirectly. This increases our fundamental understanding of virus-induced metabolic alterations and may lead to improvements in the cultivation of human NoVs.

Noncontiguous operon atlas for the Staphylococcus aureus genome

Bacteria synchronize the expression of genes with related functions by organizing genes into operons so that they are cotranscribed together in a single polycistronic messenger RNA. However, some cellular processes may benefit if the simultaneous production of the operon proteins coincides with the inhibition of the expression of an antagonist gene. To coordinate such situations, bacteria have evolved noncontiguous operons (NcOs), a subtype of operons that contain one or more genes that are transcribed in the opposite direction to the other operon genes. This structure results in overlapping transcripts whose expression is mutually repressed. The presence of NcOs cannot be predicted computationally and their identification requires a detailed knowledge of the bacterial transcriptome. In this study, we used direct RNA sequencing methodology to determine the NcOs map in the Staphylococcus aureus genome. We detected the presence of 18 NcOs in the genome of S. aureus and four in the genome of the lysogenic prophage 80α. The identified NcOs comprise genes involved in energy metabolism, metal acquisition and transport, toxin-antitoxin systems, and control of the phage life cycle. Using the menaquinone operon as a proof of concept, we show that disarrangement of the NcO architecture results in a reduction of bacterial fitness due to an increase in menaquinone levels and a decrease in the rate of oxygen consumption. Our study demonstrates the significance of NcO structures in bacterial physiology and emphasizes the importance of combining operon maps with transcriptomic data to uncover previously unnoticed functional relationships between neighbouring genes.

Protein interaction network revealed by quantitative proteomic analysis links TFIIB to multiple aspects of the transcription cycle

Although TFIIB is widely regarded as an initiation factor, recent reports have implicated it in multiple aspects of eukaryotic transcription. To investigate the broader role of TFIIB in transcription, we performed quantitative proteomic analysis of yeast TFIIB. We purified two different populations of TFIIB; one from soluble cell lysate, which is not engaged in transcription, and the other from the chromatin fraction which yields the transcriptionally active form of the protein. TFIIB purified from the chromatin exhibits several interactions that explain its non-canonical roles in transcription. RNAPII, TFIIF and TFIIH were the only components of the preinitiation complex with a significant presence in chromatin TFIIB. A notable feature was enrichment of all subunits of CF1 and Rat1 3' end processing-termination complexes in chromatin-TFIIB preparation. Subunits of the CPF termination complex were also detected in both chromatin and soluble derived TFIIB preparations. These results may explain the presence of TFIIB at the 3' end of genes during transcription as well as its role in promoter-termination interaction.

Customization of the translational complex regulates mRNA-specific translation to control CNS regeneration

In the adult mammalian central nervous system (CNS), axons fail to regenerate spontaneously after injury because of a combination of extrinsic and intrinsic factors. Despite recent advances targeting the intrinsic regenerative properties of adult neurons, the molecular mechanisms underlying axon regeneration are not fully understood. Here, we uncover a regulatory mechanism that controls the expression of key proteins involved in regeneration at the translational level. Our results show that mRNA-specific translation is critical for promoting axon regeneration. Indeed, we demonstrate that specific ribosome-interacting proteins, such as the protein Huntingtin (HTT), selectively control the translation of a specific subset of mRNAs. Moreover, modulating the expression of these translationally regulated mRNAs is crucial for promoting axon regeneration. Altogether, our findings highlight that selective translation through the customization of the translational complex is a key mechanism of axon regeneration with major implications in the development of therapeutic strategies for CNS repair.

In Vitro Expression Analysis Reveals HML6-c14 to Be an Attractive Research Target

HML6-c14, a long terminal repeat (LTR)-type retrotransposon identified by expressed sequence tag (EST) database screening, was found to undergo RNA processing resembling that of placental tissue by in vitro expression analysis. Previous in situ hybridization studies using normal placental tissue showed that the transcript remained in the nucleus. However, among the transcripts forcedly expressed in cultured cells, the transcript that retained the 3.3 kb intron was observed in the nucleus, and a part of the spliced transcript was observed outside the nucleus. To verify whether this cytoplasmic transcript could be translated, we examined the coding potential of the open reading frame (ORF), consisting of 109 codons on the spliced transcript, along with two other putative ORFs detected in the intronic region. As a result, none of the ORF-derived products could be detected by Western blotting as fusion proteins tagged with the FLAG epitope, suggesting that HML6-c14 belongs to a group of long non-coding RNA (lncRNA) genes. Promoter analysis of the upstream 6.4 kb genomic region also suggested that the 5'-LTR itself potentially retains high promoter activity. Despite losing the ability to produce functional proteins, HML6-c14 continues to retain its transcriptional ability while converting to an lncRNA gene, which is an interesting subject for future research.

Coexpression Network Construction and Visualization from Transcriptomes Underlying ER Stress Responses

Dynamic gene expression changes are primary cellular reactions in response to most stresses and developmental cues in all organisms, including plants. With the ever-decreasing cost and increasing access, high-throughput transcriptome analyses have become a significant research tool to understand a wide spectrum of complex gene regulatory mechanisms. However, it is still challenging to understand the complete picture of gene responses because of the interactive and dynamic nature of gene expression in biological networks. Coexpression network analyses followed by network mapping are being increasingly applied to overcome this challenge. In this chapter, we will introduce detailed instructions for performing a weighted coexpression network analysis (WGCNA) and network visualization using a transcriptome dataset obtained during recovery from endoplasmic reticulum (ER) stress in Arabidopsis thaliana. The streamlined workflow described here allows biologists to identify and visualize coexpression interactions among genes, accessing a comprehensive landscape of dynamic gene expression changes for further downstream analyses using their datasets.

PUS1 May Be a Potential Prognostic Biomarker and Therapeutic Target for Hepatocellular Carcinoma

Objective

The mechanisms of pseudouridine synthase (PUS) are not definite in hepatocellular carcinoma (HCC), the objective of this study is to investigate the effect of PUS genes in HCC.

Materials and Methods

Differentially expressed and prognostic gene of PUS enzymes was identified based on The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC) and Gene Expression Profiling Interactive Analysis (GEPIA) databases. For the identified gene, pseudouridine synthase 1 (PUS1), was used for further research. The clinicopathological feature of PUS1 was analyzed by Student's t-test. Prognostic significance was explored by Kaplan-Meier (KM) analysis and Cox proportional hazards regression model. Receiver operating characteristic (ROC) curve was applied to appraise diagnostic and prognostic value. The Database for Annotation, Visualization, and Integrated Discovery (DAVID) and Gene Set Enrichment Analysis (GSEA) were implemented to explore mechanism of PUS1. A Guangxi cohort was applied to verify differential expression. In vitro cell experiments were implemented to investigate the influence for proliferation, reactive oxygen species (ROS) level, migration, and invasion of HCC cells after a knockdown of PUS1.

Results

PUS1 was significantly overexpressed in HCC tissues, and patients with high PUS1 were related to unpromising clinicopathological features. Survival analysis revealed high PUS1 expression was associated with a poor overall survival (OS) and 1 year-recurrence free survival (RFS), was an independent risk factor. Meanwhile, ROC curve showed that PUS1 had a diagnostic and prognostic significance to HCC. Functional enrichment analysis implied that PUS1 may be involved in metabolic pathways, mitochondrial function, non-alcoholic fatty liver disease (NAFLD), and some important carcinogenic pathways. Cell assays revealed that knockdown of PUS1 significantly constrained the migration, proliferation, invasion and improved the ROS level of HCC cells.

Conclusion

PUS1 may be a prognostic biomarker and a underlying treatment target for HCC.

scGENA: A Single-Cell Gene Coexpression Network Analysis Framework for Clustering Cell Types and Revealing Biological Mechanisms

Single-cell RNA-sequencing (scRNA-seq) is a recent high-throughput technique that can measure gene expression, reveal cell heterogeneity, rare and complex cell populations, and discover cell types and their relationships. The analysis of scRNA-seq data is challenging because of transcripts sparsity, replication noise, and outlier cell populations. A gene coexpression network (GCN) analysis effectively deciphers phenotypic differences in specific states by describing gene–gene pairwise relationships. The underlying gene modules with different coexpression patterns partially bridge the gap between genotype and phenotype. This study presents a new framework called scGENA (single-cell gene coexpression network analysis) for GCN analysis based on scRNA-seq data. Although there are several methods for scRNA-seq data analysis, we aim to build an integrative pipeline for several purposes that cover primary data preprocessing, including data exploration, quality control, normalization, imputation, and dimensionality reduction of clustering as downstream of GCN analysis. To demonstrate this integrated workflow, an scRNA-seq dataset of the human diabetic pancreas with 1600 cells and 39,851 genes was implemented to perform all these processes in practice. As a result, scGENA is demonstrated to uncover interesting gene modules behind complex diseases, which reveal biological mechanisms. scGENA provides a state-of-the-art method for gene coexpression analysis for scRNA-seq data.

Transcriptional competition shapes proteotoxic ER stress resolution

Through dynamic activities of conserved master transcription factors (mTFs), the unfolded protein response (UPR) relieves proteostasis imbalance of the endoplasmic reticulum (ER), a condition known as ER stress^1,2. Because dysregulated UPR is lethal, the competence for fate changes of the UPR mTFs must be tightly controlled^3,4. However, the molecular mechanisms underlying regulatory dynamics of mTFs remain largely elusive. Here, we identified the abscisic acid-related regulator G-class bZIP TF2 (GBF2) and the cis-regulatory element G-box as regulatory components of the plant UPR led by the mTFs, bZIP28 and bZIP60. We demonstrate that, by competing with the mTFs at G-box, GBF2 represses UPR gene expression. Conversely, a gbf2 null mutation enhances UPR gene expression and suppresses the lethality of a bzip28 bzip60 mutant in unresolved ER stress. By demonstrating that GBF2 functions as a transcriptional repressor of the UPR, we address the long-standing challenge of identifying shared signalling components for a better understanding of the dynamic nature and complexity of stress biology. Furthermore, our results identify a new layer of UPR gene regulation hinged upon an antagonistic mTFs-GFB2 competition for proteostasis and cell fate determination.

Modulation of RNA stability regulates gene expression in two opposite ways: through buffering of RNA levels upon global perturbations and by supporting adapted differential expression

The steady state levels of RNAs, often referred to as expression levels, result from a well-balanced combination of RNA transcription and decay. Alterations in RNA levels will therefore result from tight regulation of transcription rates, decay rates or both. Here, we explore the role of RNA stability in achieving balanced gene expression and present genome-wide RNA stabilities in Drosophila melanogaster male and female cells as well as male cells depleted of proteins essential for dosage compensation. We identify two distinct RNA-stability mediated responses involved in regulation of gene expression. The first of these responds to acute and global changes in transcription and thus counteracts potentially harmful gene mis-expression by shifting the RNA stability in the direction opposite to the transcriptional change. The second response enhances inter-individual differential gene expression by adjusting the RNA stability in the same direction as a transcriptional change. Both mechanisms are global, act on housekeeping as well as non-housekeeping genes and were observed in both flies and mammals. Additionally, we show that, in contrast to mammals, modulation of RNA stability does not detectably contribute to dosage compensation of the sex-chromosomes in D. melanogaster.

The Unfolded Protein Response at the Tumor-Immune Interface

The tumor-immune interface has surged to primary relevance in an effort to understand the hurdles facing immune surveillance and cancer immunotherapy. Reports over the past decades have indicated a role for the unfolded protein response (UPR) in modulating not only tumor cell fitness and drug resistance, but also local immunity, with emphasis on the phenotype and altered function of immune cells such as myeloid cells and T cells. Emerging evidence also suggests that aneuploidy correlates with local immune dysregulation. Recently, we reported that the UPR serves as a link between aneuploidy and immune cell dysregulation in a cell nonautonomous way. These new findings add considerable complexity to the organization of the tumor microenvironment (TME) and the origin of its altered function. In this review, we summarize these data and also discuss the role of aneuploidy as a negative regulator of local immunity.

Human transcription factor protein interaction networks

Transcription factors (TFs) interact with several other proteins in the process of transcriptional regulation. Here, we identify 6703 and 1536 protein–protein interactions for 109 different human TFs through proximity-dependent biotinylation (BioID) and affinity purification mass spectrometry (AP-MS), respectively. The BioID analysis identifies more high-confidence interactions, highlighting the transient and dynamic nature of many of the TF interactions. By performing clustering and correlation analyses, we identify subgroups of TFs associated with specific biological functions, such as RNA splicing or chromatin remodeling. We also observe 202 TF-TF interactions, of which 118 are interactions with nuclear factor 1 (NFI) family members, indicating uncharacterized cross-talk between NFI signaling and other TF signaling pathways. Moreover, TF interactions with basal transcription machinery are mainly observed through TFIID and SAGA complexes. This study provides a rich resource of human TF interactions and also act as a starting point for future studies aimed at understanding TF-mediated transcription.

Transcription factors (TFs) interact with several other proteins in the process of transcriptional regulation. Here the authors identify 6703 and 1536 protein–protein interactions for 109 different human TFs through BioID and AP-MS analyses, respectively.

Anota2seq Analysis for Transcriptome-Wide Studies of mRNA Translation

mRNA translation plays a critical role in determining proteome composition. In health, regulation of mRNA translation facilitates rapid gene expression responses to intra- and extracellular signals. Moreover, dysregulated mRNA translation is a common feature in disease states, including neurological disorders and cancer. Yet, most studies of gene expression focus on analysis of mRNA levels, leaving variations in translational efficiencies largely uncharacterized. Here, we outline procedures to identify mRNA-selective alterations in translational efficiencies on a transcriptome-wide scale using the anota2seq package. Anota2seq compares expression data originating from translated mRNA to data from matched total mRNA to identify changes in translated mRNA not paralleled by corresponding changes in total mRNA (interpreted as changes in translational efficiencies impacting protein levels), congruent changes in total and translated mRNA (interpreted as changes in transcription and/or mRNA stability), and changes in total mRNA not paralleled by corresponding alterations in translated mRNA (interpreted as translational buffering). To illustrate the functionality of the anota2seq analysis package, we demonstrate a detailed analysis using a polysome-profiling data set quantified by RNA sequencing, revealing that estrogen receptor α modulates gene expression via a type of translational buffering termed offsetting. Notably, this anota2seq analysis procedure is also applicable to ribosome-profiling (RiboSeq) data sets and can be adapted to a variety of other data types and experimental contexts. Finally, we provide guidance for extending anota2seq analysis to examine associations between untranslated regions and altered translational efficiencies as well as targeted cellular functions to gain insights into mechanisms and phenotypic consequences of altered mRNA translation. Thus, this step-by-step manual allows users to interrogate selective changes in mRNA translation on a transcriptome-wide scale using the Bioconductor package anota2seq.

The unfolded protein response links tumor aneuploidy to local immune dysregulation

Aneuploidy is a chromosomal abnormality associated with poor prognosis in many cancer types. Here, we tested the hypothesis that the unfolded protein response (UPR) mechanistically links aneuploidy and local immune dysregulation. Using a single somatic copy number alteration (SCNA) score inclusive of whole‐chromosome, chromosome arm, and focal alterations in a pan‐cancer analysis of 9,375 samples in The Cancer Genome Atlas (TCGA) database, we found an inverse correlation with a cytotoxicity (CYT) score across disease stages. Co‐expression patterns of UPR genes changed substantially between SCNA^low and SCNA^high groups. Pathway activity scores showed increased activity of multiple branches of the UPR in response to aneuploidy. The PERK branch showed the strongest association with a reduced CYT score. The conditioned medium of aneuploid cells transmitted XBP1 splicing and caused IL‐6 and arginase 1 transcription in receiver bone marrow‐derived macrophages and markedly diminished the production of IFN‐γ and granzyme B in activated human T cells. We propose the UPR as a mechanistic link between aneuploidy and immune dysregulation in the tumor microenvironment.

Combined analyses of translatome and transcriptome in Arabidopsis reveal new players responding to magnesium deficiency

Translational control of gene expression, including recruitment of ribosomes to messenger RNA (mRNA), is particularly important during the response to stress. Purification of ribosome-associated mRNAs using translating ribosome affinity purification (TRAP) followed by RNA-sequencing facilitates the study of mRNAs undergoing active transcription and better proxies the translatome, or protein response, to stimuli. To identify plant responses to Magnesium (Mg) deficiency at the translational level, we combined transcriptome and translatome analyses. Excitingly, we found 26 previously unreported Mg-responsive genes that were only regulated at the translational level and not the transcriptional level, during the early response to Mg deficiency. In addition, mutants of the transcription factor ELONGATED HYPOCOTYL 5 (HY5), the H⁺ /CATION EXCHANGER 1 and 3 (CAX1 and CAX3), and UBIQUITIN 11 (UBQ11) exhibited early chlorosis phenotype under Mg deficiency, supporting their functional involvement in ion homeostasis. Overall, our study strongly supports that TRAP-seq combined with RNA-seq followed by phenotype screening could facilitate the identification of novel players during stress responses.

Inhibition of MAVS Aggregation-Mediated Type-I Interferon Signaling by Foot-and-Mouth Disease Virus VP3

As a structural protein of the Foot-and-mouth disease virus (FMDV), VP3 plays a vital role in virus assembly and inhibiting the interferon (IFN) signal transduction to promote FMDV replication. Previous studies demonstrated that FMDV VP3 blocks the type-I IFN response by inhibiting the mRNA expression of the mitochondrial antiviral-signaling protein (MAVS); however, the underlying mechanism is poorly understood. Here, we describe the specificity of FMDV VP3 interaction with the transmembrane (TM) domain of MAVS as FMDV driven type-I IFN inhibitory mechanism for its effective replication. The TM domain of MAVS governs the mitochondria localization of MAVS, and it is a key factor in type-I IFN signaling transduction via MAVS aggregation. Thereby, the interaction of FMDV VP3 with the TM domain of MAVS leads to the inhibition of MAVS mitochondria localization, self-association, and aggregation, resulting in the suppression of type-I IFN response. Collectively, these results provide a clear understanding of a key molecular mechanism used by the FMDV VP3 for the suppression of IFN responses via targeting MAVS.

Pbp1, the yeast ortholog of human Ataxin-2, functions in the cell growth on non-fermentable carbon sources

Pbp1, the yeast ortholog of human Ataxin-2, was originally isolated as a poly(A) binding protein (Pab1)-binding protein. Pbp1 regulates the Pan2-Pan3 deadenylase complex, thereby modulating the mRNA stability and translation efficiency. However, the physiological significance of Pbp1 remains unclear since a yeast strain harboring PBP1 deletion grows similarly to wild-type strain on normal glucose-containing medium. In this study, we found that Pbp1 has a role in cell growth on the medium containing non-fermentable carbon sources. While the pbp1Δ mutant showed a similar growth compared to the wild-type cell on a normal glucose-containing medium, the pbp1Δ mutant showed a slower growth on the medium containing glycerol and lactate. Microarray analyses revealed that expressions of the genes involved in gluconeogenesis, such as PCK1 and FBP1, and of the genes involved in mitochondrial function, such as COX10 and COX11, were decreased in the pbp1Δ mutant. Pbp1 regulated the expressions of PCK1 and FBP1 via their promoters, while the expressions of COX10 and COX11 were regulated by Pbp1, not through their promoters. The decreased expressions of COX10 and COX11 in the pbp1Δ mutant were recovered by the loss of Dcp1 decapping enzyme or Xrn1 5’-3’exonuclease. Our results suggest that Pbp1 regulates the expressions of the genes involved in gluconeogenesis and mitochondrial function through multiple mechanisms.

Spatio-Temporal Coupling Coordination Analysis between Urbanization and Water Resource Carrying Capacity of the Provinces in the Yellow River Basin, China

With the rapid expansion of the Chinese economy in recent years, the urbanization process of the provinces in the Yellow River Basin (YRB) has put serious pressure on the sustainability of the water resources carrying capacity (WRCC). It is necessary to analyze and diagnose the coordination state between urbanization and the WRCC. In this study, based on the Population-Economic-Social-Spatial (PESS) framework and Pressure-State-Response (PSR) model, we developed two index systems for the urbanization and WRCC, respectively. At the basis of the two index systems, the coupling coordination degree (CCD) of the two systems is calculated by applying the improved CCD model. Based on the calculated CCD for each province, the spatio-temporal analysis was performed to analyze the characteristics of CCD in the YRB. The obstacle factor model is utilized to obtain the main obstacle factors. The results show that: (1) the coordination state between the urbanization and WRCC systems was improved to some extent in 2017, compared to 2008, but there are differences in the coordination state of the different provinces in the YRB. (2) In terms of the level of urbanization, the gap between the seven provinces’ performance levels widened because urbanization grew at different rates. The WRCC system’s performance presented a fluctuating downward trend from 2008 to 2017 in the YRB. (3) The pressure subsystem had a significant impact on the two systems’ coordination state in the YRB, while the social urbanization and response subsystem had a less significant impact on the urbanization system and the WRCC system, respectively. Due to the growth of urbanization, the imbalanced development of the WRCC and urbanization has become the principal contradiction that must be solved in order to achieve sustainability in the YRB. The analysis of the coupling relationship between urbanization and WRCC may guide the policy makers in planning for realistic goals. The results provide a guide for high-quality development in the YRB.

Transcriptomic Changes in Endothelial Cells Triggered by Na,K-ATPase Inhibition: A Search for Upstream Na+i/K+i Sensitive Genes

Stimulus-dependent elevation of intracellular Ca²⁺ affects gene expression via well-documented calmodulin-mediated signaling pathways. Recently, we found that the addition of extra- and intracellular Ca²⁺ chelators increased, rather than decreased, the number of genes expressed, and that this is affected by the elevation of [Na⁺]_i/[K⁺]_i-ratio. This assumes the existence of a novel Na⁺_i/K⁺_i-mediated Ca²⁺_i-independent mechanism of excitation-transcription coupling. To identify upstream Na⁺_i/K⁺_i-sensitive genes, we examined the kinetics of transcriptomic changes in human umbilical vein endothelial cells (HUVEC) subjected to Na,K-ATPase inhibition by ouabain or K⁺-free medium. According to our data, microRNAs, transcription factors, and proteins involved in immune response and inflammation might be considered as key components of Na⁺_i/K⁺_i-mediated excitation-transcription coupling. Special attention was focused on the FOS gene and the possible mechanism of transcription regulation via G-quadruplexes, non-canonical secondary structures of nucleic acids, whose stability depends on [Na⁺]_i/[K⁺]_i-ratio. Verification of the [Na⁺]_i/[K⁺]_i-sensitive transcription regulation mechanism should be continued in forthcoming studies.

Assessing the coordination between economic growth and urban climate change in China from 2000 to 2015

The balance between economic growth and environmental protection has been a critical concern for sustainable urban development. However, among the multiple research efforts exploring the coordination between the two aspects, the widespread urban climate change has rarely been considered. This study encompasses urban climate change into the cross-system coupling analysis framework to assess its coordination with economic growth using the Coupling Coordination Degree (CCD) model. The two aspects are respectively represented using indicators of Surface Urban Heat Island Intensity (SUHII) and Gross Domestic Product (GDP). Specifically, China is selected as case study, and a total of 259 cities from the 2000-2015 period are analyzed. The spatio-temporal patterns of CCD are investigated through time series clustering to uncover its performance under diversified economic and climatic contexts. The regional inequality and spatial agglomeration effects are also examined. Results reveal significant spatio-temporal heterogeneity of CCD. Geographically, CCD varies from uncoordinated to high-level coordination. Wealthier cities in the eastern coastal region are significantly better coordinated than their inland counterparts. Temporally, the uptrend of CCD is not significant for most cities due to the relatively insufficient emphasis on urban heat island (UHI) mitigation in previous efforts. Evident spatial inequality and agglomeration patterns are also observed with slight downtrends. The spatio-temporal patterns of CCD revealed in this study indicate great necessity for the central government to develop policies suiting cities' special characteristics and contexts, and more efforts should be targeted on reducing regional imbalance. Hence, a nation-city-community policy skeleton is last outlined to enhance the pursuit of a more climate-friendly urban environment under rapid economic development. Overall, this study advances the understanding of economy-urban climate interactions and facilitates the future pursuit of better sustainable cities. The proposed workflow can be utilized for other countries with diversified urbanization processes and potentially used for comparison among different countries.

Proteogenomic Characterization Reveals Therapeutic Vulnerabilities in Lung Adenocarcinoma

To explore the biology of lung adenocarcinoma (LUAD) and identify new therapeutic opportunities, we performed comprehensive proteogenomic characterization of 110 tumors and 101 matched normal adjacent tissues (NATs) incorporating genomics, epigenomics, deep-scale proteomics, phosphoproteomics, and acetylproteomics. Multi-omics clustering revealed four subgroups defined by key driver mutations, country, and gender. Proteomic and phosphoproteomic data illuminated biology downstream of copy number aberrations, somatic mutations, and fusions and identified therapeutic vulnerabilities associated with driver events involving KRAS, EGFR, and ALK. Immune subtyping revealed a complex landscape, reinforced the association of STK11 with immune-cold behavior, and underscored a potential immunosuppressive role of neutrophil degranulation. Smoking-associated LUADs showed correlation with other environmental exposure signatures and a field effect in NATs. Matched NATs allowed identification of differentially expressed proteins with potential diagnostic and therapeutic utility. This proteogenomics dataset represents a unique public resource for researchers and clinicians seeking to better understand and treat lung adenocarcinomas.

Generally applicable transcriptome-wide analysis of translation using anota2seq

mRNA translation plays an evolutionarily conserved role in homeostasis and when dysregulated contributes to various disorders including metabolic and neurological diseases and cancer. Notwithstanding that optimal and universally applicable methods are critical for understanding the complex role of translational control under physiological and pathological conditions, approaches to analyze translatomes are largely underdeveloped. To address this, we developed the anota2seq algorithm which outperforms current methods for statistical identification of changes in translation. Notably, in contrast to available analytical methods, anota2seq also allows specific identification of an underappreciated mode of gene expression regulation whereby translation acts as a buffering mechanism which maintains protein levels despite fluctuations in corresponding mRNA abundance (‘translational buffering’). Thus, the universal anota2seq algorithm allows efficient and hitherto unprecedented interrogation of translatomes which is anticipated to advance knowledge regarding the role of translation in homeostasis and disease.

URH49 exports mRNA by remodeling complex formation and mediating the NXF1-dependent pathway

The TREX complex integrates information from nuclear mRNA processing events to ensure the timely export of mRNA to the cytoplasm. In humans, UAP56 and its paralog URH49 form distinct complexes, the TREX complex and the AREX complex, respectively, which cooperatively regulate the expression of a specific set of mRNA species on a genome wide scale. The difference in the complex formation between UAP56 and URH49 are thought to play a critical role in the regulation of target mRNAs. To date, the underlying mechanism remains poorly understood. Here we characterize the formation of the TREX complex and the AREX complex. In the ATP depleted condition, UAP56 formed an Apo-TREX complex containing the THO subcomplex but not ALYREF and CIP29. URH49 formed an Apo-AREX complex containing CIP29 but not ALYREF and the THO subcomplex. However, with the addition of ATP, both the Apo-TREX complex and the Apo-AREX complex were remodeled to highly similar ATP-TREX complex containing the THO subcomplex, ALYREF and CIP29. The knockdown of URH49 caused a reduction in its target mRNAs and a cytokinesis failure. Similarly, cytokinesis abnormality was observed in CIP29 knockdown cells, suggesting that CIP29 belongs to the URH49 regulated mRNA export pathway. Lastly, we confirmed that the export of mRNA in URH49-dependent pathway is achieved by NXF1, which is also observed in UAP56-dependent pathway. Our studies propose an mRNA export model that the mRNA selectivity depends on the Apo-form TREX/AREX complex, which is remodeled to the highly similar ATP-form complex upon ATP loading, and integrated to NXF1.

Effects of β-HgS on cell viability and intracellular oxidative stress in PC-12 cells

Traditional Tibetan medicines containing β-HgS have been used to treat chronic ailments for thousands of years. The effects were studied of β-HgS on cell viability and intracellular oxidative stress in PC-12 cells.

XAB2 depletion induces intron retention in POLR2A to impair global transcription and promote cellular senescence

XAB2 is a multi-functional protein participating processes including transcription, splicing, DNA repair and mRNA export. Here, we report POLR2A, the largest catalytic subunit of RNA polymerase II, as a major target gene down-regulated after XAB2 depletion. XAB2 depletion led to severe splicing defects of POLR2A with significant intron retention. Such defects resulted in substantial loss of POLR2A at RNA and protein levels, which further impaired global transcription. Treatment of splicing inhibitor madrasin induced similar reduction of POLR2A. Screen using TMT-based quantitative proteomics identified several proteins involved in mRNA surveillance including Dom34 with elevated expression. Inhibition of translation or depletion of Dom34 rescued the expression of POLR2A by stabilizing its mRNA. Immuno-precipitation further confirmed that XAB2 associated with spliceosome components important to POLR2A expression. Domain mapping revealed that TPR motifs 2–4 and 11 of XAB2 were critical for POLR2A expression by interacting with SNW1. Finally, we showed POLR2A mediated cell senescence caused by XAB2 deficiency. Depletion of XAB2 or POLR2A induced cell senescence by up-regulation of p53 and p21, re-expression of POLR2A after XAB2 depletion alleviated cellular senescence. These data together support that XAB2 serves as a guardian of POLR2A expression to ensure global gene expression and antagonize cell senescence.

Integrated Proteogenomic Characterization of Clear Cell Renal Cell Carcinoma

To elucidate the deregulated functional modules that drive clear cell renal cell carcinoma (ccRCC), we performed comprehensive genomic, epigenomic, transcriptomic, proteomic, and phosphoproteomic characterization of treatment-naive ccRCC and paired normal adjacent tissue samples. Genomic analyses identified a distinct molecular subgroup associated with genomic instability. Integration of proteogenomic measurements uniquely identified protein dysregulation of cellular mechanisms impacted by genomic alterations, including oxidative phosphorylation-related metabolism, protein translation processes, and phospho-signaling modules. To assess the degree of immune infiltration in individual tumors, we identified microenvironment cell signatures that delineated four immune-based ccRCC subtypes characterized by distinct cellular pathways. This study reports a large-scale proteogenomic analysis of ccRCC to discern the functional impact of genomic alterations and provides evidence for rational treatment selection stemming from ccRCC pathobiology.

The expression tractability of biological traits shaped by natural selection

Understanding how gene expression is translated to phenotype is central to modern molecular biology, and the success is contingent on the intrinsic tractability of the specific traits under examination. However, an a priori estimate of trait tractability from the perspective of gene expression is unavailable. Motivated by the concept of entropy in a thermodynamic system, we here propose such an estimate (S_T) by gauging the number (N) of expression states that underlie the same trait abnormality, with large S_T corresponding to large N. By analyzing over 200 yeast morphological traits, we show that S_T predicts the tractability of an expression-trait relationship. We further show that S_T is ultimately determined by natural selection, which builds co-regulated gene modules to minimize possible expression states.

Transmission of dynamic supercoiling in linear and multi-way branched DNAs and its regulation revealed by a fluorescent G-quadruplex torsion sensor

DNA supercoiling is an important regulator of gene activity. The transmission of transcription-generated supercoiling wave along a DNA helix provides a way for a gene being transcribed to communicate with and regulate its neighboring genes. Currently, the dynamic behavior of supercoiling transmission remains unclear owing to the lack of a suitable tool for detecting the dynamics of supercoiling transmission. In this work, we established a torsion sensor that quantitatively monitors supercoiling transmission in real time in DNA. Using this sensor, we studied the transmission of transcriptionally generated negative supercoiling in linear and multi-way DNA duplexes. We found that transcription-generated dynamic supercoiling not only transmits along linear DNA duplex but also equally diverges at and proceeds through multi-way DNA junctions. We also show that such a process is regulated by DNA–protein interactions and non-canonical DNA structures in the path of supercoiling transmission. These results imply a transcription-coupled mechanism of dynamic supercoiling-mediated intra- and inter-chromosomal signal transduction pathway and their regulation in DNA.

Coordination Analysis of Gene Expression Points to the Relative Impact of Different Regulators During Endoplasmic Reticulum Stress

Analysis of gene expression can be challenging, especially if it involves genetically diverse populations that exhibit high variation in their individual expression profile. Despite this variation, it is conceivable that in the same individuals a high degree of coordination is maintained between transcripts that belong to the same signaling modules and are associated with related biological functions. To explore this further, we calculated the correlation in the expression levels between each of ATF4, CHOP (DDIT3), GRP94, DNAJB9 (ERdj4), DNAJ3C (P58^IPK), and HSPA5 (BiP/GRP78) with the whole transcriptome in primary fibroblasts from deer mice following induction of endoplasmic reticulum (ER) stress. Since these genes are associated with different transducers of the unfolded protein response (UPR), we postulated that their profile, in terms of correlation of transcripts, reflects distinct UPR branches engaged, and therefore different biological processes. Standard gene ontology analysis was able to predict major functions associated with the corresponding transcript, and of the UPR arm related to that, namely regulation of the apoptotic response by ATF4 (PERK arm) and the ER stress-associated degradation for GRP94 (IRE1). BiP, being a global regulator of the UPR, was associated with activation of ER stress in a rather global manner. Pairwise comparison in the correlation coefficients for these genes' associated transcriptome showed the relevance of selected genes in terms of expression profiles. Conventional assessment of differential gene expression was incapable of providing meaningful information and pointed only to a generic association with stress. Collectively, this approach suggests that by evaluating the degree of coordination in gene expression, in genetically diverse biological specimens, may be useful in assigning genes in transcriptome networks, and more importantly in linking signaling nodules to specific biological functions and processes.

Estimation of Transcription Factor Activity in Knockdown Studies

Numerous methods have been developed trying to infer actual regulatory events in a sample. A prominent class of methods model genome-wide gene expression as linear equations derived from a transcription factor (TF) – gene network and optimizes parameters to fit the measured expression intensities. We apply four such methods on experiments with a TF-knockdown (KD) in human and E. coli. The transcriptome data provides clear expression signals and thus represents an extremely favorable test setting. The methods estimate activity changes of all TFs, which we expect to be highest in the KD TF. However, only in 15 out of 54 cases, the KD TFs ranked in the top 5%. We show that this poor overall performance cannot be attributed to a low effectiveness of the knockdown or the specific regulatory network provided as background knowledge. Further, the ranks of regulators related to the KD TF by the network or pathway are not significantly different from a random selection. In general, the result overlaps of different methods are small, indicating that they draw very different conclusions when presented with the same, presumably simple, inference problem. These results show that the investigated methods cannot yield robust TF activity estimates in knockdown schemes.

A Combined Transcriptomic and Genomic Analysis Identifies a Gene Signature Associated With the Response to Anti-TNF Therapy in Rheumatoid Arthritis

Background: Rheumatoid arthritis (RA) is the most frequent autoimmune disease involving the joints. Although anti-TNF therapies have proven effective in the management of RA, approximately one third of patients do not show a significant clinical response. The objective of this study was to identify new genetic variation associated with the clinical response to anti-TNF therapy in RA.

Methods: We performed a sequential multi-omic analysis integrating different sources of molecular information. First, we extracted the RNA from synovial biopsies of 11 RA patients starting anti-TNF therapy to identify gene coexpression modules (GCMs) in the RA synovium. Second, we analyzed the transcriptomic association between each GCM and the clinical response to anti-TNF therapy. The clinical response was determined at week 14 using the EULAR criteria. Third, we analyzed the association between the GCMs and anti-TNF response at the genetic level. For this objective, we used genome-wide data from a cohort of 348 anti-TNF treated patients from Spain. The GCMs that were significantly associated with the anti-TNF response were then tested for validation in an independent cohort of 2,706 anti-TNF treated patients. Finally, the functional implication of the validated GCMs was evaluated via pathway and cell type epigenetic enrichment analyses.

Results: A total of 149 GCMs were identified in the RA synovium. From these, 13 GCMs were found to be significantly associated with anti-TNF response (P < 0.05). At the genetic level, we detected two of the 13 GCMs to be significantly associated with the response to adalimumab (P = 0.0015) and infliximab (P = 0.021) in the Spain cohort. Using the independent cohort of RA patients, we replicated the association of the GCM associated with the response to adalimumab (P = 0.0019). The validated module was found to be significantly enriched for genes involved in the nucleotide metabolism (P = 2.41e-5) and epigenetic marks from immune cells, including CD4+ regulatory T cells (P = 0.041).

Conclusions: These findings show the existence of a drug-specific genetic basis for anti-TNF response, thereby supporting treatment stratification in the search for response biomarkers in RA.

Multiple roles of CTDK-I throughout the cell

The heterotrimeric carboxy-terminal domain kinase I (CTDK-I) in yeast is a cyclin-dependent kinase complex that is evolutionally conserved throughout eukaryotes and phosphorylates the C-terminal domain of the largest subunit of RNA polymerase II (RNApII) on the second-position serine (Ser2) residue of YSPTSPS heptapeptide repeats. CTDK-I plays indispensable roles in transcription elongation and transcription-coupled processing, such as the 3'-end processing of nascent mRNA transcripts. However, recent studies have revealed additional roles of CTDK-I beyond its primary effect on transcription by RNApII. Here, we describe recent advances in the regulation of genomic stability and rDNA integrity by CTDK-I and highlight the previously underappreciated cellular roles of CTDK-I in rRNA synthesis by RNA polymerase I and translational initiation and elongation. These multiple roles of CTDK-I throughout the cell expand our understanding of how this complex functions to coordinate diverse cellular processes through gene expression and how the human orthologue exerts its roles in diseased states such as tumorigenesis.

LncRNA EPR controls epithelial proliferation by coordinating Cdkn1a transcription and mRNA decay response to TGF-β

Long noncoding RNAs (lncRNAs) are emerging as regulators of fundamental biological processes. Here we report on the characterization of an intergenic lncRNA expressed in epithelial tissues which we termed EPR (Epithelial cell Program Regulator). EPR is rapidly downregulated by TGF-β and its sustained expression largely reshapes the transcriptome, favors the acquisition of epithelial traits, and reduces cell proliferation in cultured mammary gland cells as well as in an animal model of orthotopic transplantation. EPR generates a small peptide that localizes at epithelial cell junctions but the RNA molecule per se accounts for the vast majority of EPR-induced gene expression changes. Mechanistically, EPR interacts with chromatin and regulates Cdkn1a gene expression by affecting both its transcription and mRNA decay through its association with SMAD3 and the mRNA decay-promoting factor KHSRP, respectively. We propose that EPR enables epithelial cells to control proliferation by modulating waves of gene expression in response to TGF-β.

Several lncRNAs are regulated by TGF-β. Here the authors report that an intergenic lncRNA —EPR— is a component of the TGF-β signaling pathway and controls epithelial cell proliferation by altering transcription and mRNA decay of Cdkn1a. EPR overexpression restrains tumor growth of orthotopically transplanted mice.

miR-193a targets MLL1 mRNA and drastically decreases MLL1 protein production: Ectopic expression of the miRNA aberrantly lowers H3K4me3 content of the chromatin and hampers cell proliferation and viability

Mixed-lineage leukaemia 1 (MLL1) enzyme plays major role in regulating genes associated with vertebrate development. Cell physiology and homeostasis is regulated by microRNAs in diverse microenvironment. In this investigation we have identified conserved miR-193a target sites within the 3'-UTR of MLL1 gene transcript. Utilizing wild type and mutated 3'-UTR constructs and luciferase reporter assays we have clearly demonstrated that miR-193a directly targets the 3'-UTR region of the MLL1 mRNA. Ectopic expression of miR-193a modulated global H3K4 mono-, di- and tri-methylation levels and affects the expression of CAV1, a gene which is specifically modulated by H3K4me3. To determine the implications of this in vitro finding in aberrant physiological conditions we analyzed prostate cancer tissue samples. In this context miR-193a RNA was undetectable and MLL1 was highly expressed with concomitantly high levels of H3K4me, H3K4me2, and H3K4me3 enrichment in the promoters of MLL1 responsive genes. Finally, we showed that prolonged ectopic expression of miR-193a inhibits growth and cell migration, and induces apoptosis. Thus, while our study unveils amplitude of the epigenome, including miRnome it establishes that; (i) miR-193a directly target MLL1 mRNA, (ii) miR-193a impair MLL1 protein production, (iii) miR-193a reduces the overall methylation marks of the genome.

Post-transcriptional regulatory patterns revealed by protein-RNA interactions

The coordination of the synthesis of functionally-related proteins can be achieved at the post-transcriptional level by the action of common regulatory molecules, such as RNA–binding proteins (RBPs). Despite advances in the genome-wide identification of RBPs and their binding transcripts, the protein–RNA interaction space is still largely unexplored, thus hindering a broader understanding of the extent of the post-transcriptional regulation of related coding RNAs. Here, we propose a computational approach that combines protein–mRNA interaction networks and statistical analyses to provide an inferred regulatory landscape for more than 800 human RBPs and identify the cellular processes that can be regulated at the post-transcriptional level. We show that 10% of the tested sets of functionally-related mRNAs can be post-transcriptionally regulated. Moreover, we propose a classification of (i) the RBPs and (ii) the functionally-related mRNAs, based on their distinct behaviors in the functional landscape, hinting towards mechanistic regulatory hypotheses. In addition, we demonstrate the usefulness of the inferred functional landscape to investigate the cellular role of both well-characterized and novel RBPs in the context of human diseases.

O-GlcNAcylation alters the selection of mRNAs for translation and promotes 4E-BP1-dependent mitochondrial dysfunction in the retina

Diabetes promotes the posttranslational modification of proteins by O-linked addition of GlcNAc (O-GlcNAcylation) to Ser/Thr residues of proteins and thereby contributes to diabetic complications. In the retina of diabetic mice, the repressor of mRNA translation, eIF4E-binding protein 1 (4E-BP1), is O-GlcNAcylated, and sequestration of the cap-binding protein eukaryotic translation initiation factor (eIF4E) is enhanced. O-GlcNAcylation has also been detected on several eukaryotic translation initiation factors and ribosomal proteins. However, the functional consequence of this modification is unknown. Here, using ribosome profiling, we evaluated the effect of enhanced O-GlcNAcylation on retinal gene expression. Mice receiving thiamet G (TMG), an inhibitor of the O-GlcNAc hydrolase O-GlcNAcase, exhibited enhanced retinal protein O-GlcNAcylation. The principal effect of TMG on retinal gene expression was observed in ribosome-associated mRNAs (i.e. mRNAs undergoing translation), as less than 1% of mRNAs exhibited changes in abundance. Remarkably, ∼19% of the transcriptome exhibited TMG-induced changes in ribosome occupancy, with 1912 mRNAs having reduced and 1683 mRNAs having increased translational rates. In the retina, the effect of O-GlcNAcase inhibition on translation of specific mitochondrial proteins, including superoxide dismutase 2 (SOD2), depended on 4E-BP1/2. O-GlcNAcylation enhanced cellular respiration and promoted mitochondrial superoxide levels in WT cells, and 4E-BP1/2 deletion prevented O-GlcNAcylation-induced mitochondrial superoxide in cells in culture and in the retina. The retina of diabetic WT mice exhibited increased reactive oxygen species levels, an effect not observed in diabetic 4E-BP1/2-deficient mice. These findings provide evidence for a mechanism whereby diabetes-induced O-GlcNAcylation promotes oxidative stress in the retina by altering the selection of mRNAs for translation.

Transcriptionally induced enhancers in the macrophage immune response to Mycobacterium tuberculosis infection

Background

Tuberculosis is a life-threatening infectious disease caused by Mycobacterium tuberculosis (M.tb). M.tb subverts host immune responses to build a favourable niche and survive inside of host macrophages. Macrophages can control or eliminate the infection, if acquire appropriate functional phenotypes. Transcriptional regulation is a key process that governs the activation and maintenance of these phenotypes. Among the factors orchestrating transcriptional regulation during M.tb infection, transcriptional enhancers still remain unexplored.

Results

We analysed transcribed enhancers in M.tb-infected mouse bone marrow-derived macrophages. We established a link between known M.tb-responsive transcription factors and transcriptional activation of enhancers and their target genes. Our data suggest that enhancers might drive macrophage response via transcriptional activation of key immune genes, such as Tnf, Tnfrsf1b, Irg1, Hilpda, Ccl3, and Ccl4. We report enhancers acquiring transcription de novo upon infection. Finally, we link highly transcriptionally induced enhancers to activation of genes with previously unappreciated roles in M.tb infection, such as Fbxl3, Tapt1, Edn1, and Hivep1.

Conclusions

Our findings suggest the importance of macrophage host transcriptional enhancers during M.tb infection. Our study extends current knowledge of the regulation of macrophage responses to M.tb infection and provides a basis for future functional studies on enhancer-gene interactions in this process.

Electronic supplementary material

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

The Human Cell Atlas: Technical approaches and challenges

The Human Cell Atlas is a large, international consortium that aims to identify and describe every cell type in the human body. The comprehensive cellular maps that arise from this ambitious effort have the potential to transform many aspects of fundamental biology and clinical practice. Here, we discuss the technical approaches that could be used today to generate such a resource and also the technical challenges that will be encountered.

Modeling Transcriptional Rewiring in Neutrophils Through the Course of Treated Juvenile Idiopathic Arthritis

Neutrophils in children with the polyarticular form of juvenile idiopathic arthritis (JIA) display abnormal transcriptional patterns linked to fundamental metabolic derangements. In this study, we sought to determine the effects of therapy on mRNA and miRNA expression networks in polyarticular JIA. Using exon and miRNA microarrays, we studied children with untreated active JIA (ADU, n = 35), children with active disease on therapy with methotrexate ± etanercept (ADT, n = 26), and children with inactive disease also on therapy (ID, n = 14). We compared the results to findings from healthy control children (HC, n = 35). We found substantial re-ordering of mRNA and miRNA expression networks after the initiation of therapy. Each disease state was associated with a distinct transcriptional profile, with the ADT state differing the most from HC, and ID more strongly resembling HC. Changes at the mRNA level were mirrored in changes in miRNA expression patterns. The analysis of the expression dynamics from differentially expressed genes across three disease states indicated that therapeutic response is a complex process. This process does not simply involve genes slowly correcting in a linear fashion over time. Computational modeling of miRNA and transcription factor (TF) co-regulatory networks demonstrated that combinational regulation of miRNA and TF might play an important role in dynamic transcriptome changes.

Beyond Transcription: Roles of Transcription Factors in Pre-mRNA Splicing

Whereas individual steps of protein-coding gene expression in eukaryotes can be studied in isolation in vitro, it has become clear that these steps are intimately connected within cells. Connections not only ensure quality control but also fine-tune the gene expression process, which must adapt to environmental changes while remaining robust. In this review, we systematically present proven and potential mechanisms by which sequence-specific DNA-binding transcription factors can alter gene expression beyond transcription initiation and regulate pre-mRNA splicing, and thereby mRNA isoform production, by (i) influencing transcription elongation rates, (ii) binding to pre-mRNA to recruit splicing factors, and/or (iii) blocking the association of splicing factors with pre-mRNA. We propose various mechanistic models throughout the review, in some cases without explicit supportive evidence, in hopes of providing fertile ground for future studies.

Crosstalk between Receptor and Non-receptor Mediated Chemical Modes of Action in Rat Livers Converges through a Dysregulated Gene Expression Network at Tumor Suppressor Tp53

Chemicals, toxicants, and environmental stressors mediate their biologic effect through specific modes of action (MOAs). These encompass key molecular events that lead to changes in the expression of genes within regulatory pathways. Elucidating shared biologic processes and overlapping gene networks will help to better understand the toxicologic effects on biological systems. In this study we used a weighted network analysis of gene expression data from the livers of male Sprague-Dawley rats exposed to chemicals that elicit their effects through receptor-mediated MOAs (aryl hydrocarbon receptor, orphan nuclear hormone receptor, or peroxisome proliferator-activated receptor-α) or non-receptor-mediated MOAs (cytotoxicity or DNA damage). Four gene networks were highly preserved and statistically significant in each of the two MOA classes. Three of the four networks contain genes that enrich for immunity and defense. However, many canonical pathways related to an immune response were activated from exposure to the non-receptor-mediated MOA chemicals and deactivated from exposure to the receptor-mediated MOA chemicals. The top gene network contains a module with 33 genes including tumor suppressor TP53 as the central hub which was slightly up-regulated in gene expression compared to control. Although, there is crosstalk between the two MOA classes of chemicals at the TP53 gene network, more than half of the genes are dysregulated in opposite directions. For example, Thromboxane A Synthase 1 (Tbxas1), a cytochrome P450 protein coding gene regulated by Tp53, is down-regulated by exposure to the receptor-mediated chemicals but up-regulated by the non-receptor-mediated chemicals. The regulation of gene expression by the chemicals within MOA classes was consistent despite varying alanine transaminase and aspartate aminotransferase liver enzyme measurements. These results suggest that overlap in toxicologic pathways by chemicals with different MOAs provides common mechanisms for discordant regulation of gene expression within molecular networks.

Genome-wide profiling of transcribed enhancers during macrophage activation

BACKGROUND

Macrophages are sentinel cells essential for tissue homeostasis and host defence. Owing to their plasticity, macrophages acquire a range of functional phenotypes in response to microenvironmental stimuli, of which M(IFN-γ) and M(IL-4/IL-13) are well known for their opposing pro- and anti-inflammatory roles. Enhancers have emerged as regulatory DNA elements crucial for transcriptional activation of gene expression.

RESULTS

Using cap analysis of gene expression and epigenetic data, we identify on large-scale transcribed enhancers in bone marrow-derived mouse macrophages, their time kinetics, and target protein-coding genes. We observe an increase in target gene expression, concomitant with increasing numbers of associated enhancers, and find that genes associated with many enhancers show a shift towards stronger enrichment for macrophage-specific biological processes. We infer enhancers that drive transcriptional responses of genes upon M(IFN-γ) and M(IL-4/IL-13) macrophage activation and demonstrate stimuli specificity of regulatory associations. Finally, we show that enhancer regions are enriched for binding sites of inflammation-related transcription factors, suggesting a link between stimuli response and enhancer transcriptional control.

CONCLUSIONS

Our study provides new insights into genome-wide enhancer-mediated transcriptional control of macrophage genes, including those implicated in macrophage activation, and offers a detailed genome-wide catalogue of transcribed enhancers in bone marrow-derived mouse macrophages.

Transcriptome analysis of IFM-specific actin and myosin nulls in Drosophila melanogaster unravels lesion-specific expression blueprints across muscle mutations

Muscle contraction is a highly fine-tuned process that requires the precise and timely construction of large protein sub-assemblies to form sarcomeres. Mutations in many genes encoding constituent proteins of this macromolecular machine result in defective functioning of the muscle tissue. However, the pathways underlying muscle degeneration, and manifestation of myopathy phenotypes are not well understood. In this study, we explored transcriptional alterations that ensue from the absence of the two major muscle proteins - myosin and actin - using the Drosophila indirect flight muscles. Our aim was to understand how the muscle tissue responds as a whole to the absence of either of the major scaffold proteins, whether the responses are generic to the tissue; or unique to the thick versus thin filament systems. Our results indicated that muscles respond by altering gene transcriptional levels in multiple systems active in muscle remodelling, protein degradation and heat shock responses. However, there were some responses that were filament-specific signatures of muscle degeneration, like immune responses, metabolic alterations and alterations in expression of muscle structural genes and mitochondrial ribosomal genes. These general and filament-specific changes in gene expression may be of relevance to human myopathies.