The evolution of posttranscriptional regulation

Multi-omics analysis of aggregative multicellularity

All organisms have to carefully regulate their gene expression, not least during development. mRNA levels are often used as proxy for protein output; however, this approach ignores post-transcriptional effects. In particular, mRNA-protein correlation remains elusive for organisms that exhibit aggregative rather than clonal multicellularity. We addressed this issue by generating a paired transcriptomics and proteomics time series during the transition from uni-to multicellular stage in the social ameba Dictyostelium discoideum. Our data reveals that mRNA and protein levels correlate highly during unicellular growth, but decrease when multicellular development is initiated. This accentuates that transcripts alone cannot accurately predict protein levels. The dataset provides a useful resource to study gene expression during aggregative multicellular development. Additionally, our study provides an example of how to analyze and visualize mRNA and protein levels, which should be broadly applicable to other organisms and conditions.

Intake of imidacloprid in lethal and sublethal doses alters gene expression in Apis mellifera bees

Bees are important pollinators for ecosystems and agriculture; however, populations have suffered a decline that may be associated with several factors, including habitat loss, climate change, increased vulnerability to diseases and parasites and use of pesticides. The extensive use of neonicotinoids, including imidacloprid, as agricultural pesticides, leads to their persistence in the environment and accumulation in bees, pollen, nectar, and honey, thereby inducing deleterious effects. Forager honey bees face significant exposure to pesticide residues while searching for resources outside the hive, particularly systemic pesticides like imidacloprid. In this study, 360 Apis mellifera bees, twenty-one days old (supposed to be in the forager phase) previously marked were fed syrup (honey and water, 1:1 m/v) containing a lethal dose (0.081 μg/bee) or sublethal dose (0.00081 μg/bee) of imidacloprid. The syrup was provided in plastic troughs, with 250 μL added per trough onto each plastic Petri dish containing 5 bees (50 μL per bee). The bees were kept in the plastic Petri dishes inside an incubator, and after 1 and 4 h of ingestion, the bees were euthanised and stored in an ultra-freezer (-80 °C) for transcriptome analysis. Following the 1-h ingestion of imidacloprid, 1516 genes (73 from lethal dose; 1509 from sublethal dose) showed differential expression compared to the control, while after 4 h, 758 genes (733 from lethal dose; 25 from sublethal) exhibited differential expression compared to the control. All differentially expressed genes found in the brain tissue transcripts of forager bees were categorised based on gene ontology into functional groups encompassing biological processes, molecular functions, and cellular components. These analyses revealed that sublethal doses might be capable of altering more genes than lethal doses, potentially associated with a phenomenon known as insecticide-induced hormesis. Alterations in genes related to areas such as the immune system, nutritional metabolism, detoxification system, circadian rhythm, odour detection, foraging activity, and memory in bees were present after exposure to the pesticide. These findings underscore the detrimental effects of both lethal and sublethal doses of imidacloprid, thereby providing valuable insights for establishing public policies regarding the use of neonicotinoids, which are directly implicated in the compromised health of Apis mellifera bees.

Hydrolyzed Chicken Meat Extract and Its Bioactive Cyclopeptides Protect Neural Function by Attenuating Inflammation and Apoptosis via PI3K/AKT and AMPK Pathways

Cognitive decline is inevitable with age, and due to the lack of well-established pharmacotherapies for neurodegenerative disorders, dietary supplements have become important alternatives to ameliorate brain deterioration. Hydrolyzed chicken meat extract (HCE) and its bioactive components were previously found to improve neuroinflammation and cognitive decline by regulating microglia polarization. However, the effects and mechanisms of these bioactives on neurons remain unclear. Here, the most potent bioactive component on neural function in HCE was screened out, and the detailed mechanism was clarified through in vivo and in vitro experiments. We found that HCE, cyclo(Val-Pro), cyclo(Phe-Phe), cyclo(His-Pro), cyclo(Leu-Lys), and arginine exerted stronger anti-inflammatory and antioxidant effects among the 12 bioactives in amyloid β (Aβ)-treated HT-22 cells. Further transcriptome sequencing and polymerase chain reaction (PCR) array analysis showed that these bioactives participated in different signaling pathways, and cyclo(Val-Pro) was identified as the most potent cyclic dipeptide. In addition, the antiapoptotic and neuroprotective effect of cyclo(Val-Pro) was partly regulated by the activation of PI3K/AKT and AMPK pathways, and the inhibition of these pathways abolished the effect of cyclo(Val-Pro). Moreover, cyclo(Val-Pro) enhanced cognitive function and neurogenesis and alleviated neuroinflammation and oxidative stress in middle-aged mice, with an effect similar to HCE. Hippocampal transcriptome analysis further revealed that HCE and cyclo(Val-Pro) significantly enriched the neuroactive ligand-receptor interaction pathway, verified by enhanced neurotransmitter levels and upregulated neurotransmitter receptor-related gene expression. Therefore, the mechanism of cyclo(Val-Pro) on neural function might be associated with PI3K/AKT and AMPK pathway-mediated antiapoptotic effect and neurogenesis and the activation of the neurotransmitter-receptor pathway.

