Synonymous edits in the Escherichia coli genome have substantial and condition-dependent effects on fitness.

CRISPR-Cas-based genome editing is widely used in bacteria at scales ranging from construction of individual mutants to massively parallel libraries. This procedure relies on guide RNA-directed cleavage of the genome followed by repair with a template that introduces a desired mutation along with synonymous "immunizing" mutations to prevent re-cleavage of the genome after editing. Because the immunizing mutations do not change the protein sequence, they are often assumed to be neutral. However, synonymous mutations can change mRNA structures in ways that alter levels of the encoded proteins. We have tested the assumption that immunizing mutations are neutral by constructing a library of over 50,000 edits that consist of only synonymous mutations in Escherichia coli. Thousands of edits had substantial effects on fitness during growth of E. coli on acetate, a poor carbon source that is toxic at high concentrations. The percentage of high-impact edits varied considerably between genes and at different positions within genes. We reconstructed clones with high-impact edits and found that 69% indeed had significant effects on growth in acetate. Interestingly, fewer edits affected fitness during growth in glucose, a preferred carbon source, suggesting that changes in protein expression caused by synonymous mutations may be most important when an organism encounters challenging conditions. Finally, we showed that synonymous edits can have widespread effects; a synonymous edit at the 5' end of ptsI altered expression of hundreds of genes. Our results suggest that the synonymous immunizing edits introduced during CRISPR-Cas-based genome editing should not be assumed to be innocuous.

Recurrent emergence of Klebsiella pneumoniae carbapenem resistance mediated by an inhibitory ompK36 mRNA secondary structure.

Outer membrane porins in Gram-negative bacteria facilitate antibiotic influx. In Klebsiella pneumoniae, modifications in the porin OmpK36 are implicated in increasing resistance to carbapenems. An analysis of large K. pneumoniae genome collections, encompassing major healthcare-associated clones, revealed the recurrent emergence of a synonymous cytosine-to-thymine transition at position 25 (25c > t) in ompK36. We show that the 25c > t transition increases carbapenem resistance through depletion of OmpK36 from the outer membrane. The mutation attenuates K. pneumoniae in a murine pneumonia model, which accounts for its limited clonal expansion observed by phylogenetic analysis. However, in the context of carbapenem treatment, the 25c > t transition tips the balance toward treatment failure, thus accounting for its recurrent emergence. Mechanistically, the 25c > t transition mediates an intramolecular messenger RNA (mRNA) interaction between a uracil encoded by 25t and the first adenine within the Shine-Dalgarno sequence. This specific interaction leads to the formation of an RNA stem structure, which obscures the ribosomal binding site thus disrupting translation. While mutations reducing OmpK36 expression via transcriptional silencing are known, we uniquely demonstrate the repeated selection of a synonymous ompK36 mutation mediating translational suppression in response to antibiotic pressure.

The DNA damage sensor ATM kinase interacts with the p53 mRNA and guides the DNA damage response pathway.

BACKGROUND: The ATM kinase constitutes a master regulatory hub of DNA damage and activates the p53 response pathway by phosphorylating the MDM2 protein, which develops an affinity for the p53 mRNA secondary structure. Disruption of this interaction prevents the activation of the nascent p53. The link of the MDM2 protein-p53 mRNA interaction with the upstream DNA damage sensor ATM kinase and the role of the p53 mRNA in the DNA damage sensing mechanism, are still highly anticipated. METHODS: The proximity ligation assay (PLA) has been extensively used to reveal the sub-cellular localisation of the protein-mRNA and protein-protein interactions. ELISA and co-immunoprecipitation confirmed the interactions in vitro and in cells. RESULTS: This study provides a novel mechanism whereby the p53 mRNA interacts with the ATM kinase enzyme and shows that the L22L synonymous mutant, known to alter the secondary structure of the p53 mRNA, prevents the interaction. The relevant mechanistic roles in the DNA Damage Sensing pathway, which is linked to downstream DNA damage response, are explored. Following DNA damage (double-stranded DNA breaks activating ATM), activated MDMX protein competes the ATM-p53 mRNA interaction and prevents the association of the p53 mRNA with NBS1 (MRN complex). These data also reveal the binding domains and the phosphorylation events on ATM that regulate the interaction and the trafficking of the complex to the cytoplasm. CONCLUSION: The presented model shows a novel interaction of ATM with the p53 mRNA and describes the link between DNA Damage Sensing with the downstream p53 activation pathways; supporting the rising functional implications of synonymous mutations altering secondary mRNA structures.

