Simultaneous assessment of the effects of exonic mutations on RNA splicing and protein functions

Code inside the codon: The role of synonymous mutations in regulating splicing machinery and its impact on disease

In eukaryotes, precise pre-mRNA processing, including alternative splicing, is essential to carry out the intricate protein translation process. Both point mutations (that alter the translated protein sequence) and synonymous mutations (that do not alter the translated protein sequence) are capable of affecting the splicing process. Synonymous mutations are known to affect gene expression via altering mRNA stability, mRNA secondary structure, splicing processes, and translational kinetics. In higher eukaryotes, precise splicing is regulated by three weakly conserved cis-elements, 5' and 3' splice sites and the branch site. Many other cis-acting elements (exonic/intronic splicing enhancers and silencers) and trans-acting splicing factors (serine and arginine-rich proteins and heterogeneous nuclear ribonucleoproteins) have also been found to enhance or suppress the splicing process. The appearance of synonymous mutations in cis-acting elements can alter the splicing process by changing the binding pattern of splicing factors to exonic splicing enhancers or silencer motifs. This results in exon skipping, intron retention, and various other forms of alternative splicing, eventually leading to the emergence of a wide range of diseases. The focus of this review is to elucidate the role of synonymous mutations and their impact on abnormal splicing mechanisms. Further, this study highlights the function of synonymous mutation in mediating abnormal splicing in cancer and development of X-linked, and autosomal inherited diseases.

CUBAP: an interactive web portal for analyzing codon usage biases across populations

Synonymous codon usage significantly impacts translational and transcriptional efficiency, gene expression, the secondary structure of both mRNA and proteins, and has been implicated in various diseases. However, population-specific differences in codon usage biases remain largely unexplored. Here, we present a web server, https://cubap.byu.edu, to facilitate analyses of codon usage biases across populations (CUBAP). Using the 1000 Genomes Project, we calculated and visually depict population-specific differences in codon frequencies, codon aversion, identical codon pairing, co-tRNA codon pairing, ramp sequences, and nucleotide composition in 17,634 genes. We found that codon pairing significantly differs between populations in 35.8% of genes, allowing us to successfully predict the place of origin for African and East Asian individuals with 98.8% and 100% accuracy, respectively. We also used CUBAP to identify a significant bias toward decreased CTG pairing in the immunity related GTPase M (IRGM) gene in East Asian and African populations, which may contribute to the decreased association of rs10065172 with Crohn's disease in those populations. CUBAP facilitates in-depth gene-specific and codon-specific visualization that will aid in analyzing candidate genes identified in genome-wide association studies, identifying functional implications of synonymous variants, predicting population-specific impacts of synonymous variants and categorizing genetic biases unique to certain populations.

Comparative analysis of functional assay evidence use by ClinGen Variant Curation Expert Panels

BACKGROUND

The 2015 American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) guidelines for clinical sequence variant interpretation state that "well-established" functional studies can be used as evidence in variant classification. These guidelines articulated key attributes of functional data, including that assays should reflect the biological environment and be analytically sound; however, details of how to evaluate these attributes were left to expert judgment. The Clinical Genome Resource (ClinGen) designates Variant Curation Expert Panels (VCEPs) in specific disease areas to make gene-centric specifications to the ACMG/AMP guidelines, including more specific definitions of appropriate functional assays. We set out to evaluate the existing VCEP guidelines for functional assays.

METHODS

We evaluated the functional criteria (PS3/BS3) of six VCEPs (CDH1, Hearing Loss, Inherited Cardiomyopathy-MYH7, PAH, PTEN, RASopathy). We then established criteria for evaluating functional studies based on disease mechanism, general class of assay, and the characteristics of specific assay instances described in the primary literature. Using these criteria, we extensively curated assay instances cited by each VCEP in their pilot variant classification to analyze VCEP recommendations and their use in the interpretation of functional studies.

RESULTS

Unsurprisingly, our analysis highlighted the breadth of VCEP-approved assays, reflecting the diversity of disease mechanisms among VCEPs. We also noted substantial variability between VCEPs in the method used to select these assays and in the approach used to specify strength modifications, as well as differences in suggested validation parameters. Importantly, we observed discrepancies between the parameters VCEPs specified as required for approved assay instances and the fulfillment of these requirements in the individual assays cited in pilot variant interpretation.

CONCLUSIONS

Interpretation of the intricacies of functional assays often requires expert-level knowledge of the gene and disease, and current VCEP recommendations for functional assay evidence are a useful tool to improve the accessibility of functional data by providing a starting point for curators to identify approved functional assays and key metrics. However, our analysis suggests that further guidance is needed to standardize this process and ensure consistency in the application of functional evidence.

Decoding mechanisms by which silent codon changes influence protein biogenesis and function

SCOPE

Synonymous codon usage has been a focus of investigation since the discovery of the genetic code and its redundancy. The occurrences of synonymous codons vary between species and within genes of the same genome, known as codon usage bias. Today, bioinformatics and experimental data allow us to compose a global view of the mechanisms by which the redundancy of the genetic code contributes to the complexity of biological systems from affecting survival in prokaryotes, to fine tuning the structure and function of proteins in higher eukaryotes. Studies analyzing the consequences of synonymous codon changes in different organisms have revealed that they impact nucleic acid stability, protein levels, structure and function without altering amino acid sequence. As such, synonymous mutations inevitably contribute to the pathogenesis of complex human diseases. Yet, fundamental questions remain unresolved regarding the impact of silent mutations in human disorders. In the present review we describe developments in this area concentrating on mechanisms by which synonymous mutations may affect protein function and human health.

