Tubular cell and keratinocyte single-cell transcriptomics applied to lupus nephritis reveal type I IFN and fibrosis relevant pathways.

The molecular and cellular processes that lead to renal damage and to the heterogeneity of lupus nephritis (LN) are not well understood. We applied single-cell RNA sequencing (scRNA-seq) to renal biopsies from patients with LN and evaluated skin biopsies as a potential source of diagnostic and prognostic markers of renal disease. Type I interferon (IFN)-response signatures in tubular cells and keratinocytes distinguished patients with LN from healthy control subjects. Moreover, a high IFN-response signature and fibrotic signature in tubular cells were each associated with failure to respond to treatment. Analysis of tubular cells from patients with proliferative, membranous and mixed LN indicated pathways relevant to inflammation and fibrosis, which offer insight into their histologic differences. In summary, we applied scRNA-seq to LN to deconstruct its heterogeneity and identify novel targets for personalized approaches to therapy.

Author Correction: Tubular cell and keratinocyte single-cell transcriptomics applied to lupus nephritis reveal type I IFN and fibrosis relevant pathways.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Human CLP1 mutations alter tRNA biogenesis, affecting both peripheral and central nervous system function.

CLP1 is a RNA kinase involved in tRNA splicing. Recently, CLP1 kinase-dead mice were shown to display a neuromuscular disorder with loss of motor neurons and muscle paralysis. Human genome analyses now identified a CLP1 homozygous missense mutation (p.R140H) in five unrelated families, leading to a loss of CLP1 interaction with the tRNA splicing endonuclease (TSEN) complex, largely reduced pre-tRNA cleavage activity, and accumulation of linear tRNA introns. The affected individuals develop severe motor-sensory defects, cortical dysgenesis, and microcephaly. Mice carrying kinase-dead CLP1 also displayed microcephaly and reduced cortical brain volume due to the enhanced cell death of neuronal progenitors that is associated with reduced numbers of cortical neurons. Our data elucidate a neurological syndrome defined by CLP1 mutations that impair tRNA splicing. Reduction of a founder mutation to homozygosity illustrates the importance of rare variations in disease and supports the clan genomics hypothesis.

Convergence of mammalian RQC and C-end rule proteolytic pathways via alanine tailing.

Incompletely synthesized nascent chains obstructing large ribosomal subunits are targeted for degradation by ribosome-associated quality control (RQC). In bacterial RQC, RqcH marks the nascent chains with C-terminal alanine (Ala) tails that are directly recognized by proteasome-like proteases, whereas in eukaryotes, RqcH orthologs (Rqc2/NEMF [nuclear export mediator factor]) assist the Ltn1/Listerin E3 ligase in nascent chain ubiquitylation. Here, we study RQC-mediated proteolytic targeting of ribosome stalling products in mammalian cells. We show that mammalian NEMF has an additional, Listerin-independent proteolytic role, which, as in bacteria, is mediated by tRNA-Ala binding and Ala tailing. However, in mammalian cells Ala tails signal proteolysis indirectly, through a pathway that recognizes C-terminal degrons; we identify the CRL2KLHDC10 E3 ligase complex and the novel C-end rule E3, Pirh2/Rchy1, as bona fide RQC pathway components that directly bind to Ala-tailed ribosome stalling products and target them for degradation. As Listerin mutation causes neurodegeneration in mice, functionally redundant E3s may likewise be implicated in molecular mechanisms of neurodegeneration.

Cyclic [G(2',5')pA(3',5')p] is the metazoan second messenger produced by DNA-activated cyclic GMP-AMP synthase.

