Splice Site Research Articles

The Effect of Alternative Splicing Sites on Mirtron Formation and Arm Selection of Precursor microRNAs.

Mirtrons represent a subclass of microRNAs (miRNAs) that rely on the splicing machinery for their maturation. However, the molecular details of this Drosha-independent processing are still not fully understood; as an example, the Microprocessor complex cannot process the mirtronic pre-miRNA from the transcript even if splice site mutations are present. To investigate the influence of alternative splicing sites on mirtron formation, we generated Enhanced Green Fluorescent Protein (EGFP) reporters containing artificial introns to compare the processing of canonical miRNAs and mirtrons. Although mutations of both splice sites generated a complex pattern of alternative transcripts, mirtron formation was always severely affected as opposed to the normal processing of the canonical hsa-mir-33b miRNA. However, we also detected that while its formation was also hindered, the mirtron-derived hsa-mir-877-3p miRNA was less affected by certain mutations than the hsa-mir-877-5p species. By knocking down Drosha, we showed that this phenomenon is not dependent on Microprocessor activity but rather points toward the potential stability difference between the miRNAs from the different arms. Our results indicate that when the major splice sites are mutated, mirtron formation cannot be rescued by nearby alternative splice sites, and stability differences between 5p and 3p species should also be considered for functional studies of mirtrons.

U1 snRNA interactions with deep intronic sequences regulate splicing of multiple exons of spinal muscular atrophy genes.

The U1 small nuclear RNA (snRNA) forms ribonucleoprotein particles (RNPs) such as U1 snRNP and U1-TAF15 snRNP. U1 snRNP is one of the most studied RNPs due to its critical role in pre-mRNA splicing in defining the 5' splice site (5'ss) of every exon through direct interactions with sequences at exon/intron junctions. Recent reports support the role of U1 snRNP in all steps of transcription, namely initiation, elongation, and termination. Functions of U1-TAF15 snRNP are less understood, though it associates with the transcription machinery and may modulate pre-mRNA splicing by interacting with the 5'ss and/or 5'ss-like sequences within the pre-mRNA. An anti-U1 antisense oligonucleotide (ASO) that sequesters the 5' end of U1 snRNA inhibits the functions of U1 snRNP, including transcription and splicing. However, it is not known if the inhibition of U1 snRNP influences post-transcriptional regulation of pre-mRNA splicing through deep intronic sequences. We examined the effect of an anti-U1 ASO that sequesters the 5' end of U1 snRNA on transcription and splicing of all internal exons of the spinal muscular atrophy (SMA) genes, SMN1 and SMN2. Our study was enabled by the employment of a multi-exon-skipping detection assay (MESDA) that discriminates against prematurely terminated transcripts. We employed an SMN2 super minigene to determine if anti-U1 ASO differently affects splicing in the context of truncated introns. We observed substantial skipping of multiple internal exons of SMN1 and SMN2 triggered by anti-U1 treatment. Suggesting a role for U1 snRNP in interacting with deep intronic sequences, early exons of the SMN2 super minigene with truncated introns were resistant to anti-U1 induced skipping. Consistently, overexpression of engineered U1 snRNAs targeting the 5'ss of early SMN1 and SMN2 exons did not prevent exon skipping caused by anti-U1 treatment. Our results uncover a unique role of the U1 snRNA-associated RNPs in splicing regulation executed through deep intronic sequences. Findings are significant for developing novel therapies for SMA based on deep intronic targets.

Modulation of prion protein expression through cryptic splice site manipulation

Lowering expression of prion protein (PrP) is a well-validated therapeutic strategy in prion disease, but additional modalities are urgently needed. In other diseases, small molecules have proven capable of modulating pre-mRNA splicing, sometimes by forcing inclusion of cryptic exons that reduce gene expression. Here, we characterize a cryptic exon located in human PRNP’s sole intron and evaluate its potential to reduce PrP expression through incorporation into the 5’ untranslated region (5’UTR). This exon is homologous to exon 2 in non-primate species, but contains a start codon that would yield an upstream open reading frame (uORF) with a stop codon prior to a splice site if included in PRNP mRNA, potentially downregulating PrP expression through translational repression or nonsense-mediated decay. We establish a minigene transfection system and test a panel of splice site alterations, identifying mutants that reduce PrP expression by as much as 78%. Our findings nominate a new therapeutic target for lowering PrP.

