Search

During intestinal organogenesis, equipotent epithelial progenitors mature into phenotypically distinct stem cells that are responsible for life-long 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 3D chromatin conformation landscapes defining fetal and adult epithelial cells. We observed prominent differences in gene expression and enhancer activity, accompanied by changes in 3D organization and local changes in DNA accessibility and methylation, between the two cellular states. Using integrative analyses, we identified sustained 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. Altogether, our work highlights the value of unbiased profiling of regulatory landscapes for the identification of key mechanisms underlying tissue maturation.

Inter-patient variability and the similarity of healthy and leukemic stem cells have impeded the characterization of leukemic stem cells (LSCs) in acute myeloid leukemia (AML), and their differentiation landscape. Here, we introduce CloneTracer, a novel method that adds clonal resolution to single-cell RNA-seq datasets. Applied to samples from 19 AML patients, CloneTracer revealed routes of leukemic differentiation. While residual healthy cells dominated the dormant stem cell compartment, active leukemic stem cells resembled their healthy counterpart and retained erythroid capacity. By contrast, downstream myeloid progenitors were highly aberrant and constituted the disease-defining compartment: Their gene expression and differentiation state determined both chemotherapy response and the leukemia’s ability to differentiate to transcriptomically normal monocytes. Finally, we demonstrated the potential of CloneTracer to identify surface markers mis-regulated specifically in leukemic cells by intra-patient comparisons. Taken together, CloneTracer revealed a differentiation landscape that mimics its healthy counterpart and determines biology and therapy response in AML.

BackgroundAxial Spondyloarthritis can lead to significant disability and impairment in quality of life. TNF inhibitors are recommended to patients enduring active disease despite conventional treatment. Nonetheless, up to 40% of patients of patients fail to respond to TNF inhibitors. In this context, it is important to identify as early as possible patients highly likely to respond. This study aims at identifying, among axial spondyloarthritis patients undergoing treatment with the TNF inhibitor adalimumab, early molecular biomarkers differentiating good responders from non-responders after 14 weeks of treatment, as measured by ASAS20.MethodsPeripheral blood RNA sequencing and serum proteins measured by mass spectrometry were evaluated in a cohort of biologic naïve axial spondyloarthritis patients (n = 35), before (baseline) and after (3-5 days, 2 weeks and 14 weeks) treatment with adalimumab. Results from differential expression analysis were used in combination with clinical data to build logistic regression models and random forest models to predict response to adalimumab at baseline.ResultsResponders to adalimumab presented higher levels of markers of innate immunity at baseline, mostly related with neutrophils, and lower levels of adaptive immunity markers, particularly B-cells. A logistic regression model incorporating ASDAS-CRP and AFF3, the top differentially expressed gene between responders and non-responders at baseline, enabled an accurate prediction of response to adalimumab in our cohort (AUC=0.96), with random forest models suggesting 80% predictive accuracy. A treatment-associated signature suggests a reduction in inflammatory activity, with C-reactive protein and Haptoglobin showing strong and early decrease in the serum of axial spondyloarthritis patients, while a cluster of apolipoproteins showed increased expression at week 14.ConclusionsDifferences in disease activity and/or blood innate/adaptive immune cell type composition at baseline may be a major contributor to response to adalimumab in axial spondyloarthritis, where a model including clinical and blood gene expression variables shows high predictive power. Our results suggest novel molecular biomarkers of response to adalimumab at baseline.Trial registrationAxial spondyloarthritis patients were selected from participants of the Bioefficacy study - Biomarkers Identification of Anti-TNFα Agent’s Efficacy in Ankylosing Spondylitis Patients Using a Transcriptome Analysis and Mass Spectrometry (clinical trials.gov identifier NCT02492217).

