The rules and impact of nonsense-mediated mRNA decay in human cancers

Structure of the Nmd4-Upf1 complex supports conservation of the nonsense-mediated mRNA decay pathway between yeast and humans

The nonsense-mediated mRNA decay (NMD) pathway clears eukaryotic cells of mRNAs containing premature termination codons (PTCs) or normal stop codons located in specific contexts. It therefore plays an important role in gene expression regulation. The precise molecular mechanism of the NMD pathway has long been considered to differ substantially from yeast to metazoa, despite the involvement of universally conserved factors such as the central ATP-dependent RNA-helicase Upf1. Here, we describe the crystal structure of the yeast Upf1 bound to its recently identified but yet uncharacterized partner Nmd4, show that Nmd4 stimulates Upf1 ATPase activity and that this interaction contributes to the elimination of NMD substrates. We also demonstrate that a region of Nmd4 critical for the interaction with Upf1 in yeast is conserved in the metazoan SMG6 protein, another major NMD factor. We show that this conserved region is involved in the interaction of SMG6 with UPF1 and that mutations in this region affect the levels of endogenous human NMD substrates. Our results support the universal conservation of the NMD mechanism in eukaryotes.

Investigating the genetic basis of hereditary spastic paraplegia and cerebellar Ataxia in Pakistani families

BACKGROUND

Hereditary Spastic Paraplegias (HSPs) and Hereditary Cerebellar Ataxias (HCAs) are progressive neurodegenerative disorders encompassing a spectrum of neurogenetic conditions with significant overlaps of clinical features. Spastic ataxias are a group of conditions that have features of both cerebellar ataxia and spasticity, and these conditions are frequently clinically challenging to distinguish. Accurate genetic diagnosis is crucial but challenging, particularly in resource-limited settings. This study aims to investigate the genetic basis of HSPs and HCAs in Pakistani families.

METHODS

Families from Khyber Pakhtunkhwa with at least two members showing HSP or HCA phenotypes, and who had not previously been analyzed genetically, were included. Families were referred for genetic analysis by local neurologists based on the proband's clinical features and signs of a potential genetic neurodegenerative disorder. Whole Exome Sequencing (WES) and Sanger sequencing were then used to identify and validate genetic variants, and to analyze variant segregation within families to determine inheritance patterns. The mean age of onset and standard deviation were calculated to assess variability among affected individuals, and the success rate was compared with literature reports using differences in proportions and Cohen's h.

RESULTS

Pathogenic variants associated with these conditions were identified in five of eight families, segregating according to autosomal recessive inheritance. These variants included previously reported SACS c.2182 C > T, p.(Arg728*), FA2H c.159_176del, p.(Arg53_Ile58del) and SPG11 c.2146 C > T, p.(Gln716*) variants, and two previously unreported variants in SACS c.2229del, p.(Phe743Leufs*8) and ZFYVE26 c.1926_1941del, p.(Tyr643Metfs*2). Additionally, FA2H and SPG11 variants were found to have recurrent occurrences, suggesting a potential founder effect within the Pakistani population. Onset age among affected individuals ranged from 1 to 14 years (M = 6.23, SD = 3.96). The diagnostic success rate was 62.5%, with moderate effect sizes compared to previous studies.

CONCLUSIONS

The findings of this study expand the genotypic and phenotypic spectrum of HSPs and HCAs in Pakistan and emphasize the importance of utilizing exome/genome sequencing for accurate diagnosis or support accurate differential diagnosis. This approach can improve genetic counseling and clinical management, addressing the challenges of diagnosing neurodegenerative disorders in resource-limited settings.

Saturation mutagenesis-reinforced functional assays for disease-related genes

Interpretation of disease-causing genetic variants remains a challenge in human genetics. Current costs and complexity of deep mutational scanning methods are obstacles for achieving genome-wide resolution of variants in disease-related genes. Our framework, saturation mutagenesis-reinforced functional assays (SMuRF), offers simple and cost-effective saturation mutagenesis paired with streamlined functional assays to enhance the interpretation of unresolved variants. Applying SMuRF to neuromuscular disease genes FKRP and LARGE1, we generated functional scores for all possible coding single-nucleotide variants, which aid in resolving clinically reported variants of uncertain significance. SMuRF also demonstrates utility in predicting disease severity, resolving critical structural regions, and providing training datasets for the development of computational predictors. Overall, our approach enables variant-to-function insights for disease genes in a cost-effective manner that can be broadly implemented by standard research laboratories.

NMDtxDB: data-driven identification and annotation of human NMD target transcripts

The nonsense-mediated RNA decay (NMD) pathway is a crucial mechanism of mRNA quality control. Current annotations of NMD substrate RNAs are rarely data-driven, but use generally established rules. We present a data set with four cell lines and combinations for SMG5, SMG6, and SMG7 knockdowns or SMG7 knockout. Based on this data set, we implemented a workflow that combines Nanopore and Illumina sequencing to assemble a transcriptome, which is enriched for NMD target transcripts. Moreover, we use coding sequence information (CDS) from Ensembl, Gencode consensus Ribo-seq ORFs, and OpenProt to enhance the CDS annotation of novel transcript isoforms. In summary, 302,889 transcripts were obtained from the transcriptome assembly process, out of which 24% are absent from Ensembl database annotations, 48,213 contain a premature stop codon, and 6433 are significantly upregulated in three or more comparisons of NMD active versus deficient cell lines. We present an in-depth view of these results through the NMDtxDB database, which is available at https://shiny.dieterichlab.org/app/NMDtxDB, and supports the study of NMD-sensitive transcripts. We open sourced our implementation of the respective web-application and analysis workflow at https://github.com/dieterich-lab/NMDtxDB and https://github.com/dieterich-lab/nmd-wf.

Novel autosomal recessive SINO syndrome-associated KIDINS220 variants provide insight into the genotype-phenotype correlation

Background

KIDINS220 encodes a transmembrane scaffold protein, kinase D-interacting substrate of 220 kDa, that regulates neurotrophin signaling. Variants in KIDINS220 have been linked to spastic paraplegia, intellectual disability, nystagmus, and obesity (SINO) syndrome or prenatal fatal cerebral ventriculomegaly and arthrogryposis (VENARG). This study aimed to investigate the genotype-phenotype correlation of pathogenic KIDINS220 variants.

Methods

We performed whole-exome sequencing on a patient with SINO syndrome and epilepsy. Identified pathogenic variants were confirmed using Sanger sequencing and evaluated with in silico tools. A comprehensive literature review was conducted to analyze the genetic and phenotypic data of both the newly diagnosed patient and previously reported cases with KIDINS220 variants.

Results

We identified novel compound heterozygous variants in KIDINS220, c.1556C > T (p.Thr519Met) and c.2374C > T (p.Arg792*), in the patient. Our analysis revealed that biallelic loss-of-function variants in KIDINS220 are associated with VENARG or autosomal recessive SINO (AR-SINO), whereas carboxy-terminal truncated variants that escape nonsense-mediated mRNA decay and lack amino acid residues 1507-1529 are linked to autosomal dominant SINO (AD-SINO). Patients with AR-SINO exhibit more severe clinical features compared to those with AD-SINO.

Conclusions

Our study expands the spectrum of KIDINS220 variants associated with AR-SINO and provides a valuable genotype-phenotype correlation for pathogenic KIDINS220 variants.

Nonsense-Mediated mRNA Decay in Human Health and Diseases: Current Understanding, Regulatory Mechanisms and Future Perspectives

Nonsense-mediated mRNA decay (NMD) is a surveillance mechanism that is conserved across all eukaryotes ensuring the quality of transcripts by targeting messenger RNA (mRNA) harbouring premature stop codons. It regulates the gene expression by targeting aberrant mRNA carrying pre-termination codons (PTCs) and eliminates C-terminal truncated proteins. NMD distinguishes aberrant and non-aberrant transcript by looking after long 3' UTRs and exon-junction complex (EJC) downstream of stop codon that indicate the presence of PTC. Therefore, NMD modulates cellular surveillance and eliminates the truncated proteins but if the PTC escapes the surveillance pathway it can lead to potential negative phenotype resulting in genetic diseases. The alternative splicing also contributes in formation of NMD-sensitive isoforms by introducing PTC. NMD plays a complex role in cancer, it can either aggravate or downregulates the tumour. Some tumours agitate NMD to deteriorate mRNAs encoding tumour suppressor proteins, stress response proteins and neoantigens. In other case, tumours suppress the NMD to encourage the expression of oncoproteins for tumour growth and survival. This mechanism augmented in the development of new therapeutics by PTC read-through mechanism and personalized medicine. Detailed studies on NMD surveillance will possibly lead towards development of strategies for improving human health aligning with United Nations sustainable development goals (SDG 3: Good health and well-being). The potential therapeutic applications of NMD pose a challenge in terms of safe and effective modulation. Understanding the complexities of NMD regulation and its interaction with other cellular processes can lead to the development of new interventions for various diseases.

