Enteric infection induces Lark-mediated intron retention at the 5 ' end of Drosophila genes

Background RNA splicing is a key post-transcriptional mechanism that generates protein diversity and contributes to the fine-tuning of gene expression, which may facilitate adaptation to environmental challenges. Here, we employ a systems approach to study alternative splicing changes upon enteric infection in females from classical Drosophila melanogaster strains as well as 38 inbred lines. Results We find that infection leads to extensive differences in isoform ratios, which results in a more diverse transcriptome with longer 5 ' untranslated regions (5 ' UTRs). We establish a role for genetic variation in mediating inter-individual splicing differences, with local splicing quantitative trait loci (local-sQTLs) being preferentially located at the 5 ' end of transcripts and directly upstream of splice donor sites. Moreover, local-sQTLs are more numerous in the infected state, indicating that acute stress unmasks a substantial number of silent genetic variants. We observe a general increase in intron retention concentrated at the 5 ' end of transcripts across multiple strains, whose prevalence scales with the degree of pathogen virulence. The length, GC content, and RNA polymerase II occupancy of these introns with increased retention suggest that they have exon-like characteristics. We further uncover that retained intron sequences are enriched for the Lark/RBM4 RNA binding motif. Interestingly, we find that lark is induced by infection in wild-type flies, its overexpression and knockdown alter survival, and tissue-specific overexpression mimics infection-induced intron retention. Conclusion Our collective findings point to pervasive and consistent RNA splicing changes, partly mediated by Lark/RBM4, as being an important aspect of the gut response to infection.

Diverse cell-specific patterns of alternative polyadenylation in Drosophila

Bart Deplancke, Sandy Lee

Most genes in higher eukaryotes express isoforms with distinct 3' untranslated regions (3' UTRs), generated by alternative polyadenylation (APA). Since 3' UTRs are predominant locations of post-transcriptional regulation, APA can render such programs conditional, and can also alter protein sequences via alternative last exon (ALE) isoforms. We previously used 3'-sequencing from diverse Drosophila samples to define multiple tissue-specific APA landscapes. Here, we exploit comprehensive single nucleus RNA-sequencing data (Fly Cell Atlas) to elucidate cell-type expression of 3' UTRs across >250 adult Drosophila cell types. We reveal the cellular bases of multiple tissue-specific APA/ALE programs, such as 3' UTR lengthening in differentiated neurons and 3' UTR shortening in spermatocytes and spermatids. We trace dynamic 3' UTR patterns across cell lineages, including in the male germline, and discover new APA patterns in the intestinal stem cell lineage. Finally, we correlate expression of RNA binding proteins (RBPs), miRNAs and global levels of cleavage and polyadenylation (CPA) factors in several cell types that exhibit characteristic APA landscapes, yielding candidate regulators of transcriptome complexity. These analyses provide a comprehensive foundation for future investigations of mechanisms and biological impacts of alternative 3' isoforms across the major cell types of this widely-studied model organism.

NATURE PORTFOLIO2022

Phase specific transcriptional regulation of circadian clock and metabolism in mouse liver

Jonathan Aryeh Sobel