PURPOSE

This synopsis illustrates the significance of synonymous mutations in disease pathogenesis. We review the different steps of gene expression affected by silent mutations, and assess the benefits and possible harmful effects of codon optimization applied in the development of therapeutic biologics.

PHYSIOLOGICAL AND MEDICAL RELEVANCE

Understanding mechanisms by which synonymous mutations contribute to complex diseases such as cancer, neurodegeneration and genetic disorders, including the limitations of codon-optimized biologics, provides insight concerning interpretation of silent variants and future molecular therapies.

Exonic mutations in the SLC12A3 gene cause exon skipping and premature termination in Gitelman syndrome

A variety of genetic backgrounds cause the loss of function of thiazide-sensitive sodium chloride cotransporter, encoded by SLC12A3, responsible for the phenotypes in Gitelman syndrome. Recently, the phenomenon of exon skipping, in which exonic mutations result in abnormal splicing, has been associated with various diseases. Specifically, mutations in exonic splicing enhancer (ESE) sequences can promote exon skipping. Here, we used a bioinformatics program to analyze 88 missense mutations in the SLC12A3 gene and identify candidate mutations that may induce exon skipping. The three candidate mutations that reduced ESE scores the most were further investigated by minigene assay, and two (p.A356V and p.M672I) caused abnormal splicing in vitro. Furthermore, we identified the p.M672I (c.2016G>A) mutation in a patient with Gitelman syndrome and found that this single nucleotide mutation causes exclusion of exon 16 in the SLC12A3 mRNA transcript. Functional analyses revealed that the protein encoded by the aberrant SLC12A3 transcript does not transport sodium. These results suggest that aberrant exon skipping is one previously unrecognized mechanism by which missense mutations in SLC12A3 can lead to Gitelman syndrome.

The mutation spectrum of the phenylalanine hydroxylase (PAH) gene and associated haplotypes reveal ethnic heterogeneity in the Taiwanese population

Phenylalanine hydroxylase (PAH) deficiency is responsible for most cases of phenylketonuria (PKU). In this study of the PAH mutation spectrum in the Taiwanese population, 139 alleles were identified including 34 different mutations. The V190G, Q267R and F392I mutations are first reported in this study. The most common mutations, R241C, R408Q and Ex6-96A>G, account for 23.2%, 12.0% and 9.2%, of the mutant alleles, respectively. Haplotype analysis shows that R241C and Ex6-96A>G are exclusively associated with haplotype 4.3 to suggest founder effects. On the other hand, R408Q is found on two distinct haplotypes suggesting recurrent mutations. The spectrum of PAH mutations in Taiwan shows various links to those of other Asian regions, yet remarkable differences exist. Notably, R408Q, E286K and -4173_-407del, accounting for 21% of all mutant alleles in Taiwan, are very rare or are undetected among PKU cohorts of other Asian regions to suggest local founder effects. Moreover, the low homozygosity value of 0.092 hints at a high degree of ethnic heterogeneity within the Taiwanese population. Our study of PAH mutation spectrum and the associated haplotypes is useful for subsequent study on the origin and migration pattern via Taiwan, an island at the historical crossroad of migration of ancient populations.

A novel mutation (c.951C>T) in an exonic splicing enhancer results in exon 10 skipping in the human mitochondrial acetoacetyl-CoA thiolase gene

Mitochondrial acetoacetyl-CoA thiolase (T2) deficiency is an inherited disorder affecting isoleucine catabolism and ketone body metabolism. A Japanese female developed a severe ketoacidotic attack at the age of 7 months. Urinary organic acid analysis showed elevated excretion of 2-methyl-3-hydroxybutyrate but not tiglylglycine. She was diagnosed as having T2 deficiency by enzyme assay using fibroblasts. Mutation analysis revealed a compound heterozygote of c.556G>T(D186Y) and c.951C>T(D317D). Since c.951C>T does not cause amino acid change, we performed cDNA analysis and found that exon 10 skipping had occurred in the c.951C>T allele. A computer search using an ESE finder showed that an exonic splicing enhancer sequence, SF2/ASF, was located in CTGA(951)CGC. We hypothesized that the exonic splicing enhancer is necessary for accurate splicing since the first nucleotide of exon 10 is C, which weakens the splice acceptor site of intron 9. We made a mini gene construct including exon 9-truncated intron 9-exon 10-truncated intron 10-exon 11 for a splicing experiment. We also made three mutant constructs which alter the SF2/ASF site (947C>T, 951C>T, 952G>A). An min-gene splicing experiment clearly showed that exon 10 skipping was induced in all three mutant constructs. Moreover, additional substitution of G for C at the first nucleotide of exon 10 resulted in normal splicing in these three mutants. These results confirmed that c.951C>T diminished the effect of the exonic splicing enhancer and caused exon 10 skipping.