Genetic variation mining of the Chinese mitten crab (Eriocheir sinensis) based on transcriptome data from public databases.

At present, public databases house an extensive repository of transcriptome data, with the volume continuing to grow at an accelerated pace. Utilizing these data effectively is a shared interest within the scientific community. In this study, we introduced a novel strategy that harnesses SNPs and InDels identified from transcriptome data, combined with sample metadata from databases, to effectively screen for molecular markers correlated with traits. We utilized 228 transcriptome datasets of Eriocheir sinensis from the NCBI database and employed the Genome Analysis Toolkit software to identify 96 388 SNPs and 20 645 InDels. Employing the genome-wide association study analysis, in conjunction with the gender information from databases, we identified 3456 sex-biased SNPs and 639 sex-biased InDels. The KOG and KEGG annotations of the sex-biased SNPs and InDels revealed that these genes were primarily involved in the metabolic processes of E. sinensis. Combined with SnpEff annotation and PCR experimental validation, a highly sex-biased SNP located in the Kelch domain containing 4 (Klhdc4) gene, CHR67-6415071, was found to alter the splicing sites of Klhdc4, generating two splice variants, Klhdc4_a and Klhdc4_b. Additionally, Klhdc4 exhibited robust expression across the ovaries, testes, and accessory glands. The sex-biased SNPs and InDels identified in this study are conducive to the development of unisexual cultivation methods for E. sinensis, and the alternative splicing event caused by the sex-biased SNP in Klhdc4 may serve as a potential mechanism for sex regulation in E. sinensis. The analysis strategy employed in this study represents a new direction for the rational exploitation and utilization of transcriptome data in public databases.

Loss of Elp1 in cerebellar granule cell progenitors models ataxia phenotype of Familial Dysautonomia

Familial Dysautonomia (FD) is an autosomal recessive disorder caused by a splice site mutation in the gene ELP1, which disproportionally affects neurons. While classically characterized by deficits in sensory and autonomic neurons, neuronal defects in the central nervous system have also been described. Although ELP1 expression remains high in the normal developing and adult cerebellum, its role in cerebellar development is unknown. To explore the role of Elp1 in the cerebellum, we knocked out Elp1 in cerebellar granule cell progenitors (GCPs) and examined the outcome on animal behavior and cellular composition. We found that GCP-specific conditional knockout of Elp1 (Elp1cKO) resulted in ataxia by 8 weeks of age. Cellular characterization showed that the animals had smaller cerebella with fewer granule cells. This defect was already apparent as early as 7 days after birth, when Elp1cKO animals also had fewer mitotic GCPs and shorter Purkinje dendrites. Through molecular characterization, we found that loss of Elp1 was associated with an increase in apoptotic cell death and cell stress pathways in GCPs. Our study demonstrates the importance of ELP1 in the developing cerebellum, and suggests that loss of Elp1 in the GC lineage may also play a role in the progressive ataxia phenotypes of FD patients.

Novel genotypes and phenotypes in Snijders Blok-Campeau syndrome caused by CHD3 mutations.