A new nonsense pathogenic variant in exon 1 of PHOX2B leads to the diagnosis of congenital central hypoventilation syndrome with intra-familial variability

Congenital central hypoventilation syndrome (CCHS) is a rare genetic disorder of the autonomic nervous system resulting in decreased brain sensitivity to hypercapnia and hypoxia characterized by a genetic abnormality in the pair-like homeobox 2B (PHOX2B) gene. Most patients have a heterozygous expansion of the polyalanine repeat in exon 3 (PARM), while 10 % of patients have non-PARM (NPARM) mutations that can span the entire gene. The majority of pathogenic variants are de novo, but variants with incomplete penetrance can be identified in the heterozygous state. In the present report, CCHS was diagnosed in a symptomatic 3-month-old infant with neonatal respiratory distress. Genetic analysis revealed a new mutation in exon 1 of the PHOX2B gene - p.Ser28* (c.83C>G) - which was further identified in two family members, one minimally symptomatic and one asymptomatic. The identification of this new mutation supports the importance of sequencing the entire gene even when the classic PARM mutation is not found and highlights the phenotypic variability of CCHS.

RBM8A, a new target of TEAD4, promotes breast cancer progression by regulating IGF1R and IRS-2

BACKGROUND

Breast cancer (BC) is the most common malignant tumor in women worldwide, and further elucidation of the molecular mechanisms involved in BC pathogenesis is essential to improve the prognosis of BC patients. RNA Binding Motif Protein 8 A (RBM8A), with high affinity to a myriad of RNA transcripts, has been shown to play a crucial role in genesis and progression of multiple cancers. We attempted to explore its functional significance and molecular mechanisms in BC.

METHODS

Bioinformatics analysis was performed on publicly available BC datasets. qRT-PCR was used to determine the expression of RBM8A in BC tissues. MTT assay, clone formation assay and flow cytometry were employed to examine BC cell proliferation and apoptosis in vitro. RNA immunoprecipitation (RIP) and RIP-seq were used to investigate the binding of RBM8A/EIF4A3 to the mRNA of IGF1R/IRS-2. RBM8A and EIF4A3 interactions were determined by co-immunoprecipitation (Co-IP) and immunofluorescence. Chromatin immunoprecipitation (Ch-IP) and dual-luciferase reporter assay were carried out to investigate the transcriptional regulation of RBM8A by TEAD4. Xenograft model was used to explore the effects of RBM8A and TEAD4 on BC cell growth in vivo.

RESULTS

In this study, we showed that RBM8A is abnormally highly expressed in BC and knockdown of RBM8A inhibits BC cell proliferation and induces apoptosis in vitro. EIF4A3, which phenocopy RBM8A in BC, forms a complex with RBM8A in BC. Moreover, EIF4A3 and RBM8A complex regulate the expression of IGF1R and IRS-2 to activate the PI3K/AKT signaling pathway, thereby promoting BC progression. In addition, we identified TEAD4 as a transcriptional activator of RBM8A by Ch-IP, dual luciferase reporter gene and a series of functional rescue assays. Furthermore, we demonstrated the in vivo pro-carcinogenic effects of TEAD4 and RBM8A by xenograft tumor experiments in nude mice.

CONCLUSION

Collectively, these findings suggest that TEAD4 novel transcriptional target RBM8A interacts with EIF4A3 to increase IGF1R and IRS-2 expression and activate PI3K/AKT signaling pathway, thereby further promoting the malignant phenotype of BC cells.

Mutations in the U2 snRNA gene RNU2-2P cause a severe neurodevelopmental disorder with prominent epilepsy

The major spliceosome comprises the five snRNAs U1, U2, U4, U5 and U6. We recently showed that mutations in RNU4- 2, which encodes U4 snRNA, cause one of the most prevalent monogenic neurodevelopmental disorders. Here, we report that recurrent germline mutations in RNU2-2P , a 191bp gene encoding U2 snRNA, are responsible for a related disorder. By genetic association, we implicated recurrent de novo single nucleotide mutations at nucleotide positions 4 and 35 of RNU2-2P among nine cases. We replicated this finding in six additional cases, bringing the total to 15. The disorder is characterized by intellectual disability, neurodevelopmental delay, autistic behavior, microcephaly, hypotonia, epilepsy and hyperventilation. All cases display a severe and complex seizure phenotype. Our findings cement the role of major spliceosomal snRNAs in the etiologies of neurodevelopmental disorders.

Knockout, Knockdown, and the Schrödinger Paradox: Genetic Immunity to Phenotypic Recapitulation in Zebrafish

Previous research has highlighted significant phenotypic discrepancies between knockout and knockdown approaches in zebrafish, raising concerns about the reliability of these methods. However, our study suggests that these differences are not as pronounced as was once believed. By carefully examining the roles of maternal and zygotic gene contributions, we demonstrate that these factors significantly influence phenotypic outcomes, often accounting for the observed discrepancies. Our findings emphasize that morpholinos, despite their potential off-target effects, can be effective tools when used with rigorous controls. We introduce the concept of graded maternal contribution, which explains how the uneven distribution of maternal mRNA and proteins during gametogenesis impacts phenotypic variability. Our research categorizes genes into three types-susceptible, immune, and "Schrödinger" (conditional)-based on their phenotypic expression and interaction with genetic compensation mechanisms. This distinction provides new insights into the paradoxical outcomes observed in genetic studies. Ultimately, our work underscores the importance of considering both maternal and zygotic contributions, alongside rigorous experimental controls, to accurately interpret gene function and the mechanisms underlying disease. This study advocates for the continued use of morpholinos in conjunction with advanced genetic tools like CRISPR/Cas9, stressing the need for a meticulous experimental design to optimize the utility of zebrafish in genetic research and therapeutic development.

SERPINA11 related novel serpinopathy - A perinatal lethal disorder

SERPINA11 is a hitherto poorly characterised gene belonging to Clade A of the SERPIN superfamily, with unknown expression pattern and functional significance. We report a perinatal lethal phenotype in two foetuses from the same family associated with a biallelic loss of function variant in SERPINA11, and provide functional evidence to support its candidature as a Mendelian disorder. The SERPINA11 variant-associated foetal phenotype is characterised by gross and histopathological features of extracellular matrix disruption. Western blot and immunofluorescence analyses revealed SERPINA11 expression in multiple mouse tissues, with pronounced expression in the bronchiolar epithelium. We observed a significant decrease in SERPINA11 immunofluorescence in the affected foetal lung compared with a healthy gestation-matched foetus. Protein expression data from HEK293T cell lines following site-directed mutagenesis support the loss of function nature of the variant. Transcriptome analysis from the affected foetal liver indicated the possibility of reduced SERPINA11 transcript abundance. This novel serpinopathy appears to be a consequence of the loss of inhibition of serine proteases involved in extracellular matrix remodelling, revealing SERPINA11 as a protease inhibitor critical for embryonic development.

Prognostic genome and transcriptome signatures in colorectal cancers

Colorectal cancer is caused by a sequence of somatic genomic alterations affecting driver genes in core cancer pathways1. Here, to understand the functional and prognostic impact of cancer-causing somatic mutations, we analysed the whole genomes and transcriptomes of 1,063 primary colorectal cancers in a population-based cohort with long-term follow-up. From the 96 mutated driver genes, 9 were not previously implicated in colorectal cancer and 24 had not been linked to any cancer. Two distinct patterns of pathway co-mutations were observed, timing analyses identified nine early and three late driver gene mutations, and several signatures of colorectal-cancer-specific mutational processes were identified. Mutations in WNT, EGFR and TGFβ pathway genes, the mitochondrial CYB gene and 3 regulatory elements along with 21 copy-number variations and the COSMIC SBS44 signature correlated with survival. Gene expression classification yielded five prognostic subtypes with distinct molecular features, in part explained by underlying genomic alterations. Microsatellite-instable tumours divided into two classes with different levels of hypoxia and infiltration of immune and stromal cells. To our knowledge, this study constitutes the largest integrated genome and transcriptome analysis of colorectal cancer, and interlinks mutations, gene expression and patient outcomes. The identification of prognostic mutations and expression subtypes can guide future efforts to individualize colorectal cancer therapy.

Joint, multifaceted genomic analysis enables diagnosis of diverse, ultra-rare monogenic presentations

