Mitochondrial mutational spectrum is associated with mammalian longevity: a novel signature of oxidative damage.

The mutational spectrum of the mitochondrial DNA (mtDNA) does not resemble any of the known mutational signatures of the nuclear genome and variation in mtDNA mutational spectra between different tissues and organisms is still incomprehensible. Since mitochondria is tightly involved in aerobic energy production, it is expected that mtDNA mutational spectra may be affected by the oxidative damage which is increasing with organismal aging. However, the well-documented mutational signature of the oxidative damage, G>T substitutions, is typical only for the nuclear genome while it is extremely rare in mtDNA. Thus it is still unclear if there is a mitochondria-specific mutational signature of the oxidative damage. Here, reconstructing mtDNA mutational spectra for 424 mammalian species with variable generation length which is a proxy for oocyte age, we observed that the frequency of AH>GH substitutions (H - heavy chain notation) is positively correlated with organismal longevity. This mutational bias from AH to GH significantly affected the nucleotide content of analyzed 650 complete mammalian mitochondrial genomes, where fourfold degenerative synonymous positions of long-lived species become more AH poor and GH rich. Because (i) A>G is a substitution, typical for mtDNA; (ii) it is characterized by very strong asymmetry: A>G is several-fold more frequent on a heavy chain as compared to the light one; (iii) it is sensitive to the time being single-stranded during mtDNA asynchronous replication; (iv) it is associated with oxidative damage of single-stranded DNA in recent experimental studies we propose that A>G is a novel mutational signature of age-associated oxidative damage of single-stranded mtDNA. The described association of the mtDNA mutational spectra with a species-specific life-history trait can significantly affect general patterns of molecular evolution of mtDNA.

Changes in the distribution of fitness effects and adaptive mutational spectra following a single first step towards adaptation

Historical contingency and diminishing returns epistasis have been typically studied for relatively divergent genotypes and/or over long evolutionary timescales. Here, we use Saccharomyces cerevisiae to study the extent of diminishing returns and the changes in the adaptive mutational spectra following a single first adaptive mutational step. We further evolve three clones that arose under identical conditions from a common ancestor. We follow their evolutionary dynamics by lineage tracking and determine adaptive outcomes using fitness assays and whole genome sequencing. We find that diminishing returns manifests as smaller fitness gains during the 2nd step of adaptation compared to the 1st step, mainly due to a compressed distribution of fitness effects. We also find that the beneficial mutational spectra for the 2nd adaptive step are contingent on the 1st step, as we see both shared and diverging adaptive strategies. Finally, we find that adaptive loss-of-function mutations, such as nonsense and frameshift mutations, are less common in the second step of adaptation than in the first step.

Sequence Dependent Repair of 1,N6-Ethenoadenine by DNA Repair Enzymes ALKBH2, ALKBH3, and AlkB

Mutation patterns of DNA adducts, such as mutational spectra and signatures, are useful tools for diagnostic and prognostic purposes. Mutational spectra of carcinogens derive from three sources: adduct formation, replication bypass, and repair. Here, we consider the repair aspect of 1,N6-ethenoadenine (εA) by the 2-oxoglutarate/Fe(II)-dependent AlkB family enzymes. Specifically, we investigated εA repair across 16 possible sequence contexts (5′/3′ flanking base to εA varied as G/A/T/C). The results revealed that repair efficiency is altered according to sequence, enzyme, and strand context (ss- versus ds-DNA). The methods can be used to study other aspects of mutational spectra or other pathways of repair.

Mutability of mononucleotide repeats, not oxidative stress, explains the discrepancy between laboratory-accumulated mutations and the natural allele-frequency spectrum in C. elegans

