Multiple Conserved Heteroplasmic Sites in tRNA Genes in the Mitochondrial Genomes of Terrestrial Isopods (Oniscidea)

Long-read sequencing reveals atypical mitochondrial genome structure in a New Zealand marine isopod

Most animal mitochondrial genomes are small, circular and structurally conserved. However, recent work indicates that diverse taxa possess unusual mitochondrial genomes. In Isopoda, species in multiple lineages have atypical and rearranged mitochondrial genomes. However, more species of this speciose taxon need to be evaluated to understand the evolutionary origins of atypical mitochondrial genomes in this group. In this study, we report the presence of an atypical mitochondrial structure in the New Zealand endemic marine isopod, Isocladus armatus. Data from long- and short-read DNA sequencing suggest that I. armatus has two mitochondrial chromosomes. The first chromosome consists of two mitochondrial genomes that have been inverted and fused together in a circular form, and the second chromosome consists of a single mitochondrial genome in a linearized form. This atypical mitochondrial structure has been detected in other isopod lineages, and our data from an additional divergent isopod lineage (Sphaeromatidae) lends support to the hypothesis that atypical structure evolved early in the evolution of Isopoda. Additionally, we find that an asymmetrical site previously observed across many species within Isopoda is absent in I. armatus, but confirm the presence of two asymmetrical sites recently reported in two other isopod species.

Disrupted architecture and fast evolution of the mitochondrial genome of Argeia pugettensis (Isopoda): implications for speciation and fitness

BACKGROUND

Argeia pugettensis is an isopod species that parasitizes other crustaceans. Its huge native geographic range spans the Pacific from China to California, but molecular data are available only for a handful of specimens from North-American populations. We sequenced and characterised the complete mitogenome of a specimen collected in the Yellow Sea.

RESULTS

It exhibited a barcode (cox1) similarity level of only 87-89% with North-American populations, which is unusually low for conspecifics. Its mitogenome is among the largest in isopods (≈16.5 Kbp), mostly due to a large duplicated palindromic genomic segment (2 Kbp) comprising three genes. However, it lost a segment comprising three genes, nad4L-trnP-nad6, and many genes exhibited highly divergent sequences in comparison to isopod orthologues, including numerous mutations, deletions and insertions. Phylogenetic and selection analyses corroborated that this is one of the handful of most rapidly evolving available isopod mitogenomes, and that it evolves under highly relaxed selection constraints (as opposed to positive selection). However, its nuclear 18S gene is highly conserved, which suggests that rapid evolution is limited to its mitochondrial genome. The cox1 sequence analysis indicates that elevated mitogenomic evolutionary rates are not shared by North-American conspecifics, which suggests a breakdown of cox1 barcoding in this species.

CONCLUSIONS

A highly architecturally disrupted mitogenome and decoupling of mitochondrial and nuclear rates would normally be expected to have strong negative impacts on the fitness of the organism, so the existence of this lineage is a puzzling evolutionary question. Additional studies are needed to assess the phylogenetic breadth of this disrupted mitochondrial architecture and its impact on fitness.

The complete mitochondrial genome of Cymothoa indica has a highly rearranged gene order and clusters at the very base of the Isopoda clade

As a result of great diversity in life histories and a large number of described species, taxonomic and phylogenetic uncertainty permeates the entire crustacean order of Isopoda. Large molecular datasets capable of providing sufficiently high phylogenetic resolution, such as mitochondrial genomes (mitogenomes), are needed to infer their evolutionary history with confidence, but isopod mitogenomes remain remarkably poorly represented in public databases. We sequenced the complete mitogenome of Cymothoa indica, a species belonging to a family from which no mitochondrial genome was sequenced yet, Cymothoidae. The mitogenome (circular, 14484 bp, A+T = 63.8%) is highly compact, appears to be missing two tRNA genes (trnI and trnE), and exhibits a unique gene order with a large number of rearrangements. High compactness and the existence of palindromes indicate that the mechanism behind these rearrangements might be associated with linearization events in its evolutionary history, similar to those proposed for isopods from the Armadillidium genus (Oniscidea). Isopods might present an important model system to study the proposed discontinuity in the dynamics of mitochondrial genomic architecture evolution. Phylogenetic analyses (Bayesian Inference and Maximum Likelihood) conducted using nucleotide sequences of all mitochondrial genes resolved Oniscidea and Cymothoida suborders as paraphyletic. Cymothoa indica was resolved as a sister group (basal) to all remaining isopods, which challenges the accepted isopod phylogeny, where Cymothoida are the most derived, and Phreatoicidea the most basal isopod group. There is growing evidence that Cymothoida suborder might be split into two evolutionary distant clades, with parasitic species being the most basal split in the Isopoda clade, but a much larger amount of molecular resources carrying a high phylogenetic resolution will be needed to infer the remarkably complex evolutionary history of this group of animals with confidence.

Investigating the Molecular Genetic Basis of Cytoplasmic Sex Determination Caused by Wolbachia Endosymbionts in Terrestrial Isopods