Current technical advancements in plant epitranscriptomic studies

The growth and development of plants are the result of the interplay between the internal developmental programming and plant-environment interactions. Gene expression regulations in plants are made up of multi-level networks. In the past few years, many studies were carried out on co- and post-transcriptional RNA modifications, which, together with the RNA community, are collectively known as the "epitranscriptome." The epitranscriptomic machineries were identified and their functional impacts characterized in a broad range of physiological processes in diverse plant species. There is mounting evidence to suggest that the epitranscriptome provides an additional layer in the gene regulatory network for plant development and stress responses. In the present review, we summarized the epitranscriptomic modifications found so far in plants, including chemical modifications, RNA editing, and transcript isoforms. The various approaches to RNA modification detection were described, with special emphasis on the recent development and application potential of third-generation sequencing. The roles of epitranscriptomic changes in gene regulation during plant-environment interactions were discussed in case studies. This review aims to highlight the importance of epitranscriptomics in the study of gene regulatory networks in plants and to encourage multi-omics investigations using the recent technical advancements.

The Genome-Wide Characterization of Alternative Splicing and RNA Editing in the Development of Coprinopsis cinerea

Coprinopsis cinerea is one of the model species used in fungal developmental studies. This mushroom-forming Basidiomycetes fungus has several developmental destinies in response to changing environments, with dynamic developmental regulations of the organism. Although the gene expression in C. cinerea development has already been profiled broadly, previous studies have only focused on a specific stage or process of fungal development. A comprehensive perspective across different developmental paths is lacking, and a global view on the dynamic transcriptional regulations in the life cycle and the developmental paths is far from complete. In addition, knowledge on co- and post-transcriptional modifications in this fungus remains rare. In this study, we investigated the transcriptional changes and modifications in C. cinerea during the processes of spore germination, vegetative growth, oidiation, sclerotia formation, and fruiting body formation by inducing different developmental paths of the organism and profiling the transcriptomes using the high-throughput sequencing method. Transition in the identity and abundance of expressed genes drive the physiological and morphological alterations of the organism, including metabolism and multicellularity construction. Moreover, stage- and tissue-specific alternative splicing and RNA editing took place and functioned in C. cinerea. These modifications were negatively correlated to the conservation features of genes and could provide extra plasticity to the transcriptome during fungal development. We suggest that C. cinerea applies different molecular strategies in its developmental regulation, including shifts in expressed gene sets, diversifications of genetic information, and reversible diversifications of RNA molecules. Such features would increase the fungal adaptability in the rapidly changing environment, especially in the transition of developmental programs and the maintenance and balance of genetic and transcriptomic divergence. The multi-layer regulatory network of gene expression serves as the molecular basis of the functioning of developmental regulation.

Post-Transcriptional and Epigenetic Regulation of Estrogen Signaling

Post-translational and epigenetic regulation are important mechanisms controlling functions of genes and proteins. Although the "classic" estrogen receptors (ERs) have been acknowledged to function in mediating estrogen effects via transcriptional mechanisms, estrogenic agents modulate the turnover of several proteins via post-transcriptional and post-translational pathways including epigenetics. For instance, the metabolic and angiogenic action of G-protein coupled estrogen receptor (GPER) in vascular endothelial cells has been recently elucidated. By interacting with GPER, 17β-estradiol and the GPER agonist G1 enhance endothelial stability of 6-phosphofructo-2-kinase/fructose-2,6-biphosphatase 3 (PFKFB3) and capillary tube formation by increasing ubiquitin-specific peptidase 19 levels, thereby reducing PFKFB3 ubiquitination and proteasomal degradation. In addition to ligands, the functional expression and trafficking of ERs can be modulated by post-translational modification, including palmitoylation. MicroRNAs (miRNAs), the most abundant form of endogenous small RNAs in humans, regulate multiple target genes and are at the center of the multi-target regulatory network. This review also discusses the emerging evidence of how miRNAs affect glycolytic metabolism in cancer, as well as their regulation by estrogens. Restoring dysregulated miRNA expression represents a promising strategy to counteract the progression of cancer and other disease conditions. Accordingly, estrogen post-transcriptional regulatory and epigenetic mechanisms represent novel targets for pharmacological and nonpharmacological intervention for the treatment and prevention of hormone-sensitive noncommunicable diseases, including estrogen-sensitive cancers of the reproductive system in women. SIGNIFICANCE STATEMENT: The effects of estrogen are mediated by several mechanisms that are not limited to the transcriptional regulation of target genes. Slowing down the turnover of master regulators of metabolism by estrogens allows cells to rapidly adapt to environmental cues. Identification of estrogen-targeted microRNAs may lead to the development of novel RNA therapeutics that disrupt pathological angiogenesis in estrogen-dependent cancers.