Genetic analysis and molecular basis of G6PD deficiency among malaria patients in Thailand: implications for safe use of 8-aminoquinolines.

BACKGROUND: It was hypothesized that glucose-6-phosphate dehydrogenase (G6PD) deficiency confers a protective effect against malaria infection, however, safety concerns have been raised regarding haemolytic toxicity caused by radical cure with 8-aminoquinolines in G6PD-deficient individuals. Malaria elimination and control are also complicated by the high prevalence of G6PD deficiency in malaria-endemic areas. Hence, accurate identification of G6PD deficiency is required to identify those who are eligible for malaria treatment using 8-aminoquinolines. METHODS: The prevalence of G6PD deficiency among 408 Thai participants diagnosed with malaria by microscopy (71), and malaria-negative controls (337), was assessed using a phenotypic test based on water-soluble tetrazolium salts. High-resolution melting (HRM) curve analysis was developed from a previous study to enable the detection of 15 common missense, synonymous and intronic G6PD mutations in Asian populations. The identified mutations were subjected to biochemical and structural characterisation to understand the molecular mechanisms underlying enzyme deficiency. RESULTS: Based on phenotypic testing, the prevalence of G6PD deficiency (< 30% activity) was 6.13% (25/408) and intermediate deficiency (30-70% activity) was found in 15.20% (62/408) of participants. Several G6PD genotypes with newly discovered double missense variants were identified by HRM assays, including G6PD Gaohe + Viangchan, G6PD Valladolid + Viangchan and G6PD Canton + Viangchan. A significantly high frequency of synonymous (c.1311C>T) and intronic (c.1365-13T>C and c.486-34delT) mutations was detected with intermediate to normal enzyme activity. The double missense mutations were less catalytically active than their corresponding single missense mutations, resulting in severe enzyme deficiency. While the mutations had a minor effect on binding affinity, structural instability was a key contributor to the enzyme deficiency observed in G6PD-deficient individuals. CONCLUSIONS: With varying degrees of enzyme deficiency, G6PD genotyping can be used as a complement to phenotypic screening to identify those who are eligible for 8-aminoquinolines. The information gained from this study could be useful for management and treatment of malaria, as well as for the prevention of unanticipated reactions to certain medications and foods in the studied population.

Non-synonymous, synonymous, and non-coding nucleotide variants contribute to recurrently altered biological processes during retinoblastoma progression.

Retinoblastomas form in response to biallelic RB1 mutations or MYCN amplification and progress to more aggressive and therapy-resistant phenotypes through accumulation of secondary genomic changes. Progression-related changes include recurrent somatic copy number alterations and typically non-recurrent nucleotide variants, including synonymous and non-coding variants, whose significance has been unclear. To determine if nucleotide variants recurrently affect specific biological processes, we identified altered genes and over-represented variant gene ontologies in 168 exome or whole-genome-sequenced retinoblastomas and 12 tumor-matched cell lines. In addition to RB1 mutations, MYCN amplification, and established retinoblastoma somatic copy number alterations, the analyses revealed enrichment of variant genes related to diverse biological processes including histone monoubiquitination, mRNA processing (P) body assembly, and mitotic sister chromatid segregation and cytokinesis. Importantly, non-coding and synonymous variants increased the enrichment significance of each over-represented biological process term. To assess the effects of such mutations, we examined the consequences of a 3' UTR variant of PCGF3 (a BCOR-binding component of Polycomb repressive complex I), dual 3' UTR variants of CDC14B (a regulator of sister chromatid segregation), and a synonymous variant of DYNC1H1 (a regulator of P-body assembly). One PCGF3 and one of two CDC14B 3' UTR variants impaired gene expression whereas a base-edited DYNC1H1 synonymous variant altered protease sensitivity and stability. Retinoblastoma cell lines retained only ~50% of variants detected in tumors and enriched for new variants affecting p53 signaling. These findings reveal potentially important differences in retinoblastoma cell lines and tumors and implicate synonymous and non-coding variants, along with non-synonymous variants, in retinoblastoma oncogenesis.

