A role for the C. elegans Argonaute protein CSR-1 in small nuclear RNA 3' processing.

The Integrator is a multi-subunit protein complex that catalyzes the maturation of snRNA transcripts via 3' cleavage, a step required for snRNA incorporation with snRNP for spliceosome biogenesis. Here we developed a GFP based in vivo snRNA misprocessing reporter as a readout of Integrator function and performed a genome-wide RNAi screen for Integrator regulators. We found that loss of the Argonaute encoding csr-1 gene resulted in widespread 3' misprocessing of snRNA transcripts that is accompanied by a significant increase in alternative splicing. Loss of the csr-1 gene down-regulates the germline expression of Integrator subunits 4 and 6 and is accompanied by a reduced protein translation efficiency of multiple Integrator catalytic and non-catalytic subunits. Through isoform and motif mutant analysis, we determined that CSR-1's effect on snRNA processing is dependent on its catalytic slicer activity but does not involve the CSR-1a isoform. Moreover, mRNA-sequencing revealed high similarity in the transcriptome profile between csr-1 and Integrator subunit knockdown via RNAi. Together, our findings reveal CSR-1 as a new regulator of the Integrator complex and implicate a novel role of this Argonaute protein in snRNA 3' processing.

PRPF8-mediated dysregulation of hBrr2 helicase disrupts human spliceosome kinetics and 5´-splice-site selection causing tissue-specific defects.

The carboxy-terminus of the spliceosomal protein PRPF8, which regulates the RNA helicase Brr2, is a hotspot for mutations causing retinitis pigmentosa-type 13, with unclear role in human splicing and tissue-specificity mechanism. We used patient induced pluripotent stem cells-derived cells, carrying the heterozygous PRPF8 c.6926 A > C (p.H2309P) mutation to demonstrate retinal-specific endophenotypes comprising photoreceptor loss, apical-basal polarity and ciliary defects. Comprehensive molecular, transcriptomic, and proteomic analyses revealed a role of the PRPF8/Brr2 regulation in 5'-splice site (5'SS) selection by spliceosomes, for which disruption impaired alternative splicing and weak/suboptimal 5'SS selection, and enhanced cryptic splicing, predominantly in ciliary and retinal-specific transcripts. Altered splicing efficiency, nuclear speckles organisation, and PRPF8 interaction with U6 snRNA, caused accumulation of active spliceosomes and poly(A)+ mRNAs in unique splicing clusters located at the nuclear periphery of photoreceptors. Collectively these elucidate the role of PRPF8/Brr2 regulatory mechanisms in splicing and the molecular basis of retinal disease, informing therapeutic approaches.

Structural insights into human exon-defined spliceosome prior to activation.

Spliceosome is often assembled across an exon and undergoes rearrangement to span a neighboring intron. Most states of the intron-defined spliceosome have been structurally characterized. However, the structure of a fully assembled exon-defined spliceosome remains at large. During spliceosome assembly, the pre-catalytic state (B complex) is converted from its precursor (pre-B complex). Here we report atomic structures of the exon-defined human spliceosome in four sequential states: mature pre-B, late pre-B, early B, and mature B. In the previously unknown late pre-B state, U1 snRNP is already released but the remaining proteins are still in the pre-B state; unexpectedly, the RNAs are in the B state, with U6 snRNA forming a duplex with 5'-splice site and U5 snRNA recognizing the 3'-end of the exon. In the early and mature B complexes, the B-specific factors are stepwise recruited and specifically recognize the exon 3'-region. Our study reveals key insights into the assembly of the exon-defined spliceosomes and identifies mechanistic steps of the pre-B-to-B transition.

Distinct functions for the paralogous RBM41 and U11/U12-65K proteins in the minor spliceosome.

