Genetic barriers to gene flow separate divergent substitution rates across a butterfly hybrid zone

Substitution rate defines the fundamental timescale of molecular evolution which often varies in a species-specific manner. However, it is unknown under what conditions lineage-specific rates can be preserved between natural populations with frequent hybridization. Here, we show in a hybrid zone between two butterflies, Papilio syfanius and Papilio maackii, that genome-wide barriers to gene flow can effectively separate different rates of molecular evolution in linked regions. The increased substitution rate in the lowland lineage can be largely explained by temperature-induced changes to the spontaneous mutation rate. A novel method based on entropy is developed to test for the existence of barrier loci using a minimal number of samples from the hybrid zone, a robust framework when system complexity far exceeds sample information. Overall, our results suggest that during the process of speciation, the separation of substitution rates can occur locally in the genome in parallel to the separation of gene pools.

Compromised base excision repair pathway in Mycobacterium tuberculosis imparts superior adaptability in the host

Tuberculosis caused by Mycobacterium tuberculosis (Mtb) is a significant public health concern, exacerbated by the emergence of drug-resistant TB. To combat the host’s dynamic environment, Mtb encodes multiple DNA repair enzymes that play a critical role in maintaining genomic integrity. Mtb possesses a GC-rich genome, rendering it highly susceptible to cytosine deaminations, resulting in the occurrence of uracils in the DNA. UDGs encoded by ung and udgB initiate the repair; hence we investigated the biological impact of deleting UDGs in the adaptation of pathogen. We generated gene replacement mutants of uracil DNA glycosylases, individually (RvΔung, RvΔudgB) or together (RvΔdKO). The double KO mutant, RvΔdKO exhibited remarkably higher spontaneous mutation rate, in the presence of antibiotics. Interestingly, RvΔdKO showed higher survival rates in guinea pigs and accumulated large number of SNPs as revealed by whole-genome sequence analysis. Competition assays revealed the superior fitness of RvΔdKO over Rv, both in ex vivo and in vivo conditions. We propose that compromised DNA repair results in the accumulation of mutations, and a subset of these drives adaptation in the host. Importantly, this property allowed us to utilize RvΔdKO for the facile identification of drug targets.

Functional Analysis of the NucS/EndoMS of the Hyperthermophilic Archaeon Sulfolobus islandicus REY15A

EndoMS is a recently identified mismatch specific endonuclease in Thermococcales of Archaea and Mycobacteria of Bacteria. The homologs of EndoMS are conserved in Archaea and Actinobacteria, where classic MutS-MutL-mediated DNA mismatch repair pathway is absent or non-functional. Here, we report a study on the in vitro mismatch cleavage activity and in vivo function of an EndoMS homolog (SisEndoMS) from Sulfolobus islandicus REY15A, the model archaeon belonging to Crenarchaeota. SisEndoMS is highly active on duplex DNA containing G/T, G/G, and T/T mismatches. Interestingly, the cleavage activity of SisEndoMS is stimulated by the heterotrimeric PCNAs, and when Mn2+ was used as the co-factor instead of Mg2+, SisEndoMS was also active on DNA substrates containing C/T or A/G mismatches, suggesting that the endonuclease activity can be regulated by ion co-factors and accessory proteins. We compared the spontaneous mutation rate of the wild type strain REY15A and ∆endoMS by counter selection against 5-fluoroorotic acid (5-FOA). The endoMS knockout mutant had much higher spontaneous mutation rate (5.06 × 10−3) than that of the wild type (4.6 × 10−6). A mutation accumulation analysis also showed that the deletion mutant had a higher mutation occurrence than the wild type, with transition mutation being the dominant, suggesting that SisEndoMS is responsible for mutation avoidance in this hyperthermophilic archaeon. Overexpression of the wild type SisEndoMS in S. islandicus resulted in retarded growth and abnormal cell morphology, similar to strains overexpressing Hje and Hjc, the Holliday junction endonucleases. Transcriptomic analysis revealed that SisEndoMS overexpression led to upregulation of distinct gene including the CRISPR-Cas IIIB system, methyltransferases, and glycosyltransferases, which are mainly localized to specific regions in the chromosome. Collectively, our results support that EndoMS proteins represent a noncanonical DNA repair pathway in Archaea. The mechanism of the mismatch repair pathway in Sulfolobus which have a single chromosome is discussed.

