Genomic Analyses of Globodera pallida, A Quarantine Agricultural Pathogen in Idaho.

Globodera pallida is among the most significant plant-parasitic nematodes worldwide, causing major damage to potato production. Since it was discovered in Idaho in 2006, eradication efforts have aimed to contain and eradicate G. pallida through phytosanitary action and soil fumigation. In this study, we investigated genome-wide patterns of G. pallida genetic variation across Idaho fields to evaluate whether the infestation resulted from a single or multiple introduction(s) and to investigate potential evolutionary responses since the time of infestation. A total of 53 G. pallida samples (~1,042,000 individuals) were collected and analyzed, representing five different fields in Idaho, a greenhouse population, and a field in Scotland that was used for external comparison. According to genome-wide allele frequency and fixation index (Fst) analyses, most of the genetic variation was shared among the G. pallida populations in Idaho fields pre-fumigation, indicating that the infestation likely resulted from a single introduction. Temporal patterns of genome-wide polymorphisms involving (1) pre-fumigation field samples collected in 2007 and 2014 and (2) pre- and post-fumigation samples revealed nucleotide variants (SNPs, single-nucleotide polymorphisms) with significantly differentiated allele frequencies indicating genetic differentiation. This study provides insights into the genetic origins and adaptive potential of G. pallida invading new environments.

Phylogenetic evidence for the invasion of a commercialized European Phasmarhabditis hermaphrodita lineage into North America and New Zealand.

Biological control (biocontrol) as a component of pest management strategies reduces reliance on synthetic chemicals, and seemingly offers a natural approach that minimizes environmental impact. However, introducing a new organism to new environments as a classical biocontrol agent can have broad and unanticipated biodiversity effects and conservation consequences. Nematodes are currently used in a variety of commercial biocontrol applications, including the use of Phasmarhabditis hermaphrodita as an agent targeting pest slug and snail species. This species was originally discovered in Germany, and is generally thought to have European origins. P. hermaphrodita is sold under the trade name Nemaslug®, and is available only in European markets. However, this nematode species was discovered in New Zealand and the western United States, though its specific origins remained unclear. In this study, we analyzed 45 nematode strains representing eight different Phasmarhabditis species, collected from nine countries around the world. A segment of nematode mitochondrial DNA (mtDNA) was sequenced and subjected to phylogenetic analyses. Our mtDNA phylogenies were overall consistent with previous analyses based on nuclear ribosomal RNA (rRNA) loci. The recently discovered P. hermaphrodita strains in New Zealand and the United States had mtDNA haplotypes nearly identical to that of Nemaslug®, and these were placed together in an intraspecific monophyletic clade with high support in maximum likelihood and Bayesian analyses. We also examined bacteria that co-cultured with the nematode strains isolated in Oregon, USA, by analyzing 16S rRNA sequences. Eight different bacterial genera were found to associate with these nematodes, though Moraxella osloensis, the bacteria species used in the Nemaslug® formulation, was not detected. This study provided evidence that nematodes deriving from the Nemaslug® biocontrol product have invaded countries where its use is prohibited by regulatory agencies and not commercially available.

First report of molecular characterization and phylogeny of Trichuris fossor Hall, 1916 (Nematoda: Trichuridae).

Because species of Trichuris are morphologically similar and ranges of host preference are variable, using molecular data to evaluate species delineations is essential for properly quantifying biodiversity of and relationships within Trichuridae. Trichuris fossor has been reported from Thomomys spp. (Rodentia: Geomyidae, 'pocket gophers') hosts based on morphological features alone. Partial 18S rRNA sequences for specimens identified as T. fossor based on morphology, along with sequences from 26 additional taxa, were used for a phylogenetic analysis. Evolutionary histories were constructed using maximum likelihood and Bayesian inference. In both analyses, the specimens fell within the Trichuris clade with 100% support and formed a distinct subclade with 100% support. These results confirm that T. fossor is a distinct species and represent the first molecular report for it. Relatedness among species within the family were well resolved in the BI tree. This study represents an initial effort to obtain a more comprehensive view of Trichuridae by including a new clade member, T. fossor. A better understanding of Trichuridae phylogeny could contribute to further characterization of host-associations, including species that infect livestock and humans.Because species of Trichuris are morphologically similar and ranges of host preference are variable, using molecular data to evaluate species delineations is essential for properly quantifying biodiversity of and relationships within Trichuridae. Trichuris fossor has been reported from Thomomys spp. (Rodentia: Geomyidae, 'pocket gophers') hosts based on morphological features alone. Partial 18S rRNA sequences for specimens identified as T. fossor based on morphology, along with sequences from 26 additional taxa, were used for a phylogenetic analysis. Evolutionary histories were constructed using maximum likelihood and Bayesian inference. In both analyses, the specimens fell within the Trichuris clade with 100% support and formed a distinct subclade with 100% support. These results confirm that T. fossor is a distinct species and represent the first molecular report for it. Relatedness among species within the family were well resolved in the BI tree. This study represents an initial effort to obtain a more comprehensive view of Trichuridae by including a new clade member, T. fossor. A better understanding of Trichuridae phylogeny could contribute to further characterization of host-associations, including species that infect livestock and humans.

