Origins of spontaneous mutations: specificity and directionality of base-substitution, frameshift, and sequence-substitution mutageneses

Global and local genomic features together modulate the spontaneous single nucleotide mutation rate

Spontaneous mutations are evolutionary engines as they generate variants for the evolutionary downstream processes that give rise to speciation and adaptation. Single nucleotide mutations (SNM) are the most abundant type of mutations among them. Here, we perform a meta-analysis to quantify the influence of selected global genomic parameters (genome size, genomic GC content, genomic repeat fraction, number of coding genes, gene count, and strand bias in prokaryotes) and local genomic features (local GC content, repeat content, CpG content and the number of SNM at CpG islands) on spontaneous SNM rates across the tree of life (prokaryotes, unicellular eukaryotes, multicellular eukaryotes) using wild-type sequence data in two different taxon classification systems. We find that the spontaneous SNM rates in our data are correlated with many genomic features in prokaryotes and unicellular eukaryotes irrespective of their sample sizes. On the other hand, only the number of coding genes was correlated with the spontaneous SNM rates in multicellular eukaryotes primarily contributed by vertebrates data. Considering local features, we notice that local GC content and CpG content significantly were correlated with the spontaneous SNM rates in the unicellular eukaryotes, while local repeat fraction is an important feature in prokaryotes and certain specific uni- and multi-cellular eukaryotes. Such predictive features of the spontaneous SNM rates often support non-linear models as the best fit compared to the linear model. We also observe that the strand asymmetry in prokaryotes plays an important role in determining the spontaneous SNM rates but the SNM spectrum does not.

The drought-induced plasticity of mineral nutrients contributes to drought tolerance discrimination in durum wheat

Drought is a major challenge for the cultivation of durum wheat, a crucial crop for global food security. Plants respond to drought by adjusting their mineral nutrient profiles to cope with water scarcity, showing the importance of nutrient plasticity for plant acclimation and adaptation to diverse environments. Therefore, it is essential to understand the genetic basis of mineral nutrient profile plasticity in durum wheat under drought stress to select drought-tolerant varieties. The research study investigated the responses of different durum wheat genotypes to severe drought stress at the seedling stage. The study employed an ionomic, molecular, biochemical and physiological approach to shed light on distinct behaviors among different genotypes. The drought tolerance of SVEMS16, SVEVO, and BULEL was related to their capacity of maintaining or increasing nutrient's accumulation, while the limited nutrient acquisition capability of CRESO and S.CAP likely resulted in their susceptibility to drought. The study highlighted the importance of macronutrients such as SO42-, NO3-, PO43-, and K+ in stress resilience and identified variant-containing genes potentially influencing nutritional variations under drought. These findings provide valuable insights for further field studies to assess the drought tolerance of durum wheat genotypes across various growth stages, ultimately ensuring food security and sustainable production in the face of changing environmental conditions.

Properties and Mechanisms of Deletions, Insertions, and Substitutions in the Evolutionary History of SARS-CoV-2

SARS-CoV-2 has accumulated many mutations since its emergence in late 2019. Nucleotide substitutions leading to amino acid replacements constitute the primary material for natural selection. Insertions, deletions, and substitutions appear to be critical for coronavirus's macro- and microevolution. Understanding the molecular mechanisms of mutations in the mutational hotspots (positions, loci with recurrent mutations, and nucleotide context) is important for disentangling roles of mutagenesis and selection. In the SARS-CoV-2 genome, deletions and insertions are frequently associated with repetitive sequences, whereas C>U substitutions are often surrounded by nucleotides resembling the APOBEC mutable motifs. We describe various approaches to mutation spectra analyses, including the context features of RNAs that are likely to be involved in the generation of recurrent mutations. We also discuss the interplay between mutations and natural selection as a complex evolutionary trend. The substantial variability and complexity of pipelines for the reconstruction of mutations and the huge number of genomic sequences are major problems for the analyses of mutations in the SARS-CoV-2 genome. As a solution, we advocate for the development of a centralized database of predicted mutations, which needs to be updated on a regular basis.

DNA polymerase ε and δ variants drive mutagenesis in polypurine tracts in human tumors

Alterations in the exonuclease domain of DNA polymerase ε cause ultramutated cancers. These cancers accumulate AGA>ATA transversions; however, their genomic features beyond the trinucleotide motifs are obscure. We analyze the extended DNA context of ultramutation using whole-exome sequencing data from 524 endometrial and 395 colorectal tumors. We find that G>T transversions in POLE-mutant tumors predominantly affect sequences containing at least six consecutive purines, with a striking preference for certain positions within polypurine tracts. Using this signature, we develop a machine-learning classifier to identify tumors with hitherto unknown POLE drivers and validate two drivers, POLE-E978G and POLE-S461L, by functional assays in yeast. Unlike other pathogenic variants, the E978G substitution affects the polymerase domain of Pol ε. We further show that tumors with POLD1 drivers share the extended signature of POLE ultramutation. These findings expand the understanding of ultramutation mechanisms and highlight peculiar mutagenic properties of polypurine tracts in the human genome.

Innovative approach for high-throughput exploiting sex-specific markers in Japanese parrotfish Oplegnathus fasciatus

BACKGROUND

The use of sex-specific molecular markers has become a prominent method in enhancing fish production and economic value, as well as providing a foundation for understanding the complex molecular mechanisms involved in fish sex determination. Over the past decades, research on male and female sex identification has predominantly employed molecular biology methodologies such as restriction fragment length polymorphism, random amplification of polymorphic DNA, simple sequence repeat, and amplified fragment length polymorphism. The emergence of high-throughput sequencing technologies, particularly Illumina, has led to the utilization of single nucleotide polymorphism and insertion/deletion variants as significant molecular markers for investigating sex identification in fish. The advancement of sex-controlled breeding encounters numerous challenges, including the inefficiency of current methods, intricate experimental protocols, high costs of development, elevated rates of false positives, marker instability, and cumbersome field-testing procedures. Nevertheless, the emergence and swift progress of PacBio high-throughput sequencing technology, characterized by its long-read output capabilities, offers novel opportunities to overcome these obstacles.

FINDINGS

Utilizing male/female assembled genome information in conjunction with short-read sequencing data survey and long-read PacBio sequencing data, a catalog of large-segment (>100 bp) insertion/deletion genetic variants was generated through a genome-wide variant site-scanning approach with bidirectional comparisons. The sequence tagging sites were ranked based on the long-read depth of the insertion/deletion site, with markers exhibiting lower long-read depth being considered more effective for large-segment deletion variants. Subsequently, a catalog of bulk primers and simulated PCR for the male/female variant loci was developed, incorporating primer design for the target region and electronic PCR (e-PCR) technology. The Japanese parrotfish (Oplegnathus fasciatus), belonging to the Oplegnathidae family within the Centrarchiformes order, holds significant economic value as a rocky reef fish indigenous to East Asia. The criteria for rapid identification of male and female differences in Japanese parrotfish were established through agarose gel electrophoresis, which revealed 2 amplified bands for males and 1 amplified band for females. A high-throughput identification catalog of sex-specific markers was then constructed using this method, resulting in the identification of 3,639 (2,786 INS/853 DEL, ♀ as reference) and 3,672 (2,876 INS/833 DEL, ♂ as reference) markers in conjunction with 1,021 and 894 high-quality genetic sex identification markers, respectively. Sixteen differential loci were randomly chosen from the catalog for validation, with 11 of them meeting the criteria for male/female distinctions. The implementation of cost-effective and efficient technological processes would facilitate the rapid advancement of genetic breeding through expediting the high-throughput development of sex genetic markers for various species.

