Synonymous mutations make dramatic contributions to fitness when growth is limited by a weak-link enzyme

Cetuximab chemotherapy resistance: Insight into the homeostatic evolution of head and neck cancer (Review)

The complex evolution of genetic alterations in cancer that occurs in vivo is a selective process involving numerous factors and mechanisms. Chemotherapeutic agents that prevent the growth and spread of cancer cells induce selective pressure, leading to rapid artificial selection of resistant subclones. This rapid evolution is possible because antineoplastic drugs promote alterations in tumor‑cell metabolism, thus creating a bottleneck event. The few resistant cells that survive in this new environment obtain differential reproductive success that enables them to pass down the newly selected resistant gene pool. The present review aims to summarize key findings of tumor evolution, epithelial‑mesenchymal transition and resistance to cetuximab therapy in head and neck squamous cell carcinoma.

CRISPR-Cas9-assisted genome editing in E. coli elevates the frequency of unintended mutations

Cas-assisted lambda Red recombineering techniques have rapidly become a mainstay of bacterial genome editing. Such techniques have been used to construct both individual mutants and massive libraries to assess the effects of genomic changes. We have found that a commonly used Cas9-assisted editing method results in unintended mutations elsewhere in the genome in 26% of edited clones. The unintended mutations are frequently found over 200 kb from the intended edit site and even over 10 kb from potential off-target sites. We attribute the high frequency of unintended mutations to error-prone polymerases expressed in response to dsDNA breaks introduced at the edit site. Most unintended mutations occur in regulatory or coding regions and thus may have phenotypic effects. Our findings highlight the risks associated with genome editing techniques involving dsDNA breaks in E. coli and likely other bacteria and emphasize the importance of sequencing the genomes of edited cells to ensure the absence of unintended mutations.

Fitness difference between two synonymous mutations of Phytophthora infestans ATP6 gene

BACKGROUND

Sequence variation produced by mutation provides the ultimate source of natural selection for species adaptation. Unlike nonsynonymous mutation, synonymous mutations are generally considered to be selectively neutral but accumulating evidence suggests they also contribute to species adaptation by regulating the flow of genetic information and the development of functional traits. In this study, we analysed sequence characteristics of ATP6, a housekeeping gene from 139 Phytophthora infestans isolates, and compared the fitness components including metabolic rate, temperature sensitivity, aggressiveness, and fungicide tolerance among synonymous mutations.

RESULTS

We found that the housekeeping gene exhibited low genetic variation and was represented by two major synonymous mutants at similar frequency (0.496 and 0.468, respectively). The two synonymous mutants were generated by a single nucleotide substitution but differed significantly in fitness as well as temperature-mediated spatial distribution and expression. The synonymous mutant ending in AT was more common in cold regions and was more expressed at lower experimental temperature than the synonymous mutant ending in GC and vice versa.

CONCLUSION

Our results are consistent with the argument that synonymous mutations can modulate the adaptive evolution of species including pathogens and have important implications for sustainable disease management, especially under climate change.

Fitness benefits of a synonymous substitution in an ancient EF-Tu gene depend on the genetic background

Synonymous mutations are changes to DNA sequence, which occur within translated genes but which do not affect the protein sequence. Although often referred to as silent mutations, evidence suggests that synonymous mutations can affect gene expression, mRNA stability, and even translation efficiency. A collection of both experimental and bioinformatic data has shown that synonymous mutations can impact cell phenotype, yet less is known about the molecular mechanisms and potential of beneficial or adaptive effects of such changes within evolved populations. Here, we report a beneficial synonymous mutation acquired via experimental evolution in an essential gene variant encoding the translation elongation factor protein EF-Tu. We demonstrate that this particular synonymous mutation increases EF-Tu mRNA and protein levels as well as global polysome abundance on RNA transcripts. Although presence of the synonymous mutation is clearly causative of such changes, we also demonstrate that fitness benefits are highly contingent on other potentiating mutations present within the genetic background in which the mutation arose. Our results underscore the importance of beneficial synonymous mutations, especially those that affect levels of proteins that are key for cellular processes.IMPORTANCEThis study explores the degree to which synonymous mutations in essential genes can influence adaptation in bacteria. An experimental system whereby an Escherichia coli strain harboring an engineered translation protein elongation factor-Tu (EF-Tu) was subjected to laboratory evolution. We find that a synonymous mutation acquired on the gene encoding for EF-Tu is conditionally beneficial for bacterial fitness. Our findings provide insight into the importance of the genetic background when a synonymous substitution is favored by natural selection and how such changes have the potential to impact evolution when critical cellular processes are involved.

Laboratory domestication of Lactiplantibacillus plantarum alters some phenotypic traits but causes non-novel genomic impact

AIMS

Laboratory domestication has been negligibly examined in lactic acid bacteria (LAB). Lactiplantibacillus plantarum is a highly studied and industrially relevant LAB. Here, we passaged L. plantarum JGR2 in a complex medium to study the effects of domestication on the phenotypic properties and the acquisition of mutations.

METHODS AND RESULTS

Lactiplantibacillus plantarum JGR2 was passaged in mMRS medium (deMan Rogossa Sharpe supplemented with 0.05% w/v L-cysteine) in three parallel populations for 70 days. One pure culture from each population was studied for various phenotypic properties and genomic alterations. Auto-aggregation of the evolved strains was significantly reduced, and lactic acid production and ethanol tolerance were increased. Other probiotic properties and antibiotic sensitivity were not altered. Conserved synonymous and non-synonymous mutations were observed in mobile element proteins (transposases), β-galactosidase, and phosphoketolases in all three isolates. The evolved strains lost all the repeat regions and some of the functions associated with them. Most of the conserved mutations were found in the genomes of other wild-type strains available in a public database, indicating the non-novel genomic impact of laboratory passaging.

CONCLUSIONS

Laboratory domestication can affect the phenotypic and genotypic traits of L. plantarum and similar studies are necessary for other important species of LAB.

Genetic determinants of SARS-CoV-2 and the clinical outcome of COVID-19 in Southern Bangladesh

BACKGROUND

The coronavirus disease 2019 (COVID-19) pandemic has had a severe impact on population health. The genetic determinants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in southern Bangladesh are not well understood.

METHODS

This study aimed to determine the genomic variation in SARS-CoV-2 genomes that have evolved over 2 years of the pandemic in southern Bangladesh and their association with disease outcomes and virulence of this virus. We investigated demographic variables, disease outcomes of COVID-19 patients and genomic features of SARS-CoV-2.

RESULTS

We observed that the disease severity was significantly higher in adults (85.3%) than in children (14.7%), because the expression of angiotensin-converting enzyme-2 (ACE-2) diminishes with ageing that causes differences in innate and adaptive immunity. The clade GK (n = 66) was remarkable between June 2021 and January 2022. Because of the mutation burden, another clade, GRA started a newly separated clustering in December 2021. The burden was significantly higher in GRA (1.5-fold) highlighted in mild symptoms of COVID-19 patients than in other clades (GH, GK, and GR). Mutations were accumulated mainly in S (22.15 mutations per segment) and ORF1ab segments. Missense (67.5%) and synonymous (18.31%) mutations were highly noticed in adult patients with mild cases rather than severe cases, especially in ORF1ab segments. Moreover, we observed many unique mutations in S protein in mild cases compared to severe, and homology modeling revealed that those might cause more folding in the protein's alpha helix and beta sheets.

CONCLUSION

Our study identifies some risk factors such as age comorbidities (diabetes, hypertension, and renal disease) that are associated with severe COVID-19, providing valuable insight regarding prioritizing vaccination for high-risk individuals and allocating health care and resources. The findings of this work outlined the knowledge and mutational basis of SARS-CoV-2 for the next treatment steps. Further studies are needed to confirm the effects of structural and functional proteins of SARS-CoV-2 in detail for monitoring the emergence of new variants in future.

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.

