Explore the vast range of titles available.

The mechanisms of sequence divergence in angiosperm mitochondrial genomes have long been enigmatic. In particular, it is difficult to reconcile the rapid divergence of intergenic regions that can make non-coding sequences almost unrecognizable even among close relatives with the unusually high levels of sequence conservation found in genic regions. It has been hypothesized that different mutation and repair mechanisms act on genic and intergenic sequences or alternatively that mutational input is relatively constant but that selection has strikingly different effects on these respective regions. To test these alternative possibilities, we analyzed mtDNA divergence within Arabidopsis thaliana, including variants from the 1001 Genomes Project and changes accrued in published mutation accumulation (MA) lines. We found that base-substitution frequencies are relatively similar for intergenic regions and synonymous sites in coding regions, whereas indel and nonsynonymous substitutions rates are greatly depressed in coding regions, supporting a conventional model in which mutation/repair mechanisms are consistent throughout the genome but differentially filtered by selection. Most types of sequence and structural changes were undetectable in 10-generation MA lines, but we found significant shifts in relative copy number across mtDNA regions for lines grown under stressed vs. benign conditions. We confirmed quantitative variation in copy number across the A. thaliana mitogenome using both whole-genome sequencing and droplet digital PCR, further undermining the classic but oversimplified model of a circular angiosperm mtDNA structure. Our results suggest that copy number variation is one of the most fluid features of angiosperm mitochondrial genomes.

Mitochondrial and plastid genomes in land plants exhibit some of the slowest rates of sequence evolution observed in any eukaryotic genome, suggesting an exceptional ability to prevent or correct mutations. However, the mechanisms responsible for this extreme fidelity remain unclear. We tested seven candidate genes involved in cytoplasmic DNA replication, recombination, and repair (POLIA, POLIB, MSH1, RECA3, UNG, FPG, and OGG1) for effects on mutation rates in the model angiosperm Arabidopsis thaliana by applying a highly accurate DNA sequencing technique (duplex sequencing) that can detect newly arisenmitochondrial and plastidmutations even at low heteroplasmic frequencies. We find that disrupting MSH1 (but not the other candidate genes) leads to massive increases in the frequency of point mutations and small indels and changes to the mutation spectrum in mitochondrial and plastid DNA. We also used droplet digital PCR to show transmission of de novo heteroplasmies across generations in msh1 mutants, confirming a contribution to heritable mutation rates. This dual-targeted gene is part of an enigmatic lineagewithin the mutS mismatch repair family thatwe find is also present outside of green plants in multiple eukaryotic groups (stramenopiles, alveolates, haptophytes, and cryptomonads), as well as certain bacteria and viruses. MSH1 has previously been shown to limit ectopic recombination in plant cytoplasmic genomes. Our results point to a broader role in recognition and correction of errors in plant mitochondrial and plastid DNA sequence, leading to greatly suppressed mutation rates perhaps via initiation of doublestranded breaks and repair pathways based on faithful homologous recombination.