Genomics for rare disease diagnosis has advanced at a rapid pace due to our ability to perform "N-of-1" analyses on individual patients with ultra-rare diseases. The increasing sizes of ultra-rare disease cohorts internationally newly enables cohort-wide analyses for new discoveries, but well-calibrated statistical genetics approaches for jointly analyzing these patients are still under development.1,2 The Undiagnosed Diseases Network (UDN) brings multiple clinical, research and experimental centers under the same umbrella across the United States to facilitate and scale N-of-1 analyses. Here, we present the first joint analysis of whole genome sequencing data of UDN patients across the network. We introduce new, well-calibrated statistical methods for prioritizing disease genes with de novo recurrence and compound heterozygosity. We also detect pathways enriched with candidate and known diagnostic genes. Our computational analysis, coupled with a systematic clinical review, recapitulated known diagnoses and revealed new disease associations. We further release a software package, RaMeDiES, enabling automated cross-analysis of deidentified sequenced cohorts for new diagnostic and research discoveries. Gene-level findings and variant-level information across the cohort are available in a public-facing browser (https://dbmi-bgm.github.io/udn-browser/). These results show that N-of-1 efforts should be supplemented by a joint genomic analysis across cohorts.

RNA variant assessment using transactivation and transdifferentiation

Understanding the impact of splicing and nonsense variants on RNA is crucial for the resolution of variant classification as well as their suitability for precision medicine interventions. This is primarily enabled through RNA studies involving transcriptomics followed by targeted assays using RNA isolated from clinically accessible tissues (CATs) such as blood or skin of affected individuals. Insufficient disease gene expression in CATs does however pose a major barrier to RNA based investigations, which we show is relevant to 1,436 Mendelian disease genes. We term these "silent" Mendelian genes (SMGs), the largest portion (36%) of which are associated with neurological disorders. We developed two approaches to induce SMG expression in human dermal fibroblasts (HDFs) to overcome this limitation, including CRISPR-activation-based gene transactivation and fibroblast-to-neuron transdifferentiation. Initial transactivation screens involving 40 SMGs stimulated our development of a highly multiplexed transactivation system culminating in the 6- to 90,000-fold induction of expression of 20/20 (100%) SMGs tested in HDFs. Transdifferentiation of HDFs directly to neurons led to expression of 193/516 (37.4%) of SMGs implicated in neurological disease. The magnitude and isoform diversity of SMG expression following either transactivation or transdifferentiation was comparable to clinically relevant tissues. We apply transdifferentiation and/or gene transactivation combined with short- and long-read RNA sequencing to investigate the impact that variants in USH2A, SCN1A, DMD, and PAK3 have on RNA using HDFs derived from affected individuals. Transactivation and transdifferentiation represent rapid, scalable functional genomic solutions to investigate variants impacting SMGs in the patient cell and genomic context.

Ornaments for efficient allele-specific expression estimation with bias correction

Allele-specific expression plays a crucial role in unraveling various biological mechanisms, including genomic imprinting and gene expression controlled by cis-regulatory variants. However, existing methods for quantification from RNA-sequencing (RNA-seq) reads do not adequately and efficiently remove various allele-specific read mapping biases, such as reference bias arising from reads containing the alternative allele that do not map to the reference transcriptome or ambiguous mapping bias caused by reads containing the reference allele that map differently from reads containing the alternative allele. We present Ornaments, a computational tool for rapid and accurate estimation of allele-specific transcript expression at unphased heterozygous loci from RNA-seq reads while correcting for allele-specific read mapping biases. Ornaments removes reference bias by mapping reads to a personalized transcriptome and ambiguous mapping bias by probabilistically assigning reads to multiple transcripts and variant loci they map to. Ornaments is a lightweight extension of kallisto, a popular tool for fast RNA-seq quantification, that improves the efficiency and accuracy of WASP, a popular tool for bias correction in allele-specific read mapping. In experiments with simulated and human lymphoblastoid cell-line RNA-seq reads with the genomes of the 1000 Genomes Project, we demonstrate that Ornaments improves the accuracy of WASP and kallisto, is nearly as efficient as kallisto, and is an order of magnitude faster than WASP per sample, with the additional cost of constructing a personalized index for multiple samples. Additionally, we show that Ornaments finds imprinted transcripts with higher sensitivity than WASP, which detects imprinted signals only at gene level.

SMG5 Inhibition Restrains Hepatocellular Carcinoma Growth and Enhances Sorafenib Sensitivity

Hepatocellular carcinoma (HCC) has a pathogenesis that remains elusive with restricted therapeutic strategies and efficacy. This study aimed to investigate the role of SMG5, a crucial component in nonsense-mediated mRNA decay (NMD) that degrades mRNA containing a premature termination codon, in HCC pathogenesis and therapeutic resistance. We demonstrated an elevated expression of SMG5 in HCC and scrutinized its potential as a therapeutic target. Our findings revealed that SMG5 knockdown not only inhibited the migration, invasion, and proliferation of HCC cells but also influenced sorafenib resistance. Differential gene expression analysis between the control and SMG5 knockdown groups showed an upregulation of methionine adenosyltransferase 1A in the latter. High expression of methionine adenosyltransferase 1A, a catalyst for S-adenosylmethionine (SAM) production, as suggested by The Cancer Genome Atlas data, was indicative of a better prognosis for HCC. Further, an ELISA showed a higher concentration of SAM in SMG5 knockdown cell supernatants. Furthermore, we found that exogenous SAM supplementation enhanced the sensitivity of HCC cells to sorafenib alongside changes in the expression of Bax and Bcl-2, apoptosis-related proteins. Our findings underscore the important role of SMG5 in HCC development and its involvement in sorafenib resistance, highlighting it as a potential target for HCC treatment.

Genebe.net: Implementation and validation of an automatic ACMG variant pathogenicity criteria assignment

We present GeneBe, an online platform streamlining the automated application of American College of Medical Genetics and Genomics (ACMG), Association for Molecular Pathology (AMP), and the College of American Pathologists (CAP) criteria for assessment of pathogenicity of genetic variants. GeneBe utilizes automated algorithms that evaluate 17 criteria from 28, closely aligning with current guidelines and leveraging data from diverse sources, including ClinVar. The user-friendly web interface enables manual refinement of assignments for specific criteria based on site-collected data. Our algorithm demonstrates a high correlation (r = 0.90) of assigned pathogenicity scores compared to expert assessments from the ClinGen Evidence Repository and substantial concordance with ClinVar verdict assignments (κ = 0.69). Comparative analysis with other published tools reveals that GeneBe performs similarly to VarSome while being superior over TAPES and InterVar. In contrast to some other tools, GeneBe's web implementation is tracker-free and third-party request-free, safeguarding user privacy. Additionally, GeneBe offers an Application Programming Interface (API) for enhanced flexibility and integration into existing workflows and is provided free of charge for research purposes. GeneBe is available at https://genebe.net.

Exploring the Impact of Genetics in a Large Cohort of Moebius Patients by Trio Whole Exome Sequencing

Moebius syndrome (MBS) is a rare congenital disorder characterized by non-progressive facial palsy and ocular abduction paralysis. Most cases are sporadic, but also rare familial cases with autosomal dominant transmission and incomplete penetrance/variable expressivity have been described. The genetic etiology of MBS is still unclear: de novo pathogenic variants in REV3L and PLXND1 are reported in only a minority of cases, suggesting the involvement of additional causative genes. With the aim to uncover the molecular causative defect and identify a potential genetic basis of this condition, we performed trio-WES on a cohort of 37 MBS and MBS-like patients. No de novo variants emerged in REV3L and PLXND1. We then proceeded with a cohort analysis to identify possible common causative genes among all patients and a trio-based analysis using an in silico panel of candidate genes. However, identified variants emerging from both approaches were considered unlikely to be causative of MBS, mainly due to the lack of clinical overlap. In conclusion, despite this large cohort, WES failed to identify mutations possibly associated with MBS, further supporting the heterogeneity of this syndrome, and suggesting the need for integrated omics approaches to identify the molecular causes underlying MBS development.

Copy number losses of oncogenes and gains of tumor suppressor genes generate common driver mutations

Cancer driver genes can undergo positive selection for various types of genetic alterations, including gain-of-function or loss-of-function mutations and copy number alterations (CNA). We investigated the landscape of different types of alterations affecting driver genes in 17,644 cancer exomes and genomes. We find that oncogenes may simultaneously exhibit signatures of positive selection and also negative selection in different gene segments, suggesting a method to identify additional tumor types where an oncogene is a driver or a vulnerability. Next, we characterize the landscape of CNA-dependent selection effects, revealing a general trend of increased positive selection on oncogene mutations not only upon CNA gains but also upon CNA deletions. Similarly, we observe a positive interaction between mutations and CNA gains in tumor suppressor genes. Thus, two-hit events involving point mutations and CNA are universally observed regardless of the type of CNA and may signal new therapeutic opportunities. An analysis with focus on the somatic CNA two-hit events can help identify additional driver genes relevant to a tumor type. By a global inference of point mutation and CNA selection signatures and interactions thereof across genes and tissues, we identify 9 evolutionary archetypes of driver genes, representing different mechanisms of (in)activation by genetic alterations.

ZSCAN10 deficiency causes a neurodevelopmental disorder with characteristic oto-facial malformations

Neurodevelopmental disorders are major indications for genetic referral and have been linked to more than 1500 loci including genes encoding transcriptional regulators. The dysfunction of transcription factors often results in characteristic syndromic presentations; however, at least half of these patients lack a genetic diagnosis. The implementation of machine learning approaches has the potential to aid in the identification of new disease genes and delineate associated phenotypes. Next generation sequencing was performed in seven affected individuals with neurodevelopmental delay and dysmorphic features. Clinical characterization included reanalysis of available neuroimaging datasets and 2D portrait image analysis with GestaltMatcher. The functional consequences of ZSCAN10 loss were modelled in mouse embryonic stem cells (mESCs), including a knockout and a representative ZSCAN10 protein truncating variant. These models were characterized by gene expression and western blot analyses, chromatin immunoprecipitation and quantitative PCR (ChIP-qPCR) and immunofluorescence staining. Zscan10 knockout mouse embryos were generated and phenotyped. We prioritized bi-allelic ZSCAN10 loss-of-function variants in seven affected individuals from five unrelated families as the underlying molecular cause. RNA-sequencing analyses in Zscan10-/- mESCs indicated dysregulation of genes related to stem cell pluripotency. In addition, we established in mESCs the loss-of-function mechanism for a representative human ZSCAN10 protein truncating variant by showing alteration of its expression levels and subcellular localization, interfering with its binding to DNA enhancer targets. Deep phenotyping revealed global developmental delay, facial asymmetry and malformations of the outer ear as consistent clinical features. Cerebral MRI showed dysplasia of the semicircular canals as an anatomical correlate of sensorineural hearing loss. Facial asymmetry was confirmed as a clinical feature by GestaltMatcher and was recapitulated in the Zscan10 mouse model along with inner and outer ear malformations. Our findings provide evidence of a novel syndromic neurodevelopmental disorder caused by bi-allelic loss-of-function variants in ZSCAN10.

Cell-Specific Single Viral Vector CRISPR/Cas9 Editing and Genetically Encoded Tool Delivery in the Central and Peripheral Nervous Systems

CRISPR/Cas9 gene editing represents an exciting avenue to study genes of unknown function and can be combined with genetically encoded tools such as fluorescent proteins, channelrhodopsins, DREADDs, and various biosensors to more deeply probe the function of these genes in different cell types. However, current strategies to also manipulate or visualize edited cells are challenging due to the large size of Cas9 proteins and the limited packaging capacity of adeno-associated viruses (AAVs). To overcome these constraints, we developed an alternative gene editing strategy using a single AAV vector and mouse lines that express Cre-dependent Cas9 to achieve efficient cell-type specific editing across the nervous system. Expressing Cre-dependent Cas9 from a genomic locus affords space to package guide RNAs for gene editing together with Cre-dependent, genetically encoded tools to manipulate, map, or monitor neurons using a single virus. We validated this strategy with three common tools in neuroscience: ChRonos, a channelrhodopsin, for studying synaptic transmission using optogenetics, GCaMP8f for recording Ca2+ transients using photometry, and mCherry for tracing axonal projections. We tested these tools in multiple brain regions and cell types, including GABAergic neurons in the nucleus accumbens, glutamatergic neurons projecting from the ventral pallidum to the lateral habenula, dopaminergic neurons in the ventral tegmental area, and proprioceptive neurons in the periphery. This flexible approach could help identify and test the function of novel genes affecting synaptic transmission, circuit activity, or morphology with a single viral injection.

Mechanisms and Delivery of tRNA Therapeutics

Transfer ribonucleic acid (tRNA) therapeutics will provide personalized and mutation specific medicines to treat human genetic diseases for which no cures currently exist. The tRNAs are a family of adaptor molecules that interpret the nucleic acid sequences in our genes into the amino acid sequences of proteins that dictate cell function. Humans encode more than 600 tRNA genes. Interestingly, even healthy individuals contain some mutant tRNAs that make mistakes. Missense suppressor tRNAs insert the wrong amino acid in proteins, and nonsense suppressor tRNAs read through premature stop signals to generate full length proteins. Mutations that underlie many human diseases, including neurodegenerative diseases, cancers, and diverse rare genetic disorders, result from missense or nonsense mutations. Thus, specific tRNA variants can be strategically deployed as therapeutic agents to correct genetic defects. We review the mechanisms of tRNA therapeutic activity, the nature of the therapeutic window for nonsense and missense suppression as well as wild-type tRNA supplementation. We discuss the challenges and promises of delivering tRNAs as synthetic RNAs or as gene therapies. Together, tRNA medicines will provide novel treatments for common and rare genetic diseases in humans.

Deep Mutational Scanning in Disease-related Genes with Saturation Mutagenesis-Reinforced Functional Assays (SMuRF)

Interpretation of disease-causing genetic variants remains a challenge in human genetics. Current costs and complexity of deep mutational scanning methods hamper crowd-sourcing approaches toward genome-wide resolution of variants in disease-related genes. Our framework, Saturation Mutagenesis-Reinforced Functional assays (SMuRF), addresses these issues by offering simple and cost-effective saturation mutagenesis, as well as streamlining functional assays to enhance the interpretation of unresolved variants. Applying SMuRF to neuromuscular disease genes FKRP and LARGE1, we generated functional scores for all possible coding single nucleotide variants, which aid in resolving clinically reported variants of uncertain significance. SMuRF also demonstrates utility in predicting disease severity, resolving critical structural regions, and providing training datasets for the development of computational predictors. Our approach opens new directions for enabling variant-to-function insights for disease genes in a manner that is broadly useful for crowd-sourcing implementation across standard research laboratories.

Identification of nonsense-mediated decay inhibitors that alter the tumor immune landscape

Despite exciting developments in cancer immunotherapy, its broad application is limited by the paucity of targetable antigens on the tumor cell surface. As an intrinsic cellular pathway, nonsense-mediated decay (NMD) conceals neoantigens through the destruction of the RNA products from genes harboring truncating mutations. We developed and conducted a high throughput screen, based on the ratiometric analysis of transcripts, to identify critical mediators of NMD. This screen implicated disruption of kinase SMG1's phosphorylation of UPF1 as a potential disruptor of NMD. This led us to design a novel SMG1 inhibitor, KVS0001, that elevates the expression of transcripts and proteins resulting from truncating mutations in vivo and in vitro . Most importantly, KVS0001 concomitantly increased the presentation of immune-targetable HLA class I-associated peptides from NMD-downregulated proteins on the surface of cancer cells. KVS0001 provides new opportunities for studying NMD and the diseases in which NMD plays a role, including cancer and inherited diseases.

One Sentence Summary

Disruption of the nonsense-mediated decay pathway with a newly developed SMG1 inhibitor with in-vivo activity increases the expression of T-cell targetable cancer neoantigens resulting from truncating mutations.

A novel IRF6 gene mutation impacting the regulation of TGFβ2-AS1 in the TGFβ pathway: A mechanism in the development of Van der Woude syndrome

Several mutations in the IRF6 gene have been identified as a causative link to VWS. In this investigation, whole-exome sequencing (WES) and Sanger sequencing of a three-generation pedigree with an autosomal-dominant inheritance pattern affected by VWS identified a unique stop-gain mutation-c.748C>T:p.R250X-in the IRF6 gene that co-segregated exclusively with the disease phenotype. Immunofluorescence analysis revealed that the IRF6-p.R250X mutation predominantly shifted its localization from the nucleus to the cytoplasm. WES and protein interaction analyses were conducted to understand this mutation's role in the pathogenesis of VWS. Using LC-MS/MS, we found that this mutation led to a reduction in the binding of IRF6 to histone modification-associated proteins (NAA10, SNRPN, NAP1L1). Furthermore, RNA-seq results show that the mutation resulted in a downregulation of TGFβ2-AS1 expression. The findings highlight the mutation's influence on TGFβ2-AS1 and its subsequent effects on the phosphorylation of SMAD2/3, which are critical in maxillofacial development, particularly the palate. These insights contribute to a deeper understanding of VWS's molecular underpinnings and might inform future therapeutic strategies.

Mutations causing premature termination codons discriminate and generate cellular and clinical variability in HHT

ABSTRACT

For monogenic diseases caused by pathogenic loss-of-function DNA variants, attention focuses on dysregulated gene-specific pathways, usually considering molecular subtypes together within causal genes. To better understand phenotypic variability in hereditary hemorrhagic telangiectasia (HHT), we subcategorized pathogenic DNA variants in ENG/endoglin, ACVRL1/ALK1, and SMAD4 if they generated premature termination codons (PTCs) subject to nonsense-mediated decay. In 3 patient cohorts, a PTC-based classification system explained some previously puzzling hemorrhage variability. In blood outgrowth endothelial cells (BOECs) derived from patients with ACVRL1+/PTC, ENG+/PTC, and SMAD4+/PTC genotypes, PTC-containing RNA transcripts persisted at low levels (8%-23% expected, varying between replicate cultures); genes differentially expressed to Bonferroni P < .05 in HHT+/PTC BOECs clustered significantly only to generic protein terms (isopeptide-bond/ubiquitin-like conjugation) and pulse-chase experiments detected subtle protein maturation differences but no evidence for PTC-truncated protein. BOECs displaying highest PTC persistence were discriminated in unsupervised hierarchical clustering of near-invariant housekeeper genes, with patterns compatible with higher cellular stress in BOECs with >11% PTC persistence. To test directionality, we used a HeLa reporter system to detect induction of activating transcription factor 4 (ATF4), which controls expression of stress-adaptive genes, and showed that ENG Q436X but not ENG R93X directly induced ATF4. AlphaFold accurately modeled relevant ENG domains, with AlphaMissense suggesting that readthrough substitutions would be benign for ENG R93X and other less rare ENG nonsense variants but more damaging for Q436X. We conclude that PTCs should be distinguished from other loss-of-function variants, PTC transcript levels increase in stressed cells, and readthrough proteins and mechanisms provide promising research avenues.

Identification of a Novel Frameshift Variant in MYF5 Leading to External Ophthalmoplegia with Rib and Vertebral Anomalies

Myogenic transcription factors with a basic helix-loop-helix (bHLH) such as MYOD, myogenin, MRF4, and MYF5 contribute to muscle differentiation and regulation. The MYF5 gene located on chromosome 12 encodes for myogenic factor 5 (MYF5), which has a role in skeletal and extraocular muscle development and rib formation. Variants in MYF5 were found to cause external ophthalmoplegia with rib and vertebral anomalies (EORVA), a rare recessive condition. To date, three homozygous variants in MYF5 have been reported to cause EORVA in six members of four unrelated families. Here, we present a novel homozygous MYF5 frameshift variant, c.596dupA p. (Asn199Lysfs*49), causing premature protein termination and presenting with external ophthalmoplegia, ptosis, and scoliosis in three siblings from a consanguineous family of Pakistani origin. With four MYF5 variants now discovered, genetic testing and paediatric assessment for extra-ocular features should be considered in all cases of congenital ophthalmoplegia.

Genomic analysis of a spontaneous unifoliate mutant reveals gene candidates associated with compound leaf development in Vigna unguiculata [L] Walp

Molecular mechanisms which underpin compound leaf development in some legumes have been reported, but there is no previous study on the molecular genetic control of compound leaf formation in Vigna unguiculata (cowpea), an important dryland legume of African origin. In most studied species with compound leaves, class 1 KNOTTED-LIKE HOMEOBOX genes expressed in developing leaf primordia sustain morphogenetic activity, allowing leaf dissection and the development of leaflets. Other genes, such as, SINGLE LEAFLET1 in Medicago truncatula and Trifoliate in Solanum lycopersicum, are also implicated in regulating compound leaf patterning. To set the pace for an in-depth understanding of the genetics of compound leaf development in cowpea, we applied RNA-seq and whole genome shotgun sequence datasets of a spontaneous cowpea unifoliate mutant and its trifoliate wild-type cultivar to conduct comparative reference-based gene expression, de novo genome-wide isoform switch, and genome variant analyses between the two genotypes. Our results suggest that genomic variants upstream of LATE ELONGATED HYPOCOTYL and down-stream of REVEILLE4, BRASSINOSTERIOD INSENSITIVE1 and LATERAL ORGAN BOUNDARIES result in down-regulation of key components of cowpea circadian rhythm central oscillator and brassinosteroid signaling, resulting in unifoliate leaves and brassinosteroid-deficient-like phenotypes. We have stated hypotheses that will guide follow-up studies expected to provide more insights.

A Novel Four-Gene Signature Based on Nonsense-Mediated RNA Decay for Predicting Prognosis in Hepatocellular Carcinoma: Bioinformatics Analysis and Functional Validation

Purpose

Nonsense-mediated RNA decay (NMD), a surveillance pathway for selective degradation of aberrant mRNAs, is associated with cancer progression. Its potential as a predictor for aggressive hepatocellular carcinoma (HCC) is unclear. Here, we present an innovative NMD risk model for predicting HCC prognosis.

Methods

The Cancer Genome Atlas (TCGA) data of 374 liver HCC (LIHC) and 50 normal liver samples were extracted. A risk model based on NMD-related genes was developed through least absolute shrinkage and selection operator Cox (LASSO-Cox) regression of the LIHC-TCGA data. Prognostic validation was done using GSE54236, GSE116174, and GSE76427 data. Univariate and multivariate Cox regression analyses were conducted to assess the prognostic value of the model. We also constructed nomograms for survival prediction. Tumor immune infiltration was evaluated using the CIBERSORT algorithm, and the tumor cell phenotype was assessed. Finally, mouse experiments verified UPF3B knockdown effects on HCC tumor characteristics.

Results

We developed a risk model based on four NMD-related genes (PABPC1, RPL8, SMG5, and UPF3B) and validated it using GSE54236, GSE116174, and GSE76427 data. The model effectively distinguished high- and low-risk groups corresponding to unfavorable and favorable HCC outcomes. Its prognostic prediction accuracy was confirmed through time-dependent ROC analysis, and clinical-use nomograms with calibration curves were developed. Single-cell RNA sequencing results indicated significantly higher expression of SMG5 and UPF3B in tumor cells. Knockdown of SMG5 and UPF3B inhibited HCC cell proliferation, invasion, and migration, while affecting cell-cycle progression and apoptosis. In vivo, UPF3B knockdown delayed tumor growth and increased immune cell infiltration.

Conclusion

Our NMD-related gene-based risk model can help identify therapeutic targets and biomarkers for HCC. Additionally, it assists clinicians in predicting the prognosis of HCC patients.

Characterization of NEB pathogenic variants in patients reveals novel nemaline myopathy disease mechanisms and omecamtiv mecarbil force effects

Nebulin, a critical protein of the skeletal muscle thin filament, plays important roles in physiological processes such as regulating thin filament length (TFL), cross-bridge cycling, and myofibril alignment. Pathogenic variants in the nebulin gene (NEB) cause NEB-based nemaline myopathy (NEM2), a genetically heterogeneous disorder characterized by hypotonia and muscle weakness, currently lacking curative therapies. In this study, we examined a cohort of ten NEM2 patients, each with unique pathogenic variants, aiming to understand their impact on mRNA, protein, and functional levels. Results show that pathogenic truncation variants affect NEB mRNA stability and lead to nonsense-mediated decay of the mutated transcript. Moreover, a high incidence of cryptic splice site activation was found in patients with pathogenic splicing variants that are expected to disrupt the actin-binding sites of nebulin. Determination of protein levels revealed patients with either relatively normal or markedly reduced nebulin. We observed a positive relation between the reduction in nebulin and a reduction in TFL, or reduction in tension (both maximal and submaximal tension). Interestingly, our study revealed a pathogenic duplication variant in nebulin that resulted in a four-copy gain in the triplicate region of NEB and a much larger nebulin protein and longer TFL. Additionally, we investigated the effect of Omecamtiv mecarbil (OM), a small-molecule activator of cardiac myosin, on force production of type 1 muscle fibers of NEM2 patients. OM treatment substantially increased submaximal tension across all NEM2 patients ranging from 87 to 318%, with the largest effects in patients with the lowest level of nebulin. In summary, this study indicates that post-transcriptional or post-translational mechanisms regulate nebulin expression. Moreover, we propose that the pathomechanism of NEM2 involves not only shortened but also elongated thin filaments, along with the disruption of actin-binding sites resulting from pathogenic splicing variants. Significantly, our findings highlight the potential of OM treatment to improve skeletal muscle function in NEM2 patients, especially those with large reductions in nebulin levels.

Genotype-phenotype associations in 1018 individuals with SCN1A-related epilepsies

OBJECTIVE

SCN1A variants are associated with epilepsy syndromes ranging from mild genetic epilepsy with febrile seizures plus (GEFS+) to severe Dravet syndrome (DS). Many variants are de novo, making early phenotype prediction difficult, and genotype-phenotype associations remain poorly understood.

METHODS

We assessed data from a retrospective cohort of 1018 individuals with SCN1A-related epilepsies. We explored relationships between variant characteristics (position, in silico prediction scores: Combined Annotation Dependent Depletion (CADD), Rare Exome Variant Ensemble Learner (REVEL), SCN1A genetic score), seizure characteristics, and epilepsy phenotype.

RESULTS

DS had earlier seizure onset than other GEFS+ phenotypes (5.3 vs. 12.0 months, p < .001). In silico variant scores were higher in DS versus GEFS+ (p < .001). Patients with missense variants in functionally important regions (conserved N-terminus, S4-S6) exhibited earlier seizure onset (6.0 vs. 7.0 months, p = .003) and were more likely to have DS (280/340); those with missense variants in nonconserved regions had later onset (10.0 vs. 7.0 months, p = .036) and were more likely to have GEFS+ (15/29, χ2 = 19.16, p < .001). A minority of protein-truncating variants were associated with GEFS+ (10/393) and more likely to be located in the proximal first and last exon coding regions than elsewhere in the gene (9.7% vs. 1.0%, p < .001). Carriers of the same missense variant exhibited less variability in age at seizure onset compared with carriers of different missense variants for both DS (1.9 vs. 2.9 months, p = .001) and GEFS+ (8.0 vs. 11.0 months, p = .043). Status epilepticus as presenting seizure type is a highly specific (95.2%) but nonsensitive (32.7%) feature of DS.

SIGNIFICANCE

Understanding genotype-phenotype associations in SCN1A-related epilepsies is critical for early diagnosis and management. We demonstrate an earlier disease onset in patients with missense variants in important functional regions, the occurrence of GEFS+ truncating variants, and the value of in silico prediction scores. Status epilepticus as initial seizure type is a highly specific, but not sensitive, early feature of DS.

Functional requirement of alternative splicing in epithelial-mesenchymal transition of pancreatic circulating tumor

Circulating tumor cells (CTCs) that undergo epithelial-to-mesenchymal transition (EMT) can provide valuable information regarding metastasis and potential therapies. However, current studies on the EMT overlook alternative splicing. Here, we used single-cell full-length transcriptome data and mRNA sequencing of CTCs to identify stage-specific alternative splicing of partial EMT and mesenchymal states during pancreatic cancer metastasis. We classified definitive tumor and normal epithelial cells via genetic aberrations and demonstrated dynamic changes in the epithelial-mesenchymal continuum in both epithelial cancer cells and CTCs. We provide the landscape of alternative splicing in CTCs at different stages of EMT, uncovering cell-type-specific splicing patterns and splicing events in cell surface proteins suitable for therapies. We show that MBNL1 governs cell fate through alternative splicing independently of changes in gene expression and affects the splicing pattern during EMT. We found a high frequency of events that contained multiple premature termination codons and were enriched with C and G nucleotides in close proximity, which influence the likelihood of stop codon readthrough and expand the range of potential therapeutic targets. Our study provides insights into the EMT transcriptome's dynamic changes and identifies potential diagnostic and therapeutic targets in pancreatic cancer.

Heterozygous loss-of-function variants in DOCK4 cause neurodevelopmental delay and microcephaly

Neurons form the basic anatomical and functional structure of the nervous system, and defects in neuronal differentiation or formation of neurites are associated with various psychiatric and neurodevelopmental disorders. Dynamic changes in the cytoskeleton are essential for this process, which is, inter alia, controlled by the dedicator of cytokinesis 4 (DOCK4) through the activation of RAC1. Here, we clinically describe 7 individuals (6 males and one female) with variants in DOCK4 and overlapping phenotype of mild to severe global developmental delay. Additional symptoms include coordination or gait abnormalities, microcephaly, nonspecific brain malformations, hypotonia and seizures. Four individuals carry missense variants (three of them detected de novo) and three individuals carry null variants (two of them maternally inherited). Molecular modeling of the heterozygous missense variants suggests that the majority of them affect the globular structure of DOCK4. In vitro functional expression studies in transfected Neuro-2A cells showed that all missense variants impaired neurite outgrowth. Furthermore, Dock4 knockout Neuro-2A cells also exhibited defects in promoting neurite outgrowth. Our results, including clinical, molecular and functional data, suggest that loss-of-function variants in DOCK4 probable cause a variable spectrum of a novel neurodevelopmental disorder with microcephaly.

Messenger RNA Surveillance: Current Understanding, Regulatory Mechanisms, and Future Implications

Nonsense-mediated mRNA decay (NMD) is an evolutionarily conserved surveillance mechanism in eukaryotes primarily deployed to ensure RNA quality control by eliminating aberrant transcripts and also involved in modulating the expression of several physiological transcripts. NMD, the mRNA surveillance pathway, is a major form of gene regulation in eukaryotes. NMD serves as one of the most significant quality control mechanisms as it primarily scans the newly synthesized transcripts and differentiates the aberrant and non-aberrant transcripts. The synthesis of truncated proteins is restricted, which would otherwise lead to cellular dysfunctions. The up-frameshift factors (UPFs) play a central role in executing the NMD event, largely by recognizing and recruiting multiple protein factors that result in the decay of non-physiological mRNAs. NMD exhibits astounding variability in its ability across eukaryotes in an array of pathological and physiological contexts. The detailed understanding of NMD and the underlying molecular mechanisms remains blurred. This review outlines our current understanding of NMD, in regulating multifaceted cellular events during development and disease. It also attempts to identify unanswered questions that deserve further investigation.

Exome sequencing improves the molecular diagnostics of paediatric unexplained neurodevelopmental disorders

BACKGROUND

Neurodevelopmental disorders (NDDs) and/or associated multiple congenital abnormalities (MCAs) represent a genetically heterogeneous group of conditions with an adverse prognosis for the quality of intellectual and social abilities and common daily functioning. The rapid development of exome sequencing (ES) techniques, together with trio-based analysis, nowadays leads to up to 50% diagnostic yield. Therefore, it is considered as the state-of-the-art approach in these diagnoses.

RESULTS

In our study, we present the results of ES in a cohort of 85 families with 90 children with severe NDDs and MCAs. The interconnection of the in-house bioinformatic pipeline and a unique algorithm for variant prioritization resulted in a diagnostic yield of up to 48.9% (44/90), including rare and novel causative variants (41/90) and intragenic copy-number variations (CNVs) (3/90). Of the total number of 47 causative variants, 53.2% (25/47) were novel, highlighting the clinical benefit of ES for unexplained NDDs. Moreover, trio-based ES was verified as a reliable tool for the detection of rare CNVs, ranging from intragenic exon deletions (GRIN2A, ZC4H2 genes) to a 6-Mb duplication. The functional analysis using PANTHER Gene Ontology confirmed the involvement of genes with causative variants in a wide spectrum of developmental processes and molecular pathways, which form essential structural and functional components of the central nervous system.

CONCLUSION

Taken together, we present one of the first ES studies of this scale from the central European region. Based on the high diagnostic yield for paediatric NDDs in this study, 48.9%, we confirm trio-based ES as an effective and reliable first-tier diagnostic test in the genetic evaluation of children with NDDs.

Phase separation as a possible mechanism for dosage sensitivity

BACKGROUND

Deletion of haploinsufficient genes or duplication of triplosensitive ones results in phenotypic effects in a concentration-dependent manner, and the mechanisms underlying these dosage-sensitive effects remain elusive. Phase separation drives functional compartmentalization of biomolecules in a concentration-dependent manner as well, which suggests a potential link between these two processes, and warrants further systematic investigation.

RESULTS

Here we provide bioinformatic and experimental evidence to show a close link between phase separation and dosage sensitivity. We first demonstrate that haploinsufficient or triplosensitive gene products exhibit a higher tendency to undergo phase separation. Assessing the well-established dosage-sensitive genes HNRNPK, PAX6, and PQBP1 with experiments, we show that these proteins undergo phase separation. Critically, pathogenic variations in dosage-sensitive genes disturb the phase separation process either through reduced protein levels, or loss of phase-separation-prone regions. Analysis of multi-omics data further demonstrates that loss-of-function genetic perturbations on phase-separating genes cause similar dysfunction phenotypes as dosage-sensitive gene perturbations. In addition, dosage-sensitive scores derived from population genetics data predict phase-separating proteins with much better performance than available sequence-based predictors, further illustrating close ties between these two parameters.

CONCLUSIONS

Together, our study shows that phase separation is functionally linked to dosage sensitivity and provides novel insights for phase-separating protein prediction from the perspective of population genetics data.

Systematic analysis of variants escaping nonsense-mediated decay uncovers candidate Mendelian diseases

Protein-truncating variants (PTVs) near the 3' end of genes may escape nonsense-mediated decay (NMD). PTVs in the NMD-escape region (PTVescs) can cause Mendelian disease but are difficult to interpret given their varying impact on protein function. Previously, PTVesc burden was assessed in an epilepsy cohort, but no large-scale analysis has systematically evaluated these variants in rare disease. We performed a retrospective analysis of 29,031 neurodevelopmental disorder (NDD) parent-offspring trios referred for clinical exome sequencing to identify PTVesc de novo mutations (DNMs). We identified 1,376 PTVesc DNMs and 133 genes that were significantly enriched (binomial p < 0.001). The PTVesc-enriched genes included those with PTVescs previously described to cause dominant Mendelian disease (e.g., SEMA6B, PPM1D, and DAGLA). We annotated ClinVar variants for PTVescs and identified 948 genes with at least one high-confidence pathogenic variant. Twenty-two known Mendelian PTVesc-enriched genes had no prior evidence of PTVesc-associated disease. We found 22 additional PTVesc-enriched genes that are not well established to be associated with Mendelian disease, several of which showed phenotypic similarity between individuals harboring PTVesc variants in the same gene. Four individuals with PTVesc mutations in RAB1A had similar phenotypes including NDD and spasticity. PTVesc mutations in IRF2BP1 were found in two individuals who each had severe immunodeficiency manifesting in NDD. Three individuals with PTVesc mutations in LDB1 all had NDD and multiple congenital anomalies. Using a large-scale, systematic analysis of DNMs, we extend the mutation spectrum for known Mendelian disease-associated genes and identify potentially novel disease-associated genes.

Post-transcriptional Regulation of Gene Expression via Unproductive Splicing

Unproductive splicing is a mechanism of post-transcriptional gene expression control in which premature stop codons are inserted into protein-coding transcripts as a result of regulated alternative splicing, leading to their degradation via the nonsense-mediated decay pathway. This mechanism is especially characteristic of RNA-binding proteins, which regulate each other's expression levels and those of other genes in multiple auto- and cross-regulatory loops. Deregulation of unproductive splicing is a cause of serious human diseases, including cancers, and is increasingly being considered as a prominent therapeutic target. This review discusses the types of unproductive splicing events, the mechanisms of auto- and cross-regulation, nonsense-mediated decay escape, and problems in identifying unproductive splice isoforms. It also provides examples of deregulation of unproductive splicing in human diseases and discusses therapeutic strategies for its correction using antisense oligonucleotides and small molecules.

En Route to Targeted Ribosome Editing to Replenish Skin Anchor Protein LAMB3 in Junctional Epidermolysis Bullosa

Severe junctional epidermolysis bullosa is a rare genetic, postpartum lethal skin disease, predominantly caused by nonsense/premature termination codon (PTC) sequence variants in LAMB3 gene. LAMB3 encodes LAMB3, the β subunit of epidermal-dermal skin anchor laminin 332. Most translational reads of a PTC mRNA deliver truncated, nonfunctional proteins, whereas an endogenous PTC readthrough mechanism produces full-length protein at minimal and insufficient levels. Conventional translational readthrough-inducing drugs amplify endogenous PTC readthrough; however, translational readthrough-inducing drugs are either proteotoxic or nonselective. Ribosome editing is a more selective and less toxic strategy. This technique identified ribosomal protein L35/uL29 (ie, RpL35) and RpL35-ligands repurposable drugs artesunate and atazanavir as molecular tools to increase production levels of full-length LAMB3. To evaluate ligand activity in living cells, we monitored artesunate and atazanavir treatment by dual luciferase reporter assays. Production levels of full-length LAMB3 increased up to 200% upon artesunate treatment, up to 150% upon atazanavir treatment, and up to 170% upon combinatorial treatment of RpL35 ligands at reduced drug dosage, with an unrelated PTC reporter being nonresponsive. Proof of bioactivity of RpL35 ligands in selective increase of full-length LAMB3 provides the basis for an alternative, targeted therapeutic route to replenish LAMB3 in severe junctional epidermolysis bullosa.

Characterization of NEB mutations in patients reveals novel nemaline myopathy disease mechanisms and omecamtiv mecarbil force effects

Nebulin, a critical protein of the skeletal muscle thin filament, plays important roles in physiological processes such as regulating thin filament length (TFL), cross-bridge cycling, and myofibril alignment. Mutations in the nebulin gene ( NEB ) cause NEB-based nemaline myopathy (NEM2), a genetically heterogeneous disorder characterized by hypotonia and muscle weakness, currently lacking therapies targeting the underlying pathological mechanisms. In this study, we examined a cohort of ten NEM2 patients, each with unique mutations, aiming to understand their impact on mRNA, protein, and functional levels. Results show that truncation mutations affect NEB mRNA stability and lead to nonsense-mediated decay of the mutated transcript. Moreover, a high incidence of cryptic splice site activation was found in patients with splicing mutations which is expected to disrupt the actin-binding sites of nebulin. Determination of protein levels revealed patients with relatively normal nebulin levels and others with markedly reduced nebulin. We observed a positive relation between the reduction in nebulin and a reduction in TFL, and a positive relation between the reduction in nebulin level and the reduction in tension (both maximal and submaximal tension). Interestingly, our study revealed a duplication mutation in nebulin that resulted in a larger nebulin protein and longer TFL. Additionally, we investigated the effect of Omecamtiv mecarbil (OM), a small-molecule activator of cardiac myosin, on force production of type I muscle fibers of NEM2 patients. OM treatment substantially increased submaximal tension across all NEM2 patients ranging from 87-318%, with the largest effects in patients with the lowest level of nebulin. In summary, this study indicates that post-transcriptional or post-translational mechanisms regulate nebulin expression. Moreover, we propose that the pathomechanism of NEM2 involves not only shortened but also elongated thin filaments, along with the disruption of actin-binding sites resulting from splicing mutations. Significantly, our findings highlight the potential of OM treatment to improve skeletal muscle function in NEM2 patients, especially those with large reductions in nebulin levels.

Genotype-phenotype associations in Alström syndrome: a systematic review and meta-analysis

BACKGROUND

Alström syndrome (ALMS; #203800) is an ultrarare monogenic recessive disease. This syndrome is associated with variants in the ALMS1 gene, which encodes a centrosome-associated protein involved in the regulation of several ciliary and extraciliary processes, such as centrosome cohesion, apoptosis, cell cycle control and receptor trafficking. The type of variant associated with ALMS is mostly complete loss-of-function variants (97%) and they are mainly located in exons 8, 10 and 16 of the gene. Other studies in the literature have tried to establish a genotype-phenotype correlation in this syndrome with limited success. The difficulty in recruiting a large cohort in rare diseases is the main barrier to conducting this type of study.

METHODS

In this study we collected all cases of ALMS published to date. We created a database of patients who had a genetic diagnosis and an individualised clinical history. Lastly, we attempted to establish a genotype-phenotype correlation using the truncation site of the patient's longest allele as a grouping criteria.

RESULTS

We collected a total of 357 patients, of whom 227 had complete clinical information, complete genetic diagnosis and meta-information on sex and age. We have seen that there are five variants with high frequency, with p.(Arg2722Ter) being the most common variant, with 28 alleles. No gender differences in disease progression were detected. Finally, truncating variants in exon 10 seem to be correlated with a higher prevalence of liver disorders in patients with ALMS.

CONCLUSION

Pathogenic variants in exon 10 of the ALMS1 gene were associated with a higher prevalence of liver disease. However, the location of the variant in the ALMS1 gene does not have a major impact on the phenotype developed by the patient.

Exploring the relationship between IGHMBP2 gene mutations and spinal muscular atrophy with respiratory distress type 1 and Charcot-Marie-Tooth disease type 2S: a systematic review

Background

IGHMBP2 is a crucial gene for the development and maintenance of the nervous system, especially in the survival of motor neurons. Mutations in this gene have been associated with spinal muscular atrophy with respiratory distress type 1 (SMARD1) and Charcot-Marie-Tooth disease type 2S (CMT2S).

Methods

We conducted a systematic literature search using the PubMed database to identify studies published up to April 1st, 2023, that investigated the association between IGHMBP2 mutations and SMARD1 or CMT2S. We compared the non-truncating mutations and truncating mutations of the IGHMBP2 gene and selected high-frequency mutations of the IGHMBP2 gene.

Results

We identified 52 articles that investigated the association between IGHMBP2 mutations and SMARD1/CMT2S. We found 6 hotspot mutations of the IGHMBP2 gene. The truncating mutations in trans were all associated with SMARD1.

Conclusion

This study provides evidence that the complete LOF mechanism of the IGHMBP2 gene defect may be an important cause of SMARD1.

Mutation Impact on mRNA Versus Protein Expression across Human Cancers

Cancer mutations are often assumed to alter proteins, thus promoting tumorigenesis. However, how mutations affect protein expression has rarely been systematically investigated. We conduct a comprehensive analysis of mutation impacts on mRNA- and protein-level expressions of 953 cancer cases with paired genomics and global proteomic profiling across six cancer types. Protein-level impacts are validated for 47.2% of the somatic expression quantitative trait loci (seQTLs), including mutations from likely "long-tail" driver genes. Devising a statistical pipeline for identifying somatic protein-specific QTLs (spsQTLs), we reveal several gene mutations, including NF1 and MAP2K4 truncations and TP53 missenses showing disproportional influence on protein abundance not readily explained by transcriptomics. Cross-validating with data from massively parallel assays of variant effects (MAVE), TP53 missenses associated with high tumor TP53 proteins were experimentally confirmed as functional. Our study demonstrates the importance of considering protein-level expression to validate mutation impacts and identify functional genes and mutations.

Nonsense-Mediated mRNA Decay: Mechanistic Insights and Physiological Significance

Nonsense-mediated mRNA decay (NMD) is an evolutionarily conserved surveillance mechanism across eukaryotes and also regulates the expression of physiological transcripts, thus involved in gene regulation. It essentially ensures recognition and removal of aberrant transcripts. Therefore, the NMD protects the cellular system by restricting the synthesis of truncated proteins, potentially by eliminating the faulty mRNAs. NMD is an evolutionarily conserved surveillance mechanism across eukaryotes and also regulates the expression of physiological transcripts, thus involved in gene regulation as well. Primarily, the NMD machinery scans and differentiates the aberrant and non-aberrant transcripts. A myriad of cellular dysfunctions arise due to production of truncated proteins, so the NMD core proteins, the up-frameshift factors (UPFs) recognizes the faulty mRNAs and further recruits factors resulting in the mRNA degradation. NMD exhibits astounding variability in its ability in regulating cellular mechanisms including both pathological and physiological events. But, the detailed underlying molecular mechanisms in NMD remains blurred and require extensive investigation to gain insights on cellular homeostasis. The complexity in understanding of NMD pathway arises due to the involvement of numerous proteins, molecular interactions and their functioning in different steps of this process. Moreover methods such as alternative splicing generates numerous isoforms of mRNA, so it makes difficulties in understanding the impact of alternative splicing on the efficiency of NMD functioning. Role of NMD in cancer development is very complex. Studies have shown that in some cases cancer cells use NMD pathway as a tool to exploit the NMD mechanism to maintain tumor microenvironment. A greater level of understanding about the intricate mechanism of how tumor used NMD pathway for their benefits, a strategy can be developed for targeting and inhibiting NMD factors involved in pro-tumor activity. There are very little amount of information available about the NMD pathway, how it discriminate mRNAs that are targeted by NMD from those that are not. This review highlights our current understanding of NMD, specifically the regulatory mechanisms and attempts to outline less explored questions that warrant further investigations. Taken as a whole, a detailed molecular understanding of the NMD mechanism could lead to wide-ranging applications for improving cellular homeostasis and paving out strategies in combating pathological disorders leaping forward toward achieving United Nations sustainable development goals (SDG 3: Good health and well-being).

CUX1-related neurodevelopmental disorder: deep insights into phenotype-genotype spectrum and underlying pathology

Heterozygous, pathogenic CUX1 variants are associated with global developmental delay or intellectual disability. This study delineates the clinical presentation in an extended cohort and investigates the molecular mechanism underlying the disorder in a Cux1+/- mouse model. Through international collaboration, we assembled the phenotypic and molecular information for 34 individuals (23 unpublished individuals). We analyze brain CUX1 expression and susceptibility to epilepsy in Cux1+/- mice. We describe 34 individuals, from which 30 were unrelated, with 26 different null and four missense variants. The leading symptoms were mild to moderate delayed speech and motor development and borderline to moderate intellectual disability. Additional symptoms were muscular hypotonia, seizures, joint laxity, and abnormalities of the forehead. In Cux1+/- mice, we found delayed growth, histologically normal brains, and increased susceptibility to seizures. In Cux1+/- brains, the expression of Cux1 transcripts was half of WT animals. Expression of CUX1 proteins was reduced, although in early postnatal animals significantly more than in adults. In summary, disease-causing CUX1 variants result in a non-syndromic phenotype of developmental delay and intellectual disability. In some individuals, this phenotype ameliorates with age, resulting in a clinical catch-up and normal IQ in adulthood. The post-transcriptional balance of CUX1 expression in the heterozygous brain at late developmental stages appears important for this favorable clinical course.

Reclassification of variants of tumor suppressor genes based on Sanger RNA sequencing without NMD inhibition

Introduction: RNA sequence analysis can be effectively used to identify aberrant splicing, and tumor suppressor genes are adequate targets considering their loss-of-function mechanisms. Sanger sequencing is the simplest method for RNA sequence analysis; however, because of its insufficient sensitivity in cases with nonsense-mediated mRNA decay (NMD), the use of cultured specimens with NMD inhibition has been recommended, hindering its wide adoption. Method: The results of Sanger sequencing of peripheral blood RNA without NMD inhibition performed on potential splicing variants of tumor suppressor genes were retrospectively reviewed. For negative cases, in which no change was identified in the transcript, the possibility of false negativity caused by NMD was assessed through a review of the up-to-date literature. Results: Eleven potential splice variants of various tumor suppressor genes were reviewed. Six variants were classified as pathogenic or likely pathogenic based on the nullifying effect identified by Sanger RNA sequencing. Four variants remained as variants of uncertain significance because of identified in-frame changes or normal expression of both alleles. The result of one variant was suspected to be a false negative caused by NMD after reviewing a recent study that reported the same variant as causing a nullifying effect on the affected transcript. Conclusion: Although RNA changes found in the majority of cases were expected to undergo NMD by canonical rules, most cases (10/11) were interpretable by Sanger RNA sequencing without NMD inhibition due to incomplete NMD efficiency or allele-specific expression despite highly efficient NMD.

Translation-coupled mRNA quality control mechanisms

mRNA surveillance pathways are essential for accurate gene expression and to maintain translation homeostasis, ensuring the production of fully functional proteins. Future insights into mRNA quality control pathways will enable us to understand how cellular mRNA levels are controlled, how defective or unwanted mRNAs can be eliminated, and how dysregulation of these can contribute to human disease. Here we review translation-coupled mRNA quality control mechanisms, including the non-stop and no-go mRNA decay pathways, describing their mechanisms, shared trans-acting factors, and differences. We also describe advances in our understanding of the nonsense-mediated mRNA decay (NMD) pathway, highlighting recent mechanistic findings, the discovery of novel factors, as well as the role of NMD in cellular physiology and its impact on human disease.

The penetrance of rare variants in cardiomyopathy-associated genes: A cross-sectional approach to estimating penetrance for secondary findings

Understanding the penetrance of pathogenic variants identified as secondary findings (SFs) is of paramount importance with the growing availability of genetic testing. We estimated penetrance through large-scale analyses of individuals referred for diagnostic sequencing for hypertrophic cardiomyopathy (HCM; 10,400 affected individuals, 1,332 variants) and dilated cardiomyopathy (DCM; 2,564 affected individuals, 663 variants), using a cross-sectional approach comparing allele frequencies against reference populations (293,226 participants from UK Biobank and gnomAD). We generated updated prevalence estimates for HCM (1:543) and DCM (1:220). In aggregate, the penetrance by late adulthood of rare, pathogenic variants (23% for HCM, 35% for DCM) and likely pathogenic variants (7% for HCM, 10% for DCM) was substantial for dominant cardiomyopathy (CM). Penetrance was significantly higher for variant subgroups annotated as loss of function or ultra-rare and for males compared to females for variants in HCM-associated genes. We estimated variant-specific penetrance for 316 recurrent variants most likely to be identified as SFs (found in 51% of HCM- and 17% of DCM-affected individuals). 49 variants were observed at least ten times (14% of affected individuals) in HCM-associated genes. Median penetrance was 14.6% (±14.4% SD). We explore estimates of penetrance by age, sex, and ancestry and simulate the impact of including future cohorts. This dataset reports penetrance of individual variants at scale and will inform the management of individuals undergoing genetic screening for SFs. While most variants had low penetrance and the costs and harms of screening are unclear, some individuals with highly penetrant variants may benefit from SFs.

aenmd: annotating escape from nonsense-mediated decay for transcripts with protein-truncating variants

SUMMARY

DNA changes that cause premature termination codons (PTCs) represent a large fraction of clinically relevant pathogenic genomic variation. Typically, PTCs induce transcript degradation by nonsense-mediated mRNA decay (NMD) and render such changes loss-of-function alleles. However, certain PTC-containing transcripts escape NMD and can exert dominant-negative or gain-of-function (DN/GOF) effects. Therefore, systematic identification of human PTC-causing variants and their susceptibility to NMD contributes to the investigation of the role of DN/GOF alleles in human disease. Here we present aenmd, a software for annotating PTC-containing transcript-variant pairs for predicted escape from NMD. aenmd is user-friendly and self-contained. It offers functionality not currently available in other methods and is based on established and experimentally validated rules for NMD escape; the software is designed to work at scale, and to integrate seamlessly with existing analysis workflows. We applied aenmd to variants in the gnomAD, Clinvar, and GWAS catalog databases and report the prevalence of human PTC-causing variants in these databases, and the subset of these variants that could exert DN/GOF effects via NMD escape.

AVAILABILITY AND IMPLEMENTATION

aenmd is implemented in the R programming language. Code is available on GitHub as an R-package (github.com/kostkalab/aenmd.git), and as a containerized command-line interface (github.com/kostkalab/aenmd_cli.git).