CONCLUSIONS

Our study utilized assembled genome information from male and female individuals obtained from PacBio, in addition to data from short-read sequencing data survey and long-read PacBio sequencing data. We extensively employed genome-wide variant site scanning and identification, high-throughput primer design of target regions, and e-PCR batch amplification, along with statistical analysis and ranking of the long-read depth of the variant sites. Through this integrated approach, we successfully compiled a catalog of large insertion/deletion sites (>100 bp) in both male and female Japanese parrotfish.

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.

Revisiting the Role of Genetic Variation in Adaptation

AbstractEvolutionary biologists have thought about the role of genetic variation during adaptation for a very long time-before we understood the organization of the genetic code, the provenance of genetic variation, and how such variation influenced the phenotypes on which natural selection acts. Half a century after the discovery of the structure of DNA and the unraveling of the genetic code, we have a rich understanding of these problems and the means to both delve deeper and widen our perspective across organisms and natural populations. The 2022 Vice Presidential Symposium of the American Society of Naturalists highlighted examples of recent insights into the role of genetic variation in adaptive processes, which are compiled in this special section. The work was conducted in different parts of the world, included theoretical and empirical studies with diverse organisms, and addressed distinct aspects of how genetic variation influences adaptation. In our introductory article to the special section, we discuss some important recent insights about the generation and maintenance of genetic variation, its impacts on phenotype and fitness, its fate in natural populations, and its role in driving adaptation. By placing the special section articles in the broader context of recent developments, we hope that this overview will also serve as a useful introduction to the field.

Genome-wide probing of eukaryotic nascent RNA structure elucidates cotranscriptional folding and its antimutagenic effect

The transcriptional intermediates of RNAs fold into secondary structures with multiple regulatory roles, yet the details of such cotranscriptional RNA folding are largely unresolved in eukaryotes. Here, we present eSPET-seq (Structural Probing of Elongating Transcripts in eukaryotes), a method to assess the cotranscriptional RNA folding in Saccharomyces cerevisiae. Our study reveals pervasive structural transitions during cotranscriptional folding and overall structural similarities between nascent and mature RNAs. Furthermore, a combined analysis with genome-wide R-loop and mutation rate approximations provides quantitative evidence for the antimutator effect of nascent RNA folding through competitive inhibition of the R-loops, known to facilitate transcription-associated mutagenesis. Taken together, we present an experimental evaluation of cotranscriptional folding in eukaryotes and demonstrate the antimutator effect of nascent RNA folding. These results suggest genome-wide coupling between the processing and transmission of genetic information through RNA folding.

Oxidative stress drives mutagenesis through transcription-coupled repair in bacteria

In bacteria, mutations lead to the evolution of antibiotic resistance, which is one of the main public health problems of the twenty-first century. Therefore, determining which cellular processes most frequently contribute to mutagenesis, especially in cells that have not been exposed to exogenous DNA damage, is critical. Here, we show that endogenous oxidative stress is a key driver of mutagenesis and the subsequent development of antibiotic resistance. This is the case for all classes of antibiotics and highly divergent species tested, including patient-derived strains. We show that the transcription-coupled repair pathway, which uses the nucleotide excision repair proteins (TC-NER), is responsible for endogenous oxidative stress-dependent mutagenesis and subsequent evolution. This suggests that a majority of mutations arise through transcription-associated processes rather than the replication fork. In addition to determining that the NER proteins play a critical role in mutagenesis and evolution, we also identify the DNA polymerases responsible for this process. Our data strongly suggest that cooperation between three different mutagenic DNA polymerases, likely at the last step of TC-NER, is responsible for mutagenesis and evolution. Overall, our work identifies a highly conserved pathway that drives mutagenesis due to endogenous oxidative stress, which has broad implications for all diseases of evolution, including antibiotic resistance development.

Rebuilding ships while at sea-Character individuality, homology, and evolutionary innovation

How novel traits originate in evolution is still one of the most perplexing questions in Evolutionary Biology. Building on a previous account of evolutionary innovation, I here propose that evolutionary novelties are those individualized characters that are not homologous to any characters in the ancestor. To clarify this definition, I here provide a detailed analysis of the concepts of "character individuality" and "homology" first, before addressing their role for our understanding of evolutionary innovation. I will argue (1) that functional as well as structural considerations are important for character individualization; and (2) that compositional (structural) and positional homology need to be clearly distinguished to properly describe the evolutionary transformations of hierarchically structured characters. My account will therefore integrate functional and structural perspectives and put forward a new multi-level view of character identity and transformation.

Establishing Linkages Among DNA Damage, Mutagenesis, and Genetic Diseases

DNA damage by chemicals, radiation, or oxidative stress leads to a mutational spectrum, which is complex because it is determined in part by lesion structure, the DNA sequence context of the lesion, lesion repair kinetics, and the type of cells in which the lesion is replicated. Accumulation of mutations may give rise to genetic diseases such as cancer and therefore understanding the process underlying mutagenesis is of immense importance to preserve human health. Chemical or physical agents that cause cancer often leave their mutational fingerprints, which can be used to back-calculate the molecular events that led to disease. To make a clear link between DNA lesion structure and the mutations a given lesion induces, the field of single-lesion mutagenesis was developed. In the last three decades this area of research has seen much growth in several directions, which we attempt to describe in this Perspective.

Differential Susceptibility of Retinal Neurons to the Loss of Mitochondrial Biogenesis Factor Nrf1

The retina, the accessible part of the central nervous system, has served as a model system to study the relationship between energy utilization and metabolite supply. When the metabolite supply cannot match the energy demand, retinal neurons are at risk of death. As the powerhouse of eukaryotic cells, mitochondria play a pivotal role in generating ATP, produce precursors for macromolecules, maintain the redox homeostasis, and function as waste management centers for various types of metabolic intermediates. Mitochondrial dysfunction has been implicated in the pathologies of a number of degenerative retinal diseases. It is well known that photoreceptors are particularly vulnerable to mutations affecting mitochondrial function due to their high energy demand and susceptibility to oxidative stress. However, it is unclear how defective mitochondria affect other retinal neurons. Nuclear respiratory factor 1 (Nrf1) is the major transcriptional regulator of mitochondrial biogenesis, and loss of Nrf1 leads to defective mitochondria biogenesis and eventually cell death. Here, we investigated how different retinal neurons respond to the loss of Nrf1. We provide in vivo evidence that the disruption of Nrf1-mediated mitochondrial biogenesis results in a slow, progressive degeneration of all retinal cell types examined, although they present different sensitivity to the deletion of Nrf1, which implicates differential energy demand and utilization, as well as tolerance to mitochondria defects in different neuronal cells. Furthermore, transcriptome analysis on rod-specific Nrf1 deletion uncovered a previously unknown role of Nrf1 in maintaining genome stability.

Seed banks alter the molecular evolutionary dynamics of Bacillus subtilis