Mutations are changes to DNA sequences which drive evolution by supplying raw genetic variation for natural selection to act upon. At the same time, mutations tend to have negative fitness consequences and are the source of genetic diseases. Such costs and benefits of mutation create opposing forces of selection on mutation rate modifiers, which are alleles (typically in DNA repair genes) responsible for increases or decreases to mutation rates. Essentially all eukaryotes possess at least two genomes: the nuclear genome (nucDNA) and the endosymbiotically derived mitochondrial genome (mtDNA). Photosynthetic plants and algae additionally possess endosymbiotically derived plastid genomes (cpDNAs). Together, the mtDNA and cpDNA are referred to as organellar genomes. Chapter 1 of this dissertation provides a framework for understanding how DNA repair machinery and mutation rates have evolved in complex eukaryotic cells. Chapters 2 and 3 focus on specific repair pathways active in organellar genomes. Finally, Chapter 4 shifts focus to understand how environmental perturbations in expression level impact mutation in plant nuclear genomes.Repair of organellar genomes is conducted by nuclear-encoded genes that are translated in the cytosol and targeted to the organelles. In terms of evolutionary history, organellar repair machinery is a mosaic network of bacterial-like repair genes, which came into the cell with the organelles, and nuclear repair genes that have been co-opted for organellar function. In some cases, repair proteins are targeted to both the nucDNA and mtDNA (and/or cpDNA in plants) to perform similar functions. This is the case as for many base excision repair (BER) proteins, which identify and remove of chemically damaged bases. In contrast, organellar repair arsenals are thought to lack canonical mismatch-repair (MMR) and nucleotide excision repair (NER), which are both important repair pathways in nuclear genomes. Instead, the diverse eukaryotic lineages have adopted unique strategies for organellar genome maintenance. As a result, there is a tremendous diversity in mtDNA mutation rates, which show over a 4000-fold variation across eukaryotes. Interestingly, much of this variation is driven by the extremely low point mutation rates plant mtDNAs. In addition, plant organellar genomes are more recombinationally active and plant mtDNAs are structurally unstable compared to the mtDNAs of other eukaryotes.Chapter 2 of this dissertation explores the mechanistic basis of low point mutation rates and recombination-mediated repair in plant organellar genomes. We performed high-fidelity Duplex Sequencing on a panel of Arabidopsis thaliana lines lacking specific organellar genome repair genes. We report large point mutation increases in mutant lines lacking MSH1, a mutS homolog that has been proposed to induce double-stranded breaks at the site of DNA mismatches, effectively shuttling such lesions into homologous recombination (HR) pathways that play important roles in plant organellar genome replication and repair. We see smaller point mutation increases in other mutant lines lacking RADA, RECA1 and RECA3. In addition, we generated long-read Oxford Nanopore sequencing to characterize repeat-mediated recombination in several of the mutants in the panel. Our findings provide valuable insights into the mechanisms driving the fascinating patterns of organellar genome evolution in land plants. The aforementioned lack of NER in organellar genomes is surprising given the importance of NER as a ‘catchall’ for repair of a variety of bulky DNA lesions in nuclear genomes. Chapter 3 focuses in on the fate of photodamage (an important type of bulky DNA damage) in organellar DNA. To do so, we leverage publicly available XR-seq datasets, which were generated to quantify and map active NER excision products in nuclear genomes following UV exposure. The taxonomic scope of chapter is expanded from plants to also include fungi (the brewer’s yeast Saccharomyces cerevisiae) and animals (the model fruit fly Drosophila melanogaster). We find that mtDNA-derived XR-seq reads in A. thaliana and S. cerevisiae have distinct and repeatable patterns in terms of length and internal positioning of pyrimidine dimers (known targets of photodamage formation). These data mirror established patterns of NER-derived reads originating from the nuclear genomes, raising the exciting possibility that NER-like repair pathways may exist in for repair of photodamage in organellar genomes. The focus of chapter 4 shifts to understanding how environmental changes impact mutation in plant nuclear genomes. The textbook view of mutation and adaptation is that mutations occur randomly with respect to their environment-specific fitness consequences. However, this view of random mutation is challenged as evidence increasingly establishes a correlation between increased expression and decreased mutation via the coordination of transcription and DNA repair machinery at the molecular level. As a result of this correlation, intragenomic mutation rates likely vary with changing environments given that expression levels are environmentally labile. Therefore, certain genes may be predisposed to higher or lower mutation rates depending on the environment, though the magnitude and importance of this effect remains largely untested. A technical challenge in addressing these questions is that large scale plant mutation datasets are time and resource intensive to generate. A recent plant study relied on low frequency somatic calls from Illumina based shotgun libraries to generate a large number of mutations, but others report that most of these inferred mutations are sequencing errors. To overcome these challenges, we took a novel approach to measuring somatic mutations by using Duplex Sequencing to quantify mutations in targeted regions of the A. thaliana nucDNA. We identified a set of differentially expressed genes in plants grown at different temperatures, which we then targeted for mutation detection using hybrid capture. In addition to wild type (WT) lines, we also studied mutant lines deficient in BER and MMR to test if either of these pathways are responsible for the correlation between expression and mutation in plants. We found large point mutation increases in the MMR mutants compared to WT plants, which displayed surprisingly few mutations at either temperature. Though the small number of WT mutations precluded a meaningful comparison of expression and mutation in a WT background, this result if nonetheless valuable for establishing that the true frequency of somatic mutations in plants is indeed very low suggesting that previous estimates likely conflated Illumina based sequencing artifacts with mutations. Mutation rates vary by over three orders of magnitude across the tree of life. Much of this variation is captured in mitochondrial mutation rates. The chapters of this dissertation provide valuable insights into the molecular processes that drive mutation rate variation in eukaryotic genomes. Such mechanistic understandings are critical for advancing the broader field of mutation rate evolution.

