Dana K. Howe

A genome-wide view of Caenorhabditis elegans base-substitution mutation processes

Knowledge of mutation processes is central to understanding virtually all evolutionary phenomena and the underlying nature of genetic disorders and cancers. However, the limitations of standard molecular mutation detection methods have historically precluded a genome-wide understanding of mutation rates and spectra in the nuclear genomes of multicellular organisms. We applied two high-throughput DNA sequencing technologies to identify and characterize hundreds of spontaneously arising base-substitution mutations in 10 Caenorhabditis elegans mutation-accumulation (MA)-line nuclear genomes. C. elegans mutation rate estimates were similar to previous calculations based on smaller numbers of mutations. Mutations were distributed uniformly within and among chromosomes and were not associated with recombination rate variation in the MA lines, suggesting that intragenomic variation in genetic hitchhiking and/or background selection are primarily responsible for the chromosomal distribution patterns of polymorphic nucleotides in C. elegans natural populations. A strong mutational bias from G/C to A/T nucleotides was detected in the MA lines, implicating oxidative DNA damage as a major endogenous mutagenic force in C. elegans. The observed mutational bias also suggests that the C. elegans nuclear genome cannot be at equilibrium because of mutation alone. Transversions dominate the spectrum of spontaneous mutations observed here, whereas transitions dominate patterns of allegedly neutral polymorphism in natural populations of C. elegans and many other animal species; this observation challenges the assumption that natural patterns of molecular variation in noncoding regions of the nuclear genome accurately reflect underlying mutation processes.

Variation in base-substitution mutation in experimental and natural lineages of Caenorhabditis nematodes

Variation among lineages in the mutation process has the potential to impact diverse biological processes ranging from susceptibilities to genetic disease to the mode and tempo of molecular evolution. The combination of high-throughput DNA sequencing (HTS) with mutation-accumulation (MA) experiments has provided a powerful approach to genome-wide mutation analysis, though insights into mutational variation have been limited by the vast evolutionary distances among the few species analyzed. We performed a HTS analysis of MA lines derived from four Caenorhabditis nematode natural genotypes: C. elegans N2 and PB306 and C. briggsae HK104 and PB800. Total mutation rates did not differ among the four sets of MA lines. A mutational bias toward G:C→A:T transitions and G:C→T:A transversions was observed in all four sets of MA lines. Chromosome-specific rates were mostly stable, though there was some evidence for a slightly elevated X chromosome mutation rate in PB306. Rates were homogeneous among functional coding sequence types and across autosomal cores, arms, and tips. Mutation spectra were similar among the four MA line sets but differed significantly when compared with patterns of natural base-substitution polymorphism for 13/14 comparisons performed. Our findings show that base-substitution mutation processes in these closely related animal lineages are mostly stable but differ from natural polymorphism patterns in these two species.

High Rate of Large Deletions in Caenorhabditis briggsae Mitochondrial Genome Mutation Processes

Mitochondrial DNA (mtDNA) mutations underlie a variety of human genetic disorders and are associated with the aging process. mtDNA polymorphisms are widely used in a variety of evolutionary applications. Although mtDNA mutation spectra are known to differ between distantly related model organisms, the extent to which mtDNA mutation processes vary between more closely related species and within species remains enigmatic. We analyzed mtDNA divergence in two sets of 250-generation Caenorhabditis briggsae mutation-accumulation (MA) lines, each derived from a different natural isolate progenitor: strain HK104 from Okayama, Japan, and strain PB800 from Ohio, United States. Both sets of C. briggsae MA lines accumulated numerous large heteroplasmic mtDNA deletions, whereas only one similar event was observed in a previous analysis of Caenorhabditis elegans MA line mtDNA. Homopolymer length change mutations were frequent in both sets of C. briggsae MA lines and occurred in both intergenic and protein-coding gene regions. The spectrum of C. briggsae mtDNA base substitution mutations differed from the spectrum previously observed in C. elegans. In C. briggsae, the HK104 MA lines experienced many different base substitution types, whereas the PB800 lines displayed only C:G → T:A transitions, although the difference was not significant. Over half of the mtDNA base substitutions detected in the C. briggsae MA lines were in a heteroplasmic state, whereas all those previously characterized in C. elegans MA line mtDNA were fixed changes, indicating a narrower mtDNA bottleneck in C. elegans as compared with C. briggsae. Our results show that C. briggsae mtDNA is highly susceptible to large deletions and that the mitochondrial mutation process varies between Caenorhabditis nematode species.

