EBNA2-EBF1 complexes promote MYC expression and metabolic processes driving S-phase progression of Epstein-Barr virus-infected B cells.

Epstein-Barr virus (EBV) is a human tumor virus which preferentially infects resting human B cells. Upon infection in vitro, EBV activates and immortalizes these cells. The viral latent protein EBV nuclear antigen 2 (EBNA2) is essential for B cell activation and immortalization; it targets and binds the cellular and ubiquitously expressed DNA-binding protein CBF1, thereby transactivating a plethora of viral and cellular genes. In addition, EBNA2 uses its N-terminal dimerization (END) domain to bind early B cell factor 1 (EBF1), a pioneer transcription factor specifying the B cell lineage. We found that EBNA2 exploits EBF1 to support key metabolic processes and to foster cell cycle progression of infected B cells in their first cell cycles upon activation. The α1-helix within the END domain was found to promote EBF1 binding. EBV mutants lacking the α1-helix in EBNA2 can infect and activate B cells efficiently, but activated cells fail to complete the early S phase of their initial cell cycle. Expression of MYC, target genes of MYC and E2F, as well as multiple metabolic processes linked to cell cycle progression are impaired in EBVΔα1-infected B cells. Our findings indicate that EBF1 controls B cell activation via EBNA2 and, thus, has a critical role in regulating the cell cycle of EBV-infected B cells. This is a function of EBF1 going beyond its well-known contribution to B cell lineage specification.

Tissue-specific regulation of lipid polyester synthesis genes controlling oxygen permeation into Lotus japonicus nodules.

Legumes establish endosymbiotic associations with nitrogen-fixing rhizobia, which they host inside root nodules. Here, specific physiological and morphological adaptations, such as the production of oxygen-binding leghemoglobin proteins and the formation of an oxygen diffusion barrier in the nodule periphery, are essential to protect the oxygen-labile bacterial nitrogenase enzyme. The molecular basis of the latter process remains elusive as the identification of required genes is limited by the epistatic effect of nodule organogenesis over nodule infection and rhizobia accommodation. We overcame this by exploring the phenotypic diversity of Lotus japonicus accessions that uncouple nodule organogenesis from nodule infection when inoculated with a subcompatible Rhizobium strain. Using comparative transcriptomics, we identified genes with functions associated with oxygen homeostasis and deposition of lipid polyesters on cell walls to be specifically up-regulated in infected compared to noninfected nodules. As hydrophobic modification of cell walls is pivotal for creating diffusion barriers like the root endodermis, we focused on two Fatty acyl-CoA Reductase genes that were specifically activated in the root and/or in the nodule endodermis. Mutant lines in a Fatty acyl-CoA Reductase gene expressed exclusively in the nodule endodermis had decreased deposition of polyesters on this cell layer and increased nodule permeability compared to wild-type plants. Oxygen concentrations were significantly increased in the inner cortex of mutant nodules, which correlated with reduced nitrogenase activity, and impaired shoot growth. These results provide the first genetic evidence for the formation of the nodule oxygen diffusion barrier, a key adaptation enabling nitrogen fixation in legume nodules.

TET1 regulates gene expression and repression of endogenous retroviruses independent of DNA demethylation.

DNA methylation (5-methylcytosine (5mC)) is critical for genome stability and transcriptional regulation in mammals. The discovery that ten-eleven translocation (TET) proteins catalyze the oxidation of 5mC to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC) revolutionized our perspective on the complexity and regulation of DNA modifications. However, to what extent the regulatory functions of TET1 can be attributed to its catalytic activity remains unclear. Here, we use genome engineering and quantitative multi-omics approaches to dissect the precise catalytic vs. non-catalytic functions of TET1 in murine embryonic stem cells (mESCs). Our study identifies TET1 as an essential interaction hub for multiple chromatin modifying complexes and a global regulator of histone modifications. Strikingly, we find that the majority of transcriptional regulation depends on non-catalytic functions of TET1. In particular, we show that TET1 is critical for the establishment of H3K9me3 and H4K20me3 at endogenous retroviral elements (ERVs) and their silencing that is independent of its canonical role in DNA demethylation. Furthermore, we provide evidence that this repression of ERVs depends on the interaction between TET1 and SIN3A. In summary, we demonstrate that the non-catalytic functions of TET1 are critical for regulation of gene expression and the silencing of endogenous retroviruses in mESCs.

