Diversity and evolution of sex determination systems in terrestrial isopods

Wolbachia: Advancing into a Second Century

Wolbachia pipientis had its scientific debut nearly a century ago and has recently emerged as a target for therapeutic treatment of filarial infections and an attractive tool for control of arthropod pests. Wolbachia was known as a biological entity before DNA was recognized as the molecule that carries the genetic information on which life depends, and before arthropods and nematodes were grouped in the Ecdysozoa. Today, some investigators consider Wolbachia the most abundant endosymbiont on earth, given the numbers of its hosts and its diverse mutualistic, commensal, and parasitic roles in their life histories. Recent advances in molecular technologies have revolutionized our understanding of Wolbachia and its associated reproductive phenotypes. New models have emerged for its investigation, and substantial progress has been made towards Wolbachia-based interventions in medicine and agriculture. Here I introduce Wolbachia, with a focus on aspects of its biology that are covered in greater detail in subsequent chapters.

The Scorpionfly (Panorpa cognata) Genome Highlights Conserved and Derived Features of the Peculiar Dipteran X Chromosome

Many insects carry an ancient X chromosome-the Drosophila Muller element F-that likely predates their origin. Interestingly, the X has undergone turnover in multiple fly species (Diptera) after being conserved for more than 450 My. The long evolutionary distance between Diptera and other sequenced insect clades makes it difficult to infer what could have contributed to this sudden increase in rate of turnover. Here, we produce the first genome and transcriptome of a long overlooked sister-order to Diptera: Mecoptera. We compare the scorpionfly Panorpa cognata X-chromosome gene content, expression, and structure to that of several dipteran species as well as more distantly related insect orders (Orthoptera and Blattodea). We find high conservation of gene content between the mecopteran X and the dipteran Muller F element, as well as several shared biological features, such as the presence of dosage compensation and a low amount of genetic diversity, consistent with a low recombination rate. However, the 2 homologous X chromosomes differ strikingly in their size and number of genes they carry. Our results therefore support a common ancestry of the mecopteran and ancestral dipteran X chromosomes, and suggest that Muller element F shrank in size and gene content after the split of Diptera and Mecoptera, which may have contributed to its turnover in dipteran insects.

Evolution of sex determination in crustaceans

Sex determination (SD) involves mechanisms that determine whether an individual will develop into a male, female, or in rare cases, hermaphrodite. Crustaceans harbor extremely diverse SD systems, including hermaphroditism, environmental sex determination (ESD), genetic sex determination (GSD), and cytoplasmic sex determination (e.g., Wolbachia controlled SD systems). Such diversity lays the groundwork for researching the evolution of SD in crustaceans, i.e., transitions among different SD systems. However, most previous research has focused on understanding the mechanism of SD within a single lineage or species, overlooking the transition across different SD systems. To help bridge this gap, we summarize the understanding of SD in various clades of crustaceans, and discuss how different SD systems might evolve from one another. Furthermore, we review the genetic basis for transitions between different SD systems (i.e., Dmrt genes) and propose the microcrustacean Daphnia (clade Branchiopoda) as a model to study the transition from ESD to GSD.

Heterogeneous distribution of sex ratio distorters in natural populations of the isopod Armadillidium vulgare

In the isopod Armadillidium vulgare, many females produce progenies with female-biased sex ratios, owing to two feminizing sex ratio distorters (SRD): Wolbachia endosymbionts and the f element. We investigated the distribution and population dynamics of these SRD and mitochondrial DNA variation in 16 populations from Europe and Japan. Confirming and extending results from the 1990s, we found that the SRD are present at variable frequencies in populations and that the f element is overall more frequent than Wolbachia. The two SRD never co-occur at high frequency in any population, suggesting an apparent mutual exclusion. We also detected Wolbachia or the f element in some males, which probably reflects insufficient titer to induce feminization or presence of masculinizing alleles. Our results are consistent with a single integration event of a Wolbachia genome in the A. vulgare genome at the origin of the f element, which contradicts an earlier hypothesis of frequent losses and gains. We identified strong linkage between Wolbachia strains and mitochondrial haplotypes, but no association between the f element and mitochondrial background. Our results open new perspectives on SRD evolutionary dynamics in A. vulgare, the evolution of genetic conflicts and their impact on the variability of sex determination systems.

