Search

MicroRNA (miRNA) loaded Argonaute (AGO) complexes regulate gene expression via direct base pairing with their mRNA targets. Previous works suggest that up to 60% of mammalian transcripts might be subject to miRNA-mediated regulation, but it remains largely unknown which fraction of these interactions are functional in a specific cellular context. Here, we integrate transcriptome data from a set of miRNA-depleted mouse embryonic stem cell (mESC) lines with published miRNA interaction predictions and AGO-binding profiles. Using this integrative approach, combined with molecular validation data, we present evidence that < 10% of expressed genes are functionally and directly regulated by miRNAs in mESCs. In addition, analyses of the stem cell-specific miR-290-295 cluster target genes identify TFAP4 as an important transcription factor for early development. The extensive datasets developed in this study will support the development of improved predictive models for miRNA-mRNA functional interactions., EMBO Reports, 23 (9), ISSN:1469-221X, ISSN:1469-3178

Transcriptional and translational control are key determinants of gene expression, however, to what extent these two processes can be collectively coordinated is still poorly understood. Here, we use Nanopore long-read sequencing and cap analysis of gene expression (CAGE-seq) to document the landscape of 5 & PRIME; and 3 & PRIME; untranslated region (UTR) isoforms and transcription start sites of epidermal stem cells, wild-type keratinocytes and squamous cell carcinomas. Focusing on squamous cell carcinomas, we show that a small cohort of genes with alternative 5 & PRIME;UTR isoforms exhibit overall increased translational efficiencies and are enriched in ribosomal proteins and splicing factors. By combining polysome fractionations and CAGE-seq, we further characterize two of these UTR isoform genes with identical coding sequences and demonstrate that the underlying transcription start site heterogeneity frequently results in 5 & PRIME; terminal oligopyrimidine (TOP) and pyrimidine-rich translational element (PRTE) motif switches to drive mTORC1-dependent translation of the mRNA. Genome-wide, we show that highly translated squamous cell carcinoma transcripts switch towards increased use of 5 & PRIME;TOP and PRTE motifs, have generally shorter 5 & PRIME;UTRs and expose decreased RNA secondary structures. Notably, we found that the two 5 & PRIME;TOP motif-containing, but not the TOP-less, RPL21 transcript isoforms strongly correlated with overall survival in human head and neck squamous cell carcinoma patients. Our findings warrant isoform-specific analyses in human cancer datasets and suggest that switching between 5 & PRIME;UTR isoforms is an elegant and simple way to alter protein synthesis rates, set their sensitivity to the mTORC1-dependent nutrient-sensing pathway and direct the translational potential of an mRNA by the precise 5 & PRIME;UTR sequence., Oncogene, 24 (9), ISSN:0950-9232, ISSN:1476-5594

Low flows can impact water use and instream ecology. Therefore, reliable predictions of low-flow metrics are crucial. In this study, we assess which catchment characteristics (climate, topography, geology and landcover) can explain the spatial variability of low-flow metrics at two different scales: the regional scale and the small headwater catchment scale. For the regional-scale analysis, we calculated the mean 7-day annual minimum flow (qmin), the mean of the flow that is exceeded 95% of the year (q95), and the master recession constant (C) for 280 independent gauging stations across the Swiss Plateau and the Swiss Alps for the 2000–2018 period. We assessed the relation between 44 catchment characteristics and the three low-flow metrics based on correlation analysis and a random forest model. Low-flow magnitudes across the Swiss Plateau were positively correlated with the fraction of the area covered by sandstone bedrock or alluvium, and with the area that has a slope between 10° and 30°. Across the Swiss Alps, low-flow magnitudes were positively correlated with the fraction of area with slopes between 30° and 60°, and the area with glacial deposits and debris cover. There was good agreement between observations and predictions by the random forest regression model with the top 11 catchment characteristics for both regions: for 80% of the Swiss Plateau catchments and 60% of the Swiss Alpine catchments, we could predict the three low-flow metrics within an error of 30%. The residuals of the regression model, however, varied across short distances, suggesting that local catchment characteristics affect the variability of low-flow metrics. For the local-scale headwater catchments, we conducted 1-day snapshot field campaigns in 16 catchments during low-flow periods in 2015 and 2016. The measurements in these sub-catchments also showed that areas with sandstone bedrock and a good storage-to-river connectivity had above average low-flow magnitudes. Including knowledge on local catchment characteristics may help to improve regional low-flow predictions, however, not all local catchment characteristics were useful descriptors at larger scales., Hydrological Processes, 36 (1), ISSN:0885-6087, ISSN:1099-1085

