Patterns of selection on synonymous and nonsynonymous variants in Drosophila miranda

We have investigated patterns of within-species polymorphism and between-species divergence for synonymous and nonsynonymous variants at a set of autosomal and X-linked loci of Drosophila miranda. D. pseudoobscura and D. affinis were used for the between-species comparisons. The results suggest the action of purifying selection on nonsynonymous, polymorphic variants. Among synonymous polymorphisms, there is a significant excess of synonymous mutations from preferred to unpreferred codons and of GC to AT mutations. There was no excess of GC to AT mutations among polymorphisms at noncoding sites. This suggests that selection is acting to maintain the use of preferred codons. Indirect evidence suggests that biased gene conversion in favor of GC base pairs may also be operating. The joint intensity of selection and biased gene conversion, in terms of the product of effective population size and the sum of the selection and conversion coefficients, was estimated to be approximately 0.65.     

Decoding the epigenetic effects of chromatin

A report on the Alan Wolffe EMBO Workshop on 'Chromatin and Epigenetics', Heidelberg, Germany, 19-22 June 2003.     

Species-specific codon context rules unveil non-neutrality effects of synonymous mutations

BackgroundCodon pair usage (codon context) is a species specific gene primary structure feature whose evolutionary and functional roles are poorly understood. The data available show that codon-context has direct impact on both translation accuracy and efficiency, but one does not yet understand how it affects these two translation variables or whether context biases shape gene evolution.ConclusionsSince in vivo studies provide evidence for a role of codon context on decoding fidelity in E. coli and for decoding efficiency in mammalian cells, our data support the hypothesis that, like codon usage, codon context modulates the evolution of gene primary structure and fine tunes the structure of open reading frames for high genome translational fidelity and efficiency in the 3 domains of life.     

Nucleotide composition and synonymous codon usage of open reading frames in Norovirus GII.4 variants

Abstract

Norovirus GII.4 variants, a genotype in genogroup II belonging to the genus Norovirus, is a single-strand positive sense RNA containing three open reading frames (ORF1, ORF2 and ORF3) and is the most important pathogen causing nonbacterial gastroenteritis outbreaks. By using bioinformatic softwares such as Codon W, SPSS and so on, a total of 292 strains of the viruses isolated from 1974 to 2016 were analyzed for nucleotide composition and synonymous codon usage in each ORF. The result shows that it is enriched for A over the other bases in nucleotide composition, G behind the other bases in the 3^rd site of all synonymous codons in the three ORFs. The patterns of nucleotide composition and codon bias of ORF2 are similar to those of ORF3 and different from those of ORF1. There are generally UpA motif and CpG motif in the codons with the lowest proportion. Correspondence analysis indicates that the codon usage may be changing over a certain time period for ORF1 in 2006 and 2012, ORF2 in 2012, and ORF3 in 2013. ENC (effective number of codons) plot and other analyses indicate that both natural selection and mutational pressure play partly roles in the ORFs, but natural selection is more important for ORF2 and ORF3. Besides, we also found all optimal codons in the ORFs. The study provides a basic understanding of the mechanism for norovirus GII.4 codon usage bias. Abbreviations ORF Open Reading Frame

ENC Effective Number of Codons

COA correspondence analysis

RSCU Relative Synonymous Codon Usage

CAI Codon Adaptation Index

CBI Codon Bias Index

Fop frequency of optimal codons

L_sym number of synonymous codons

L_aa length amino acids

GRAVY grand average of hydropathicity

Aroma aromaticity

Communicated by Ramaswamy H. Sarma     

The fitness effects of synonymous mutations in DNA and RNA viruses

Despite being silent with respect to protein sequence, synonymous nucleotide substitutions can be targeted by natural selection directly at the DNA or RNA level. However, there has been no systematic assessment of how frequent this type of selection is. Here, we have constructed 53 single random synonymous substitution mutants of the bacteriophages Qβ and ΦX174 by site-directed mutagenesis and assayed their fitness. Analysis of this mutant collection and of previous studies undertaken with a variety of single-stranded (ss) viruses demonstrates that selection at synonymous sites is stronger in RNA viruses than in DNA viruses. We estimate that this type of selection contributes approximately 18% of the overall mutational fitness effects in ssRNA viruses under our assay conditions and that random synonymous substitutions have a 5% chance of being lethal to the virus, whereas in ssDNA viruses, these figures drop to 1.4% and 0%, respectively. In contrast, the effects of nonsynonymous substitutions appear to be similar in ssRNA and ssDNA viruses.     

Rare missense and synonymous variants in UBE1 are associated with X-linked infantile spinal muscular atrophy