Regulatory and coding sequences of TRNP1 co-evolve with brain size and cortical folding in mammals.

Brain size and cortical folding have increased and decreased recurrently during mammalian evolution. Identifying genetic elements whose sequence or functional properties co-evolve with these traits can provide unique information on evolutionary and developmental mechanisms. A good candidate for such a comparative approach is TRNP1, as it controls proliferation of neural progenitors in mice and ferrets. Here, we investigate the contribution of both regulatory and coding sequences of TRNP1 to brain size and cortical folding in over 30 mammals. We find that the rate of TRNP1 protein evolution (ω) significantly correlates with brain size, slightly less with cortical folding and much less with body size. This brain correlation is stronger than for >95% of random control proteins. This co-evolution is likely affecting TRNP1 activity, as we find that TRNP1 from species with larger brains and more cortical folding induce higher proliferation rates in neural stem cells. Furthermore, we compare the activity of putative cis-regulatory elements (CREs) of TRNP1 in a massively parallel reporter assay and identify one CRE that likely co-evolves with cortical folding in Old World monkeys and apes. Our analyses indicate that coding and regulatory changes that increased TRNP1 activity were positively selected either as a cause or a consequence of increases in brain size and cortical folding. They also provide an example how phylogenetic approaches can inform biological mechanisms, especially when combined with molecular phenotypes across several species.

Prime-seq, efficient and powerful bulk RNA sequencing.

Cost-efficient library generation by early barcoding has been central in propelling single-cell RNA sequencing. Here, we optimize and validate prime-seq, an early barcoding bulk RNA-seq method. We show that it performs equivalently to TruSeq, a standard bulk RNA-seq method, but is fourfold more cost-efficient due to almost 50-fold cheaper library costs. We also validate a direct RNA isolation step, show that intronic reads are derived from RNA, and compare cost-efficiencies of available protocols. We conclude that prime-seq is currently one of the best options to set up an early barcoding bulk RNA-seq protocol from which many labs would profit.

In vivo PDX CRISPR/Cas9 screens reveal mutual therapeutic targets to overcome heterogeneous acquired chemo-resistance.

Resistance towards cancer treatment represents a major clinical obstacle, preventing cure of cancer patients. To gain mechanistic insights, we developed a model for acquired resistance to chemotherapy by treating mice carrying patient derived xenografts (PDX) of acute lymphoblastic leukemia with widely-used cytotoxic drugs for 18 consecutive weeks. In two distinct PDX samples, tumors initially responded to treatment, until stable disease and eventually tumor re-growth evolved under therapy, at highly similar kinetics between replicate mice. Notably, replicate tumors developed different mutations in TP53 and individual sets of chromosomal alterations, suggesting independent parallel clonal evolution rather than selection, driven by a combination of stochastic and deterministic processes. Transcriptome and proteome showed shared dysregulations between replicate tumors providing putative targets to overcome resistance. In vivo CRISPR/Cas9 dropout screens in PDX revealed broad dependency on BCL2, BRIP1 and COPS2. Accordingly, venetoclax re-sensitized derivative tumors towards chemotherapy, despite genomic heterogeneity, demonstrating direct translatability of the approach. Hence, despite the presence of multiple resistance-associated genomic alterations, effective rescue treatment for polychemotherapy-resistant tumors can be identified using functional testing in preclinical models.

PARP14 is a novel target in STAT6 mutant follicular lymphoma.

The variable clinical course of follicular lymphoma (FL) is determined by the molecular heterogeneity of tumor cells and complex interactions within the tumor microenvironment (TME). IL-4 producing follicular helper T cells (TFH) are critical components of the FL TME. Binding of IL-4 to IL-4R on FL cells activates JAK/STAT signaling. We identified STAT6 mutations (STAT6MUT) in 13% of FL (N = 33/258), all clustered within the DNA binding domain. Gene expression data and immunohistochemistry showed upregulation of IL-4/STAT6 target genes in STAT6MUT FL, including CCL17, CCL22, and FCER2 (CD23). Functionally, STAT6MUT was gain-of-function by serial replating phenotype in pre-B CFU assays. Expression of STAT6MUT enhanced IL-4 induced FCER2/CD23, CCL17 and CCL22 expression and was associated with nuclear accumulation of pSTAT6. RNA sequencing identified PARP14 -a transcriptional switch and co-activator of STAT6- among the top differentially upregulated genes in IL-4 stimulated STAT6MUT lymphoma cells and in STAT6MUT primary FL cells. Quantitative chromatin immunoprecipitation (qChIP) demonstrated binding of STAT6MUT but not STAT6WT to the PARP14 promotor. Reporter assays showed increased IL-4 induced transactivation activity of STAT6MUT at the PARP14 promotor, suggesting a self-reinforcing regulatory circuit. Knock-down of PARP14 or PARP-inhibition abrogated the STAT6MUT gain-of-function phenotype. Thus, our results identify PARP14 as a novel therapeutic target in STAT6MUT FL.

