Detecting host responses to microbial stimulation using primary epithelial organoids.

The intestinal epithelium is constantly exposed to microbes residing in the lumen. Traditionally, the response to microbial interactions has been studied in cell lines derived from cancerous tissues, e.g. Caco-2. It is, however, unclear how the responses in these cancer cell lines reflect the responses of a normal epithelium and whether there might be microbial strain-specific effects. To address these questions, we derived organoids from the small intestine from a cohort of healthy individuals. Culturing intestinal epithelium on a flat laminin matrix induced their differentiation, facilitating analysis of microbial responses via the apical membrane normally exposed to the luminal content. Here, it was evident that the healthy epithelium across multiple individuals (n = 9) demonstrates robust acute both common and strain-specific responses to a range of probiotic bacterial strains (BB-12Ⓡ, LGGⓇ, DSM33361, and Bif195). Importantly, parallel experiments using the Caco-2 cell line provide no acute response. Collectively, we demonstrate that primary epithelial cells maintained as organoids represent a valuable resource for assessing interactions between the epithelium and luminal microbes across individuals, and that these models are likely to contribute to a better understanding of host microbe interactions.

Transcriptional and epigenomic profiling identifies YAP signaling as a key regulator of intestinal epithelium maturation.

During intestinal organogenesis, equipotent epithelial progenitors mature into phenotypically distinct stem cells that are responsible for lifelong maintenance of the tissue. While the morphological changes associated with the transition are well characterized, the molecular mechanisms underpinning the maturation process are not fully understood. Here, we leverage intestinal organoid cultures to profile transcriptional, chromatin accessibility, DNA methylation, and three-dimensional (3D) chromatin conformation landscapes in fetal and adult epithelial cells. We observed prominent differences in gene expression and enhancer activity, which are accompanied by local changes in 3D organization, DNA accessibility, and methylation between the two cellular states. Using integrative analyses, we identified sustained Yes-Associated Protein (YAP) transcriptional activity as a major gatekeeper of the immature fetal state. We found the YAP-associated transcriptional network to be regulated at various levels of chromatin organization and likely to be coordinated by changes in extracellular matrix composition. Together, our work highlights the value of unbiased profiling of regulatory landscapes for the identification of key mechanisms underlying tissue maturation.

Promoter sequence and architecture determine expression variability and confer robustness to genetic variants.

Genetic and environmental exposures cause variability in gene expression. Although most genes are affected in a population, their effect sizes vary greatly, indicating the existence of regulatory mechanisms that could amplify or attenuate expression variability. Here, we investigate the relationship between the sequence and transcription start site architectures of promoters and their expression variability across human individuals. We find that expression variability can be largely explained by a promoter's DNA sequence and its binding sites for specific transcription factors. We show that promoter expression variability reflects the biological process of a gene, demonstrating a selective trade-off between stability for metabolic genes and plasticity for responsive genes and those involved in signaling. Promoters with a rigid transcription start site architecture are more prone to have variable expression and to be associated with genetic variants with large effect sizes, while a flexible usage of transcription start sites within a promoter attenuates expression variability and limits genotypic effects. Our work provides insights into the variable nature of responsive genes and reveals a novel mechanism for supplying transcriptional and mutational robustness to essential genes through multiple transcription start site regions within a promoter.

PAMP-triggered genetic reprogramming involves widespread alternative transcription initiation and an immediate transcription factor wave.

Immune responses triggered by pathogen-associated molecular patterns (PAMPs) are key to pathogen defense, but drivers and stabilizers of the growth-to-defense genetic reprogramming remain incompletely understood in plants. Here, we report a time-course study of the establishment of PAMP-triggered immunity (PTI) using cap analysis of gene expression. We show that around 15% of all transcription start sites (TSSs) rapidly induced during PTI define alternative transcription initiation events. From these, we identify clear examples of regulatory TSS change via alternative inclusion of target peptides or domains in encoded proteins, or of upstream open reading frames in mRNA leader sequences. We also find that 60% of PAMP response genes respond earlier than previously thought. In particular, a cluster of rapidly and transiently PAMP-induced genes is enriched in transcription factors (TFs) whose functions, previously associated with biological processes as diverse as abiotic stress adaptation and stem cell activity, appear to converge on growth restriction. Furthermore, examples of known potentiators of PTI, in one case under direct mitogen-activated protein kinase control, support the notion that the rapidly induced TFs could constitute direct links to PTI signaling pathways and drive gene expression changes underlying establishment of the immune state.

