A synonymous single nucleotide polymorphism in DeltaF508 CFTR alters the secondary structure of the mRNA and the expression of the mutant protein

Mechanistic Approaches to Improve Correction of the Most Common Disease-Causing Mutation in Cystic Fibrosis

The most common mutation in the cystic fibrosis transmembrane conductance regulator (CFTR) gene leads to deletion of the phenylalanine at position 508 (ΔF508) in the CFTR protein and causes multiple folding and functional defects. Contrary to large-scale efforts by industry and academia, no significant therapeutic benefit has been achieved with a single “corrector”. Therefore, investigations concentrate on drug combinations. Orkambi (Vertex Pharmaceuticals), the first FDA-approved drug for treatment of cystic fibrosis (CF) caused by this mutation, is a combination of a corrector (VX-809) that facilitates ΔF508 CFTR biogenesis and a potentiator (VX-770), which improves its function. Yet, clinical trials utilizing this combination showed only modest therapeutic benefit. The low efficacy Orkambi has been attributed to VX-770-mediated destabilization of VX-809-rescued ΔF508 CFTR. Here we report that the negative effects of VX-770 can be reversed by increasing the half-life of the endoplasmic reticulum (ER) form (band B) of ΔF508 CFTR with another corrector (Corr-4a.) Although Corr-4a alone has only minimal effects on ΔF508 CFTR rescue, it increases the half-life of ΔF508 CFTR band B when it is present during half-life measurements. Our data shows that stabilization of band B ΔF508 CFTR with Corr-4a and simultaneous rescue with VX-809, leads to a >2-fold increase in cAMP-activated, CFTR_inh-172-inhibited currents compared to VX-809 alone, or VX-809+VX-770. The negative effects of VX-770 and the Corr-4a protection are specific to the native I507-ATT ΔF508 CFTR without affecting the inherently more stable, synonymous variant I507-ATC ΔF508 CFTR. Our studies emphasize that stabilization of ΔF508 CFTR band B in the ER might improve its functional rescue by Orkambi.

The dark matter of the cancer genome: aberrations in regulatory elements, untranslated regions, splice sites, non-coding RNA and synonymous mutations

Cancer is a disease of the genome caused by oncogene activation and tumor suppressor gene inhibition. Deep sequencing studies including large consortia such as TCGA and ICGC identified numerous tumor‐specific mutations not only in protein‐coding sequences but also in non‐coding sequences. Although 98% of the genome is not translated into proteins, most studies have neglected the information hidden in this “dark matter” of the genome. Malignancy‐driving mutations can occur in all genetic elements outside the coding region, namely in enhancer, silencer, insulator, and promoter as well as in 5′‐UTR and 3′‐UTR. Intron or splice site mutations can alter the splicing pattern. Moreover, cancer genomes contain mutations within non‐coding RNA, such as microRNA, lncRNA, and lincRNA. A synonymous mutation changes the coding region in the DNA and RNA but not the protein sequence. Importantly, oncogenes such as TERT or miR‐21 as well as tumor suppressor genes such as TP53/p53,APC,BRCA1, or RB1 can be affected by these alterations. In summary, coding‐independent mutations can affect gene regulation from transcription, splicing, mRNA stability to translation, and hence, this largely neglected area needs functional studies to elucidate the mechanisms underlying tumorigenesis. This review will focus on the important role and novel mechanisms of these non‐coding or allegedly silent mutations in tumorigenesis.

The Rare Codon AGA Is Involved in Regulation of Pyoluteorin Biosynthesis in Pseudomonas protegens Pf-5

The soil bacterium Pseudomonas protegens Pf-5 can colonize root and seed surfaces of many plants, protecting them from infection by plant pathogenic fungi and oomycetes. The capacity to suppress disease is attributed to Pf-5's production of a large spectrum of antibiotics, which is controlled by complex regulatory circuits operating at the transcriptional and post-transcriptional levels. In this study, we analyzed the genomic sequence of Pf-5 for codon usage patterns and observed that the six rarest codons in the genome are present in all seven known antibiotic biosynthesis gene clusters. In particular, there is an abundance of rare codons in pltR, which encodes a member of the LysR transcriptional regulator family that controls the expression of pyoluteorin biosynthetic genes. To test the hypothesis that rare codons in pltR influence pyoluteorin production, we generated a derivative of Pf-5 in which 23 types of rare codons in pltR were substituted with synonymous preferred codons. The resultant mutant produced pyoluteorin at levels 15 times higher than that of the wild-type Pf-5. Accordingly, the promoter activity of the pyoluteorin biosynthetic gene pltL was 20 times higher in the codon-modified stain than in the wild-type. pltR has six AGA codons, which is the rarest codon in the Pf-5 genome. Substitution of all six AGA codons with preferred Arg codons resulted in a variant of pltR that conferred increased pyoluteorin production and pltL promoter activity. Furthermore, overexpression of tRNAUCUArg, the cognate tRNA for the AGA codon, significantly increased pyoluteorin production by Pf-5. A bias in codon usage has been linked to the regulation of many phenotypes in eukaryotes and prokaryotes but, to our knowledge, this is the first example of the role of a rare codon in the regulation of antibiotic production by a Gram-negative bacterium.

A Variant of GJD2, Encoding for Connexin 36, Alters the Function of Insulin Producing β-Cells

Signalling through gap junctions contributes to control insulin secretion and, thus, blood glucose levels. Gap junctions of the insulin-producing β-cells are made of connexin 36 (Cx36), which is encoded by the GJD2 gene. Cx36-null mice feature alterations mimicking those observed in type 2 diabetes (T2D). GJD2 is also expressed in neurons, which share a number of common features with pancreatic β-cells. Given that a synonymous exonic single nucleotide polymorphism of human Cx36 (SNP rs3743123) associates with altered function of central neurons in a subset of epileptic patients, we investigated whether this SNP also caused alterations of β-cell function. Transfection of rs3743123 cDNA in connexin-lacking HeLa cells resulted in altered formation of gap junction plaques and cell coupling, as compared to those induced by wild type (WT) GJD2 cDNA. Transgenic mice expressing the very same cDNAs under an insulin promoter revealed that SNP rs3743123 expression consistently lead to a post-natal reduction of islet Cx36 levels and β-cell survival, resulting in hyperglycemia in selected lines. These changes were not observed in sex- and age-matched controls expressing WT hCx36. The variant GJD2 only marginally associated to heterogeneous populations of diabetic patients. The data document that a silent polymorphism of GJD2 is associated with altered β-cell function, presumably contributing to T2D pathogenesis.

Whisper mutations: cryptic messages within the genetic code

Recent years have seen a great expansion in our understandings of how silent mutations can drive a disease and that mRNAs are not only mere messengers between the genome and the encoded proteins but also encompass regulatory activities. This review focuses on how silent mutations within open reading frames can affect the functional properties of the encoded protein. We describe how mRNAs exert control of cell biological processes governed by the encoded proteins via translation kinetics, protein folding, mRNA stability, spatio-temporal protein expression and by direct interactions with cellular factors. These examples illustrate how additional levels of information lie within the coding sequences and that the degenerative genetic code is not redundant and have co-evolved with the encoded proteins. Hence, so called synonymous mutations are not always silent but 'whisper'.