Snijders Blok-Campeau syndrome (SNIBCPS) is a rare genetic disorder characterized by facial abnormalities, hypotonia, macrocephaly, and global developmental delay (GDD) caused by mutations in CHD3 gene. There is limited information on SNIBCPS and few studies on its pathogenic gene CHD3. We utilized whole-exome sequencing, in vitro minigene splicing assay analysis, and construction of protein models to validate the suspected pathogenic mutation. In addition, the PubMed database was searched using the keywords "Snijders Blok-Campeau syndrome," "CHD3," or "SNIBCPS" to summarize the gene mutations and clinical phenotypic characteristics of children with SNIBCPS. We identified a non-frameshift variant c.3592_c.3606delGCCAAGAGAAAGATG, a splice site variant c.1708-1G>T, and two missense variants, c. 2954G>C (p.Arg985Pro) and c.3371C>T (p.A1124V), in CHD3 variants with SNIBCPS. Importantly, the c.3592_c.3606delGCCAAGAGAAAGATG, c.1708-1G>T, and c.3371C > T (p.A1124V) loci were not reported, and the children in this study also had phenotypic features of unibrow, transverse palmar creases, tracheal bronchus, and hypomelanosis of Ito (HI). The c.1708-1G>T classical splicing mutation leads to abnormal shearing of mRNA, forming a truncated protein that ultimately affects gene function. Our findings have expanded the spectrum of genetic variants and clinical features in children with SNIBCPS. Splicing analysis of CHD3 is an important method to understand the pathogenesis of spliced cells.

COQ7 splice site variant causing a spastic paraparesis phenotype in siblings

COQ7 splice site variant causing a spastic paraparesis phenotype in siblings

Alternative splicing responses to salt stress in Glycyrrhiza uralensis revealed by global profiling of transcriptome RNA-seq datasets.

Excessive reactive oxygen species stress due to salinity poses a significant threat to the growth of Glycyrrhiza uralensis Fisch. To adapt to salt stress, G. uralensis engages in alternative splicing (AS) to generate a variety of proteins that help it withstand the effects of salt stress. While several studies have investigated the impact of alternative splicing on plants stress responses, the mechanisms by which AS interacts with transcriptional regulation to modulate the salt stress response in G. uralensis remain poorly understood. In this study, we utilized high-throughput RNA sequencing data to perform a comprehensive analysis of AS events at various time points in G. uralensis under salt stress, with exon skipping (SE) being the predominant AS type. KEGG enrichment analysis was performed on the different splicing genes (DSG), and pathways associated with AS were significantly enriched, including RNA transport, mRNA surveillance, and spliceosome. This indicated splicing regulation of genes, resulting in AS events under salt stress conditions. Moreover, plant response to salt stress pathways were also enriched, such as mitogen-activated protein kinase signaling pathway - plant, flavonoid biosynthesis, and oxidative phosphorylation. We focused on four differentially significant genes in the MAPK pathway by AS and qRT-PCR analysis. The alternative splicing type of MPK4 and SnRK2 was skipped exon (SE). ETR2 and RbohD were retained intron (RI) and alternative 5'splice site (A5SS), respectively. The expression levels of isoform1 of these four genes displayed different but significant increases in different tissue sites and salt stress treatment times. These findings suggest that MPK4, SnRK2, ETR2, and RbohD in G. uralensis activate the expression of isoform1, leading to the production of more isoform1 protein and thereby enhancing resistance to salt stress. These findings suggest that salt-responsive AS directly and indirectly governs G. uralensis salt response. Further investigations into AS function and mechanism during abiotic stresses may offer novel references for bolstering plant stress tolerance.

Application of whole exome sequencing technology in the diagnosis of hereditary eye diseases

Objective: To evaluate the application of whole exome sequencing (WES) in the diagnosis of hereditary eye diseases. Methods: A total of 24 patients who came to the Obstetrics and Gynecology Hospital of Fudan University for reproductive genetic counseling from December 2020 to December 2023 with the main complaint of congenital eye disorders were included in this study. All cases had no known infections or exposure to known teratogenic drugs, karyotype and chromosome microarray analysis (CMA) abnormalities. Genomic DNA was extracted from the peripheral blood of the probands and their family members and tested for WES. Among them, three individual WES and 21 Trio WES were performed. Potential pathogenic sites were screened and analyzed by Sanger sequencing. For RPGRIP1:c.1611+26G>A site, minigene vector was constructed and RT-qPCR was performed to detect the effect of mRNA splicing. Results: A total of 24 families were collected in this study, of which 20 yielded positive results, achieving a diagnosis rate of 83.3% (20/24). The results involved 21 genes and identified 30 distinct variants, 19 of which were new variants reported. Prenatal diagnostic analysis of family 3 revealed that the fetus carried a c.6970G>T heterozygous nonsense mutation in the PRPF8 gene. The results of RT-PCR with the minigene vector at the non-classical splice site in family 24 indicated that the transcription product of the mutant plasmid was partially retained 104 bp in intron 12, resulting in a p.Glu538Valfs*12 alteration of the protein. Conclusions: The high detection rate of WES in the diagnosis of hereditary eye diseases further supports the advantages of its application as an important molecular detection tool for determining the etiology of hereditary eye diseases.