Tire microplastics exposure in soil induces changes in expression profile of immune-related genes in terrestrial crustacean Porcellio scaber

Tire particles pose a potential threat to terrestrial organisms because they are deposited in large quantities in the soil by tire wear abrasion, and moreover their chemical complexity poses an additional risk. Microplastics can affect several physiological processes in organisms, including those related to immunity. Therefore, we investigated the expression profile of selected immune-related genes (MnSod, Manganese Superoxide dismutase; Cat, Catalase; CypG, Cyclophilin G; Nos, Nitric oxide synthase; Ppae2a, Prophenoloxidase-activating enzyme 2a; Dscam, Down syndrome cell adhesion molecule; Myd88, Myeloid-differentiation factor 88; Toll4, Toll-like receptor 4; Mas-like, Masquerade-like protein) in haemocytes and the digestive gland hepatopancreas of terrestrial crustacean Porcellio scaber after two different time exposures (4 and 14 days) to tire particles in soil. Our results reveal for the first time the response of P. scaber after microplastic exposure at the transcriptome level. We observed time- and tissue-dependent changes in the expression of the analysed genes, with more pronounced alterations in haemocytes after 14 days of exposure. Some minor changes were also observed in hepatopancreas after 4 days. Changes in the expression profile of the analysed genes are a direct indication of a modulated immune status of the test organism, which, however, does not represent an adverse effect on the test organism under the given conditions. Nevertheless, the question remains whether the observed change in immune status affects the immunocompetence of the test organism.

Evolutionary transition to XY sex chromosomes associated with Y-linked duplication of a male hormone gene in a terrestrial isopod

Sex chromosomes are highly variable in some taxonomic groups, but the evolutionary mechanisms underlying this diversity are not well understood. In terrestrial isopod crustaceans, evolutionary turnovers in sex chromosomes are frequent, possibly caused by Wolbachia, a vertically-transmitted endosymbiont causing male-to-female sex reversal. Here, we use surgical manipulations and genetic crosses, plus genome sequencing, to examine sex chromosomes in the terrestrial isopod Trachelipus rathkei. Although an earlier cytogenetics study suggested a ZZ/ZW sex chromosome system in this species, we surprisingly find multiple lines of evidence that in our study population, sex is determined by an XX/XY system. Consistent with a recent evolutionary origin for this XX/XY system, the putative male-specific region of the genome is small. The genome shows evidence of Y-linked duplications of the gene encoding the androgenic gland hormone, a major component of male sexual differentiation in isopods. Our analyses also uncover sequences horizontally acquired from past Wolbachia infections, consistent with the hypothesis that Wolbachia may have interfered with the evolution of sex determination in T. rathkei. Overall, these results provide evidence for the co-occurrence of multiple sex chromosome systems within T. rathkei, further highlighting the relevance of terrestrial isopods as models for the study of sex chromosome evolution.

Insulin-like androgenic gland hormone from the shrimp Fenneropenaeus merguiensis: Expression, gene organization and transcript variants