Here, we identify RBM41 as a novel unique protein component of the minor spliceosome. RBM41 has no previously recognized cellular function but has been identified as a paralog of U11/U12-65K, a known unique component of the U11/U12 di-snRNP. Both proteins use their highly similar C-terminal RRMs to bind to 3'-terminal stem-loops in U12 and U6atac snRNAs with comparable affinity. Our BioID data indicate that the unique N-terminal domain of RBM41 is necessary for its association with complexes containing DHX8, an RNA helicase, which in the major spliceosome drives the release of mature mRNA from the spliceosome. Consistently, we show that RBM41 associates with excised U12-type intron lariats, is present in the U12 mono-snRNP, and is enriched in Cajal bodies, together suggesting that RBM41 functions in the post-splicing steps of the minor spliceosome assembly/disassembly cycle. This contrasts with U11/U12-65K, which uses its N-terminal region to interact with U11 snRNP during intron recognition. Finally, while RBM41 knockout cells are viable, they show alterations in U12-type 3' splice site usage. Together, our results highlight the role of the 3'-terminal stem-loop of U12 snRNA as a dynamic binding platform for the U11/U12-65K and RBM41 proteins, which function at distinct stages of the assembly/disassembly cycle.

2'-O-methylation (Nm) in RNA: progress, challenges, and future directions.

RNA 2'-O-methylation (Nm) is highly abundant in noncoding RNAs including ribosomal RNA (rRNA), transfer RNA (tRNA), and small nuclear RNA (snRNA), and occurs in the 5' cap of virtually all messenger RNAs (mRNAs) in higher eukaryotes. More recently, Nm has also been reported to occur at internal sites in mRNA. High-throughput methods have been developed for the transcriptome-wide detection of Nm. However, these methods have mostly been applied to abundant RNAs such as rRNA, and the validity of the internal mRNA Nm sites detected with these approaches remains controversial. Nonetheless, Nm in both coding and noncoding RNAs has been demonstrated to impact cellular processes, including translation and splicing. In addition, Nm modifications at the 5' cap and possibly at internal sites in mRNA serve to prevent the binding of nucleic acid sensors, thus preventing the activation of the innate immune response by self-mRNAs. Finally, Nm has been implicated in a variety of diseases including cancer, cardiovascular diseases, and neurologic syndromes. In this review, we discuss current challenges in determining the distribution, regulation, function, and disease relevance of Nm, as well as potential future directions for the field.

Unconventional features in the transcription and processing of spliceosomal small nuclear RNAs in the protozoan parasite Trichomonas vaginalis.

Trichomonas vaginalis is a medically important protozoan parasite, and a deep-branching, evolutionarily divergent unicellular eukaryote that has conserved several key features of eukaryotic gene expression. Trichomonas vaginalis possesses a metazoan/plant-like capping apparatus, mRNAs with a cap 1 structure and spliceosomes containing the five small nuclear RNAs (snRNAs). However, in contrast to metazoan and plant snRNAs, the structurally conserved T. vaginalis snRNAs were initially identified as lacking the canonical guanosine cap nucleotide. To explain this unusual condition, we sought to investigate transcriptional and processing features of the spliceosomal snRNAs in this protist. Here, we show that T. vaginalis spliceosomal snRNA genes mostly lack typical eukaryotic promoters. In contrast to other eukaryotes, the putative TATA box in the T. vaginalis U6 snRNA gene was found to be dispensable for transcription or RNA polymerase selectivity. Moreover, U6 transcription in T. vaginalis was virtually insensitive to tagetitoxin compared with other cellular transcripts produced by the same RNA polymerase III. Most important and unexpected, snRNA transcription in T. vaginalis appears to bypass capping as we show that these transcripts retain their original 5'-triphosphate groups. In conclusion, transcription and processing of spliceosomal snRNAs in T. vaginalis deviate considerably from the conventional rules of other eukaryotes.

Studying plant vascular development using single-cell approaches.

Vascular cells form a highly complex and heterogeneous tissue. Its composition, function, shape, and arrangement vary with the developmental stage and between organs and species. Understanding the transcriptional regulation underpinning this complexity thus requires a high-resolution technique that is capable of capturing rapid events during vascular cell formation. Single-cell and single-nucleus RNA sequencing (sc/snRNA-seq) approaches provide powerful tools to extract transcriptional information from these lowly abundant and dynamically changing cell types, which allows the reconstruction of developmental trajectories. Here, we summarize and reflect on recent studies using single-cell transcriptomics to study vascular cell types and discuss current and future implementations of sc/snRNA-seq approaches in the field of vascular development.

A role for SNU66 in maintaining 5' splice site identity during spliceosome assembly.

