A missense mutation in SNRPE linked to non-syndromal microcephaly interferes with U snRNP assembly and pre-mRNA splicing

Machine learning-based gene expression biomarkers to distinguish Zika and Dengue virus infections: implications for diagnosis

Zika virus (ZIKV) and Dengue virus (DENV) infections cause severe disease in humans and are significant socio-economic burden worldwide. These flavivirus infections are difficult to diagnose serologically due to antigenic overlap. The phylogenetic analysis shows that ZIKV clusters with DENVs at a higher node of the phylogenetic tree with significant genomic and structural similarity. Our study aims to identify gene biomarkers for the classification of Dengue and Zika viral infections using machine learning algorithms and bioinformatics analysis. The gene expression count matrix for single-cell RNA sequencing dataset GSE110496 was analyzed using binary classifiers, namely Logistic regression, Support Vector Machines, Random Forest, and Decision trees. The GSE110496 dataset represents a unique study of the transcriptional and translational dynamics of DENV and ZIKV infections at 4-, 12-, 24-, and 48-h time points for human hepatoma (Huh7) cells. Out of which 24-h time point has been analyzed in this study, at the optimal threshold of viral molecules. Feature selection was performed using two different approaches Random Forest Classifier (RFC) for gene ranking and Recursive Feature Elimination (RFE). Out of which RFE, showed more accuracy and precision. The classification accuracy of 89.4% and the precision of 90% were obtained using selected 10 gene features. SCY1 Like Pseudokinase 3 (SCYL3), Chromosome 1 Open Reading Frame 112 (C1orf112), Complement factor H (CFH), Heme-binding protein 1 (HEBP1), Cadherin 1 (CDH1), Nibrin (NBN), Histone deacetylase 5 (HDAC5), nuclear receptor subfamily 0, group B, member 2 (NR0B2), Annexin A9 (ANXA9) and Alcohol dehydrogenase 6 (ADH6) are the proposed gene biomarkers in this study. The functional analysis of the reported biomarkers was performed using KEGG and GO with the WEB-based Gene SeT AnaLysis Toolkit (WebGestalt). The relationship of the selected biomarkers with DENV and ZIKV infections analyzed using a gene-gene interaction network showed important interactions for viral entry, replication, translation, and metabolic pathways. These biomarkers are potential diagnostic markers for DENV and ZIKV infections based on machine learning analysis and need further experimental validation.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s13337-024-00885-8.

Two High-Quality Cygnus Genome Assemblies Reveal Genomic Variations Associated with Plumage Color

As an exemplary model for examining molecular mechanisms responsible for extreme phenotypic variations, plumage color has garnered significant interest. The Cygnus genus features two species, Cygnus olor and Cygnus atratus, that exhibit striking disparities in plumage color. However, the molecular foundation for this differentiation has remained elusive. Herein, we present two high-quality genomes for C. olor and C. atratus, procured using the Illumina and Nanopore technologies. The assembled genome of C. olor was 1.12 Gb in size with a contig N50 of 26.82 Mb, while its counterpart was 1.13 Gb in size with a contig N50 of 21.91 Mb. A comparative analysis unveiled three genes (TYR, SLC45A2, and SLC7A11) with structural variants in the melanogenic pathway. Notably, we also identified a novel gene, PWWP domain containing 2A (PWWP2A), that is related to plumage color, for the first time. Using targeted gene modification analysis, we demonstrated the potential genetic effect of the PWWP2A variant on pigment gene expression and melanin production. Finally, our findings offer insight into the intricate pattern of pigmentation and the role of polygenes in birds. Furthermore, these two high-quality genome references provide a comprehensive resource and perspective for comparative functional and genetic studies of evolution within the Cygnus genus.

mRNA isoform balance in neuronal development and disease

Balanced mRNA isoform diversity and abundance are spatially and temporally regulated throughout cellular differentiation. The proportion of expressed isoforms contributes to cell type specification and determines key properties of the differentiated cells. Neurons are unique cell types with intricate developmental programs, characteristic cellular morphologies, and electrophysiological potential. Neuron-specific gene expression programs establish these distinctive cellular characteristics and drive diversity among neuronal subtypes. Genes with neuron-specific alternative processing are enriched in key neuronal functions, including synaptic proteins, adhesion molecules, and scaffold proteins. Despite the similarity of neuronal gene expression programs, each neuronal subclass can be distinguished by unique alternative mRNA processing events. Alternative processing of developmentally important transcripts alters coding and regulatory information, including interaction domains, transcript stability, subcellular localization, and targeting by RNA binding proteins. Fine-tuning of mRNA processing is essential for neuronal activity and maintenance. Thus, the focus of neuronal RNA biology research is to dissect the transcriptomic mechanisms that underlie neuronal homeostasis, and consequently, predispose neuronal subtypes to disease. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA in Disease and Development > RNA in Development.