Patients with chronic obstructive pulmonary disease (COPD) are still waiting for curative treatments. Considering the environmental cause of COPD (e.g., cigarette smoke) and disease phenotypes, including stem-cell senescence and impaired differentiation, we hypothesized that COPD will be associated with altered epigenetic signaling in lung cells. We generated genome-wide DNA methylation maps at single CpG resolution of primary human lung fibroblasts (HLFs) isolated from distal parenchyma of ex-smoker controls and COPD patients, with both mild and severe disease. The epigenetic landscape is markedly changed in lung fibroblasts across COPD stages, with DNA methylation changes occurring predominantly in regulatory regions, including promoters and enhancers. RNA sequencing of matched fibroblasts demonstrated dysregulation of genes involved in proliferation, DNA repair, and extracellular matrix organization. Notably, we identified epigenetic and transcriptional dysregulation already in mild COPD patients, providing unique insights into early disease. Integration of profiling data identified 110 candidate regulators of disease phenotypes, including epigenetic factors. Using phenotypic screens, we verified the regulator capacity of multiple candidates and linked them to repair processes in the human lung.Our study provides first integrative high-resolution epigenetic and transcriptomic maps of human lung fibroblasts across stages of COPD. We reveal novel transcriptomic and epigenetic signatures associated with COPD onset and progression and identify new candidate regulators involved in the pathogenesis of chronic respiratory diseases. The presence of various epigenetic factors among the candidates demonstrates that epigenetic regulation in COPD is an exciting research field that holds promise for novel therapeutic avenues for patients.

Control of gene expression in specific tissues and/or at certain stages of development allows the study and manipulation of gene function with high precision. Site-specific genome recombination by the Flippase (FLP) and Cre enzymes has proven particularly relevant. Joint efforts of many research groups have led to the creation of efficient FLP and Cre drivers to regulate gene expression in a variety of tissues in Caenorhabditis elegans. Here, we extend this toolkit by the addition of FLP lines that drive recombination specifically in distal tip cells, the somatic gonad, coelomocytes and the epithelial P lineage. In some cases, recombination-mediated gene knockouts do not completely deplete protein levels due to persistence of long-lived proteins. To overcome this, we developed a spatiotemporally regulated degradation system for GFP fusion proteins (GFPdeg) based on FLP-mediated recombination. Using two stable nuclear pore proteins, MEL-28/ELYS and NPP-2/NUP85 as examples, we report the benefit of combining tissue-specific gene knockout and protein degradation to achieve complete protein depletion. We also demonstrate that FLP-mediated recombination can be utilized to identify transcriptomes in a C. elegans tissue of interest. We have adapted RNA polymerase DamID (RAPID) for the FLP toolbox and by focusing on a well-characterized tissue, the hypodermis, we show that the vast majority of genes identified by RAPID are known to be expressed in this tissue. These tools allow combining FLP activity for simultaneous gene inactivation and transcriptomic profiling, thus enabling the inquiry of gene function in various complex biological processes.

Anti-cancer therapies often exhibit only short-term effects. Tumors typically develop drug resistance causing relapses that might be tackled with drug combinations. Identification of the right combination is challenging and would benefit from high-content, high-throughput combinatorial screens directly on patient biopsies. However, such screens require a large amount of material, normally not available from patients. To address these challenges, we developed a scalable microfluidic workflow to screen hundreds of drug combinations in picoliter-size droplets using transcriptome changes as a readout for drug effects. We devised a deterministic combinatorial DNA barcoding approach to encode treatment conditions, enabling the gene expression-based readout of drug effects in a highly multiplexed fashion. We applied our method to screen the effect of 420 drug combinations on the transcriptome of K562 cells using only ∼250 single cell droplets per condition, to successfully predict synergistic and antagonistic drug pairs, as well as their pathway activities.