In animals, sexual differences between males and females are usually determined by sex chromosomes. Alternatively, sex may also be determined by vertically transmitted intracellular microbial endosymbionts. The best known cytoplasmic sex manipulative endosymbiont is Wolbachia which can, for instance, feminize genetic males into phenotypic females in the terrestrial isopod Armadillidium vulgare. However, the molecular genetic basis of cytoplasmic sex determination is unknown. To identify candidate genes of feminization induced by Wolbachia strain wVulC from A. vulgare, we sequenced the genome of Wolbachia strain wCon from Cylisticus convexus, the most closely related known Wolbachia strain to wVulC that does not induce feminization, and compared it to the wVulC genome. Then, we performed gene expression profiling of the 216 resulting wVulC candidate genes throughout host developmental stages in A. vulgare and the heterologous host C. convexus. We identified a set of 35 feminization candidate genes showing differential expression during host sexual development. Interestingly, 27 of the 35 genes are present in the f element, which is a piece of a feminizing Wolbachia genome horizontally transferred into the nuclear genome of A. vulgare and involved in female sex determination. Assuming that the molecular genetic basis of feminization by Wolbachia and the f element is the same, the 27 genes are candidates for acting as master sex determination genes in A. vulgare females carrying the f element.

The first complete mitochondrial genome of a parasitic isopod supports Epicaridea Latreille, 1825 as a suborder and reveals the less conservative genome of isopods

The complete mitochondrial genome sequence of the holoparasitic isopod Gyge ovalis (Shiino, 1939) has been determined. The mitogenome is 14,268 bp in length and contains 34 genes: 13 protein-coding genes, two ribosomal RNA, 19 tRNA and a control region. Three tRNA genes (trnE, trnI and trnS1) are missing. Most of the tRNA genes show secondary structures which derive from the usual cloverleaf pattern except for trnC which is characterised by the loss of the DHU-arm. Compared to the isopod ground pattern and Eurydice pulchra Leach, 1815 (suborder Cymothoida Wägele, 1989), the genome of G. ovalis shows few differences, with changes only around the control region. However, the genome of G. ovalis is very different from that of non-cymothoidan isopods and reveals that the gene order evolution in isopods is less conservative compared to other crustaceans. Phylogenic trees were constructed using maxiumum likelihood and Bayesian inference analyses based on 13 protein-coding genes. The results do not support the placement of G. ovalis with E. pulchra and Bathynomus sp. in the same suborder; rather, G. ovalis appears to have a closer relationship to Ligia oceanica (Linnaeus, 1767), but this result suggests a need for more data and further analysis. Nevertheless, these results cast doubt that Epicaridea Latreille, 1825 can be placed as an infraorder within the suborder Cymothoida, and Epicaridea appears to also deserve subordinal rank. Further development of robust phylogenetic relationships across Isopoda Latreille, 1817 will require more genetic data from a greater diversity of taxa belonging to all isopod suborders.

Untangling Heteroplasmy, Structure, and Evolution of an Atypical Mitochondrial Genome by PacBio Sequencing

The highly compact mitochondrial (mt) genome of terrestrial isopods (Oniscidae) presents two unusual features. First, several loci can individually encode two tRNAs, thanks to single nucleotide polymorphisms at anticodon sites. Within-individual variation (heteroplasmy) at these loci is thought to have been maintained for millions of years because individuals that do not carry all tRNA genes die, resulting in strong balancing selection. Second, the oniscid mtDNA genome comes in two conformations: a ∼14 kb linear monomer and a ∼28 kb circular dimer comprising two monomer units fused in palindrome. We hypothesized that heteroplasmy actually results from two genome units of the same dimeric molecule carrying different tRNA genes at mirrored loci. This hypothesis, however, contradicts the earlier proposition that dimeric molecules result from the replication of linear monomers-a process that should yield totally identical genome units within a dimer. To solve this contradiction, we used the SMRT (PacBio) technology to sequence mirrored tRNA loci in single dimeric molecules. We show that dimers do present different tRNA genes at mirrored loci; thus covalent linkage, rather than balancing selection, maintains vital variation at anticodons. We also leveraged unique features of the SMRT technology to detect linear monomers closed by hairpins and carrying noncomplementary bases at anticodons. These molecules contain the necessary information to encode two tRNAs at the same locus, and suggest new mechanisms of transition between linear and circular mtDNA. Overall, our analyses clarify the evolution of an atypical mt genome where dimerization counterintuitively enabled further mtDNA compaction.

Diversity and evolution of sex determination systems in terrestrial isopods