Fluctuations in the availability of resources constrain the growth and reproduction of individuals, which subsequently affects the evolution of their respective populations. Many organisms contend with such fluctuations by entering a reversible state of reduced metabolic activity, a phenomenon known as dormancy. This pool of dormant individuals (i.e. a seed bank) does not reproduce and is expected to act as an evolutionary buffer, though it is difficult to observe this effect directly over an extended evolutionary timescale. Through genetic manipulation, we analyze the molecular evolutionary dynamics of Bacillus subtilis populations in the presence and absence of a seed bank over 700 days. The ability of these bacteria to enter a dormant state increased the accumulation of genetic diversity over time and altered the trajectory of mutations, findings that were recapitulated using simulations based on a mathematical model of evolutionary dynamics. While the ability to form a seed bank did not alter the degree of negative selection, we found that it consistently altered the direction of molecular evolution across genes. Together, these results show that the ability to form a seed bank can affect the direction and rate of molecular evolution over an extended evolutionary timescale.

De Novo L509P Mutation of the TGFBI Gene Associated with Slit-Lamp Findings of Lattice Corneal Dystrophy Type IIIA

Background: Mutations of the transforming growth factor-β-induced (TGFBI) gene produce various types of corneal dystrophy. Here, we report a novel de novo L509P mutation not located in a known hot spot of the transforming growth factor-β-induced (TGFBI) gene and its clinical phenotype, which resembles that of lattice corneal dystrophy type IIIA (LCD IIIA). Case presentation: A 36-year-old man (proband) visited our clinic due to decreased visual acuity with intermittent ocular irritation in conjunction with painful recurrent erosions in both eyes for 10 years. Molecular genetic analyses revealed a TGFBI L509P mutation (c.1526T>C) in the proband and one of his sons. Interestingly, neither TGFBI mutations nor corneal abnormalities were detected in either of the proband’s biological parents, indicating the occurrence of a de novo L509P mutation. Clinical examinations, including slit-lamp retro-illumination and Fourier-domain anterior segment optical coherence tomography (FD-OCT), revealed gray deposits in the anterior stroma and deeper refractile lines extending from limbus to limbus in both corneas of the proband, consistent with a diagnosis of LCD IIIA. Superficial diffuse haze and surface irregularity were observed in conjunction with corneal erosions and visual impairment, necessitating phototherapeutic keratectomy (PTK). A 60 μm PTK of the Bowman layer and anterior stroma of the proband’s left eye was performed following the removal of the epithelium in order to remove superficial corneal opacities. His BCVA improved from 20/400 to 20/50 at postoperative week 8 and was maintained for 45 months. Pinhole-corrected VA was 20/20 at the last visit, and corneal opacities had not recurred. Conclusions: An inheritable de novo mutation of L509P in the TGFBI gene can produce severe LCD IIIA, which can be successfully treated with OCT-guided PRK.

A study on UHPLC-MS/MS analyses of DNA and RNA oxidative damage metabolites in patients with cervical carcinoma: 8-oxoG in urine as a potential biomarker of cervical carcinoma

Objective

Methods

Results

Conclusions

The MDM2 Single-Nucleotide Polymorphism T309G Is Associated With the Development of Epimacular Membranes

Purpose: To investigate the role of the mouse double minute 2 (MDM2) gene single-nucleotide polymorphism (SNP) T309G in the development of epimacular membranes (EMMs) by analyzing the genotype distribution and consistency of the polymorphism in paired membrane-blood samples.

Methods: This was a cross–sectional genetic association study of patients with proliferative vitreoretinopathy (PVR) or EMMs. PVR membranes (PVRMs), internal limiting membranes (ILMs) (PVR-ILMs) and blood samples (PVR-blood) from patients with PVR, and EMMs, EMM-ILMs and EMM-blood from patients with EMMs were collected. The genotype of all samples was determined by Sanger sequencing. Sex composition, mean age, the genotype distribution of MDM2 T309G, the allelic frequency of the MDM2 SNP309 G allele (% G) and the somatic mutation rate at the MDM2 T309G locus (% M) were analyzed and compared. The PVR and healthy Chinese donor groups were used as controls for different comparisons.

Results: The EMM group of 62 patients was older than the PVR group of 61 patients by an average of 8.87 years (p < 0.0001), but the two groups were statistically similar in the sex composition (p = 0.1754). Importantly, G allele carriers were at a higher risk of developing EMMs than non-G allele carriers (p = 0.0479; OR = 2.047). Moreover, EMM-blood exhibited a significantly higher % G than blood samples from healthy Chinese donors (EMM-blood: 56.78%, donors: 45.61%; p = 0.0256; OR = 1.567). Regarding membrane-blood consistency, % M was significantly different between PVRMs and EMMs (PVRMs: 2.63%, EMMs: 21.57%; p = 0.0097; OR = 10.18) but not between different types of ILMs (PVR-ILMs: 18.18%, EMM-ILMs: 29.17%; p = 0.6855). Furthermore, EMMs (p = 0.0053; OR = 8.250) and EMM-ILMs (p = 0.0233; OR = 14.40) from patients with preoperative macular holes were more predisposed toward somatic mutations at the MDM2 T309G locus than those from patients without preoperative macular holes.

Conclusions:MDM2 T309G is associated with the development of EMMs. Herein, the MDM2 SNP309 G allele is first reported as an associated factor of EMMs in a Chinese population. In addition, EMMs and ILMs are genetically unstable at the MDM2 T309G locus, especially when complicated with preoperative macular holes.

Sono-Controllable and ROS-Sensitive CRISPR-Cas9 Genome Editing for Augmented/Synergistic Ultrasound Tumor Nanotherapy

The potential of the cluster regularly interspaced short palindromic repeat (CRISPR)-associated protein 9 (Cas9)-based therapeutic genome editing is severely hampered by the difficulties in precise regulation of the in vivo activity of the CRISPR-Cas9 system. Herein, sono-controllable and reactive oxygen species (ROS)-sensitive sonosensitizer-integrated metal-organic frameworks (MOFs), denoted as P/M@CasMTH1, are developed for augmented sonodynamic therapy (SDT) efficacy using the genome-editing technology. P/M@CasMTH1 nanoparticles comprise singlet oxygen (1 O2 )-generating MOF structures anchored with CRISPR-Cas9 systems via 1 O2 -cleavable linkers, which serve not only as a delivery vector of CRISPR-Cas9 targeting MTH1, but also as a sonoregulator to spatiotemporally activate the genome editing. P/M@CasMTH1 escapes from the lysosomes, harvests the ultrasound (US) energy and converts it into abundant 1 O2 to induce SDT. The generated ROS subsequently trigger cleavage of ROS-responsive thioether bonds, thus inducing controllable release of the CRISPR-Cas9 system and initiation of genome editing. The genomic disruption of MTH1 conspicuously augments the therapeutic efficacy of SDT by destroying the self-defense system in tumor cells, thereby causing cellular apoptosis and tumor suppression. This therapeutic strategy for synergistic MTH1 disruption and abundant 1 O2 generation provides a paradigm for augmenting SDT efficacy based on the emerging nanomedicine-enabled genome-editing technology.

mTOR-Activating Mutations in RRAGD Are Causative for Kidney Tubulopathy and Cardiomyopathy

BACKGROUND

Over the last decade, advances in genetic techniques have resulted in the identification of rare hereditary disorders of renal magnesium and salt handling. Nevertheless, approximately 20% of all patients with tubulopathy lack a genetic diagnosis.

METHODS

We performed whole-exome and -genome sequencing of a patient cohort with a novel, inherited, salt-losing tubulopathy; hypomagnesemia; and dilated cardiomyopathy. We also conducted subsequent in vitro functional analyses of identified variants of RRAGD, a gene that encodes a small Rag guanosine triphosphatase (GTPase).

