In recent years, tetracyclines, such as doxycycline, have become broadly used to control gene expression by virtue of the Tet-on/Tet-off systems. However, the wide range of direct effects of tetracycline use has not been fully appreciated. We show here that these antibiotics induce a mitonuclear protein imbalance through their effects on mitochondrial translation, an effect that likely reflects the evolutionary relationship between mitochondria and proteobacteria. Even at low concentrations, tetracyclines induce mitochondrial proteotoxic stress, leading to changes in nuclear gene expression and altered mitochondrial dynamics and function in commonly used cell types, as well as worms, flies, mice, and plants. Given that tetracyclines are so widely applied in research, scientists should be aware of their potentially confounding effects on experimental results. Furthermore, these results caution against extensive use of tetracyclines in livestock due to potential downstream impacts on the environment and human health.

The Drosophila melanogaster Genetic Reference Panel (DGRP) is a community resource of 205 sequenced inbred lines, derived to improve our understanding of the effects of naturally occurring genetic variation on molecular and organismal phenotypes. We used an integrated genotyping strategy to identify 4,853,802 single nucleotide polymorphisms (SNPs) and 1,296,080 non-SNP variants. Our molecular population genomic analyses show higher deletion than insertion mutation rates and stronger purifying selection on deletions. Weaker selection on insertions than deletions is consistent with our observed distribution of genome size determined by flow cytometry, which is skewed toward larger genomes. Insertion/deletion and single nucleotide polymorphisms are positively correlated with each other and with local recombination, suggesting that their nonrandom distributions are due to hitchhiking and background selection. Our cytogenetic analysis identified 16 polymorphic inversions in the DGRP. Common inverted and standard karyotypes are genetically divergent and account for most of the variation in relatedness among the DGRP lines. Intriguingly, variation in genome size and many quantitative traits are significantly associated with inversions. Approximately 50% of the DGRP lines are infected with Wolbachia, and four lines have germline insertions of Wolbachia sequences, but effects of Wolbachia infection on quantitative traits are rarely significant. The DGRP complements ongoing efforts to functionally annotate the Drosophila genome. Indeed, 15% of all D. melanogaster genes segregate for potentially damaged proteins in the DGRP, and genome-wide analyses of quantitative traits identify novel candidate genes. The DGRP lines, sequence data, genotypes, quality scores, phenotypes, and analysis and visualization tools are publicly available.

The Anaplastic Lymphoma Kinase (Alk) receptor tyrosine kinase (RTK) plays a critical role in the specification of founder cells (FCs) in the Drosophila visceral mesoderm (VM) during embryogenesis. Reporter gene and CRISPR/Cas9 deletion analysis reveals enhancer regions in and upstream of the Alk locus that influence tissue-specific expression in the amnioserosa (AS), the VM and the epidermis. By performing high throughput yeast one-hybrid screens (Y1H) with a library of Drosophila transcription factors (TFs) we identify Odd-paired (Opa), the Drosophila homologue of the vertebrate Zic family of TFs, as a novel regulator of embryonic Alk expression. Further characterization identifies evolutionarily conserved Opa-binding cis-regulatory motifs in one of the Alk associated enhancer elements. Employing Alk reporter lines as well as CRISPR/Cas9-mediated removal of regulatory elements in the Alk locus, we show modulation of Alk expression by Opa in the embryonic AS, epidermis and VM. In addition, we identify enhancer elements that integrate input from additional TFs, such as Binou (Bin) and Bagpipe (Bap), to regulate VM expression of Alk in a combinatorial manner. Taken together, our data show that the Opa zinc finger TF is a novel regulator of embryonic Alk expression.

Single-cell RNA-sequencing (scRNA-seq) allows whole transcriptome profiling of thousands of individual cells, enabling the molecular exploration of tissues at the cellular level. Such analytical capacity is of great interest to many research groups in the world, yet these groups often lack the expertise to handle complex scRNA-seq datasets.

Gene families expand by gene duplication, and resulting paralogs diverge through mutation. Functional diversification can include neofunctionalization as well as subfunctionalization of ancestral functions. In addition, redundancy in which multiple genes fulfill overlapping functions is often maintained. Here, we use the family of 40 Caenorhabditis elegans insulins to gain insight into the balance between specificity and redundancy. The insulin/insulin-like growth factor (IIS) pathway comprises a single receptor, DAF-2. To date, no single insulin-like peptide recapitulates all DAF-2-associated phenotypes, likely due to redundancy between insulin-like genes. To provide a first-level annotation of potential patterns of redundancy, we comprehensively delineate the spatiotemporal and conditional expression of all 40 insulins in living animals. We observe extensive dynamics in expression that can explain the lack of simple patterns of pairwise redundancy. We propose a model in which gene families evolve to attain differential alliances in different tissues and in response to a range of environmental stresses.