N-Glycanase 1 (NGLY1) deficiency is a rare and complex genetic disorder. Although recent studies have shed light on the molecular underpinnings of NGLY1 deficiency, a systematic characterization of gene and protein expression changes in patient-derived cells has been lacking. Here, we performed RNA-sequencing and mass spectrometry to determine the transcriptomes and proteomes of 66 cell lines representing 4 different cell types derived from 14 NGLY1 deficient patients and 17 controls. While gene and protein expression levels agreed well with each other, expression differences were more pronounced at the protein level. Although NGLY1 protein levels were up to 9.5-fold downregulated in patients compared to parent controls, depending on the genotype, NGLY1 protein was still detectable in all patient-derived lymphoblastoid cell lines. Consistent with the role of NGLY1 as a regulator of the transcription factor Nrf1, we observed a cell type-independent downregulation of proteasomal genes in NGLY1 deficient cells. In contrast, genes involved in ribosomal mRNA processing were upregulated in multiple cell types. In addition, we observed cell type-specific effects. For example, genes and proteins involved in glutathione synthesis, such as the glutamate-cystein ligase subunits GCLC and GCLM, were downregulated specifically in lymphoblastoid cells. We provide a web application that enables access to all results generated in this study at https://apps.embl.de/ngly1browser. This resource will guide future studies of NGLY1 deficiency in directions that are most relevant to patients.

SummaryVariants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are replacing the initial wild-type strain, jeopardizing current efforts to contain the pandemic. Amino acid exchanges in the spike protein are of particular concern as they can render the virus more transmissible or reduce vaccine efficacy. Here, we conducted whole genome sequencing of SARS-CoV-2 positive samples from the Rhine-Neckar district in Germany during January-March 2021. We detected a total of 166 samples positive for a variant with a distinct mutational pattern in the spike gene comprising L18F, L452R, N501Y, A653V, H655Y, D796Y and G1219V with a later gain of A222V. This variant was designated A.27.RN according to its phylogenetic clade classification. It emerged in parallel with the B.1.1.7 variant, increased to >50% of all SARS-CoV-2 variants by week five. Subsequently it decreased to

Single-cell genomics has transformed our understanding of complex cellular systems. However, excessive costs and a lack of strategies for the purification of newly identified cell types impede their functional characterization and large-scale profiling. Here, we have generated high content single-cell proteo-genomic reference maps of human blood and bone marrow that quantitatively link the expression of up to 197 surface markers to cellular identities and biological processes across all major hematopoietic cell types in healthy aging and leukemia. These reference maps enable the automatic design of cost-effective high-throughput cytometry schemes that outperform state-of-the-art approaches, accurately reflect complex topologies of cellular systems, and permit the purification of precisely defined cell states. The systematic integration of cytometry and proteo-genomic data enables measuring the functional capacities of precisely mapped cell states at the single-cell level. Our study serves as an accessible resource and paves the way for a data-driven era in cytometry.

McQ is a SARS-CoV-2 quantification assay that couples early-stage barcoding with high-throughput sequencing to enable multiplexed processing of thousands of samples. McQ is based on homemade enzymes to enable low-cost testing of large sample pools, circumventing supply chain shortages.Implementation of cost-efficient high-throughput methods for detection of RNA viruses such as SARS-CoV-2 is a potent strategy to curb ongoing and future pandemics. Here we describe Multiplexed SARS-CoV-2 Quantification platform (McQ), an in-expensive scalable framework for SARS-CoV-2 quantification in saliva samples. McQ is based on the parallel sequencing of barcoded amplicons generated from SARS- CoV-2 genomic RNA. McQ uses indexed, target-specific reverse transcription (RT) to generate barcoded cDNA for amplifying viral- and human-specific regions. The barcoding system enables early sample pooling to reduce hands-on time and makes the ap-proach scalable to thousands of samples per sequencing run. Robust and accurate quantification of viral load is achieved by measuring the abundance of Unique Molecular Identifiers (UMIs) introduced during reverse transcription. The use of homemade reverse transcriptase and polymerase enzymes and non-proprietary buffers reduces RNA to library reagent costs to 92 cents/sample and circumvents potential supply chain short-ages. We demonstrate the ability of McQ to robustly quantify various levels of viral RNA in 838 clinical samples and accu-rately diagnose positive and negative control samples in a test-ing workflow entailing self-sampling and automated RNA ex-traction from saliva. The implementation of McQ is modular, scalable and could be extended to other pathogenic targets in future.