Recent studies identified cyclic GMP-AMP (cGAMP) as a metazoan second messenger triggering an interferon response. cGAMP is generated from GTP and ATP by cytoplasmic dsDNA sensor cGAMP synthase (cGAS). We combined structural, chemical, biochemical, and cellular assays to demonstrate that this second messenger contains G(2',5')pA and A(3',5')pG phosphodiester linkages, designated c[G(2',5')pA(3',5')p]. We show that, upon dsDNA binding, cGAS is activated through conformational transitions, resulting in formation of a catalytically competent and accessible nucleotide-binding pocket for generation of c[G(2',5')pA(3',5')p]. We demonstrate that cyclization occurs in a stepwise manner through initial generation of 5'-pppG(2',5')pA prior to cyclization to c[G(2',5')pA(3',5')p], with the latter positioned precisely in the catalytic pocket. Mutants of cGAS dsDNA-binding or catalytic pocket residues exhibit reduced or abrogated activity. Our studies have identified c[G(2',5')pA(3',5')p] as a founding member of a family of metazoan 2',5'-containing cyclic heterodinucleotide second messengers distinct from bacterial 3',5' cyclic dinucleotides.

Structure-function analysis of STING activation by c[G(2',5')pA(3',5')p] and targeting by antiviral DMXAA.

Binding of dsDNA by cyclic GMP-AMP (cGAMP) synthase (cGAS) triggers formation of the metazoan second messenger c[G(2',5')pA(3',5')p], which binds the signaling protein STING with subsequent activation of the interferon (IFN) pathway. We show that human hSTING(H232) adopts a "closed" conformation upon binding c[G(2',5')pA(3',5')p] and its linkage isomer c[G(2',5')pA(2',5')p], as does mouse mSting(R231) on binding c[G(2',5')pA(3',5')p], c[G(3',5')pA(3',5')p] and the antiviral agent DMXAA, leading to similar "closed" conformations. Comparing hSTING to mSting, 2',5'-linkage-containing cGAMP isomers were more specific triggers of the IFN pathway compared to the all-3',5'-linkage isomer. Guided by structural information, we identified a unique point mutation (S162A) placed within the cyclic-dinucleotide-binding site of hSTING that rendered it sensitive to the otherwise mouse-specific drug DMXAA, a conclusion validated by binding studies. Our structural and functional analysis highlights the unexpected versatility of STING in the recognition of natural and synthetic ligands within a small-molecule pocket created by the dimerization of STING.

The G3BP1-Family-USP10 Deubiquitinase Complex Rescues Ubiquitinated 40S Subunits of Ribosomes Stalled in Translation from Lysosomal Degradation.

Ribosome-associated quality control (RQC) purges aberrant mRNAs and nascent polypeptides in a multi-step molecular process initiated by the E3 ligase ZNF598 through sensing of ribosomes collided at aberrant mRNAs and monoubiquitination of distinct small ribosomal subunit proteins. We show that G3BP1-family-USP10 complexes are required for deubiquitination of RPS2, RPS3, and RPS10 to rescue modified 40S subunits from programmed degradation. Knockout of USP10 or G3BP1 family proteins increased lysosomal ribosomal degradation and perturbed ribosomal subunit stoichiometry, both of which were rescued by a single K214R substitution of RPS3. While the majority of RPS2 and RPS3 monoubiquitination resulted from ZNF598-dependent sensing of ribosome collisions initiating RQC, another minor pathway contributed to their monoubiquitination. G3BP1 family proteins have long been considered RNA-binding proteins, however, our results identified 40S subunits and associated mRNAs as their predominant targets, a feature shared by stress granules to which G3BP1 family proteins localize under stress.

A role for neuronal piRNAs in the epigenetic control of memory-related synaptic plasticity.

Small RNA-mediated gene regulation during development causes long-lasting changes in cellular phenotypes. To determine whether small RNAs of the adult brain can regulate memory storage, a process that requires stable and long-lasting changes in the functional state of neurons, we generated small RNA libraries from the Aplysia CNS. In these libraries, we discovered an unexpectedly abundant expression of a 28 nucleotide sized class of piRNAs in brain, which had been thought to be germline specific. These piRNAs have unique biogenesis patterns, predominant nuclear localization, and robust sensitivity to serotonin, a modulatory transmitter that is important for memory. We find that the Piwi/piRNA complex facilitates serotonin-dependent methylation of a conserved CpG island in the promoter of CREB2, the major inhibitory constraint of memory in Aplysia, leading to enhanced long-term synaptic facilitation. These findings provide a small RNA-mediated gene regulatory mechanism for establishing stable long-term changes in neurons for the persistence of memory.