Mutation of Arabidopsis SME1 and Sm core assembly improves oxidative stress resilience

Alternative splicing is a key posttranscriptional gene regulatory process, acting in diverse adaptive and basal plant processes. Splicing of precursor-messenger RNA (pre-mRNA) is catalyzed by a dynamic ribonucleoprotein complex, designated the spliceosome. In a suppressor screen, we identified a nonsense mutation in the Smith (Sm) antigen protein SME1 to alleviate photorespiratory H₂O₂-dependent cell death in catalase deficient plants. Similar attenuation of cell death was observed upon chemical inhibition of the spliceosome, suggesting pre-mRNA splicing inhibition to be responsible for the observed cell death alleviation. Furthermore, the sme1-2 mutants showed increased tolerance to the reactive oxygen species inducing herbicide methyl viologen. Both an mRNA-seq and shotgun proteomic analysis in sme1-2 mutants displayed a constitutive molecular stress response, together with extensive alterations in pre-mRNA splicing of transcripts encoding metabolic enzymes and RNA binding proteins, even under unstressed conditions. Using SME1 as a bait to identify protein interactors, we provide experimental evidence for almost 50 homologs of the mammalian spliceosome-associated protein to reside in the Arabidopsis thaliana spliceosome complexes and propose roles in pre-mRNA splicing for four uncharacterized plant proteins. Furthermore, as for sme1-2, a mutant in the Sm core assembly protein ICLN resulted in a decreased sensitivity to methyl viologen. Taken together, these data show that both a perturbed Sm core composition and assembly results in the activation of a defense response and in enhanced resilience to oxidative stress.

Novel blended SNRPE-related spliceosomopathy phenotype characterized by microcephaly and congenital atrichia

Variants in genes encoding core components of the spliceosomes are associated with craniofacial syndromes, collectively called craniofacial spliceosomopathies. SNRPE encodes a core component of pre-mRNA processing U-rich small nuclear ribonuclear proteins (UsnRNPs). Heterozygous variants in SNRPE have been reported in six families with isolated hypotrichosis simplex in addition to one case of isolated non syndromic congenital microcephaly. Here, we report a patient with a novel blended phenotype of microcephaly and congenital atrichia with multiple congenital anomalies due to a de novo intronic SNRPE deletion, c.82-28_82-16del, which results in exon skipping. As discussed within, this phenotype, which we propose be named SNRPE-related syndromic microcephaly and hypotrichosis, overlaps other craniofacial splicesosomopathies.

Congenital Microcephaly: A Debate on Diagnostic Challenges and Etiological Paradigm of the Shift from Isolated/Non-Syndromic to Syndromic Microcephaly

Congenital microcephaly (CM) exhibits broad clinical and genetic heterogeneity and is thus categorized into several subtypes. However, the recent bloom of disease-gene discoveries has revealed more overlaps than differences in the underlying genetic architecture for these clinical sub-categories, complicating the differential diagnosis. Moreover, the mechanism of the paradigm shift from a brain-restricted to a multi-organ phenotype is only vaguely understood. This review article highlights the critical factors considered while defining CM subtypes. It also presents possible arguments on long-standing questions of the brain-specific nature of CM caused by a dysfunction of the ubiquitously expressed proteins. We argue that brain-specific splicing events and organ-restricted protein expression may contribute in part to disparate clinical manifestations. We also highlight the role of genetic modifiers and de novo variants in the multi-organ phenotype of CM and emphasize their consideration in molecular characterization. This review thus attempts to expand our understanding of the phenotypic and etiological variability in CM and invites the development of more comprehensive guidelines.

The emerging significance of splicing in vertebrate development

Splicing is a crucial regulatory node of gene expression that has been leveraged to expand the proteome from a limited number of genes. Indeed, the vast increase in intron number that accompanied vertebrate emergence might have aided the evolution of developmental and organismal complexity. Here, we review how animal models for core spliceosome components have provided insights into the role of splicing in vertebrate development, with a specific focus on neuronal, neural crest and skeletal development. To this end, we also discuss relevant spliceosomopathies, which are developmental disorders linked to mutations in spliceosome subunits. Finally, we discuss potential mechanisms that could underlie the tissue-specific phenotypes often observed upon spliceosome inhibition and identify gaps in our knowledge that, we hope, will inspire further research.

The comprehensive expression and functional analysis of m6A modification "readers" in hepatocellular carcinoma

