Inheritable Silencing of Endogenous Genes by Hit-and-Run Targeted Epigenetic Editing.

Gene silencing is instrumental to interrogate gene function and holds promise for therapeutic applications. Here, we repurpose the endogenous retroviruses' silencing machinery of embryonic stem cells to stably silence three highly expressed genes in somatic cells by epigenetics. This was achieved by transiently expressing combinations of engineered transcriptional repressors that bind to and synergize at the target locus to instruct repressive histone marks and de novo DNA methylation, thus ensuring long-term memory of the repressive epigenetic state. Silencing was highly specific, as shown by genome-wide analyses, sharply confined to the targeted locus without spreading to nearby genes, resistant to activation induced by cytokine stimulation, and relieved only by targeted DNA demethylation. We demonstrate the portability of this technology by multiplex gene silencing, adopting different DNA binding platforms and interrogating thousands of genomic loci in different cell types, including primary T lymphocytes. Targeted epigenome editing might have broad application in research and medicine.

Interferon-inducible phospholipids govern IFITM3-dependent endosomal antiviral immunity.

The interferon-induced transmembrane proteins (IFITM) are implicated in several biological processes, including antiviral defense, but their modes of action remain debated. Here, taking advantage of pseudotyped viral entry assays and replicating viruses, we uncover the requirement of host co-factors for endosomal antiviral inhibition through high-throughput proteomics and lipidomics in cellular models of IFITM restriction. Unlike plasma membrane (PM)-localized IFITM restriction that targets infectious SARS-CoV2 and other PM-fusing viral envelopes, inhibition of endosomal viral entry depends on lysines within the conserved IFITM intracellular loop. These residues recruit Phosphatidylinositol 3,4,5-trisphosphate (PIP3) that we show here to be required for endosomal IFITM activity. We identify PIP3 as an interferon-inducible phospholipid that acts as a rheostat for endosomal antiviral immunity. PIP3 levels correlated with the potency of endosomal IFITM restriction and exogenous PIP3 enhanced inhibition of endocytic viruses, including the recent SARS-CoV2 Omicron variant. Together, our results identify PIP3 as a critical regulator of endosomal IFITM restriction linking it to the Pi3K/Akt/mTORC pathway and elucidate cell-compartment-specific antiviral mechanisms with potential relevance for the development of broadly acting antiviral strategies.

Loss of ribonuclease DIS3 hampers genome integrity in myeloma by disrupting DNA:RNA hybrid metabolism.

The ribonuclease DIS3 is one of the most frequently mutated genes in the hematological cancer multiple myeloma, yet the basis of its tumor suppressor function in this disease remains unclear. Herein, exploiting the TCGA dataset, we found that DIS3 plays a prominent role in the DNA damage response. DIS3 inactivation causes genomic instability by increasing mutational load, and a pervasive accumulation of DNA:RNA hybrids that induces genomic DNA double-strand breaks (DSBs). DNA:RNA hybrid accumulation also prevents binding of the homologous recombination (HR) machinery to double-strand breaks, hampering DSB repair. DIS3-inactivated cells become sensitive to PARP inhibitors, suggestive of a defect in homologous recombination repair. Accordingly, multiple myeloma patient cells mutated for DIS3 harbor an increased mutational burden and a pervasive overexpression of pro-inflammatory interferon, correlating with the accumulation of DNA:RNA hybrids. We propose DIS3 loss in myeloma to be a driving force for tumorigenesis via DNA:RNA hybrid-dependent enhanced genome instability and increased mutational rate. At the same time, DIS3 loss represents a liability that might be therapeutically exploited in patients whose cancer cells harbor DIS3 mutations.

Scalable integration of multiomic single-cell data using generative adversarial networks.

MOTIVATION: Single-cell profiling has become a common practice to investigate the complexity of tissues, organs, and organisms. Recent technological advances are expanding our capabilities to profile various molecular layers beyond the transcriptome such as, but not limited to, the genome, the epigenome, and the proteome. Depending on the experimental procedure, these data can be obtained from separate assays or the very same cells. Yet, integration of more than two assays is currently not supported by the majority of the computational frameworks avaiable. RESULTS: We here propose a Multi-Omic data integration framework based on Wasserstein Generative Adversarial Networks suitable for the analysis of paired or unpaired data with a high number of modalities (>2). At the core of our strategy is a single network trained on all modalities together, limiting the computational burden when many molecular layers are evaluated. AVAILABILITY AND IMPLEMENTATION: Source code of our framework is available at https://github.com/vgiansanti/MOWGAN.