X-linked infantile spinal muscular atrophy (XL-SMA) is an X-linked disorder presenting with the clinical features hypotonia, areflexia, and multiple congenital contractures (arthrogryposis) associated with loss of anterior horn cells and infantile death. To identify the XL-SMA disease gene, we performed large-scale mutation analysis in genes located between markers DXS8080 and DXS7132 (Xp11.3–Xq11.1). This resulted in detection of three rare novel variants in exon 15 of UBE1 that segregate with disease: two missense mutations (c.1617 G→T, p.Met539Ile; c.1639 A→G, p.Ser547Gly) present each in one XL-SMA family, and one synonymous C→T substitution (c.1731 C→T, p.Asn577Asn) identified in another three unrelated families. Absence of the missense mutations was demonstrated for 3550 and absence of the synonymous mutation was shown in 7914 control X chromosomes; therefore, these results yielded statistical significant evidence for the association of the synonymous substitution and the two missense mutations with XL-SMA (p = 2.416 × 10⁻¹⁰, p = 0.001815). We also demonstrated that the synonymous C→T substitution leads to significant reduction of UBE1 expression and alters the methylation pattern of exon 15, implying a plausible role of this DNA element in developmental UBE1 expression in humans. Our observations indicate first that XL-SMA is part of a growing list of neurodegenerative disorders associated with defects in the ubiquitin-proteasome pathway and second that synonymous C→T transitions might have the potential to affect gene expression.     

Characterization of coding synonymous and non-synonymous variants in ADAMTS13 using ex vivo and in silico approaches

Synonymous variations, which are defined as codon substitutions that do not change the encoded amino acid, were previously thought to have no effect on the properties of the synthesized protein(s). However, mounting evidence shows that these "silent" variations can have a significant impact on protein expression and function and should no longer be considered "silent". Here, the effects of six synonymous and six non-synonymous variations, previously found in the gene of ADAMTS13, the von Willebrand Factor (VWF) cleaving hemostatic protease, have been investigated using a variety of approaches. The ADAMTS13 mRNA and protein expression levels, as well as the conformation and activity of the variants have been compared to that of wild-type ADAMTS13. Interestingly, not only the non-synonymous variants but also the synonymous variants have been found to change the protein expression levels, conformation and function. Bioinformatic analysis of ADAMTS13 mRNA structure, amino acid conservation and codon usage allowed us to establish correlations between mRNA stability, RSCU, and intracellular protein expression. This study demonstrates that variants and more specifically, synonymous variants can have a substantial and definite effect on ADAMTS13 function and that bioinformatic analysis may allow development of predictive tools to identify variants that will have significant effects on the encoded protein.     

Decoding variants in drug-metabolizing enzymes and transporters in solid tumor patients by whole-exome sequencing

BackgroundPharmacogenetics is involved in customizing therapy according to the genetic makeup of an individual, and is applicable for chemotherapy, radiotherapy as well as targeted therapy. Drug metabolizing enzymes (DMEs) involving both phase I, and phase II reactions are widely studied. Our study was involved in whole exome sequencing (WES) of cancer patients, followed by analysis for identifying key variations in DMEs, and associated transporters that have a potential impact on treatment outcome.MethodologyA total of 181 solid tumor patients at stage >/= III were subjected to WES by the SureSelect^XT Human All Exon V6 + UTR library preparation kit, and sequencing in the Illumina NextSeq 550 system. Bioinformatics analysis involved use of GATK pipeline, and the variants were further assessed for population frequency, functional impact with annovar insilico algorithms. Further variant information from significant DMEs, and transporters were extracted and analyzed with PharmGKB to assess level of evidence and infer their impact on the pathways involved in drug response.ResultsThe total study cohort of 181 solid tumor patients included 60 males, and 121 females respectively. Among DMEs, deleterious mutation in dihydropyrimidine dehydrogenase (DPYD; rs67376798), solute carrier organic anion transporter family member 1B1 (SLCO1B1*5), and cytochrome P450 2D6 (CYP2D6*10) associated with metabolism of anticancer drugs was detected to be in high frequency of 26%, 21% and 25% respectively.ConclusionOur analysis detected variations in both phase I and phase II DMEs, as well as associated transporter genes which has been documented to reduce drug efficacy, as well as cause grade 3 and 4 toxicity. Our study reiterates the significance of pharmacogenomics in stratifying patients for appropriate therapy regimen focused at better treatment outcome and quality of life.     

Variation in synonymous substitution rates among mammalian genes and the correlation between synonymous and nonsynonymous divergences

Using mammalian gene sequences, the variances in the numbers of synonymous and nonsynonymous substitutions among genes were estimated together with the correlation coefficient between the two. The expected correlation coefficient can be obtained under the neutral theory using these estimated values of the variances. The expected coefficient is found to often be one-half to two-thirds of the observed value. Possible causes for the disagreement were discussed, such as correlated selective constraints on the two types of substitutions and excess doublet mutations. The variance of mutation rate and that of selective constraint were also estimated. The results show that the coefficient of variation of the former is 0.2–0.3, whereas that of the latter is 0.7–0.9. Correspondence to: T. Ohta     