Personalized B cell response to the Lactobacillus rhamnosus GG probiotic in healthy human subjects: a randomized trial.

The specific effects of administering live probiotics in the human gut are not well characterized. To this end, we investigated the immediate effect of Lactobacillus rhamnosus GG (LGG) in the jejunum of 27 healthy volunteers 2 h after ingestion using a combination of global RNA sequencing of human biopsies and bacterial DNA sequencing in a multi-visit, randomized, cross-over design (ClinicalTrials.gov number NCT03140878). While LGG was detectable in jejunum after 2 h in treated subjects, the gene expression response vs. placebo was subtle if assessed across all subjects. However, clustering analysis revealed that one-third of subjects exhibited a strong and consistent LGG response involving hundreds of genes, where genes related to B cell activation were upregulated, consistent with prior results in mice. Immunohistochemistry and single cell-based deconvolution analyses showed that this B cell signature likely is due to activation and proliferation of existing B cells rather than B cell immigration to the tissue. Our results indicate that the LGG strain has an immediate effect in the human gut in a subpopulation of individuals. In extension, our data strongly suggest that studies on in vivo probiotic effects in humans require large cohorts and must take individual variation into account.

Characterization of Arabidopsis thaliana Promoter Bidirectionality and Antisense RNAs by Inactivation of Nuclear RNA Decay Pathways.

In animals, RNA polymerase II initiates transcription bidirectionally from gene promoters to produce pre-mRNAs on the forward strand and promoter upstream transcripts (PROMPTs) on the reverse strand. PROMPTs are degraded by the nuclear exosome. Previous studies based on nascent RNA approaches concluded that Arabidopsis (Arabidopsis thaliana) does not produce PROMPTs. Here, we used steady-state RNA sequencing in mutants defective in nuclear RNA decay including the exosome to reassess the existence of Arabidopsis PROMPTs. While they are rare, we identified ∼100 cases of exosome-sensitive PROMPTs in Arabidopsis. Such PROMPTs are sources of small interfering RNAs in exosome-deficient mutants, perhaps explaining why plants have evolved mechanisms to suppress PROMPTs. In addition, we found ∼200 long, unspliced and exosome-sensitive antisense RNAs that arise from transcription start sites within parts of the genome encoding 3'-untranslated regions on the sense strand. The previously characterized noncoding RNA that regulates expression of the key seed dormancy regulator, DELAY OF GERMINATION1, is a typical representative of this class of RNAs. Transcription factor genes are overrepresented among loci with exosome-sensitive antisense RNAs, suggesting a potential for widespread control of gene expression via this class of noncoding RNAs. Lastly, we assess the use of alternative promoters in Arabidopsis and compare the accuracy of existing TSS annotations.

Comprehensive profiling of the fission yeast transcription start site activity during stress and media response.

Fission yeast, Schizosaccharomyces pombe, is an attractive model organism for transcriptional and chromatin biology research. Such research is contingent on accurate annotation of transcription start sites (TSSs). However, comprehensive genome-wide maps of TSSs and their usage across commonly applied laboratory conditions and treatments for S. pombe are lacking. To this end, we profiled TSS activity genome-wide in S. pombe cultures exposed to heat shock, nitrogen starvation, hydrogen peroxide and two commonly applied media, YES and EMM2, using Cap Analysis of Gene Expression (CAGE). CAGE-based annotation of TSSs is substantially more accurate than existing PomBase annotation; on average, CAGE TSSs fall 50-75 bp downstream of PomBase TSSs and co-localize with nucleosome boundaries. In contrast to higher eukaryotes, dispersed TSS distributions are not common in S. pombe. Our data recapitulate known S. pombe stress expression response patterns and identify stress- and media-responsive alternative TSSs. Notably, alteration of growth medium induces changes of similar magnitude as some stressors. We show a link between nucleosome occupancy and genetic variation, and that the proximal promoter region is genetically diverse between S. pombe strains. Our detailed TSS map constitutes a central resource for S. pombe gene regulation research.