Phenotypic and mutational spectrum of 17 Chinese patients with Menkes Disease.

Menkes Disease (MD) is a fatal X-linked recessive disorder caused by mutations in the ATP7A gene. Severe cases typically die before the age of three. Mild MD and occipital horn syndrome are variants of MD characterized by a less severe phenotype and longer survival. This case series aims to validate previous findings, expand the clinical phenotype, identify novel ATP7A mutations of MD patients. Observational data with follow-up were collected from 17 genetically diagnosed Chinese MD patients. All 17 patients exhibited neurological symptoms, including delayed motor milestones (100%) and seizures (58.8%). Unspecific pregnancy or delivery complications occurred in 9 patients (52.9%). The most prevalent connective tissue problems were abnormal hair (76.5%), followed by skeletal and dental abnormalities (52.9%), skin problems (41.2%) and hernia (35.3%). Sensorineural hearing loss (17.6%) was previously unreported. Coronary artery aneurysm and patent foramen ovale (5.9%) were infrequent. One 16-year-old boy carries pathological exon 3-4 deletion, presents novel mild phenotype including short stature and cerebellar ataxia. Out of 13 patients with follow-up (median: 24 months), 7 patients (53.8%) died with median survival of 40 months (range: 21-48 months), 3 patients (23.1%) show severe motor development delay and 2 (15.4%) have refractory epilepsy, only the mild MD patient shows improved cerebellar ataxia. Sixteen ATP7A mutations were identified including 6 small indels (37.5%), 5 nonsense mutations (31.2%), 2 missense mutations (12.5%), 2 exon deletions (12.5%), and 1 splice site mutation (6.25%). Fourteen mutations were novel. Our study further broadens the phenotypic and genotypic spectrums of Menkes disease.

Long-read sequencing identifies an SVA_D retrotransposon insertion deep within the intron of ATP7A as a novel cause of occipital horn syndrome

BackgroundSINE-VNTR-Alu (SVA) retrotransposons move from one genomic location to another in a ‘copy-and-paste’ manner. They continue to move actively and cause monogenic diseases through various mechanisms. Currently, disease-causing SVA retrotransposons...

Relevance and regulation of alternative splicing in plant secondary metabolism: current understanding and future directions.

The secondary metabolism of plants is an essential life process enabling organisms to navigate various stages of plant development and cope with ever-changing environmental stresses. Secondary metabolites, abundantly found in nature, possess significant medicinal value. Among the regulatory mechanisms governing these metabolic processes, alternative splicing stands out as a widely observed post-transcriptional mechanism present in multicellular organisms. It facilitates the generation of multiple mRNA transcripts from a single gene by selecting different splicing sites. Selective splicing events in plants are widely induced by various signals, including external environmental stress and hormone signals. These events ultimately regulate the secondary metabolic processes and the accumulation of essential secondary metabolites in plants by influencing the synthesis of primary metabolites, hormone metabolism, biomass accumulation, and capillary density. Simultaneously, alternative splicing plays a crucial role in enhancing protein diversity and the abundance of the transcriptome. This paper provides a summary of the factors inducing alternative splicing events in plants and systematically describes the progress in regulating alternative splicing with respect to different secondary metabolites, including terpenoid, phenolic compounds, and nitrogen-containing compounds. Such elucidation offers critical foundational insights for understanding the role of alternative splicing in regulating plant metabolism and presents novel avenues and perspectives for bioengineering.

Latent Dirichlet allocation mixture models for nucleotide sequence analysis.