Tumor suppressor PNRC1 blocks rRNA maturation by recruiting the decapping complex to the nucleolus.

Focal deletions occur frequently in the cancer genome. However, the putative tumor-suppressive genes residing within these regions have been difficult to pinpoint. To robustly identify these genes, we implemented a computational approach based on non-negative matrix factorization, NMF, and interrogated the TCGA dataset. This analysis revealed a metagene signature including a small subset of genes showing pervasive hemizygous deletions, reduced expression in cancer patient samples, and nucleolar function. Amid the genes belonging to this signature, we have identified PNRC1, a nuclear receptor coactivator. We found that PNRC1 interacts with the cytoplasmic DCP1α/DCP2 decapping machinery and hauls it inside the nucleolus. PNRC1-dependent nucleolar translocation of the decapping complex is associated with a decrease in the 5'-capped U3 and U8 snoRNA fractions, hampering ribosomal RNA maturation. As a result, PNRC1 ablates the enhanced proliferation triggered by established oncogenes such as RAS and MYC These observations uncover a previously undescribed mechanism of tumor suppression, whereby the cytoplasmic decapping machinery is hauled within nucleoli, tightly regulating ribosomal RNA maturation.

Nested Stochastic Block Models applied to the analysis of single cell data.

Single cell profiling has been proven to be a powerful tool in molecular biology to understand the complex behaviours of heterogeneous system. The definition of the properties of single cells is the primary endpoint of such analysis, cells are typically clustered to underpin the common determinants that can be used to describe functional properties of the cell mixture under investigation. Several approaches have been proposed to identify cell clusters; while this is matter of active research, one popular approach is based on community detection in neighbourhood graphs by optimisation of modularity. In this paper we propose an alternative and principled solution to this problem, based on Stochastic Block Models. We show that such approach not only is suitable for identification of cell groups, it also provides a solid framework to perform other relevant tasks in single cell analysis, such as label transfer. To encourage the use of Stochastic Block Models, we developed a python library, schist, that is compatible with the popular scanpy framework.

Gene Expression Analysis in Patients with Cocaine-Induced Midline Destructive Lesions.

BACKGROUND AND OBJECTIVES: Cocaine users may present with positive antineutrophil cytoplasmic antibodies (ANCA) and severe midline destructive lesions (CIMDL) which are histologically characterized by massive apoptosis. However, histopathological and laboratory studies suggest that autoimmunity may not be the main pathogenic driver. We analyzed gene expression both in cell lines of nasal mucosa exposed to cocaine and in CIMDL patients to determine whether genetic predisposition might cause such lesions, which are observed in a minority of cocaine abusers. MATERIALS AND METHODS: The genetic expression profile of nasal mucosa exposed to cocaine was analyzed. Rare variants of expressed genes were searched in patients with CIMDL using exome sequencing and bio-informatics. RESULTS: We identified 462 genes that were induced by cocaine, mainly related to apoptosis and autophagy in response to oxidative stress. Under the hypothesis that genes linked to the phenotype are also induced by cocaine itself, a rare variants burden test was performed to select genes that were significantly enriched in rare mutations. Next, 11 cocaine abusers with CIMDL and no other relevant medical comorbidities underwent exome sequencing, and 12 genes that were significantly enriched in the burden test and present in at least 10 patients were identified. An in-depth analysis of these genes revealed their involvement in apoptosis, tissue homeostasis, autophagy, and response to oxidative stress. CONCLUSIONS: Oxidative stress and rare genetic alterations in the response to reactive oxygen species, apoptosis, autophagy, and tissue regeneration are plausible drivers of damage affecting nasal mucosa exposed to cocaine crystals and, consequently, the pathogenic mechanism behind CIMDL.

Identification of novel genetic variants associated with short stature in a Baka Pygmies population.

Human growth is a complex trait determined by genetic factors in combination with external stimuli, including environment, nutrition and hormonal status. In the past, several genome-wide association studies (GWAS) have collectively identified hundreds of genetic variants having a putative effect on determining adult height in different worldwide populations. Theoretically, a valuable approach to better understand the mechanisms of complex traits as adult height is to study a population exhibiting extreme stature phenotypes, such as African Baka Pygmies. After phenotypic characterization, we sequenced the whole exomes of a cohort of Baka Pygmies and their non-Pygmies Bantu neighbors to highlight genetic variants associated with the reduced stature. Whole exome data analysis revealed 29 single nucleotide polymorphisms (SNPs) significantly associated with the reduced height in the Baka group. Among these variants, we focused on SNP rs7629425, located in the 5'-UTR of the Hyaluronidase-2 (HYAL2) gene. The frequency of the alternative allele was significantly increased compared to African and non-African populations. In vitro luciferase assay showed significant differences in transcription modulation by rs7629425 C/T alleles. In conclusion, our results suggested that the HYAL2 gene variants may play a role in the etiology of short stature in Baka Pygmies population.