Selfish Little Circles: Transmission Bias and Evolution of Large Deletion-Bearing Mitochondrial DNA in Caenorhabditis briggsae Nematodes

Selfish DNA poses a significant challenge to genome stability and organismal fitness in diverse eukaryotic lineages. Although selfish mitochondrial DNA (mtDNA) has known associations with cytoplasmic male sterility in numerous gynodioecious plant species and is manifested as petite mutants in experimental yeast lab populations, examples of selfish mtDNA in animals are less common. We analyzed the inheritance and evolution of mitochondrial DNA bearing large heteroplasmic deletions including nad5 gene sequences (nad5Δ mtDNA), in the nematode Caenorhabditis briggsae. The deletion is widespread in C. briggsae natural populations and is associated with deleterious organismal effects. We studied the inheritance patterns of nad5Δ mtDNA using eight sets of C. briggsae mutation-accumulation (MA) lines, each initiated from a different natural strain progenitor and bottlenecked as single hermaphrodites across generations. We observed a consistent and strong drive toward higher levels of deletion-bearing molecules in the heteroplasmic pool of mtDNA after ten generations of bottlenecking. Our results demonstrate a uniform transmission bias whereby nad5Δ mtDNA accumulates to higher levels relative to intact mtDNA in multiple genetically diverse natural strains of C. briggsae. We calculated an average 1% per-generation transmission bias for deletion-bearing mtDNA relative to intact genomes. Our study, coupled with known deleterious phenotypes associated with high deletion levels, shows that nad5Δ mtDNA are selfish genetic elements that have evolved in natural populations of C. briggsae, offering a powerful new system to study selfish mtDNA dynamics in metazoans.

Selective sweeps and parallel mutation in the adaptive recovery from deleterious mutation in Caenorhabditis elegans.

Deleterious mutation poses a serious threat to human health and the persistence of small populations. Although adaptive recovery from deleterious mutation has been well-characterized in prokaryotes, the evolutionary mechanisms by which multicellular eukaryotes recover from deleterious mutation remain unknown. We applied high-throughput DNA sequencing to characterize genomic divergence patterns associated with the adaptive recovery from deleterious mutation using a Caenorhabditis elegans recovery-line system. The C. elegans recovery lines were initiated from a low-fitness mutation-accumulation (MA) line progenitor and allowed to independently evolve in large populations (N ∼ 1000) for 60 generations. All lines rapidly regained levels of fitness similar to the wild-type (N2) MA line progenitor. Although there was a near-zero probability of a single mutation fixing due to genetic drift during the recovery experiment, we observed 28 fixed mutations. Cross-generational analysis showed that all mutations went from undetectable population-level frequencies to a fixed state in 10-20 generations. Many recovery-line mutations fixed at identical timepoints, suggesting that the mutations, if not beneficial, hitchhiked to fixation during selective sweep events observed in the recovery lines. No MA line mutation reversions were detected. Parallel mutation fixation was observed for two sites in two independent recovery lines. Analysis using a C. elegans interactome map revealed many predicted interactions between genes with recovery line-specific mutations and genes with previously accumulated MA line mutations. Our study suggests that recovery-line mutations identified in both coding and noncoding genomic regions might have beneficial effects associated with compensatory epistatic interactions.

In Vivo Quantification Reveals Extensive Natural Variation in Mitochondrial Form and Function in Caenorhabditis briggsae