Analysis of subcellular RNA fractions demonstrates significant genetic regulation of gene expression in human brain post-transcriptionally

Gaining insight into the genetic regulation of gene expression in human brain is key to the interpretation of genome-wide association studies for major neurological and neuropsychiatric diseases. Expression quantitative trait loci (eQTL) analyses have largely been used to achieve this, providing valuable insights into the genetic regulation of steady-state RNA in human brain, but not distinguishing between molecular processes regulating transcription and stability. RNA quantification within cellular fractions can disentangle these processes in cell types and tissues which are challenging to model in vitro. We investigated the underlying molecular processes driving the genetic regulation of gene expression specific to a cellular fraction using allele-specific expression (ASE). Applying ASE analysis to genomic and transcriptomic data from paired nuclear and cytoplasmic fractions of anterior prefrontal cortex, cerebellar cortex and putamen tissues from 4 post-mortem neuropathologically-confirmed control human brains, we demonstrate that a significant proportion of genetic regulation of gene expression occurs post-transcriptionally in the cytoplasm, with genes undergoing this form of regulation more likely to be synaptic. These findings have implications for understanding the structure of gene expression regulation in human brain, and importantly the interpretation of rapidly growing single-nucleus brain RNA-sequencing and eQTL datasets, where cytoplasm-specific regulatory events could be missed.

The Glomerulus Multiomics Analysis Provides Deeper Insights into Diabetic Nephropathy

Although diabetic nephropathy (DN) is the leading cause of the end-stage renal disease, the exact regulation mechanisms remain unknown. In this study, we integrated the transcriptomics and proteomics profiles of glomeruli isolated from 50 biopsy-proven DN patients and 25 controls to investigate the latest findings about DN pathogenesis. First, 1152 genes exhibited differential expression at the mRNA or protein level, and 364 showed significant association. These strong correlated genes were divided into four different functional modules. Moreover, a regulatory network of the transcription factors (TFs)-target genes (TGs) was constructed, with 30 TFs upregulated at the protein levels and 265 downstream TGs differentially expressed at the mRNA levels. These TFs are the integration centers of several signal transduction pathways and have tremendous therapeutic potential for regulating the aberrant production of TGs and the pathological process of DN. Furthermore, 29 new DN-specific splice-junction peptides were discovered with high confidence; these peptides may play novel functions in the pathological course of DN. So, our in-depth integrative transcriptomics-proteomics analysis provided deeper insights into the pathogenesis of DN and opened the potential avenue for finding new therapeutic interventions. MS raw files were deposited into the proteomeXchange with the dataset identifier PXD040617.

Multidimensional plasticity jointly contributes to rapid acclimation to environmental challenges during biological invasions

Rapid plastic response to environmental changes, which involves extremely complex underlying mechanisms, is crucial for organismal survival during many ecological and evolutionary processes such as those in global change and biological invasions. Gene expression is among the most studied molecular plasticity, while co- or posttranscriptional mechanisms are still largely unexplored. Using a model invasive ascidian Ciona savignyi, we studied multidimensional short-term plasticity in response to hyper- and hyposalinity stresses, covering the physiological adjustment, gene expression, alternative splicing (AS), and alternative polyadenylation (APA) regulations. Our results demonstrated that rapid plastic response varied with environmental context, timescales, and molecular regulatory levels. Gene expression, AS, and APA regulations independently acted on different gene sets and corresponding biological functions, highlighting their nonredundant roles in rapid environmental adaptation. Stress-induced gene expression changes illustrated the use of a strategy of accumulating free amino acids under high salinity and losing/reducing them during low salinity to maintain the osmotic homoeostasis. Genes with more exons were inclined to use AS regulations, and isoform switches in functional genes such as SLC2a5 and Cyb5r3 resulted in enhanced transporting activities by up-regulating the isoforms with more transmembrane regions. The extensive 3'-untranslated region (3'UTR) shortening through APA was induced by both salinity stresses, and APA regulation predominated transcriptomic changes at some stages of stress response. The findings here provide evidence for complex plastic mechanisms to environmental changes, and thereby highlight the importance of systemically integrating different levels of regulatory mechanisms in studying initial plasticity in evolutionary trajectories.

Not1 and Not4 inversely determine mRNA solubility that sets the dynamics of co-translational events

BACKGROUND

The Ccr4-Not complex is mostly known as the major eukaryotic deadenylase. However, several studies have uncovered roles of the complex, in particular of the Not subunits, unrelated to deadenylation and relevant for translation. In particular, the existence of Not condensates that regulate translation elongation dynamics has been reported. Typical studies that evaluate translation efficiency rely on soluble extracts obtained after the disruption of cells and ribosome profiling. Yet cellular mRNAs in condensates can be actively translated and may not be present in such extracts.

RESULTS