RESULTS

In eight children from unrelated families with a tubulopathy characterized by hypomagnesemia, hypokalemia, salt wasting, and nephrocalcinosis, we identified heterozygous missense variants in RRAGD that mostly occurred de novo. Six of these patients also had dilated cardiomyopathy and three underwent heart transplantation. We identified a heterozygous variant in RRAGD that segregated with the phenotype in eight members of a large family with similar kidney manifestations. The GTPase RagD, encoded by RRAGD, plays a role in mediating amino acid signaling to the mechanistic target of rapamycin complex 1 (mTORC1). RagD expression along the mammalian nephron included the thick ascending limb and the distal convoluted tubule. The identified RRAGD variants were shown to induce a constitutive activation of mTOR signaling in vitro.

CONCLUSIONS

Our findings establish a novel disease, which we call autosomal dominant kidney hypomagnesemia (ADKH-RRAGD), that combines an electrolyte-losing tubulopathy and dilated cardiomyopathy. The condition is caused by variants in the RRAGD gene, which encodes Rag GTPase D; these variants lead to an activation of mTOR signaling, suggesting a critical role of Rag GTPase D for renal electrolyte handling and cardiac function.

Modeling a cataract disorder in mice with prime editing

Prime editing enables efficient introduction of targeted transversions, insertions, and deletions in mammalian cells and several organisms. However, genetic disease models with base deletions by prime editing have not yet been reported in mice. Here, we successfully generate a mouse model with a cataract disorder through microinjection of prime editor 3 (PE3) plasmids to efficiently induce targeted single-base deletion. Notably, a generated mouse with a high G-deletion rate (38.2%) displays a nuclear cataract phenotype; the PE3-induced deletions in mutant mice achieve high rates of germline transmission to their progenies, with phenotypic inheritance of cataract. Our data propose that modeling a genetic disease with a single nucleotide deletion in mice can be achieved with prime genome editing in vivo.

Principles of seed banks and the emergence of complexity from dormancy

Across the tree of life, populations have evolved the capacity to contend with suboptimal conditions by engaging in dormancy, whereby individuals enter a reversible state of reduced metabolic activity. The resulting seed banks are complex, storing information and imparting memory that gives rise to multi-scale structures and networks spanning collections of cells to entire ecosystems. We outline the fundamental attributes and emergent phenomena associated with dormancy and seed banks, with the vision for a unifying and mathematically based framework that can address problems in the life sciences, ranging from global change to cancer biology.

The Mutational Robustness of the Genetic Code and Codon Usage in Environmental Context: A Non-Extremophilic Preference?

The genetic code was evolved, to some extent, to minimize the effects of mutations. The effects of mutations depend on the amino acid repertoire, the structure of the genetic code and frequencies of amino acids in proteomes. The amino acid compositions of proteins and corresponding codon usages are still under selection, which allows us to ask what kind of environment the standard genetic code is adapted to. Using simple computational models and comprehensive datasets comprising genomic and environmental data from all three domains of Life, we estimate the expected severity of non-synonymous genomic mutations in proteins, measured by the change in amino acid physicochemical properties. We show that the fidelity in these physicochemical properties is expected to deteriorate with extremophilic codon usages, especially in thermophiles. These findings suggest that the genetic code performs better under non-extremophilic conditions, which not only explains the low substitution rates encountered in halophiles and thermophiles but the revealed relationship between the genetic code and habitat allows us to ponder on earlier phases in the history of Life.

Copper Ions Induce DNA Sequence Variation in Zygotic Embryo Culture-Derived Barley Regenerants

In vitro tissue culture could be exploited to study cellular mechanisms that induce sequence variation. Altering the metal ion composition of tissue culture medium affects biochemical pathways involved in tissue culture-induced variation. Copper ions are involved in the mitochondrial respiratory chain and Yang cycle. Copper ions may participate in oxidative mutations, which may contribute to DNA sequence variation. Silver ions compete with copper ions to bind to the complex IV subunit of the respiratory chain, thus affecting the Yang cycle and DNA methylation. The mechanisms underlying somaclonal variation are unknown. In this study, we evaluated embryo-derived barley regenerants obtained from a single double-haploid plant via embryo culture under varying copper and silver ion concentrations and different durations of in vitro culture. Morphological variation among regenerants and the donor plant was not evaluated. Methylation-sensitive Amplified Fragment Length Polymorphism analysis of DNA samples showed DNA methylation pattern variation in CG and CHG (H = A, C, or T) sequence contexts. Furthermore, modification of in vitro culture conditions explained DNA sequence variation, demethylation, and de novo methylation in the CHG context, as indicated by analysis of variance. Linear regression indicated that DNA sequence variation was related to de novo DNA methylation in the CHG context. Mediation analysis showed the role of copper ions as a mediator of sequence variation in the CHG context. No other contexts showed a significant sequence variation in mediation analysis. Silver ions did not act as a mediator between any methylation contexts and sequence variation. Thus, incorporating copper ions in the induction medium should be treated with caution.

Pseudouridines of tRNA Anticodon Stem-Loop Have Unexpected Role in Mutagenesis in Pseudomonas sp

Pseudouridines are known to be important for optimal translation. In this study we demonstrate an unexpected link between pseudouridylation of tRNA and mutation frequency in Pseudomonas species. We observed that the lack of pseudouridylation activity of pseudouridine synthases TruA or RluA elevates the mutation frequency in Pseudomonas putida 3 to 5-fold. The absence of TruA but not RluA elevates mutation frequency also in Pseudomonas aeruginosa. Based on the results of genetic studies and analysis of proteome data, the mutagenic effect of the pseudouridylation deficiency cannot be ascribed to the involvement of error-prone DNA polymerases or malfunctioning of DNA repair pathways. In addition, although the deficiency in TruA-dependent pseudouridylation made P. putida cells more sensitive to antimicrobial compounds that may cause oxidative stress and DNA damage, cultivation of bacteria in the presence of reactive oxygen species (ROS)-scavenging compounds did not eliminate the mutator phenotype. Thus, the elevated mutation frequency in the absence of tRNA pseudouridylation could be the result of a more specific response or, alternatively, of a cumulative effect of several small effects disturbing distinct cellular functions, which remain undetected when studied independently. This work suggests that pseudouridines link the translation machinery to mutation frequency.

Exploring the Biochemical Origin of DNA Sequence Variation in Barley Plants Regenerated via in Vitro Anther Culture

Tissue culture is an essential tool for the regeneration of uniform plant material. However, tissue culture conditions can be a source of abiotic stress for plants, leading to changes in the DNA sequence and methylation patterns. Despite the growing evidence on biochemical processes affected by abiotic stresses, how these altered biochemical processes affect DNA sequence and methylation patterns remains largely unknown. In this study, the methylation-sensitive Amplified Fragment Length Polymorphism (metAFLP) approach was used to investigate de novo methylation, demethylation, and sequence variation in barley regenerants derived by anther culture. Additionally, we used Attenuated Total Reflectance Fourier Transform Infrared (ATR-FTIR) spectroscopy to identify the spectral features of regenerants, which were then analyzed by mediation analysis. The infrared spectrum ranges (710-690 and 1010-940 cm-1) identified as significant in the mediation analysis were most likely related to β-glucans, cellulose, and S-adenosyl-L-methionine (SAM). Additionally, the identified compounds participated as predictors in moderated mediation analysis, explaining the role of demethylation of CHG sites (CHG_DMV) in in vitro tissue culture-induced sequence variation, depending on the duration of tissue culture. The data demonstrate that ATR-FTIR spectroscopy is a useful tool for studying the biochemical compounds that may affect DNA methylation patterns and sequence variation, if combined with quantitative characteristics determined using metAFLP molecular markers and mediation analysis. The role of β-glucans, cellulose, and SAM in DNA methylation, and in cell wall, mitochondria, and signaling, are discussed to highlight the putative cellular mechanisms involved in sequence variation.