Participation of RecJ in the base excision repair pathway of Deinococcus radiodurans

RecJ reportedly participates in the base excision repair (BER) pathway, but structural and functional data are scarce. Herein, the Deinococcus radiodurans RecJ (drRecJ) deletion strain exhibited extreme sensitivity to hydrogen peroxide and methyl-methanesulphonate, as well as a high spontaneous mutation rate and an accumulation of unrepaired abasic sites in vivo, indicating the involvement of drRecJ in the BER pathway. The binding affinity and nuclease activity preference of drRecJ toward DNA substrates containing a 5′-P-dSpacer group, a 5′-deoxyribose-phosphate (dRP) mimic, were established. A 1.9 Å structure of drRecJ in complex with 5′-P-dSpacer-modified single-stranded DNA (ssDNA) revealed a 5′-monophosphate binding pocket and occupancy of 5′-dRP in the drRecJ nuclease core. The mechanism for RecJ 5′-dRP catalysis was explored using structural and biochemical data, and the results implied that drRecJ is not a canonical 5′-dRP lyase. Furthermore, in vitro reconstitution assays indicated that drRecJ tends to participate in the long-patch BER pathway rather than the short-patch BER pathway.

Evolution of Mutation Rate in Astronomically Large Phytoplankton Populations

Genetic diversity is expected to be proportional to population size, yet, there is a well-known, but unexplained lack of genetic diversity in large populations—the “Lewontin’s paradox.” Larger populations are expected to evolve lower mutation rates, which may help to explain this paradox. Here, we test this conjecture by measuring the spontaneous mutation rate in a ubiquitous unicellular marine phytoplankton species Emiliania huxleyi (Haptophyta) that has modest genetic diversity despite an astronomically large population size. Genome sequencing of E. huxleyi mutation accumulation lines revealed 455 mutations, with an unusual GC-biased mutation spectrum. This yielded an estimate of the per site mutation rate µ = 5.55×10−10 (CI 95%: 5.05×10−10 – 6.09×10−10), which corresponds to an effective population size Ne ∼ 2.7×106. Such a modest Ne is surprising for a ubiquitous and abundant species that accounts for up to 10% of global primary productivity in the oceans. Our results indicate that even exceptionally large populations do not evolve mutation rates lower than ∼10−10 per nucleotide per cell division. Consequently, the extreme disparity between modest genetic diversity and astronomically large population size in the plankton species cannot be explained by an unusually low mutation rate.

Low Spontaneous Mutation Rate and Pleistocene Radiation of Pea Aphids

Accurate estimates of divergence times are essential to understand the evolutionary history of species. It allows linking evolutionary histories of the diverging lineages with past geological, climatic, and other changes in environment and shed light on the processes involved in speciation. The pea aphid radiation includes multiple host races adapted to different legume host plants. It is thought that diversification in this system occurred very recently, over the past 8,000–16,000 years. This young age estimate was used to link diversification in pea aphids to the onset of human agriculture, and led to the establishment of the pea aphid radiation as a model system in the study of speciation with gene flow. Here, we re-examine the age of the pea aphid radiation, by combining a mutation accumulation experiment with a genome-wide estimate of divergence between distantly related pea aphid host races. We estimate the spontaneous mutation rate for pea aphids as 2.7×10-10 per haploid genome per parthenogenic generation. Using this estimate of mutation rate and the genome-wide genetic differentiation observed between pea aphid host races, we show that the pea aphid radiation is much more ancient than assumed previously, predating Neolithic agriculture by several hundreds of thousands of years. Our results rule out human agriculture as the driver of diversification of the pea aphid radiation, and call for re-assessment of the role of allopatric isolation during Pleistocene climatic oscillations in divergence of the pea aphid complex.