The TIA1 RNA-Binding Protein Family Regulates EIF2AK2-Mediated Stress Response and Cell Cycle Progression.

TIA1 and TIAL1 encode a family of U-rich element mRNA-binding proteins ubiquitously expressed and conserved in metazoans. Using PAR-CLIP, we determined that both proteins bind target sites with identical specificity in 3' UTRs and introns proximal to 5' as well as 3' splice sites. Double knockout (DKO) of TIA1 and TIAL1 increased target mRNA abundance proportional to the number of binding sites and also caused accumulation of aberrantly spliced mRNAs, most of which are subject to nonsense-mediated decay. Loss of PRKRA by mis-splicing triggered the activation of the double-stranded RNA (dsRNA)-activated protein kinase EIF2AK2/PKR and stress granule formation. Ectopic expression of PRKRA cDNA or knockout of EIF2AK2 in DKO cells rescued this phenotype. Perturbation of maturation and/or stability of additional targets further compromised cell cycle progression. Our study reveals the essential contributions of the TIA1 protein family to the fidelity of mRNA maturation, translation, and RNA-stress-sensing pathways in human cells.

A comprehensive analysis of the effects of the deaminase AID on the transcriptome and methylome of activated B cells.

Beyond its well-characterized functions in antibody diversification, the cytidine deaminase AID can catalyze off-target DNA damage and has been hypothesized to edit RNA and mediate DNA demethylation. To comprehensively examine the effects of AID on the transcriptome and the pattern of DNA methylation ('methylome'), we analyzed AID-deficient (Aicda(-/-)), wild-type and AID-overexpressing activated B cells by high-throughput RNA sequencing (RNA-Seq) and reduced-representation bisulfite sequencing (RRBS). These analyses confirmed the known role of AID in immunoglobulin isotype switching and also demonstrated few other effects of AID on gene expression. Additionally, we detected no evidence of AID-dependent editing of mRNA or microRNA. Finally, the RRBS data did not support the proposed role for AID in regulating DNA methylation. Thus, despite evidence of its additional activities in other systems, antibody diversification seems to be the sole physiological function of AID in activated B cells.

Transcriptome-wide identification of RNA-binding protein and microRNA target sites by PAR-CLIP.

RNA transcripts are subject to posttranscriptional gene regulation involving hundreds of RNA-binding proteins (RBPs) and microRNA-containing ribonucleoprotein complexes (miRNPs) expressed in a cell-type dependent fashion. We developed a cell-based crosslinking approach to determine at high resolution and transcriptome-wide the binding sites of cellular RBPs and miRNPs. The crosslinked sites are revealed by thymidine to cytidine transitions in the cDNAs prepared from immunopurified RNPs of 4-thiouridine-treated cells. We determined the binding sites and regulatory consequences for several intensely studied RBPs and miRNPs, including PUM2, QKI, IGF2BP1-3, AGO/EIF2C1-4 and TNRC6A-C. Our study revealed that these factors bind thousands of sites containing defined sequence motifs and have distinct preferences for exonic versus intronic or coding versus untranslated transcript regions. The precise mapping of binding sites across the transcriptome will be critical to the interpretation of the rapidly emerging data on genetic variation between individuals and how these variations contribute to complex genetic diseases.

Discriminating Neoplastic from Nonneoplastic Tissues Using an miRNA-Based Deep Cancer Classifier.