Cytokine signaling is responsible for coordinating conserved epithelial regeneration and immune responses in the digestive tract. In the Drosophila midgut, Upd3 is a major cytokine, which is induced in enterocytes (EC) and enteroblasts (EB) upon oral infection, and initiates intestinal stem cell (ISC) dependent tissue repair. To date, the genetic network directing upd3 transcription remains largely uncharacterized. Here, we have identified the key infection-responsive enhancers of the upd3 gene and show that distinct enhancers respond to various stresses. Furthermore, through functional genetic screening, bioinformatic analyses and yeast one-hybrid screening, we determined that the transcription factors Scalloped (Sd), Mothers against dpp (Mad), and D-Fos are principal regulators of upd3 expression. Our study demonstrates that upd3 transcription in the gut is regulated by the activation of multiple pathways, including the Hippo, TGF-β/Dpp, and Src, as well as p38-dependent MAPK pathways. Thus, these essential pathways, which are known to control ISC proliferation cell-autonomously, are also activated in ECs to promote tissue turnover the regulation of upd3 transcription.

Resistance to enteric pathogens is a complex trait at the crossroads of multiple biological processes. We have previously shown in the Drosophila Genetic Reference Panel (DGRP) that resistance to infection is highly heritable, but our understanding of how the effects of genetic variants affect different molecular mechanisms to determine gut immunocompetence is still limited.

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.

Non-coding variants coordinate transcription factor (TF) binding and chromatin mark enrichment changes over regions spanning >100 kb. These molecularly coordinated regions are named "variable chromatin modules" (VCMs), providing a conceptual framework of how regulatory variation might shape complex traits. To better understand the molecular mechanisms underlying VCM formation, here, we mechanistically dissect a VCM-modulating noncoding variant that is associated with reduced chronic lymphocytic leukemia (CLL) predisposition and disease progression. This common, germline variant constitutes a 5-bp indel that controls the activity of an AXIN2 gene-linked VCM by creating a MEF2 binding site, which, upon binding, activates a super-enhancer-like regulatory element. This triggers a large change in TF binding activity and chromatin state at an enhancer cluster spanning >150 kb, coinciding with subtle, long-range chromatin compaction and robust AXIN2 up-regulation. Our results support a model in which the indel acts as an AXIN2 VCM-activating TF nucleation event, which modulates CLL pathology.

Single-cell transcriptomics (scRNA-seq) has greatly advanced our ability to characterize cellular heterogeneity1. However, scRNA-seq requires lysing cells, which impedes further molecular or functional analyses on the same cells. Here, we established Live-seq, a single-cell transcriptome profiling approach that preserves cell viability during RNA extraction using fluidic force microscopy2,3, thus allowing to couple a cell's ground-state transcriptome to its downstream molecular or phenotypic behaviour. To benchmark Live-seq, we used cell growth, functional responses and whole-cell transcriptome read-outs to demonstrate that Live-seq can accurately stratify diverse cell types and states without inducing major cellular perturbations. As a proof of concept, we show that Live-seq can be used to directly map a cell's trajectory by sequentially profiling the transcriptomes of individual macrophages before and after lipopolysaccharide (LPS) stimulation, and of adipose stromal cells pre- and post-differentiation. In addition, we demonstrate that Live-seq can function as a transcriptomic recorder by preregistering the transcriptomes of individual macrophages that were subsequently monitored by time-lapse imaging after LPS exposure. This enabled the unsupervised, genome-wide ranking of genes on the basis of their ability to affect macrophage LPS response heterogeneity, revealing basal Nfkbia expression level and cell cycle state as important phenotypic determinants, which we experimentally validated. Thus, Live-seq can address a broad range of biological questions by transforming scRNA-seq from an end-point to a temporal analysis approach.