The RNA exosome contributes to gene expression regulation during stem cell differentiation.

Gene expression programs change during cellular transitions. It is well established that a network of transcription factors and chromatin modifiers regulate RNA levels during embryonic stem cell (ESC) differentiation, but the full impact of post-transcriptional processes remains elusive. While cytoplasmic RNA turnover mechanisms have been implicated in differentiation, the contribution of nuclear RNA decay has not been investigated. Here, we differentiate mouse ESCs, depleted for the ribonucleolytic RNA exosome, into embryoid bodies to determine to which degree RNA abundance in the two states can be attributed to changes in transcription versus RNA decay by the exosome. As a general observation, we find that exosome depletion mainly leads to the stabilization of RNAs from lowly transcribed loci, including several protein-coding genes. Depletion of the nuclear exosome cofactor RBM7 leads to similar effects. In particular, transcripts that are differentially expressed between states tend to be more exosome sensitive in the state where expression is low. We conclude that the RNA exosome contributes to down-regulation of transcripts with disparate expression, often in conjunction with transcriptional down-regulation.

Characterization of the enhancer and promoter landscape of inflammatory bowel disease from human colon biopsies.

Inflammatory bowel disease (IBD) is a chronic intestinal disorder, with two main types: Crohn's disease (CD) and ulcerative colitis (UC), whose molecular pathology is not well understood. The majority of IBD-associated SNPs are located in non-coding regions and are hard to characterize since regulatory regions in IBD are not known. Here we profile transcription start sites (TSSs) and enhancers in the descending colon of 94 IBD patients and controls. IBD-upregulated promoters and enhancers are highly enriched for IBD-associated SNPs and are bound by the same transcription factors. IBD-specific TSSs are associated to genes with roles in both inflammatory cascades and gut epithelia while TSSs distinguishing UC and CD are associated to gut epithelia functions. We find that as few as 35 TSSs can distinguish active CD, UC, and controls with 85% accuracy in an independent cohort. Our data constitute a foundation for understanding the molecular pathology, gene regulation, and genetics of IBD.

A new class of genic nuclear RNA species in Arabidopsis.

Targeting of ArabidopsisPHABULOSA (PHB) mRNA by miR166 has been implicated in gene body methylation at the PHB locus. We report that the PHB locus produces an array of stable nuclear RNA species that are neither polyadenylated nor capped. Their biogenesis requires neither RNA polymerases IV/V nor miR166-guided cleavage. The PHB RNAs are insensitive to mutation of nuclear RNA decay pathways and are conserved in several Brassicaceae species, suggesting functional relevance. Similar RNA species are also produced by another body-methylated locus encoding the miR414 target eIF2. Our data reveal the existence of a new class of genic nuclear RNA species.

Relation between NOD2 genotype and changes in innate signaling in Crohn's disease on mRNA and miRNA levels.