Male sex differentiation in the crustacean is best known to be controlled by the insulin-like androgenic gland hormone (IAG). In this report, the cDNA and gene of the shrimp Fenneropenaeus merguiensis FmIAG were studied and characterized. FmIAG gene shares a high sequence identity in the coding region as well as the promoter region with that of F. chinensis. FmIAG gene is most likely consists of 5 exons and 4 introns. The cDNA reported here is the most abundant transcript that retained cryptic intron 4. The use of different splicing acceptor sites in exon 2 can produce a long-form FmIAG transcript variant with 6 additional amino acids inserted. Splicing of cryptic intron 4 would produce a transcript variant with a different C-terminal end. Therefore 4 different FmIAG transcripts can be produced from the FmIAG gene. During the molt cycle, the expression level of FmIAG was low in the early intermolt, increase steadily towards the late premolt and decreased rapidly in the early postmolt. In addition to the androgenic gland, FmIAG is also expressed in the hepatopancreas and ovary of adult females. Unilateral eyestalk ablation caused a significant increase in FmIAG transcript suggesting that the eyestalk consists of inhibiting factor(s) that suppressesFmIAGexpression. To explore the function of FmIAG in males, injection of FmIAG dsRNA knock-down the expression of FmIAG and up-regulated the expression of the vitellogenin gene in the testis and hepatopancreas. Interestingly a CHH-like gene identified in the androgenic gland was down-regulated. CHH-like gene knock-down resulted in altered expression of FmIAG in males suggesting that the CHH-like may be involved in FmIAG's regulation. RT-PCR with specific primers to the different transcript variant were used to determine if there is an association of different sizes of male and the type of IAG transcript. Results indicated that a high percentage of the large male shrimp expressed the long-form of FmIAG. The results suggested that FmIAG may be useful as a size marker for male shrimp aquaculture. In summary, the results of this study have expanded our knowledge of shrimp insulin-like androgenic gland hormone in male sex development and its potential role as a marker gene for growth regulation in shrimp.

Sex Determination and Differentiation in Decapod and Cladoceran Crustaceans: An Overview of Endocrine Regulation

Mechanisms underlying sex determination and differentiation in animals are known to encompass a diverse array of molecular clues. Recent innovations in high-throughput sequencing and mass spectrometry technologies have been widely applied in non-model organisms without reference genomes. Crustaceans are no exception. They are particularly diverse among the Arthropoda and contain a wide variety of commercially important fishery species such as shrimps, lobsters and crabs (Order Decapoda), and keystone species of aquatic ecosystems such as water fleas (Order Branchiopoda). In terms of decapod sex determination and differentiation, previous approaches have attempted to elucidate their molecular components, to establish mono-sex breeding technology. Here, we overview reports describing the physiological functions of sex hormones regulating masculinization and feminization, and gene discovery by transcriptomics in decapod species. Moreover, this review summarizes the recent progresses of studies on the juvenile hormone-driven sex determination system of the branchiopod genus Daphnia, and then compares sex determination and endocrine systems between decapods and branchiopods. This review provides not only substantial insights for aquaculture research, but also the opportunity to re-organize the current and future trends of this field.

Shedding Light on the Antimicrobial Peptide Arsenal of Terrestrial Isopods: Focus on Armadillidins, a New Crustacean AMP Family

In crustaceans, antimicrobial peptides (AMPs) are clustered into four major groups according to their amino acid composition and structure: (1) single-domain peptides containing cysteine residues such as anti-lipopolysaccharide-factor (ALF), (2) multi-domain or chimeric AMPs such as crustins, (3) non-conventional AMPs, and (4) linear single-domain AMPs. The majority of AMPs has been described in commercially exploited crustaceans, particularly decapods living in aquatic environments (crab, shrimp, lobster, and crayfish). Here, we aimed at establishing the AMPs repertoire of terrestrial isopods (Oniscidea), an original suborder of crustaceans adapted to life outside of the aquatic environment. Using transcriptomic data from 21 species, we identified 110 ALF and 73 crustin sequences. We also characterized the full-length sequence of armadillidins from 17 species, similar to the AMP previously described in the terrestrial isopod Armadillidium vulgare. Furthermore, we tested the antimicrobial activity of three armadillidin peptides characterized from three distantly related species. This analysis revealed similar activity spectra against pathogens, despite extensive structural variation among the tested peptides. In addition to conventional crustacean AMPs, our work highlights armadillidins as a new and independent family of AMPs specific to the Oniscidea, thus opening new perspectives concerning the study of the immune system of terrestrial isopods.

How to identify sex chromosomes and their turnover