Land plant organellar genomes have extremely low rates of point mutation yet also experience high rates of recombination and genome instability. Characterizing the molecular machinery responsible for these patterns is critical for understanding the evolution of these genomes. While much progress has been made toward understanding recombination activity in land plant organellar genomes, the relationship between recombination pathways and point mutation rates remains uncertain. The organellar-targeted mutS homolog MSH1 has previously been shown to suppress point mutations as well as non-allelic recombination between short repeats in Arabidopsis thaliana. We therefore implemented high-fidelity Duplex Sequencing to test if other genes that function in recombination and maintenance of genome stability also affect point mutation rates. We found small to moderate increases in the frequency of single nucleotide variants (SNVs) and indels in mitochondrial and/or plastid genomes of A. thaliana mutant lines lacking radA, recA1, or recA3. In contrast, osb2 and why2 mutants did not exhibit an increase in point mutations compared to wild-type (WT) controls. In addition, we analyzed the distribution of SNVs in previously generated Duplex Sequencing data from A. thaliana organellar genomes and found unexpected strand asymmetries and large effects of flanking nucleotides on mutation rates in WT plants and msh1 mutants. Finally, using long-read Oxford Nanopore sequencing, we characterized structural variants in organellar genomes of the mutant lines and show that different short repeat sequences become recombinationally active in different mutant backgrounds. Together, these complementary sequencing approaches shed light on how recombination may impact the extraordinarily low point mutation rates in plant organellar genomes.

Compared to free-living bacteria, endosymbionts of sap-feeding insects have tiny and rapidly evolving genomes. Increased genetic drift, high mutation rates, and relaxed selection associated with host control of key cellular functions all likely contribute to genome decay. Phylogenetic comparisons have revealed massive variation in endosymbiont evolutionary rate, but such methods make it difficult to partition the effects of mutation vs. selection. For example, the ancestor of auchenorrhynchan insects contained two obligate endosymbionts, Sulcia and a betaproteobacterium (BetaSymb; called Nasuia in leafhoppers) that exhibit divergent rates of sequence evolution and different propensities for loss and replacement in the ensuing ~300 Ma. Here, we use the auchenorrhynchan leafhopper Macrosteles sp. nr. severini, which retains both of the ancestral endosymbionts, to test the hypothesis that differences in evolutionary rate are driven by differential mutagenesis. We used a high-fidelity technique known as duplex sequencing to measure and compare low-frequency variants in each endosymbiont. Our direct detection of de novo mutations reveals that the rapidly evolving endosymbiont (Nasuia) has a much higher frequency of single-nucleotide variants than the more stable endosymbiont (Sulcia) and a mutation spectrum that is even more AT-biased than implied by the 83.1% AT content of its genome. We show that indels are common in both endosymbionts but differ substantially in length and distribution around repetitive regions. Our results suggest that differences in long-term rates of sequence evolution in Sulcia vs. BetaSymb, and perhaps the contrasting degrees of stability of their relationships with the host, are driven by differences in mutagenesis. Competing Interest Statement The authors have declared no competing interest.