The BXD genetic reference population is a recombinant inbred panel descended from crosses between the C57BL/6 (B6) and DBA/2 (D2) strains of mice, which segregate for about 5 million sequence variants. Recently, some of these variants have been established with effects on general metabolic phenotypes such as glucose response and bone strength. Here we phenotype 43 BXD strains and observe they have large variation (-5-fold) in their spontaneous activity during waking hours. QTL analyses indicate that -40% of this variance is attributable to a narrow locus containing the aryl hydrocarbon receptor (Ahr), a basic helix-loop-helix transcription factor with well-established roles in development and xenobiotic metabolism. Strains with the D2 allele of Ahr have reduced gene expression compared to those with the B6 allele, and have significantly higher spontaneous activity. This effect was also observed in B6 mice with a congenic D2 Ahr interval, and in B6 mice with a humanized AHR allele which, like the D2 allele, is expressed much less and has less enzymatic activity than the B6 allele. Ahr is highly conserved in invertebrates, and strikingly inhibition of its orthologs in D. melanogaster and C. elegans (spineless and ahr-1) leads to marked increases in basal activity. In mammals, Ahr has numerous ligands, but most are either non-selective (e.g. resveratrol) or highly toxic (e.g., 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)). Thus, we chose to examine a major environmental influence--long term feeding with high fat diet (HFD)--to see if the effects of Ahr are dependent on major metabolic differences. Interestingly, while HFD robustly halved movement across all strains, the QTL position and effects of Ahr remained unchanged, indicating that the effects are independent. The highly consistent effects of Ahr on movement indicate that changes in its constitutive activity have a role on spontaneous movement and may influence human behavior.

Alterations in glucose metabolism and epithelial-mesenchymal transition (EMT) constitute two important characteristics of carcinoma progression toward invasive cancer. Despite an extensive characterization of each of them separately, the links between EMT and glucose metabolism of tumor cells remain elusive. Here we show that the neuronal glucose transporter GLUT3 contributes to glucose uptake and proliferation of lung tumor cells that have undergone an EMT.

Gut immunocompetence involves immune, stress and regenerative processes. To investigate the determinants underlying inter-individual variation in gut immunocompetence, we perform enteric infection of 140 Drosophila lines with the entomopathogenic bacterium Pseudomonas entomophila and observe extensive variation in survival. Using genome-wide association analysis, we identify several novel immune modulators. Transcriptional profiling further shows that the intestinal molecular state differs between resistant and susceptible lines, already before infection, with one transcriptional module involving genes linked to reactive oxygen species (ROS) metabolism contributing to this difference. This genetic and molecular variation is physiologically manifested in lower ROS activity, lower susceptibility to ROS-inducing agent, faster pathogen clearance and higher stem cell activity in resistant versus susceptible lines. This study provides novel insights into the determinants underlying population-level variability in gut immunocompetence, revealing how relatively minor, but systematic genetic and transcriptional variation can mediate overt physiological differences that determine enteric infection susceptibility.

Development of eye tissue is initiated by a conserved set of transcription factors termed retinal determination network (RDN). In the fruit fly Drosophila melanogaster, the zinc-finger transcription factor Glass acts directly downstream of the RDN to control identity of photoreceptor as well as non-photoreceptor cells. Tight control of spatial and temporal gene expression is a critical feature during development, cell-fate determination as well as maintenance of differentiated tissues. The molecular mechanisms that control expression of glass, however, remain largely unknown. We here identify complex regulatory mechanisms controlling expression of the glass locus. All information to recapitulate glass expression are contained in a compact 5.2 kb cis-acting genomic element by combining different cell-type specific and general enhancers with repressor elements. Moreover, the immature RNA of the locus contains an alternative small open reading frame (smORF) upstream of the actual glass translation start, resulting in a small peptide instead of the three possible Glass protein isoforms. CRISPR/Cas9-based mutagenesis shows that the smORF is not required for the formation of functioning photoreceptors, but is able to attenuate effects of glass misexpression. Furthermore, editing the genome to generate glass loci eliminating either one or two isoforms shows that only one of the three proteins is critical for formation of functioning photoreceptors, while removing the two other isoforms did not cause defects in developmental or photoreceptor function. Our results show that eye development and function is largely unaffected by targeted manipulations of critical features of the glass transcript, suggesting a strong selection pressure to allow the formation of a functioning eye.

Combinatorial interactions among transcription factors (TFs) play essential roles in generating gene expression specificity and diversity in metazoans. Using yeast 2-hybrid (Y2H) assays on nearly all sequence-specific Drosophila TFs, we identified 1,983 protein-protein interactions (PPIs), more than doubling the number of currently known PPIs among Drosophila TFs. For quality assessment, we validated a subset of our interactions using MITOMI and bimolecular fluorescence complementation assays. We combined our interactome with prior PPI data to generate an integrated Drosophila TF-TF binary interaction network. Our analysis of ChIP-seq data, integrating PPI and gene expression information, uncovered different modes by which interacting TFs are recruited to DNA. We further demonstrate the utility of our Drosophila interactome in shedding light on human TF-TF interactions. This study reveals how TFs interact to bind regulatory elements in vivo and serves as a resource of Drosophila TF-TF binary PPIs for understanding tissue-specific gene regulation.

