Functional analysis of synonymous substitutions predicted to affect splicing of the CFTR gene

Moderating the neutralist-selectionist debate: exactly which propositions are we debating, and which arguments are valid?

Half a century after its foundation, the neutral theory of molecular evolution continues to attract controversy. The debate has been hampered by the coexistence of different interpretations of the core proposition of the neutral theory, the 'neutral mutation-random drift' hypothesis. In this review, we trace the origins of these ambiguities and suggest potential solutions. We highlight the difference between the original, the revised and the nearly neutral hypothesis, and re-emphasise that none of them equates to the null hypothesis of strict neutrality. We distinguish the neutral hypothesis of protein evolution, the main focus of the ongoing debate, from the neutral hypotheses of genomic and functional DNA evolution, which for many species are generally accepted. We advocate a further distinction between a narrow and an extended neutral hypothesis (of which the latter posits that random non-conservative amino acid substitutions can cause non-ecological phenotypic divergence), and we discuss the implications for evolutionary biology beyond the domain of molecular evolution. We furthermore point out that the debate has widened from its initial focus on point mutations, and also concerns the fitness effects of large-scale mutations, which can alter the dosage of genes and regulatory sequences. We evaluate the validity of neutralist and selectionist arguments and find that the tested predictions, apart from being sensitive to violation of underlying assumptions, are often derived from the null hypothesis of strict neutrality, or equally consistent with the opposing selectionist hypothesis, except when assuming molecular panselectionism. Our review aims to facilitate a constructive neutralist-selectionist debate, and thereby to contribute to answering a key question of evolutionary biology: what proportions of amino acid and nucleotide substitutions and polymorphisms are adaptive?

Principles and Practical Considerations for the Analysis of Disease-Associated Alternative Splicing Events Using the Gateway Cloning-Based Minigene Vectors pDESTsplice and pSpliceExpress

Splicing is an important RNA processing step. Genetic variations can alter the splicing process and thereby contribute to the development of various diseases. Alterations of the splicing pattern can be examined by gene expression analyses, by computational tools for predicting the effects of genetic variants on splicing, and by splicing reporter minigene assays for studying alternative splicing events under defined conditions. The minigene assay is based on transient transfection of cells with a vector containing a genomic region of interest cloned between two constitutive exons. Cloning can be accomplished by the use of restriction enzymes or by site-specific recombination using Gateway cloning. The vectors pDESTsplice and pSpliceExpress represent two minigene systems based on Gateway cloning, which are available through the Addgene plasmid repository. In this review, we describe the features of these two splicing reporter minigene systems. Moreover, we provide an overview of studies in which determinants of alternative splicing were investigated by using pDESTsplice or pSpliceExpress. The studies were reviewed with regard to the investigated splicing regulatory events and the experimental strategy to construct and perform a splicing reporter minigene assay. We further elaborate on how analyses on the regulation of RNA splicing offer promising prospects for gaining important insights into disease mechanisms.

Evaluation of both exonic and intronic variants for effects on RNA splicing allows for accurate assessment of the effectiveness of precision therapies