Mutations are the source of novel genetic diversity but can also lead to disease and maladaptation. The conventional view is that mutations occur randomly with respect to their environment-specific fitness consequences. However, intragenomic mutation rates can vary dramatically due to transcription coupled repair and based on local epigenomic modifications, which are non-uniformly distributed across genomes. One sequence feature associated with decreased mutation is higher expression level, which can vary depending on environmental cues. To understand whether the association between expression level and mutation rate creates a systematic relationship with environment-specific fitness effects, we perturbed expression through a heat treatment in . We quantified gene expression to identify differentially expressed genes, which we then targeted for mutation detection using Duplex Sequencing. This approach provided a highly accurate measurement of the frequency of rare somatic mutations in vegetative plant tissues, which has been a recent source of uncertainty in plant mutation research. We included mutant lines lacking mismatch repair (MMR) and base excision repair (BER) capabilities to understand how repair mechanisms may drive biased mutation accumulation. We found wild type (WT) and BER mutant mutation frequencies to be very low (mean variant frequency 1.8×10 and 2.6×10 , respectively), while MMR mutant frequencies were significantly elevated (1.13×10 ). These results show that somatic variant frequencies are extremely low in WT plants, indicating that larger datasets will be needed to address the fundamental evolutionary question as to whether environmental change leads to gene-specific changes in mutation rate.

Land plant organellar genomes have extremely low rates of point mutation yet also experience high rates of recombination and genome instability. Characterizing the molecular machinery responsible for these patterns is critical for understanding the evolution of these genomes. While much progress has been made towards understanding recombination activity in land plant organellar genomes, the relationship between recombination pathways and point mutation rates remains uncertain. The organellar targeted homolog MSH1 has previously been shown to suppress point mutations as well as non-allelic recombination between short repeats in . We therefore implemented high-fidelity Duplex Sequencing to test if other genes that function in recombination and maintenance of genome stability also affect point mutation rates. We found small to moderate increases in the frequency of single nucleotide variants (SNVs) and indels in mitochondrial and/or plastid genomes of mutant lines lacking , , or . In contrast, and mutants did not exhibit an increase in point mutations compared to wild type (WT) controls. In addition, we analyzed the distribution of SNVs in previously generated Duplex Sequencing data from organellar genomes and found unexpected strand asymmetries and large effects of flanking nucleotides on mutation rates in WT plants and mutants. Finally, using long-read Oxford Nanopore sequencing, we characterized structural variants in organellar genomes of the mutant lines and show that different short repeat sequences become recombinationally active in different mutant backgrounds. Together, these complementary sequencing approaches shed light on how recombination may impact the extraordinarily low point mutation rates in plant organellar genomes.

UV light is a potent mutagen that induces bulky DNA damage in the form of cyclobutane pyrimidine dimers (CPDs). In eukaryotic cells, photodamage and other bulky lesions occurring in nuclear genomes (nucDNAs) can be repaired through nucleotide excision repair (NER), where dual incisions on both sides of a damaged site precede the removal of a single-stranded oligonucleotide containing the damage. Mitochondrial genomes (mtDNAs) are also susceptible to damage from UV light, but current views hold that the only way to eliminate bulky DNA damage in mtDNAs is through mtDNA degradation. Damage-containing oligonucleotides excised during NER can be captured with anti-damage antibodies and sequenced (XR-seq) to produce high resolution maps of active repair locations following UV exposure. We analyzed previously published datasets from , and to identify reads originating from the mtDNA (and plastid genome in ). In and , the mtDNA-mapping reads have unique length distributions compared to the nuclear-mapping reads. The dominant fragment size was 26 nt in and 28 nt in with distinct secondary peaks occurring in 2-nt ( ) or 4-nt ( ) intervals. These reads also show a nonrandom distribution of di-pyrimidines (the substrate for CPD formation) with TT enrichment at positions 7-8 of the reads. Therefore, UV damage to mtDNA appears to result in production of DNA fragments of characteristic lengths and positions relative to the damaged location. We hypothesize that these fragments may reflect the outcome of a previously uncharacterized mechanism of NER-like repair in mitochondria or a programmed mtDNA degradation pathway.