N6-methyladenosine (m6A) modification regulators are essential for the diagnosis and treatment of various cancers. However, the comprehensive analysis about roles of m6A "readers" in hepatocellular carcinoma (HCC) remains unclear. UALCAN, GEPIA2, HPA, Kaplan Meier plotter, cBioPortal, STRING WebGestalt, Metascape and TIMER 2.0 database and Cytoscape software were used to comprehensively analyze the bioinformatic data. We found that m6A "readers" were upregulated at the mRNA level and protein level in HCC patients. Highly expressed YTHDF1, IGF2BP3 and NKAP were positively correlated with advanced HCC stage and had a poor prognosis in OS and PFS. The gene alterations of m6A "readers" happened frequently, and YTHDF3 had the highest mutation rate. The function of m6A "readers" on HCC may be closely correlated with splicing related proteins (including HNRNP family, SNRP family, and SR family), metabolic process, protein binding and RNA splicing related signaling pathways. Moreover, although the correlation of YTHDF3 and CD8+ T cell infiltration, and the correlation of IGF2BP3 and infiltration of mast cells and CAF are negative, most m6A "readers" had a positive correlation with immune cells (including CD8+ T cell, CD4+ T cell, Tregs, B cell, neutrophil, monocyte, macrophage, myeloid dendritic cell, nature killer cell, mast cell, and CAF). Macrophages, CD4+ T cell, Treg, B cell, monocyte, and myeloid dendritic cell had a positively strong correlation (Rho>0.4) with most m6A "readers" (such as YTHDC1, YTHDC2, YTHDF1, IGF2BP3, HNRNPA2B1 and HNRNPC). In conclusion, by comprehensive analysis of m6A "readers", we found that they were involved in the prognosis of HCC, and m6A "readers" might regulate the development and progression of HCC by participating in some metabolism-related and RNA splicing-related signaling pathways as well as immune cell infiltration.

Pan-cancer pervasive upregulation of 3' UTR splicing drives tumourigenesis

Most mammalian genes generate messenger RNAs with variable untranslated regions (UTRs) that are important post-transcriptional regulators. In cancer, shortening at 3' UTR ends via alternative polyadenylation can activate oncogenes. However, internal 3' UTR splicing remains poorly understood as splicing studies have traditionally focused on protein-coding alterations. Here we systematically map the pan-cancer landscape of 3' UTR splicing and present this in SpUR ( http://www.cbrc.kaust.edu.sa/spur/home/ ). 3' UTR splicing is widespread, upregulated in cancers, correlated with poor prognosis and more prevalent in oncogenes. We show that antisense oligonucleotide-mediated inhibition of 3' UTR splicing efficiently reduces oncogene expression and impedes tumour progression. Notably, CTNNB1 3' UTR splicing is the most consistently dysregulated event across cancers. We validate its upregulation in hepatocellular carcinoma and colon adenocarcinoma, and show that the spliced 3' UTR variant is the predominant contributor to its oncogenic functions. Overall, our study highlights the importance of 3' UTR splicing in cancer and may launch new avenues for RNA-based anti-cancer therapeutics.

Genome-wide analysis of runs of homozygosity in Italian Mediterranean buffalo

Runs of homozygosity (ROH) are a powerful tool to explore patterns of genomic inbreeding in animal populations and detect signatures of selection. The present study used ROH analysis to evaluate the genome-wide patterns of homozygosity, inbreeding levels, and distribution of ROH islands using the SNP data sets from 899 Mediterranean buffaloes. A total of 42,433 ROH segments were identified, with an average of 47.20 segments per individual. The ROH comprising mostly shorter segments (1-4 Mb) accounted for approximately 72.29% of all ROH. In contrast, the larger ROH (>8 Mb) class accounted for only 7.97% of all ROH segments. Estimated inbreeding coefficients from ROH (F_ROH) ranged from 0.0201 to 0.0371. Pearson correlations between F_ROH and genomic relationship matrix increased with the increase of ROH length. We identified ROH hotspots in 12 genomic regions, located on chromosomes 1, 2, 3, 5, 17, and 19, harboring a total of 122 genes. Protein-protein interaction (PPI) analysis revealed the clustering of these genes into 7 PPI networks. Many genes located in these regions were associated with different production traits. In addition, 5 ROH islands overlapped with cattle quantitative trait loci that were mainly associated with milk traits. These findings revealed the genome-wide autozygosity patterns and inbreeding levels in Mediterranean buffalo. Our study identified many candidate genes related to production traits that could be used to assist in selective breeding for genetic improvement of buffalo.

Mammalian splicing divergence is shaped by drift, buffering in trans, and a scaling law