Functional effects of variants of the RNA chaperone Hfq

The ring-shaped RNA chaperone Hfq has recently received much attention owing to its multiple roles in RNA metabolism. In this study we have performed a mutational analysis of the Escherichia coli hfq gene, and have studied the effects of amino acid substitutions at several positions in the Hfq protein as well as of C-terminal truncations on its role in phage Qbeta replication, in repression of a target mRNA, and on the stability of the small regulatory RNA DsrA. These functional studies provided insights into the interaction of Hfq with RNA and suggested a role for the C-terminus of Hfq in DsrA stability.     

Compensating effects on the cytotoxicity of ribonuclease A variants

Ribonuclease (RNase) A can be endowed with cytotoxic activity by enabling it to evade the cytosolic ribonuclease inhibitor protein (RI). Enhancing its conformational stability can increase further its cytotoxicity. Herein, the A4C/K41R/G88R/V118C variant of RNase A was created to integrate four individual changes that greatly decrease RI affinity (K41R/G88R) and increase conformational stability (A4C/V118C). Yet, the variant suffers a decrease in ribonucleolytic activity and is only as potent a cytotoxin as its precursors. Thus, individual changes that increase cytotoxicity can have offsetting consequences. Overall, cytotoxicity correlates well with the maintenance of ribonucleolytic activity in the presence of RI. The parameter (k(cat)/K(m))(cyto), which reports on the ability of a ribonuclease to manifest its ribonucleolytic activity in the cytosol, is especially useful in predicting the cytotoxicity of an RNase A variant.     

Negative feedback buffers effects of regulatory variants

Mechanisms conferring robustness against regulatory variants have been controversial. Previous studies suggested widespread buffering of RNA misexpression on protein levels during translation. We do not find evidence that translational buffering is common. Instead, we find extensive buffering at the level of RNA expression, exerted through negative feedback regulation acting in trans, which reduces the effect of regulatory variants on gene expression. Our approach is based on a novel experimental design in which allelic differential expression in a yeast hybrid strain is compared to allelic differential expression in a pool of its spores. Allelic differential expression in the hybrid is due to cis-regulatory differences only. Instead, in the pool of spores allelic differential expression is not only due to cis-regulatory differences but also due to local trans effects that include negative feedback. We found that buffering through such local trans regulation is widespread, typically compensating for about 15% of cis-regulatory effects on individual genes. Negative feedback is stronger not only for essential genes, indicating its functional relevance, but also for genes with low to middle levels of expression, for which tight regulation matters most. We suggest that negative feedback is one mechanism of Waddington''s canalization, facilitating the accumulation of genetic variants that might give selective advantage in different environments.     

Site-to-site variation of synonymous substitution rates

We develop a new model for studying the molecular evolution of protein-coding DNA sequences. In contrast to existing models, we incorporate the potential for site-to-site heterogeneity of both synonymous and nonsynonymous substitution rates. We demonstrate that within-gene heterogeneity of synonymous substitution rates appears to be common. Using the new family of models, we investigate the utility of a variety of new statistical inference procedures, and we pay particular attention to issues surrounding the detection of sites undergoing positive selection. We discuss how failure to model synonymous rate variation in the model can lead to misidentification of sites as positively selected.     

Understanding the contribution of synonymous mutations to human disease

Synonymous mutations - sometimes called 'silent' mutations - are now widely acknowledged to be able to cause changes in protein expression, conformation and function. The recent increase in knowledge about the association of genetic variants with disease, particularly through genome-wide association studies, has revealed a substantial contribution of synonymous SNPs to human disease risk and other complex traits. Here we review current understanding of the extent to which synonymous mutations influence disease, the various molecular mechanisms that underlie these effects and the implications for future research and biomedical applications.     

Forbidden synonymous substitutions in coding regions

In the evolution of highly conserved genes, a few "synonymous" substitutions at third bases that would not alter the protein sequence are forbidden or very rare, presumably as a result of functional requirements of the gene or the messenger RNA. Another 10% or 20% of codons are significantly less variable by synonymous substitution than are the majority of codons. The changes that occur at the majority of third bases are subject to codon usage restrictions. These usage restrictions control sequence similarities between very distant genes. For example, 70% of third bases are identical in calmodulin genes of man and trypanosome. Third-base similarities of distant genes for conserved proteins are mathematically predicted, on the basis of the G+C composition of third bases. These observations indicate the need for reexamination of methods used to calculate synonymous substitutions.