Limb-Girdle Muscular Dystrophy Type 2A Resulting From c.C479G and c.G1818A Mutations in the Calpain-3 Gene

Limb-Girdle Muscular Dystrophy type 2A (LGMD2A) is an autosomal recessive disorder characterized by progressive weakness of proximal muscles. Here, we describe a patient with clinical features consistent with LGMD2A who harbors 2 rare changes in the CAPN3 gene sequence of unknown clinical significance. Mechanisms by which these 2 mutations could affect the protein are discussed. The c.C479G mutation seems to affect the proteolytic domain of calpain-3. Whereas the novel mutation c.G1818A seems to affect mRNA translation of the protein region involved in titin binding. We strongly believe that these genomic variants in CAPN3 are indeed deleterious and thus are currently misclassified. Since LGMD2 is considered a disorder of autosomal recessive inheritance, further population studies involving the molecular characterization of symptomatic patients must be performed as well as in vitro studies to ascertain the functional effects of these specific variants.

Breakthrough therapies: Cystic fibrosis (CF) potentiators and correctors

Cystic Fibrosis is caused by mutations in the Cystic Fibrosis Transmembrane conductance Regulator (CFTR) gene resulting in abnormal protein function. Recent advances of targeted molecular therapies and high throughput screening have resulted in multiple drug therapies that target many important mutations in the CFTR protein. In this review, we provide the latest results and current progress of CFTR modulators for the treatment of cystic fibrosis, focusing on potentiators of CFTR channel gating and Phe508del processing correctors for the Phe508del CFTR mutation. Special emphasis is placed on the molecular basis underlying these new therapies and emerging results from the latest clinical trials. The future directions for augmenting the rescue of Phe508del with CFTR modulators are also emphasized.

A synonymous codon change alters the drug sensitivity of ΔF508 cystic fibrosis transmembrane conductance regulator

Synonymous mutations, such as I507-ATC→ATT, in deletion of Phe508 in cystic fibrosis transmembrane conductance regulator (ΔF508 CFTR), the most frequent disease-associated mutant of CFTR, may affect protein biogenesis, structure, and function and contribute to an altered disease phenotype. Small-molecule drugs are being developed to correct ΔF508 CFTR. To understand correction mechanisms and the consequences of synonymous mutations, we analyzed the effect of mechanistically distinct correctors, corrector 4a (C4) and lumacaftor (VX-809), on I507-ATT and I507-ATC ΔF508 CFTR biogenesis and function. C4 stabilized I507-ATT ΔF508 CFTR band B, but without considerable biochemical and functional correction. VX-809 biochemically corrected ∼10% of both of the variants, leading to stable, forskolin+3-isobutyl-1-methylxanthine (IBMX)-activated whole-cell currents in the presence of the corrector. Omitting VX-809 during whole-cell recordings led to a spontaneous decline of the currents, suggesting posttranslational stabilization by VX-809. Treatment of cells with the C4+VX-809 combination resulted in enhanced rescue and 2-fold higher forskolin+IBMX-activated currents of both I507-ATT and I507-ATC ΔF508 CFTR, compared with VX-809 treatment alone. The lack of an effect of C4 on I507-ATC ΔF508 CFTR, but its additive effect in combination with VX-809, implies that C4 acted on VX-809-modified I507-ATC ΔF508 CFTR. Our results suggest that binding of C4 and VX-809 to ΔF508 CFTR is conformation specific and provide evidence that synonymous mutations can alter the drug sensitivity of proteins.

Influenza virus M2 targets cystic fibrosis transmembrane conductance regulator for lysosomal degradation during viral infection

We sought to determine the mechanisms by which influenza infection of human epithelial cells decreases cystic fibrosis transmembrane conductance regulator (CFTR) expression and function. We infected human bronchial epithelial (NHBE) cells and murine nasal epithelial (MNE) cells with various strains of influenza A virus. Influenza infection significantly reduced CFTR short circuit currents (Isc) and protein levels at 8 hours postinfection. We then infected CFTR expressing human embryonic kidney (HEK)-293 cells (HEK-293 CFTRwt) with influenza virus encoding a green fluorescent protein (GFP) tag and performed whole-cell and cell-attached patch clamp recordings. Forskolin-stimulated, GlyH-101-sensitive CFTR conductances, and CFTR open probabilities were reduced by 80% in GFP-positive cells; Western blots also showed significant reduction in total and plasma membrane CFTR levels. Knockdown of the influenza matrix protein 2 (M2) with siRNA, or inhibition of its activity by amantadine, prevented the decrease in CFTR expression and function. Lysosome inhibition (bafilomycin-A1), but not proteasome inhibition (lactacystin), prevented the reduction in CFTR levels. Western blots of immunoprecipitated CFTR from influenza-infected cells, treated with BafA1, and probed with antibodies against lysine 63-linked (K-63) or lysine 48-linked (K-48) polyubiquitin chains supported lysosomal targeting. These results highlight CFTR damage, leading to early degradation as an important contributing factor to influenza infection-associated ion transport defects.

Optimization of translation profiles enhances protein expression and solubility

mRNA is translated with a non-uniform speed that actively coordinates co-translational folding of protein domains. Using structure-based homology we identified the structural domains in epoxide hydrolases (EHs) and introduced slow-translating codons to delineate the translation of single domains. These changes in translation speed dramatically improved the solubility of two EHs of metagenomic origin in Escherichia coli. Conversely, the importance of transient attenuation for the folding, and consequently solubility, of EH was evidenced with a member of the EH family from Agrobacterium radiobacter, which partitions in the soluble fraction when expressed in E. coli. Synonymous substitutions of codons shaping the slow-transiting regions to fast-translating codons render this protein insoluble. Furthermore, we show that low protein yield can be enhanced by decreasing the free folding energy of the initial 5’-coding region, which can disrupt mRNA secondary structure and enhance ribosomal loading. This study provides direct experimental evidence that mRNA is not a mere messenger for translation of codons into amino acids but bears an additional layer of information for folding, solubility and expression level of the encoded protein. Furthermore, it provides a general frame on how to modulate and fine-tune gene expression of a target protein.

The nicotinic acetylcholine receptor alpha 4 subunit contains a functionally relevant SNP Haplotype

BACKGROUND

Non-coding single nucleotide polymorphisms within the nicotinic acetylcholine receptor alpha 4 subunit gene (CHRNA4) are robustly associated with various neurological and behavioral phenotypes including schizophrenia, cognition and smoking. The most commonly associated polymorphisms are located in exon 5 and segregate as part of a haplotype. So far it is unknown if this haplotype is indeed functional, or if the observed associations are an indirect effect caused by linkage disequilibrium with not yet identified adjacent functional variants. We therefore analyzed the functional relevance of the exon 5 haplotype alleles.