We have analyzed natural variation in mitochondrial form and function among a set of Caenorhabditis briggsae isolates known to harbor mitochondrial DNA structural variation in the form of a heteroplasmic nad5 gene deletion (nad5Δ) that correlates negatively with organismal fitness. We performed in vivo quantification of 24 mitochondrial phenotypes including reactive oxygen species level, membrane potential, and aspects of organelle morphology, and observed significant among-isolate variation in 18 traits. Although several mitochondrial phenotypes were non-linearly associated with nad5Δ levels, most of the among-isolate phenotypic variation could be accounted for by phylogeographic clade membership. In particular, isolate-specific mitochondrial membrane potential was an excellent predictor of clade membership. We interpret this result in light of recent evidence for local adaptation to temperature in C. briggsae. Analysis of mitochondrial-nuclear hybrid strains provided support for both mtDNA and nuclear genetic variation as drivers of natural mitochondrial phenotype variation. This study demonstrates that multicellular eukaryotic species are capable of extensive natural variation in organellar phenotypes and highlights the potential of integrating evolutionary and cell biology perspectives.

Using Mitogenomic and Nuclear Ribosomal Sequence Data to Investigate the Phylogeny of the Xiphinema americanum Species Complex

Nematodes within the Xiphinema americanum species complex are economically important because they vector nepoviruses which cause considerable damage to a variety of agricultural crops. The taxonomy of X. americanum species complex is controversial, with the number of putative species being the subject of debate. Accurate phylogenetic knowledge of this group is highly desirable as it may ultimately reveal genetic differences between species. For this study, nematodes belonging to the X. americanum species complex, including potentially mixed species populations, were collected from 12 geographically disparate locations across the U.S. from different crops and in varying association with nepoviruses. At least four individuals from each population were analyzed. A portion of the 18S nuclear ribosomal DNA (rDNA) gene was sequenced for all individuals while the internal transcribed spacer region 1 (ITS1) of rDNA was cloned and 2 to 6 clones per individual were sequenced. Mitochondrial genomes for numerous individuals were sequenced in parallel using high-throughput DNA sequencing (HTS) technology. Phylogenetic analysis of the 18S rDNA revealed virtually identical sequences across all populations. Analysis of ITS1 rDNA sequences revealed several well-supported clades, with some degree of congruence with geographic location and viral transmission, but also numerous presumably paralogous sequences that failed to form clades with other sequences from the same population. Analysis of mitochondrial DNA (mtDNA) indicated the presence of three distinct monophyletic clades of X. americanum species complex nematodes. Two clades contained nematodes found in association with nepovirus and the third contained divergent mtDNA sequences from three nematode populations from the western U.S. where nepovirus was absent. The inherent heterogeneity in ITS1 rDNA sequence data and lack of informative sites in 18S rDNA analysis suggests that mtDNA may be more useful in sorting out the taxonomic confusion of the X. americanum species complex.

Selfish Mitochondrial DNA Proliferates and Diversifies in Small, but not Large, Experimental Populations of Caenorhabditis briggsae

Evolutionary interactions across levels of biological organization contribute to a variety of fundamental processes including genome evolution, reproductive mode transitions, species diversification, and extinction. Evolutionary theory predicts that so-called “selfish” genetic elements will proliferate when the host effective population size (Ne) is small, but direct tests of this prediction remain few. We analyzed the evolutionary dynamics of deletion-containing mitochondrial DNA (ΔmtDNA) molecules, previously characterized as selfish elements, in six different natural strains of the nematode Caenorhabditis briggsae allowed to undergo experimental evolution in a range of population sizes (N = 1, 10, 100, and 1,000) for a maximum of 50 generations. Mitochondrial DNA (mtDNA) was analyzed for replicate lineages at each five-generation time point. Ten different ΔmtDNA molecule types were observed and characterized across generations in the experimental populations. Consistent with predictions from evolutionary theory, lab lines evolved in small-population sizes (e.g., nematode N = 1) were more susceptible to accumulation of high levels of preexisting ΔmtDNA compared with those evolved in larger populations. New ΔmtDNA elements were observed to increase in frequency and persist across time points, but almost exclusively at small population sizes. In some cases, ΔmtDNA levels decreased across generations when population size was large (nematode N = 1,000). Different natural strains of C. briggsae varied in their susceptibilities to ΔmtDNA accumulation, owing in part to preexisting compensatory mtDNA alleles in some strains that prevent deletion formation. This analysis directly demonstrates that the evolutionary trajectories of ΔmtDNA elements depend upon the population-genetic environments and molecular-genetic features of their hosts.