Nutritional conditions and oxygen concentration affect spontaneous occurrence of homologous recombination events but not spontaneous mutagenesis in Escherichia coli

Effects of environmental factors for growth of Escherichia coli on spontaneous mutagenesis and homologous recombination events are described. By analyzing rifampicin-resistant (Rif^r) mutation frequencies in an E. coli strain lacking MutM and MutY repair enzymes, which suppress base substitution mutations caused by 8-oxoguanine (7,8 dihydro-8-oxoguanine; 8-oxoG) in DNA, we examined levels of oxidative DNA damage produced in normally growing cells. The level of 8-oxoG DNA damage was about 9- and 63-fold higher in cells grown in M9-glucose and M9-glycerol media, respectively, than in those grown in LB medium. We also found that about 14-fold more 8-oxoG DNA damage was produced in cells grown in about 0.1% oxygen than in those grown in the normal atmosphere. However, Rif^r mutation frequency in wild-type cells was unchanged in such different growth conditions, suggesting that the capacity of repair mechanisms is sufficient to suppress mutations caused by 8-oxoG even at very high levels in cells growing in the particular conditions. On the other hand, the frequency of spontaneous homologous recombination events in wild-type E. coli cells varied with different growth conditions. When cells were grown in M9-glucose and M9-glycerol media, the spontaneous recombination frequency increased to about 4.3- and 7.3-fold, respectively, higher than that in LB medium. Likewise, the spontaneous recombination frequency was about 3.5-fold higher in cells growing in the hypoxic condition than in cells growing in the atmosphere. When cells were grown in anaerobic conditions, the recombination frequency decreased to half of that in the atmosphere. These data indicated that spontaneous homologous recombination is highly responsive to environmental factors such as nutrition and oxygen concentration.

Frameshifting preserves key physicochemical properties of proteins

Frameshifts in protein coding sequences are widely perceived as resulting in either nonfunctional or even deleterious protein products. Indeed, frameshifts typically lead to markedly altered protein sequences and premature stop codons. By analyzing complete proteomes from all three domains of life, we demonstrate that, in contrast, several key physicochemical properties of protein sequences exhibit significant robustness against +1 and -1 frameshifts. In particular, we show that hydrophobicity profiles of many protein sequences remain largely invariant upon frameshifting. For example, over 2,900 human proteins exhibit a Pearson's correlation coefficient R between the hydrophobicity profiles of the original and the +1-frameshifted variants greater than 0.7, despite an average sequence identity between the two of only 6.5% in this group. We observe a similar effect for protein sequence profiles of affinity for certain nucleobases as well as protein sequence profiles of intrinsic disorder. Finally, analysis of significance and optimality demonstrates that frameshift stability is embedded in the structure of the universal genetic code and may have contributed to shaping it. Our results suggest that frameshifting may be a powerful evolutionary mechanism for creating new proteins with vastly different sequences, yet similar physicochemical properties to the proteins from which they originate.

Mutation and Recombination Rates Vary Across Bacterial Chromosome

Bacteria evolve as a result of mutations and acquisition of foreign DNA by recombination processes. A growing body of evidence suggests that mutation and recombination rates are not constant across the bacterial chromosome. Bacterial chromosomal DNA is organized into a compact nucleoid structure which is established by binding of the nucleoid-associated proteins (NAPs) and other proteins. This review gives an overview of recent findings indicating that the mutagenic and recombination processes in bacteria vary at different chromosomal positions. Involvement of NAPs and other possible mechanisms in these regional differences are discussed. Variations in mutation and recombination rates across the bacterial chromosome may have implications in the evolution of bacteria.

Augment of Oxidative Damage with Enhanced Photodynamic Process and MTH1 Inhibition for Tumor Therapy

Tumor cells adapt to reactive oxygen species (ROS) attacking by launching DNA damage repairing mechanisms such as nucleotide pool sanitizing enzyme mutt homologue 1 (MTH1) to mitigate the oxidatively induced DNA lesions, which could greatly limit the therapeutic efficiency of current oxidation therapy. Here, an amplified oxidative damage strategy for tumor therapy was proposed that was focused not only on the enhancement of ROS generation but also the inhibition of subsequent MTH1 enzyme activity simultaneously. In our formulation, mesoporous silica-coated Prussian blue nanoplatforms (PB@MSN) with excellent catalase-like activity and drug loading capability were employed to encapsulate MTH1 inhibitor TH287, followed by the modification of tetraphenylporphrin zinc (Zn-Por) via metallo-supramolecular coordination (PMPT), where Zn-Por behaved as photodynamic and fluorescence imaging agents, as well as acid-responsive gatekeepers. The intelligent PMPT nanosystems could induce the decomposition of H₂O₂ to relieve the hypoxic tumor environment, thus elevating the generation of singlet oxygen for improved oxidative damage. In the meantime, controllable-released TH287 from pores could hinder MTH1-mediated damage repairing process and aggravate oxidative damage, thereby resulting in cellular toxicity as well as tumor growth inhibition.

Identification of spontaneous mutation for broad-spectrum brown planthopper resistance in a large, long-term fast neutron mutagenized rice population

BACKGROUND

The development of rice varieties with broad-spectrum resistance to insect pests is the most promising approach for controlling a fast evolving insect pest such as the brown planthopper (BPH). To cope with rapid evolution, discovering new sources of broad-spectrum resistance genes is the ultimate goal.

RESULTS

We used a forward genetics approach to identify BPH resistance genes in rice (Oryza sativa L.) using double digest restriction site-associated DNA sequencing (ddRADseq) for quantitative trait loci (QTL)-seq of the backcross inbred lines (BILs) derived from a cross between the BPH-susceptible cultivar KDML105 and BPH-resistant cultivar Rathu Heenati (RH). Two major genomic regions, located between 5.78-7.78 Mb (QBPH4.1) and 15.22-17.22 Mb (QBPH4.2) on rice chromosome 4, showed association with BPH resistance in both pooled BILs and individual highly resistant and susceptible BILs. The two most significant candidate resistance genes located within the QBPH4.1 and QBPH4.2 windows were lectin receptor kinase 3 (OsLecRK3) and sesquiterpene synthase 2 (OsSTPS2), respectively. Functional markers identified in these two genes were used for reverse screening 9323 lines of the fast neutron (FN)-mutagenized population developed from the BPH-susceptible, purple-pigmented, indica cultivar Jao Hom Nin (JHN). Nineteen FN-mutagenized lines (0.24%) carried mutations in the OsLecRK3 and/or OsSTPS2 gene. Among these mutant lines, only one highly resistant line (JHN4) and three moderately resistant lines (JHN09962, JHN12005, and JHN19525) were identified using three active, local BPH populations. The 19 mutant lines together with three randomly selected mutant lines, which did not harbor mutations in the two target genes, were screened further for mutations in six known BPH resistance genes including BPH9, BPH14, BPH18, BPH26, BPH29, and BPH32. Multiple single nucleotide polymorphisms (SNPs) and insertion-deletion (Indel) mutations were identified, which formed gene-specific haplotype patterns (HPs) essential for broad-spectrum resistance to BPH in both BILs and JHN mutant populations.