Incorporating mutational heterogeneity to identify genes that are enriched for synonymous mutations in cancer.

BACKGROUND: Synonymous mutations, which change the DNA sequence but not the encoded protein sequence, can affect protein structure and function, mRNA maturation, and mRNA half-lives. The possibility that synonymous mutations might be enriched in cancer has been explored in several recent studies. However, none of these studies control for all three types of mutational heterogeneity (patient, histology, and gene) that are known to affect the accurate identification of non-synonymous cancer-associated genes. Our goal is to adopt the current standard for non-synonymous mutations in an investigation of synonymous mutations. RESULTS: Here, we create an algorithm, MutSigCVsyn, an adaptation of MutSigCV, to identify cancer-associated genes that are enriched for synonymous mutations based on a non-coding background model that takes into account the mutational heterogeneity across these levels. Using MutSigCVsyn, we first analyzed 2572 cancer whole-genome samples from the Pan-cancer Analysis of Whole Genomes (PCAWG) to identify non-synonymous cancer drivers as a quality control. Indicative of the algorithm accuracy we find that 58.6% of these candidate genes were also found in Cancer Census Gene (CGC) list, and 66.2% were found within the PCAWG cancer driver list. We then applied it to identify 30 putative cancer-associated genes that are enriched for synonymous mutations within the same samples. One of the promising gene candidates is the B cell lymphoma 2 (BCL-2) gene. BCL-2 regulates apoptosis by antagonizing the action of proapoptotic BCL-2 family member proteins. The synonymous mutations in BCL2 are enriched in its anti-apoptotic domain and likely play a role in cancer cell proliferation. CONCLUSION: Our study introduces MutSigCVsyn, an algorithm that accounts for mutational heterogeneity at patient, histology, and gene levels, to identify cancer-associated genes that are enriched for synonymous mutations using whole genome sequencing data. We identified 30 putative candidate genes that will benefit from future experimental studies on the role of synonymous mutations in cancer biology.

Natural Selection Shapes Codon Usage in the Human Genome.

Synonymous codon usage has been identified as a determinant of translational efficiency and mRNA stability in model organisms and human cell lines. However, whether natural selection shapes human codon content to optimize translation efficiency is unclear. Furthermore, aside from those that affect splicing, synonymous mutations are typically ignored as potential contributors to disease. Using genetic sequencing data from nearly 200,000 individuals, we uncover clear evidence that natural selection optimizes codon content in the human genome. In deriving intolerance metrics to quantify gene-level constraint on synonymous variation, we discover that dosage-sensitive genes, DNA-damage-response genes, and cell-cycle-regulated genes are particularly intolerant to synonymous variation. Notably, we illustrate that reductions in codon optimality in BRCA1 can attenuate its function. Our results reveal that synonymous mutations most likely play an underappreciated role in human variation.

Non-random Codon Usage of Synonymous and Non-synonymous Mutations in the Human HLA-A Gene.

The structure and function of human leucocyte antigen (HLA-A) is well known and is an extremely variable protein. From the public HLA-A database, we chose 26 high frequency HLA-A alleles (45% of sequenced alleles). Using five arbitrary references from these alleles, we analyzed synonymous mutations at the third codon position (sSNP3) and non-synonymous mutations (NSM). Both mutation types showed non-random locations of 29 sSNP3 codons and 71 NSM codons in the five reference lists. Most sSNP3 codons show identical mutation types with many mutations resulting from cytosine deamination. We proposed 23 ancestral parents of sSNP3 in five reference sequences using conserved parents in five unidirectional codons and 18 majority parents in reciprocal codons. These 23 proposed ancestral parents show exclusive codon usage of G3 or C3 parents located on both DNA strands that mutate to A3 or T3 variants mostly (76%) by cytosine deamination The sSNP3 and NSM show clear separation of the two variant types with most sSNP3 located in conserved areas in exons 2, 3 and 4, compared to most NSM appearing in two Variable Areas with no sSNP3 in the latter parts of exons 2 (α1) and 3 (α2). The Variable Areas contain NSM (polymorphic) residues at the center of the groove that bind the foreign peptide. We find distinctly different mutation patterns in NSM codons from those of sSNP3. Namely, G-C to A-T mutation frequency was much smaller, suggesting that evolutional pressures of deamination and other mechanisms applied to the two areas are significantly different.