Next-generation sequencing has enabled the collection of large biological data sets, allowing novel molecular-based classification methods to be developed for increased understanding of disease. miRNAs are small regulatory RNA molecules that can be quantified using next-generation sequencing and are excellent classificatory markers. Herein, a deep cancer classifier (DCC) was adapted to differentiate neoplastic from nonneoplastic samples using comprehensive miRNA expression profiles from 1031 human breast and skin tissue samples. The classifier was fine-tuned and evaluated using 750 neoplastic and 281 nonneoplastic breast and skin tissue samples. Performance of the DCC was compared with two machine-learning classifiers: support vector machine and random forests. In addition, performance of feature extraction through the DCC was also compared with a developed feature selection algorithm, cancer specificity. The DCC had the highest performance of area under the receiver operating curve and high performance in both sensitivity and specificity, unlike machine-learning and feature selection models, which often performed well in one metric compared with the other. In particular, deep learning had noticeable advantages with highly heterogeneous data sets. In addition, our cancer specificity algorithm identified candidate biomarkers for differentiating neoplastic and nonneoplastic tissue samples (eg, miR-144 and miR-375 in breast cancer and miR-375 and miR-451 in skin cancer).

DND1 maintains germline stem cells via recruitment of the CCR4-NOT complex to target mRNAs.

The vertebrate-conserved RNA-binding protein DND1 is required for the survival of primordial germ cells (PGCs), as well as the suppression of germ cell tumours in mice. Here we show that in mice DND1 binds a UU(A/U) trinucleotide motif predominantly in the 3' untranslated regions of mRNA, and destabilizes target mRNAs through direct recruitment of the CCR4-NOT deadenylase complex. Transcriptomic analysis reveals that the extent of suppression is dependent on the number of DND1-binding sites. This DND1-dependent mRNA destabilization is required for the survival of mouse PGCs and spermatogonial stem cells by suppressing apoptosis. The spectrum of target RNAs includes positive regulators of apoptosis and inflammation, and modulators of signalling pathways that regulate stem-cell pluripotency, including the TGFß superfamily, all of which are aberrantly elevated in DND1-deficient PGCs. We propose that the induction of the post-transcriptional suppressor DND1 synergizes with concurrent transcriptional changes to ensure precise developmental transitions during cellular differentiation and maintenance of the germ line.

AUF1 promotes let-7b loading on Argonaute 2.

Eukaryotic gene expression is tightly regulated post-transcriptionally by RNA-binding proteins (RBPs) and microRNAs. The RBP AU-rich-binding factor 1 (AUF1) isoform p37 was found to have high affinity for the microRNA let-7b in vitro (Kd = ∼ 6 nM) in cells. Ribonucleoprotein immunoprecipitation, in vitro association, and single-molecule-binding analyses revealed that AUF1 promoted let-7b loading onto Argonaute 2 (AGO2), the catalytic component of the RNA-induced silencing complex (RISC). In turn, AGO2-let-7 triggered target mRNA decay. Our findings uncover a novel mechanism by which AUF1 binding and transfer of microRNA let-7 to AGO2 facilitates let-7-elicited gene silencing.

Topology preserving stratification of tissue neoplasticity using Deep Neural Maps and microRNA signatures.

BACKGROUND: Accurate cancer classification is essential for correct treatment selection and better prognostication. microRNAs (miRNAs) are small RNA molecules that negatively regulate gene expression, and their dyresgulation is a common disease mechanism in many cancers. Through a clearer understanding of miRNA dysregulation in cancer, improved mechanistic knowledge and better treatments can be sought. RESULTS: We present a topology-preserving deep learning framework to study miRNA dysregulation in cancer. Our study comprises miRNA expression profiles from 3685 cancer and non-cancer tissue samples and hierarchical annotations on organ and neoplasticity status. Using unsupervised learning, a two-dimensional topological map is trained to cluster similar tissue samples. Labelled samples are used after training to identify clustering accuracy in terms of tissue-of-origin and neoplasticity status. In addition, an approach using activation gradients is developed to determine the attention of the networks to miRNAs that drive the clustering. Using this deep learning framework, we classify the neoplasticity status of held-out test samples with an accuracy of 91.07%, the tissue-of-origin with 86.36%, and combined neoplasticity status and tissue-of-origin with an accuracy of 84.28%. The topological maps display the ability of miRNAs to recognize tissue types and neoplasticity status. Importantly, when our approach identifies samples that do not cluster well with their respective classes, activation gradients provide further insight in cancer subtypes or grades. CONCLUSIONS: An unsupervised deep learning approach is developed for cancer classification and interpretation. This work provides an intuitive approach for understanding molecular properties of cancer and has significant potential for cancer classification and treatment selection.