CONCLUSION

On the one hand, HPs of OsLekRK2-3, OsSTPS2, and BPH32 determined broad-spectrum resistance to BPH among RH-derived BILs. On the other hand, in the JHN mutant population, BPH9 together with seven significant genes on chromosome 4 played a crucial role in BPH resistance.

The non-random patterns of genetic variation induced by asymmetric somatic hybridization in wheat

BACKGROUND

Asymmetric somatic hybridization is an efficient crop breeding approach by introducing several exogenous chromatin fragments, which leads to genomic shock and therefore induces genome-wide genetic variation. However, the fundamental question concerning the genetic variation such as whether it occurs randomly and suffers from selection pressure remains unknown.

RESULTS

Here, we explored this issue by comparing expressed sequence tags of a common wheat cultivar and its asymmetric somatic hybrid line. Both nucleotide substitutions and indels (insertions and deletions) had lower frequencies in coding sequences than in un-translated regions. The frequencies of nucleotide substitutions and indels were both comparable between chromosomes with and without introgressed fragments. Nucleotide substitutions distributed unevenly and were preferential to indel-flanking sequences, and the frequency of nucleotide substitutions at 5'-flanking sequences of indels was obviously higher in chromosomes with introgressed fragments than in those without exogenous fragment. Nucleotide substitutions and indels both had various frequencies among seven groups of allelic chromosomes, and the frequencies of nucleotide substitutions were strongly negatively correlative to those of indels. Among three sets of genomes, the frequencies of nucleotide substitutions and indels were both heterogeneous, and the frequencies of nucleotide substitutions exhibited drastically positive correlation to those of indels.

CONCLUSIONS

Our work demonstrates that the genetic variation induced by asymmetric somatic hybridization is attributed to both whole genomic shock and local chromosomal shock, which is a predetermined and non-random genetic event being closely associated with selection pressure. Asymmetric somatic hybrids provide a worthwhile model to further investigate the nature of genomic shock induced genetic variation.

Reduced intrinsic DNA curvature leads to increased mutation rate

BACKGROUND

Mutation rates vary across the genome. Many trans factors that influence mutation rates have been identified, as have specific sequence motifs at the 1-7-bp scale, but cis elements remain poorly characterized. The lack of understanding regarding why different sequences have different mutation rates hampers our ability to identify positive selection in evolution and to identify driver mutations in tumorigenesis.

RESULTS

Here, we use a combination of synthetic genes and sequences of thousands of isolated yeast colonies to show that intrinsic DNA curvature is a major cis determinant of mutation rate. Mutation rate negatively correlates with DNA curvature within genes, and a 10% decrease in curvature results in a 70% increase in mutation rate. Consistently, both yeast and humans accumulate mutations in regions with small curvature. We further show that this effect is due to differences in the intrinsic mutation rate, likely due to differences in mutagen sensitivity and not due to differences in the local activity of DNA repair.

CONCLUSIONS

Our study establishes a framework for understanding the cis properties of DNA sequence in modulating the local mutation rate and identifies a novel causal source of non-uniform mutation rates across the genome.

GC content elevates mutation and recombination rates in the yeast Saccharomyces cerevisiae

The chromosomes of many eukaryotes have regions of high GC content interspersed with regions of low GC content. In the yeast Saccharomyces cerevisiae, high-GC regions are often associated with high levels of meiotic recombination. In this study, we constructed URA3 genes that differ substantially in their base composition [URA3-AT (31% GC), URA3-WT (43% GC), and URA3-GC (63% GC)] but encode proteins with the same amino acid sequence. The strain with URA3-GC had an approximately sevenfold elevated rate of ura3 mutations compared with the strains with URA3-WT or URA3-AT About half of these mutations were single-base substitutions and were dependent on the error-prone DNA polymerase ζ. About 30% were deletions or duplications between short (5-10 base) direct repeats resulting from DNA polymerase slippage. The URA3-GC gene also had elevated rates of meiotic and mitotic recombination relative to the URA3-AT or URA3-WT genes. Thus, base composition has a substantial effect on the basic parameters of genome stability and evolution.

DNA damage and epigenetic alteration in soybean farmers exposed to complex mixture of pesticides

Exposure to pesticides can trigger genotoxic and mutagenic processes through different pathways. However, epidemiological studies are scarce, and further work is needed to find biomarkers sensitive to the health of exposed populations. Considering that there are few evaluations of soybean farmers, the aim of this study was to assess the effects of human exposure to complex mixtures of pesticides. The alkaline comet assay modified with restriction enzyme (hOGG1: human 8-oxoguanine DNA glycosylase) was used to detect oxidised guanine, and compared with the buccal micronucleus cytome assay, global methylation, haematological parameters, biochemical analyses (serum glutamic oxaloacetic transaminase, serum glutamic pyruvic transaminase, gamma-glutamyl-transferase and butyrylcholinesterase), and particle-induced X-ray emission (PIXE) for the analysis of inorganic elements. Farm workers (n = 137) exposed to different types of pesticides were compared with a non-exposed reference group (control; n = 83). Results of the enzyme-modified comet assay suggest oxidation of guanine in DNA generated by pesticides exposure. It was observed that DNA damage (comet assay and micronucleus test) was significantly increased in exposed individuals compared to the unexposed group. The micronucleus test demonstrated elimination of nuclear material by budding, defective cytokinesis and dead cells. Occupationally exposed individuals also showed genomic hypermethylation of DNA, which correlated with micronucleus frequency. No differences were detected regarding the haematological and biochemical parameters. Finally, significantly higher concentrations of Al and P were observed in the urine of the soybean farmers. DNA damage could be a consequence of the ability of the complex mixture, including Al and P, to cause oxidative damage. These data indicate that persistent genetic instability associated with hypermethylation of DNA in soybean workers after long-term exposure to a low-level to pesticides mixtures may be critical for the development of adverse health effects such as cancer.

The Emergence of Canalization and Evolvability in an Open-Ended, Interactive Evolutionary System

Many believe that an essential component for the discovery of the tremendous diversity in natural organisms was the evolution of evolvability, whereby evolution speeds up its ability to innovate by generating a more adaptive pool of offspring. One hypothesized mechanism for evolvability is developmental canalization, wherein certain dimensions of variation become more likely to be traversed and others are prevented from being explored (e.g., offspring tend to have similar-size legs, and mutations affect the length of both legs, not each leg individually). While ubiquitous in nature, canalization is rarely reported in computational simulations of evolution, which deprives us of in silico examples of canalization to study and raises the question of which conditions give rise to this form of evolvability. Answering this question would shed light on why such evolvability emerged naturally, and it could accelerate engineering efforts to harness evolution to solve important engineering challenges. In this article, we reveal a unique system in which canalization did emerge in computational evolution. We document that genomes entrench certain dimensions of variation that were frequently explored during their evolutionary history. The genetic representation of these organisms also evolved to be more modular and hierarchical than expected by chance, and we show that these organizational properties correlate with increased fitness. Interestingly, the type of computational evolutionary experiment that produced this evolvability was very different from traditional digital evolution in that there was no objective, suggesting that open-ended, divergent evolutionary processes may be necessary for the evolution of evolvability.