Evidence for Strong Mutation Bias toward, and Selection against, U Content in SARS-CoV-2: Implications for Vaccine Design.

Large-scale re-engineering of synonymous sites is a promising strategy to generate vaccines either through synthesis of attenuated viruses or via codon-optimized genes in DNA vaccines. Attenuation typically relies on deoptimization of codon pairs and maximization of CpG dinucleotide frequencies. So as to formulate evolutionarily informed attenuation strategies that aim to force nucleotide usage against the direction favored by selection, here, we examine available whole-genome sequences of SARS-CoV-2 to infer patterns of mutation and selection on synonymous sites. Analysis of mutational profiles indicates a strong mutation bias toward U. In turn, analysis of observed synonymous site composition implicates selection against U. Accounting for dinucleotide effects reinforces this conclusion, observed UU content being a quarter of that expected under neutrality. Possible mechanisms of selection against U mutations include selection for higher expression, for high mRNA stability or lower immunogenicity of viral genes. Consistent with gene-specific selection against CpG dinucleotides, we observe systematic differences of CpG content between SARS-CoV-2 genes. We propose an evolutionarily informed approach to attenuation that, unusually, seeks to increase usage of the already most common synonymous codons. Comparable analysis of H1N1 and Ebola finds that GC3 deviated from neutral equilibrium is not a universal feature, cautioning against generalization of results.

The Heterogeneous Landscape and Early Evolution of Pathogen-Associated CpG Dinucleotides in SARS-CoV-2.

COVID-19 can lead to acute respiratory syndrome, which can be due to dysregulated immune signaling. We analyze the distribution of CpG dinucleotides, a pathogen-associated molecular pattern, in the SARS-CoV-2 genome. We characterize CpG content by a CpG force that accounts for statistical constraints acting on the genome at the nucleotidic and amino acid levels. The CpG force, as the CpG content, is overall low compared with other pathogenic betacoronaviruses; however, it widely fluctuates along the genome, with a particularly low value, comparable with the circulating seasonal HKU1, in the spike coding region and a greater value, comparable with SARS and MERS, in the highly expressed nucleocapside coding region (N ORF), whose transcripts are relatively abundant in the cytoplasm of infected cells and present in the 3'UTRs of all subgenomic RNA. This dual nature of CpG content could confer to SARS-CoV-2 the ability to avoid triggering pattern recognition receptors upon entry, while eliciting a stronger response during replication. We then investigate the evolution of synonymous mutations since the outbreak of the COVID-19 pandemic, finding a signature of CpG loss in regions with a greater CpG force. Sequence motifs preceding the CpG-loss-associated loci in the N ORF match recently identified binding patterns of the zinc finger antiviral protein. Using a model of the viral gene evolution under human host pressure, we find that synonymous mutations seem driven in the SARS-CoV-2 genome, and particularly in the N ORF, by the viral codon bias, the transition-transversion bias, and the pressure to lower CpG content.

Synonymous mutations that regulate translation speed might play a non-negligible role in liver cancer development.

BACKGROUND: Synonymous mutations do not change the protein sequences. Automatically, they have been regarded as neutral events and are ignored in the mutation-based cancer studies. However, synonymous mutations will change the codon optimality, resulting in altered translational velocity. METHODS: We fully utilized the transcriptome and translatome of liver cancer and normal tissue from ten patients. We profiled the mutation spectrum and examined the effect of synonymous mutations on translational velocity. RESULTS: Synonymous mutations that increase the codon optimality significantly enhanced the translational velocity, and were enriched in oncogenes. Meanwhile, synonymous mutations decreasing codon optimality slowed down translation, and were enriched in tumor suppressor genes. These synonymous mutations significantly contributed to the translational changes in tumor samples compared to normal samples. CONCLUSIONS: Synonymous mutations might play a role in liver cancer development by altering codon optimality and translational velocity. Synonymous mutations should no longer be ignored in the genome-wide studies.