In this work, by analyzing soluble and insoluble mRNA decay intermediates in yeast, we determine that insoluble mRNAs are enriched for ribosomes dwelling at non-optimal codons compared to soluble mRNAs. mRNA decay is higher for soluble RNAs, but the proportion of co-translational degradation relative to the overall mRNA decay is higher for insoluble mRNAs. We show that depletion of Not1 and Not4 inversely impacts mRNA solubilities and, for soluble mRNAs, ribosome dwelling according to codon optimality. Depletion of Not4 solubilizes mRNAs with lower non-optimal codon content and higher expression that are rendered insoluble by Not1 depletion. By contrast, depletion of Not1 solubilizes mitochondrial mRNAs, which are rendered insoluble upon Not4 depletion.

CONCLUSIONS

Our results reveal that mRNA solubility defines the dynamics of co-translation events and is oppositely regulated by Not1 and Not4, a mechanism that we additionally determine may already be set by Not1 promoter association in the nucleus.

The formation of HuR/YB1 complex is required for the stabilization of target mRNA to promote myogenesis

mRNA stability is the mechanism by which cells protect transcripts allowing their expression to execute various functions that affect cell metabolism and fate. It is well-established that RNA binding proteins (RBPs) such as HuR use their ability to stabilize mRNA targets to modulate vital processes such as muscle fiber formation (myogenesis). However, the machinery and the mechanisms regulating mRNA stabilization are still elusive. Here, we identified Y-Box binding protein 1 (YB1) as an indispensable HuR binding partner for mRNA stabilization and promotion of myogenesis. Both HuR and YB1 bind to 409 common mRNA targets, 147 of which contain a U-rich consensus motif in their 3' untranslated region (3'UTR) that can also be found in mRNA targets in other cell systems. YB1 and HuR form a heterodimer that associates with the U-rich consensus motif to stabilize key promyogenic mRNAs. The formation of this complex involves a small domain in HuR (227-234) that if mutated prevents HuR from reestablishing myogenesis in siHuR-treated muscle cells. Together our data uncover that YB1 is a key player in HuR-mediated stabilization of pro-myogenic mRNAs and provide the first indication that the mRNA stability mechanism is as complex as other key cellular processes such as mRNA decay and translation.

The missing link between genetic association and regulatory function

The genetic basis of most traits is highly polygenic and dominated by non-coding alleles. It is widely assumed that such alleles exert small regulatory effects on the expression of cis-linked genes. However, despite the availability of gene expression and epigenomic datasets, few variant-to-gene links have emerged. It is unclear whether these sparse results are due to limitations in available data and methods, or to deficiencies in the underlying assumed model. To better distinguish between these possibilities, we identified 220 gene-trait pairs in which protein-coding variants influence a complex trait or its Mendelian cognate. Despite the presence of expression quantitative trait loci near most GWAS associations, by applying a gene-based approach we found limited evidence that the baseline expression of trait-related genes explains GWAS associations, whether using colocalization methods (8% of genes implicated), transcription-wide association (2% of genes implicated), or a combination of regulatory annotations and distance (4% of genes implicated). These results contradict the hypothesis that most complex trait-associated variants coincide with homeostatic expression QTLs, suggesting that better models are needed. The field must confront this deficit and pursue this 'missing regulation.'

Full-length transcriptome sequencing reveals the molecular mechanism of potato seedlings responding to low-temperature

BACKGROUND

Potato (Solanum tuberosum L.) is one of the world's most important crops, the cultivated potato is frost-sensitive, and low-temperature severely influences potato production. However, the mechanism by which potato responds to low-temperature stress is unclear. In this research, we apply a combination of second-generation sequencing and third-generation sequencing technologies to sequence full-length transcriptomes in low-temperature-sensitive cultivars to identify the important genes and main pathways related to low-temperature resistance.

RESULTS

In this study, we obtained 41,016 high-quality transcripts, which included 15,189 putative new transcripts. Amongst them, we identified 11,665 open reading frames, 6085 simple sequence repeats out of the potato dataset. We used public available genomic contigs to analyze the gene features, simple sequence repeat, and alternative splicing event of 24,658 non-redundant transcript sequences, predicted the coding sequence and identified the alternative polyadenylation. We performed cluster analysis, GO, and KEGG functional analysis of 4518 genes that were differentially expressed between the different low-temperature treatments. We examined 36 transcription factor families and identified 542 transcription factors in the differentially expressed genes, and 64 transcription factors were found in the AP2 transcription factor family which was the most. We measured the malondialdehyde, soluble sugar, and proline contents and the expression genes changed associated with low temperature resistance in the low-temperature treated leaves. We also tentatively speculate that StLPIN10369.5 and StCDPK16 may play a central coordinating role in the response of potatoes to low temperature stress.

CONCLUSIONS

Overall, this study provided the first large-scale full-length transcriptome sequencing of potato and will facilitate structure-function genetic and comparative genomics studies of this important crop.

Mammalian splicing divergence is shaped by drift, buffering in trans, and a scaling law