In spliceosome assembly, the 5' splice site is initially recognized by U1 snRNA. U1 leaves the spliceosome during the assembly process, therefore other factors contribute to the maintenance of 5' splice site identity as it is loaded into the catalytic site. Recent structural data suggest that human tri-snRNP 27K (SNRP27) M141 and SNU66 H734 interact to stabilize the U4/U6 quasi-pseudo knot at the base of the U6 snRNA ACAGAGA box in pre-B complex. Previously, we found that mutations in Caenorhabditis elegans at SNRP-27 M141 promote changes in alternative 5'ss usage. We tested whether the potential interaction between SNRP-27 M141 and SNU-66 H765 (the C. elegans equivalent position to human SNU66 H734) contributes to maintaining 5' splice site identity during spliceosome assembly. We find that SNU-66 H765 mutants promote alternative 5' splice site usage. Many of the alternative 5' splicing events affected by SNU-66(H765G) overlap with those affected SNRP-27(M141T). Double mutants of snrp-27(M141T) and snu-66(H765G) are homozygous lethal. We hypothesize that mutations at either SNRP-27 M141 or SNU-66 H765 allow the spliceosome to load alternative 5' splice sites into the active site. Tests with mutant U1 snRNA and swapped 5' splice sites indicate that the ability of SNRP-27 M141 and SNU-66 H765 mutants to affect a particular 5' splice alternative splicing event is dependent on both the presence of a weaker consensus 5'ss nearby and potentially nearby splicing factor binding sites. Our findings confirm a new role for the C terminus of SNU-66 in maintenance of 5' splice site identity during spliceosome assembly.

TBX3 reciprocally controls key trophoblast lineage decisions in villi during human placenta development in the first trimester.

Human trophoblastic lineage development is intertwined with placental development and pregnancy outcomes, but the regulatory mechanisms underpinning this process remain inadequately understood. In this study, based on single-nuclei RNA sequencing (snRNA-seq) analysis of the human early maternal-fetal interface, we compared the gene expression pattern of trophoblast at different developmental stages. Our findings reveal a predominant upregulation of TBX3 during the transition from villous cytotrophoblast (VCT) to syncytiotrophoblast (SCT), but downregulation of TBX3 as VCT progresses into extravillous trophoblast cells (EVT). Immunofluorescence analysis verified the primary expression of TBX3 in SCT, partial expression in MKi67-positive VCT, and absence in HLA-G-positive EVT, consistent with our snRNA-seq results. Using immortalized trophoblastic cell lines (BeWo and HTR8/SVneo) and human primary trophoblast stem cells (hTSCs), we observed that TBX3 knockdown impedes SCT formation through RAS-MAPK signaling, while TBX3 overexpression disrupts the cytoskeleton structure of EVT and hinders EVT differentiation by suppressing FAK signaling. In conclusion, our study suggests that the spatiotemporal expression of TBX3 plays a critical role in regulating trophoblastic lineage development via distinct signaling pathways. This underscores TBX3 as a key determinant during hemochorial placental development.

Disease-Associated Neurotoxic Astrocyte Markers in Alzheimer Disease Based on Integrative Single-Nucleus RNA Sequencing.

Alzheimer disease (AD) is an irreversible neurodegenerative disease, and astrocytes play a key role in its onset and progression. The aim of this study is to analyze the characteristics of neurotoxic astrocytes and identify novel molecular targets for slowing down the progression of AD. Single-nucleus RNA sequencing (snRNA-seq) data were analyzed from various AD cohorts comprising about 210,654 cells from 53 brain tissue. By integrating snRNA-seq data with bulk RNA-seq data, crucial astrocyte types and genes associated with the prognosis of patients with AD were identified. The expression of neurotoxic astrocyte markers was validated using 5 × FAD and wild-type (WT) mouse models, combined with experiments such as western blot, quantitative real-time PCR (qRT-PCR), and immunofluorescence. A group of neurotoxic astrocytes closely related to AD pathology was identified, which were involved in inflammatory responses and pathways related to neuron survival. Combining snRNA and bulk tissue data, ZEP36L, AEBP1, WWTR1, PHYHD1, DST and RASL12 were identified as toxic astrocyte markers closely related to disease severity, significantly elevated in brain tissues of 5 × FAD mice and primary astrocytes treated with Aß. Among them, WWTR1 was significantly increased in astrocytes of 5 × FAD mice, driving astrocyte inflammatory responses, and has been identified as an important marker of neurotoxic astrocytes. snRNA-seq analysis reveals the biological functions of neurotoxic astrocytes. Six genes related to AD pathology were identified and validated, among which WWTR1 may be a novel marker of neurotoxic astrocytes.