Molecular bases for strong phenotypic effects of single synonymous codon substitutions in the E. coli ccdB toxin gene

BACKGROUND

Single synonymous codon mutations typically have only minor or no effects on gene function. Here, we estimate the effects on cell growth of ~ 200 single synonymous codon mutations in an operonic context by mutating almost all positions of ccdB, the 101-residue long cytotoxin of the ccdAB Toxin-Antitoxin (TA) operon to most degenerate codons. Phenotypes were assayed by transforming the mutant library into CcdB sensitive and resistant E. coli strains, isolating plasmid pools, and subjecting them to deep sequencing. Since autoregulation is a hallmark of TA operons, phenotypes obtained for ccdB synonymous mutants after transformation in a RelE toxin reporter strain followed by deep sequencing provided information on the amount of CcdAB complex formed.

RESULTS

Synonymous mutations in the N-terminal region involved in translation initiation showed the strongest non-neutral phenotypic effects. We observe an interplay of numerous factors, namely, location of the codon, codon usage, t-RNA abundance, formation of anti-Shine Dalgarno sequences, predicted transcript secondary structure, and evolutionary conservation in determining phenotypic effects of ccdB synonymous mutations. Incorporation of an N-terminal, hyperactive synonymous mutation, in the background of the single synonymous codon mutant library sufficiently increased translation initiation, such that mutational effects on either folding or termination of translation became more apparent. Introduction of putative pause sites not only affects the translational rate, but might also alter the folding kinetics of the protein in vivo.

CONCLUSION

In summary, the study provides novel insights into diverse mechanisms by which synonymous mutations modulate gene function. This information is useful in optimizing heterologous gene expression in E. coli and understanding the molecular bases for alteration in gene expression that arise due to synonymous mutations.

Polymorphism in ovine ADCY8 gene and its association with residual feed intake in Hu sheep

Feed efficiency makes up a large part of sheep production, which also has a crucial impact on the economic benefits of producers. This study explores the polymorphism of Adenylyl cyclase 8 gene associated with feed efficiency and detects the expression characteristics of ADCY8 in ten tissues of Hu sheep. The polymorphism of ADCY8 was recognized by using PCR amplification and Sanger sequencing, KASPar technology was used for genotyping subsequently, and the relationship between SNP and RFI is also studied. The results indicated that an intronic mutation g.24799148 C > T (rs 423395741) was identified in ADCY8, and association analysis showed that the SNP g.24799148 C > T (rs 423395741) was significantly associated with RFI at 100-120, 100-140, 100-160, and 100-180 days (p < 0.05). The quantitative real-time PCR (qRT-PCR) result showed that ADCY8 was expressed in ten tissues, and the expression of ADCY8 gene in rumen tissue was significantly higher than in the other tested tissues. Therefore, these results indicated that the ADCY8 mutation locus may be used as a candidate molecular marker for evaluating the feed efficiency of Hu sheep.

Association between IRS-1, PPAR-γ and LEP genes polymorphisms and growth traits in rabbits

Single nucleotide polymorphisms are commonly associated with changes in quantitative traits, and have been considered useful markers for improving different traits in livestock. The current study aimed to explore the effect of three SNPs located in Insulin receptor substrate (IRS-1), Peroxisome proliferator-activated receptor γ (PPAR-γ), and Leptin (LEP) genes on the growth traits of rabbits. Individuals from three rabbit breeds were genotyped using RFLP-PCR. The IRS-1 variant (c.189T > G) was associated with post-weaning body weight, and body weight gains, However, the effect on growth rates was insignificant in Baladi Red and V-line rabbits. The PPAR-γ variant (c.207A > C) was significantly associated with 8-wk body weights in V-line rabbits, 10-wk body weights, and growth rates from 8 to 10 weeks of age in New Zealand rabbits. However, the differences between genotypes were insignificant for body weight gains and average daily gain. The LEP gene mutation (g.16079636C > G) had significant effects on body weights at 6 and 8 weeks of age in New Zealand White rabbits and 8 weeks of age in Baladi Red rabbits were associated with the presence of the C allele. Concludingly, the results stressed the importance of the IRS-1 gene in post-weaning growth and suggested the existence of breed-specific effects for PPAR-γ and LEP.

Mutation bias and adaptation in bacteria

Genetic mutation, which provides the raw material for evolutionary adaptation, is largely a stochastic force. However, there is ample evidence showing that mutations can also exhibit strong biases, with some mutation types and certain genomic positions mutating more often than others. It is becoming increasingly clear that mutational bias can play a role in determining adaptive outcomes in bacteria in both the laboratory and the clinic. As such, understanding the causes and consequences of mutation bias can help microbiologists to anticipate and predict adaptive outcomes. In this review, we provide an overview of the mechanisms and features of the bacterial genome that cause mutational biases to occur. We then describe the environmental triggers that drive these mechanisms to be more potent and outline the adaptive scenarios where mutation bias can synergize with natural selection to define evolutionary outcomes. We conclude by describing how understanding mutagenic genomic features can help microbiologists predict areas sensitive to mutational bias, and finish by outlining future work that will help us achieve more accurate evolutionary forecasts.

Whole-Genome Sequencing Reveals That Regulatory and Low Pleiotropy Variants Underlie Local Adaptation to Environmental Variability in Purple Sea Urchins

AbstractOrganisms experience environments that vary across both space and time. Such environmental heterogeneity shapes standing genetic variation and may influence species' capacity to adapt to rapid environmental change. However, we know little about the kind of genetic variation that is involved in local adaptation to environmental variability. To address this gap, we sequenced the whole genomes of 140 purple sea urchins (Strongylocentrotus purpuratus) from seven populations that vary in their degree of pH variability. Despite no evidence of global population structure, we found a suite of single-nucleotide polymorphisms (SNPs) tightly correlated with local pH variability (outlier SNPs), which were overrepresented in regions putatively involved in gene regulation (long noncoding RNA and enhancers), supporting the idea that variation in regulatory regions is important for local adaptation to variability. In addition, outliers in genes were found to be (i) enriched for biomineralization and ion homeostasis functions related to low pH response, (ii) less central to the protein-protein interaction network, and (iii) underrepresented among genes highly expressed during early development. Taken together, these results suggest that loci that underlie local adaptation to pH variability in purple sea urchins fall in regions with potentially low pleiotropic effects (based on analyses involving regulatory regions, network centrality, and expression time) involved in low pH response (based on functional enrichment).

A beneficial synonymous substitution in EF-Tu is contingent on genetic background

Synonymous mutations are changes to DNA sequence that occur within translated genes but which do not affect the protein sequence. Although often referred to as silent mutations, evidence suggests that synonymous mutations can affect gene expression, mRNA stability, and even translation efficiency. A collection of both experimental and bioinformatic data has shown that synonymous mutations can impact cell phenotype, yet less is known about the molecular mechanisms and potential of beneficial or adaptive effects of such changes within evolved populations. Here, we report a beneficial synonymous mutation acquired via experimental evolution in an essential gene variant encoding the translation Elongation Factor protein EF-Tu. We demonstrate that this particular synonymous mutation increases EF-Tu mRNA and protein levels, as well as the polysome abundance on global transcripts. Although presence of the synonymous mutation is clearly causative of such changes, we also demonstrate that fitness benefits are highly contingent on other potentiating mutations present within the genetic background in which the mutation arose. Our results underscore the importance of beneficial synonymous mutations, especially those that affect levels of proteins that are key for cellular processes.

Development of a Streptomyces-based system for facile thioholgamide library generation and analysis