SARS-CoV2 envelope protein: non-synonymous mutations and its consequences.

In the NCBI database, as on June 6, 2020, total number of available complete genome sequences of SARS-CoV2 across the world is 3617. The envelope (E) protein of SARS-CoV2 possesses several non-synonymous mutations over the transmembrane and C-terminus domains in 15 (0.414%) genomes among 3617 SARS-CoV2 genomes, analyzed. More precisely, 10(0.386%) out of 2588 genomes from the USA, 3(0.806%) from Asia, 1 (0.348%) from Europe and 1 (0.274%) from Oceania contained the missense mutations over the E-protein of SARS-CoV2 genomes. The C-terminus motif DLLV has been to DFLV and YLLV in the proteins from QJR88103 (Australia: Victoria) and QKI36831 (China: Guangzhou) respectively, which might affect the binding of this motif with the host protein PALS1.

Trends of mutation accumulation across global SARS-CoV-2 genomes: Implications for the evolution of the novel coronavirus.

To understand SARS-CoV-2 microevolution, this study explored the genome-wide frequency, gene-wise distribution, and molecular nature of all point-mutations detected across its 71,703 RNA-genomes deposited in GISAID till 21 August 2020. Globally, nsp1/nsp2 and orf7a/orf3a were the most mutation-ridden non-structural and structural genes respectively. Phylogeny of 4618 spatiotemporally-representative genomes revealed that entities belonging to the early lineages are mostly spread over Asian countries, including India, whereas the recently-derived lineages are more globally distributed. Of the total 20,163 instances of polymorphism detected across global genomes, 12,594 and 7569 involved transitions and transversions, predominated by cytidine-to-uridine and guanosine-to-uridine conversions, respectively. Positive selection of nonsynonymous mutations (dN/dS >1) in most of the structural, but not the non-structural, genes indicated that SARS-CoV-2 has already harmonized its replication/transcription machineries with the host metabolism, while it is still redefining virulence/transmissibility strategies at the molecular level. Mechanistic bases and evolutionary/pathogenicity-related implications are discussed for the predominant mutation-types.

Protein Engineering Approaches to Enhance Fungal Laccase Production in S. cerevisiae.

Laccases secreted by saprotrophic basidiomycete fungi are versatile biocatalysts able to oxidize a wide range of aromatic compounds using oxygen as the sole requirement. Saccharomyces cerevisiae is a preferred host for engineering fungal laccases. To assist the difficult secretion of active enzymes by yeast, the native signal peptide is usually replaced by the preproleader of S. cerevisiae alfa mating factor (MFα1). However, in most cases, only basal enzyme levels are obtained. During directed evolution in S. cerevisiae of laccases fused to the α-factor preproleader, we demonstrated that mutations accumulated in the signal peptide notably raised enzyme secretion. Here we describe different protein engineering approaches carried out to enhance the laccase activity detected in the liquid extracts of S. cerevisiae cultures. We demonstrate the improved secretion of native and engineered laccases by using the fittest mutated α-factor preproleader obtained through successive laccase evolution campaigns in our lab. Special attention is also paid to the role of protein N-glycosylation in laccase production and properties, and to the introduction of conserved amino acids through consensus design enabling the expression of certain laccases otherwise not produced by the yeast. Finally, we revise the contribution of mutations accumulated in laccase coding sequence (CDS) during previous directed evolution campaigns that facilitate enzyme production.

Genome-wide analysis on the maize genome reveals weak selection on synonymous mutations.