Elucidating the functional consequence of molecular defects underlying genetic diseases enables appropriate design of therapeutic options. Treatment of cystic fibrosis (CF) is an exemplar of this paradigm as the development of CFTR modulator therapies has allowed for targeted and effective treatment of individuals harboring specific genetic variants. However, the mechanism of these drugs limits effectiveness to particular classes of variants that allow production of CFTR protein. Thus, assessment of the molecular mechanism of individual variants is imperative for proper assignment of these precision therapies. This is particularly important when considering variants that affect pre-mRNA splicing, thus limiting success of the existing protein-targeted therapies. Variants affecting splicing can occur throughout exons and introns and the complexity of the process of splicing lends itself to a variety of outcomes, both at the RNA and protein levels, further complicating assessment of disease liability and modulator response. To investigate the scope of this challenge, we evaluated splicing and downstream effects of 52 naturally occurring CFTR variants (exonic = 15, intronic = 37). Expression of constructs containing select CFTR intronic sequences and complete CFTR exonic sequences in cell line models allowed for assessment of RNA and protein-level effects on an allele by allele basis. Characterization of primary nasal epithelial cells obtained from individuals harboring splice variants corroborated in vitro data. Notably, we identified exonic variants that result in complete missplicing and thus a lack of modulator response (e.g. c.2908G>A, c.523A>G), as well as intronic variants that respond to modulators due to the presence of residual normally spliced transcript (e.g. c.4242+2T>C, c.3717+40A>G). Overall, our data reveals diverse molecular outcomes amongst both exonic and intronic variants emphasizing the need to delineate RNA, protein, and functional effects of each variant in order to accurately assign precision therapies.

Genetic variants that impact pre-mRNA splicing are a common cause of genetic disease and have varying downstream molecular consequences. As a result, precision therapies that function at the protein level are not always effective for these variants and thus careful assessment is necessary. Here we evaluate RNA-level effects of 52 variants in the cystic fibrosis transmembrane conductance regulator (CFTR) gene and show that study of splicing and its consequences allows for more accurate assignment of precision therapies.

Antisense oligonucleotide-mediated correction of CFTR splicing improves chloride secretion in cystic fibrosis patient-derived bronchial epithelial cells

Cystic fibrosis (CF) is an autosomal recessive disorder caused by mutations in the CF transmembrane conductance regulator (CFTR) gene, encoding an anion channel that conducts chloride and bicarbonate across epithelial membranes. Mutations that disrupt pre-mRNA splicing occur in >15% of CF cases. One common CFTR splicing mutation is CFTR c.3718-2477C>T (3849+10 kb C>T), which creates a new 5′ splice site, resulting in splicing to a cryptic exon with a premature termination codon. Splice-switching antisense oligonucleotides (ASOs) have emerged as an effective therapeutic strategy to block aberrant splicing. We test an ASO targeting the CFTR c.3718-2477C>T mutation and show that it effectively blocks aberrant splicing in primary bronchial epithelial (hBE) cells from CF patients with the mutation. ASO treatment results in long-term improvement in CFTR activity in hBE cells, as demonstrated by a recovery of chloride secretion and apical membrane conductance. We also show that the ASO is more effective at recovering chloride secretion in our assay than ivacaftor, the potentiator treatment currently available to these patients. Our findings demonstrate the utility of ASOs in correcting CFTR expression and channel activity in a manner expected to be therapeutic in patients.

Synonymous variants associated with Alzheimer disease in multiplex families

OBJECTIVE

Synonymous variants can lead to disease; nevertheless, the majority of sequencing studies conducted in Alzheimer disease (AD) only assessed coding variation.

METHODS

To detect synonymous variants modulating AD risk, we conducted a whole-genome sequencing study on 67 Caribbean Hispanic (CH) families multiply affected by AD. Identified disease-associated variants were further assessed in an independent cohort of CHs, expression quantitative trait locus (eQTL) data, brain autopsy data, and functional experiments.

RESULTS