This study globally investigates the allelic splicing pattern in multiple tissues of an F1 hybrid mouse and reveals the underlying driving forces shaping such tissue-dependent splicing divergence.

Alternative splicing is ubiquitous, but the mechanisms underlying its pattern of evolutionary divergence across mammalian tissues are still underexplored. Here, we investigated the cis-regulatory divergences and their relationship with tissue-dependent trans-regulation in multiple tissues of an F1 hybrid between two mouse species. Large splicing changes between tissues are highly conserved and likely reflect functional tissue-dependent regulation. In particular, micro-exons frequently exhibit this pattern with high inclusion levels in the brain. Cis-divergence of splicing appears to be largely non-adaptive. Although divergence is in general associated with higher densities of sequence variants in regulatory regions, events with high usage of the dominant isoform apparently tolerate more mutations, explaining why their exon sequences are highly conserved but their intronic splicing site flanking regions are not. Moreover, we demonstrate that non-adaptive mutations are often masked in tissues where accurate splicing likely is more important, and experimentally attribute such buffering effect to trans-regulatory splicing efficiency.

Update of recent findings in genetic hair disorders

Genetic hair disorders, although unusual, are not very rare, and dermatologists often have opportunities to see patients. Significant advances in molecular genetics have led to identifying many causative genes for genetic hair disorders, including the recently identified causative genes, such as LSS and C3ORF52. Many patients have been detected with autosomal recessive woolly hair/hypotrichosis in the Japanese population caused by founder mutations in the LIPH gene. Additionally, many patients with genetic hair disorders caused by other genes have been reported in East Asia including Japan. Understanding genetic hair disorders is essential for dermatologists, and the findings obtained from analyzing these diseases will contribute to revealing the mechanisms of hair follicle morphogenesis and development in humans.

Do Fragile X Syndrome and Other Intellectual Disorders Converge at Aberrant Pre-mRNA Splicing?

Fragile X Syndrome is a neuro-developmental disorder caused by the silencing of the FMR1 gene, resulting in the loss of its protein product, FMRP. FMRP binds mRNA and represses general translation in the brain. Transcriptome analysis of the Fmr1-deficient mouse hippocampus reveals widespread dysregulation of alternative splicing of pre-mRNAs. Many of these aberrant splicing changes coincide with those found in post-mortem brain tissue from individuals with autism spectrum disorders (ASDs) as well as in mouse models of intellectual disability such as PTEN hamartoma syndrome (PHTS) and Rett Syndrome (RTT). These splicing changes could result from chromatin modifications (e.g., in FXS, RTT) and/or splicing factor alterations (e.g., PTEN, autism). Based on the identities of the RNAs that are mis-spliced in these disorders, it may be that they are at least partly responsible for some shared pathophysiological conditions. The convergence of splicing aberrations among these autism spectrum disorders might be crucial to understanding their underlying cognitive impairments.

The Identification of RNA-Binding Proteins Functionally Associated with Tumor Progression in Gastrointestinal Cancer

Simple Summary

Previous investigations described bioinformatic analyses based on the mRNA expression and somatic mutation as useful strategies for identifying cancer-associated molecules that were potential candidates of therapeutic targets. However, these data included secondary changes and non-functional alterations that do not influence tumor progression. Investigations, including our own studies, have shown that some RBPs shuttle cytoplasm and nuclei, and their affinity to RNAs is regulated by posttranslational modifications, such as phosphorylation. Therefore, the functional assessment of individual molecules is the most suitable strategy for identifying cancer-associated genes with or without expressional changes and mutations. This report showed for the first time that a functional assessment using an siRNA library was useful for identifying therapeutic targets from molecular groups, including RBPs, that had not been identified by expressional and mutational analyses.

Abstract

Previous investigations have indicated that RNA-binding proteins (RBPs) are key molecules for the development of organs, differentiation, cell growth and apoptosis in cancer cells as well as normal cells. A bioinformatics analysis based on the mRNA expression and a somatic mutational database revealed the association between aberrant expression/mutations of RBPs and cancer progression. However, this method failed to detect functional alterations in RBPs without changes in the expression, thus leading to false negatives. To identify major tumor-associated RBPs, we constructed an siRNA library based on the database of RBPs and assessed the influence on the growth of colorectal, pancreatic and esophageal cancer cells. A comprehensive analysis of siRNA functional screening findings using 1198 siRNAs targeting 416 RBPs identified 41 RBPs in which 50% inhibition of cell growth was observed in cancer cells. Among these RBPs, 12 showed no change in the mRNA expression and no growth suppression in non-cancerous cells when downregulated by specific siRNAs. We herein report for the first time cancer-promotive RBPs identified by a novel functional assessment using an siRNA library of RBPs combined with expressional and mutational analyses.