This study globally investigates the allelic splicing pattern in multiple tissues of an F1 hybrid mouse and reveals the underlying driving forces shaping such tissue-dependent splicing divergence.

Alternative splicing is ubiquitous, but the mechanisms underlying its pattern of evolutionary divergence across mammalian tissues are still underexplored. Here, we investigated the cis-regulatory divergences and their relationship with tissue-dependent trans-regulation in multiple tissues of an F1 hybrid between two mouse species. Large splicing changes between tissues are highly conserved and likely reflect functional tissue-dependent regulation. In particular, micro-exons frequently exhibit this pattern with high inclusion levels in the brain. Cis-divergence of splicing appears to be largely non-adaptive. Although divergence is in general associated with higher densities of sequence variants in regulatory regions, events with high usage of the dominant isoform apparently tolerate more mutations, explaining why their exon sequences are highly conserved but their intronic splicing site flanking regions are not. Moreover, we demonstrate that non-adaptive mutations are often masked in tissues where accurate splicing likely is more important, and experimentally attribute such buffering effect to trans-regulatory splicing efficiency.

Impact of medical imaging on the epigenome - low-dose exposure in the course of computed tomography does not induce detectable changes of DNA-methylation profiles in peripheral blood cells

BACKGROUND

Computed tomography (CT) is a main contributor to artificial low-dose exposure. Understanding the biological effects induced by CT exposure and their dependency on the characteristics of photon spectra is essential for knowledge-driven risk assessment. In a previous gene expression study, we have identified upregulation of AEN, BAX, DDB2, EDA2R and FDXR after ex vivo exposure with single-energy CT and dual-energy CT (DECT). In this study, we focused on CT-induced changes of DNA methylation. This epigenetic modification of DNA is a central regulator of gene expression and instrumental in preserving genome integrity. Previous studies reported focal hypermethylation and global hypomethylation after exposure with doses above 100 mSv, however, the effect of low dose exposure on DNA methylation is hardly explored.

MATERIALS AND METHODS

DNA was isolated from peripheral blood of three healthy individuals 6 h after ex vivo exposition to single-energy (80 kV and 150 kV) and DECT (80 kV/Sn150 kV) with a calculated effective dose of 7.0 ± 0.08 mSv. The experimental setting was identical to the one used in our previous gene expression study enabling a direct comparison of gene expression results with changes of DNA methylation identified in this study. DNA methylation was analyzed by high-throughput sequencing of bisulfite-treated DNA targeted methylation sequencing.

RESULTS

Unsupervised hierarchical clustering based on DNA methylation profiles of all samples created three distinct clusters. Formation of these three clusters was solely determined by the origin of samples, indicating the absence of prominent irradiation-associated changes of DNA methylation. In line with this observation, inter-individual comparison of non-irradiated samples revealed 1163, 1224 and 4550 significant differentially methylated regions (DMRs), respectively, whereas the pairwise comparison of irradiated and non-irradiated samples failed to identify irradiation-induced DMRs in any of the three probands. This even applied to the genomic regions harboring AEN, BAX, DDB2, EDA2R and FDXR, the five genes known to be upregulated by CT exposure.

CONCLUSIONS

CT exposure with various photon spectra did not result in detectable changes of DNA methylation. However, minor effects in a subpopulation of irradiated cells cannot be ruled out. Thus, future studies with extended observation intervals are needed to investigate DNA methylation changes that are induced by indirect effects at later points of time or become detectable by clonal expansion of affected cells. Moreover, our data suggest that DNA methylation analysis is less sensitive in detecting immediate effects of low-dose irradiation when compared to gene expression analysis.

Resolving the microalgal gene landscape at the strain level: A novel hybrid transcriptome of Emiliania huxleyi CCMP3266

Microalgae are key ecological players with a complex evolutionary history. Genomic diversity, in addition to limited availability of high-quality genomes, challenge studies that aim to elucidate molecular mechanisms underlying microalgal ecophysiology. Here, we present a novel and comprehensive transcriptomic hybrid approach to generate a reference for genetic analyses, and resolve the microalgal gene landscape at the strain level. The approach is demonstrated for a strain of the coccolithophore microalga Emiliania huxleyi, which is a species complex with considerable genome variability. The investigated strain is commonly studied as a model for algal-bacterial interactions, and was therefore sequenced in the presence of bacteria to elicit the expression of interaction-relevant genes. We applied complementary PacBio Iso-Seq full-length cDNA, and poly(A)-independent Illumina total RNA sequencing, which resulted in a de novo assembled, near complete hybrid transcriptome. In particular, hybrid sequencing improved the reconstruction of long transcripts and increased the recovery of full-length transcript isoforms. To use the resulting hybrid transcriptome as a reference for genetic analyses, we demonstrate a method that collapses the transcriptome into a genome-like dataset, termed "synthetic genome" (sGenome). We used the sGenome as a reference to visually confirm the robustness of the CCMP3266 gene assembly, to conduct differential gene expression analysis, and to characterize novel E. huxleyi genes. The newly-identified genes contribute to our understanding of E. huxleyi genome diversification, and are predicted to play a role in microbial interactions. Our transcriptomic toolkit can be implemented in various microalgae to facilitate mechanistic studies on microalgal diversity and ecology. Importance Microalgae are key players in the ecology and biogeochemistry of our oceans. Efforts to implement genomic and transcriptomic tools in laboratory studies involving microalgae suffer from the lack of published genomes. In the case of coccolithophore microalgae, the problem has long been recognized; the model species Emiliania huxleyi is a species complex with genomes composed of a core, and a large variable portion. To study the role of the variable portion in niche adaptation, and specifically in microbial interactions, strain-specific genetic information is required. Here we present a novel transcriptomic hybrid approach, and generated strain-specific genome-like information. We demonstrate our approach on an E. huxleyi strain that is co-cultivated with bacteria. By constructing a "synthetic genome", we generated comprehensive gene annotations that enabled accurate analyses of gene expression patterns. Importantly, we unveiled novel genes in the variable portion of E. huxleyi that play putative roles in microbial interactions.

