Altered Pre-mRNA Splicing Caused by a Novel Intronic Mutation c.1443+5G>A in the Dihydropyrimidinase (DPYS) Gene

Dihydropyrimidinase (DHP) deficiency is an autosomal recessive disease caused by mutations in the DPYS gene. Patients present with highly elevated levels of dihydrouracil and dihydrothymine in their urine, blood and cerebrospinal fluid. The analysis of the effect of mutations in DPYS on pre-mRNA splicing is hampered by the fact that DHP is primarily expressed in liver and kidney cells. The minigene approach can detect mRNA splicing aberrations using cells that do not express the endogenous mRNA. We have used a minigene-based approach to analyze the effects of a presumptive pre-mRNA splicing mutation in two newly identified Chinese pediatric patients with DHP deficiency. Mutation analysis of DPYS showed that both patients were compound heterozygous for a novel intronic mutation c.1443+5G>A in intron 8 and a previously described missense mutation c.1001A>G (p.Q334R) in exon 6. Wild-type and the mutated minigene constructs, containing exons 7, 8 and 9 of DPYS, yielded different splicing products after expression in HEK293 cells. The c.1443+5G>A mutation resulted in altered pre-mRNA splicing of the DPYS minigene construct with full skipping of exon 8. Analysis of the DHP crystal structure showed that the deletion of exon 8 severely affects folding, stability and homooligomerization of the enzyme as well as disruption of the catalytic site. Thus, the analysis suggests that the c.1443+5G>A mutation results in aberrant splicing of the pre-mRNA encoding DHP, underlying the DHP deficiency in two unrelated Chinese patients.     

Polymorphisms of Leptin-b Gene Associated with Growth Traits in Orange-Spotted Grouper (Epinephelus coioides)

In mammals, leptin has been demonstrated to perform important roles in many physiological activities and to influence development, growth, metabolism and reproduction. However, in fish, its function is still unclear. Duplicate leptin genes, leptin-a and leptin-b, have been identified in the orange-spotted grouper. In the present study, the polymorphisms in the leptin-b gene of the orange-spotted grouper were detected, and the relation between these polymorphisms and 12 growth traits were analyzed. Six polymorphisms (including 3 single nucleotide polymorphisms (c.14G>A, c.93A>G, c.149G>A) in exon 1, 2 SNPs (c.181A>G, c.193G>A) in intron 1, and 1 SNP (c.360C>T) in exon 2) were identified and genotyped from 200 different individuals. The results revealed that the SNP c.149G>A was significantly associated with growth traits, that the heterozygous mutation genotype GA having negative effects on growth traits. However, the other five SNPs (c.14G>A, c.93A>G, c.181A>G, c.193G>A, c.360C>T) did not show significant associations with all the growth traits. Compared with our findings in leptin-a gene, the results suggested that the leptin-a hormone has more important physiological effects in fish bodies than the leptin-b type. Moreover, leptin genes were supposed to be one class of major candidate genes of regulating growth traits in the orange-spotted grouper.     

Mechanistic Insights of Aberrant Splicing with Splicing Factor Mutations Found in Myelodysplastic Syndromes

Pre-mRNA splicing is an essential process for gene expression in higher eukaryotes, which requires a high order of accuracy. Mutations in splicing factors or regulatory elements in pre-mRNAs often result in many human diseases. Myelodysplastic syndrome (MDS) is a heterogeneous group of chronic myeloid neoplasms characterized by many symptoms and a high risk of progression to acute myeloid leukemia. Recent findings indicate that mutations in splicing factors represent a novel class of driver mutations in human cancers and affect about 50% of Myelodysplastic syndrome (MDS) patients. Somatic mutations in MDS patients are frequently found in genes SF3B1, SRSF2, U2AF1, and ZRSR2. Interestingly, they are involved in the recognition of 3′ splice sites and exons. It has been reported that mutations in these splicing regulators result in aberrant splicing of many genes. In this review article, we first describe molecular mechanism of pre-mRNA splicing as an introduction and mainly focus on those four splicing factors to describe their mutations and their associated aberrant splicing patterns.     

Regulation of Survival Motor Neuron Gene Expression by Calcium Signaling