Strings of nucleotides carrying biological information are typically described as sequence motifs represented by weight matrices or consensus sequences. However, many signals in DNA or RNA are recognized by multiple factors in temporal sequence, consist of distinct alternative motifs, or are best described by base composition. Here we apply the latent Dirichlet allocation (LDA) mixture model to nucleotide sequences. Using positions in an alignment of human or Drosophila splice sites as samples, we show that LDA readily identifies motifs, including such elusive cases as the intron branch site. Using whole sequences with positional k-mers as features, LDA can identify sequence subtypes enriched in long vs. short introns. LDA with bulk k-mers can reliably distinguish reading frame and species of origin in coding sequences from humans and Drosophila. We find that LDA is a useful model for describing heterogeneous signals, for assigning individual sequences to subtypes, and for identifying and characterizing sequences that do not fit recognized subtypes. Because LDA topic models are interpretable, they also aid the discovery of new motifs, even those present in a small fraction of samples. In summary, LDA can identify and characterize signals in nucleotide sequences, including candidate regulatory factors involved in biological processes.

Reference-informed prediction of alternative splicing and splicing-altering mutations from sequences.

Alternative splicing plays a crucial role in protein diversity and gene expression regulation in higher eukaryotes, and mutations causing dysregulated splicing underlie a range of genetic diseases. Computational prediction of alternative splicing from genomic sequences not only provides insight into gene-regulatory mechanisms but also helps identify disease-causing mutations and drug targets. However, the current methods for the quantitative prediction of splice site usage still have limited accuracy. Here, we present DeltaSplice, a deep neural network model optimized to learn the impact of mutations on quantitative changes in alternative splicing from the comparative analysis of homologous genes. The model architecture enables DeltaSplice to perform "reference-informed prediction" by incorporating the known splice site usage of a reference gene sequence to improve its prediction on splicing-altering mutations. We benchmarked DeltaSplice and several other state-of-the-art methods on various prediction tasks, including evolutionary sequence divergence on lineage-specific splicing and splicing-altering mutations in human populations and neurodevelopmental disorders, and demonstrated that DeltaSplice outperformed consistently. DeltaSplice predicted ∼15% of splicing quantitative trait loci (sQTLs) in the human brain as causal splicing-altering variants. It also predicted splicing-altering de novo mutations outside the splice sites in a subset of patients affected by autism and other neurodevelopmental disorders (NDDs), including 19 genes with recurrent splicing-altering mutations. Integration of splicing-altering mutations with other types of de novo mutation burdens allowed the prediction of eight novel NDD-risk genes. Our work expanded the capacity of in silico splicing models with potential applications in genetic diagnosis and the development of splicing-based precision medicine.

Association of ACTN4 Gene Mutation with Primary Nephrotic Syndrome in Children in Guangxi Autonomous Region, China.

Objective: To investigate the association between ACTN4 gene mutation and primary nephrotic syndrome (PNS) in children in Guangxi Autonomous Region, China. Methods: The high-throughput sequencing technology was used to sequence ACTN4 gene in 155 children with PNS in Guangxi Autonomous Region in China, with 98 healthy children serving as controls. Twenty-three exon-specific capture probes targeting ACTN4 were designed and used to hybridize with the genomic DNA library. The targeted genomic region DNA fragments were enriched and sequenced. The protein levels of ACTN4 in both case and control groups were quantified using ELISA method. Results: Bioinformatics analysis revealed five unique ACTN4 mutations exclusively in patients with PNS, including c.1516G>A (p.G506S) on one exon in 2 patients, c.1442 + 10G>A at the splice site in 1 patient, c.1649A>G (p.D550G) on exon in 1 patient, c.2191-4G>A at the cleavage site in 2 patients, and c.2315C>T (p.A772V) on one exon in 1 patient. The c.1649A>G (p.D550G) and c.2315C>T (p.A772V) were identified from the same patient. Notably, c.1649A>G (p.D550G) represents a novel mutation in ACTN4. In addition, three other ACTN4 polymorphisms occurred in both case and control groups, including c.162 + 6C>T (1 patient in case group and 2 patients in control group), c.572 + 11G>A (1 patient in case group and 2 patients in control group), and c.2191-5C>T (4 patients in the case group and 3 patients in control group). The serum ACTN4 concentration in the case group was markedly higher, averaging 544.7 ng/mL (range: 264.6-952.6 ng/mL), compared with 241.20 ng/mL (range: 110.75-542.35 ng/mL) in the control group. Conclusion: Five ACTN4 polymorphisms were identified among children with PNS in Guangxi Autonomous Region, China, including the novel mutation c.1649A>G. The lower serum levels of α-actinin-4 in the case group suggest that this protein might play a protective role in PNS.