Copy Number Gains (CNGs) lead to genetic heterogeneity, driving evolution and carcinogenesis. The mechanisms promoting CNG formation however remain poorly characterized. We show that abnormal expression of the replication licensing factor Cdc18 in fission yeast, which leads to genome-wide re-replication, drives the formation of CNGs at different genomic loci, promoting the acquisition of new selectable traits. Whole genome sequencing reveals Mb long, primarily extrachromosomal amplicons. Genetic analysis shows that homology-mediated repair is required to resolve re-replication intermediates into heritable CNGs. Consistently, we show that in mammalian cells overexpression of CDC6 and/or CDT1 leads to CNGs and promotes drug resistance. In human cells, multiple repair pathways are activated upon rereplication and act antagonistically, with RAD52 promoting and 53BP1 inhibiting CNG formation. In tumours, CDT1 and/or CDC6 overexpression correlates with copy number gains genome-wide. We propose re-replication as an evolutionary-conserved driver of CNGs, highlighting a link between aberrant licensing, CNGs and cancer.

Human intestinal epithelial cells form a primary barrier protecting us from pathogens, yet only limited knowledge is available about individual contribution of each cell type to mounting an immune response against infection. Here, we developed a pipeline combining single-cell RNA-Seq and highly-multiplex RNA imaging and applied it to human intestinal organoids infected with human astrovirus, a model human enteric virus. We found that interferon controls the infection and that astrovirus infects all major cell types and lineages with a preferential infection of proliferating cells. Intriguingly, each intestinal epithelial cell lineage has a unique basal expression of interferon-stimulated genes and, upon astrovirus infection, undergoes an antiviral transcriptional reprogramming by upregulating distinct sets of interferon-stimulated genes. These findings suggest that in the human intestinal epithelium, each cell lineage plays a unique role in resolving virus infection. Our pipeline can be applicable to other organoids and viruses, opening new avenues to unravel roles of individual cell types in viral pathogenesis.

Gene activation requires the cooperative activity of multiple transcription factors at cis-regulatory elements. Yet, most transcription factors have short residence time, questioning the requirement of their physical co-occupancy on DNA to achieve cooperativity. Here, we advance Single Molecule Footprinting to detect individual molecular interactions of transcription factors and nucleosomes with DNA at mouse cis-regulatory elements. We apply this strategy to quantify the simultaneous binding of multiple transcription factors on single DNA molecules. Analysis of the binary occupancy patterns at thousands of motif combinations reveals that for most types of transcription factors high DNA co-occupancy can occur in absence of direct physical interaction, at sites of competition with nucleosomes. Perturbation of pairwise interactions demonstrates the function of molecular co-occupancy for binding cooperativity. These findings elucidate the binding cooperativity mechanism used by transcription factors in absence of strict organisation of their binding motifs, a characteristic feature of most of enhancers.