RESULTS

Using voltage clamp experiments we were able to show that the CHRNA4 haplotype alleles differ with respect to their functional effects on receptor sensitivity including reversal of receptor sensitivity between low and high acetylcholine concentrations. The results indicate that underlying mechanisms might include differences in codon usage bias and changes in mRNA stability.

CONCLUSIONS

Our data demonstrate that the complementary alleles of the CHRNA4 exon 5 haplotype are functionally relevant, and might therefore be causative for the above mentioned associations.

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.

Novel Splice Site Mutation in MAMLD1 in a Patient with Hypospadias

MAMLD1 is a causative gene for disorders of sex development. Several MAMLD1 mutations have been shown to cause hypospadias by generating dysfunctional proteins and/or unstable mRNAs. Here, we identified an intronic mutation of MAMLD1 (g.IVS4-2A>G) in 1 of 180 hypospadias patients. RT-PCR of the patient's skin sample showed normal expression of full-length MAMLD1 and markedly reduced expression of a known splice variant lacking exon 4. A hitherto unreported splice variant that lacks exon 5 was similarly identified in samples of the patient and control individuals. The full-length transcript of the patient contained mutant mRNA lacking the first 10 nucleotides of exon 5 (c.1822_1831delACTCATGTAG, p.K609fsX1070). In vitro assays using cells expressing the full-length wild-type and mutant proteins revealed reduced expression of the mutant. The expression of the wild-type and mutant MAMLD1 showed parallel changes upon treatment with a proteasome inhibitor and a translation inhibitor. The mutant-expressing cells exerted low transactivation activity for the Hes3 promoter, which reflected limited expression of the mutant protein. These results imply that the pathogenic events resulting from MAMLD1 mutations include splice errors. Furthermore, this study raises the possibility of translation failure of MAMLD1 mutants, which deserves further investigation.

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.

Why Selection Might Be Stronger When Populations Are Small: Intron Size and Density Predict within and between-Species Usage of Exonic Splice Associated cis-Motifs

The nearly neutral theory predicts that small effective population size provides the conditions for weakened selection. This is postulated to explain why our genome is more “bloated” than that of, for example, yeast, ours having large introns and large intergene spacer. If a bloated genome is also an error prone genome might it, however, be the case that selection for error-mitigating properties is stronger in our genome? We examine this notion using splicing as an exemplar, not least because large introns can predispose to noisy splicing. We thus ask whether, owing to genomic decay, selection for splice error-control mechanisms is stronger, not weaker, in species with large introns and small populations. In humans much information defining splice sites is in cis-exonic motifs, most notably exonic splice enhancers (ESEs). These act as splice-error control elements. Here then we ask whether within and between-species intron size is a predictor of the commonality of exonic cis-splicing motifs. We show that, as predicted, the proportion of synonymous sites that are ESE-associated and under selection in humans is weakly positively correlated with the size of the flanking intron. In a phylogenetically controlled framework, we observe, also as expected, that mean intron size is both predicted by N_e.μ and is a good predictor of cis-motif usage across species, this usage coevolving with splice site definition. Unexpectedly, however, across taxa intron density is a better predictor of cis-motif usage than intron size. We propose that selection for splice-related motifs is driven by a need to avoid decoy splice sites that will be more common in genes with many and large introns. That intron number and density predict ESE usage within human genes is consistent with this, as is the finding of intragenic heterogeneity in ESE density. As intronic content and splice site usage across species is also well predicted by N_e.μ, the result also suggests an unusual circumstance in which selection (for cis-modifiers of splicing) might be stronger when population sizes are smaller, as here splicing is noisier, resulting in a greater need to control error-prone splicing.

Roles for Synonymous Codon Usage in Protein Biogenesis

Owing to the degeneracy of the genetic code, a protein sequence can be encoded by many different synonymous mRNA coding sequences. Synonymous codon usage was once thought to be functionally neutral, but evidence now indicates it is shaped by evolutionary selection and affects other aspects of protein biogenesis beyond specifying the amino acid sequence of the protein. Synonymous rare codons, once thought to have only negative impacts on the speed and accuracy of translation, are now known to play an important role in diverse functions, including regulation of cotranslational folding, covalent modifications, secretion, and expression level. Mutations altering synonymous codon usage are linked to human diseases. However, much remains unknown about the molecular mechanisms connecting synonymous codon usage to efficient protein biogenesis and proper cell physiology. Here we review recent literature on the functional effects of codon usage, including bioinformatics approaches aimed at identifying general roles for synonymous codon usage.

Synonymous codon usage affects the expression of wild type and F508del CFTR

The cystic fibrosis transmembrane conductance regulator (CFTR) is an anion channel composed of 1480 amino acids. The major mutation responsible for cystic fibrosis results in loss of amino acid residue, F508 (F508del). Loss of F508 in CFTR alters the folding pathway resulting in endoplasmic-reticulum-associated degradation. This study investigates the role of synonymous codon in the expression of CFTR and CFTR F508del in human HEK293 cells. DNA encoding the open reading frame (ORF) for CFTR containing synonymous codon replacements was expressed using a heterologous vector integrated into the genome. The results indicate that the codon usage greatly affects the expression of CFTR. While the promoter strength driving expression of the ORFs was largely unchanged and the mRNA half-lives were unchanged, the steady-state levels of the mRNA varied by as much as 30-fold. Experiments support that this apparent inconsistency is attributed to nonsense mediated decay independent of exon junction complex. The ratio of CFTR/mRNA indicates that mRNA containing native codons was more efficient in expressing mature CFTR as compared to mRNA containing synonymous high-expression codons. However, when F508del CFTR was expressed after codon optimization, a greater percentage of the protein escaped endoplasmic-reticulum-associated degradation resulting in considerable levels of mature F508del CFTR on the plasma membrane, which showed channel activity. These results indicate that codon usage has an effect on mRNA levels and protein expression, for CFTR, and likely on chaperone-assisted folding pathway, for F508del CFTR.

Association of COMT and TPH-2 genes with DSM-5 based PTSD symptoms

BACKGROUND

Dopaminergic and serotonergic systems have been implicated in PTSD. The present study evaluated the association of four catechol-O-methyltransferase (COMT) gene loci, and the joint effect of COMT and tryptophan hydroxylase 2 (TPH2) genes on PTSD symptoms.

METHODS

Subjects included 200 Caucasian Armenian adults exposed to the 1988 Spitak earthquake from 12 multigenerational (3-5 generations) families. Instruments used included the UCLA PTSD Reaction Index based on DSM-5 criteria, and the Beck Depression Inventory.

RESULTS