Derivatizing natural products (NPs) is essential in structure-activity relationship (SAR) studies, compound optimization, and drug development. Ribosomally synthesized and post-translationally modified peptides (RiPPs) represent one of the major classes of natural products. Thioholgamide represents thioamitide - a recently emerged family of RiPPs with unique structures and great potential in anticancer drug development. Although the method for generating the RiPP library by codon substitutions in the precursor peptide gene is straightforward, the techniques to perform RiPP derivatization in Actinobacteria remain limited and time-consuming. Here, we report a facile system for producing a library of randomized thioholgamide derivatives utilizing an optimized Streptomyces host. This technique enabled us to access all possible amino acid substitutions of the thioholgamide molecule, one position at a time. Out of 152 potential derivatives, 85 were successfully detected, revealing the impact of amino acid substitutions on thioholgamide post-translational modifications (PTMs). Moreover, new PTMs were observed among thioholgamide derivatives: thiazoline heterocycles, which have not yet been reported for thioamitides, and S-methylmethionine, which is very rare in nature. The obtained library was subsequently used for thioholgamide SAR studies and stability assays.

Cancer genomes tolerate deleterious coding mutations through somatic copy number amplifications of wild-type regions

Cancers evolve under the accumulation of thousands of somatic mutations and chromosomal aberrations. While most coding mutations are deleterious, almost all protein-coding genes lack detectable signals of negative selection. This raises the question of how tumors tolerate such large amounts of deleterious mutations. Using 8,690 tumor samples from The Cancer Genome Atlas, we demonstrate that copy number amplifications frequently cover haploinsufficient genes in mutation-prone regions. This could increase tolerance towards the deleterious impact of mutations by creating safe copies of wild-type regions and, hence, protecting the genes therein. Our findings demonstrate that these potential buffering events are highly influenced by gene functions, essentiality, and mutation impact and that they occur early during tumor evolution. We show how cancer type-specific mutation landscapes drive copy number alteration patterns across cancer types. Ultimately, our work paves the way for the detection of novel cancer vulnerabilities by revealing genes that fall within amplifications likely selected during evolution to mitigate the effect of mutations.

Whole-genome sequencing analysis of two heat-evolved Escherichia coli strains

BACKGROUND

High temperatures cause a suite of problems for cells, including protein unfolding and aggregation; increased membrane fluidity; and changes in DNA supercoiling, RNA stability, transcription and translation. Consequently, enhanced thermotolerance can evolve through an unknown number of genetic mechanisms even in the simple model bacterium Escherichia coli. To date, each E. coli study exploring this question resulted in a different set of mutations. To understand the changes that can arise when an organism evolves to grow at higher temperatures, we sequenced and analyzed two previously described E. coli strains, BM28 and BM28 ΔlysU, that have been laboratory adapted to the highest E. coli growth temperature reported to date.

RESULTS

We found three large deletions in the BM28 and BM28 ΔlysU strains of 123, 15 and 8.5 kb in length and an expansion of IS10 elements. We found that BM28 and BM28 ΔlysU have considerably different genomes, suggesting that the BM28 culture that gave rise to BM28 and BM28 ΔlysU was a mixed population of genetically different cells. Consistent with published findings of high GroESL expression in BM28, we found that BM28 inexplicitly carries the groESL bearing plasmid pOF39 that was maintained simply by high-temperature selection pressure. We identified over 200 smaller insertions, deletions, single nucleotide polymorphisms and other mutations, including changes in master regulators such as the RNA polymerase and the transcriptional termination factor Rho. Importantly, this genome analysis demonstrates that the commonly cited findings that LysU plays a crucial role in thermotolerance and that GroESL hyper-expression is brought about by chromosomal mutations are based on a previous misinterpretation of the genotype of BM28.

CONCLUSIONS

This whole-genome sequencing study describes genetically distinct mechanisms of thermotolerance evolution from those found in other heat-evolved E. coli strains. Studying adaptive laboratory evolution to heat in simple model organisms is important in the context of climate change. It is important to better understand genetic mechanisms of enhancing thermotolerance in bacteria and other organisms, both in terms of optimizing laboratory evolution methods for various organisms and in terms of potential genetic engineering of organisms most at risk or most important to our societies and ecosystems.

Most synonymous allelic variants in HIV tat are not silent

The genetic code has degenerate codons that produce no change in the translated protein sequence and are generally thought to be silent. However, some synonymous variants are clearly not silent. Herein, we questioned the frequency of non-silent synonymous variants. We tested how random synonymous variants in the HIV Tat transcription factor effect transcription of an LTR-GFP reporter. Our model system has the advantage of directly measuring the function of the gene in human cells. Approximately, 67% of synonymous variants in Tat were non-silent, either having reduced activity or were full loss-of-function alleles. Eight mutant codons had higher codon usage than wild type, accompanied by reduced transcriptional activity. These were clustered on a loop in the Tat structure. We conclude that most synonymous Tat variants are not silent in human cells, and 25% are associated with changes in codon usage, likely effecting protein folding.

Effects of synonymous mutations on kinetic properties and structure of firefly luciferase: Molecular dynamics simulation, molecular docking, RNA folding, and experimental study

Although synonymous mutations have long been thought to lack striking results, a growing body of research shows these mutations have highly variable effects. In this study, the impact of synonymous mutations in the development of thermostable luciferase was investigated using a combination of experimental and theoretical approaches. Using bioinformatics analysis, the codon usage features in the Lampyridae family's luciferases were studied and four synonymous mutations of Arg in luciferase were created. An exciting result was that the analysis of kinetic parameters showed a slight increase in the thermal stability of the mutant luciferase. AutoDock Vina, %MinMax algorithm, and UNAFold Server were used to perform molecular docking, folding rate, and RNA folding, respectively. Here, it was assumed that in the region (Arg337) with a moderate propensity for coil, synonymous mutation altered the rate of translation, which in turn may lead to a slight change in the structure of the enzyme. According to the molecular dynamics simulation data, local minor global flexibility is observed in the context of the protein conformation. A plausible explanation is that this flexibility may strengthen hydrophobic interactions due to its sensitivity to a molecular collision. Accordingly, thermostability originated mainly from hydrophobic interaction.

Molecular evidence of parallel evolution in a cyanophage

Antagonistic interactions between bacteriophage (phage) and its bacterial host drives the continual selection for resistance and counter-defence. To date, much remains unknown about the genomic evolution that occurs as part of the underlying mechanisms. Such is the case for the marine cyanobacteria Synechococcus and viruses (cyanophages) that infect them. Here, we monitored host and phage abundances, alongside genomic changes to the phage populations, in a 500-day (~55 bacterial generations) infection experiment between Synechococcus sp. WH7803 and the T4-type cyanophage S-PM2d, run parallel in three replicate chemostats (plus one control chemostat). Flow cytometric count of total abundances revealed relatively similar host-phage population dynamics across the chemostats, starting with a cycle of host population collapse and recovery that led to phases of host-phage coexistence. Whole-genome analysis of the S-PM2d populations detected an assemblage of strongly selected and repeatable genomic changes, and therefore parallel evolution in the phage populations, early in the experiment (sampled on day 39). These consisted mostly of non-synonymous single-nucleotide-polymorphisms and a few instances of indel, altogether affecting 18 open-reading-frames, the majority of which were predicted to encode virion structures including those involved in phage adsorption onto host (i.e., baseplate wedge, short tail fibre, adhesin component). Mutations that emerged later (sampled on day 500), on the other hand, were found at a larger range of frequencies, with many lacking repeatability across the chemostats. This is indicative of some degree of between-population divergence in the phage evolutionary trajectory over time. A few of the early and late mutations were detected within putative auxiliary metabolic genes, but these generally occurred in only one or two of the chemostats. Less repeatable mutations may have higher fitness costs, thus drawing our attention onto the role of trade-offs in modulating the trajectory of a host-phage coevolution.

Parallel Evolution in Predatory Bdellovibrio sp. NC01 during Long-Term Coculture with a Single Prey Strain