An RNA-binding protein complex regulates the purine-dependent expression of a nucleobase transporter in trypanosomes

Post-transcriptional regulation of gene expression is particularly important in trypanosomatid protozoa. RNA-binding proteins (RBPs) regulate mRNA stability and translation, yet information about how RBPs are able to link environmental cues to post-transcriptional control is scarce. In Trypanosoma brucei, we have previously characterized a short RNA stem-loop cis-element (PuRE, Purine Responsive Element) within the 3'-UTR of the NT8 nucleobase transporter mRNA that is necessary and sufficient to confer a strong repression of gene expression in response to purines. In this study, we have identified a protein complex composed of two RNA-binding proteins (PuREBP1 and PuREBP2) that binds to the PuRE in vitro and to NT8 mRNA in vivo. Depletion of PuREBP1 by RNA interference results in the upregulation of just NT8 and the mRNAs encoding the amino acid transporter AATP6 paralogues. Moreover, we found that the PuREBP1/2 complex is associated with only a handful of mRNAs, and that it is responsible for the observed purine-dependent regulation of NT8 expression.

Limited expression of non-integrating CpG-free plasmid is associated with increased nucleosome enrichment

CpG-free pDNA was reported to facilitate sustained transgene expression with minimal inflammation in vivo as compared to CpG-containing pDNA. However, the expression potential and impact of CpG-free pDNA in in vitro model have never been described. Hence, in this study, we analyzed the transgene expression profiles of CpG-free pDNA in vitro to determine the influence of CpG depletion from the transgene. We found that in contrast to the published in vivo studies, CpG-free pDNA expressed a significantly lower level of luciferase than CpG-rich pDNA in several human cell lines. By comparing novel CpG-free pDNA carrying CpG-free GFP (pZGFP: 0 CpG) to CpG-rich GFP (pRGFP: 60 CpGs), we further showed that the discrepancy was not influenced by external factors such as gene transfer agent, cell species, cell type, and cytotoxicity. Moreover, pZGFP exhibited reduced expression despite having equal gene dosage as pRGFP. Analysis of mRNA distribution revealed that the mRNA export of pZGFP and pRGFP was similar; however, the steady state mRNA level of pZGFP was significantly lower. Upon further investigation, we found that the CpG-free transgene in non-integrating CpG-free pDNA backbone acquired increased nucleosome enrichment as compared with CpG-rich transgene, which may explain the observed reduced level of steady state mRNA. Our findings suggest that nucleosome enrichment could regulate non-integrating CpG-free pDNA expression and has implications on pDNA design.

Molecular and evolutionary processes generating variation in gene expression

Heritable variation in gene expression is common within and between species. This variation arises from mutations that alter the form or function of molecular gene regulatory networks that are then filtered by natural selection. High-throughput methods for introducing mutations and characterizing their cis- and trans-regulatory effects on gene expression (particularly, transcription) are revealing how different molecular mechanisms generate regulatory variation, and studies comparing these mutational effects with variation seen in the wild are teasing apart the role of neutral and non-neutral evolutionary processes. This integration of molecular and evolutionary biology allows us to understand how the variation in gene expression we see today came to be and to predict how it is most likely to evolve in the future.

Transcriptome and translatome co-evolution in mammals