The adjusted heritabilitiy of vulnerability to DSM-5 based PTSD symptoms was 0.60 (p<10(-4)). There was a significant association of the COMT allele rs4633C with total PTSD (p<0.03), and D category (p<0.04) (negative alterations in cognitions and mood) severity scores, but not with C category (avoidance) scores. There was no genetic correlation between C and D category severity scores. COMT allele rs4633C and the TPH-2 allele rs11178997T together accounted for 7% of the variance in PTSD severity scores (p<0.001). None of the COMT alleles were associated with depression.

LIMITATIONS

The ratings of earthquake exposure and prior trauma may have been subject to recall bias. The findings may not be generalizable to other ethnic/racial populations.

CONCLUSION

COMT allele rs4633C may be causally related and/or is in linkage disequilibrium with gene(s) that are causally related to PTSD symptoms. Carriers of these COMT and the TPH-2 alleles may be at increased risk for PTSD. The findings provide biological support for dividing DSM-IV category C symptoms into DSM-5 categories C and D.

Timing is everything: unifying codon translation rates and nascent proteome behavior

Experiments have demonstrated that changing the rate at which the ribosome translates a codon position in an mRNA molecule's open reading frame can alter the behavior of the newly synthesized protein. That is, codon translation rates can govern nascent proteome behavior. We emphasize that this phenomenon is a manifestation of the nonequilibrium nature of cotranslational processes, and as such, there exist theoretical tools that offer a potential means to quantitatively predict the influence of codon translation rates on the broad spectrum of nascent protein behaviors including cotranslational folding, aggregation, and translocation. We provide a review of the experimental evidence for the impact that codon translation rates can have, followed by a discussion of theoretical methods that can describe this phenomenon. The development and application of these tools are likely to provide fundamental insights into protein maturation and homeostasis, codon usage bias in organisms, the origins of translation related diseases, and new rational design methods for biotechnology and biopharmaceutical applications.

Exposing synonymous mutations

Synonymous codon changes, which do not alter protein sequence, were previously thought to have no functional consequence. Although this concept has been overturned in recent years, there is no unique mechanism by which these changes exert biological effects. A large repertoire of both experimental and bioinformatic methods has been developed to understand the effects of synonymous variants. Results from this body of work have provided global insights into how biological systems exploit the degeneracy of the genetic code to control gene expression, protein folding efficiency, and the coordinated expression of functionally related gene families. Although it is now clear that synonymous variants are important in a variety of contexts, from human disease to the safety and efficacy of therapeutic proteins, there is no clear consensus on the approaches to identify and validate these changes. Here, we review the diverse methods to understand the effects of synonymous mutations.

RNAdigest: a web-based tool for the analysis and prediction of structure-specific RNAse digestion results

Despite recent developments in analyzing RNA secondary structures, relatively few RNA structures have been determined. To date, many investigators have relied on the traditional method of using structure-specific RNAse enzymes to probe RNA secondary structures. However, if these data were combined with novel computational approaches, investigators would have an informative and valuable tool for RNA structural analysis. To this end, we created the web server “RNAdigest.” RNAdigest uses mfold RNA structural models in order to predict the results of RNAse digestion experiments. Furthermore, RNAdigest also utilizes both RNA sequence and the experimental digestion patterns to formulate the constraints for predicting secondary structures of the RNA. Thus, RNAdigest allows for the structural interpretation of RNAse digestion experiments. Overall, RNAdigest simplifies RNAse digestion result analyses while allowing for the identification of unique fragments. These unique fragments can then be used for testing predicted mfold structures and for designing structural-specific DNA/RNA probes.

The impact of RNA structure on coding sequence evolution in both bacteria and eukaryotes

Background

Many studies have found functional RNA secondary structures are selectively conserved among species. But, the effect of RNA structure selection on coding sequence evolution remains unknown. To address this problem, we systematically investigated the relationship between nucleotide conservation level and its structural sensitivity in four model organisms, Escherichia coli, yeast, fly, and mouse.

Results

We define structurally sensitive sites as those with putative local structure-disruptive mutations. Using both the Mantel-Haenszel procedure and association test, we found structurally sensitive nucleotide sites evolved more slowly than non-sensitive sites in all four organisms. Furthermore, we observed that this association is more obvious in highly expressed genes and region near the start codon.

Conclusion

We conclude that structurally sensitive sites in mRNA sequences normally have less nucleotide divergence in all species we analyzed. This study extends our understanding of the impact of RNA structure on coding sequence evolution, and is helpful to the development of a codon model with RNA structure information.

Association of CREB1 gene polymorphism with drug seeking behaviour in eastern Indian addicts

cAMP response element binding protein (CREB) is a major transcription factor which plays an important role in a wide array of cellular functions. CREB also has a significant function in developing substance abuse. A study was undertaken to identify the single nucleotide polymorphisms (SNP) at selective areas of CREB1 gene in heroin as well as in alcohol addicts in comparison with control population. One hundred and forty control subjects, 112 heroin and 102 alcoholics, all male and residing in Kolkata, a city in eastern India participated in the study. SNPs from several exonic regions of CREB1 gene were assessed to investigate possible associations with addiction. One SNP in exon 3, rs35349697, demonstrated a significant correlation with opioid addiction as well as with alcohol addiction. A novel SNP, also located in exon 3, was identified which showed epistatic interaction with rs35349697 to decrease susceptibility to narcotic addiction in the population. The study is the first report on the identification of a role of CREB1 gene polymorphism with addiction.

Post-transcriptional regulation of cystic fibrosis transmembrane conductance regulator expression and function by microRNAs

MicroRNAs (miRNAs) are increasingly recognized as important posttranscriptional regulators of gene expression, and changes in their actions can contribute to disease states. Little is understood regarding miRNA functions in the airway epithelium under normal or diseased conditions. We profiled miRNA expression in well-differentiated primary cultures of human cystic fibrosis (CF) and non-CF airway epithelia, and discovered that miR-509-3p and miR-494 concentrations were increased in CF epithelia. Human non-CF airway epithelia, transfected with the mimics of miR-509-3p or miR-494, showed decreased cystic fibrosis transmembrane conductance regulator (CFTR) expression, whereas their respective anti-miRs exerted the opposite effect. Interestingly, the two miRNAs acted cooperatively in regulating CFTR expression. Upon infecting non-CF airway epithelial cells with Staphylococcus aureus, or upon stimulating them with the proinflammatory cytokines TNF-α or IL-1β, we observed an increased expression of both miRNAs and a concurrent decrease in CFTR expression and function, suggesting that inflammatory mediators may regulate these miRNAs. Transfecting epithelia with anti-miRs for miR-509-3p and miR-494, or inhibiting NF-κB signaling before stimulating cells with TNFα or IL-1β, suppressed these responses, suggesting that the expression of both miRNAs was responsive to NF-κB signaling. Thus, miR-509-3p and miR-494 are dynamic regulators of CFTR abundance and function in normal, non-CF airway epithelia.

Associations between variants of FADS genes and omega-3 and omega-6 milk fatty acids of Canadian Holstein cows

BACKGROUND

