Fragmentation in mitochondrial genomes in relation to elevated sequence divergence and extreme rearrangements

Background A single circular mitochondrial (mt) genome is a common feature across most metazoans. The mt-genome includes protein-coding genes involved in oxidative phosphorylation, as well as RNAs necessary for translation of mt-RNAs, whose order and number are highly conserved across animal clades, with few known exceptions of alternative mt-gene order or mt-genome architectures. One such exception consists of the fragmented mitochondrial genome, a type of genome architecture where mt-genes are split across two or more mt-chromosomes. However, the origins of mt-genome fragmentation and its effects on mt-genome evolution are unknown. Here, we investigate these origin and potential mechanisms underlying mt-genome fragmentation, focusing on a genus of booklice, Liposcelis, which exhibits elevated sequence divergence, frequent rearrangement of mt-gene order, and fragmentation of the mt genome, and compare them to other Metazoan clades. Results We found this genus Liposcelis exhibits very low conservation of mt-gene order across species, relative to other metazoans. Levels of gene order rearrangement were, however, unrelated to whether or not mt-genomes were fragmented or intact, suggesting mitochondrial genome fragmentation is not affecting mt-gene order directly. We further investigated possible mechanisms underpinning these patterns and revealed very high conservation of non-coding sequences at the edges of multiple recombination regions across populations of one particular Liposcelis species, supportive of a hypothesis that mt-fragmentation arises from recombination errors between mt-genome copies. We propose these errors may arise as a consequence of a heightened mutation rate in clades exhibiting mt-fragmentation. Consistent with this, we observed a striking pattern across three Metazoan phyla (Arthropoda, Nematoda, Cnidaria) characterised by members exhibiting high levels of mt-gene order rearrangement and cases of mt-fragmentation, whereby the mt-genomes of species more closely related to species with fragmented mt-genomes diverge more rapidly despite experiencing strong purifying selection. Conclusions We showed that contrary to expectations, mt-genome fragmentation is not correlated with the increase in mt-genome rearrangements. Furthermore, we present evidence that fragmentation of the mt-genome may be part of a general relaxation of a natural selection on the mt-genome, thus providing new insights into the origins of mt-genome fragmentation and evolution.

All ANIs are not created equal: implications for prokaryotic species boundaries and integration of ANIs into polyphasic taxonomy

In prokaryotic taxonomy, a set of criteria is commonly used to delineate species. These criteria are generally based on cohesion at the phylogenetic, phenotypic and genomic levels. One such criterion shown to have promise in the genomic era is average nucleotide identity (ANI), which provides an average measure of similarity across homologous regions shared by a pair of genomes. However, despite the popularity and relative ease of using this metric, ANI has undergone numerous refinements, with variations in genome fragmentation, homologue detection parameters and search algorithms. To test the robustness of a 95–96 % species cut-off range across all the commonly used ANI approaches, seven different methods were used to calculate ANI values for intra- and interspecies datasets representing three classes in the Proteobacteria . As a reference point, these methods were all compared to the widely used blast-based ANI (i.e. ANIb as implemented in JSpecies), and regression analyses were performed to investigate the correlation of these methods to ANIb with more than 130000 individual data points. From these analyses, it was clear that ANI methods did not provide consistent results regarding the conspecificity of isolates. Most of the methods investigated did not correlate perfectly with ANIb, particularly between 90 and 100% identity, which includes the proposed species boundary. There was also a difference in the correlation of methods for the different taxon sets. Our study thus suggests that the specific approach employed needs to be considered when ANI is used to delineate prokaryotic species. We furthermore suggest that one would first need to determine an appropriate cut-off value for a specific taxon set, based on the intraspecific diversity of that group, before conclusions on conspecificity of isolates can be made, and that the resulting species hypotheses be confirmed with analyses based on evolutionary history as part of the polyphasic approach to taxonomy.

A single-cell genome reveals diplonemid-like ancestry of kinetoplastid mitochondrial gene structure

Euglenozoa comprises euglenids, kinetoplastids, and diplonemids, with each group exhibiting different and highly unusual mitochondrial genome organizations. Although they are sister groups, kinetoplastids and diplonemids have very distinct mitochondrial genome architectures, requiring widespread insertion/deletion RNA editing and extensive trans -splicing, respectively, in order to generate functional transcripts. The evolutionary history by which these differing processes arose remains unclear. Using single-cell genomics, followed by small sub unit ribosomal DNA and multigene phylogenies, we identified an isolated marine cell that branches on phylogenetic trees as a sister to known kinetoplastids. Analysis of single-cell amplified genomic material identified multiple mitochondrial genome contigs. These revealed a gene architecture resembling that of diplonemid mitochondria, with small fragments of genes encoded out of order and or on different contigs, indicating that these genes require extensive trans -splicing. Conversely, no requirement for kinetoplastid-like insertion/deletion RNA-editing was detected. Additionally, while we identified some proteins so far only found in kinetoplastids, we could not unequivocally identify mitochondrial RNA editing proteins. These data invite the hypothesis that extensive genome fragmentation and trans -splicing were the ancestral states for the kinetoplastid-diplonemid clade but were lost during the kinetoplastid radiation. This study demonstrates that single-cell approaches can successfully retrieve lineages that represent important new branches on the tree of life, and thus can illuminate major evolutionary and functional transitions in eukaryotes. This article is part of a discussion meeting issue ‘Single cell ecology’.

Targeted genome fragmentation with CRISPR/Cas9 improves hybridization capture, reduces PCR bias, and enables efficient high-accuracy sequencing of small targets

ABSTRACTCurrent next-generation sequencing techniques suffer from inefficient target enrichment and frequent errors. To address these issues, we have developed a targeted genome fragmentation approach based on CRISPR/Cas9 digestion. By designing all fragments to similar lengths, regions of interest can be size-selected prior to library preparation, increasing hybridization capture efficiency. Additionally, homogenous length fragments reduce PCR bias and maximize read usability. We combine this novel target enrichment approach with ultra-accurate Duplex Sequencing. The result, termed CRISPR-DS, is a robust targeted sequencing technique that overcomes the inherent challenges of small target enrichment and enables the detection of ultra-low frequency mutations with small DNA inputs.

Genome expansion via lineage splitting and genome reduction in the cicada endosymbiont Hodgkinia

Comparative genomics from mitochondria, plastids, and mutualistic endosymbiotic bacteria has shown that the stable establishment of a bacterium in a host cell results in genome reduction. Although many highly reduced genomes from endosymbiotic bacteria are stable in gene content and genome structure, organelle genomes are sometimes characterized by dramatic structural diversity. Previous results from Candidatus Hodgkinia cicadicola, an endosymbiont of cicadas, revealed that some lineages of this bacterium had split into two new cytologically distinct yet genetically interdependent species. It was hypothesized that the long life cycle of cicadas in part enabled this unusual lineage-splitting event. Here we test this hypothesis by investigating the structure of the Ca. Hodgkinia genome in one of the longest-lived cicadas, Magicicada tredecim. We show that the Ca. Hodgkinia genome from M. tredecim has fragmented into multiple new chromosomes or genomes, with at least some remaining partitioned into discrete cells. We also show that this lineage-splitting process has resulted in a complex of Ca. Hodgkinia genomes that are 1.1-Mb pairs in length when considered together, an almost 10-fold increase in size from the hypothetical single-genome ancestor. These results parallel some examples of genome fragmentation and expansion in organelles, although the mechanisms that give rise to these extreme genome instabilities are likely different.