Experimental evolution provides a powerful tool for examining how Bdellovibrio evolves in response to unique selective pressures associated with its predatory lifestyle. We tested how Bdellovibrio sp. NC01 adapts to long-term coculture with Pseudomonas sp. NC02, which is less susceptible to predation compared to other Gram-negative bacteria. Analyzing six replicate Bdellovibrio populations across six time points spanning 40 passages and 2,880 h of coculture, we detected 30 to 40 new mutations in each population that exceeded a frequency of 5%. Nonsynonymous substitutions were the most abundant type of new mutation, followed by small indels and synonymous substitutions. After completing the final passage, we detected 20 high-frequency (>75%) mutations across all six evolved Bdellovibrio populations. Eighteen of these alter protein sequences, and most increased in frequency rapidly. Four genes acquired a high-frequency mutation in two or more evolved Bdellovibrio populations, reflecting parallel evolution and positive selection. The genes encode a sodium/phosphate cotransporter family protein (Bd2221), a metallophosphoesterase (Bd0054), a TonB family protein (Bd0396), and a hypothetical protein (Bd1601). Tested prey range and predation efficiency phenotypes did not differ significantly between evolved Bdellovibrio populations and the ancestor; however, all six evolved Bdellovibrio populations demonstrated enhanced starvation survival compared to the ancestor. These results suggest that, instead of evolving improved killing of Pseudomonas sp. NC02, Bdellovibrio evolved to better withstand nutrient limitation in the presence of this prey strain. The mutations identified here point to genes and functions that may be important for Bdellovibrio adaptation to the different selective pressures of long-term coculture with Pseudomonas. IMPORTANCE Bdellovibrio attack and kill Gram-negative bacteria, including drug-resistant pathogens of animals and plants. This lifestyle is unusual among bacteria, and it imposes unique selective pressures on Bdellovibrio. Determining how Bdellovibrio evolve in response to these pressures is valuable for understanding the mechanisms that govern predation. We applied experimental evolution to test how Bdellovibrio sp. NC01 evolved in response to long-term coculture with a single Pseudomonas strain, which NC01 can kill, but with low efficiency. Our experimental design imposed different selective pressures on the predatory bacteria and tracked the evolutionary trajectories of replicate Bdellovibrio populations. Using genome sequencing, we identified Bdellovibrio genes that acquired high-frequency mutations in two or more populations. Using phenotype assays, we determined that evolved Bdellovibrio populations did not improve their ability to kill Pseudomonas, but rather are better able to survive starvation. Overall, our results point to functions that may be important for Bdellovibrio adaptation.

GWAS and functional studies suggest a role for altered DNA repair in the evolution of drug resistance in Mycobacterium tuberculosis

The emergence of drug resistance in Mycobacterium tuberculosis (Mtb) is alarming and demands in-depth knowledge for timely diagnosis. We performed genome-wide association analysis using 2237 clinical strains of Mtb to identify novel genetic factors that evoke drug resistance. In addition to the known direct targets, we identified for the first time, a strong association between mutations in DNA repair genes and the multidrug-resistant phenotype. To evaluate the impact of variants identified in the clinical samples in the evolution of drug resistance, we utilized knockouts and complemented strains in Mycobacterium smegmatis and Mtb. Results show that variant mutations compromised the functions of MutY and UvrB. MutY variant showed enhanced survival compared with wild-type (Rv) when the Mtb strains were subjected to multiple rounds of ex vivo antibiotic stress. In an in vivo guinea pig infection model, the MutY variant outcompeted the wild-type strain. We show that novel variant mutations in the DNA repair genes collectively compromise their functions and contribute to better survival under antibiotic/host stress conditions.

How to Recruit a Promiscuous Enzyme to Serve a New Function

Promiscuous enzymes can be recruited to serve new functions when a genetic or environmental change makes catalysis of a novel reaction important for fitness or even survival. Subsequently, gene duplication and divergence can lead to evolution of an efficient and specialized new enzyme. Every organism likely has thousands of promiscuous enzyme activities that provide a vast reservoir of catalytic potential. However, much of this potential may not be accessible. We compiled kinetic parameters for promiscuous reactions catalyzed by 108 enzymes. The median value of k_cat/K_M is a very modest 31 M^-1 s^-1. Based upon the fluxes through metabolic pathways in E. coli, we estimate that many, if not most, promiscuous activities are too inefficient to impact fitness. However, mutations can elevate the level of an insufficient promiscuous activity by increasing enzyme expression, improving k_cat/K_M, or altering concentrations of the promiscuous and native substrates and allosteric regulators. Particularly in large bacterial populations, stochastic mutations may provide a viable pathway for recruitment of even inefficient promiscuous activities.

Long-Term Adaptation to Galactose as a Sole Carbon Source Selects for Mutations Outside the Canonical GAL Pathway

Galactose is a secondary fermentable sugar that requires specific regulatory and structural genes for its assimilation, which are under catabolite repression by glucose. When glucose is absent, the catabolic repression is attenuated, and the structural GAL genes are fully activated. In Saccharomyces cerevisiae, the GAL pathway is under selection in environments where galactose is present. However, it is unclear the adaptive strategies in response to long-term propagation in galactose as a sole carbon source in laboratory evolution experiments. Here, we performed a 4,000-generation evolution experiment using 48 diploid Saccharomyces cerevisiae populations to study adaptation in galactose. We show that fitness gains were greater in the galactose-evolved population than in identically evolved populations with glucose as a sole carbon source. Whole-genome sequencing of 96 evolved clones revealed recurrent de novo single nucleotide mutations in candidate targets of selection, copy number variations, and ploidy changes. We find that most mutations that improve fitness in galactose lie outside of the canonical GAL pathway. Reconstruction of specific evolved alleles in candidate target of selection, SEC23 and IRA1, showed a significant increase in fitness in galactose compared to glucose. In addition, most of our evolved populations (28/46; 61%) fixed aneuploidies on Chromosome VIII, suggesting a parallel adaptive amplification. Finally, we show greater loss of extrachromosomal elements in our glucose-evolved lineages compared with previous glucose evolution. Broadly, these data further our understanding of the evolutionary pressures that drive adaptation to less-preferred carbon sources.

Comparative Genomics of Streptococcus oralis Identifies Large Scale Homologous Recombination and a Genetic Variant Associated with Infection

The viridans group streptococci (VGS) are a large consortium of commensal streptococci that colonize the human body. Many species within this group are opportunistic pathogens causing bacteremia and infective endocarditis (IE), yet little is known about why some strains cause invasive disease. Identification of virulence determinants is complicated by the difficulty of distinguishing between the closely related species of this group. Here, we analyzed genomic data from VGS that were isolated from blood cultures in patients with invasive infections and oral swabs of healthy volunteers and then determined the best-performing methods for species identification. Using whole-genome sequence data, we characterized the population structure of a diverse sample of Streptococcus oralis isolates and found evidence of frequent recombination. We used multiple genome-wide association study tools to identify candidate determinants of invasiveness. These tools gave consistent results, leading to the discovery of a single synonymous single nucleotide polymorphism (SNP) that was significantly associated with invasiveness. This SNP was within a previously undescribed gene that was conserved across the majority of VGS species. Using the growth in the presence of human serum and a simulated infective endocarditis vegetation model, we were unable to identify a phenotype for the enriched allele in laboratory assays, suggesting a phenotype may be specific to natural infection. These data highlighted the power of analyzing natural populations for gaining insight into pathogenicity, particularly for organisms with complex population structures like the VGS. IMPORTANCE The viridians group streptococci (VGS) are a large collection of closely related commensal streptococci, with many being opportunistic pathogens causing invasive diseases, such as bacteremia and infective endocarditis. Little is known about virulence determinants in these species, and there is a distinct lack of genomic information available for the VGS. In this study, we collected VGS isolates from invasive infections and healthy volunteers and performed whole-genome sequencing for a suite of downstream analyses. We focused on a diverse sample of Streptococcus oralis genomes and identified high rates of recombination in the population as well as a single genome variant highly enriched in invasive isolates. The variant lies within a previously uncharacterized gene, nrdM, which shared homology with the anaerobic ribonucleoside triphosphate reductase, nrdD, and was highly conserved among VGS. This work increased our knowledge of VGS genomics and indicated that differences in virulence potential among S. oralis isolates were, at least in part, genetically determined.