Mitochondrial DNA Variation and Selfish Propagation Following Experimental Bottlenecking in Two Distantly Related Caenorhabditis briggsae Isolates.

Understanding mitochondrial DNA (mtDNA) evolution and inheritance has broad implications for animal speciation and human disease models. However, few natural models exist that can simultaneously represent mtDNA transmission bias, mutation, and copy number variation. Certain isolates of the nematode Caenorhabditis briggsae harbor large, naturally-occurring mtDNA deletions of several hundred basepairs affecting the NADH dehydrogenase subunit 5 (nduo-5) gene that can be functionally detrimental. These deletion variants can behave as selfish DNA elements under genetic drift conditions, but whether all of these large deletion variants are transmitted in the same preferential manner remains unclear. In addition, the degree to which transgenerational mtDNA evolution profiles are shared between isolates that differ in their propensity to accumulate the nduo-5 deletion is also unclear. We address these knowledge gaps by experimentally bottlenecking two isolates of C. briggsae with different nduo-5 deletion frequencies for up to 50 generations and performing total DNA sequencing to identify mtDNA variation. We observed multiple mutation profile differences and similarities between C. briggsae isolates, a potentially species-specific pattern of copy number dysregulation, and some evidence for genetic hitchhiking in the deletion-bearing isolate. Our results further support C. briggsae as a practical model for characterizing naturally-occurring mtgenome variation and contribute to the understanding of how mtgenome variation persists in animal populations and how it presents in mitochondrial disease states.

Complex Transmission Patterns and Age-Related Dynamics of a Selfish mtDNA Deletion.

Despite wide-ranging implications of selfish mitochondrial DNA (mtDNA) elements for human disease and topics in evolutionary biology (e.g., speciation), the forces controlling their formation, age-related accumulation, and offspring transmission remain largely unknown. Selfish mtDNA poses a significant challenge to genome integrity, mitochondrial function, and organismal fitness. For instance, numerous human diseases are associated with mtDNA mutations; however, few genetic systems can simultaneously represent pathogenic mitochondrial genome evolution and inheritance. The nematode Caenorhabditis briggsae is one such system. Natural C. briggsae isolates harbor varying levels of a large-scale deletion affecting the mitochondrial nduo-5 gene, termed nad5Δ. A subset of these isolates contains putative compensatory mutations that may reduce the risk of deletion formation. We studied the dynamics of nad5Δ heteroplasmy levels during animal development and transmission from mothers to offspring in genetically diverse C. briggsae natural isolates. Results support previous work demonstrating that nad5Δ is a selfish element and that heteroplasmy levels of this deletion can be quite plastic, exhibiting high degrees of inter-family variability and divergence between generations. The latter is consistent with a mitochondrial bottleneck effect, and contrasts with previous findings from a laboratory-derived model uaDf5 mtDNA deletion in C. elegans. However, we also found evidence for among-isolate differences in the ability to limit nad5Δ accumulation, the pattern of which suggested that forces other than the compensatory mutations are important in protecting individuals and populations from rampant mtDNA deletion expansion over short time scales.

Variable Abundance and Distribution of Wolbachia and Cardinium Endosymbionts in Plant-Parasitic Nematode Field Populations.