Fatty acid desaturase 1 (FADS1) and 2 (FADS2) genes code respectively for the enzymes delta-5 and delta-6 desaturases which are rate limiting enzymes in the synthesis of polyunsaturated omega-3 and omega-6 fatty acids (FAs). Omega-3 and-6 FAs as well as conjugated linoleic acid (CLA) are present in bovine milk and have demonstrated positive health effects in humans. Studies in humans have shown significant relationships between genetic variants in FADS1 and 2 genes with plasma and tissue concentrations of omega-3 and-6 FAs. The aim of this study was to evaluate the extent of sequence variations within these two genes in Canadian Holstein cows as well as the association between sequence variants and health promoting FAs in milk.

RESULTS

Thirty three SNPs were detected within the studied regions of genes including a synonymous mutation (FADS1-07, rs42187261, 306Tyr > Tyr) in exon 8 of FADS1, a non-synonymous mutation (FADS2-14, rs211580559, 294Ala > Val) within FADS2 exon 7, a splice site SNP (FADS2-05, rs211263660), a 3'UTR SNP (FADS2-23, rs109772589), and another 3'UTR SNP with an effect on a microRNA binding site within FADS2 gene (FADS2-19, rs210169303). Association analyses showed significant relations between three out of seven tested SNPs and several FAs. Significant associations (FDR P < 0.05) were recorded between FADS2-23 (rs109772589) and two omega-6 FAs (dihomogamma linolenic acid [C20:3n6] and arachidonic acid [C20:4n6]), FADS1-07 (rs42187261) and one omega-3 FA (eicosapentaenoic acid, C20:5n3) and tricosanoic acid (C23:0), and one intronic SNP, FADS1-01 (rs136261927) and C20:3n6.

CONCLUSION