BACKGROUND: Synonymous mutations are able to change the tAI (tRNA adaptation index) of a codon and consequently affect the local translation rate. Intuitively, one may hypothesize that those synonymous mutations which increase the tAI values are favored by natural selection. RESULTS: We use the maize (Zea mays) genome to test our assumption. The first supporting evidence is that the tAI-increasing synonymous mutations have higher fixed-to-polymorphic ratios than the tAI-decreasing ones. Next, the DAF (derived allele frequency) or MAF (minor allele frequency) of the former is significantly higher than the latter. Moreover, similar results are obtained when we investigate CAI (codon adaptation index) instead of tAI. CONCLUSION: The synonymous mutations in the maize genome are not strictly neutral. The tAI-increasing mutations are positively selected while those tAI-decreasing ones undergo purifying selection. This selection force might be weak but should not be automatically ignored.

Causative and common PHOX2B variants define a broad phenotypic spectrum.

Paired Like homeobox 2B (PHOX2B) is a gene crucial for the differentiation of the neural lineages of the autonomic nervous system (ANS), whose coding mutations cause congenital central hypoventilation syndrome (CCHS). The vast majority of PHOX2B mutations in CCHS is represented by expansions of a polyalanine region in exon 3, collectively defined PARMs (PolyAlanine Repeat Mutations), the minority being frameshift, missense and nonsense mutations, defined as NPARMs (Non-PARMs). While PARMs are nearly exclusively associated with isolated CCHS, most of NPARMs is detected in syndromic CCHS, presenting with neuroblastoma and/or Hirschsprung disease. More recently, evidence of a complex role of PHOX2B in the pathogenesis of a wider spectrum of ANS disorders has emerged. Indeed, common and hypomorphic PHOX2B variants, including synonymous, polyalanine-contractions, gene deletions may influence the occurrence of either apparent life-threatening event (ALTE), Sudden Infant Death Syndrome (SIDS), neuroblastoma, or isolated HSCR, likely through small effects on PHOX2B expression levels. After an introduction to the role of PHOX2B in the ANS development, causative mutations, common variants, and gene expression deregulation of the PHOX2B gene are discussed, though the involvement of synonymous variants and contractions requires further confirmations with respect to ANS disorders and molecular mechanisms underlying the PHOX2B phenotypic heterogeneity.

Parsing the synonymous mutations in the maize genome: isoaccepting mutations are more advantageous in regions with codon co-occurrence bias.

BACKGROUND: Synonymous mutations do not change amino acids but do sometimes change the tRNAs (anticodons) that decode a particular codon. An isoaccepting codon is a synonymous codon that shares the same tRNA. If a mutated codon could base pair with the same anticodon as the original, the mutation is termed an isoaccepting mutation. An interesting but less-studied type of codon bias is codon co-occurrence bias. There is a trend to cluster the isoaccepting codons in the genome. The proposed advantage of codon co-occurrence bias is that the tRNA released from the ribosome E site could be quickly recharged and subsequently decode the following isoaccepting codons. This advantage would enhance translation efficiency. In plant species, whether there are signals of positive selection on isoaccepting mutations in the codon co-occurred regions has not been studied. RESULTS: We termed polymorphic mutations in coding regions using publicly available RNA-seq data in maize (Zea mays). Next, we classified all synonymous mutations into three categories according to the context, i.e., the relationship between the focal codon and the previous codon, as follows: isoaccepting, nonisoaccepting and nonsynonymous. We observed higher fractions of isoaccepting mutations in the isoaccepting context. If we looked at the minor allele frequency (MAF) spectrum, the isoaccepting mutations have a higher MAF in the isoaccepting context than that in other regions, and accordingly, the nonisoaccepting mutations have a higher MAF in the nonisoaccepting context. CONCLUSION: Our results indicate that in regions with codon co-occurrence bias, natural selection maintains this pattern by suppressing the nonisoaccepting mutations. However, if the consecutive codons are nonisoaccepting, mutations tend to switch these codons to become isoaccepting. Our study demonstrates that the codon co-occurrence bias in the maize genome is selectively maintained by natural selection and that the advantage of this trend could potentially be the rapid recharging and reuse of tRNAs to increase translation efficiency.

An analysis of mutational signatures of synonymous mutations across 15 cancer types.