Important clues about natural selection can be gleaned from discrepancies between the properties of segregating genetic variants and of mutations accumulated experimentally under minimal selection, provided the mutational process is the same in the laboratory as in nature. The base-substitution spectrum differs between C. elegans laboratory mutation accumulation (MA) experiments and the standing site-frequency spectrum, which has been argued to be in part owing to increased oxidative stress in the laboratory environment. Using genome sequence data from C. elegans MA lines carrying a mutation (mev-1) that increases the cellular titer of reactive oxygen species (ROS), leading to increased oxidative stress, we find the base-substitution spectrum is similar between mev-1, its wild-type progenitor (N2), and another set of MA lines derived from a different wild strain (PB306). Conversely, the rate of short insertions is greater in mev-1, consistent with studies in other organisms in which environmental stress increased the rate of insertion–deletion mutations. Further, the mutational properties of mononucleotide repeats in all strains are different from those of nonmononucleotide sequence, both for indels and base-substitutions, and whereas the nonmononucleotide spectra are fairly similar between MA lines and wild isolates, the mononucleotide spectra are very different, with a greater frequency of A:T → T:A transversions and an increased proportion of ±1-bp indels. The discrepancy in mutational spectra between laboratory MA experiments and natural variation is likely owing to a consistent (but unknown) effect of the laboratory environment that manifests itself via different modes of mutability and/or repair at mononucleotide loci.

Patterns of within-host genetic diversity in SARS-CoV-2

Monitoring the spread of SARS-CoV-2 and reconstructing transmission chains has become a major public health focus for many governments around the world. The modest mutation rate and rapid transmission of SARS-CoV-2 prevents the reconstruction of transmission chains from consensus genome sequences, but within-host genetic diversity could theoretically help identify close contacts. Here we describe the patterns of within-host diversity in 1181 SARS-CoV-2 samples sequenced to high depth in duplicate. 95.1% of samples show within-host mutations at detectable allele frequencies. Analyses of the mutational spectra revealed strong strand asymmetries suggestive of damage or RNA editing of the plus strand, rather than replication errors, dominating the accumulation of mutations during the SARS-CoV-2 pandemic. Within- and between-host diversity show strong purifying selection, particularly against nonsense mutations. Recurrent within-host mutations, many of which coincide with known phylogenetic homoplasies, display a spectrum and patterns of purifying selection more suggestive of mutational hotspots than recombination or convergent evolution. While allele frequencies suggest that most samples result from infection by a single lineage, we identify multiple putative examples of co-infection. Integrating these results into an epidemiological inference framework, we find that while sharing of within-host variants between samples could help the reconstruction of transmission chains, mutational hotspots and rare cases of superinfection can confound these analyses.

In vivo evolution of an emerging zoonotic bacterial pathogen in an immunocompromised human host

Zoonotic transfer of animal pathogens to human hosts can generate novel agents, but the genetic events following such host jumps are not well studied. Here we characterize the mechanisms driving adaptive evolution of the emerging zoonotic pathogen Bordetella hinzii in a patient with interleukin-12 receptor β1 deficiency. Genomic sequencing of 24 B. hinzii isolates cultured from blood and stool over 45 months revealed a clonal lineage that had undergone extensive within-host genetic and phenotypic diversification. Twenty of 24 isolates shared an E9G substitution in the DNA polymerase III ε-subunit active site, resulting in a proofreading deficiency. Within this proofreading-deficient clade, multiple lineages with mutations in DNA repair genes and altered mutational spectra emerged and dominated clinical cultures for more than 12 months. Multiple enzymes of the tricarboxylic acid cycle and gluconeogenesis pathways were repeatedly mutated, suggesting rapid metabolic adaptation to the human environment. Furthermore, an excess of G:C > T:A transversions suggested that oxidative stress shaped genetic diversification during adaptation. We propose that inactivation of DNA proofreading activity in combination with prolonged, but sub-lethal, oxidative attack resulting from the underlying host immunodeficiency facilitated rapid genomic adaptation. These findings suggest a fundamental role for host immune phenotype in shaping pathogen evolution following zoonotic infection.

Modeling sequence-space exploration and emergence of epistatic signals in protein evolution

During their evolution, proteins explore sequence space via an interplay between random mutations and phenotypic selection. Here we build upon recent progress in reconstructing data-driven fitness landscapes for families of homologous proteins, to propose stochastic models of experimental protein evolution. These models predict quantitatively important features of experimentally evolved sequence libraries, like fitness distributions and position-specific mutational spectra. They also allow us to efficiently simulate sequence libraries for a vast array of combinations of experimental parameters like sequence divergence, selection strength and library size. We showcase the potential of the approach in re-analyzing two recent experiments to determine protein structure from signals of epistasis emerging in experimental sequence libraries. To be detectable, these signals require sufficiently large and sufficiently diverged libraries. Our modeling framework offers a quantitative explanation for the variable success of recently published experiments. Furthermore, we can forecast the outcome of time- and resource-intensive evolution experiments, opening thereby a way to computationally optimize experimental protocols.