Exome variants associated with asthma and allergy

The mutational spectrum of asthma and allergy associated genes is not known although recent biobank based exome sequencing studies included these traits. We therefore conducted a secondary analysis of exome data from 281,104 UK Biobank samples for association of mostly rare variants with asthma, allergic rhinitis and atopic dermatitis. Variants of interest (VOI) were tabulated, shared genes annotated and compared to earlier genome-wide SNP association studies (GWAS), whole genome sequencing, exome and bisulfit sequencing studies. 354 VOI were significantly associated with asthma, allergic rhinitis and atopic dermatitis. They cluster mainly in two large regions on chromosome 6 and 17. After exclusion of the variants associated with atopic dermatitis and redundant variants, 321 unique VOI remain in 122 unique genes. 30 genes are shared among the 87 genes with increased and the 65 genes with decreased risk for allergic disease. 85% of genes identified earlier by common GWAS SNPs are not replicated here. Most identified genes are located in interferon ɣ and IL33 signaling pathway. These genes include already known but also new pharmacological targets, including the IL33 receptor ST2/IL1RL1, as well as TLR1, ALOX15, GSDMA, BTNL2, IL13 and IKZF3. Future pharmacological studies will need to included these VOI for stratification of the study population paving the way to individualized treatment.

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 High Mutational Sensitivity of ccdA Antitoxin Is Linked to Codon Optimality

Deep mutational scanning studies suggest that synonymous mutations are typically silent and that most exposed, nonactive-site residues are tolerant to mutations. Here, we show that the ccdA antitoxin component of the Escherichia coli ccdAB toxin-antitoxin system is unusually sensitive to mutations when studied in the operonic context. A large fraction (∼80%) of single-codon mutations, including many synonymous mutations in the ccdA gene shows inactive phenotype, but they retain native-like binding affinity towards cognate toxin, CcdB. Therefore, the observed phenotypic effects are largely not due to alterations in protein structure/stability, consistent with a large region of CcdA being intrinsically disordered. E. coli codon preference and strength of ribosome-binding associated with translation of downstream ccdB gene are found to be major contributors of the observed ccdA mutant phenotypes. In select cases, proteomics studies reveal altered ratios of CcdA:CcdB protein levels in vivo, suggesting that the ccdA mutations likely alter relative translation efficiencies of the two genes in the operon. We extend these results by studying single-site synonymous mutations that lead to loss of function phenotypes in the relBE operon upon introduction of rarer codons. Thus, in their operonic context, genes are likely to be more sensitive to both synonymous and nonsynonymous point mutations than inferred previously.

Genetic variants in CYP11B1 influence the susceptibility to coronary heart disease

BACKGROUND

Genetic factors are important risk factors to develop coronary heart disease (CHD). In this study, we mainly explored whether CYP11B1 mutations influence CHD risk among Chinese Han population.

METHODS

Six variants were genotyped using Agena MassARRAY system from 509 CHD patients and 509 healthy controls. The correlations between CYP11B1 mutations and CHD risk were assessed using odds ratio (OR) and 95% confidence interval (95% CI) by logistic regression. The haplotype analysis and were ultifactor dimensionality reduction (MDR) were conducted.

RESULTS

In the overall analysis, CYP11B1 polymorphisms were not correlated with CHD susceptibility. In the stratified analysis, we found that rs5283, rs6410, and rs4534 are significantly associated with susceptibility to CHD dependent on age and gender (p < 0.05). Moreover, we also observed that rs5283 and rs4534 could affect diabetes/hypertension risk among CHD patients (p < 0.05). In addition, the C_rs4736312A_rs5017238C_rs5301G_rs5283T_rs6410C_rs4534 haplotype of CYP11B1 reduce the susceptibility to CHD (p < 0.05).

CONCLUSIONS

We found that rs4534, rs6410 and rs5283 in CYP11B1 gene influence the susceptibility to CHD, which depend on age and gender.

Synonymous mutations in representative yeast genes are mostly strongly non-neutral

Synonymous mutations in protein-coding genes do not alter protein sequences and are thus generally presumed to be neutral or nearly neutral^1-5. Here, to experimentally verify this presumption, we constructed 8,341 yeast mutants each carrying a synonymous, nonsynonymous or nonsense mutation in one of 21 endogenous genes with diverse functions and expression levels and measured their fitness relative to the wild type in a rich medium. Three-quarters of synonymous mutations resulted in a significant reduction in fitness, and the distribution of fitness effects was overall similar-albeit nonidentical-between synonymous and nonsynonymous mutations. Both synonymous and nonsynonymous mutations frequently disturbed the level of mRNA expression of the mutated gene, and the extent of the disturbance partially predicted the fitness effect. Investigations in additional environments revealed greater across-environment fitness variations for nonsynonymous mutants than for synonymous mutants despite their similar fitness distributions in each environment, suggesting that a smaller proportion of nonsynonymous mutants than synonymous mutants are always non-deleterious in a changing environment to permit fixation, potentially explaining the common observation of substantially lower nonsynonymous than synonymous substitution rates. The strong non-neutrality of most synonymous mutations, if it holds true for other genes and in other organisms, would require re-examination of numerous biological conclusions about mutation, selection, effective population size, divergence time and disease mechanisms that rely on the assumption that synoymous mutations are neutral.

Genomics of diversification of Pseudomonas aeruginosa in cystic fibrosis lung-like conditions

Pseudomonas aeruginosa is among the most problematic opportunistic pathogens for adults with cystic fibrosis (CF), causing repeated and resilient infections in the lung and surrounding airways. Evidence suggests that long-term infections are associated with diversification into specialized types but the underlying cause of that diversification and the effect it has on the persistence of infections remains poorly understood. Here, we use evolve-and-resequence experiments to investigate the genetic changes accompanying rapid, de novo phenotypic diversification in lab environments designed to mimic two aspects of human lung ecology: spatial structure and complex nutritional content. After ∼220 generations of evolution, we find extensive genetic variation present in all environments, including those that most closely resemble the CF lung. We use the abundance and frequency of nonsynonymous and synonymous mutations to estimate the ratio of mutations that are selectively neutral (hitchhikers) to those that are under positive selection (drivers). A significantly lower proportion of driver mutations in spatially structured populations suggests that reduced dispersal generates subpopulations with reduced effective population size, decreasing the supply of beneficial mutations and causing more divergent evolutionary trajectories. In addition, we find mutations in a handful of genes typically associated with chronic infection in the CF lung, including one gene associated with antibiotic resistance. This demonstrates that many of the genetic changes considered to be hallmarks of CF lung adaptation can arise as a result of adaptation to a novel environment and do not necessarily require antimicrobial treatment, immune system suppression, or competition from other microbial species to occur.

Evolutionary Analyses of Emerging GII.2[P16] and GII.4 Sydney [P16] Noroviruses