Although sex is a fundamental component of eukaryotic reproduction, the genetic systems that control sex determination are highly variable. In many organisms the presence of sex chromosomes is associated with female or male development. Although certain groups possess stable and conserved sex chromosomes, others exhibit rapid sex chromosome evolution, including transitions between male and female heterogamety, and turnover in the chromosome pair recruited to determine sex. These turnover events have important consequences for multiple facets of evolution, as sex chromosomes are predicted to play a central role in adaptation, sexual dimorphism, and speciation. However, our understanding of the processes driving the formation and turnover of sex chromosome systems is limited, in part because we lack a complete understanding of interspecific variation in the mechanisms by which sex is determined. New bioinformatic methods are making it possible to identify and characterize sex chromosomes in a diverse array of non-model species, rapidly filling in the numerous gaps in our knowledge of sex chromosome systems across the tree of life. In turn, this growing data set is facilitating and fueling efforts to address many of the unanswered questions in sex chromosome evolution. Here, we synthesize the available bioinformatic approaches to produce a guide for characterizing sex chromosome system and identity simultaneously across clades of organisms. Furthermore, we survey our current understanding of the processes driving sex chromosome turnover, and highlight important avenues for future research.

Sex chromosomes control vertical transmission of feminizing Wolbachia symbionts in an isopod

Microbial endosymbiosis is widespread in animals, with major ecological and evolutionary implications. Successful symbiosis relies on efficient vertical transmission through host generations. However, when symbionts negatively affect host fitness, hosts are expected to evolve suppression of symbiont effects or transmission. Here, we show that sex chromosomes control vertical transmission of feminizing Wolbachia endosymbionts in the isopod Armadillidium nasatum. Theory predicts that the invasion of an XY/XX species by cytoplasmic sex ratio distorters is unlikely because it leads to fixation of the unusual (and often lethal or infertile) YY genotype. We demonstrate that A. nasatum X and Y sex chromosomes are genetically highly similar and that YY individuals are viable and fertile, thereby enabling Wolbachia spread in this XY-XX species. Nevertheless, we show that Wolbachia cannot drive fixation of YY individuals, because infected YY females do not transmit Wolbachia to their offspring, unlike XX and XY females. The genetic basis fits the model of a Y-linked recessive allele (associated with an X-linked dominant allele), in which the homozygous state suppresses Wolbachia transmission. Moreover, production of all-male progenies by infected YY females restores a balanced sex ratio at the host population level. This suggests that blocking of Wolbachia transmission by YY females may have evolved to suppress feminization, thereby offering a whole new perspective on the evolutionary interplay between microbial symbionts and host sex chromosomes.

Unravelling the molecular mechanisms of nickel in woodlice

During the last few years, there has been an alarming increase in the amount of nickel (Ni) being released into the environment, primarily due to its use in the production of stainless steel but also from other sources such as batteries manufacturing and consequent disposal. The established biotic ligand models provide precise estimates for Ni bioavailability, in contrast, studies describing the mechanisms underpinning toxicological effect of Ni are scarce. This study exploits RNA-seq to determine the transcriptomic responses of isopods using Porcellionides pruinosus as an example of a terrestrial metal-resistant woodlouse. Furthermore, the recently proposed model for Ni adverse outcome pathways (Ni-AOP) presents an unprecedented opportunity to fit isopod responses to Ni toxicity and define Porcellionides pruinosus as a metalomic model. Prior to this study, P. pruinosus represented an important environmental sentinel, though lacking genetic/omic data. The reference transcriptome generated here thus represents a major advance and a novel resource. A detailed annotation of the transcripts obtained is presented together with the homology to genes/gene products from Metazoan and Arthropoda phylum, Gene Ontology (GO) classification, clusters of orthologous groups (COG) and assignment to KEGG metabolic pathways. The differential gene expression comparison was determined in response to nickel (Ni) exposure and used to derive the enriched pathways and processes. It revealed a significant impact on ion trafficking and storage, oxidative stress, neurotoxicity, reproduction impairment, genetics and epigenetics. Many of the processes observed support the current Ni-AOP although the data highlights that the current model can be improved by including epigenetic endpoints, which represents key chronic risks under a scenario of Ni toxicity.