The bacterial endosymbiont Wolbachia interacts with different invertebrate hosts, engaging in diverse symbiotic relationships. Wolbachia is often a reproductive parasite in arthropods, but an obligate mutualist in filarial nematodes. Wolbachia was recently discovered in plant-parasitic nematodes, and, is thus far known in just two genera Pratylenchus and Radopholus, yet the symbiont's function remains unknown. The occurrence of Wolbachia in these economically important plant pests offers an unexplored biocontrol strategy. However, development of Wolbachia-based biocontrol requires an improved understanding of symbiont-host functional interactions and the symbiont's prevalence among nematode field populations. This study used a molecular-genetic approach to assess the prevalence of a Wolbachia lineage (wPpe) in 32 field populations of Pratylenchus penetrans. Populations were examined from eight different plant species in Washington, Oregon, and California. Nematodes were also screened for the endosymbiotic bacterium Cardinium (cPpe) that was recently shown to co-infect P. penetrans. Results identified wPpe in 9/32 and cPpe in 1/32 of P. penetrans field populations analyzed. No co-infection was observed in field populations. Wolbachia was detected in nematodes from 4/8 plant-hosts examined (raspberry, strawberry, clover, and lily), and in all three states surveyed. Cardinium was detected in nematodes from mint in Washington. In the wPpe-infected P. penetrans populations collected from raspberry, the prevalence of wPpe infection ranged from 11 to 58%. This pattern is unlike that in filarial nematodes where Wolbachia is an obligate mutualist and occurs in 100% of the host. Further analysis of wPpe-infected populations revealed female-skewed sex ratios (up to 96%), with the degree of skew positively correlating with wPpe prevalence. Uninfected nematode populations had approximately equal numbers of males and females. Comparisons of 54 wPpe 16S ribosomal RNA sequences revealed high similarity across the geographic isolates, with 45 of 54 isolates being identical at this locus. The complete absence of wPpe among some populations and low prevalence in others suggest that this endosymbiont is not an obligate mutualist of P. penetrans. The observed sex ratio bias in wPpe-infected nematode populations is similar to that observed in arthropods where Wolbachia acts as a reproductive manipulator, raising the question of a similar role in plant-parasitic nematodes.

Sex and Mitonuclear Adaptation in Experimental Caenorhabditis elegans Populations.

To reveal phenotypic and functional genomic patterns of mitonuclear adaptation, a laboratory adaptation study with Caenorhabditis elegans nematodes was conducted in which independently evolving lines were initiated from a low-fitness mitochondrial electron transport chain (ETC) mutant, gas-1 Following 60 generations of evolution in large population sizes with competition for food resources, two distinct classes of lines representing different degrees of adaptive response emerged: a low-fitness class that exhibited minimal or no improvement compared to the gas-1 mutant ancestor, and a high-fitness class containing lines that exhibited partial recovery of wild-type fitness. Many lines that achieved higher reproductive and competitive fitness levels were also noted to evolve high frequencies of males during the experiment, consistent with adaptation in these lines having been facilitated by outcrossing. Whole-genome sequencing and analysis revealed an enrichment of mutations in loci that occur in a gas-1-centric region of the C. elegans interactome and could be classified into a small number of functional genomic categories. A highly nonrandom pattern of mitochondrial DNA mutation was observed within high-fitness gas-1 lines, with parallel fixations of nonsynonymous base substitutions within genes encoding NADH dehydrogenase subunits I and VI. These mitochondrial gene products reside within ETC complex I alongside the nuclear-encoded GAS-1 protein, suggesting that rapid adaptation of select gas-1 recovery lines was driven by fixation of compensatory mitochondrial mutations.

Genome Announcement: The Draft Genomes of Two Radopholus similis populations from Costa Rica.

Radopholus similis is an economically important pest of both banana and citrus in tropical regions. Here we present draft genomes from two populations of R. similis from Costa Rica that were created and assembled using short read libraries from Illumina HiSeq technology.

Comparative Genomics of Wolbachia-Cardinium Dual Endosymbiosis in a Plant-Parasitic Nematode.