GII.2[P16] and GII.4 Sydney [P16] are currently the two predominant norovirus genotypes. This study sought to clarify their evolutionary patterns by analyzing the major capsid VP1 and RNA-dependent RNA polymerase (RdRp) genes. Sequence diversities were analyzed at both nucleotide and amino acid levels. Selective pressures were evaluated with the Hyphy package in different models. Phylogenetic trees were constructed by the maximum likelihood method from full VP1 sequences, and evolutionary rates were estimated by the Bayesian Markov Chain Monte Carlo approach. The results showed that (1) several groups of tightly linked mutations between the RdRp and VP1 genes were detected in the GII.2[P16] and GII.4[P16] noroviruses, and most of these mutations were synonymous, which may lead to a better viral fitness to the host; (2) although the pattern of having new GII.4 variants every 2–4 years has been broken, both the pre- and the post-2015 Sydney VP1 had comparable evolutionary rates to previously epidemic GII.4 variants, and half of the major antigenic sites on GII.4 Sydney had residue substitutions and several caused obvious changes in the carbohydrate-binding surface that may potentially alter the property of the virus; and (3) GII.4 Sydney variants during 2018–21 showed geographical specificity in East Asia, South Asia, and North America; the antigenic sites of GII.2 are strictly conserved, but the GII.2 VP1 chronologically evolved into nine different sublineages over time, with sublineage IX being the most prevalent one since 2018. This study suggested that both VP1 and RdRp of the GII.2[P16] and GII.4 Sydney [P16] noroviruses exhibited different evolutionary directions. GII.4[P16] is likely to generate potential novel epidemic variants by accumulating mutations in the P2 domain, similar to previously epidemic GII.4 variants, while GII.2[P16] has conserved predicted antigenicity and may evolve by changing the properties of nonstructural proteins, such as polymerase replicational fidelity and efficiency. This study expands the understanding of the evolutionary dynamics of GII.2[P16] and GII.4[P16] noroviruses and may predict the emergence of new variants.

Dominant clade-featured SARS-CoV-2 co-occurring mutations reveal plausible epistasis: An in silico based hypothetical model

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved into eight fundamental clades with four of these clades (G, GH, GR, and GV) globally prevalent in 2020. To explain plausible epistatic effects of the signature co-occurring mutations of these circulating clades on viral replication and transmission fitness, we proposed a hypothetical model using in silico approach. Molecular docking and dynamics analyses showed the higher infectiousness of a spike mutant through more favorable binding of G614 with the elastase-2. RdRp mutation p.P323L significantly increased genome-wide mutations (p < 0.0001), allowing for more flexible RdRp (mutated)-NSP8 interaction that may accelerate replication. Superior RNA stability and structural variation at NSP3:C241T might impact protein, RNA interactions, or both. Another silent 5'-UTR:C241T mutation might affect translational efficiency and viral packaging. These four G-clade-featured co-occurring mutations might increase viral replication. Sentinel GH-clade ORF3a:p.Q57H variants constricted the ion-channel through intertransmembrane-domain interaction of cysteine(C81)-histidine(H57). The GR-clade N:p.RG203-204KR would stabilize RNA interaction by a more flexible and hypo-phosphorylated SR-rich region. GV-clade viruses seemingly gained the evolutionary advantage of the confounding factors; nevertheless, N:p.A220V might modulate RNA binding with no phenotypic effect. Our hypothetical model needs further retrospective and prospective studies to understand detailed molecular events and their relationship to the fitness of SARS-CoV-2.

QTL mapping of a Brazilian bioethanol strain links the cell wall protein-encoding gene GAS1 to low pH tolerance in S. cerevisiae

BACKGROUND

Saccharomyces cerevisiae is largely applied in many biotechnological processes, from traditional food and beverage industries to modern biofuel and biochemicals factories. During the fermentation process, yeast cells are usually challenged in different harsh conditions, which often impact productivity. Regarding bioethanol production, cell exposure to acidic environments is related to productivity loss on both first- and second-generation ethanol. In this scenario, indigenous strains traditionally used in fermentation stand out as a source of complex genetic architecture, mainly due to their highly robust background-including low pH tolerance.

RESULTS

In this work, we pioneer the use of QTL mapping to uncover the genetic basis that confers to the industrial strain Pedra-2 (PE-2) acidic tolerance during growth at low pH. First, we developed a fluorescence-based high-throughput approach to collect a large number of haploid cells using flow cytometry. Then, we were able to apply a bulk segregant analysis to solve the genetic basis of low pH resistance in PE-2, which uncovered a region in chromosome X as the major QTL associated with the evaluated phenotype. A reciprocal hemizygosity analysis revealed the allele GAS1, encoding a β-1,3-glucanosyltransferase, as the casual variant in this region. The GAS1 sequence alignment of distinct S. cerevisiae strains pointed out a non-synonymous mutation (A631G) prevalence in wild-type isolates, which is absent in laboratory strains. We further showcase that GAS1 allele swap between PE-2 and a low pH-susceptible strain can improve cell viability on the latter of up to 12% after a sulfuric acid wash process.

CONCLUSION

This work revealed GAS1 as one of the main causative genes associated with tolerance to growth at low pH in PE-2. We also showcase how GAS1PE-2 can improve acid resistance of a susceptible strain, suggesting that these findings can be a powerful foundation for the development of more robust and acid-tolerant strains. Our results collectively show the importance of tailored industrial isolated strains in discovering the genetic architecture of relevant traits and its implications over productivity.

Genome-wide association study of gastric cancer- and duodenal ulcer-derived Helicobacter pylori strains reveals discriminatory genetic variations and novel oncoprotein candidates

Genome-wide association studies (GWASs) can reveal genetic variations associated with a phenotype in the absence of any hypothesis of candidate genes. The problem of false-positive sites linked with the responsible site might be bypassed in bacteria with a high homologous recombination rate, such as Helicobacter pylori, which causes gastric cancer. We conducted a small-sample GWAS (125 gastric cancer cases and 115 controls) followed by prediction of gastric cancer and control (duodenal ulcer) H. pylori strains. We identified 11 single nucleotide polymorphisms (eight amino acid changes) and three DNA motifs that, combined, allowed effective disease discrimination. They were often informative of the underlying molecular mechanisms, such as electric charge alteration at the ligand-binding pocket, alteration in subunit interaction, and mode-switching of DNA methylation. We also identified three novel virulence factors/oncoprotein candidates. These results provide both defined targets for further informatic and experimental analyses to gain insights into gastric cancer pathogenesis and a basis for identifying a set of biomarkers for distinguishing these H. pylori-related diseases.

A domain damage index to prioritizing the pathogenicity of missense variants

Prioritizing causal variants is one major challenge for the clinical application of sequencing data. Prompted by the observation that 74.3% of missense pathogenic variants locate in protein domains, we developed an approach named domain damage index (DDI). DDI identifies protein domains depleted of rare missense variations in the general population, which can be further used as a metric to prioritize variants. DDI is significantly correlated with phylogenetic conservation, variant-level metrics, and reported pathogenicity. DDI achieved great performance for distinguishing pathogenic variants from benign ones in three benchmark datasets. The combination of DDI with the other two best approaches improved the performance of each individual method considerably, suggesting DDI provides a powerful and complementary way of variant prioritization.

Molecular signatures of Janthinobacterium lividum from Trinidad support high potential for crude oil metabolism

Background

Janthinobacterium lividum is considered to be a psychrotrophic bacterial species. For the first time in the literature, J. lividum strains were isolated from Trinidad presenting with atypical features - hydrocarbonoclastic and able to survive in a tropical environment.

Methods

Identification of the Trinidad strains was carried out through 16S rRNA phylogenetic analysis. Gene-specific primers were designed to target the VioA which encodes violacein pigment and the EstA/B gene which encodes secreted extracellular lipase. Bioinformatics analyses were carried out on the nucleotide and amino acid sequences of VioA and EstA/B genes of the Trinidad Janthinobacterium strains to assess functionality and phylogenetic relatedness to other Janthinobacterium sequences specifically and more broadly, to other members of the Oxalobacteraceae family of betaproteobacteria.

Results