Gene expression programs define shared and species-specific phenotypes, but their evolution remains largely uncharacterized beyond the transcriptome layer¹. Here we report an analysis of the co-evolution of translatomes and transcriptomes using ribosome-profling and matched RNA-sequencing data for three organs (brain, liver and testis) in fve mammals (human, macaque, mouse, opossum and platypus) and a bird (chicken). Our within-species analyses reveal that translational regulation is widespread in the diferent organs, in particular across the spermatogenic cell types of the testis. The between-species divergence in gene expression is around 20% lower at the translatome layer than at the transcriptome layer owing to extensive buffering between the expression layers, which especially preserved old, essential and housekeeping genes. Translational upregulation specifcally counterbalanced global dosage reductions during the evolution of sex chromosomes and the efects of meiotic sex-chromosome inactivation during spermatogenesis. Despite the overall prevalence of bufering, some genes evolved faster at the translatome layer—potentially indicating adaptive changes in expression; testis tissue shows the highest fraction of such genes. Further analyses incorporating mass spectrometry proteomics data establish that the co-evolution of transcriptomes and translatomes is refected at the proteome layer. Together, our work uncovers co-evolutionary patterns and associated selective forces across the expression layers, and provides a resource for understanding their interplay in mammalian organs.

Highly dynamic transcriptional reprogramming and shorter isoform shifts under acute stresses during biological invasions

Phenotypic plasticity has been increasingly recognized for its importance in adaptation to novel environments, and initial rapid plastic response to acute stresses usually serves as the stepping stone for future adaptation. Differential gene expression and alternative splicing have been proposed as two underlying mechanisms for rapid plastic response to environmental stresses. Here, we used an invasive model species, Ciona savignyi, to investigate the temporary plastic changes under temperature stresses on gene expression and alternative splicing. Our results revealed rapid and highly dynamic gene expression reprogramming and alternative splicing switch under acute stresses. Distinct transcriptional response profiles were triggered by two types of temperature stresses, showing resilience recovery and increasing divergence under heat and cold challenges, respectively. Interestingly, alternative exons were more inclined to be skipped under both heat and cold stresses, leading to shorter isoforms but with maintained Open Reading Frames (ORFs). Although similar response patterns were observed between differential gene expression and alternative splicing, low overlap between Differentially Expressed Genes (DEGs) and Differentially Alternative Spliced Genes (DASGs) suggests that distinct gene sets and associated functions should be involved in temperature challenges. Thus, alternative splicing should offer an additional layer of plastic response to environmental challenges. Finally, we identified key plastic genes involved in both gene expression regulation and alternative splicing. The results obtained here shed light on adaptation and accommodation mechanisms during biological invasions, particularly for acute environmental changes at early stages of biological invasions such as transport and introduction.

GC-AG Introns Features in Long Non-coding and Protein-Coding Genes Suggest Their Role in Gene Expression Regulation

Long non-coding RNAs (lncRNAs) are recognized as an important class of regulatory molecules involved in a variety of biological functions. However, the regulatory mechanisms of long non-coding genes expression are still poorly understood. The characterization of the genomic features of lncRNAs is crucial to get insight into their function. In this study, we exploited recent annotations by GENCODE to characterize the genomic and splicing features of long non-coding genes in comparison with protein-coding ones, both in human and mouse. Our analysis highlighted differences between the two classes of genes in terms of their gene architecture. Significant differences in the splice sites usage were observed between long non-coding and protein-coding genes (PCG). While the frequency of non-canonical GC-AG splice junctions represents about 0.8% of total splice sites in PCGs, we identified a significant enrichment of the GC-AG splice sites in long non-coding genes, both in human (3.0%) and mouse (1.9%). In addition, we found a positional bias of GC-AG splice sites being enriched in the first intron in both classes of genes. Moreover, a significant shorter length and weaker donor and acceptor sites were found comparing GC-AG introns to GT-AG introns. Genes containing at least one GC-AG intron were found conserved in many species, more prone to alternative splicing and a functional analysis pointed toward their enrichment in specific biological processes such as DNA repair. Our study shows for the first time that GC-AG introns are mainly associated with lncRNAs and are preferentially located in the first intron. Additionally, we discovered their regulatory potential indicating the existence of a new mechanism of non-coding and PCGs expression regulation.

The complexity of alternative splicing and landscape of tissue-specific expression in lotus (Nelumbo nucifera) unveiled by Illumina- and single-molecule real-time-based RNA-sequencing

Alternative splicing (AS) plays a critical role in regulating different physiological and developmental processes in eukaryotes, by dramatically increasing the diversity of the transcriptome and the proteome. However, the saturation and complexity of AS remain unclear in lotus due to its limitation of rare obtainment of full-length multiple-splice isoforms. In this study, we apply a hybrid assembly strategy by combining single-molecule real-time sequencing and Illumina RNA-seq to get a comprehensive insight into the lotus transcriptomic landscape. We identified 211,802 high-quality full-length non-chimeric reads, with 192,690 non-redundant isoforms, and updated the lotus reference gene model. Moreover, our analysis identified a total of 104,288 AS events from 16,543 genes, with alternative 3ʹ splice-site being the predominant model, following by intron retention. By exploring tissue datasets, 370 tissue-specific AS events were identified among 12 tissues. Both the tissue-specific genes and isoforms might play important roles in tissue or organ development, and are suitable for ‘ABCE’ model partly in floral tissues. A large number of AS events and isoform variants identified in our study enhance the understanding of transcriptional diversity in lotus, and provide valuable resource for further functional genomic studies.