Structural insights into branch site proofreading by human spliceosome.

Selection of the pre-mRNA branch site (BS) by the U2 small nuclear ribonucleoprotein (snRNP) is crucial to prespliceosome (A complex) assembly. The RNA helicase PRP5 proofreads BS selection but the underlying mechanism remains unclear. Here we report the atomic structures of two sequential complexes leading to prespliceosome assembly: human 17S U2 snRNP and a cross-exon pre-A complex. PRP5 is anchored on 17S U2 snRNP mainly through occupation of the RNA path of SF3B1 by an acidic loop of PRP5; the helicase domain of PRP5 associates with U2 snRNA; the BS-interacting stem-loop (BSL) of U2 snRNA is shielded by TAT-SF1, unable to engage the BS. In the pre-A complex, an initial U2-BS duplex is formed; the translocated helicase domain of PRP5 stays with U2 snRNA and the acidic loop still occupies the RNA path. The pre-A conformation is specifically stabilized by the splicing factors SF1, DNAJC8 and SF3A2. Cancer-derived mutations in SF3B1 damage its association with PRP5, compromising BS proofreading. Together, these findings reveal key insights into prespliceosome assembly and BS selection or proofreading by PRP5.

THUMPD2 catalyzes the N2-methylation of U6 snRNA of the spliceosome catalytic center and regulates pre-mRNA splicing and retinal degeneration.

The mechanisms by which the relatively conserved spliceosome manages the enormously large number of splicing events that occur in humans (∼200 000 versus ∼300 in yeast) are poorly understood. Here, we show deposition of one RNA modification-N2-methylguanosine (m2G) on the G72 of U6 snRNA (the catalytic center of the spliceosome) promotes efficient pre-mRNA splicing activity in human cells. This modification was identified to be conserved among vertebrates. Further, THUMPD2 was demonstrated as the methyltransferase responsible for U6 m2G72 by explicitly recognizing the U6-specific sequences and structural elements. The knock-out of THUMPD2 eliminated U6 m2G72 and impaired the pre-mRNA splicing activity, resulting in thousands of changed alternative splicing events of endogenous pre-mRNAs in human cells. Notably, the aberrantly spliced pre-mRNA population elicited the nonsense-mediated mRNA decay pathway. We further show that THUMPD2 was associated with age-related macular degeneration and retinal function. Our study thus demonstrates how an RNA epigenetic modification of the major spliceosome regulates global pre-mRNA splicing and impacts physiology and disease.

Examining the capacity of human U1 snRNA variants to facilitate pre-mRNA splicing.

The human U1 snRNA is encoded by a multigene family consisting of transcribed variants and defective pseudogenes. Many variant U1 (vU1) snRNAs have been demonstrated to not only be transcribed but also processed by the addition of a trimethylated guanosine cap, packaged into snRNPs, and assembled into spliceosomes; however, their capacity to facilitate pre-mRNA splicing has, so far, not been tested. A recent systematic analysis of the human snRNA genes identified 178 U1 snRNA genes that are present in the genome as either tandem arrays or single genes on multiple chromosomes. Of these, 15 were found to be expressed in human tissues and cell lines, although at significantly low levels from their endogenous loci, <0.001% of the canonical U1 snRNA. In this study, we found that placing the variants in the context of the regulatory elements of the RNU1-1 gene improves the expression of many variants to levels comparable to the canonical U1 snRNA. Application of a previously established HeLa cell-based minigene reporter assay to examine the capacity of the vU1 snRNAs to support pre-mRNA splicing revealed that even though the exogenously expressed variant snRNAs were enriched in the nucleus, only a few had a measurable effect on splicing.