Wolbachia and Cardinium are among the most important and widespread of all endosymbionts, occurring in nematodes and more than half of insect and arachnid species, sometimes as coinfections. These symbionts are of significant interest as potential biocontrol agents due to their abilities to cause major effects on host biology and reproduction through cytoplasmic incompatibility, sex ratio distortion, or obligate mutualism. The ecological and metabolic effects of coinfections are not well understood. This study examined a Wolbachia-Cardinium coinfection in the plant-parasitic nematode (PPN), Pratylenchus penetrans, producing the first detailed study of such a coinfection using fluorescence in situ hybridization (FISH), polymerase chain reaction (PCR), and comparative genomic analysis. Results from FISH and single-nematode PCR showed 123/127 individuals in a focal population carried Cardinium (denoted strain cPpe), and 48% were coinfected with Wolbachia strain wPpe. Both endosymbionts showed dispersed tissue distribution with highest densities in the anterior intestinal walls and gonads. Phylogenomic analyses confirmed an early place of cPpe and long distance from a sister strain in another PPN, Heterodera glycines, supporting a long history of both Cardinium and Wolbachia in PPNs. The genome of cPpe was 1.36 Mbp with 35.8% GC content, 1,131 predicted genes, 41% having no known function, and missing biotin and lipoate synthetic capacity and a plasmid present in other strains, despite having a slightly larger genome compared to other sequenced Cardinium. The larger genome revealed expansions of gene families likely involved in host-cellular interactions. More than 2% of the genes of cPpe and wPpe were identified as candidate horizontally transferred genes, with some of these from eukaryotes, including nematodes. A model of the possible Wolbachia-Cardinium interaction is proposed with possible complementation in function for pathways such as methionine and fatty acid biosynthesis and biotin transport.

Comparative Genomic Analysis of 130 Bacteriophages Infecting Bacteria in the Genus Pseudomonas.

Bacteria of the genus Pseudomonas are genetically diverse and ubiquitous in the environment. Like other bacteria, those of the genus Pseudomonas are susceptible to bacteriophages which can significantly affect their host in many ways, ranging from cell lysis to major changes in morphology and virulence. Insights into phage genomes, evolution, and functional relationships with their hosts have the potential to contribute to a broader understanding of Pseudomonas biology, and the development of novel phage therapy strategies. Here we provide a broad-based comparative and evolutionary analysis of 130 complete Pseudomonas phage genome sequences available in online databases. We discovered extensive variation in genome size (ranging from 3 to 316 kb), G + C percentage (ranging from 37 to 66%), and overall gene content (ranging from 81-96% of genome space). Based on overall nucleotide similarity and the numbers of shared gene products, 100 out of 130 genome sequences were grouped into 12 different clusters; 30 were characterized as singletons, which do not have close relationships with other phage genomes. For 5/12 clusters, constituent phage members originated from two or more different Pseudomonas host species, suggesting that phage in these clusters can traverse bacterial species boundaries. An analysis of CRISPR spacers in Pseudomonas bacterial genome sequences supported this finding. Substantial diversity was revealed in analyses of phage gene families; out of 4,462 total families, the largest had only 39 members and there were 2,992 families with only one member. An evolutionary analysis of 72 phage gene families, based on patterns of nucleotide diversity at non-synonymous and synonymous sites, revealed strong and consistent signals for purifying selection. Our study revealed highly diverse and dynamic Pseudomonas phage genomes, and evidence for a dominant role of purifying selection in shaping the evolution of genes encoded in them.

Molecular signatures of host specificity linked to habitat specialization in Exaiptasia sea anemones.

Rising ocean temperatures associated with global climate change induce breakdown of the symbiosis between coelenterates and photosynthetic microalgae of the genus Symbiodinium. Association with more thermotolerant partners could contribute to resilience, but the genetic mechanisms controlling specificity of hosts for particular Symbiodinium types are poorly known. Here, we characterize wild populations of a sea anemone laboratory model system for anthozoan symbiosis, from contrasting environments in Caribbean Panama. Patterns of anemone abundance and symbiont diversity were consistent with specialization of holobionts for particular habitats, with Exaiptasia pallida/S. minutum (ITS2 type B1) abundant on vertical substrate in thermally stable, shaded environments but E. brasiliensis/Symbiodinium sp. (ITS2 clade A) more common in shallow areas subject to high temperature and irradiance. Population genomic sequencing revealed a novel E. pallida population from the Bocas del Toro Archipelago that only harbors S. minutum. Loci most strongly associated with divergence of the Bocas-specific population were enriched in genes with putative roles in cnidarian symbiosis, including activators of the complement pathway of the innate immune system, thrombospondin-type-1 repeat domain proteins, and coordinators of endocytic recycling. Our findings underscore the importance of unmasking cryptic diversity in natural populations and the role of long-term evolutionary history in mediating interactions with Symbiodinium.