Mutation hotspots, geographical and temporal distribution of SARS-CoV-2 lineages in Brazil, February 2020 to February 2021: insights and limitations from uneven sequencing efforts

The COVID-19 pandemic has already reached approximately 110 million people and it is associated with 2.5 million deaths worldwide. Brazil is the third worst-hit country, with approximately 10.2 million cases and 250 thousand deaths. International efforts have been established to share information about SARS-CoV-2 epidemiology and evolution. However, sequencing facilities and research investments are very heterogeneous across different regions and countries. The understanding of the SARS-CoV-2 evolution plays a significant role in the development of effective strategies for public health and disease management. We aimed to analyze the available and high-quality genome sequences from Brazil between February 2020 and February 2021 to identify mutation hotspots, geographical and temporal distribution of SARS-CoV-2 lineages by using phylogenetics and phylodynamics analyses. We describe heterogeneous and episodic sequencing efforts, the progression of the different lineages along time, evaluating mutational spectra and frequency oscillations derived from the prevalence of novel and specific lineages across different Brazilian regions. We found at least seven major (1-7) and two minor clades (4.2 and 5.3) related to the six most prevalent lineages in the country and described its spatial distribution and dynamics. The emergence and recent frequency shift of lineages (P.1 and P.2) containing mutations of concern in the spike protein (e. g., E484K, N501Y) draws attention due to their association with immune evasion and enhanced receptor binding affinity. Improvements in genomic surveillance are of paramount importance and should be extended in Brazil to better inform policy makers and enable evidence-based decisions to fight the COVID-19 pandemic.

Case Report: Detection of a Novel Germline PALB2 Deletion in a Young Woman With Hereditary Breast Cancer: When the Patient's Phenotype History Doesn't Lie

The partner and localizer of BRCA2 (PALB2) is a major BRCA2 binding partner that participates in homologous recombination repair in response to DNA double-strand breaks. Germline alterations of the PALB2 gene have recently been associated with a high risk of developing breast cancer. We investigated a 37-year-old Caucasian woman with breast cancer and family history of breast cancer using targeted next generation sequencing. A novel heterozygous deletion involving exons 5 and 6 was found in the PALB2 gene, and resulted in the production of a truncated PALB2 protein. These findings expand the mutational spectra of PALB2-associated breast cancer, and may improve the mutation-based screening and genetic diagnosis of breast cancer.

Mutability of mononucleotide repeats, not oxidative stress, explains the discrepancy between laboratory-accumulated mutations and the natural allele-frequency spectrum in C. elegans

Important clues about natural selection can be gleaned from discrepancies between the properties of segregating genetic variants and of mutations accumulated experimentally under minimal selection, provided the mutational process is the same in the lab as in nature. The ratio of transitions to transversions (Ts/Tv) is consistently lower in C. elegans mutation accumulation (MA) experiments than in nature, which has been argued to be in part due to increased oxidative stress in the lab environment. Using whole-genome sequence data from a set of C. elegans MA lines carrying a mutation (mev-1) that increases the cellular titer of reactive oxygen species (ROS), leading to increased endogenous oxidative stress, we find that the base-substitution spectrum is similar between mev-1 lines, its wild-type progenitor (N2), and another set of MA lines derived from a different wild strain (PB306). By contrast, the rate of short insertions is greater in the mev-1 lines, consistent with studies in other organisms in which environmental stress led to an increase in the rate of insertion-deletion mutations. Further, the mutational properties of mononucleotide repeats in all strains are qualitatively different from those of non-mononucleotide sequence, both for indels and base-substitutions, and whereas the non-mononucleotide spectra are fairly similar between MA lines and wild isolates, the mononucleotide spectra are very different. The discrepancy in mutational spectra between lab MA experiments and natural variation is likely due to a consistent (but unknown) effect of the lab environment that manifests itself via different modes of mutability and/or repair at mononucleotide loci.