Rare synonymous variants in 4 genes (CDH23, SLC9A3R1, RHBDD2, and ITIH2) segregated with AD status in multiplex families and had a significantly higher frequency in these families compared with reference populations of similar ancestry. In comparison to subjects without dementia, expression of CDH23 (β = 0.53, p = 0.006) and SLC9A3R1 (β = 0.50, p = 0.02) was increased, and expression of RHBDD2 (β = -0.70, p = 0.02) decreased in individuals with AD at death. In line with this finding, increased expression of CDH23 (β = 0.26 ± 0.08, p = 4.9E-4) and decreased expression of RHBDD2 (β = -0.60 ± 0.12, p = 5.5E-7) were related to brain amyloid load (p = 0.0025). SLC9A3R1 expression was associated with burden of TDP43 pathology (β = 0.58 ± 0.17, p = 5.9E-4). Using eQTL data, the CDH23 variant was in linkage disequilibrium with variants modulating CDH23 expression levels (top single nucleotide polymorphism: rs11000035, p = 4.85E-6, D' = 1.0). Using minigene splicing assays, the CDH23 and SLC9A3R1 variants affected splicing efficiency.

CONCLUSIONS

These findings suggest that CDH23, SLC9A3R1, RHBDD2, and possibly ITIH2, which are involved in synaptic function, the glutamatergic system, and innate immunity, contribute to AD etiology. In addition, this study supports the notion that synonymous variants contribute to AD risk and that comprehensive scrutinization of this type of genetic variation is warranted and critical.

CAGI experiments: Modeling sequence variant impact on gene splicing using predictions from computational tools

Improving predictions of phenotypic consequences for genomic variants is part of ongoing efforts in the scientific community to gain meaningful insights into genomic function. Within the framework of the critical assessment of genome interpretation experiments, we participated in the Vex-seq challenge, which required predicting the change in the percent spliced in measure (ΔΨ) for 58 exons caused by more than 1,000 genomic variants. Experimentally determined through the Vex-seq assay, the Ψ quantifies the fraction of reads that include an exon of interest. Predicting the change in Ψ associated with specific genomic variants implies determining the sequence changes relevant for splicing regulators, such as splicing enhancers and silencers. Here we took advantage of two computational tools, SplicePort and SPANR, that incorporate relevant sequence features in their models of splice sites and exon-inclusion level, respectively. Specifically, we used the SplicePort and SPANR outputs to build mathematical models of the experimental data obtained for the variants in the training set, which we then used to predict the ΔΨ associated with the mutations in the test set. We show that the sequence changes captured by these computational tools provide a reasonable foundation for modeling the impact on splicing associated with genomic variants.

regSNPs-splicing: a tool for prioritizing synonymous single-nucleotide substitution

While synonymous single-nucleotide variants (sSNVs) have largely been unstudied, since they do not alter protein sequence, mounting evidence suggests that they may affect RNA conformation, splicing, and the stability of nascent-mRNAs to promote various diseases. Accurately prioritizing deleterious sSNVs from a pool of neutral ones can significantly improve our ability of selecting functional genetic variants identified from various genome-sequencing projects, and, therefore, advance our understanding of disease etiology. In this study, we develop a computational algorithm to prioritize sSNVs based on their impact on mRNA splicing and protein function. In addition to genomic features that potentially affect splicing regulation, our proposed algorithm also includes dozens structural features that characterize the functions of alternatively spliced exons on protein function. Our systematical evaluation on thousands of sSNVs suggests that several structural features, including intrinsic disorder protein scores, solvent accessible surface areas, protein secondary structures, and known and predicted protein family domains, show significant differences between disease-causing and neutral sSNVs. Our result suggests that the protein structure features offer an added dimension of information while distinguishing disease-causing and neutral synonymous variants. The inclusion of structural features increases the predictive accuracy for functional sSNV prioritization.

The functional relevance of somatic synonymous mutations in melanoma and other cancers

Recent technological advances in sequencing have flooded the field of cancer research with knowledge about somatic mutations for many different cancer types. Most cancer genomics studies focus on mutations that alter the amino acid sequence, ignoring the potential impact of synonymous mutations. However, accumulating experimental evidence has demonstrated clear consequences for gene function, leading to a widespread recognition of the functional role of synonymous mutations and their causal connection to various diseases. Here, we review the evidence supporting the direct impact of synonymous mutations on gene function via gene splicing; mRNA stability, folding, and translation; protein folding; and miRNA-based regulation of expression. These results highlight the functional contribution of synonymous mutations to oncogenesis and the need to further investigate their detection and prioritization for experimental assessment.

Analysis of long-range interactions in primary human cells identifies cooperative CFTR regulatory elements

A mechanism by which control DNA elements regulate transcription over large linear genomic distances is by achieving close physical proximity with genes, and looping of the intervening chromatin paths. Alterations of such regulatory 'chromatin looping' systems are likely to play a critical role in human genetic disease at large. Here, we studied the spatial organization of a ≈790 kb locus encompassing the cystic fibrosis transmembrane conductance regulator (CFTR) gene. Dysregulation of CFTR is responsible for cystic fibrosis, which is the most common lethal genetic disorder in Caucasian populations. CFTR is a relatively large gene of 189 kb with a rather complex tissue-specific and temporal expression profile. We used chromatin conformation at the CFTR locus to identify new DNA sequences that regulate its transcription. By comparing 5C chromatin interaction maps of the CFTR locus in expressing and non-expressing human primary cells, we identified several new contact points between the CFTR promoter and its surroundings, in addition to regions featuring previously described regulatory elements. We demonstrate that two of these novel interacting regions cooperatively increase CFTR expression, and suggest that the new enhancer elements located on either side of the gene are brought together through chromatin looping via CTCF.

Conservation of CFTR codon frequency through primates suggests synonymous mutations could have a functional effect

Cystic fibrosis is an inherited chronic disease that affects the lungs and digestive system, with a prevalence of about 1:3000 people. Cystic fibrosis is caused by mutations in CFTR gene, which lead to a defective function of the chloride channel, the cystic fibrosis transmembrane conductance regulator (CFTR). Up-to-date, more than 1900 mutations have been reported in CFTR. However for an important proportion of them, their functional effects and the relation to disease are still not understood. Many of these mutations are silent (or synonymous), namely they do not alter the encoded amino acid. These synonymous mutations have been considered as neutral to protein function. However, more recent evidence in bacterial and human proteins has put this concept under revision. With the aim of understanding possible functional effects of synonymous mutations in CFTR, we analyzed human and primates CFTR codon usage and divergence patterns. We report the presence of regions enriched in rare and frequent codons. This spatial pattern of codon preferences is conserved in primates, but this cannot be explained by sequence conservation alone. In sum, the results presented herein suggest a functional implication of these regions of the gene that may be maintained by purifying selection acting to preserve a particular codon usage pattern along the sequence. Overall these results support the idea that several synonymous mutations in CFTR may have functional importance, and could be involved in the disease.

A rigorous approach for selection of optimal variant sets for carrier screening with demonstration of clinical utility

Carrier screening for certain diseases is recommended by major medical and Ashkenazi Jewish (AJ) societies. Most carrier screening panels test only for common, ethnic-specific variants. However, with formerly isolated ethnic groups becoming increasingly intermixed, this approach is becoming inadequate. Our objective was to develop a rigorous process to curate all variants, for relevant genes, into a database and then apply stringent clinical validity classification criteria to each in order to retain only those with clear evidence for pathogenicity. The resulting variant set, in conjunction with next-generation DNA sequencing (NGS), then affords the capability for an ethnically diverse, comprehensive, highly specific carrier-screening assay. The clinical utility of our approach was demonstrated by screening a pan-ethnic population of 22,864 individuals for Bloom syndrome carrier status using a BLM variant panel comprised of 50 pathogenic variants. In addition to carriers of the common AJ founder variant, we identified 57 carriers of other pathogenic BLM variants. All variants reported had previously been curated and their clinical validity documented, or were of a type that met our stringent, preassigned validity criteria. Thus, it was possible to confidently report an increased number of Bloom’s syndrome carriers compared to traditional, ethnicity-based screening, while not reducing the specificity of the screening due to reporting variants of unknown clinical significance.

Experimental assessment of splicing variants using expression minigenes and comparison with in silico predictions

Assessment of the functional consequences of variants near splice sites is a major challenge in the diagnostic laboratory. To address this issue, we created expression minigenes (EMGs) to determine the RNA and protein products generated by splice site variants (n = 10) implicated in cystic fibrosis (CF). Experimental results were compared with the splicing predictions of eight in silico tools. EMGs containing the full-length Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) coding sequence and flanking intron sequences generated wild-type transcript and fully processed protein in Human Embryonic Kidney (HEK293) and CF bronchial epithelial (CFBE41o-) cells. Quantification of variant induced aberrant mRNA isoforms was concordant using fragment analysis and pyrosequencing. The splicing patterns of c.1585-1G>A and c.2657+5G>A were comparable to those reported in primary cells from individuals bearing these variants. Bioinformatics predictions were consistent with experimental results for 9/10 variants (MES), 8/10 variants (NNSplice), and 7/10 variants (SSAT and Sroogle). Programs that estimate the consequences of mis-splicing predicted 11/16 (HSF and ASSEDA) and 10/16 (Fsplice and SplicePort) experimentally observed mRNA isoforms. EMGs provide a robust experimental approach for clinical interpretation of splice site variants and refinement of in silico tools.