Adaptive Evolution under Extreme Genetic Drift in Oxidatively Stressed Caenorhabditis elegans.

A mutation-accumulation (MA) experiment with Caenorhabditis elegans nematodes was conducted in which replicate, independently evolving lines were initiated from a low-fitness mitochondrial electron transport chain mutant, gas-1. The original intent of the study was to assess the effect of electron transport chain dysfunction involving elevated reactive oxygen species production on patterns of spontaneous germline mutation. In contrast to results of standard MA experiments, gas-1 MA lines evolved slightly higher mean fitness alongside reduced among-line genetic variance compared with their ancestor. Likewise, the gas-1 MA lines experienced partial recovery to wildtype reactive oxygen species levels. Whole-genome sequencing and analysis revealed that the molecular spectrum but not the overall rate of nuclear DNA mutation differed from wildtype patterns. Further analysis revealed an enrichment of mutations in loci that occur in a gas-1-centric region of the C. elegans interactome, and could be classified into a small number of functional-genomic categories. Characterization of a backcrossed four-mutation set isolated from one gas-1 MA line revealed this combination to be beneficial on both gas-1 mutant and wildtype genetic backgrounds. Our combined results suggest that selection favoring beneficial mutations can be powerful even under unfavorable population genetic conditions, and agree with fitness landscape theory predicting an inverse relationship between population fitness and the likelihood of adaptation.

The Draft Genome of Globodera ellingtonae.

Globodera ellingtonae is a newly described potato cyst nematode (PCN) found in Idaho, Oregon, and Argentina. Here, we present a genome assembly for G. ellingtonae, a relative of the quarantine nematodes G. pallida and G. rostochiensis, produced using data from Illumina and Pacific Biosciences DNA sequencing technologies.

Transcription through the eye of a needle: daily and annual cyclic gene expression variation in Douglas-fir needles.

BACKGROUND: Perennial growth in plants is the product of interdependent cycles of daily and annual stimuli that induce cycles of growth and dormancy. In conifers, needles are the key perennial organ that integrates daily and seasonal signals from light, temperature, and water availability. To understand the relationship between seasonal cycles and seasonal gene expression responses in conifers, we examined diurnal and circannual needle mRNA accumulation in Douglas-fir (Pseudotsuga menziesii) needles at diurnal and circannual scales. Using mRNA sequencing, we sampled 6.1 × 109 reads from 19 trees and constructed a de novo pan-transcriptome reference that includes 173,882 tree-derived transcripts. Using this reference, we mapped RNA-Seq reads from 179 samples that capture daily and annual variation. RESULTS: We identified 12,042 diurnally-cyclic transcripts, 9299 of which showed homology to annotated genes from other plant genomes, including angiosperm core clock genes. Annual analysis revealed 21,225 circannual transcripts, 17,335 of which showed homology to annotated genes from other plant genomes. The timing of maximum gene expression is associated with light intensity at diurnal scales and photoperiod at annual scales, with approximately half of transcripts reaching maximum expression +/- 2 h from sunrise and sunset, and +/- 20 days from winter and summer solstices. Comparisons with published studies from other conifers shows congruent behavior in clock genes with Japanese cedar (Cryptomeria), and a significant preservation of gene expression patterns for 2278 putative orthologs from Douglas-fir during the summer growing season, and 760 putative orthologs from spruce (Picea) during the transition from fall to winter. CONCLUSIONS: Our study highlight the extensive diurnal and circannual transcriptome variability demonstrated in conifer needles. At these temporal scales, 29% of expressed transcripts show a significant diurnal cycle, and 58.7% show a significant circannual cycle. Remarkably, thousands of genes reach their annual peak activity during winter dormancy. Our study establishes the fine-scale timing of daily and annual maximum gene expression for diverse needle genes in Douglas-fir, and it highlights the potential for using this information for evaluating hypotheses concerning the daily or seasonal timing of gene activity in temperate-zone conifers, and for identifying cyclic transcriptome components in other conifer species.

Natural Variation in Responses to Acute Heat and Cold Stress in a Sea Anemone Model System for Coral Bleaching.