Spinal muscular atrophy (SMA) is caused by homozygous survival of motor neurons 1 (SMN1) gene deletion, leaving a duplicate gene, SMN2, as the sole source of SMN protein. However, a defect in SMN2 splicing, involving exon 7 skipping, results in a low level of functional SMN protein. Therefore, the upregulation of SMN protein expression from the SMN2 gene is generally considered to be one of the best therapeutic strategies to treat SMA. Most of the SMA drug discovery is based on synthetic compounds, and very few natural compounds have been explored thus far. Here, we performed an unbiased mechanism-independent and image-based screen of a library of microbial metabolites in SMA fibroblasts using an SMN-specific immunoassay. In doing so, we identified brefeldin A (BFA), a well-known inhibitor of ER-Golgi protein trafficking, as a strong inducer of SMN protein. The profound increase in SMN protein was attributed to, in part, the rescue of the SMN2 pre-mRNA splicing defect. Intriguingly, BFA increased the intracellular calcium concentration, and the BFA-induced exon 7 inclusion of SMN2 splicing, was abrogated by the depletion of intracellular calcium and by the pharmacological inhibition of calcium/calmodulin-dependent kinases (CaMKs). Moreover, BFA considerably reduced the expression of Tra2-β and SRSF9 proteins in SMA fibroblasts and enhanced the binding of PSF and hnRNP M to an exonic splicing enhancer (ESE) of exon 7. Together, our results demonstrate a significant role for calcium and its signaling on the regulation of SMN splicing, probably through modulating the expression/activity of splicing factors.     

Polymorphism of the Flap Endonuclease 1 Gene in Keratoconus and Fuchs Endothelial Corneal Dystrophy

Oxidative stress is implicated in the pathogenesis of many diseases, including serious ocular diseases, keratoconus (KC) and Fuchs endothelial corneal dystrophy (FECD). Flap endonuclease 1 (FEN1) plays an important role in the repair of oxidative DNA damage in the base excision repair pathway. We determined the association between two single nucleotide polymorphisms (SNPs), c.–441G>A (rs174538) and g.61564299G>T (rs4246215), in the FEN1 gene and the occurrence of KC and FECD. This study involved 279 patients with KC, 225 patients with FECD and 322 control individuals. Polymerase chain reaction (PCR) and length polymorphism restriction fragment analysis (RFLP) were applied. The T/T genotype of the g.61564299G>T polymorphism was associated with an increased occurrence of KC and FECD. There was no association between the c.–441G>A polymorphism and either disease. However, the GG haplotype of both polymorphisms was observed more frequently and the GT haplotype less frequently in the KC group than the control. The AG haplotype was associated with increased FECD occurrence. Our findings suggest that the g.61564299G>T and c.–441G>A polymorphisms in the FEN1 gene may modulate the risk of keratoconus and Fuchs endothelial corneal dystrophy.     

The GAUGAA Motif Is Responsible for the Binding between circSMARCA5 and SRSF1 and Related Downstream Effects on Glioblastoma Multiforme Cell Migration and Angiogenic Potential

Circular RNAs (circRNAs) are a large class of RNAs with regulatory functions within cells. We recently showed that circSMARCA5 is a tumor suppressor in glioblastoma multiforme (GBM) and acts as a decoy for Serine and Arginine Rich Splicing Factor 1 (SRSF1) through six predicted binding sites (BSs). Here we characterized RNA motifs functionally involved in the interaction between circSMARCA5 and SRSF1. Three different circSMARCA5 molecules (Mut1, Mut2, Mut3), each mutated in two predicted SRSF1 BSs at once, were obtained through PCR-based replacement of wild-type (WT) BS sequences and cloned in three independent pcDNA3 vectors. Mut1 significantly decreased its capability to interact with SRSF1 as compared to WT, based on the RNA immunoprecipitation assay. In silico analysis through the “Find Individual Motif Occurrences” (FIMO) algorithm showed GAUGAA as an experimentally validated SRSF1 binding motif significantly overrepresented within both predicted SRSF1 BSs mutated in Mut1 (q-value = 0.0011). U87MG and CAS-1, transfected with Mut1, significantly increased their migration with respect to controls transfected with WT, as revealed by the cell exclusion zone assay. Immortalized human brain microvascular endothelial cells (IM-HBMEC) exposed to conditioned medium (CM) harvested from U87MG and CAS-1 transfected with Mut1 significantly sprouted more than those treated with CM harvested from U87MG and CAS-1 transfected with WT, as shown by the tube formation assay. qRT-PCR showed that the intracellular pro- to anti-angiogenic Vascular Endothelial Growth Factor A (VEGFA) mRNA isoform ratio and the amount of total VEGFA mRNA secreted in CM significantly increased in Mut1-transfected CAS-1 as compared to controls transfected with WT. Our data suggest that GAUGAA is the RNA motif responsible for the interaction between circSMARCA5 and SRSF1 as well as for the circSMARCA5-mediated control of GBM cell migration and angiogenic potential.     

