Around 60% of individuals with neurodevelopmental disorders (NDD) remain undiagnosed after comprehensive genetic testing, primarily of protein-coding genes¹. Large genome-sequenced cohorts are improving our ability to discover new diagnoses in the non-coding genome. Here we identify the non-coding RNA RNU4-2 as a syndromic NDD gene. RNU4-2 encodes the U4 small nuclear RNA (snRNA), which is a critical component of the U4/U6.U5 tri-snRNP complex of the major spliceosome². We identify an 18 base pair region of RNU4-2 mapping to two structural elements in the U4/U6 snRNA duplex (the T-loop and stem III) that is severely depleted of variation in the general population, but in which we identify heterozygous variants in 115 individuals with NDD. Most individuals (77.4%) have the same highly recurrent single base insertion (n.64_65insT). In 54 individuals in whom it could be determined, the de novo variants were all on the maternal allele. We demonstrate that RNU4-2 is highly expressed in the developing human brain, in contrast to RNU4-1 and other U4 homologues. Using RNA sequencing, we show how 5′ splice-site use is systematically disrupted in individuals with RNU4-2 variants, consistent with the known role of this region during spliceosome activation. Finally, we estimate that variants in this 18 base pair region explain 0.4% of individuals with NDD. This work underscores the importance of non-coding genes in rare disorders and will provide a diagnosis to thousands of individuals with NDD worldwide.

The major spliceosome contains five small nuclear RNAs (snRNAs; U1, U2, U4, U5 and U6) essential for splicing. Variants in RNU4-2, encoding U4, cause a neurodevelopmental disorder called ReNU syndrome. We investigated de novo variants in 50 snRNA-encoding genes in a French cohort of 23,649 individuals with rare disorders and gathered additional cases through international collaborations. Altogether, we identified 145 previously unreported probands with (likely) pathogenic variants in RNU4-2 and 21 individuals with de novo and/or recurrent variants in RNU5B-1 and RNU5A-1, encoding U5. Pathogenic variants typically arose de novo on the maternal allele and cluster in regions critical for splicing. RNU4-2 variants mainly localize to two structures, the stem III and T-loop/quasi-pseudoknot, which position the U6 ACAGAGA box for 5′ splice site recognition and associate with different phenotypic severity. RNU4-2 variants result in specific defects in alternative 5′ splice site usage and methylation patterns (episignatures) that correlate with variant location and clinical severity. This study establishes RNU5B-1 as a neurodevelopmental disorder gene, suggests RNU5A-1 as a strong candidate and highlights the role of de novo variants in snRNAs.

Mutations in the mitochondrial or nuclear genomes are associated with a diverse group of human disorders characterized by impaired mitochondrial respiration. Within this group, an increasing number of mutations have been identified in nuclear genes involved in mitochondrial RNA biology. The TEFM gene encodes the mitochondrial transcription elongation factor responsible for enhancing the processivity of mitochondrial RNA polymerase, POLRMT. We report for the first time that TEFM variants are associated with mitochondrial respiratory chain deficiency and a wide range of clinical presentations including mitochondrial myopathy with a treatable neuromuscular transmission defect. Mechanistically, we show muscle and primary fibroblasts from the affected individuals have reduced levels of promoter distal mitochondrial RNA transcripts. Finally, tefm knockdown in zebrafish embryos resulted in neuromuscular junction abnormalities and abnormal mitochondrial function, strengthening the genotype-phenotype correlation. Our study highlights that TEFM regulates mitochondrial transcription elongation and its defect results in variable, tissue-specific neurological and neuromuscular symptoms.

PNPT1 (PNPase—polynucleotide phosphorylase) is involved in multiple RNA processing functions in the mitochondria. Bi-allelic pathogenic PNPT1 variants cause heterogeneous clinical phenotypes affecting multiple organs without any established genotype–phenotype correlations. Defects in PNPase can cause variable combined respiratory chain complex defects. Recently, it has been suggested that PNPase can lead to activation of an innate immune response. To better understand the clinical and molecular spectrum of patients with bi-allelic PNPT1 variants, we captured detailed clinical and molecular phenotypes of all 17 patients reported in the literature, plus seven new patients, including a 78-year-old male with the longest reported survival. A functional follow-up of genomic sequencing by cDNA studies confirmed a splicing defect in a novel, apparently synonymous, variant. Patient fibroblasts showed an accumulation of mitochondrial unprocessed PNPT1 transcripts, while blood showed an increased interferon response. Our findings suggest that functional analyses of the RNA processing function of PNPase are more sensitive than testing downstream defects in oxidative phosphorylation (OXPHPOS) enzyme activities. This research extends our knowledge of the clinical and functional consequences of bi-allelic pathogenic PNPT1 variants that may guide management and further efforts into understanding the pathophysiological mechanisms for therapeutic development.