Protective immune trajectories in early viral containment of non-pneumonic SARS-CoV-2 infection.

The antiviral immune response to SARS-CoV-2 infection can limit viral spread and prevent development of pneumonic COVID-19. However, the protective immunological response associated with successful viral containment in the upper airways remains unclear. Here, we combine a multi-omics approach with longitudinal sampling to reveal temporally resolved protective immune signatures in non-pneumonic and ambulatory SARS-CoV-2 infected patients and associate specific immune trajectories with upper airway viral containment. We see a distinct systemic rather than local immune state associated with viral containment, characterized by interferon stimulated gene (ISG) upregulation across circulating immune cell subsets in non-pneumonic SARS-CoV2 infection. We report reduced cytotoxic potential of Natural Killer (NK) and T cells, and an immune-modulatory monocyte phenotype associated with protective immunity in COVID-19. Together, we show protective immune trajectories in SARS-CoV2 infection, which have important implications for patient prognosis and the development of immunomodulatory therapies.

Retrotransposons as pathogenicity factors of the plant pathogenic fungus Botrytis cinerea.

BACKGROUND: Retrotransposons are genetic elements inducing mutations in all domains of life. Despite their detrimental effect, retrotransposons can become temporarily active during epigenetic reprogramming and cellular stress response, which may accelerate host genome evolution. In fungal pathogens, a positive role has been attributed to retrotransposons when shaping genome architecture and expression of genes encoding pathogenicity factors; thus, retrotransposons are known to influence pathogenicity. RESULTS: We uncover a hitherto unknown role of fungal retrotransposons as being pathogenicity factors, themselves. The aggressive fungal plant pathogen, Botrytis cinerea, is known to deliver some long-terminal repeat (LTR) deriving regulatory trans-species small RNAs (BcsRNAs) into plant cells to suppress host gene expression for infection. We find that naturally occurring, less aggressive B. cinerea strains possess considerably lower copy numbers of LTR retrotransposons and had lost retrotransposon BcsRNA production. Using a transgenic proof-of-concept approach, we reconstitute retrotransposon expression in a BcsRNA-lacking B. cinerea strain, which results in enhanced aggressiveness in a retrotransposon and BcsRNA expression-dependent manner. Moreover, retrotransposon expression in B. cinerea leads to suppression of plant defence-related genes during infection. CONCLUSIONS: We propose that retrotransposons are pathogenicity factors that manipulate host plant gene expression by encoding trans-species BcsRNAs. Taken together, the novelty that retrotransposons are pathogenicity factors will have a broad impact on studies of host-microbe interactions and pathology.

A non-invasive method to generate induced pluripotent stem cells from primate urine.

Comparing the molecular and cellular properties among primates is crucial to better understand human evolution and biology. However, it is difficult or ethically impossible to collect matched tissues from many primates, especially during development. An alternative is to model different cell types and their development using induced pluripotent stem cells (iPSCs). These can be generated from many tissue sources, but non-invasive sampling would decisively broaden the spectrum of non-human primates that can be investigated. Here, we report the generation of primate iPSCs from urine samples. We first validate and optimize the procedure using human urine samples and show that suspension- Sendai Virus transduction of reprogramming factors into urinary cells efficiently generates integration-free iPSCs, which maintain their pluripotency under feeder-free culture conditions. We demonstrate that this method is also applicable to gorilla and orangutan urinary cells isolated from a non-sterile zoo floor. We characterize the urinary cells, iPSCs and derived neural progenitor cells using karyotyping, immunohistochemistry, differentiation assays and RNA-sequencing. We show that the urine-derived human iPSCs are indistinguishable from well characterized PBMC-derived human iPSCs and that the gorilla and orangutan iPSCs are well comparable to the human iPSCs. In summary, this study introduces a novel and efficient approach to non-invasively generate iPSCs from primate urine. This will extend the zoo of species available for a comparative approach to molecular and cellular phenotypes.