Functional assessment of potential splice site variants in arrhythmogenic right ventricular dysplasia/cardiomyopathy

BACKGROUND

Interpretation of genetic screening results in arrhythmogenic right ventricular dysplasia/cardiomyopathy (ARVD/C) often is difficult. Pathogenicity of variants with uncertain clinical significance may be predicted by software algorithms. However, functional assessment can unambiguously demonstrate the effect of such variants.

OBJECTIVE

The purpose of this study was to perform functional analysis of potential splice site variants in ARVD/C patients.

METHODS

Nine variants in desmosomal (PKP2, JUP, DSG2, DSC2) genes with potential RNA splicing effect were analyzed. The variants were found in patients who fulfilled 2010 ARVD/C Task Force Criteria (n = 7) or had suspected ARVD/C (n = 2). Total RNA was isolated from fresh blood samples and subjected to reverse transcriptase polymerase chain reaction.

RESULTS

An effect on splicing was predicted by software algorithms for all variants. Of the 9 variants, 5 were intronic and 4 exonic. RNA analysis showed a functional effect on mRNA splicing by exon skipping, generation of new splice sites, or activation of cryptic sites in 6 variants. All 5 intronic variants tested severely impaired splicing. Only 1 of 4 exonic potential splice site variants was shown to have a deleterious effect on splicing. The remaining 3 exonic variants had no detectable effect on splicing, and heterozygous presence in mRNA confirmed biallelic expression.