Altered HIV-1 mRNA Splicing Due to Drug-Resistance-Associated Mutations in Exon 2/2b

The underlying molecular mechanism and their general effect on the replication capacity of HIV 1 drug-resistance-associated mutations is often poorly understood. To elucidate the effect of two such mutations located in a region with a high density of spicing regulatory elements on the HIV-1-splicing outcome, bioinformatic predictions were combined with transfection and infection experiments. Results show that the previously described R263K drug-resistance-associated integrase mutation has additionally a severe effect on the ESE2b splicing regulatory element (SRE) in exon 2b, which causes loss of SD2b recognition. This was confirmed by an R263R silent mutation with a similar predicted effect on the exon 2b SRE. In contrast, a V260I mutation and its silent counterpart with a lower effect on ESS2b did not exhibit any differences in the splicing pattern. Since HIV-1 highly relies on a balanced splicing reaction, changes in the splicing outcome can contribute to changes in viral replication and might add to the effect of escape mutations toward antiviral drugs. Thus, a classification of mutations purely addressing proteins is insufficient.     

A Novel Genetic Variant in the WFS1 Gene in a Patient with Partial Uniparental Mero-Isodisomy of Chromosome 4

Wolfram syndrome is a rare autosomal recessive disorder characterized by optic atrophy and diabetes mellitus. Wolfram syndrome type 1 (WFS1) is caused by bi-allelic pathogenic variations in the wolframin gene. We described the first case of WFS1 due to a maternal inherited mutation with uniparental mero-isodisomy of chromosome 4. Diabetes mellitus was diagnosed at 11 years of age, with negative anti-beta cells antibodies. Blood glucose control was optimal with low insulin requirement. No pathogenic variations in the most frequent gene causative of maturity-onset diabetes of the young subtypes were detected. At 17.8 years old, a rapid reduction in visual acuity occurred. Genetic testing revealed the novel homozygous variant c.1369A>G; p.Arg457Gly in the exon 8 of wolframin gene. It was detected in a heterozygous state only in the mother while the father showed a wild type sequence. In silico disease causing predictions performed by Polyphen2 classified it as “likely damaging”, while Mutation Tester and Sift suggested it was “polymorphism” and “tolerated”, respectively. High resolution SNP-array analysis was suggestive of segmental uniparental disomy on chromosome 4. In conclusion, to the best of our knowledge, we describe the first patient with partial uniparental mero-isodisomy of chromosome 4 carrying a novel mutation in the wolframin gene. The clinical phenotype observed in the patient and the analysis performed suggest that the genetic variant detected is pathogenetic.     

Alternative Splicing Regulation of Low-Frequency Genetic Variants in Exon 2 of TREM2 in Alzheimer’s Disease by Splicing-Based Aggregation