Non-reversible tissue fixation retains extracellular vesicles for in situ imaging.

Extracellular vesicles (EVs) are secreted nanosized particles with many biological functions and pathological associations. The inability to image EVs in fixed tissues has been a major limitation to understanding their role in healthy and diseased tissue microenvironments. Here, we show that crosslinking mammalian tissues with formaldehyde results in significant EV loss, which can be prevented by additional fixation with 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide (EDC) for visualization of EVs in a range of normal and cancer tissues.

Co-repressor CBFA2T2 regulates pluripotency and germline development.

Developmental specification of germ cells lies at the heart of inheritance, as germ cells contain all of the genetic and epigenetic information transmitted between generations. The critical developmental event distinguishing germline from somatic lineages is the differentiation of primordial germ cells (PGCs), precursors of sex-specific gametes that produce an entire organism upon fertilization. Germ cells toggle between uni- and pluripotent states as they exhibit their own 'latent' form of pluripotency. For example, PGCs express a number of transcription factors in common with embryonic stem (ES) cells, including OCT4 (encoded by Pou5f1), SOX2, NANOG and PRDM14 (refs 2, 3, 4). A biochemical mechanism by which these transcription factors converge on chromatin to produce the dramatic rearrangements underlying ES-cell- and PGC-specific transcriptional programs remains poorly understood. Here we identify a novel co-repressor protein, CBFA2T2, that regulates pluripotency and germline specification in mice. Cbfa2t2(-/-) mice display severe defects in PGC maturation and epigenetic reprogramming. CBFA2T2 forms a biochemical complex with PRDM14, a germline-specific transcription factor. Mechanistically, CBFA2T2 oligomerizes to form a scaffold upon which PRDM14 and OCT4 are stabilized on chromatin. Thus, in contrast to the traditional 'passenger' role of a co-repressor, CBFA2T2 functions synergistically with transcription factors at the crossroads of the fundamental developmental plasticity between uni- and pluripotency.

Human cGAS catalytic domain has an additional DNA-binding interface that enhances enzymatic activity and liquid-phase condensation.

The cyclic GMP-AMP synthase (cGAS)-cGAMP-STING pathway plays a key role in innate immunity, with cGAS sensing both pathogenic and mislocalized DNA in the cytoplasm. Human cGAS (h-cGAS) constitutes an important drug target for control of antiinflammatory responses that can contribute to the onset of autoimmune diseases. Recent studies have established that the positively charged N-terminal segment of cGAS contributes to enhancement of cGAS enzymatic activity as a result of DNA-induced liquid-phase condensation. We have identified an additional cGASCD-DNA interface (labeled site-C; CD, catalytic domain) in the crystal structure of a human SRY.cGASCD-DNA complex, with mutations along this basic site-C cGAS interface disrupting liquid-phase condensation, as monitored by cGAMP formation, gel shift, spin-down, and turbidity assays, as well as time-lapse imaging of liquid droplet formation. We expand on an earlier ladder model of cGAS dimers bound to a pair of parallel-aligned DNAs to propose a multivalent interaction-mediated cluster model to account for DNA-mediated condensation involving both the N-terminal domain of cGAS and the site-C cGAS-DNA interface. We also report the crystal structure of the h-cGASCD-DNA complex containing a triple mutant that disrupts the site-C interface, with this complex serving as a future platform for guiding cGAS inhibitor development at the DNA-bound h-cGAS level. Finally, we solved the structure of RU.521 bound in two alternate alignments to apo h-cGASCD, thereby occupying more of the catalytic pocket and providing insights into further optimization of active-site-binding inhibitors.