Rising ocean temperatures disrupt the symbiosis between corals and their microalgae, accelerating global decline of coral reef ecosystems. Because of the difficulty of performing laboratory experiments with corals, the sea anemone Aiptasia has emerged as an important model system for molecular studies of coral bleaching and symbiosis. Here, we investigate natural variation in bleaching responses among different genetic lineages of Aiptasia. Both heat- and cold-induced paths to symbiosis breakdown were analyzed. Significant genetic variation in response to acute heat stress was observed, with severe bleaching of two Aiptasia strains from Hawaii but minimal bleaching of strains from the U.S. South Atlantic, including the strain used to generate the Aiptasia reference genome. Both strains from Hawaii hosted Symbiodinium type B1, whereas strains from the U.S. South Atlantic hosted type A4 or B2. In contrast to the results from exposures to acute heat stress, negligible variation was observed in response to a pulsed cold shock despite moderate bleaching across all strains. These results support our hypothesis that bleaching responses to distinct stressors are independent. Our findings emphasize the role of stress regime when predicting adaptive responses of symbiotic cnidarians to changing climates, because genetic variation may exist for some forms of stress-induced bleaching but not others.

Genomic evidence for plant-parasitic nematodes as the earliest Wolbachia hosts.

Wolbachia, one of the most widespread endosymbionts, is a target for biological control of mosquito-borne diseases (malaria and dengue virus), and antibiotic elimination of infectious filarial nematodes. We sequenced and analyzed the genome of a new Wolbachia strain (wPpe) in the plant-parasitic nematode Pratylenchus penetrans. Phylogenomic analyses placed wPpe as the earliest diverging Wolbachia, suggesting two evolutionary invasions into nematodes. The next branches comprised strains in sap-feeding insects, suggesting Wolbachia may have first evolved as a nutritional mutualist. Genome size, protein content, %GC, and repetitive DNA allied wPpe with mutualistic Wolbachia, whereas gene repertoire analyses placed it between parasite (A, B) and mutualist (C, D, F) groups. Conservation of iron metabolism genes across Wolbachia suggests iron homeostasis as a potential factor in its success. This study enhances our understanding of this globally pandemic endosymbiont, highlighting genetic patterns associated with host changes. Combined with future work on this strain, these genomic data could help provide potential new targets for plant-parasitic nematode control.

Paternal Mitochondrial Transmission in Intra-Species Caenorhabditis briggsae Hybrids.

To study mitochondrial-nuclear genetic interactions in the nematode Caenorhabditis briggsae, our three laboratories independently created 38 intra-species cytoplasmic-nuclear hybrid (cybrid) lines. Although the cross design combines maternal mitotypes with paternal nuclear genotypes, eight lines (21%) unexpectedly contained paternal mitotypes. All eight share in common ancestry of one of two genetically related strains. This unexpected parallel observation of paternal mitochondrial transmission, undesirable given our intent of creating cybrids, provides a serendipitous experimental model and framework to study the molecular and evolutionary basis of uniparental mitochondrial inheritance.

The mitochondrial genome of Globodera ellingtonae is composed of two circles with segregated gene content and differential copy numbers.

BACKGROUND: The evolution of animal mitochondrial (mt) genomes has resulted in a highly conserved structure: a single compact circular chromosome approximately 14 to 20 kb long. Within the last two decades exceptions to this conserved structure, such as the division of the genome into multiple chromosomes, have been reported in a diverse set of metazoans. We report on the two circle multipartite mt genome of a newly described cyst nematode, Globodera ellingtonae. RESULTS: The G. ellingtonae mt genome was found to be comprised of two circles, each larger than any other multipartite circular mt chromosome yet reported, and both were larger than the single mt circle of the model nematode Caenorhabditis elegans. The genetic content of the genome was disproportionately divided between the two circles, although they shared a ~6.5 kb non-coding region. The 17.8 kb circle (mtDNA-I) contained ten protein-coding genes and two tRNA genes, whereas the 14.4 kb circle (mtDNA-II) contained two protein-coding genes, 20 tRNA genes and both rRNA genes. Perhaps correlated with this division of genetic content, the copy number of mtDNA-II was more than four-fold that of mtDNA-I in individual nematodes. The difference in copy number increased between second-stage and fourth-stage juveniles. CONCLUSIONS: The segregation of gene types to different mt circles in G. ellingtonae could provide benefit by localizing gene functional types to independent transcriptional units. This is the first report of both two-circle and several-circle mt genomes within a single genus. The differential copy number associated with this multipartite mt organization could provide a model system for deconstructing mechanisms regulating mtDNA copy number both in somatic cells and during germline development.