16S rRNA confirmed the identity of the Trinidad strains as J. lividum and resolved three of the Trinidad strains at the intra-specific level. Typical motility patterns of this species were recorded. VioAp sequences were highly conserved, however, synonymous substitutions located outside of the critical sites for enzyme function were detected for the Trinidad strains. Comparisons with PDB 6g2p model from aa231 to aa406 further indicated no functional disruption of the VioA gene of the Trinidad strains. Phylogeny of the VioA protein sequences inferred placement of all J. lividum taxa into a highly supported species-specific clade (bs = 98%). EstA/Bp sequences were highly conserved, however, synonymous substitutions were detected that were unique to the Trinidad strains. Phylogenetic inference positioned the Trinidad consensus VioA and EstA protein sequences in a clearly distinct branch.

Conclusions

The findings showed that the primary sequence of VioAp and EstA/Bp were unique to the Trinidad strains and these molecular signatures were reflected in phylogenetic inference. Our results supported chemotaxis, possible elective inactivation of VioA gene expression and secreted lipase activity as survival mechanisms of the Trinidad strains in petrogenic conditions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02346-4.

A mutational hotspot that determines highly repeatable evolution can be built and broken by silent genetic changes

Mutational hotspots can determine evolutionary outcomes and make evolution repeatable. Hotspots are products of multiple evolutionary forces including mutation rate heterogeneity, but this variable is often hard to identify. In this work, we reveal that a near-deterministic genetic hotspot can be built and broken by a handful of silent mutations. We observe this when studying homologous immotile variants of the bacteria Pseudomonas fluorescens, AR2 and Pf0-2x. AR2 resurrects motility through highly repeatable de novo mutation of the same nucleotide in >95% lines in minimal media (ntrB A289C). Pf0-2x, however, evolves via a number of mutations meaning the two strains diverge significantly during adaptation. We determine that this evolutionary disparity is owed to just 6 synonymous variations within the ntrB locus, which we demonstrate by swapping the sites and observing that we are able to both break (>95% to 0%) and build (0% to 80%) a deterministic mutational hotspot. Our work reveals a key role for silent genetic variation in determining adaptive outcomes.

Amplicon remodeling and genomic mutations drive population dynamics after segmental amplification

New enzymes often evolve by duplication and divergence of genes encoding enzymes with promiscuous activities that have become important in the face of environmental opportunities or challenges. Amplifications that increase the copy number of the gene under selection commonly amplify many surrounding genes. Extra copies of these coamplified genes must be removed, either during or after evolution of a new enzyme. Here we report that amplicon remodeling can begin even before mutations occur in the gene under selection. Amplicon remodeling and mutations elsewhere in the genome that indirectly increase fitness result in complex population dynamics, leading to emergence of clones that have improved fitness by different mechanisms. In this work, one of the two most successful clones had undergone two episodes of amplicon remodeling, leaving only four coamplified genes surrounding the gene under selection. Amplicon remodeling in the other clone resulted in removal of 111 genes from the genome, an acceptable solution under these selection conditions, but one that would certainly impair fitness under other environmental conditions.

Superspreading in the emergence of COVID-19 variants

Superspreading and variants of concern (VOC) of the human pathogen severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are the main catalyzers of the coronavirus disease 2019 (COVID-19) pandemic. However, measuring their individual impact is challenging. By examining the largest database of SARS-CoV-2 genomes The Global Initiative on Sharing Avian Influenza Data [GISAID; n >1.2 million high-quality (HQ) sequences], we present evidence suggesting that superspreading has had a key role in the epidemiological predominance of VOC. There are clear signatures in the database compatible with large superspreading events (SSEs) coinciding chronologically with the worst epidemiological scenarios triggered by VOC. The data suggest that, without the randomness effect of the genetic drift facilitated by superspreading, new VOC of SARS-CoV-2 would have had more limited chance of success.

Relationship of Postoperative Pain and PONV after Minimally Invasive Surgery with the Serotonin Concentrations and Receptors' Gene Polymorphisms

(1) Background: there is a steady increase in the number of procedures performed via minimally invasive surgery, which have many benefits, but post-operative nausea and vomiting (PONV) and significant pain are still a common problem (2) Methods: 300 infertile women (18–40 years old) undergoing minimal invasive surgery. Interventions: laparoscopy and hysteroscopy performing, evaluation of postoperative symptoms, serotonin concentrations assessment, identify genetic polymorphisms. (3) Results: serotonin concentrations were significantly lower among women who required opioids (p = 0.006). The presence of the GG genotype in the rs6318 polymorphism of the 5HTR2C gene had a protective effect on PONV (OR = 0.503; C.I. = [0.300–0.841]; p = 0.008), when the GG variant of the rs11214763 polymorphism of the 5HTR3B gene, when the risk of PONV was 1.65-fold higher (OR = 1.652; C.I. = [1.003–2.723]; p = 0.048). Pain intensity was significantly higher among women with GG genotype of the rs6296 polymorphism of the 5HTR1B gene (OR = 1.660; C.I. = [1.052–2.622]; p = 0.029).; (4) Conclusions: the evaluation of serotonin concentration predicts requirement for opioid pain relief medication. The polymorphisms of the serotonin receptors affect the intensity of postoperative complaints.

Evolutionary pressures and codon bias in low complexity regions of plasmodia

The biological meaning of low complexity regions in the proteins of Plasmodium species is a topic of discussion in evolutionary biology. There is a debate between selectionists and neutralists, who either attribute or do not attribute an effect of low-complexity regions on the fitness of these parasites, respectively. In this work, we comparatively study 22 Plasmodium species to understand whether their low complexity regions undergo a neutral or, rather, a selective and species-dependent evolution. The focus is on the connection between the codon repertoire of the genetic coding sequences and the occurrence of low complexity regions in the corresponding proteins. The first part of the work concerns the correlation between the length of plasmodial proteins and their propensity at embedding low complexity regions. Relative synonymous codon usage, entropy, and other indicators reveal that the incidence of low complexity regions and their codon bias is species-specific and subject to selective evolutionary pressure. We also observed that protein length, a relaxed selective pressure, and a broad repertoire of codons in proteins, are strongly correlated with the occurrence of low complexity regions. Overall, it seems plausible that the codon bias of low-complexity regions contributes to functional innovation and codon bias enhancement of proteins on which Plasmodium species rest as successful evolutionary parasites.

Effects of synonymous mutations beyond codon bias: The evidence for adaptive synonymous substitutions from microbial evolution experiments

Synonymous mutations are often assumed to be neutral with respect to fitness because they do not alter the encoded amino acid and so cannot be 'seen' by natural selection. Yet a growing body of evidence suggests that synonymous mutations can have fitness effects that drive adaptive evolution through their impacts on gene expression and protein folding. Here, we review what microbial experiments have taught us about the contribution of synonymous mutations to adaptation. A survey of site-directed mutagenesis experiments reveals the distributions of fitness effects for nonsynonymous and synonymous mutations are more similar, especially for beneficial mutations, than expected if all synonymous mutations were neutral, suggesting they should drive adaptive evolution more often than is typically observed. A review of experimental evolution studies where synonymous mutations have contributed to adaptation shows they can impact fitness through a range of mechanisms including the creation of illicit RNA polymerase binding sites impacting transcription and changes to mRNA folding stability that modulate translation. We suggest that clonal interference in evolving microbial populations may be the reason synonymous mutations play a smaller role in adaptive evolution than expected based on their observed fitness effects. We finish by discussing the impacts of falsely assuming synonymous mutations are neutral and discuss directions for future work exploring the role of synonymous mutations in adaptive evolution.

Setting the stage for evolution of a new enzyme

The evolution of novel enzymes has fueled the diversification of life on earth for billions of years. Insights into events that set the stage for the evolution of a new enzyme can be obtained from ancestral reconstruction and laboratory evolution. Ancestral reconstruction can reveal the emergence of a promiscuous activity in a pre-existing protein and the impact of subsequent mutations that enhance a new activity. Laboratory evolution provides a more holistic view by revealing mutations elsewhere in the genome that indirectly enhance the level of a newly important enzymatic activity. This review will highlight recent studies that probe the early stages of the evolution of a new enzyme from these complementary points of view.