Global analysis of binding sites of U2AF1 and ZRSR2 reveals RNA elements required for mutually exclusive splicing by the U2- and U12-type spliceosome.

Recurring mutations in genes encoding 3' splice-site recognition proteins, U2AF1 and ZRSR2 are associated with human cancers. Here, we determined binding sites of the proteins to reveal that U2-type and U12-type splice sites are recognized by U2AF1 and ZRSR2, respectively. However, some sites are spliced by both the U2-type and U12-type spliceosomes, indicating that well-conserved consensus motifs in some U12-type introns could be recognized by the U2-type spliceosome. Nucleotides flanking splice sites of U12-type introns are different from those flanking U2-type introns. Remarkably, the AG dinucleotide at the positions -1 and -2 of 5' splice sites of U12-type introns with GT-AG termini is not present. AG next to 5' splice site introduced by a single nucleotide substitution at the -2 position could convert a U12-type splice site to a U2-type site. The class switch of introns by a single mutation and the bias against G at the -1 position of U12-type 5' splice site support the notion that the identities of nucleotides in exonic regions adjacent to splice sites are fine-tuned to avoid recognition by the U2-type spliceosome. These findings may shed light on the mechanism of selectivity in U12-type intron splicing and the mutations that affect splicing.

Disrupted hypothalamic transcriptomics and proteomics in a mouse model of type 2 diabetes exposed to recurrent hypoglycaemia.