Mismatch repair systems might facilitate the chromosomal recombination induced by N-nitrosodimethylamine, but not by N-nitrosodiethylamine, in Drosophila

Mismatch repair (MMR) systems play important roles in maintaining the high fidelity of genomic DNA. It is well documented that a lack of MMR increases the mutation rate, including base exchanges and small insertion/deletion loops; however, it is unknown whether MMR deficiency affects the frequency of chromosomal recombination in somatic cells. To investigate the effects of MMR on chromosomal recombination, we used the Drosophila wing-spot test, which efficiently detects chromosomal recombination. We prepared MMR (MutS)-deficient flies (spel1(−/−)) using a fly line generated in this study. The spontaneous mutation rate as measured by the wing-spot test was slightly higher in MutS-deficient flies than in wild-type (spel1(+/−)) flies. Previously, we showed that N-nitrosodimethylamine (NDMA)-induced chromosomal recombination more frequently than N-nitrosodiethylamine (NDEA) in Drosophila. When the wing-spot test was performed using MMR-deficient flies, unexpectedly, the rate of NDMA-induced mutation was significantly lower in spel1(−/−) flies than in spel1(+/−) flies. In contrast, the rate of mutation induced by NDEA was higher in spel1(−/−) flies than in spel1(+/−) flies. These results suggest that in Drosophila, the MutS homologue protein recognises methylated DNA lesions more efficiently than ethylated ones, and that MMR might facilitate mutational chromosomal recombination due to DNA double-strand breaks via the futile cycle induced by MutS recognition of methylated lesions.

Low spontaneous mutation rate and Pleistocene radiation of pea aphids

Accurate estimates of divergence times are essential to understand the evolutionary history of species. It allows linking evolutionary histories of the diverging lineages with past geological, climatic and other changes in environment and shed light on the processes involved in speciation. The pea aphid radiation includes multiple host races adapted to different legume host plants. It is thought that diversification in this system occurred very recently, over the past 8,000 to 16,000 years. This young age estimate was used to link diversification in pea aphids to the onset of human agriculture, and lead to the establishment of the pea aphid radiation as a model system in the study of speciation with gene flow. Here, we re-examine the age of the pea aphid radiation, by combining a mutation accumulation experiment with a genome-wide estimate of divergence between distantly related pea aphid host races. We estimate the spontaneous mutation rate for pea aphids as 2.27 × 10−10 per haploid genome per parthenogenic generation. Using this estimate of mutation rate and the genome-wide genetic differentiation observed between pea aphid host races, we show that the pea aphid radiation is much more ancient than assumed previously, predating Neolithic agriculture by several hundreds of thousands of years. Our results rule out human agriculture as the driver of diversification of the pea aphid radiation, and call for re-assessment of the role of allopatric isolation during Pleistocene climatic oscillations in divergence of the pea aphid complex.

First Estimation of the Spontaneous Mutation Rate in Diatoms

Mutations are the origin of genetic diversity, and the mutation rate is a fundamental parameter to understand all aspects of molecular evolution. The combination of mutation–accumulation experiments and high-throughput sequencing enabled the estimation of mutation rates in most model organisms, but several major eukaryotic lineages remain unexplored. Here, we report the first estimation of the spontaneous mutation rate in a model unicellular eukaryote from the Stramenopile kingdom, the diatom Phaeodactylum tricornutum (strain RCC2967). We sequenced 36 mutation accumulation lines for an average of 181 generations per line and identified 156 de novo mutations. The base substitution mutation rate per site per generation is μbs = 4.77 × 10−10 and the insertion–deletion mutation rate is μid = 1.58 × 10−11. The mutation rate varies as a function of the nucleotide context and is biased toward an excess of mutations from GC to AT, consistent with previous observations in other species. Interestingly, the mutation rates between the genomes of organelles and the nucleus differ, with a significantly higher mutation rate in the mitochondria. This confirms previous claims based on indirect estimations of the mutation rate in mitochondria of photosynthetic eukaryotes that acquired their plastid through a secondary endosymbiosis. This novel estimate enables us to infer the effective population size of P. tricornutum to be Ne∼8.72 × 106.