Our study has demonstrated positive associations between three SNPs within FADS1 and FADS2 genes (a SNP within the 3'UTR, a synonymous SNP and an intronic SNP), with three milk PUFAs of Canadian Holstein cows thus suggesting possible involvement of synonymous and non-coding region variants in FA synthesis. These SNPs may serve as potential genetic markers in breeding programs to increase milk FAs that are of benefit to human health.

A silent exonic SNP in kdm3a affects nucleic acids structure but does not regulate experimental autoimmune encephalomyelitis

Defining genetic variants that predispose for diseases is an important initiative that can improve biological understanding and focus therapeutic development. Genetic mapping in humans and animal models has defined genomic regions controlling a variety of phenotypes known as quantitative trait loci (QTL). Causative disease determinants, including single nucleotide polymorphisms (SNPs), lie within these regions and can often be identified through effects on gene expression. We previously identified a QTL on rat chromosome 4 regulating macrophage phenotypes and immune-mediated diseases including experimental autoimmune encephalomyelitis (EAE). Gene analysis and a literature search identified lysine-specific demethylase 3A (Kdm3a) as a potential regulator of these phenotypes. Genomic sequencing determined only two synonymous SNPs in Kdm3a. The silent synonymous SNP in exon 15 of Kdm3a caused problems with quantitative PCR detection in the susceptible strain through reduced amplification efficiency due to altered secondary cDNA structure. Shape Probability Shift analysis predicted that the SNP often affects RNA folding; thus, it may impact protein translation. Despite these differences in rats, genetic knockout of Kdm3a in mice resulted in no dramatic effect on immune system development and activation or EAE susceptibility and severity. These results provide support for tools that analyze causative SNPs that impact nucleic acid structures.

Resveratrol enhances airway surface liquid depth in sinonasal epithelium by increasing cystic fibrosis transmembrane conductance regulator open probability

BACKGROUND

Chronic rhinosinusitis engenders enormous morbidity in the general population, and is often refractory to medical intervention. Compounds that augment mucociliary clearance in airway epithelia represent a novel treatment strategy for diseases of mucus stasis. A dominant fluid and electrolyte secretory pathway in the nasal airways is governed by the cystic fibrosis transmembrane conductance regulator (CFTR). The objectives of the present study were to test resveratrol, a strong potentiator of CFTR channel open probability, in preparation for a clinical trial of mucociliary activators in human sinus disease.

METHODS

Primary sinonasal epithelial cells, immortalized bronchoepithelial cells (wild type and F508del CFTR), and HEK293 cells expressing exogenous human CFTR were investigated by Ussing chamber as well as patch clamp technique under non-phosphorylating conditions. Effects on airway surface liquid depth were measured using confocal laser scanning microscopy. Impact on CFTR gene expression was measured by quantitative reverse transcriptase polymerase chain reaction.

RESULTS

Resveratrol is a robust CFTR channel potentiator in numerous mammalian species. The compound also activated temperature corrected F508del CFTR and enhanced CFTR-dependent chloride secretion in human sinus epithelium ex vivo to an extent comparable to the recently approved CFTR potentiator, ivacaftor. Using inside out patches from apical membranes of murine cells, resveratrol stimulated an ~8 picosiemens chloride channel consistent with CFTR. This observation was confirmed in HEK293 cells expressing exogenous CFTR. Treatment of sinonasal epithelium resulted in a significant increase in airway surface liquid depth (in µm: 8.08+/-1.68 vs. 6.11+/-0.47,control,p<0.05). There was no increase CFTR mRNA.

CONCLUSION

Resveratrol is a potent chloride secretagogue from the mucosal surface of sinonasal epithelium, and hydrates airway surface liquid by increasing CFTR channel open probability. The foundation for a clinical trial utilizing resveratrol as a therapeutic intervention to increase mucociliary transport and airway surface liquid hydration in sinus disease is strongly supported by these findings.

The silent codon change I507-ATC->ATT contributes to the severity of the ΔF508 CFTR channel dysfunction

The most common disease-causing mutation in the cystic fibrosis transmembrane conductance regulator (CFTR) gene is the out-of-frame deletion of 3 nucleotides (CTT). This mutation leads to the loss of phenylalanine-508 (ΔF508) and a silent codon change (SCC) for isoleucine-507 (I507-ATC→ATT). ΔF508 CFTR is misfolded and degraded by endoplasmic reticulum-associated degradation (ERAD). We have demonstrated that the I507-ATC→ATT SCC alters ΔF508 CFTR mRNA structure and translation dynamics. By comparing the biochemical and functional properties of the I507-ATT and I507-ATC ΔF508 CFTR, we establish that the I507-ATC→ATT SCC contributes to the cotranslational misfolding, ERAD, and to the functional defects associated with ΔF508 CFTR. We demonstrate that the I507-ATC ΔF508 CFTR is less susceptible to the ER quality-control machinery during translation than the I507-ATT, although 27°C correction is necessary for sufficient cell-surface expression. Whole-cell patch-clamp recordings indicate sustained, thermally stable cAMP-activated Cl(-) transport through I507-ATC and unstable function of the I507-ATT ΔF508 CFTR. Single-channel recordings reveal improved gating properties of the I507-ATC compared to I507-ATT ΔF508 CFTR (NPo=0.45±0.037 vs. NPo=0.09±0.002; P<0.001). Our results signify the role of the I507-ATC→ATT SCC in the ΔF508 CFTR defects and support the importance of synonymous codon choices in determining the function of gene products.

RNAsnp: efficient detection of local RNA secondary structure changes induced by SNPs

Structural characteristics are essential for the functioning of many noncoding RNAs and cis-regulatory elements of mRNAs. SNPs may disrupt these structures, interfere with their molecular function, and hence cause a phenotypic effect. RNA folding algorithms can provide detailed insights into structural effects of SNPs. The global measures employed so far suffer from limited accuracy of folding programs on large RNAs and are computationally too demanding for genome-wide applications. Here, we present a strategy that focuses on the local regions of maximal structural change between mutant and wild-type. These local regions are approximated in a “screening mode” that is intended for genome-wide applications. Furthermore, localized regions are identified as those with maximal discrepancy. The mutation effects are quantified in terms of empirical P values. To this end, the RNAsnp software uses extensive precomputed tables of the distribution of SNP effects as function of length and GC content. RNAsnp thus achieves both a noise reduction and speed-up of several orders of magnitude over shuffling-based approaches. On a data set comprising 501 SNPs associated with human-inherited diseases, we predict 54 to have significant local structural effect in the untranslated region of mRNAs. RNAsnp is available at http://rth.dk/resources/rnasnp.

A gene-specific method for predicting hemophilia-causing point mutations

A fundamental goal of medical genetics is the accurate prediction of genotype-phenotype correlations. As an approach to develop more accurate in silico tools for prediction of disease-causing mutations of structural proteins, we present a gene- and disease-specific prediction tool based on a large systematic analysis of missense mutations from hemophilia A (HA) patients. Our HA-specific prediction tool, HApredictor, showed disease prediction accuracy comparable to other publicly available prediction software. In contrast to those methods, its performance is not limited to non-synonymous mutations. Given the role of synonymous mutations in disease and drug codon optimization, we propose that utilizing a gene- and disease-specific method can be highly useful to make functional predictions possible even for synonymous mutations. Incorporating computational metrics at both nucleotide and amino acid levels along with multiple protein sequence/structure alignment significantly improved the predictive performance of our tool. HApredictor is freely available for download at http://www.ncbi.nlm.nih.gov/CBBresearch/Przytycka/HA_Predict/index.htm.

Crucial role for phylogenetically conserved cytoplasmic loop 3 in ABCC4 protein expression

The ABC transporter ABCC4 is recognized as an ATP-dependent exporter of endogenous substances as well as an increasing variety of anionic chemotherapeutics. A loss-of-function variant of zebrafish Abcc4 was identified with a single amino acid substitution in the cytoplasmic loop T804M. Because this substituted amino acid is highly conserved among ABCC4 orthologs and is located in cytoplasmic loop 3 (CL3), we investigated the impact of this mutation on human and zebrafish Abcc4 expression. We demonstrate that zebrafish Abcc4 T804M or human ABCC4 T796M exhibit substantially reduced expression, coupled with impaired plasma membrane localization. To understand the molecular basis for the localization defect, we developed a homology model of zebrafish Abcc4. The homology model suggested that the bulky methionine substitution disrupted side-chain contacts. Molecular dynamic simulations of a fragment of human or zebrafish CL3 containing a methionine substitution indicated altered helicity coupled with reduced thermal stability. Trifluoroethanol challenge coupled with circular dichroism revealed that the methionine substitution disrupted the ability of this fragment of CL3 to readily form an α-helix. Furthermore, expression and plasma membrane localization of these mutant ABCC4/Abcc4 proteins are mostly rescued by growing cells at subphysiological temperatures. Because the cystic fibrosis transmembrane conductance regulator (ABCC7) is closely related to ABCC4, we extended this by engineering certain pathogenic CFTR-CL3 mutations, and we showed they destabilized human and zebrafish ABCC4. Altogether, our studies provide the first evidence for a conserved domain in CL3 of ABCC4 that is crucial in ensuring its proper plasma membrane localization.

Age-related macular degeneration-associated silent polymorphisms in HtrA1 impair its ability to antagonize insulin-like growth factor 1

Synonymous single nucleotide polymorphisms (SNPs) within a transcript's coding region produce no change in the amino acid sequence of the protein product and are therefore intuitively assumed to have a neutral effect on protein function. We report that two common variants of high-temperature requirement A1 (HTRA1) that increase the inherited risk of neovascular age-related macular degeneration (NvAMD) harbor synonymous SNPs within exon 1 of HTRA1 that convert common codons for Ala34 and Gly36 to less frequently used codons. The frequent-to-rare codon conversion reduced the mRNA translation rate and appeared to compromise HtrA1's conformation and function. The protein product generated from the SNP-containing cDNA displayed enhanced susceptibility to proteolysis and a reduced affinity for an anti-HtrA1 antibody. The NvAMD-associated synonymous polymorphisms lie within HtrA1's putative insulin-like growth factor 1 (IGF-1) binding domain. They reduced HtrA1's abilities to associate with IGF-1 and to ameliorate IGF-1-stimulated signaling events and cellular responses. These observations highlight the relevance of synonymous codon usage to protein function and implicate homeostatic protein quality control mechanisms that may go awry in NvAMD.

Sequence and structure-specific elements of HERG mRNA determine channel synthesis and trafficking efficiency

Human ether-á-gogo-related gene (HERG) encodes a potassium channel that is highly susceptible to deleterious mutations resulting in susceptibility to fatal cardiac arrhythmias. Most mutations adversely affect HERG channel assembly and trafficking. Why the channel is so vulnerable to missense mutations is not well understood. Since nothing is known of how mRNA structural elements factor in channel processing, we synthesized a codon-modified HERG cDNA (HERG-CM) where the codons were synonymously changed to reduce GC content, secondary structure, and rare codon usage. HERG-CM produced typical IKr-like currents; however, channel synthesis and processing were markedly different. Translation efficiency was reduced for HERG-CM, as determined by heterologous expression, in vitro translation, and polysomal profiling. Trafficking efficiency to the cell surface was greatly enhanced, as assayed by immunofluorescence, subcellular fractionation, and surface labeling. Chimeras of HERG-NT/CM indicated that trafficking efficiency was largely dependent on 5' sequences, while translation efficiency involved multiple areas. These results suggest that HERG translation and trafficking rates are independently governed by noncoding information in various regions of the mRNA molecule. Noncoding information embedded within the mRNA may play a role in the pathogenesis of hereditary arrhythmia syndromes and could provide an avenue for targeted therapeutics.

Influenza matrix protein 2 alters CFTR expression and function through its ion channel activity

The human cystic fibrosis transmembrane conductance regulator (CFTR) is a cyclic AMP-activated chloride (Cl(-)) channel in the lung epithelium that helps regulate the thickness and composition of the lung epithelial lining fluid. We investigated whether influenza M2 protein, a pH-activated proton (H(+)) channel that traffics to the plasma membrane of infected cells, altered CFTR expression and function. M2 decreased CFTR activity in 1) Xenopus oocytes injected with human CFTR, 2) epithelial cells (HEK-293) stably transfected with CFTR, and 3) human bronchial epithelial cells (16HBE14o-) expressing native CFTR. This inhibition was partially reversed by an inhibitor of the ubiquitin-activating enzyme E1. Next we investigated whether the M2 inhibition of CFTR activity was due to an increase of secretory organelle pH by M2. Incubation of Xenopus oocytes expressing CFTR with ammonium chloride or concanamycin A, two agents that alkalinize the secretory pathway, inhibited CFTR activity in a dose-dependent manner. Treatment of M2- and CFTR-expressing oocytes with the M2 ion channel inhibitor amantadine prevented the loss in CFTR expression and activity; in addition, M2 mutants, lacking the ability to transport H(+), did not alter CFTR activity in Xenopus oocytes and HEK cells. Expression of an M2 mutant retained in the endoplasmic reticulum also failed to alter CFTR activity. In summary, our data show that M2 decreases CFTR activity by increasing secretory organelle pH, which targets CFTR for destruction by the ubiquitin system. Alteration of CFTR activity has important consequences for fluid regulation and may potentially modify the immune response to viral infection.