TREM2 is among the most well-known Alzheimer’s disease (AD) risk genes; however, the functional roles of its AD-associated variants remain to be elucidated, and most known risk alleles are low-frequency variants whose investigation is challenging. Here, we utilized a splicing-guided aggregation method in which multiple low-frequency TREM2 variants were bundled together to investigate the functional impact of those variants on alternative splicing in AD. We analyzed whole genome sequencing (WGS) and RNA-seq data generated from cognitively normal elderly controls (CN) and AD patients in two independent cohorts, representing three regions in the frontal lobe of the human brain: the dorsolateral prefrontal cortex (CN = 213 and AD = 376), frontal pole (CN = 72 and AD = 175), and inferior frontal (CN = 63 and AD = 157). We observed an exon skipping event in the second exon of TREM2, with that exon tending to be more frequently skipped (p = 0.0012) in individuals having at least one low-frequency variant that caused loss-of-function for a splicing regulatory element. In addition, genes differentially expressed between AD patients with high vs. low skipping of the second exon (i.e., loss of a TREM2 functional domain) were significantly enriched in immune-related pathways. Our splicing-guided aggregation method thus provides new insight into the regulation of alternative splicing of the second exon of TREM2 by low-frequency variants and could be a useful tool for further exploring the potential molecular mechanisms of multiple, disease-associated, low-frequency variants.     

A Novel Homozygous Variant in DYSF Gene Is Associated with Autosomal Recessive Limb Girdle Muscular Dystrophy R2/2B

Mutations in the DYSF gene, encoding dysferlin, are responsible for Limb Girdle Muscular Dystrophy type R2/2B (LGMDR2/2B), Miyoshi myopathy (MM), and Distal Myopathy with Anterior Tibialis onset (MDAT). The size of the gene and the reported inter and intra familial phenotypic variability make early diagnosis difficult. Genetic analysis was conducted using Next Gene Sequencing (NGS), with a panel of 40 Muscular Dystrophies associated genes we designed. In the present study, we report a new missense variant c.5033G>A, p.Cys1678Tyr ({"type":"entrez-nucleotide","attrs":{"text":"NM_003494","term_id":"1675063869"}}NM_003494) in the exon 45 of DYSF gene related to Limb Girdle Muscular Dystrophy type R2/2B in a 57-year-old patient affected with LGMD from a consanguineous family of south Italy. Both healthy parents carried this variant in heterozygosity. Genetic analysis extended to two moderately affected sisters of the proband, showed the presence of the variant c.5033G>A in both in homozygosity. These data indicate a probable pathological role of the variant c.5033G>A never reported before in the onset of LGMDR2/2B, pointing at the NGS as powerful tool for identifying LGMD subtypes. Moreover, the collection and the networking of genetic data will increase power of genetic-molecular investigation, the management of at-risk individuals, the development of new therapeutic targets and a personalized medicine.     

Deep Molecular Characterization of Milder Spinal Muscular Atrophy Patients Carrying the c.859G>C Variant in SMN2

Spinal muscular atrophy (SMA) is a severe neuromuscular disorder caused by biallelic loss or pathogenic variants in the SMN1 gene. Copy number and modifier intragenic variants in SMN2, an almost identical paralog gene of SMN1, are known to influence the amount of complete SMN proteins. Therefore, SMN2 is considered the main phenotypic modifier of SMA, although genotype–phenotype correlation is not absolute. We present eleven unrelated SMA patients with milder phenotypes carrying the c.859G>C-positive modifier variant in SMN2. All were studied by a specific NGS method to allow a deep characterization of the entire SMN region. Analysis of two homozygous cases for the variant allowed us to identify a specific haplotype, Smn2-859C.1, in association with c.859G>C. Two other cases with the c.859G>C variant in their two SMN2 copies showed a second haplotype, Smn2-859C.2, in cis with Smn2-859C.1, assembling a more complex allele. We also identified a previously unreported variant in intron 2a exclusively linked to the Smn2-859C.1 haplotype (c.154-1141G>A), further suggesting that this region has been ancestrally conserved. The deep molecular characterization of SMN2 in our cohort highlights the importance of testing c.859G>C, as well as accurately assessing the SMN2 region in SMA patients to gain insight into the complex genotype–phenotype correlations and improve prognostic outcomes.     

Different Role of Tumor Necrosis Factor-α Polymorphism in Non-Hodgkin Lymphomas among Caucasian and Asian Populations: A Meta-Analysis