The ubiquitous redox coenzyme nicotinamide adenine dinucleotide (NAD) acts as a non-canonical cap structure on prokaryotic and eukaryotic ribonucleic acids. Here we find that in budding yeast, NAD-RNAs are abundant (>1400 species), short (, NAD (nicotinamide adenine dinucleotide) acts as a non-canonical RNA cap structure in bacteria and eukaryotes. Here the authors demonstrate the whole landscape of budding yeast NAD-RNAs which are subject to diverse surveillance pathways, suggesting that NAD caps in budding yeast are mostly dysfunctional.

Structural variation (SV), where rearrangements delete, duplicate, invert or translocate DNA segments, is a major source of somatic cell variation. It can arise in rapid bursts, mediate genetic heterogenity, and dysregulate cancer-related pathways. The challenge to systematically discover SVs in single cells remains unsolved, with copy-neutral and complex variants typically escaping detection. We developed single cell tri-channel-processing (scTRIP), a computational framework that jointly integrates read depth, template strand and haplotype phase to comprehensively discover SVs in single cells. We surveyed SV landscapes of 565 single cell genomes, including transformed epithelial cells and patient-derived leukemic samples, and discovered abundant SV classes including inversions, translocations and large-scale genomic rearrangements mediating oncogenic dysregulation. We dissected the ‘molecular karyotype’ of the leukemic samples and examined their clonal structure. Different from prior methods, scTRIP also enabled direct detection and discrimination of SV mutational processes in individual cells, including breakage-fusion-bridge cycles. scTRIP will facilitate studies of clonal evolution, genetic mosaicism and somatic SV formation, and could improve disease classification for precision medicine.

SummarySynergy of resistance and disease tolerance mechanisms is necessary for an effective immune response leading to survival and return to homeostasis when an organism is challenged by infection. Antibiotics are used for their resistance enhancement capabilities by decreasing pathogen load, but several classes have long been known to have beneficial effects that cannot be explained strictly on the basis of their capacity to control the infectious agent. Here we report that tetracycline antibiotics, a class of ribosome-targeting drugs, robustly protects against sepsis by inducing disease tolerance, independently from their direct antibiotic properties. Mechanistically, we find that mitochondrial inhibition of protein synthesis perturbs the electron transfer chain and leads to improved damage repair in the lung and fatty acid oxidation and glucocorticoid sensitivity in the liver. Using a partial and acute deletion ofCRIF1in the liver, a critical mitoribosomal component for protein synthesis, we find that mice are protected against bacterial sepsis, an observation which is phenocopied by the transient inhibition of complex I of ETC by phenformin. Together, we demonstrate that ribosome-targeting antibiotics are beneficial beyond their antibacterial activity and that mitochondrial protein synthesis inhibition leading to ETC perturbation is a novel mechanism for the induction of disease tolerance.

During heart regeneration in the zebrafish, fibrotic tissue is replaced by newly formed cardiomyocytes derived from pre-existing ones. It is unclear whether the heart is comprised of several cardiomyocyte populations bearing different capacity to replace lost myocardium. Here, usingsox10genetic fate mapping, we identified a subset of pre-existent cardiomyocytes in the adult zebrafish heart with a distinct gene expression profile that expanded massively after cryoinjury. Genetic ablation ofsox10+cardiomyocytes severely impaired cardiac regeneration revealing that they play a crucial role for heart regeneration.

SummaryMicrobial communities in nature often feature complex compositional dynamics yet also stable coexistence of diverse species. The mechanistic underpinnings of such dynamic stability remain unclear as system-wide studies have been limited to small engineered communities or synthetic assemblies. Here we show how kefir, a natural milk-fermenting community, realizes stable coexistence through spatio-temporal orchestration of species and metabolite dynamics. During milk fermentation, kefir grains (a polysaccharide matrix synthesized by kefir microbes) grow in mass but remain unchanged in composition. In contrast, the milk is colonized in a dynamic fashion with early members opening metabolic niches for the followers. Through large-scale mapping of metabolic preferences and inter-species interactions, we show how microbes poorly suited for milk survive in, and even dominate, the community through metabolic cooperation and uneven partitioning between the grain and the liquid phase. Overall, our findings reveal how spatio-temporal dynamics promote stable coexistence and have implications for deciphering and modulating complex microbial ecosystems.

Alternative mRNA splicing patterns are determined by the combinatorial control of regulator proteins and their target RNA sequences. We have recently characterized human hnRNP L as a global regulator of alternative splicing, binding to diverse C/A-rich elements. To systematically identify hnRNP L target genes on a genome-wide level, we have combined splice-sensitive microarray analysis and an RNAi-knockdown approach. As a result, we describe 11 target genes of hnRNP L that were validated by RT-PCR and that represent several new modes of hnRNP L-dependent splicing regulation, involving both activator and repressor functions: first, intron retention; second, inclusion or skipping of cassette-type exons; third, suppression of multiple exons; and fourth, alternative poly(A) site selection. In sum, this approach revealed a surprising diversity of splicing-regulatory processes as well as poly(A) site selection in which hnRNP L is involved.