AIMS/HYPOTHESIS: Repeated exposures to insulin-induced hypoglycaemia in people with diabetes progressively impairs the counterregulatory response (CRR) that restores normoglycaemia. This defect is characterised by reduced secretion of glucagon and other counterregulatory hormones. Evidence indicates that glucose-responsive neurons located in the hypothalamus orchestrate the CRR. Here, we aimed to identify the changes in hypothalamic gene and protein expression that underlie impaired CRR in a mouse model of defective CRR. METHODS: High-fat-diet fed and low-dose streptozocin-treated C57BL/6N mice were exposed to one (acute hypoglycaemia [AH]) or multiple (recurrent hypoglycaemia [RH]) insulin-induced hypoglycaemic episodes and plasma glucagon levels were measured. Single-nuclei RNA-seq (snRNA-seq) data were obtained from the hypothalamus and cortex of mice exposed to AH and RH. Proteomic data were obtained from hypothalamic synaptosomal fractions. RESULTS: The final insulin injection resulted in similar plasma glucose levels in the RH group and AH groups, but glucagon secretion was significantly lower in the RH group (AH: 94.5±9.2 ng/l [n=33]; RH: 59.0±4.8 ng/l [n=37]; p<0.001). Analysis of snRNA-seq data revealed similar proportions of hypothalamic cell subpopulations in the AH- and RH-exposed mice. Changes in transcriptional profiles were found in all cell types analysed. In neurons from RH-exposed mice, we observed a significant decrease in expression of Avp, Pmch and Pcsk1n, and the most overexpressed gene was Kcnq1ot1, as compared with AH-exposed mice. Gene ontology analysis of differentially expressed genes (DEGs) indicated a coordinated decrease in many oxidative phosphorylation genes and reduced expression of vacuolar H+- and Na+/K+-ATPases; these observations were in large part confirmed in the proteomic analysis of synaptosomal fractions. Compared with AH-exposed mice, oligodendrocytes from RH-exposed mice had major changes in gene expression that suggested reduced myelin formation. In astrocytes from RH-exposed mice, DEGs indicated reduced capacity for neurotransmitters scavenging in tripartite synapses as compared with astrocytes from AH-exposed mice. In addition, in neurons and astrocytes, multiple changes in gene expression suggested increased amyloid beta (Aß) production and stability. The snRNA-seq analysis of the cortex showed that the adaptation to RH involved different biological processes from those seen in the hypothalamus. CONCLUSIONS/INTERPRETATION: The present study provides a model of defective counterregulation in a mouse model of type 2 diabetes. It shows that repeated hypoglycaemic episodes induce multiple defects affecting all hypothalamic cell types and their interactions, indicative of impaired neuronal network signalling and dysegulated hypoglycaemia sensing, and displaying features of neurodegenerative diseases. It also shows that repeated hypoglycaemia leads to specific molecular adaptation in the hypothalamus when compared with the cortex. DATA AVAILABILITY: The transcriptomic dataset is available via the GEO ( http://www.ncbi.nlm.nih.gov/geo/ ), using the accession no. GSE226277. The proteomic dataset is available via the ProteomeXchange data repository ( http://www.proteomexchange.org ), using the accession no. PXD040183.

Transcriptome profile of subsynaptic myonuclei at the neuromuscular junction in embryogenesis.

Skeletal muscle fiber is a large syncytium with multiple and evenly distributed nuclei. Adult subsynaptic myonuclei beneath the neuromuscular junction (NMJ) express specific genes, the products of which coordinately function in the maintenance of the pre- and post-synaptic regions. However, the gene expression profiles that promote the NMJ formation during embryogenesis remain largely unexplored. We performed single-nucleus RNA sequencing (snRNA-seq) analysis of embryonic and neonatal mouse diaphragms, and found that each myonucleus had a distinct transcriptome pattern during the NMJ formation. Among the previously reported NMJ-constituting genes, Dok7, Chrna1, and Chrnd are specifically expressed in subsynaptic myonuclei at E18.5. In the E18.5 diaphragm, ca. 10.7% of the myonuclei express genes for the NMJ formation (Dok7, Chrna1, and Chrnd) together with four representative ß-catenin regulators (Amotl2, Ptprk, Fam53b, and Tcf7l2). Additionally, the temporal gene expression patterns of these seven genes are synchronized in differentiating C2C12 myoblasts. Amotl2 and Ptprk are expressed in the sarcoplasm, where ß-catenin serves as a structural protein to organize the membrane-anchored NMJ structure. In contrast, Fam53b and Tcf7l2 are expressed in the myonucleus, where ß-catenin serves as a transcriptional coactivator in Wnt/ß-catenin signaling at the NMJ. In C2C12 myotubes, knockdown of Amotl2 or Ptprk markedly, and that of Fam53b and Tcf7l2 less efficiently, impair the clustering of acetylcholine receptors. In contrast, knockdown of Fam53b and Tcf7l2, but not of Amotl2 or Ptprk, impairs the gene expression of Slit2 encoding an axonal attractant for motor neurons, which is required for the maturation of motor nerve terminal. Thus, Amotl2 and Ptprk exert different roles at the NM compared to Fam53b and Tcf7l2. Additionally, Wnt ligands originating from the spinal motor neurons and the perichondrium/chondrocyte are likely to work remotely on the subsynaptic nuclei and the myotendinous junctional nuclei, respectively. We conclude that snRNA-seq analysis of embryonic/neonatal diaphragms reveal a novel coordinated expression profile especially in the Wnt/ß-catenin signaling that regulate the formation of the embryonic NMJ.

snRNA-seq analysis in multinucleated myogenic FSHD cells identifies heterogeneous FSHD transcriptome signatures associated with embryonic-like program activation and oxidative stress-induced apoptosis.

The sporadic nature of DUX4 expression in FSHD muscle challenges comparative transcriptome analyses between FSHD and control samples. A variety of DUX4 and FSHD-associated transcriptional changes have been identified, but bulk RNA-seq strategies prohibit comprehensive analysis of their spatiotemporal relation, interdependence and role in the disease process. In this study, we used single-nucleus RNA-sequencing of nuclei isolated from patient- and control-derived multinucleated primary myotubes to investigate the cellular heterogeneity in FSHD. Taking advantage of the increased resolution in snRNA-sequencing of fully differentiated myotubes, two distinct populations of DUX4-affected nuclei could be defined by their transcriptional profiles. Our data provides insights into the differences between these two populations and suggests heterogeneity in two well-known FSHD-associated transcriptional aberrations: increased oxidative stress and inhibition of myogenic differentiation. Additionally, we provide evidence that DUX4-affected nuclei share transcriptome features with early embryonic cells beyond the well-described cleavage stage, progressing into the 8-cell and blastocyst stages. Altogether, our data suggests that the FSHD transcriptional profile is defined by a mixture of individual and sometimes mutually exclusive DUX4-induced responses and cellular state-dependent downstream effects.

Exploring the highly reduced spliceosome of Pseudoloma neurophilia.

Spliceosomal introns evolved early in eukaryogenesis, originating from self-splicing group II introns that invaded the proto-eukaryotic genome1. Elements of these ribozymes, now called snRNAs (U1, U2, U4, U5, U6), were co-opted to excise these invasive elements. Prior to eukaryotic diversification, the spliceosome is predicted to have accumulated hundreds of proteins2. This early complexification has obscured our understanding of spliceosomal evolution. Reduced systems with few introns and tiny spliceosomes give insights into the plasticity of the splicing reaction and provide an opportunity to study the evolution of the spliceosome3,4. Microsporidia are intracellular parasites possessing extremely reduced genomes that have lost many, and in some instances all, introns5. In the purportedly intron-lacking genome of the microsporidian Pseudoloma neurophilia6, we identified two introns that are spliced at high levels. Furthermore, with only 14 predicted proteins, the P. neurophilia spliceosome could be the smallest known. Intriguingly, the few proteins retained are divergent compared to canonical orthologs. Even the central spliceosomal protein Prp8, which originated from the proteinaceous component of group II introns, is extremely divergent. This is unusual given that Prp8 is highly conserved across eukaryotes, including other microsporidia. All five P. neurophilia snRNAs are present, and all but U2 have diverged extensively, likely resulting from the loss of interacting proteins. Despite this divergence, U1 and U2 are predicted to pair with intron sequences more extensively than previously described. The P. neurophilia spliceosome is retained to splice a mere two introns and, with few proteins and reliance on RNA-RNA interactions, could function in a manner more reminiscent of presumed ancestral splicing.

Does decreased autophagy and dysregulation of LC3A in astrocytes play a role in major depressive disorder?

Astrocytic dysfunction contributes to the molecular pathogenesis of major depressive disorder (MDD). However, the astrocytic subtype that mainly contributes to MDD etiology and whether dysregulated autophagy in astrocytes is associated with MDD remain unknown. Using a single-nucleus RNA sequencing (snRNA-seq) atlas, three astrocyte subtypes were identified in MDD, while C2 State-1Q astrocytes showed aberrant changes in both cell proportion and most differentially expressed genes compared with other subtypes. Moreover, autophagy pathways were commonly inhibited in astrocytes in the prefrontal cortices (PFCs) of patients with MDD, especially in C2 State-1Q astrocytes. Furthermore, by integrating snRNA-seq and bulk transcriptomic data, we found significant reductions in LC3A expression levels in the PFC region of CUMS-induced depressed mice, as well as in postmortem PFC tissues and peripheral blood samples from patients with MDD. These results were further validated by qPCR using whole-blood samples from patients with MDD and healthy controls. Finally, LC3A expression in the whole blood of patients with MDD was negatively associated with the severity of depressive symptoms. Overall, our results underscore autophagy inhibition in PFC astrocytes as a common molecular characteristic in MDD and might reveal a novel potential diagnostic marker LC3A.

Deletions of singular U1 snRNA gene significantly interfere with transcription and 3'-end mRNA formation.

Small nuclear RNAs (snRNAs) are structural and functional cores of the spliceosome. In metazoan genomes, each snRNA has multiple copies/variants, up to hundreds in mammals. However, the expressions and functions of each copy/variant in one organism have not been systematically studied. Focus on U1 snRNA genes, we investigated all five copies in Drosophila melanogaster using two series of constructed strains. Analyses of transgenic flies that each have a U1 promoter-driven gfp revealed that U1:21D is the major and ubiquitously expressed copy, and the other four copies have specificities in developmental stages and tissues. Mutant strains that each have a precisely deleted copy of U1-gene exhibited various extents of defects in fly morphology or mobility, especially deletion of U1:82Eb. Interestingly, splicing was changed at limited levels in the deletion strains, while large amounts of differentially-expressed genes and alternative polyadenylation events were identified, showing preferences in the down-regulation of genes with 1-2 introns and selection of proximal sites for 3'-end polyadenylation. In vitro assays suggested that Drosophila U1 variants pulled down fewer SmD2 proteins compared to the canonical U1. This study demonstrates that all five U1-genes in Drosophila have physiological functions in development and play regulatory roles in transcription and 3'-end formation.