Spontaneous mutation rate is a plastic trait associated with population density across domains of life

Rates of random, spontaneous mutation can vary plastically, dependent upon the environment. Such plasticity affects evolutionary trajectories and may be adaptive. We recently identified an inverse plastic association between mutation rate and population density at 1 locus in 1 species of bacterium. It is unknown how widespread this association is, whether it varies among organisms, and what molecular mechanisms of mutagenesis or repair are required for this mutation-rate plasticity. Here, we address all 3 questions. We identify a strong negative association between mutation rate and population density across 70 years of published literature, comprising hundreds of mutation rates estimated using phenotypic markers of mutation (fluctuation tests) from all domains of life and viruses. We test this relationship experimentally, determining that there is indeed density-associated mutation-rate plasticity (DAMP) at multiple loci in both eukaryotes and bacteria, with up to 23-fold lower mutation rates at higher population densities. We find that the degree of plasticity varies, even among closely related organisms. Nonetheless, in each domain tested, DAMP requires proteins scavenging the mutagenic oxidised nucleotide 8-oxo-dGTP. This implies that phenotypic markers give a more precise view of mutation rate than previously believed: having accounted for other known factors affecting mutation rate, controlling for population density can reduce variation in mutation-rate estimates by 93%. Widespread DAMP, which we manipulate genetically in disparate organisms, also provides a novel trait to use in the fight against the evolution of antimicrobial resistance. Such a prevalent environmental association and conserved mechanism suggest that mutation has varied plastically with population density since the early origins of life.

Dynamics and genetic diversification of Escherichia coli during experimental adaptation to an anaerobic environment

Background

Many bacteria are facultative anaerobes, and can proliferate in both anoxic and oxic environments. Under anaerobic conditions, fermentation is the primary means of energy generation in contrast to respiration. Furthermore, the rates and spectra of spontaneous mutations that arise during anaerobic growth differ to those under aerobic growth. A long-term selection experiment was undertaken to investigate the genetic changes that underpin how the facultative anaerobe, Escherichia coli, adapts to anaerobic environments.

Methods

Twenty-one populations of E. coli REL4536, an aerobically evolved 10,000th generation descendent of the E. coli B strain, REL606, were established from a clonal ancestral culture. These were serially sub-cultured for 2,000 generations in a defined minimal glucose medium in strict aerobic and strict anaerobic environments, as well as in a treatment that fluctuated between the two environments. The competitive fitness of the evolving lineages was assessed at approximately 0, 1,000 and 2,000 generations, in both the environment of selection and the alternative environment. Whole genome re-sequencing was performed on random colonies from all lineages after 2,000-generations. Mutations were identified relative to the ancestral genome, and based on the extent of parallelism, traits that were likely to have contributed towards adaptation were inferred.

Results

There were increases in fitness relative to the ancestor among anaerobically evolved lineages when tested in the anaerobic environment, but no increases were found in the aerobic environment. For lineages that had evolved under the fluctuating regime, relative fitness increased significantly in the anaerobic environment, but did not increase in the aerobic environment. The aerobically-evolved lineages did not increase in fitness when tested in either the aerobic or anaerobic environments. The strictly anaerobic lineages adapted more rapidly to the anaerobic environment than did the fluctuating lineages. Two main strategies appeared to predominate during adaptation to the anaerobic environment: modification of energy generation pathways, and inactivation of non-essential functions. Fermentation pathways appeared to alter through selection for mutations in genes such as nadR, adhE, dcuS/R, and pflB. Mutations were frequently identified in genes for presumably dispensable functions such as toxin-antitoxin systems, prophages, virulence and amino acid transport. Adaptation of the fluctuating lineages to the anaerobic environments involved mutations affecting traits similar to those observed in the anaerobically evolved lineages.

Discussion

There appeared to be strong selective pressure for activities that conferred cell yield advantages during anaerobic growth, which include restoring activities that had previously been inactivated under long-term continuous aerobic evolution of the ancestor.

New insights into the generation and role of de novo mutations in health and disease

Aside from inheriting half of the genome of each of our parents, we are born with a small number of novel mutations that occurred during gametogenesis and postzygotically. Recent genome and exome sequencing studies of parent-offspring trios have provided the first insights into the number and distribution of these de novo mutations in health and disease, pointing to risk factors that increase their number in the offspring. De novo mutations have been shown to be a major cause of severe early-onset genetic disorders such as intellectual disability, autism spectrum disorder, and other developmental diseases. In fact, the occurrence of novel mutations in each generation explains why these reproductively lethal disorders continue to occur in our population. Recent studies have also shown that de novo mutations are predominantly of paternal origin and that their number increases with advanced paternal age. Here, we review the recent literature on de novo mutations, covering their detection, biological characterization, and medical impact.

Interpreting the Dependence of Mutation Rates on Age and Time

Mutations can originate from the chance misincorporation of nucleotides during DNA replication or from DNA lesions that arise between replication cycles and are not repaired correctly. We introduce a model that relates the source of mutations to their accumulation with cell divisions, providing a framework for understanding how mutation rates depend on sex, age, and cell division rate. We show that the accrual of mutations should track cell divisions not only when mutations are replicative in origin but also when they are non-replicative and repaired efficiently. One implication is that observations from diverse fields that to date have been interpreted as pointing to a replicative origin of most mutations could instead reflect the accumulation of mutations arising from endogenous reactions or exogenous mutagens. We further find that only mutations that arise from inefficiently repaired lesions will accrue according to absolute time; thus, unless life history traits co-vary, the phylogenetic “molecular clock” should not be expected to run steadily across species.

Modeling how the source of mutations relates to their rate of accumulation with age, sex, and number of cell divisions helps to explain perplexing observations about germline and somatic mutations.

We relate how mutations arise to how they accumulate in different sexes, with age and with cell division. This model provides a single framework within which to interpret emerging results from evolutionary biology, human genetics, and cancer genetics. We show that the accrual of mutations should track cell divisions not only when mutations originate during DNA replication but also when they arise through non-replicative mechanisms and are repaired efficiently. This realization means that previous observations of correlations between mutation and cell division rates actually provide little support to the commonly held belief that most germline and somatic mutations arise from replication errors. We further find that only mutations that arise from inefficiently repaired lesions will accrue according to absolute time; thus, without covariation in life history traits, the phylogenetic “molecular clock” should not be expected to run at constant rates across species.

Proton tunneling in the A∙T Watson-Crick DNA base pair: myth or reality?

The results and conclusions reached by Godbeer et al. in their recent work, that proton tunneling in the A∙T(WC) Watson-Crick (WC) DNA base pair occurs according to the Löwdin's (L) model, but with a small (~10(-9)) probability were critically analyzed. Here, it was shown that this finding overestimates the possibility of the proton tunneling at the A∙T(WC)↔A*∙T*(L) tautomerization, because this process cannot be implemented as a chemical reaction. Furthermore, it was outlined those biologically important nucleobase mispairs (A∙A*↔A*∙A, G∙G*↔G*∙G, T∙T*↔T*∙T, C∙C*↔C*∙C, H∙H*↔H*∙H (H - hypoxanthine)) - the players in the field of the spontaneous point mutagenesis - where the tunneling of protons is expected and for which the application of the model proposed by Godbeer et al. can be productive.