The Argonaute proteins (AGO) are well-known for their essential role in post-transcriptional gene silencing in the microRNA (miRNA) biogenesis pathway. Only two AGOs (AGO1 and AGO2) are expressed in mouse embryonic stem cells (mESCs). The transcriptome ofAgomutant mESCs revealed a large and specific set of misregulated genes, compared to other miRNA biogenesis factor mutant cells, suggesting additional functions for the AGOs in stem cells. In this study, we endeavored to understand miRNA-independent roles of the AGOs in gene expression regulation through the integration of multiple datasets. Correlation ofAgomutant differential gene expression with ENCODE histone modification data of WT mESCs revealed that affected genes were regulated by the repressive histone modification H3K27me3. We validated this observation by performing chromatin immunoprecipitation followed by sequencing and observed a global loss of H3K27me3 inAgomutant cells. Nevertheless, this reduction explains only a small part of the specific differential gene expression observed inAgomutant mESCs. By integrating chromatin accessibility data in conjunction with prediction of transcription factor binding sites, we identified differential binding for five transcription factors, including KLF4 as a key modulator of more than half of the specific misregulation of gene expression in the absence of AGO proteins. Our findings illustrate that in addition to chromatin state, information about transcription factor binding is more revelatory in understanding the multi-layered mechanism adopted by cells to regulate gene expression. These data also highlight the importance of an integrative approach to unravel the variety of noncanonical functions of AGOs in mESCs.

SUMMARYMicroRNA (miRNA) loaded Argonaute (AGO) complexes regulate gene expression via direct base pairing with their mRNA targets. Current prediction approaches identified that between 20 to 60% of mammalian transcriptomes are regulated by miRNAs, but it remains largely unknown which fraction of these interactions are functional in a specific cellular context. Here, we integrated transcriptome data from a set of miRNA-depleted mouse embryonic stem cell (mESC) lines with published miRNA interaction predictions and AGO-binding profiles. This integrative approach, combined with molecular validation data, identified that only 6% of expressed genes are functionally and directly regulated by miRNAs in mESCs. In addition, analyses of the stem cell-specific miR-290-295 cluster target genes identified TFAP4 as an important transcription factor for early development. The extensive datasets developed in this study will support the development of improved predictive models for miRNA-mRNA functional interactions.

Transcriptional and translational control are key determinants of gene expression, however, to what extent these two processes can be collectively coordinated is still poorly understood. Here we use long-read sequencing to document the 5’and 3’untranslated region (UTR) isoform landscape of epidermal stem cells, wild-type keratinocytes and squamous cell carcinomas. Focusing on squamous cell carcinomas, we show that a small cohort of genes with alternative 5’UTR isoforms exhibit overall increased translational efficiencies and are enriched in ribosomal proteins and splicing factors. These 5’UTR isoforms with identical coding sequences either include or exclude 5’terminal oligopyrimidine (TOP) motifs and result in vastly altered translational efficiencies of the mRNA. Our findings suggest that switching between TOP and non-TOP motif-containing 5’UTR isoforms is an elegant and simple way to alter protein synthesis rates, set their sensitivity to the mTORC1-dependent nutrient-sensing pathway and direct the translational potential of an mRNA by the precise 5’UTR sequence.