CONCLUSION

Six variants of uncertain clinical significance in the PKP2, JUP, and DSG2 genes showed a deleterious effect on mRNA splicing, indicating these are ARVD/C-related pathogenic splice site mutations. These results highlight the importance of functional assessment of potential splice site variants to improve patient care and facilitate cascade screening.

Defining the disease liability of variants in the cystic fibrosis transmembrane conductance regulator gene

Allelic heterogeneity in disease-causing genes presents a substantial challenge to the translation of genomic variation to clinical practice. Few of the almost 2,000 variants in the cystic fibrosis transmembrane conductance regulator (CFTR) gene have empirical evidence that they cause cystic fibrosis. To address this gap, we collected both genotype and phenotype data for 39,696 cystic fibrosis patients in registries and clinics in North America and Europe. Among these patients, 159 CFTR variants had an allele frequency of ≥0.01%. These variants were evaluated for both clinical severity and functional consequence with 127 (80%) meeting both clinical and functional criteria consistent with disease. Assessment of disease penetrance in 2,188 fathers of cystic fibrosis patients enabled assignment of 12 of the remaining 32 variants as neutral while the other 20 variants remained indeterminate. This study illustrates that sourcing data directly from well-phenotyped subjects can address the gap in our ability to interpret clinically-relevant genomic variation.