Brown adipocytes regulate energy expenditure via mitochondrial uncoupling, which makes them attractive therapeutic targets to tackle obesity. However, the regulatory mechanisms underlying brown adipogenesis are still poorly understood. To address this, we profiled the transcriptome and chromatin state during mouse brown fat cell differentiation, revealing extensive gene expression changes and chromatin remodeling, especially during the first day post-differentiation. To identify putatively causal regulators, we performed transcription factor binding site overrepresentation analyses in active chromatin regions and prioritized factors based on their expression correlation with the bona-fide brown adipogenic marker Ucp1 across multiple mouse and human datasets. Using loss-of-function assays, we evaluated both the phenotypic effect as well as the transcriptomic impact of several putative regulators on the differentiation process, uncovering ZFP467, HOXA4 and Nuclear Factor I A (NFIA) as novel transcriptional regulators. Of these, NFIA emerged as the regulator yielding the strongest molecular and cellular phenotypes. To examine its regulatory function, we profiled the genomic localization of NFIA, identifying it as a key early regulator of terminal brown fat cell differentiation.

The enormous variation in human lifespan is in part due to a myriad of sequence variants, only a few of which have been revealed to date. Since many life-shortening events are related to diseases, we developed a Mendelian randomization-based method combining 58 disease-related GWA studies to derive longevity priors for all HapMap SNPs. A Bayesian association scan, informed by these priors, for parental age of death in the UK Biobank study (n=116,279) revealed 16 independent SNPs with significant Bayes factor at a 5% false discovery rate (FDR). Eleven of them replicate (5% FDR) in five independent longevity studies combined; all but three are depleted of the life-shortening alleles in older Biobank participants. Further analysis revealed that brain expression levels of nearby genes (RBM6, SULT1A1 and CHRNA5) might be causally implicated in longevity. Gene expression and caloric restriction experiments in model organisms confirm the conserved role for RBM6 and SULT1A1 in modulating lifespan.

Natural genetic variation affects circadian rhythms across the evolutionary tree, but the underlying molecular mechanisms are poorly understood. We investigated population-level, molecular circadian clock variation by generating >700 tissue-specific transcriptomes of Drosophila melanogaster (w1118 ) and 141 Drosophila Genetic Reference Panel (DGRP) lines. This comprehensive circadian gene expression atlas contains >1700 cycling genes including previously unknown central circadian clock components and tissue-specific regulators. Furthermore, >30% of DGRP lines exhibited aberrant circadian gene expression, revealing abundant genetic variation-mediated, intertissue circadian expression desynchrony. Genetic analysis of one line with the strongest deviating circadian expression uncovered a novel cry mutation that, as shown by protein structural modeling and brain immunohistochemistry, disrupts the light-driven flavin adenine dinucleotide cofactor photoreduction, providing in vivo support for the importance of this conserved photoentrainment mechanism. Together, our study revealed pervasive tissue-specific circadian expression variation with genetic variants acting upon tissue-specific regulatory networks to generate local gene expression oscillations.

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.

Since its emergence in December 2019, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has spread globally and become a major public health burden. Despite its close phylogenetic relationship to SARS-CoV, SARS-CoV-2 exhibits increased human-to-human transmission dynamics, likely due to efficient early replication in the upper respiratory epithelium of infected individuals. Since different temperatures encountered in the human upper and lower respiratory tract (33°C and 37°C, respectively) have been shown to affect the replication kinetics of several respiratory viruses, as well as host innate immune response dynamics, we investigated the impact of temperature on SARS-CoV-2 and SARS-CoV infection using the primary human airway epithelial cell culture model. SARS-CoV-2, in contrast to SARS-CoV, replicated to higher titers when infections were performed at 33°C rather than 37°C. Although both viruses were highly sensitive to type I and type III interferon pretreatment, a detailed time-resolved transcriptome analysis revealed temperature-dependent interferon and pro-inflammatory responses induced by SARS-CoV-2 that were inversely proportional to its replication efficiency at 33°C or 37°C. These data provide crucial insight on pivotal virus-host interaction dynamics and are in line with characteristic clinical features of SARS-CoV-2 and SARS-CoV, as well as their respective transmission efficiencies.