Tumor necrosis factor-α (TNF-α) is an immunoregulatory cytokine involved in B- and T-cell function, and also plays an important role in inflammation and cancer. TNF-α-308G>A has been associated with constitutively elevated TNF-α expression. Several studies have reported the association between the TNF-α-308G>A polymorphism and non-Hodgkin lymphomas (NHL) risk, however, results are still inconsistent. To solve these conflicts, we conducted the first meta-analysis to assess the effect of TNF-α-308G>A polymorphism on the risk of NHL and various subtypes (additive model) including 10,619 cases and 12,977 controls in Caucasian and Asian populations. Our meta-analysis indicated that TNF-α-308G>A polymorphism is not associated with NHL risk when pooling all studies together (OR = 1.06, 95% CI: 0.92–1.23, p = 0.413). In stratified analyses, we found TNF-α-308A allele was significantly associated with higher risk of NHL, B-cell lymphomas (BCL), T-cell lymphomas (TCL) and diffuse large B-cell lymphomas (DLBCL) in Caucasians (OR = 1.22, 95% CI: 1.06–1.40, p = 0.007; OR = 1.18, 95% CI: 1.03–1.34, p = 0.014; OR = 1.20, 95% CI: 1.01–1.42, p = 0.040; OR = 1.21, 95% CI: 1.11–1.32, p < 0.001, respectively). Interestingly, it was associated with decreased risk of NHL, BCL and DLBCL in Asians (OR = 0.75, 95% CI: 0.66–0.86, p < 0.001; OR = 0.70, 95% CI: 0.52–0.94, p = 0.018; OR = 0.70, 95% CI: 0.57–0.86, p = 0.001). These findings also suggest TNF-α might play a distinct role in pathogenesis of NHL in different populations.     

Quantitative Analysis of the Potency of Equimolar Two-Drug Combinations and Combi-Molecules Involving Kinase Inhibitors In Vitro: The Concept of Balanced Targeting

The median-effect principle proposed by Chou and Talalay is the most effective approach to parameterize interactions between several agents in combination. However, this method cannot be used to evaluate the effectiveness of equimolar drug combinations, which are comparative references for dual-targeting molecular design. Here, using data acquired through the development of “combi-molecules” blocking two kinases (e.g., EGFR-c-Src and EGFR-c-Met), we established potency indices for equimolar and dual-targeted inhibitors. If the fold difference (κ) between the IC50 of the two individual kinase inhibitors was >6, the IC50 of their equimolar combination resembled that of the more potent inhibitor. Hence, the “combi-targeting” of the two kinases was considered “imbalanced” and the combination ineffective. However, if κ ≤ 6, the IC50 of the combination fell below that of each individual drug and the combi-targeting was considered “balanced” and the combination effective. We also showed that combi-molecules should be compared with equimolar combinations only under balanced conditions and propose a new parameter Ω for validating their effectiveness. A multi-targeted drug is effective if Ω < 1, where Ω is defined as the IC50 of the drug divided by that of the corresponding equimolar combination. Our study provides a methodology to determine the in vitro potency of equimolar two-drug combinations as well as combi-/hybrid molecules inhibiting two different kinase targets.     

Heterozygous and Homozygous Variants in SORL1 Gene in Alzheimer’s Disease Patients: Clinical,Neuroimaging and Neuropathological Findings

In the last few years, the SORL1 gene has been strongly implicated in the development of Alzheimer’s disease (AD). We performed whole-exome sequencing on 37 patients with early-onset dementia or family history suggestive of autosomal dominant dementia. Data analysis was based on a custom panel that included 46 genes related to AD and dementia. SORL1 variants were present in a high proportion of patients with candidate variants (15%, 3/20). We expand the clinical manifestations associated with the SORL1 gene by reporting detailed clinical and neuroimaging findings of six unrelated patients with AD and SORL1 mutations. We also present for the first time a patient with the homozygous truncating variant c.364C>T (p.R122*) in SORL1, who also had severe cerebral amyloid angiopathy. Furthermore, we report neuropathological findings and immunochemistry assays from one patient with the splicing variant c.4519+5G>A in the SORL1 gene, in which AD was confirmed by neuropathological examination. Our results highlight the heterogeneity of clinical presentation and familial dementia background of SORL1-associated AD and suggest that SORL1 might be contributing to AD development as a risk factor gene rather than as a major autosomal dominant gene.