Construction of a Genomic Bacterial Artificial Chromosome (BAC) Library for the Prawn Macrobrachium rosenbergii and Initial Analysis of ZW Chromosome-Derived BAC Inserts

Knowledge on sex determination has proven valuable for commercial production of the prawn Macrobrachium rosenbergii due to sex dimorphism of the male and female individuals. Previous studies indicated that prawn sex is determined by a ZW-ZZ chromosomal system, but no genomic information is available for the sex chromosome. Herein, we constructed a genomic bacterial artificial chromosome (BAC) library and identified the ZW-derived BAC clones for initial analysis of the sex chromosomal DNA sequence. The arrayed BAC library contains 200,448 clones with average insert size of 115.4 kb, corresponding to ∼ 4× coverage of the estimated 5.38 Gb genome. Based on a short female-specific marker, a Z- and a W-fragment were retrieved with the genomic walking method. Screening the BAC library using a ZW-specific marker as probe resulted in 12 positive clones. From these, a Z-derived (P331M17) and a W-derived (P122G2) BAC clones were randomly selected and sequenced by PacBio method. We report the construction of a large insert, deep-coverage, and high-quality BAC library for M. rosenbergii that provides a useful resource for positional cloning of target genes, genomic organization, and comparative genomics analysis. Our study not only confirmed the ZW/ZZ system but also discovered sex-linked genes on ZW chromosomes for the first time, contributing to a comprehensive understanding of the genomic structure of sex chromosomes in M. rosenbergii.

A bioinformatics approach to identifying Wolbachia infections in arthropods

Wolbachia is the most widespread endosymbiont, infecting >20% of arthropod species, and capable of drastically manipulating the host's reproductive mechanisms. Conventionally, diagnosis has relied on PCR amplification; however, PCR is not always a reliable diagnostic technique due to primer specificity, strain diversity, degree of infection and/or tissue sampled. Here, we look for evidence of Wolbachia infection across a wide array of arthropod species using a bioinformatic approach to detect the Wolbachia genes ftsZ, wsp, and the groE operon in next-generation sequencing samples available through the NCBI Sequence Read Archive. For samples showing signs of infection, we attempted to assemble entire Wolbachia genomes, and in order to better understand the relationships between hosts and symbionts, phylogenies were constructed using the assembled gene sequences. Out of the 34 species with positively identified infections, eight species of arthropod had not previously been recorded to harbor Wolbachia infection. All putative infections cluster with known representative strains belonging to supergroup A or B, which are known to only infect arthropods. This study presents an efficient bioinformatic approach for post-sequencing diagnosis and analysis of Wolbachia infection in arthropods.

Sex chromosome turnovers and genetic drift: a simulation study

The recent advances of new genomic technologies have enabled the identification and characterization of sex chromosomes in an increasing number of nonmodel species, revealing that many plants and animals undergo frequent sex chromosome turnovers. What evolutionary forces drive these turnovers remains poorly understood, but it was recently proposed that drift might play a more important role than generally assumed. We analysed the dynamics of different types of turnovers using individual-based simulations and show that when mediated by genetic drift, turnovers are usually easier to achieve than substitutions at neutral markers, but that their dynamics and relative likelihoods vary with the type of the resident and emergent sex chromosome system (XY and/or ZW) and the dominance relationships among the sex-determining factors. Focusing on turnovers driven by epistatically dominant mutations, we find that drift-mediated turnovers that preserve the heterogamety pattern are 2-4× more likely than those along which the heterogametic sex changes. This ratio nevertheless decreases along with effective population size and can even reverse in case of extreme polygyny. This can be attributed to a 'drift-induced' selective force, known to influence transitions between male and female heterogamety, but which according to our study does not affect turnovers that preserve the heterogametic sex.

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.