The cystic fibrosis gene: a molecular genetic perspective

The positional cloning of the gene responsible for cystic fibrosis (CF) was the important first step in understanding the basic defect and pathophysiology of the disease. This study aims to provide a historical account of key developments as well as factors that contributed to the cystic fibrosis transmembrane conductance regulator (CFTR) gene identification work. A redefined gene structure based on the full sequence of the gene derived from the Human Genome Project is presented, along with brief reviews of the transcription regulatory sequences for the CFTR gene, the role of mRNA splicing in gene regulation and CF disease, and, various related sequences in the human genome and other species. Because CF mutations and genotype-phenotype correlations are covered by our colleagues (Ferec C, Cutting GR. 2012. Assessing the disease-liability of mutations in CFTR. Cold Spring Harb Perspect Med doi: 10.1101/cshperspect.a009480), we only attempt to provide an introduction of the CF mutation database here for reference purposes.

Genetics is a major determinant of expression of the human hepatic uptake transporter OATP1B1, but not of OATP1B3 and OATP2B1

BACKGROUND

Organic anion transporting polypeptide (OATP) 1B1, OATP1B3, and OATP2B1 (encoded by SLCO1B1, SLCO1B3, SLCO2B1) mediate the hepatic uptake of endogenous compounds like bile acids and of drugs, for example, the lipid-lowering atorvastatin, thereby influencing hepatobiliary elimination. Here we systematically elucidated the contribution of SLCO variants on expression of the three hepatic OATPs under consideration of additional important covariates.

METHODS

Expression was quantified by RT-PCR and immunoblotting in 143 Caucasian liver samples. A total of 109 rare and common variants in the SLCO1B3-SLCO1B1 genomic region and the SLCO2B1 gene were genotyped by MALDI-TOF mass spectrometry and genome-wide SNP microarray technology. SLCO1B1 haplotypes affecting hepatic OATP1B1 expression were associated with pharmacokinetic data of the OATP1B1 substrate atorvastatin (n = 82).

RESULTS

Expression of OATP1B1, OATP1B3, and OATP2B1 at the mRNA and protein levels showed marked interindividual variability. All three OATPs were expressed in a coordinated fashion. By a multivariate regression analysis adjusted for non-genetic and transcription covariates, increased OATP1B1 expression was associated with the coding SLCO1B1 variant c.388A > G (rs2306283) even after correction for multiple testing (P = 0.00034). This held true for haplotypes harboring c.388A > G but not the functional variant c.521T > C (rs4149056) associated with statin-related myopathy. c.388A > G also significantly affected atorvastatin pharmacokinetics. SLCO variants and non-genetic and regulatory covariates together accounted for 59% of variability of OATP1B1 expression.

CONCLUSIONS

Our results show that expression of OATP1B1, but not of OATP1B3 and OATP2B1, is significantly affected by genetic variants. The SLCO1B1 variant c.388A > G is the major determinant with additional consequences on atorvastatin plasma levels.

A common polymorphism in the LDL receptor gene has multiple effects on LDL receptor function

A common synonymous single nucleotide polymorphism in exon 12 of the low-density lipoprotein receptor (LDLR) gene, rs688, has been associated with increased plasma total and LDL cholesterol in several populations. Using immortalized lymphoblastoid cell lines from a healthy study population, we confirmed an earlier report that the minor allele of rs688 is associated with increased exon 12 alternative splicing (P < 0.05) and showed that this triggered nonsense-mediated decay (NMD) of the alternatively spliced LDLR mRNA. However, since synonymous single nucleotide polymorphisms may influence structure and function of the encoded proteins by co-translational effects, we sought to test whether rs688 was also functional in the full-length mRNA. In HepG2 cells expressing LDLR cDNA constructs engineered to contain the major or minor allele of rs688, the latter was associated with a smaller amount of LDLR protein at the cell surface (-21.8 ± 0.6%, P = 0.012), a higher amount in the lysosome fraction (+25.7 ± 0.3%, P = 0.037) and reduced uptake of fluorescently labeled LDL (-24.3 ± 0.7%, P < 0.01). Moreover, in the presence of exogenous proprotein convertase subtilisin/kexin type 9 (PCSK9), a protein that reduces cellular LDL uptake by promoting lysosomal degradation of LDLR, the minor allele resulted in reduced capacity of a PCSK9 monoclonal antibody to increase LDL uptake. These findings are consistent with the hypothesis that rs688, which is located in the β-propeller region of LDLR, has effects on LDLR activity beyond its role in alternative splicing due to impairment of LDLR endosomal recycling and/or PCSK9 binding, processes in which the β-propeller is critically involved.

Sounds of silence: synonymous nucleotides as a key to biological regulation and complexity

Messenger RNA is a key component of an intricate regulatory network of its own. It accommodates numerous nucleotide signals that overlap protein coding sequences and are responsible for multiple levels of regulation and generation of biological complexity. A wealth of structural and regulatory information, which mRNA carries in addition to the encoded amino acid sequence, raises the question of how these signals and overlapping codes are delineated along non-synonymous and synonymous positions in protein coding regions, especially in eukaryotes. Silent or synonymous codon positions, which do not determine amino acid sequences of the encoded proteins, define mRNA secondary structure and stability and affect the rate of translation, folding and post-translational modifications of nascent polypeptides. The RNA level selection is acting on synonymous sites in both prokaryotes and eukaryotes and is more common than previously thought. Selection pressure on the coding gene regions follows three-nucleotide periodic pattern of nucleotide base-pairing in mRNA, which is imposed by the genetic code. Synonymous positions of the coding regions have a higher level of hybridization potential relative to non-synonymous positions, and are multifunctional in their regulatory and structural roles. Recent experimental evidence and analysis of mRNA structure and interspecies conservation suggest that there is an evolutionary tradeoff between selective pressure acting at the RNA and protein levels. Here we provide a comprehensive overview of the studies that define the role of silent positions in regulating RNA structure and processing that exert downstream effects on proteins and their functions.