Crohn's disease is associated with an altered innate immune response of pathogenic importance. This altered response can be associated to loss-of-function polymorphisms in the NOD2 (nucleotide-binding oligomerization domain-containing protein 2) gene, but also changes in transcriptional and post-transcriptional regulatory layers, including microRNA activity. Here, we characterized the link between NOD2 genotype and inflammatory-mediated changes in innate signaling by studying transcriptional and post-transcriptional activity in response to NOD2-agonist muramyl dipeptide in monocytes from healthy controls, and Crohn's disease patients with and without NOD2 loss-of-function polymorphisms. We measured the expression of genes and microRNAs in monocytes from these subjects after stimulation with muramyl dipeptide. Gene expression profiles mainly distinguished the actual muramyl dipeptide response, but not the genotype. A hyper-responsive phenotype was found in Crohn's disease patients without NOD2 mutations, characterized by upregulated cytokine receptors and general downregulation of microRNA expression. Conversely, microRNA expression could identify genotype-specific differences between subject groups but exhibited little change upon muramyl dipeptide treatment. Only two microRNAs showed muramyl dipeptide-induced response, including miR-155, which was found to regulate multiple genes and whose host gene was one of the highest muramyl dipeptide responders. miR-155 was upregulated in Crohn's disease patients with NOD2 mutations following lipopolysaccharide and Escherichia coli treatment, but the upregulation was substantially reduced upon muramyl dipeptide treatment. While Crohn's disease patients with NOD2 mutations on average showed a reduced muramyl dipeptide response, the cohort exhibited large individual variance: a small subset had inflammatory responses almost comparable to wild-type patients on both gene and miR-155 regulatory levels.

Identification of Gene Transcription Start Sites and Enhancers Responding to Pulmonary Carbon Nanotube Exposure in Vivo.

Increased use of nanomaterials in industry, medicine, and consumer products has raised concerns over their toxicity. To ensure safe use of nanomaterials, understanding their biological effects at the molecular level is crucial. In particular, the regulatory mechanisms responsible for the cascade of genes activated by nanomaterial exposure are not well-characterized. To this end, we profiled the genome-wide usage of gene transcription start sites and linked active enhancer regions in lungs of C57BL/6 mice 24 h after intratracheal instillation of a single dose of the multiwalled carbon nanotube (MWCNT) Mitsui-7. Our results revealed a massive gene regulatory response, where expression of key inflammatory genes (e.g., Csf3, Il24, and Fgf23) was increased >100-fold 24 h after Mitsui-7 exposure. Many of the Mitsui-7-responsive transcription start sites were alternative transcription start sites for known genes, and the number of alternative transcription start sites used in a given gene was correlated with overall Mitsui-7 response. Strikingly, genes that were up-regulated after Mitsui-7 exposure only through their main annotated transcription start site were linked to inflammatory and defense responses, while genes up-regulated only through alternative transcription start sites were functionally heterogeneous and not inflammation-associated. Furthermore, we identified almost 12 000 active enhancers, many of which were Mitsui-7-responsive, and we identified similarly responding putative target genes. Overall, our study provides the location and activity of Mitsui-7-induced enhancers and transcription start sites, providing a useful resource for targeted experiments elucidating the biological effects of nanomaterials and the identification of biomarkers for early detection of MWCNT-induced inflammation.

On-the-fly selection of cell-specific enhancers, genes, miRNAs and proteins across the human body using SlideBase.

Genomics consortia have produced large datasets profiling the expression of genes, micro-RNAs, enhancers and more across human tissues or cells. There is a need for intuitive tools to select subsets of such data that is the most relevant for specific studies. To this end, we present SlideBase, a web tool which offers a new way of selecting genes, promoters, enhancers and microRNAs that are preferentially expressed/used in a specified set of cells/tissues, based on the use of interactive sliders. With the help of sliders, SlideBase enables users to define custom expression thresholds for individual cell types/tissues, producing sets of genes, enhancers etc. which satisfy these constraints. Changes in slider settings result in simultaneous changes in the selected sets, updated in real time. SlideBase is linked to major databases from genomics consortia, including FANTOM, GTEx, The Human Protein Atlas and BioGPS.Database URL: http://slidebase.binf.ku.dk.

The Adipose Transcriptional Response to Insulin Is Determined by Obesity, Not Insulin Sensitivity.