De novo assembly, characterization and annotation for the transcriptome of Sphaeroma terebrans and microsatellite marker discovery

Sphaeroma terebrans, an economically and ecologically important marine wood-boring isopod, is mainly distributed in tropical and subtropical mangroves. Nevertheless, available genomic and transcriptomic information for this isopod is extremely deficient. Here, we first performed the assembly of S. terebrans transcriptome by Illumina sequencing. A total of 51,092 high-quality unigenes with an average length of 641 bp were obtained and 19,915 unigenes were successfully annotated in four public databases. Only 9932 out of 19915 unigenes were commonly annotated by all four databases. In addition, 9609 unigenes were categorized into 54 function categories of Gene Ontology (GO), and 14,512 unigenes were successfully grouped into 25 functional categories of the EuKaryotic Orthologous Groups (KOG) database. Moreover, 11,507 unigenes were assigned to 228 pathways of the Kyoto Encyclopedia of Genes and Genomes (KEGG). Out of 51,092 unigenes, a total of 4257 different microsatellites with motifs range from di- to hexa-nucleotide were identified from 3324 unigene sequences. Among 64 primer pairs selected for validation, 35 were successful in PCR amplification and 13 exhibited obvious repeat polymorphisms in the wild population of S. terebrans in Dongzhaigang (Hainan Island). The transcriptome dataset and the identified microsatellite markers (SSRs) will provide abundant information for researches on the discovery of new genes, metabolic mechanism and genetic diversity of S. terebrans.

Drift-Induced Selection Between Male and Female Heterogamety

Evolutionary transitions between male and female heterogamety are common in both vertebrates and invertebrates. Theoretical studies of these transitions have found that, when all genotypes are equally fit, continuous paths of intermediate equilibria link the two sex chromosome systems. This observation has led to a belief that neutral evolution along these paths can drive transitions, and that arbitrarily small fitness differences among sex chromosome genotypes can determine the system to which evolution leads. Here, we study stochastic evolutionary dynamics along these equilibrium paths. We find non-neutrality, both in transitions retaining the ancestral pair of sex chromosomes, and in those creating a new pair. In fact, substitution rates are biased in favor of dominant sex determining chromosomes, which fix with higher probabilities than mutations of no effect. Using diffusion approximations, we show that this non-neutrality is a result of "drift-induced selection" operating at every point along the equilibrium paths: stochastic jumps off the paths return with, on average, a directional bias in favor of the dominant segregating sex chromosome. Our results offer a novel explanation for the observed preponderance of dominant sex determining genes, and hint that drift-induced selection may be a common force in standard population genetic systems.

Evolutionary Significance of Wolbachia-to-Animal Horizontal Gene Transfer: Female Sex Determination and the f Element in the Isopod Armadillidium vulgare

An increasing number of horizontal gene transfer (HGT) events from bacteria to animals have been reported in the past years, many of which involve Wolbachia bacterial endosymbionts and their invertebrate hosts. Most transferred Wolbachia genes are neutrally-evolving fossils embedded in host genomes. A remarkable case of Wolbachia HGT for which a clear evolutionary significance has been demonstrated is the “f element”, a nuclear Wolbachia insert involved in female sex determination in the terrestrial isopod Armadillidium vulgare. The f element represents an instance of bacteria-to-animal HGT that has occurred so recently that it was possible to infer the donor (feminizing Wolbachia closely related to the wVulC Wolbachia strain of A. vulgare) and the mechanism of integration (a nearly complete genome inserted by micro-homology-mediated recombination). In this review, we summarize our current knowledge of the f element and discuss arising perspectives regarding female sex determination, unstable inheritance, population dynamics and the molecular evolution of the f element. Overall, the f element unifies three major areas in evolutionary biology: symbiosis, HGT and sex determination. Its characterization highlights the tremendous impact sex ratio distorters can have on the evolution of sex determination mechanisms and sex chromosomes in animals and plants.

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.