Loss-of-function mutations in the histone methyltransferase EZH2 promote chemotherapy resistance in AML.

Chemotherapy resistance is the main impediment in the treatment of acute myeloid leukaemia (AML). Despite rapid advances, the various mechanisms inducing resistance development remain to be defined in detail. Here we report that loss-of-function mutations (LOF) in the histone methyltransferase EZH2 have the potential to confer resistance against the chemotherapeutic agent cytarabine. We identify seven distinct EZH2 mutations leading to loss of H3K27 trimethylation via multiple mechanisms. Analysis of matched diagnosis and relapse samples reveal a heterogenous regulation of EZH2 and a loss of EZH2 in 50% of patients. We confirm that loss of EZH2 induces resistance against cytarabine in the cell lines HEK293T and K562 as well as in a patient-derived xenograft model. Proteomics and transcriptomics analysis reveal that resistance is conferred by upregulation of multiple direct and indirect EZH2 target genes that are involved in apoptosis evasion, augmentation of proliferation and alteration of transmembrane transporter function. Our data indicate that loss of EZH2 results in upregulation of its target genes, providing the cell with a selective growth advantage, which mediates chemotherapy resistance.

Benchmarking single-cell RNA-sequencing protocols for cell atlas projects.

Single-cell RNA sequencing (scRNA-seq) is the leading technique for characterizing the transcriptomes of individual cells in a sample. The latest protocols are scalable to thousands of cells and are being used to compile cell atlases of tissues, organs and organisms. However, the protocols differ substantially with respect to their RNA capture efficiency, bias, scale and costs, and their relative advantages for different applications are unclear. In the present study, we generated benchmark datasets to systematically evaluate protocols in terms of their power to comprehensively describe cell types and states. We performed a multicenter study comparing 13 commonly used scRNA-seq and single-nucleus RNA-seq protocols applied to a heterogeneous reference sample resource. Comparative analysis revealed marked differences in protocol performance. The protocols differed in library complexity and their ability to detect cell-type markers, impacting their predictive value and suitability for integration into reference cell atlases. These results provide guidance both for individual researchers and for consortium projects such as the Human Cell Atlas.

Studying Cancer Heterogeneity by Single-Cell RNA Sequencing.

A major hurdle for the treatment of cancer is the incomplete understanding of its evolution through the course of its emergence, dispersal, and relapse. Genetic and epigenetic changes in combination with external cues and selective forces are the driving factors behind tumor heterogeneity. Understanding this variability within and across patients may partly explain the unpredictable outcomes of cancer treatments. Measuring the variation of gene expression levels within cells of the same tumor is a crucial part of this endeavor. Hence, the recently developed single-cell RNA-sequencing (scRNA-seq) technologies have become a valuable tool for cancer research. In practice, however, this is still challenging, especially for clinical samples. Here, we describe mcSCRB-seq (molecular crowding single-cell RNA barcoding and sequencing), a highly sensitive and powerful plate-based scRNA-seq method, which shows great capability to generate transcriptome data for cancer cells. mcSCRB-seq is not only characterized by high sensitivity due to molecular crowding and the use of unique molecular identifiers (UMIs) but also features an easy workflow and a low per-cell cost and does not require specialized equipment.

Sensitive and powerful single-cell RNA sequencing using mcSCRB-seq.

Single-cell RNA sequencing (scRNA-seq) has emerged as a central genome-wide method to characterize cellular identities and processes. Consequently, improving its sensitivity, flexibility, and cost-efficiency can advance many research questions. Among the flexible plate-based methods, single-cell RNA barcoding and sequencing (SCRB-seq) is highly sensitive and efficient. Here, we systematically evaluate experimental conditions of this protocol and find that adding polyethylene glycol considerably increases sensitivity by enhancing cDNA synthesis. Furthermore, using Terra polymerase increases efficiency due to a more even cDNA amplification that requires less sequencing of libraries. We combined these and other improvements to develop a scRNA-seq library protocol we call molecular crowding SCRB-seq (mcSCRB-seq), which we show to be one of the most sensitive, efficient, and flexible scRNA-seq methods to date.