Metabolically healthy obese subjects display preserved insulin sensitivity and a beneficial white adipose tissue gene expression pattern. However, this observation stems from fasting studies when insulin levels are low. We investigated adipose gene expression by 5'Cap-mRNA sequencing in 17 healthy non-obese (NO), 21 insulin-sensitive severely obese (ISO), and 30 insulin-resistant severely obese (IRO) subjects, before and 2 hr into a hyperinsulinemic euglycemic clamp. ISO and IRO subjects displayed a clear but globally similar transcriptional response to insulin, which differed from the small effects observed in NO subjects. In the obese, 231 genes were altered; 71 were enriched in ISO subjects (e.g., phosphorylation processes), and 52 were enriched in IRO subjects (e.g., cellular stimuli). Common cardio-metabolic risk factors and gender do not influence these findings. This study demonstrates that differences in the acute transcriptional response to insulin are primarily driven by obesity per se, challenging the notion of healthy obese adipose tissue, at least in severe obesity.

Regulating retrotransposon activity through the use of alternative transcription start sites.

Retrotransposons, the ancestors of retroviruses, have the potential for gene disruption and genomic takeover if not kept in check. Paradoxically, although host cells repress these elements by multiple mechanisms, they are transcribed and are even activated under stress conditions. Here, we describe a new mechanism of retrotransposon regulation through transcription start site (TSS) selection by altered nucleosome occupancy. We show that Fun30 chromatin remodelers cooperate to maintain a high level of nucleosome occupancy at retrotransposon-flanking long terminal repeat (LTR) elements. This enforces the use of a downstream TSS and the production of a truncated RNA incapable of reverse transcription and retrotransposition. However, in stressed cells, nucleosome occupancy at LTR elements is reduced, and the TSS shifts to allow for productive transcription. We propose that controlled retrotransposon transcription from a nonproductive TSS allows for rapid stress-induced activation, while preventing uncontrolled transposon activity in the genome.

Nuclear stability and transcriptional directionality separate functionally distinct RNA species.

Mammalian genomes are pervasively transcribed, yielding a complex transcriptome with high variability in composition and cellular abundance. Although recent efforts have identified thousands of new long non-coding (lnc) RNAs and demonstrated a complex transcriptional repertoire produced by protein-coding (pc) genes, limited progress has been made in distinguishing functional RNA from spurious transcription events. This is partly due to present RNA classification, which is typically based on technical rather than biochemical criteria. Here we devise a strategy to systematically categorize human RNAs by their sensitivity to the ribonucleolytic RNA exosome complex and by the nature of their transcription initiation. These measures are surprisingly effective at correctly classifying annotated transcripts, including lncRNAs of known function. The approach also identifies uncharacterized stable lncRNAs, hidden among a vast majority of unstable transcripts. The predictive power of the approach promises to streamline the functional analysis of known and novel RNAs.

Identification of TNF-α-responsive promoters and enhancers in the intestinal epithelial cell model Caco-2.

The Caco-2 cell line is one of the most important in vitro models for enterocytes, and is used to study drug absorption and disease, including inflammatory bowel disease and cancer. In order to use the model optimally, it is necessary to map its functional entities. In this study, we have generated genome-wide maps of active transcription start sites (TSSs), and active enhancers in Caco-2 cells with or without tumour necrosis factor (TNF)-α stimulation to mimic an inflammatory state. We found 520 promoters that significantly changed their usage level upon TNF-α stimulation; of these, 52% are not annotated. A subset of these has the potential to confer change in protein function due to protein domain exclusion. Moreover, we locate 890 transcribed enhancer candidates, where ∼50% are changing in usage after TNF-α stimulation. These enhancers share motif enrichments with similarly responding gene promoters. As a case example, we characterize an enhancer regulating the laminin-5 γ2-chain (LAMC2) gene by nuclear factor (NF)-κB binding. This report is the first to present comprehensive TSS and enhancer maps over Caco-2 cells, and highlights many novel inflammation-specific promoters and enhancers.

An atlas of active enhancers across human cell types and tissues.