Frequency of rs1051338 and rs116928232 Variants in Individuals from Northwest Mexico.

LIPA, situated on chromosome 10q23.2-q23.3, encodes the enzyme lysosomal acid lipase (LAL) (EC 3.1.1.13). Genetic alterations in LIPA lead to lysosomal acid lipase deficiency (LALD), an inborn error causing lipid metabolism anomalies and impairing cholesterol and triacylglyceride degradation. Over 40 LIPA variants have been documented, yet this study focuses on just two. The rs1051338 variant (NM_000235:c.46A>C) affects the signal peptide in Exon 2, whereas rs116928232, located in Exon 8, alters the splice site (NM_000235:c.894G>A), impacting lysosomal acid lipase activity. Considering the diverse clinical manifestations of LALD and the rising hepatic steatosis prevalence in Mexican population, mainlydue to diet, these variants were investigated within this demographic to uncover potential contributing factors. This study aimed to reveal the frequency of rs1051338 and rs116928232 among healthy mestizo individuals in Northwest Mexico, marking a significant genetic exploration in this demographic. Three hundred ten healthy mestizo individuals underwent PCR-RFLP analysis for both variants, and Sanger sequencing was performed for variant rs116928232. Bioinformatic analysis was also performed to predict protein changes. Allele frequencies for rs1051338 (FA = 0.39, p value = 0.15) and rs116928232 (FA = 0.0016, p value = 0.49) aligned with reported data, while bioinformatic analysis allowed us to identify the protein alteration observed in both variants; finally, the variants showed no linkage between them (normalized D' = 1.03, p value = 0.56). Allelic frequencies closely matched reported data, and protein structure analysis confirmed variant impacts on LAL enzyme function. Notably, this study marks the first analysis of rs1051338 and rs116928232 in a healthy Mexican mestizo population.

Taxonomy of introns and the evolution of minor introns.

Classification of introns, which is crucial to understanding their evolution and splicing, has historically been binary and has resulted in the naming of major and minor introns that are spliced by their namesake spliceosome. However, a broad range of intron consensus sequences exist, leading us to here reclassify introns as minor, minor-like, hybrid, major-like, majorand non-canonical introns in 263 species across six eukaryotic supergroups. Through intron orthology analysis, we discovered that minor-like introns are a transitory node for intron conversion across evolution. Despite close resemblance of their consensus sequences to minor introns, these introns possess an AG dinucleotide at the -1 and -2 position of the 5' splice site, a salient feature of major introns. Through combined analysis of CoLa-seq, CLIP-seq for major and minor spliceosome components, and RNAseq from samples in which the minor spliceosome is inhibited we found that minor-like introns are also an intermediate class from a splicing mechanism perspective. Importantly, this analysis has provided insight into the sequence elements that have evolved to make minor-like introns amenable to recognition by both minor and major spliceosome components. We hope that this revised intron classification provides a new framework to study intron evolution and splicing.

Alternative splicing of lncRNA LAIR fine-tunes the regulation of neighboring yield-related gene LRK1 expression in rice.

The diversity in alternative splicing of long noncoding RNAs (lncRNAs) poses a challenge for functional annotation of lncRNAs. Moreover, little is known on the effects of alternatively spliced lncRNAs on crop yield. In this study, we cloned nine isoforms resulting from the alternative splicing of the lncRNA LAIR in rice. The LAIR isoforms are generated via alternative 5'/3' splice sites and different combinations of specific introns. All LAIR isoforms activate the expression of the neighboring LRK1 gene and enhance yield-related rice traits. In addition, there are slight differences in the binding ability of LAIR isoforms to the epigenetic modification-related proteins OsMOF and OsWDR5, which affect the enrichment of H4K16ac and H3K4me3 at the LRK1 locus, and consequently fine-tune the regulation of LRK1 expression and yield-related traits. These differences in binding may be caused by polymorphic changes to the RNA secondary structure resulting from alternative splicing. It was also observed that the composition of LAIR isoforms was sensitive to abiotic stress. These findings suggest that the alternative splicing of LAIR leads to the formation of a functional transcript population that precisely regulates yield-related gene expression, which may be relevant for phenotypic polymorphism-based crop breeding under changing environmental conditions.

Geneapp: A web application for visualizing alternative splicing for biomedicine

Alternative Splicing (AS) is an essential mechanism for eukaryotes. However, the consequences of deleting a single exon can be dramatic for the organism and can lead to cancer in humans. Additionally, alternative 5′ and 3′ splice sites, which define the boundaries of exons, also play key roles to human disorders. Therefore, Investigating AS events is crucial for understanding the molecular basis of human diseases and developing therapeutic strategies. Workflow for AS event analysis can be sampling followed by data analysis with bioinformatics to identify the different AS events in the control and case samples, data visualization for curation, and selection of relevant targets for experimental validation. The raw output of the analysis software does not favor the inspection of events by bioinformaticians requiring custom scripts for data visualization. In this work, we propose the Geneapp application with three modules: GeneappScript, GeneappServer, and GeneappExplorer. GeneappScript is a wrapper that assists in identifying AS in samples compared in two different approaches, while GeneappServer integrates data from AS analysis already performed by the user. In GeneappExplorer, the user visualizes the previous dataset by exploring AS events in genes with functional annotation. This targeted screens that Geneapp allows to perform helps in the identification of targets for experimental validation to confirm the hypotheses under study. The Geneapp is freely available for non-commercial use at https://geneapp.net to advance research on AS for bioinformatics.

Molecular genetics, neuroimaging outcomes, and structural analyses of novel and recurrent variants of WDR62 gene in two consanguineous Pakistani families with autosomal recessive primary microcephaly.

Autosomal recessive primary microcephaly (MCPH) is a rare neurodevelopmental and genetically heterogeneous disorder, characterized by small cranium size (> -3 SD below mean) and often results in varying degree of intellectual disability. Thirty genes have been identified for the etiology of this disorder due to its clinical and genetic heterogeneity. Here, we report two consanguineous Pakistani families affected with MCPH exhibiting mutation in WDR62 gene. The investigation approach involved Next Generation Sequencing (NGS) gene panel sequencing coupled with linkage analysis followed by validation of identified variants through automated Sanger sequencing and Barcode-Tagged (BT) sequencing. The molecular genetic analysis revealed one novel splice site variant (NM_001083961.2(WDR62):c.1372-1del) in Family A and one known exonic variant NM_001083961.2(WDR62):c.3936dup (p.Val1313Argfs*18) in Family B. Magnetic Resonance Imaging (MRI) scans were also employed to gain insights into the structural architecture of affected individuals. Neurological assessments showed the reduced gyral and sulcal patterns along with normal corpus callosum in affected individuals harboring novel variant. In silico assessments of the identified variants were conducted using different tools to confirm the pathogenicity of these variants. Through In silico analyses, both variants were identified as disease causing and protein modeling of exonic variant indicates subtle conformational alterations in prophesied protein structure. This study identifies a novel variant (c.1372-1del) and a recurrent pathogenic variant c.3936dup (p.Val1313Argfs*18) in the WDR62 gene among the Pakistani population, expanding the mutation spectrum for MCPH. These findings emphasize the importance of genetic counseling and awareness to reduce consanguinity and address the burden of this disorder.