The angiosperm genus Silene has been the subject of extensive study in the field of ecology and evolution, but the availability of high-quality reference genome sequences has been limited for this group. Here, we report a chromosome-level assembly for the genome of Silene conica based on PacBio HiFi, Hi-C and Bionano technologies. The assembly produced 10 scaffolds (one per chromosome) with a total length of 862 Mb and only ~1% gap content. These results confirm previous observations that S. conica and its relatives have a reduced base chromosome number relative to the genus's ancestral state of 12. Silene conica has an exceptionally large mitochondrial genome (>11 Mb), predominantly consisting of sequence of unknown origins. Analysis of shared sequence content suggests that it is unlikely that transfer of nuclear DNA is the primary driver of this mitochondrial genome expansion. More generally, this assembly should provide a valuable resource for future genomic studies in Silene, including comparative analyses with related species that recently evolved sex chromosomes.The angiosperm genus Silene has been the subject of extensive study in the field of ecology and evolution, but the availability of high-quality reference genome sequences has been limited for this group. Here, we report a chromosome-level assembly for the genome of Silene conica based on PacBio HiFi, Hi-C and Bionano technologies. The assembly produced 10 scaffolds (one per chromosome) with a total length of 862 Mb and only ~1% gap content. These results confirm previous observations that S. conica and its relatives have a reduced base chromosome number relative to the genus's ancestral state of 12. Silene conica has an exceptionally large mitochondrial genome (>11 Mb), predominantly consisting of sequence of unknown origins. Analysis of shared sequence content suggests that it is unlikely that transfer of nuclear DNA is the primary driver of this mitochondrial genome expansion. More generally, this assembly should provide a valuable resource for future genomic studies in Silene, including comparative analyses with related species that recently evolved sex chromosomes.Whole-genome sequences have been largely lacking for species in the genus Silene even though these flowering plants have been used for studying ecology, evolution, and genetics for over a century. Here, we address this gap by providing a high-quality nuclear genome assembly for S. conica, a species known to have greatly accelerated rates of sequence and structural divergence in its mitochondrial and plastid genomes. This resource will be valuable in understanding the coevolutionary interactions between nuclear and cytoplasmic genomes and in comparative analyses across this highly diverse genus.SignificanceWhole-genome sequences have been largely lacking for species in the genus Silene even though these flowering plants have been used for studying ecology, evolution, and genetics for over a century. Here, we address this gap by providing a high-quality nuclear genome assembly for S. conica, a species known to have greatly accelerated rates of sequence and structural divergence in its mitochondrial and plastid genomes. This resource will be valuable in understanding the coevolutionary interactions between nuclear and cytoplasmic genomes and in comparative analyses across this highly diverse genus.

The mechanisms of sequence divergence in angiosperm mitochondrial genomes have long been enigmatic. In particular, it is difficult to reconcile the rapid divergence of intergenic regions that can make non-coding sequences almost unrecognizable even among close relatives with the unusually high levels of sequence conservation found in genic regions. It has been hypothesized that different mutation/repair mechanisms act on genic and intergenic sequences or alternatively that mutational input is relatively constant but that selection has strikingly different effects on these respective regions. To test these alternative possibilities, we analyzed mtDNA divergence within Arabidopsis thaliana, including variants from the 1001 Genomes Project and changes accrued in published mutation accumulation (MA) lines. We found that base-substitution frequencies are relatively similar for intergenic regions and synonymous sites in coding regions, whereas indel and nonsynonymous substitutions rates are greatly depressed in coding regions, supporting a conventional model in which mutation/repair mechanisms are consistent throughout the genome but differentially filtered by selection. Most types of sequence and structural changes were undetectable in 10-generation MA lines, but we found significant shifts in relative copy number across mtDNA regions for lines grown under stressed vs. benign conditions. We confirmed quantitative variation in copy number across the A. thaliana mitogenome using both whole-genome sequencing and droplet digital PCR, further undermining the classic but oversimplified model of a circular angiosperm mtDNA structure. Our results suggest that copy number variation is one of the most rapidly evolving features in angiosperm mtDNA, even outpacing rearrangements in these notoriously structurally diverse genomes.