Enhancers control the correct temporal and cell-type-specific activation of gene expression in multicellular eukaryotes. Knowing their properties, regulatory activity and targets is crucial to understand the regulation of differentiation and homeostasis. Here we use the FANTOM5 panel of samples, covering the majority of human tissues and cell types, to produce an atlas of active, in vivo-transcribed enhancers. We show that enhancers share properties with CpG-poor messenger RNA promoters but produce bidirectional, exosome-sensitive, relatively short unspliced RNAs, the generation of which is strongly related to enhancer activity. The atlas is used to compare regulatory programs between different cells at unprecedented depth, to identify disease-associated regulatory single nucleotide polymorphisms, and to classify cell-type-specific and ubiquitous enhancers. We further explore the utility of enhancer redundancy, which explains gene expression strength rather than expression patterns. The online FANTOM5 enhancer atlas represents a unique resource for studies on cell-type-specific enhancers and gene regulation.

Molecular constituents of the extracellular matrix in rat liver mounting a hepatic progenitor cell response for tissue repair.

BACKGROUND: Tissue repair in the adult mammalian liver occurs in two distinct processes, referred to as the first and second tiers of defense. We undertook to characterize the changes in molecular constituents of the extracellular matrix when hepatic progenitor cells (HPCs) respond in a second tier of defense to liver injury. RESULTS: We used transcriptional profiling on rat livers responding by a first tier (surgical removal of 70% of the liver mass (PHx protocol)) and a second tier (70% hepatectomy combined with exposure to 2-acetylaminofluorene (AAF/PHx protocol)) of defense to liver injury and compared the transcriptional signatures in untreated rat liver (control) with those from livers of day 1, day 5 and day 9 post hepatectomy in both protocols. Numerous transcripts encoding specific subunits of collagens, laminins, integrins, and various other extracellular matrix structural components were differentially up- or down-modulated (P < 0.01). The levels of a number of transcripts were significantly up-modulated, mainly in the second tier of defense (Agrn, Bgn, Fbn1, Col4a1, Col8a1, Col9a3, Lama5, Lamb1, Lamb2, Itga4, Igtb2, Itgb4, Itgb6, Nid2), and their signal intensities showed a strong or very strong correlation with Krt1-19, a well-established marker of a ductular/HPC reaction. Furthermore, a significant up-modulation and very strong correlation between the transcriptional profiles of Krt1-19 and St14 encoding matriptase, a component of a novel protease system, was found in the second tier of defense. Real-time PCR confirmed the modulation of St14 transcript levels and strong correlation to Krt-19 and also showed a significant up-modulation and strong correlation to Spint1 encoding HAI-1, a cognate inhibitor of matriptase. Immunodetection and three-dimensional reconstructions showed that laminin, Collagen1a1, agrin and nidogen1 surrounded bile ducts, proliferating cholangiocytes, and HPCs in ductular reactions regardless of the nature of defense. Similarly, matriptase and HAI-1 were expressed in cholangiocytes regardless of the tier of defense, but in the second tier of defense, a subpopulation of HPCs in ductular reactions co-expressed HAI-1 and the fetal hepatocyte marker Dlk1. CONCLUSION: Transcriptional profiling and immunodetection, including three-dimensional reconstruction, generated a detailed overview of the extracellular matrix constituents expressed in a second tier of defense to liver injury.

Detection of active matriptase using a biotinylated chloromethyl ketone peptide.

Matriptase is a member of the family of type II transmembrane serine proteases that is essential for development and maintenance of several epithelial tissues. Matriptase is synthesized as a single-chain zymogen precursor that is processed into a two-chain disulfide-linked form dependent on its own catalytic activity leading to the hypothesis that matriptase functions at the pinnacle of several protease induced signal cascades. Matriptase is usually found in either its zymogen form or in a complex with its cognate inhibitor hepatocyte growth factor activator inhibitor 1 (HAI-1), whereas the active non-inhibited form has been difficult to detect. In this study, we have developed an assay to detect enzymatically active non-inhibitor-complexed matriptase by using a biotinylated peptide substrate-based chloromethyl ketone (CMK) inhibitor. Covalently CMK peptide-bound matriptase is detected by streptavidin pull-down and subsequent analysis by Western blotting. This study presents a novel assay for detection of enzymatically active matriptase in living human and murine cells. The assay can be applied to a variety of cell systems and species.