Stratification of patients with pancreatic ductal adenocarcinoma (PDAC) remains a key challenge in the field of clinical oncology. No predictive biomarkers have yet been found for any available treatment options. Previously, we identified SERPINB7 as a putative biomarker for PDAC and thus, herein, we aimed to validate our previous findings and assessed the predictive value of SERPINB7. Patients who underwent surgery and received gemcitabine (gem) or gemcitabine plus nab-paclitaxel (gem/nab) as adjuvant therapy, between 2011 and 2017, were included in this study (n = 57). Expression level of SERPINB7 was assessed in tumor tissue by immunohistochemistry (IHC) and RNA in situ hybridization (RNA ISH). Its association with disease-free survival (DFS) and overall survival (OS) was investigated. While IHC did not show any correlation between survival and the protein level of SERPINB7, RNA ISH revealed that expression of SERPINB7 was associated with a poor DFS (P = .01) and OS (P = .002) in the gem group but not in the gem/nab. Adjusted Cox-regression analysis confirmed the independent predictive value of SERPINB7 on OS (P = .006, HR: 3.47; 95% CI: 1.49–8.09) in the gem group. In conclusion, SERPINB7 was identified as the first predictive RNA biomarker for PDAC. This study suggests that patients who expressed SERPINB7 might receive another treatment than gem alone., Translational Oncology, 12 (1)

The fibroblast growth factor (FGF) and the transforming growth factor-&beta, (TGF-&beta, ) pathways are both involved in the maintenance of human embryonic stem cells (hESCs) and regulate the onset of their differentiation. Their converging functions have suggested that these pathways might share a wide range of overlapping targets. Published studies have focused on the long-term effects (24&ndash, 48 h) of FGF and TGF-&beta, inhibition in hESCs, identifying direct and indirect target genes. In this study, we focused on the earliest transcriptome changes occurring between 3 and 9 h after FGF and TGF-&beta, inhibition to identify direct target genes only. Our analysis clearly shows that only a handful of target transcripts are common to both pathways. This is surprising in light of the previous literature, and has implications for models of cell signaling in human pluripotent cells. In addition, we identified STOX2 as a novel primary target of the TGF-&beta, signaling pathway. We show that STOX2 might act as a novel SMAD2/4 cofactor. Taken together, our results provide insights into the effect of cell signaling on the transcription profile of human pluripotent cells