Biology of the mRNA Splicing Machinery and Its Dysregulation in Cancer Providing Therapeutic Opportunities

Dysregulation of messenger RNA (mRNA) processing—in particular mRNA splicing—is a hallmark of cancer. Compared to normal cells, cancer cells frequently present aberrant mRNA splicing, which promotes cancer progression and treatment resistance. This hallmark provides opportunities for developing new targeted cancer treatments. Splicing of precursor mRNA into mature mRNA is executed by a dynamic complex of proteins and small RNAs called the spliceosome. Spliceosomes are part of the supraspliceosome, a macromolecular structure where all co-transcriptional mRNA processing activities in the cell nucleus are coordinated. Here we review the biology of the mRNA splicing machinery in the context of other mRNA processing activities in the supraspliceosome and present current knowledge of its dysregulation in lung cancer. In addition, we review investigations to discover therapeutic targets in the spliceosome and give an overview of inhibitors and modulators of the mRNA splicing process identified so far. Together, this provides insight into the value of targeting the spliceosome as a possible new treatment for lung cancer.

PSMC2, ORC5 and KRTDAP are specific biomarkers for HPV-negative head and neck squamous cell carcinoma

The prognosis of patients with human papillomavirus (HPV)-negative head and neck squamous cell carcinoma (HNSCC) is poorer than those with HPV-positive HNSCC. The present study aimed to identify novel and specific biomarkers of HPV-negative HNSCC using bioinformatics analysis and associated experiments. The gene expression profiles of HPV-negative HNSCC tissues and corresponding clinical data were downloaded from The Cancer Genome Atlas database and used in a weighted gene co-expression network analysis. Genes in clinically significant co-expression modules were used to construct a protein-protein interaction (PPI) network. The genes demonstrating a high degree score in the PPI network and a high correlation with tumor grade were considered hub genes. The diagnostic value of the hub genes associated with HPV-negative and HPV-positive HNSCC was analyzed using differential expression gene (DEG) analysis, immunohistochemical (IHC) staining and a receiver operating characteristic (ROC) curve analysis. Seven genes [Serrate RNA effector molecule (SRRT), checkpoint kinase 2 (CHEK2), small nuclear ribonucleoprotein polypeptide E (SNRPE), proteasome 26S subunit ATPase 2 (PSMC2), origin recognition complex subunit 5 (ORC5), S100 calcium binding protein A7 and keratinocyte differentiation associated protein (KRTDAP)] were demonstrated to be hub genes in clinically significant co-expression modules. DEG, IHC and ROC curve analyses revealed that SRRT, CHEK2 and SNRPE were significantly upregulated in HPV-negative and HPV-positive HNSCC tissues compared with in adjacent tissues, and these genes demonstrated a high diagnostic value for distinguishing HNSCC tissues. However, PSMC2, ORC5 and KRTDAP were the only differentially expressed genes identified in HPV-negative HNSCC tissues, and these genes demonstrated a high diagnostic value for HPV-negative HNSCC. PSMC2, ORC5 and KRTDAP may therefore serve as novel and specific biomarkers for HPV-negative HNSCC, potentially improving the diagnosis and treatment of patients with HPV-negative HNSCC.

Dissecting the Genetic and Etiological Causes of Primary Microcephaly

Autosomal recessive primary microcephaly (MCPH; “small head syndrome”) is a rare, heterogeneous disease arising from the decreased production of neurons during brain development. As of August 2020, the Online Mendelian Inheritance in Man (OMIM) database lists 25 genes (involved in molecular processes such as centriole biogenesis, microtubule dynamics, spindle positioning, DNA repair, transcriptional regulation, Wnt signaling, and cell cycle checkpoints) that are implicated in causing MCPH. Many of these 25 genes were only discovered in the last 10 years following advances in exome and genome sequencing that have improved our ability to identify disease-causing variants. Despite these advances, many patients still lack a genetic diagnosis. This demonstrates a need to understand in greater detail the molecular mechanisms and genetics underlying MCPH. Here, we briefly review the molecular functions of each MCPH gene and how their loss disrupts the neurogenesis program, ultimately demonstrating that microcephaly arises from cell cycle dysregulation. We also explore the current issues in the genetic basis and clinical presentation of MCPH as additional avenues of improving gene/variant prioritization. Ultimately, we illustrate that the detailed exploration of the etiology and inheritance of MCPH improves the predictive power in identifying previously unknown MCPH candidates and diagnosing microcephalic patients.