TBC1D24-TLDc-related epilepsy exercise-induced dystonia: rescue by antioxidants in a disease model.

Genetic mutations in TBC1D24 have been associated with multiple phenotypes, with epilepsy being the main clinical manifestation. The TBC1D24 protein consists of the unique association of a Tre2/Bub2/Cdc16 (TBC) domain and a TBC/lysin motif domain/catalytic (TLDc) domain. More than 50 missense and loss-of-function mutations have been described and are spread over the entire protein. Through whole genome/exome sequencing we identified compound heterozygous mutations, R360H and G501R, within the TLDc domain, in an index family with a Rolandic epilepsy exercise-induced dystonia phenotype (http://omim.org/entry/608105). A 20-year long clinical follow-up revealed that epilepsy was self-limited in all three affected patients, but exercise-induced dystonia persisted into adulthood in two. Furthermore, we identified three additional sporadic paediatric patients with a remarkably similar phenotype, two of whom had compound heterozygous mutations consisting of an in-frame deletion I81_K84 and an A500V mutation, and the third carried T182M and G511R missense mutations, overall revealing that all six patients harbour a missense mutation in the subdomain of TLDc between residues 500 and 511. We solved the crystal structure of the conserved Drosophila TLDc domain. This allowed us to predict destabilizing effects of the G501R and G511R mutations and, to a lesser degree, of R360H and potentially A500V. Next, we characterized the functional consequences of a strong and a weak TLDc mutation (TBC1D24G501R and TBC1D24R360H) using Drosophila, where TBC1D24/Skywalker regulates synaptic vesicle trafficking. In a Drosophila model neuronally expressing human TBC1D24, we demonstrated that the TBC1D24G501R TLDc mutation causes activity-induced locomotion and synaptic vesicle trafficking defects, while TBC1D24R360H is benign. The neuronal phenotypes of the TBC1D24G501R mutation are consistent with exacerbated oxidative stress sensitivity, which is rescued by treating TBC1D24G501R mutant animals with antioxidants N-acetylcysteine amide or α-tocopherol as indicated by restored synaptic vesicle trafficking levels and sustained behavioural activity. Our data thus show that mutations in the TLDc domain of TBC1D24 cause Rolandic-type focal motor epilepsy and exercise-induced dystonia. The humanized TBC1D24G501R fly model exhibits sustained activity and vesicle transport defects. We propose that the TBC1D24/Sky TLDc domain is a reactive oxygen species sensor mediating synaptic vesicle trafficking rates that, when dysfunctional, causes a movement disorder in patients and flies. The TLDc and TBC domain mutations' response to antioxidant treatment we observed in the animal model suggests a potential for combining antioxidant-based therapeutic approaches to TBC1D24-associated disorders with previously described lipid-altering strategies for TBC domain mutations.

miR-17∼92 family clusters control iNKT cell ontogenesis via modulation of TGF-ß signaling.

Invariant natural killer T cells (iNKT) cells are T lymphocytes displaying innate effector functions, acquired through a distinct thymic developmental program regulated by microRNAs (miRNAs). Deleting miRNAs by Dicer ablation (Dicer KO) in thymocytes selectively impairs iNKT cell survival and functional differentiation. To unravel this miRNA-dependent program, we systemically identified transcripts that were differentially expressed between WT and Dicer KO iNKT cells at different differentiation stages and predicted to be targeted by the iNKT cell-specific miRNAs. TGF-ß receptor II (TGF-ßRII), critically implicated in iNKT cell differentiation, was found up-regulated in iNKT Dicer KO cells together with enhanced TGF-ß signaling. miRNA members of the miR-17∼92 family clusters were predicted to target Tgfbr2 mRNA upon iNKT cell development. iNKT cells lacking all three miR-17∼92 family clusters (miR-17∼92, miR-106a∼363, miR-106b∼25) phenocopied both increased TGF-ßRII expression and signaling, and defective effector differentiation, displayed by iNKT Dicer KO cells. Consistently, genetic ablation of TGF-ß signaling in the absence of miRNAs rescued iNKT cell differentiation. These results elucidate the global impact of miRNAs on the iNKT cell developmental program and uncover the targeting of a lineage-specific cytokine signaling by miRNAs as a mechanism regulating innate-like T-cell development and effector differentiation.