Sex determination systems are highly variable in many taxa, sometimes even between closely related species. Yet the number and direction of transitions between these systems have seldom been characterized, and the underlying mechanisms are still poorly understood. Here we generated transcriptomes for 19 species of terrestrial isopod crustaceans, many of which are infected by Wolbachia bacterial endosymbionts. Using 88 single-copy orthologous genes, we reconstructed a fully resolved and dated phylogeny of terrestrial isopods. An original approach involving crossings of sex-reversed individuals allowed us to characterize the heterogametic systems of five species (one XY/XX and four ZW/ZZ). Mapping of these and previously known heterogametic systems onto the terrestrial isopod phylogeny revealed between 3 and 13 transitions of sex determination systems during the evolution of these taxa, most frequently from female to male heterogamety. Our results support that WW individuals are viable in many species, suggesting sex chromosomes are at an incipient stage of their evolution. Together, these data are consistent with the hypothesis that nucleo-cytoplasmic conflicts generated by Wolbachia endosymbionts triggered recurrent turnovers of sex determination systems in terrestrial isopods. They further establish terrestrial isopods as a model to study evolutionary transitions in sex determination systems and pave the way to molecularly characterize these systems.

The first complete mitogenome of the South China deep-sea giant isopod Bathynomus sp. (Crustacea: Isopoda: Cirolanidae) allows insights into the early mitogenomic evolution of isopods

In this study, the complete mitochondrial (mt) genome sequence of the South China deep‐sea giant isopod Bathynomus sp. was determined, and this study is the first to explore in detail the mt genome of a deep‐sea member of the order Isopoda. This species belongs to the genus Bathynomus, the members of which are saprophagous residents of the deep‐sea benthic environment; based on their large size, Bathynomus is included in the “supergiant group” of isopods. The mt genome of Bathynomus sp. is 14,965 bp in length and consists of 13 protein‐coding genes, two ribosomal RNA genes, only 18 transfer RNA genes, and a noncoding control region 362 bp in length, which is the smallest control region discovered in Isopoda to date. Although the overall genome organization is typical for metazoans, the mt genome of Bathynomus sp. shows a number of derived characters, such as an inversion of 10 genes when compared to the pancrustacean ground pattern. Rearrangements in some genes (e.g., cob, trnT, nad5, and trnF) are shared by nearly all isopod mt genomes analyzed thus far, and when compared to the putative isopod ground pattern, five rearrangements were found in Bathynomus sp. Two tRNAs exhibit modified secondary structures: The TΨC arm is absent from trnQ, and trnC lacks the DHU. Within the class Malacostraca, trnC arm loss is only found in other isopods. Phylogenetic analysis revealed that Bathynomus sp. (Cymothoida) and Sphaeroma serratum (Sphaeromatidea) form a single clade, although it is unclear whether Cymothoida is monophyletic or paraphyletic. Moreover, the evolutionary rate of Bathynomus sp. (dN/dS [nonsynonymous mutational rate/synonymous mutational rate] = 0.0705) is the slowest measured to date among Cymothoida, which may be associated with its relatively constant deep‐sea environment. Overall, our results may provide useful information for understanding the evolution of deep‐sea Isopoda species.

Candidate pathogenicity islands in the genome of 'Candidatus Rickettsiella isopodorum', an intracellular bacterium infecting terrestrial isopod crustaceans

The bacterial genus Rickettsiellabelongs to the order Legionellales in the Gammaproteobacteria, and consists of several described species and pathotypes, most of which are considered to be intracellular pathogens infecting arthropods. Two members of this genus, R. grylliand R. isopodorum, are known to infect terrestrial isopod crustaceans. In this study, we assembled a draft genomic sequence for R. isopodorum, and performed a comparative genomic analysis with R. grylli. We found evidence for several candidate genomic island regions in R. isopodorum, none of which appear in the previously available R. grylli genome sequence.Furthermore, one of these genomic island candidates in R. isopodorum contained a gene that encodes a cytotoxin partially homologous to those found in Photorhabdus luminescensand Xenorhabdus nematophilus (Enterobacteriaceae), suggesting that horizontal gene transfer may have played a role in the evolution of pathogenicity in Rickettsiella. These results lay the groundwork for future studies on the mechanisms underlying pathogenesis in R. isopodorum, and this system may provide a good model for studying the evolution of host-microbe interactions in nature.

Large gene overlaps and tRNA processing in the compact mitochondrial genome of the crustacean Armadillidium vulgare

A faithful expression of the mitochondrial DNA is crucial for cell survival. Animal mitochondrial DNA (mtDNA) presents a highly compact gene organization. The typical 16.5 kbp animal mtDNA encodes 13 proteins, 2 rRNAs and 22 tRNAs. In the backyard pillbug Armadillidium vulgare, the rather small 13.9 kbp mtDNA encodes the same set of proteins and rRNAs as compared to animal kingdom mtDNA, but seems to harbor an incomplete set of tRNA genes. Here, we first confirm the expression of 13 tRNA genes in this mtDNA. Then we show the extensive repair of a truncated tRNA, the expression of tRNA involved in large gene overlaps and of tRNA genes partially or fully integrated within protein-coding genes in either direct or opposite orientation. Under selective pressure, overlaps between genes have been likely favored for strong genome size reduction. Our study underlines the existence of unknown biochemical mechanisms for the complete gene expression of A. vulgare mtDNA, and of co-evolutionary processes to keep overlapping genes functional in a compacted mitochondrial genome.