RNA sequencing (RNA-seq) has become a standard procedure to investigate transcriptional changes between conditions and is routinely used in research and clinics. While standard differential expression (DE) analysis between two conditions has been extensively studied, and improved over the past decades, RNA-seq time course (TC) DE analysis algorithms are still in their early stages. In this study, we compare, for the first time, existing TC RNA-seq tools on an extensive simulation data set and validated the best performing tools on published data. Surprisingly, TC tools were outperformed by the classical pairwise comparison approach on short time series (, Briefings in Bioinformatics, 20 (1), ISSN:1467-5463, ISSN:1477-4054

Pancreatic cancer is a dismal disease with a mortality rate almost similar to its incidence rate. To date, there are neither validated predictive nor prognostic biomarkers for this lethal disease. Thus, the aim of the present study was to retrospectively investigate the capability of biochemical parameters and molecular profiles to predict survival of patients with metastatic pancreatic ductal adenocarcinoma (mPDAC) who participated in a phase II clinical trial to test the safety and efficacy of the combination treatment of capecitabine plus nab-paclitaxel. Herein, we investigated the association of 18 biochemical parameters obtained from routine diagnosis and the clinical outcome of the 30 patients enrolled in the clinical trial. Furthermore, we analysed formalin-fixed paraffin-embedded (FFPE) tumour tissue to identify molecular biomarkers via RNA seq and the Illumina TruSeq Amplicon Cancer panel which covers 48 hotspot genes. Our analysis identified SERPINB7 as a novel transcript and a DNA mutation signature that might predict a poor outcome of disease. Moreover, we identified the bilirubin basal level as an independent predictive factor for overall survival in our study cohort., Scientific Reports, 7, ISSN:2045-2322

Phenotypic differences between closely related species are thought to arise primarily from changes in gene expression due to mutations in cis-regulatory sequences (enhancers). However, it has remained unclear how frequently mutations alter enhancer activity or create functional enhancers de novo. Here we use STARR-seq, a recently developed quantitative enhancer assay, to determine genome-wide enhancer activity profiles for five Drosophila species in the constant trans-regulatory environment of Drosophila melanogaster S2 cells. We find that the functions of a large fraction of D. melanogaster enhancers are conserved for their orthologous sequences owing to selection and stabilizing turnover of transcription factor motifs. Moreover, hundreds of enhancers have been gained since the D. melanogaster–Drosophila yakuba split about 11 million years ago without apparent adaptive selection and can contribute to changes in gene expression in vivo. Our finding that enhancer activity is often deeply conserved and frequently gained provides functional insights into regulatory evolution.

Neuronal L-type voltage-gated calcium channels (LTCCs) are involved in several physiological functions, but increased activity of LTCCs has been linked to pathology. Due to the coupling of LTCC-mediated Ca2+ influx to Ca2+-dependent conductances, such as KCa or non-specific cation channels, LTCCs act as important regulators of neuronal excitability. Augmentation of after-hyperpolarizations may be one mechanism that shows how elevated LTCC activity can lead to neurological malfunctions. However, little is known about other impacts on electrical discharge activity. We used pharmacological up-regulation of LTCCs to address this issue on primary rat hippocampal neurons. Potentiation of LTCCs with Bay K8644 enhanced excitatory postsynaptic potentials to various degrees and eventually resulted in paroxysmal depolarization shifts (PDS). Under conditions of disturbed Ca2+ homeostasis, PDS were evoked frequently upon LTCC potentiation. Exposing the neurons to oxidative stress using hydrogen peroxide also induced LTCC-dependent PDS. Hence, raising LTCC activity had unidirectional effects on brief electrical signals and increased the likeliness of epileptiform events. However, long-lasting seizure-like activity induced by various pharmacological means was affected by Bay K8644 in a bimodal manner, with increases in one group of neurons and decreases in another group. In each group, isradipine exerted the opposite effect. This suggests that therapeutic reduction in LTCC activity may have little beneficial or even adverse effects on long-lasting abnormal discharge activities. However, our data identify enhanced activity of LTCCs as one precipitating cause of PDS. Because evidence is continuously accumulating that PDS represent important elements in neuropathogenesis, LTCCs may provide valuable targets for neuroprophylactic therapy. Electronic supplementary material The online version of this article (doi:10.1007/s12017-013-8234-1) contains supplementary material, which is available to authorized users.

The Argonaute proteins (AGOs) are well known for their role in post-transcriptional gene silencing in the microRNA (miRNA) pathway. Here we show that in mouse embryonic stem cells, AGO1&2 serve additional functions that go beyond the miRNA pathway. Through the combined deletion of both Agos, we identified a specific set of genes that are uniquely regulated by AGOs but not by the other miRNA biogenesis factors. Deletion of Ago2&1 caused a global reduction of the repressive histone mark H3K27me3 due to downregulation at protein levels of Polycomb repressive complex 2 components. By integrating chromatin accessibility, prediction of transcription factor binding sites, and chromatin immunoprecipitation sequencing data, we identified the pluripotency factor KLF4 as a key modulator of AGO1&2-regulated genes. Our findings revealed a novel axis of gene regulation that is mediated by noncanonical functions of AGO proteins that affect chromatin states and gene expression using mechanisms outside the miRNA pathway., Stem Cell Reports, 17 (5), ISSN:2213-6711

Analysis of gene expression has contributed to a plethora of biological and medical research studies. Microarrays have been intensively used for the profiling of gene expression during diverse developmental processes, treatments and diseases. New massively parallel sequencing methods, often named as RNA-sequencing (RNA-seq) are extensively improving our understanding of gene regulation and signaling networks. Computational methods developed originally for microarrays analysis can now be optimized and applied to genome-wide studies in order to have access to a better comprehension of the whole transcriptome. This review addresses current challenges on RNA-seq analysis and specifically focuses on new bioinformatics tools developed for time series experiments. Furthermore, possible improvements in analysis, data integration as well as future applications of differential expression analysis are discussed., Computational and Structural Biotechnology Journal, 13, ISSN:2001-0370