BACKGROUND: Synonymous mutations have been identified to play important roles in cancer development, although they do not modify the protein sequences. However, relatively little research has specifically delineated the functionality of synonymous mutations in cancer. RESULTS: We investigated the nucleotide-based and amino acid-based features of synonymous mutations across 15 cancer types from The Cancer Genome Atlas (TCGA), and revealed novel driver candidates by identifying hotspot mutations. Firstly, synonymous mutations were analyzed between TCGA and 1000 Genomes Project at nucleotide and amino acid levels. We found that C:G → T:A transitions were the most frequent single-base substitutions, and leucine underwent the largest number of synonymous mutations in TCGA due to prevalent C → T transition, which induced the transformation between optimal and non-optimal codons. Next, 97 synonymous hotspot mutations in 86 genes were nominated as candidate drivers with potential cancer risk by considering the mutational rates across different sequence contexts. We observed that non-CpG-island GC transition sequence context was positively selected across most of cancer types, and different sequence contexts under which hotspot mutations occur could be significance for genetic differences and functional features. We also found that the hotspots were more conserved than neutral mutations of hotspot-mutation-containing-genes and frequently happened at leucine. In addition, we mapped hotspots, neutral and non-hotspot mutations of hotspot-mutation-containing-genes to their respective protein domains and found ion transport domain was the most frequent one, which could mediate the cell interaction and had relevant implication for tumor therapy. And the signatures of synonymous hotspots were qualitatively similar with those of harmful missense variants. CONCLUSIONS: We illustrated the preferences of cancer associated synonymous mutations, especially hotspots, and laid the groundwork for understanding the synonymous mutations act as drivers in cancer.

Generation and comparative genomics of synthetic dengue viruses.

BACKGROUND: Synthetic virology is an important multidisciplinary scientific field, with emerging applications in biotechnology and medicine, aiming at developing methods to generate and engineer synthetic viruses. In particular, many of the RNA viruses, including among others the Dengue and Zika, are widespread pathogens of significant importance to human health. The ability to design and synthesize such viruses may contribute to exploring novel approaches for developing vaccines and virus based therapies. RESULTS: Here we develop a full multidisciplinary pipeline for generation and analysis of synthetic RNA viruses and specifically apply it to Dengue virus serotype 2 (DENV-2). The major steps of the pipeline include comparative genomics of endogenous and synthetic viral strains. Specifically, we show that although the synthetic DENV-2 viruses were found to have lower nucleotide variability, their phenotype, as reflected in the study of the AG129 mouse model morbidity, RNA levels, and neutralization antibodies, is similar or even more pathogenic in comparison to the wildtype master strain. Additionally, the highly variable positions, identified in the analyzed DENV-2 population, were found to overlap with less conserved homologous positions in Zika virus and other Dengue serotypes. These results may suggest that synthetic DENV-2 could enhance virulence if the correct sequence is selected. CONCLUSIONS: The approach reported in this study can be used to generate and analyze synthetic RNA viruses both on genotypic and on phenotypic level. It could be applied for understanding the functionality and the fitness effects of any set of mutations in viral RNA and for editing RNA viruses for various target applications.

Codon-Resolution Analysis Reveals a Direct and Context-Dependent Impact of Individual Synonymous Mutations on mRNA Level.

Codon usage bias (CUB) refers to the observation that synonymous codons are not used equally frequently in a genome. CUB is stronger in more highly expressed genes, a phenomenon commonly explained by stronger natural selection on translational accuracy and/or efficiency among these genes. Nevertheless, this phenomenon could also occur if CUB regulates gene expression at the mRNA level, a hypothesis that has not been tested until recently. Here, we attempt to quantify the impact of synonymous mutations on mRNA level in yeast using 3,556 synonymous variants of a heterologous gene encoding green fluorescent protein (GFP) and 523 synonymous variants of an endogenous gene TDH3. We found that mRNA level was positively correlated with CUB among these synonymous variants, demonstrating a direct role of CUB in regulating transcript concentration, likely via regulating mRNA degradation rate, as our additional experiments suggested. More importantly, we quantified the effects of individual synonymous mutations on mRNA level and found them dependent on 1) CUB and 2) mRNA secondary structure, both in proximal sequence contexts. Our study reveals the pleiotropic effects of synonymous codon usage and provides an additional explanation for the well-known correlation between CUB and gene expression level.