The molecular clock has been conserved from cyanobacteria to mammals and is believed to align behavioral and biochemical processes with the diurnal cycle. This cellular mechanism has been an advantage to increase the fitness of organisms through the ability to anticipate food availability or predator presence. Accumulating evidence has revealed that the circadian clock is intimately interconnected with the metabolic cycle. But, the nature of these interconnections is yet not clear at the transcriptional level, and the contribution of cis-regulatory modules has not been elucidated. Accessible chromatin regions of the genome are implicated in diverse processes, such as gene regulation through enhancers and promoters, insulation of genomic domains or alternative splicing. They allow cell-type specific programs that are controlled by tissue-specific transcription factors and chromatin modifiers, although tissue-specific regulation of the circadian clock remains unclear. Therefore, we compared genome-wide BMAL1 binding and chromatin accessibility in the liver and in NIH3T3 fibroblasts. Moreover, for the first time, our study explores the dynamics of accessible regions, histone 3 lysine 27 acetylation (H3K27ac) and RNA polymerase II (Pol II) every 4 hours over one day in mouse liver. In this study, we show that a substantial fraction of these accessible sites are oscillating during a diurnal cycle with a circadian period. We observed that these accessible regions are enriched in the proximity of actively transcribed genes, and that they are dynamically affected by the binding of transcription factors such as BMAL1. We investigated wild-type (WT) and Bmal1-/- genotypes, in night restricted feeding regimen and in Light-Dark (LD) cycle, to study the circadian clock regulatory network underlying diurnal transcription. Therefore, we applied a penalized generalized linear model to infer the activity of transcription factor binding motifs in oscillating accessible sites using Pol II loadings at the transcription start sites of nearby genes. We were able to recapitulate the known regulatory elements of the circadian clock, notably E-box, D-Box, and ROR-responsive elements (RRE) in the WT genotype. On the other hand, we found that Forkhead box (FOX), glucocorticoids responsive elements (GRE), and C-AMP Response Element (CRE) were the main contributors of the regulation of oscillating genes in Bmal1-/-. Finally, using a mixture model to detect footprints with a base pair resolution, we studied the dynamics of the accessibility overlapping E-box. Our last analysis suggested that BMAL1/CLOCK is binding on double E-boxes with a spacer of 6 or 7 bp in a hetero-tetramer configuration. A 3D structure model further supported this binding mode. In Summary, we used DNase I-seq and ChIP-seq of Pol II and H3K27ac to study the circadian chromatin landscape in a 4h time-resolved experimental design. We uncovered the underlying circadian transcriptional regulatory network and, we dissected the chromatin accessibility around BMAL1 binding sites at a base pair resolution, which led to an unappreciated mode of binding of BMAL1/CLOCK in a hetero-tetramer conformation on double E-boxes. Lastly, we found tissue-specific factors that might contribute to tissue-specific binding of BMAL1/CLOCK.

EPFL2016

Elucidating regulatory mechanisms shaping the transcriptome and proteome of the gut immune response in Drosophila melanogaster

Michael Vincent Frochaux

Work on Drosophila melanogaster paved the way to our current understanding of modern genetics. Since then, this model organism has contributed greatly to various fields such as neurobiology, development, and immunology. The discovery and analysis of the various pathways of the Drosophila immune response led the way to the in-depth characterization of the gut immune response. Furthermore, a genome-wide association study on a population of Drosophila identified resistance to infection by the entomopathogenic bacteria Pseudomonas entomophila (P.e.) as a complex phenotype, with highly resistant and susceptible lines, and highlighted genes and pathways mediating inter-individual differences. However, the molecular mechanisms mediating the resistance to enteric infection remain largely uncharacterized. This study uses the same Drosophila population to analyze what are the molecular mechanisms mediating the resistance to enteric infection. We analyzed the genetic determinants of gene expression variation among the most resistant and susceptible lines in response to infection using expression quantitative trait loci (eQTL) analysis. We also characterized the mode of action of these eQTLs into cis- and trans-acting. Furthermore, we identified the gene nutcracker (ntc) as the only differentially-expressed gene between resistance classes. We then demonstrate that loss of function ntc mutants are significantly more susceptible to infection and then identify a single nucleotide polymorphism (SNP) located in a transcription factor binding site (TFBS) which affects the binding affinity of the transcription factor Broad leading to allele-specific expression variation of the gene ntc. If the transcriptomic response to infection has been thoroughly characterized, it is not the case for the proteomic response. We thus sought to characterize for the first time the proteomic response of the Drosophila to enteric infection. We performed mass spectrometry and RNA sequencing on dissected guts from flies infected by two Gram-negative bacteria, Erwinia carotovora carotovora 15 (Ecc15) or Pseudomonas entomophila, 4h and 16h after infection. We found that a large portion of the measurable proteome (12%) varies after infection and that protein changes are strongly time- and infection-dependent. We showed the relatively poor correlation between gene expression and protein abundance in the Drosophila enteric immune response. Finally, we performed a screen of several proteins that were identified in our proteomics work but had previously not been found when assessing the gene expression response to infection using loss-of-function mutants and identifying 7 genes that modulate overall susceptibility to P.e. infection. In summary, this work analyze the genetic determinants of gene expression variation in the DGRP population upon infection and characterize them according to their mode of action. Furthermore, we describe how a non-coding variant lowers resistance to infection by modulating ntc gene expression through altered Broad repressor binding. Finally, it provides the first comprehensive characterization of the Drosophila gut proteome upon infection by Gram-negative bacteria which can be the basis of future proteomic analysis of the enteric immune response