MutT from the fish pathogen Aliivibrio salmonicida is a cold-active nucleotide-pool sanitization enzyme with unexpectedly high thermostability

Upon infection by pathogenic bacteria, production of reactive oxygen species (ROS) is part of the host organism's first line of defence. ROS damage a number of macromolecules, and in order to withstand such a harsh environment, the bacteria need to have well-functioning ROS scavenging and repair systems. Herein, MutT is an important nucleotide-pool sanitization enzyme, which degrades 8-oxo-dGTP and thus prevents it from being incorporated into DNA. In this context, we have performed a comparative biochemical and structural analysis of MutT from the fish pathogen Aliivibrio salmonicida (AsMutT) and the human pathogen Vibrio cholerae (VcMutT), in order to analyse their function as nucleotide sanitization enzymes and also determine possible cold-adapted properties of AsMutT. The biochemical characterisation revealed that both enzymes possess activity towards the 8-oxo-dGTP substrate, and that AsMutT has a higher catalytic efficiency than VcMutT at all temperatures studied. Calculations based on the biochemical data also revealed a lower activation energy (E a) for AsMutT compared to VcMutT, and differential scanning calorimetry experiments showed that AsMutT displayed an unexpected higher melting temperature (T m) value than VcMutT. A comparative analysis of the crystal structure of VcMutT, determined to 2.42 Å resolution, and homology models of AsMutT indicate that three unique Gly residues in loops of VcMutT, and additional long range ion-pairs in AsMutT could explain the difference in temperature stability of the two enzymes. We conclude that AsMutT is a stable, cold-active enzyme with high catalytic efficiency and reduced E a, compared to the mesophilic VcMutT.

A mutant ('lab strain') of the hyperthermophilic archaeon Pyrococcus furiosus, lacking flagella, has unusual growth physiology

A mutant ('lab strain') of the hyperthermophilic archaeon Pyrococcus furiosus DSM3638 exhibited an extended exponential phase and atypical cell aggregation behavior. Genomic DNA from the mutant culture was sequenced and compared to wild-type (WT) DSM3638, revealing 145 genes with one or more insertions, deletions, or substitutions (12 silent, 33 amino acid substitutions, and 100 frame shifts). Approximately, half of the mutated genes were transposases or hypothetical proteins. The WT transcriptome revealed numerous changes in amino acid and pyrimidine biosynthesis pathways coincidental with growth phase transitions, unlike the mutant whose transcriptome reflected the observed prolonged exponential phase. Targeted gene deletions, based on frame-shifted ORFs in the mutant genome, in a genetically tractable strain of P. furiosus (COM1) could not generate the extended exponential phase behavior observed for the mutant. For example, a putative radical SAM family protein (PF2064) was the most highly up-regulated ORF (>25-fold) in the WT between exponential and stationary phase, although this ORF was unresponsive in the mutant; deletion of this gene in P. furiosus COM1 resulted in no apparent phenotype. On the other hand, frame-shifting mutations in the mutant genome negatively impacted transcription of a flagellar biosynthesis operon (PF0329-PF0338).Consequently, cells in the mutant culture lacked flagella and, unlike the WT, showed minimal evidence of exopolysaccharide-based cell aggregation in post-exponential phase. Electron microscopy of PF0331-PF0337 deletions in P. furiosus COM1 showed that absence of flagella impacted normal cell aggregation behavior and, furthermore, indicated that flagella play a key role, beyond motility, in the growth physiology of P. furiosus.

MUTYH, an adenine DNA glycosylase, mediates p53 tumor suppression via PARP-dependent cell death

p53-regulated caspase-independent cell death has been implicated in suppression of tumorigenesis, however, the regulating mechanisms are poorly understood. We previously reported that 8-oxoguanine (8-oxoG) accumulation in nuclear DNA (nDNA) and mitochondrial DNA triggers two distinct caspase-independent cell death through buildup of single-strand DNA breaks by MutY homolog (MUTYH), an adenine DNA glycosylase. One pathway depends on poly-ADP-ribose polymerase (PARP) and the other depends on calpains. Deficiency of MUTYH causes MUTYH-associated familial adenomatous polyposis. MUTYH thereby suppresses tumorigenesis not only by avoiding mutagenesis, but also by inducing cell death. Here, we identified the functional p53-binding site in the human MUTYH gene and demonstrated that MUTYH is transcriptionally regulated by p53, especially in the p53/DNA mismatch repair enzyme, MLH1-proficient colorectal cancer-derived HCT116+Chr3 cells. MUTYH-small interfering RNA, an inhibitor for p53 or PARP suppressed cell death without an additive effect, thus revealing that MUTYH is a potential mediator of p53 tumor suppression, which is known to be upregulated by MLH1. Moreover, we found that the p53-proficient, mismatch repair protein, MLH1-proficient colorectal cancer cell line express substantial levels of MUTYH in nuclei but not in mitochondria, suggesting that 8-oxoG accumulation in nDNA triggers MLH1/PARP-dependent cell death. These results provide new insights on the molecular mechanism of tumorigenesis and potential new strategies for cancer therapies.

Cellular levels of 8-oxoguanine in either DNA or the nucleotide pool play pivotal roles in carcinogenesis and survival of cancer cells

8-Oxoguanine, a major oxidized base lesion formed by reactive oxygen species, causes G to T transversion mutations or leads to cell death in mammals if it accumulates in DNA. 8-Oxoguanine can originate as 8-oxo-dGTP, formed in the nucleotide pool, or by direct oxidation of the DNA guanine base. MTH1, also known as NUDT1, with 8-oxo-dGTP hydrolyzing activity, 8-oxoguanine DNA glycosylase (OGG1) an 8-oxoG DNA glycosylase, and MutY homolog (MUTYH) with adenine DNA glycosylase activity, minimize the accumulation of 8-oxoG in DNA; deficiencies in these enzymes increase spontaneous and induced tumorigenesis susceptibility. However, different tissue types have different tumorigenesis susceptibilities. These can be reversed by combined deficiencies in the defense systems, because cell death induced by accumulation of 8-oxoG in DNA is dependent on MUTYH, which can be suppressed by MTH1 and OGG1. In cancer cells encountering high oxidative stress levels, a high level of 8-oxo-dGTP accumulates in the nucleotide pool, and cells therefore express increased levels of MTH1 in order to eliminate 8-oxo-dGTP. Suppression of MTH1 may be an efficient strategy for killing cancer cells; however, because MTH1 and OGG1 protect normal tissues from oxidative-stress-induced cell death, it is important that MTH1 inhibition does not increase the risk of healthy tissue degeneration.

Engineering reduced evolutionary potential for synthetic biology

The field of synthetic biology seeks to engineer reliable and predictable behaviors in organisms from collections of standardized genetic parts. However, unlike other types of machines, genetically encoded biological systems are prone to changes in their designed sequences due to mutations in their DNA sequences after these devices are constructed and deployed. Thus, biological engineering efforts can be confounded by undesired evolution that rapidly breaks the functions of parts and systems, particularly when they are costly to the host cell to maintain. Here, we explain the fundamental properties that determine the evolvability of biological systems. Then, we use this framework to review current efforts to engineer the DNA sequences that encode synthetic biology devices and the genomes of their microbial hosts to reduce their ability to evolve and therefore increase their genetic reliability so that they maintain their intended functions over longer timescales.