In vitro antiviral activity of the anti-HCV drugs daclatasvir and sofosbuvir against SARS-CoV-2, the aetiological agent of COVID-19.

BACKGROUND: Current approaches of drug repurposing against COVID-19 have not proven overwhelmingly successful and the SARS-CoV-2 pandemic continues to cause major global mortality. SARS-CoV-2 nsp12, its RNA polymerase, shares homology in the nucleotide uptake channel with the HCV orthologue enzyme NS5B. Besides, HCV enzyme NS5A has pleiotropic activities, such as RNA binding, that are shared with various SARS-CoV-2 proteins. Thus, anti-HCV NS5B and NS5A inhibitors, like sofosbuvir and daclatasvir, respectively, could be endowed with anti-SARS-CoV-2 activity. METHODS: SARS-CoV-2-infected Vero cells, HuH-7 cells, Calu-3 cells, neural stem cells and monocytes were used to investigate the effects of daclatasvir and sofosbuvir. In silico and cell-free based assays were performed with SARS-CoV-2 RNA and nsp12 to better comprehend the mechanism of inhibition of the investigated compounds. A physiologically based pharmacokinetic model was generated to estimate daclatasvir's dose and schedule to maximize the probability of success for COVID-19. RESULTS: Daclatasvir inhibited SARS-CoV-2 replication in Vero, HuH-7 and Calu-3 cells, with potencies of 0.8, 0.6 and 1.1 µM, respectively. Although less potent than daclatasvir, sofosbuvir alone and combined with daclatasvir inhibited replication in Calu-3 cells. Sofosbuvir and daclatasvir prevented virus-induced neuronal apoptosis and release of cytokine storm-related inflammatory mediators, respectively. Sofosbuvir inhibited RNA synthesis by chain termination and daclatasvir targeted the folding of secondary RNA structures in the SARS-CoV-2 genome. Concentrations required for partial daclatasvir in vitro activity are achieved in plasma at Cmax after administration of the approved dose to humans. CONCLUSIONS: Daclatasvir, alone or in combination with sofosbuvir, at higher doses than used against HCV, may be further fostered as an anti-COVID-19 therapy.

Plasma microRNA Interindividual Variability in Healthy Individuals, Pregnant Women, and an Individual with a Stably Altered Neuroendocrine Phenotype.

BACKGROUND: Extracellular RNAs (exRNAs) in biofluids are amenable to quantitative analysis and proposed as noninvasive biomarkers for monitoring organ function. Cell-lineage-specific microRNAs (miRNAs) are present in plasma as soluble ribonucleoproteins or enclosed in exRNA carriers and transported through the vasculature. However, more extensive studies of healthy individuals are needed to gain insights into the variability of plasma miRNA abundance and composition. METHODS: The exRNA composition of platelet-depleted plasma collected twice from 236 healthy individuals was characterized by small RNA sequencing. Plasma of pregnant women featuring dramatically increased placental miRNAs and samples from subject P12 with noticeably increased epithelial- and neuroendocrine-origin miRNAs were included for comparison. The miRNA content of 10 000g and 100 000g pellet fractions of plasma generated by ultracentrifugation was also determined. Data analysis methods included Pearson correlation, differential gene expression, and unsupervised clustering. RESULTS: The abundance changes for more variable miRNAs in plasma of normal individuals correlated between coexpressed cell-lineage-specific miRNAs of the liver, neuroendocrine organs, epithelial cells, and muscle. ExRNA of pellet fractions contained <2% of total plasma miRNA with modest enrichment of lineage-specific and variable miRNAs compared to supernatant. The abundance fold changes of miRNAs observed in pregnancy and P12 compared to normal exceeded interquartile variability of healthy individuals. The neuroendocrine miRNA signature of P12 persisted for more than 4 years and was absent in other individuals. CONCLUSIONS: This study defines the framework and effect size for screening of extensive plasma collections for miRNA phenotypes and biomarker discovery.