Comparative Genomics of a Plant-Parasitic Nematode Endosymbiont Suggest a Role in Nutritional Symbiosis

Bacterial mutualists can modulate the biochemical capacity of animals. Highly coevolved nutritional mutualists do this by synthesizing nutrients missing from the host’s diet. Genomics tools have advanced the study of these partnerships. Here we examined the endosymbiont Xiphinematobacter (phylum Verrucomicrobia) from the dagger nematode Xiphinema americanum, a migratory ectoparasite of numerous crops that also vectors nepovirus. Previously, this endosymbiont was identified in the gut, ovaries, and eggs, but its role was unknown. We explored the potential role of this symbiont using fluorescence in situ hybridization, genome sequencing, and comparative functional genomics. We report the first genome of an intracellular Verrucomicrobium and the first exclusively intracellular non-Wolbachia nematode symbiont. Results revealed that Xiphinematobacter had a small 0.916-Mb genome with only 817 predicted proteins, resembling genomes of other mutualist endosymbionts. Compared with free-living relatives, conserved proteins were shorter on average, and there was large-scale loss of regulatory pathways. Despite massive gene loss, more genes were retained for biosynthesis of amino acids predicted to be essential to the host. Gene ontology enrichment tests showed enrichment for biosynthesis of arginine, histidine, and aromatic amino acids, as well as thiamine and coenzyme A, diverging from the profiles of relatives Akkermansia muciniphilia (in the human colon), Methylacidiphilum infernorum, and the mutualist Wolbachia from filarial nematodes. Together, these features and the location in the gut suggest that Xiphinematobacter functions as a nutritional mutualist, supplementing essential nutrients that are depleted in the nematode diet. This pattern points to evolutionary convergence with endosymbionts found in sap-feeding insects.

Genome Skimming: A Rapid Approach to Gaining Diverse Biological Insights into Multicellular Pathogens

Genomic data acquisition is now trivial for biologists. Yet, moving from millions of sequence reads to an assembled and annotated genome continues to pose a daunting challenge. The first animal genome sequenced arose from the free-living model nematode Caenorhabditis elegans [1]. This venture provided an unprecedented foundation for new insights into genome function and ‘omics tool development. However, the C. elegans endeavor has been tough to repeat, even with the advent of new high-throughput DNA sequencing technologies. For example, the first plant-parasitic nematode (PPN) genomes were published ten years after the C. elegans genome [2,3], and only five publication-quality PPN genomes are presently available [4–6]. Fig 1 overviews the course of a typical genome project. Millions of DNA sequences are initially collected in a matter of days, thanks to new DNA sequencing technologies. Early analytical phases (quality control and initial assembly) are also quick and usually straightforward. However, the subsequent computational stages (refining the assembly, gene prediction, and annotation) present significant bioinformatics bottlenecks. These lengthy in silico steps require multiple iterative stages of analysis, finally leading to a finished genome deemed “good enough” for publication. These latter stages often take years. Fig 1 Genome skimming schematic. The term “genome skimming” was recently coined [7–9] to describe shallow sequencing approaches aiming to uncover conserved ortholog sequences for phylogenomic studies. Here, we overview a genome skimming strategy applied to six PPN species but expand the scope beyond phylogenetics and toward diverse questions relating to pathogen function and biology. We demonstrate our strategy’s utility in rapidly revealing insights and new hypotheses relating to nematode genome structure, effector genes, and endosymbionts.