Nonsense-mediated mRNA decay: The challenge of telling right from wrong in a complex transcriptome

The nonsense-mediated mRNA decay pathway selects and degrades its targets using a dense network of RNA-protein and protein-protein interactions. Together, these interactions allow the pathway to collect copious information about the translating mRNA, including translation termination status, splice junction positions, mRNP composition, and 3'UTR length and structure. The core NMD machinery, centered on the RNA helicase UPF1, integrates this information to determine the efficiency of decay. A picture of NMD is emerging in which many factors contribute to the dynamics of decay complex assembly and disassembly, thereby influencing the probability of decay. The ability of the NMD pathway to recognize mRNP features of diverse potential substrates allows it to simultaneously perform quality control and regulatory functions. In vertebrates, increased transcriptome complexity requires balance between these two functions since high NMD efficiency is desirable for maintenance of quality control fidelity but may impair expression of normal mRNAs. NMD has adapted to this challenge by employing mechanisms to enhance identification of certain potential substrates, while using sequence-specific RNA-binding proteins to shield others from detection. These elaborations on the conserved NMD mechanism permit more sensitive post-transcriptional gene regulation but can have severe deleterious consequences, including the failure to degrade pathogenic aberrant mRNAs in many B cell lymphomas. This article is categorized under: RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Turnover and Surveillance > Turnover/Surveillance Mechanisms.

Towards Better Understanding of KSHV Life Cycle: from Transcription and Posttranscriptional Regulations to Pathogenesis

Kaposi’s sarcoma-associated herpesvirus (KSHV), also known as human herpesvirus-8 (HHV-8), is etiologically linked to the development of Kaposi’s sarcoma, primary effusion lymphoma, and multicentric Castleman’s disease. These malignancies often occur in immunosuppressed individuals, making KSHV infection-associated diseases an increasing global health concern with persistence of the AIDS epidemic. KSHV exhibits biphasic life cycles between latent and lytic infection and extensive transcriptional and posttranscriptional regulation of gene expression. As a member of the herpesvirus family, KSHV has evolved many strategies to evade the host immune response, which help the virus establish a successful lifelong infection. In this review, we summarize the current research status on the biology of latent and lytic viral infection, the regulation of viral life cycles and the related pathogenesis.

Writing a wrong: Coupled RNA polymerase II transcription and RNA quality control

Processing and maturation of precursor RNA species is coupled to RNA polymerase II transcription. Co-transcriptional RNA processing helps to ensure efficient and proper capping, splicing, and 3' end processing of different RNA species to help ensure quality control of the transcriptome. Many improperly processed transcripts are not exported from the nucleus, are restricted to the site of transcription, and are in some cases degraded, which helps to limit any possibility of aberrant RNA causing harm to cellular health. These critical quality control pathways are regulated by the highly dynamic protein-protein interaction network at the site of transcription. Recent work has further revealed the extent to which the processes of transcription and RNA processing and quality control are integrated, and how critically their coupling relies upon the dynamic protein interactions that take place co-transcriptionally. This review focuses specifically on the intricate balance between 3' end processing and RNA decay during transcription termination. This article is categorized under: RNA Turnover and Surveillance > Turnover/Surveillance Mechanisms RNA Processing > 3' End Processing RNA Processing > Splicing Mechanisms RNA Processing > Capping and 5' End Modifications.

Emerging Roles of LSM Complexes in Posttranscriptional Regulation of Plant Response to Abiotic Stress

It has long been assumed that the wide reprogramming of gene expression that modulates plant response to unfavorable environmental conditions is mainly controlled at the transcriptional level. A growing body of evidence, however, indicates that posttranscriptional regulatory mechanisms also play a relevant role in this control. Thus, the LSMs, a family of proteins involved in mRNA metabolism highly conserved in eukaryotes, have emerged as prominent regulators of plant tolerance to abiotic stress. Arabidopsis contains two main LSM ring-shaped heteroheptameric complexes, LSM1-7 and LSM2-8, with different subcellular localization and function. The LSM1-7 ring is part of the cytoplasmic decapping complex that regulates mRNA stability. On the other hand, the LSM2-8 complex accumulates in the nucleus to ensure appropriate levels of U6 snRNA and, therefore, correct pre-mRNA splicing. Recent studies reported unexpected results that led to a fundamental change in the assumed consideration that LSM complexes are mere components of the mRNA decapping and splicing cellular machineries. Indeed, these data have demonstrated that LSM1-7 and LSM2-8 rings operate in Arabidopsis by selecting specific RNA targets, depending on the environmental conditions. This specificity allows them to actively imposing particular gene expression patterns that fine-tune plant responses to abiotic stresses. In this review, we will summarize current and past knowledge on the role of LSM rings in modulating plant physiology, with special focus on their function in abiotic stress responses.