From genotype to phenotype: adaptations of Pseudomonas aeruginosa to the cystic fibrosis environment

Pseudomonas aeruginosa is one of the main microbial species colonizing the lungs of cystic fibrosis patients and is responsible for the decline in respiratory function. Despite the hostile pulmonary environment, P. aeruginosa is able to establish chronic infections thanks to its strong adaptive capacity. Various longitudinal studies have attempted to compare the strains of early infection with the adapted strains of chronic infection. Thanks to new '-omics' techniques, convergent genetic mutations, as well as transcriptomic and proteomic dysregulations have been identified. As a consequence of this evolution, the adapted strains of P. aeruginosa have particular phenotypes that promote persistent infection.

Characterizing genomic variants and mutations in SARS-CoV-2 proteins from Indian isolates

SARS-CoV-2 is mutating and creating divergent variants by altering the composition of essential constituent proteins. Pharmacologically, it is crucial to understand the diverse mechanism of mutations for stable vaccine or anti-viral drug design. Our current study concentrates on all the constituent proteins of 469 SARS-CoV-2 genome samples, derived from Indian patients. However, the study may easily be extended to the samples across the globe.

We perform clustering analysis towards identifying unique variants in each of the SARS-CoV-2 proteins. A total of 536 mutated positions within the coding regions of SARS-CoV-2 proteins are detected among the identified variants from Indian isolates. We quantify mutations by focusing on the unique variants of each SARS-CoV-2 protein. We report the average number of mutation per variant, percentage of mutated positions, synonymous and non-synonymous mutations, mutations occurring in three codon positions and so on. Our study reveals the most susceptible six (06) proteins, which are ORF1ab, Spike (S), Nucleocapsid (N), ORF3a, ORF7a, and ORF8. Several non-synonymous substitutions are observed to be unique in different SARS-CoV-2 proteins. A total of 57 possible deleterious amino acid substitutions are predicted, which may impact on the protein functions. Several mutations show a large decrease in protein stability and are observed in putative functional domains of the proteins that might have some role in disease pathogenesis. We observe a good number of physicochemical property change during above deleterious substitutions.

Polymorphisms in ovine ME1 and CA1 genes and their association with feed efficiency in Hu sheep

Improving feed efficiency would increase profitability for producers. The objective of this study was to detect the expression levels of ME1 and CA1 and the polymorphisms of ME1 and CA1 associated with the feed conversion ratio (FCR) and residual feed intake (RFI) of Hu sheep by using qRT-PCR, pooled DNA sequencing and KASPar assay. The qRT-PCR results indicated that the expression levels of ME1 and CA1 were significantly higher in the liver tissues of low-RFI sheep than in those of the high-RFI sheep (p <.01). Association analysis demonstrated that the polymorphism ME1 g.453 C > T was significantly associated with FCR and RFI (p <.05). The polymorphism CA1 g.199 C > T had a significantly associated with FCR (p <.05) and no association with RFI (p >.05). Significant differences were observed between the combined genotypes and FCR and RFI at different measurement periods (p <.05). Thus, we propose the use of these two polymorphisms as new candidate molecular markers for improving feed efficiency in sheep populations.

ΦX174 Attenuation by Whole-Genome Codon Deoptimization

Natural selection acting on synonymous mutations in protein-coding genes influences genome composition and evolution. In viruses, introducing synonymous mutations in genes encoding structural proteins can drastically reduce viral growth, providing a means to generate potent, live-attenuated vaccine candidates. However, an improved understanding of what compositional features are under selection and how combinations of synonymous mutations affect viral growth is needed to predictably attenuate viruses and make them resistant to reversion. We systematically recoded all nonoverlapping genes of the bacteriophage ΦX174 with codons rarely used in its Escherichia coli host. The fitness of recombinant viruses decreases as additional deoptimizing mutations are made to the genome, although not always linearly, and not consistently across genes. Combining deoptimizing mutations may reduce viral fitness more or less than expected from the effect size of the constituent mutations and we point out difficulties in untangling correlated compositional features. We test our model by optimizing the same genes and find that the relationship between codon usage and fitness does not hold for optimization, suggesting that wild-type ΦX174 is at a fitness optimum. This work highlights the need to better understand how selection acts on patterns of synonymous codon usage across the genome and provides a convenient system to investigate the genetic determinants of virulence.

What Has a Century of Quantitative Genetics Taught Us About Nature's Genetic Tool Kit?

The complexity of heredity has been appreciated for decades: Many traits are controlled not by a single genetic locus but instead by polymorphisms throughout the genome. The importance of complex traits in biology and medicine has motivated diverse approaches to understanding their detailed genetic bases. Here, we focus on recent systematic studies, many in budding yeast, which have revealed that large numbers of all kinds of molecular variation, from noncoding to synonymous variants, can make significant contributions to phenotype. Variants can affect different traits in opposing directions, and their contributions can be modified by both the environment and the epigenetic state of the cell. The integration of prospective (synthesizing and analyzing variants) and retrospective (examining standing variation) approaches promises to reveal how natural selection shapes quantitative traits. Only by comprehensively understanding nature's genetic tool kit can we predict how phenotypes arise from the complex ensembles of genetic variants in living organisms.

Synthetic cross-phyla gene replacement and evolutionary assimilation of major enzymes

The ability of DNA to produce a functional protein even after transfer to a foreign host is of fundamental importance in both evolutionary biology and biotechnology, enabling horizontal gene transfer in the wild and heterologous expression in the lab. However, the influence of genetic particulars on DNA functionality in a new host is poorly understood, as are the evolutionary mechanisms of assimilation and refinement. Here, we describe an automation-enabled large-scale experiment wherein Escherichia coli strains were evolved in parallel after replacement of the genes pgi or tpiA with orthologous DNA from donor species spanning all domains of life, from humans to hyperthermophilic archaea. Via analysis of hundreds of clones evolved for 50,000+ cumulative generations across dozens of independent lineages, we show that orthogene-upregulating mutations can completely mitigate fitness defects that result from initial non-functionality, with coding sequence changes unnecessary. Gene target, donor species and genomic location of the swap all influenced outcomes-both the nature of adaptive mutations (often synonymous) and the frequency with which strains successfully evolved to assimilate the foreign DNA. Additionally, time series DNA sequencing and replay evolution experiments revealed transient copy number expansions, the contingency of lineage outcome on first-step mutations and the ability for strains to escape from suboptimal local fitness maxima. Overall, this study establishes the influence of various DNA and protein features on cross-species genetic interchangeability and evolutionary outcomes, with implications for both horizontal gene transfer and rational strain design.

Synonymous genome recoding: a tool to explore microbial biology and new therapeutic strategies

Synthetic genome recoding is a new means of generating designed organisms with altered phenotypes. Synonymous mutations introduced into the protein coding region tolerate modifications in DNA or mRNA without modifying the encoded proteins. Synonymous genome-wide recoding has allowed the synthetic generation of different small-genome viruses with modified phenotypes and biological properties. Recently, a decreased cost of chemically synthesizing DNA and improved methods for assembling DNA fragments (e.g. lambda red recombination and CRISPR-based editing) have enabled the construction of an Escherichia coli variant with a 4-Mb synthetic synonymously recoded genome with a reduced number of sense codons (n = 59) encoding the 20 canonical amino acids. Synonymous genome recoding is increasing our knowledge of microbial interactions with innate immune responses, identifying functional genome structures, and strategically ameliorating cis-inhibitory signaling sequences related to splicing, replication (in eukaryotes), and complex microbe functions, unraveling the relevance of codon usage for the temporal regulation of gene expression and the microbe mutant spectrum and adaptability. New biotechnological and therapeutic applications of this methodology can easily be envisaged. In this review, we discuss how synonymous genome recoding may impact our knowledge of microbial biology and the development of new and better therapeutic methodologies.