Functions of transmembrane domain 3 of human melanocortin-4 receptor

The melanocortin-4 receptor (MC4R) is a G protein-coupled receptor critical for maintaining energy homeostasis. Transmembrane domain 3 (TM3) of MC4R contains residues that were suggested to be essential in ligand binding and signaling. Several MC4R mutations in TM3 are associated with human obesity. To gain a better understanding of the functions of TM3, we analyzed the functions of 26 residues in TM3 using alanine-scanning mutagenesis. We showed that all mutants had normal cell-surface expression. Four mutants were defective in ligand binding and signaling and six mutants had normal ligand binding but impaired cAMP production. L140A had increased basal cAMP level. To further characterize the function of L140, we generated 17 additional L140 mutants. Fifteen L140 mutants had significantly decreased cell-surface expression, with L140R and L140V expressed normally. Ten L140 mutants had increased basal cAMP activities. Four L140 mutants were defective in ligand-stimulated cAMP generation. Interestingly, with the ERK1/2 pathway, we showed that nine constitutively active mutants had similar levels of basal pERK1/2 as that of WT, and two signaling defective mutants had similar levels of pERK1/2 as that of WT upon agonist stimulation, different from their cAMP signaling properties, suggesting biased signaling in these mutant receptors. In summary, we identified 13 residues in TM3 that were essential for ligand binding and/or signaling. Moreover, L140 was critical for locking MC4R in inactive conformation and several mutants showed biased signaling in cAMP and ERK1/2 signaling pathways.

Sensitive measurement of single-nucleotide polymorphism-induced changes of RNA conformation: application to disease studies

Single-nucleotide polymorphisms (SNPs) are often linked to critical phenotypes such as diseases or responses to vaccines, medications and environmental factors. However, the specific molecular mechanisms by which a causal SNP acts is usually not obvious. Changes in RNA secondary structure emerge as a possible explanation necessitating the development of methods to measure the impact of single-nucleotide variation on RNA structure. Despite the recognition of the importance of considering the changes in Boltzmann ensemble of RNA conformers in this context, a formal method to perform directly such comparison was lacking. Here, we solved this problem and designed an efficient method to compute the relative entropy between the Boltzmann ensembles of the native and a mutant structure. On the basis of this theoretical progress, we developed a software tool, remuRNA, and investigated examples of its application. Comparing the impact of common SNPs naturally occurring in populations with the impact of random point mutations, we found that structural changes introduced by common SNPs are smaller than those introduced by random point mutations. This suggests a natural selection against mutations that significantly change RNA structure and demonstrates, surprisingly, that randomly inserted point mutations provide inadequate estimation of random mutations effects. Subsequently, we applied remuRNA to determine which of the disease-associated non-coding SNPs are potentially related to RNA structural changes.

Functional Rescue of F508del-CFTR Using Small Molecule Correctors

High-throughput screens for small molecules that are effective in “correcting” the functional expression of F508del-CFTR have yielded several promising hits. Two such compounds are currently in clinical trial. Despite this success, it is clear that further advances will be required in order to restore 50% or greater of wild-type CFTR function to the airways of patients harboring the F508del-CFTR protein. Progress will be enhanced by our better understanding of the molecular and cellular defects caused by the F508del mutation, present in 90% of CF patients. The goal of this chapter is to review the current understanding of defects caused by F508del in the CFTR protein and in CFTR-mediated interactions important for its biosynthesis, trafficking, channel function, and stability at the cell surface. Finally, we will discuss the gaps in our knowledge regarding the mechanism of action of existing correctors, the unmet need to discover compounds which restore proper CFTR structure and function in CF affected tissues and new strategies for therapy development.

Functional study of risk loci of stem cell-associated gene lin-28B and associations with disease survival outcomes in epithelial ovarian cancer

Several single-nucleotide polymorphisms (SNPs) of the stem cell-associated gene lin-28B have been identified in association with ovarian cancer and ovarian cancer-related risk factors. However, whether these SNPs are functional or might be potential biomarkers for ovarian cancer prognosis remains unknown. The purposes of this study were to investigate the functional relevance of the identified lin-28B SNPs, as well as the associations of genotype and phenotype with epithelial ovarian cancer (EOC) survival. We analyzed five SNPs and mRNA levels of lin-28B in 211 primary EOC tissues using Taqman(®) SNP genotyping assays and SYBR green-based real-time PCR, respectively. The RNA secondary structures at the region of a genome-wide association-identified intronic rs314276 were analyzed theoretically with mfold and experimentally with circular dichroism spectroscopy. We found that rs314276 was a cis-acting expression quantitative trait locus (eQTL) in both additive and dominant models, while rs7759938 and rs314277 were significant or of borderline significance in dominant models only. The rs314276 variant significantly affects RNA secondary structure. No SNPs alone were associated with patient survival. However, we found that among patients initially responding to chemotherapy, those with higher lin-28B expression had higher mortality risk (hazard ratio =3.27, 95% confidence interval: 1.63-6.56) and relapse risk (hazard ratio = 2.53, 95% confidence interval: 1.41-4.54) than those with lower expression, and these associations remained in multivariate analyses. These results suggest that rs314276 alters RNA secondary structure and thereby influences gene expression, and that lin-28B is a cancer stem cell-associated marker, which may be a pharmaceutical target in the management of EOC.

Exon 10 CFTR gene mutation in male infertility

BACKGROUND

About 10% of infertilities with obstructive azoospermia are congenital and caused by CF gene mutations. M469I mutation was observed for the first time in Taiwanese patients. This mutation not only causes CF, but also may be the origin of infertility too.

OBJECTIVE

In this study, we aimed in designing a rapid, reliable RFLP-PCR procedure for detection of M469I mutation. The correlation and association between M469I mutation with infertility was investigated in this study.

MATERIALS AND METHODS

one hundred ten patients (90 non obstructive and 20 obstructive) and 60 normal individuals were considered in this study. M469I mutation was detected using RFLP-PCR. This technique was completely designed for M469I genotyping, for the first time in our study. Amplification of the region surrounding the mutation in exon 10 of CFTR gene was then performed. RFLP analysis was carried out using the NdeI restriction enzyme.

RESULTS

All genomic DNA samples were genotyped successfully. M469I mutation was observed only in patients group. Therefore, genotype containing mutant allele (GT) has been detected only in the patients group. There was no significant correlation between GT and TT genotypes with infertility (p=0.437).

CONCLUSION

The M469I mutation has only been observed in Exon 10 CFTR gene of infertile patients, not in the control group. This mutation causes congenital bilateral absence of vaz deferens and finally infertility. This indicates a strong association between the M469I mutation and male infertility. Therefore, this is a CF-causing CFTR mutation that could be considered as a cause of infertility.