Evolution of mitochondrial genomes in Baikalian amphipods

Cryptic diversity and speciation in an endemic copepod crustacean Harpacticella inopinata within Lake Baikal

Ancient lakes are hotspots of species diversity, posing challenges and opportunities for exploration of the dynamics of endemic diversification. Lake Baikal in Siberia, the oldest lake in the world, hosts a particularly rich crustacean fauna, including the largest known species flock of harpacticoid copepods with some 70 species. Here, we focused on exploring the diversity and evolution within a single nominal species, Harpacticella inopinata Sars, 1908, using molecular markers (mitochondrial COI, nuclear ITS1 and 28S rRNA) and a set of qualitative and quantitative morphological traits. Five major mitochondrial lineages were recognized, with model-corrected COI distances of 0.20-0.37. A concordant pattern was seen in the nuclear data set, and qualitative morphological traits also distinguish a part of the lineages. All this suggests the presence of several hitherto unrecognized cryptic taxa within the baikalian H. inopinata, with long independent histories. The abundances, distributions and inferred demographic histories were different among taxa. Two taxa, H. inopinata CE and H. inopinata CW, were widespread on the eastern and western coasts, respectively, and were largely allopatric. Patterns in mitochondrial variation, that is, shallow star-like haplotype networks, suggest these taxa have spread through the lake relatively recently. Three other taxa, H. inopinata RE, RW and RW2, instead were rare and had more localized distributions on either coast, but showed deeper intraspecies genealogies, suggesting older regional presence. The rare taxa were often found in sympatry with the others and occasionally introgressed by mtDNA from the common ones. The mitochondrial divergence between and within the H. inopinata lineages is still unexpectedly deep, suggesting an unusually high molecular rate. The recognition of true systematic diversity in the evaluation and management of ecosystems is important in hotspots, as it is everywhere else, while the translation of the diversity into a formal taxonomy remains a challenge.

The complete mitochondrial genome of freshwater gammarid Gammarus nipponensis (Crustacea: Amphipoda: Gammaridae)

This study presents the complete mitochondrial genome sequence of Gammarus nipponensis, a freshwater crustacean found in the western regions of Honshu, Shikoku and Kyushu in Japan. The entire genome is 16,429 bp in length, encoding a standard set of 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes, as well as the putative control regions. The mitochondrial genome of G. nipponensis is characterized by a high concentration of A and T nucleotides (67.1%). Notably, the mitogenome contains long TATTTTA repeats in the control region 2 at 686 bp long. This newly available genome information will be useful for studying the evolutionary relationships within the genus Gammarus and for understanding diversification among G. nipponensis populations.

Comparison of the evolutionary patterns of DNA repeats in ancient and young invertebrate species flocks of Lake Baikal

DNA repeat composition of low coverage (0.1-0.5) genomic libraries of four amphipods species endemic to Lake Baikal (East Siberia) and four endemic gastropod species of the fam. Baicaliidae have been compared to each other. In order to do so, a neighbor joining tree was inferred for each quartet of species (amphipods and mollusks) based on the ratio of repeat classes shared in each pair of species. The topology of this tree was compared to the phylogenies inferred for the same species from the concatenated protein-coding mitochondrial nucleotide sequences. In all species analyzed, the fraction of DNA repeats involved circa half of the genome. In relatively more ancient amphipods (most recent common ancestor, MRCA, existed approximately sixty millions years ago), the most abundant were species-specific repeats, while in much younger Baicaliidae (MRCA equal to ca. three millions years) most of the DNA repeats were shared among all four species. If the presence/absence of a repeat is regarded as a separate independent trait, and the ratio of shared to total numbers of repeats in a species pair is used as the measure of distance, the topology of the NJ tree is the same as the quartet phylogeny inferred for the mitogenomes protein coding nucleotide sequences. Meanwhile, in each group of species, a substantial number of repeats were detected pointing to the possibility of non-neutral evolution or a horizontal transfer between species occupying the same biotope. These repeats were shared by non-sister groups while being absent in the sister genomes. On the other hand, in such cases some traits of ecological significance were also shared.

Experimental Crossing Confirms Reproductive Isolation between Cryptic Species within Eulimnogammarus verrucosus (Crustacea: Amphipoda) from Lake Baikal

Ancient lakes are known speciation hotspots. One of the most speciose groups in the ancient Lake Baikal are gammaroid amphipods (Crustacea: Amphipoda: Gammaroidea). There are over 350 morphological species and subspecies of amphipods in Baikal, but the extent of cryptic variation is still unclear. One of the most common species in the littoral zone of the lake, Eulimnogammarus verrucosus (Gerstfeldt, 1858), was recently found to comprise at least three (pseudo)cryptic species based on molecular data. Here, we further explored these species by analyzing their mitogenome-based phylogeny, genome sizes with flow cytometry, and their reproductive compatibility. We found divergent times of millions of years and different genome sizes in the three species (6.1, 6.9 and 8 pg), further confirming their genetic separation. Experimental crossing of the western and southern species, which are morphologically indistinguishable and have adjacent ranges, showed their separation with a post-zygotic reproductive barrier, as hybrid embryos stopped developing roughly at the onset of gastrulation. Thus, the previously applied barcoding approach effectively indicated the separate biological species within E. verrucosus. These results provide new data for investigating genome evolution and highlight the need for precise tracking of the sample origin in any studies in this morphospecies.

Describing novel mitochondrial genomes of Antarctic amphipods

To date, only one mitogenome from an Antarctic amphipod has been published. Here, novel complete mitochondrial genomes (mitogenomes) of two morphospecies are assembled, namely, Charcotia amundseni and Eusirus giganteus. For the latter species, we have assembled two mitogenomes from different genetic clades of this species. The lengths of Eusirus and Charcotia mitogenomes range from 15,534 to 15,619 base pairs and their mitogenomes are composed of 13 protein coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and 1 putative control region CR. Some tRNAs display aberrant structures suggesting that minimalization is also ongoing in amphipod mitogenomes. The novel mitogenomes of the two Antarctic species have features distinguishing them from other amphipod mitogenomes such as a lower AT-richness in the whole mitogenomes and a negative GC- skew in both strands of protein coding genes. The genetically most variable mitochondrial regions of amphipods are nad6 and atp8, while cox1 shows low nucleotide diversity among closely and more distantly related species. In comparison to the pancrustacean mitochondrial ground pattern, E. giganteus shows a translocation of the nad1 gene, while cytb and nad6 genes are translocated in C. amundseni. Phylogenetic analysis based on mitogenomes illustrates that Eusirus and Charcotia cluster together with other species belonging to the same amphipod superfamilies. In the absence of reference nuclear genomes, mitogenomes can be useful to develop markers for studying population genetics or evolutionary relationships at higher taxonomic levels.

The Mitochondrial Genome of a Freshwater Pelagic Amphipod Macrohectopus branickii Is among the Longest in Metazoa

There are more than 350 species of amphipods (Crustacea) in Lake Baikal, which have emerged predominantly through the course of endemic radiation. This group represents a remarkable model for studying various aspects of evolution, one of which is the evolution of mitochondrial (mt) genome architectures. We sequenced and assembled the mt genome of a pelagic Baikalian amphipod species Macrohectopus branickii. The mt genome is revealed to have an extraordinary length (42,256 bp), deviating significantly from the genomes of other amphipod species and the majority of animals. The mt genome of M. branickii has a unique gene order within amphipods, duplications of the four tRNA genes and Cox2, and a long non-coding region, that makes up about two thirds of the genome’s size. The extension of the mt genome was most likely caused by multiple duplications and inversions of regions harboring ribosomal RNA genes. In this study, we analyzed the patterns of mt genome length changes in amphipods and other animal phyla. Through a statistical analysis, we demonstrated that the variability in the mt genome length may be a characteristic of certain phyla and is primarily conferred by expansions of non-coding regions.

Mitochondrial Genomes, Phylogenetic Associations, and SNP Recovery for the Key Invasive Ponto-Caspian Amphipods in Europe

The Ponto-Caspian region is the main donor of invasive amphipods to freshwater ecosystems, with at least 13 species successfully established in European inland waters. Dikerogammarus spp. and Pontogammarus robustoides are among the most successful, due to their strong invasive impact on local biota. However, genomic knowledge about these invaders is scarce, while phylogeography and population genetics have been based on short fragments of mitochondrial markers or nuclear microsatellites. In this study, we provide: (i) a reconstruction of six mitogenomes for four invasive gammarids (D. villosus, D. haemobaphes, D. bispinosus, and P. robustoides); (ii) a comparison between the structure of the newly obtained mitogenomes and those from the literature; (iii) SNP calling rates for individual D. villosus and D. haemobaphes from different invasion sites across Europe; and (iv) the first time-calibrated full mitogenome phylogeny reconstruction of several Ponto-Caspian taxa. We found that, in comparison to other gammarids, the mitogenomes of Ponto-Caspian species show a translocation between the tRNA-E and tRNA-R positions. Phylogenetic reconstruction using the mitogenomes identified that Ponto-Caspian gammarids form a well-supported group that originated in the Miocene. Our study supports paraphyly in the family Gammaridae. These provided mitogenomes will serve as vital genetic resources for the development of new markers for PCR-based identification methods and demographic studies.

Evolutionary History of Inversions in Directional Mutational Pressures in Crustacean Mitochondrial Genomes: Implications for Evolutionary Studies

Inversions of the origin of replication (ORI) in mitochondrial genomes produce asymmetrical mutational pressures that can cause strong base composition skews. Due to skews often being overlooked, the total number of crustacean lineages that underwent ORI events remains unknown. We analysed skews, cumulative skew plots, conserved sequence motifs, and mitochondrial architecture of all 965 available crustacean mitogenomes (699 unique species). We found indications of an ORI in 159 (22.7%) species, and mapped these to 23 ORI events: 16 identified with confidence and 7 putative (13 newly proposed, and for 5 we improved the resolution). Two ORIs occurred at or above the order level: Isopoda and Copepoda. Shifts in skew plots are not a precise tool for identifying the replication mechanism. We discuss how ORIs can produce mutational bursts in mitogenomes and show how these can interfere with various types of evolutionary studies. Phylogenetic analyses were plagued by artefactual clustering, and ORI lineages exhibited longer branches, a higher number of synonymous substitutions, higher mutational saturation, and higher compositional heterogeneity. ORI events also affected codon usage and protein properties. We discuss how this may have caused erroneous interpretation of data in previous studies that did not account for skew patterns.

Low annual temperature likely prevents the Holarctic amphipod Gammarus lacustris from invading Lake Baikal

Species with effective thermal adaptation mechanisms allowing them to thrive within a wide temperature range can benefit from climatic changes as they can displace highly specialized species. Here, we studied the adaptive capabilities of the Baikal endemic amphipods Eulimnogammarus verrucosus (Gerstfeld, 1858) and Eulimnogammarus cyaneus (Dybowsky, 1874) compared to the potential Holarctic Baikal invader Gammarus lacustris Sars, 1863 at the cellular level including the energy metabolism and the antioxidant system. All species were long-term exposed to a range of temperatures (1.5 °C to mimic winter conditions and the three species-specific preferred temperatures (i.e., 6 °C for E. verrucosus, 12 °C for E. cyaneus and 15 °C for G. lacustris). At 1.5 °C, we found species-specific metabolic alterations (i.e., significantly reduced ATP content and lactate dehydrogenase activity) indicating limitations on the activity level in the Holarctic G. lacustris. Although the two Baikal endemic amphipod species largely differ in thermal tolerance, no such limitations were found at 1.5 °C. However, the cold-stenothermal Baikal endemic E. verrucosus showed changes indicating a higher involvement of anaerobic metabolism at 12 °C and 15 °C, while the metabolic responses of the more eurythermal Baikal endemic E. cyaneus may support aerobic metabolism and an active lifestyle at all exposure temperatures. Rising temperatures in summer may provide a competitive advantage for G. lacustris compared to the Baikal species but the inactive lifestyle in the cold is likely preventing G. lacustris from establishing a stable population in Lake Baikal.

The mitochondrial genomes of five spring and groundwater amphipods of the family Crangonyctidae (Crustacea: Amphipoda) from eastern North America

We sequenced the mitochondrial genomes of one spring-dwelling (Crangonyx forbesi) and four groundwater amphipods (Bactrurus brachycaudus, Stygobromus allegheniensis, S. pizzinii, and S. t. potomacus) from eastern North America using a shotgun sequencing approach on an Illumina HiSeq 4000 (Illumina, San Diego, CA). All five mitochondrial genomes encoded 13 protein-coding genes, 22 transfer RNAs (tRNAs), and two ribosomal RNAs (rRNAs) representative of subphylum Crustacea. Although the four groundwater species exhibited gene orders nearly identical to the ancestral pancrustacean gene order, the spring-dwelling species, C. forbesi, possessed a transposition of the trnH–nad4–nad4l loci downstream after nad6–cytb–trnS2. Moreover, a long nad5 locus, longer rrnL, and rrnS loci, and unconventional start codons distinguished C. forbesi from the four groundwater amphipods. Overall, our five amphipod mitogenomes add to the increasing publicly available mitogenome resources for amphipods that are not only valuable for studying the evolutionary relationships of this diverse group of crustaceans but for exploring the evolution of mitochondrial genomes in general.

Comparative Mitogenomics in Hyalella (Amphipoda: Crustacea)

We present the sequencing and comparative analysis of 17 mitochondrial genomes of Nearctic and Neotropical amphipods of the genus Hyalella, most from the Andean Altiplano. The mitogenomes obtained comprised the usual 37 gene-set of the metazoan mitochondrial genome showing a gene rearrangement (a reverse transposition and a reversal) between the North and South American Hyalella mitogenomes. Hyalella mitochondrial genomes show the typical AT-richness and strong nucleotide bias among codon sites and strands of pancrustaceans. Protein-coding sequences are biased towards AT-rich codons, with a preference for leucine and serine amino acids. Numerous base changes (539) were found in tRNA stems, with 103 classified as fully compensatory, 253 hemi-compensatory and the remaining base mismatches and indels. Most compensatory Watson–Crick switches were AU -> GC linked in the same haplotype, whereas most hemi-compensatory changes resulted in wobble GU and a few AC pairs. These results suggest a pairing fitness increase in tRNAs after crossing low fitness valleys. Branch-site level models detected positive selection for several amino acid positions in up to eight mitochondrial genes, with atp6 and nad5 as the genes displaying more sites under selection.

First draft genome for the sand-hopper Trinorchestia longiramus

Crustacean amphipods are important trophic links between primary producers and higher consumers. Although most amphipods occur in or around aquatic environments, the family Talitridae is the only family found in terrestrial and semi-terrestrial habitats. The sand-hopper Trinorchestia longiramus is a talitrid species often found in the sandy beaches of South Korea. In this study, we present the first draft genome assembly and annotation of this species. We generated ~380.3 Gb of sequencing data assembled in a 0.89 Gb draft genome. Annotation analysis estimated 26,080 protein-coding genes, with 89.9% genome completeness. Comparison with other amphipods showed that T. longiramus has 327 unique orthologous gene clusters, many of which are expanded gene families responsible for cellular transport of toxic substances, homeostatic processes, and ionic and osmotic stress tolerance. This first talitrid genome will be useful for further understanding the mechanisms of adaptation in terrestrial environments, the effects of heavy metal toxicity, as well as for studies of comparative genomic variation across amphipods.

Architectural instability, inverted skews and mitochondrial phylogenomics of Isopoda: outgroup choice affects the long-branch attraction artefacts

The majority strand of mitochondrial genomes of crustaceans usually exhibits negative GC skews. Most isopods exhibit an inversed strand asymmetry, believed to be a consequence of an inversion of the replication origin (ROI). Recently, we proposed that an additional ROI event in the common ancestor of Cymothoidae and Corallanidae families resulted in a double-inverted skew (negative GC), and that taxa with homoplastic skews cluster together in phylogenetic analyses (long-branch attraction, LBA). Herein, we further explore these hypotheses, for which we sequenced the mitogenome of Asotana magnifica (Cymothoidae), and tested whether our conclusions were biased by poor taxon sampling and inclusion of outgroups. (1) The new mitogenome also exhibits a double-inverted skew, which supports the hypothesis of an additional ROI event in the common ancestor of Cymothoidae and Corallanidae families. (2) It exhibits a unique gene order, which corroborates that isopods possess exceptionally destabilized mitogenomic architecture. (3) Improved taxonomic sampling failed to resolve skew-driven phylogenetic artefacts. (4) The use of a single outgroup exacerbated the LBA, whereas both the use of a large number of outgroups and complete exclusion of outgroups ameliorated it.

Disentangling the interplay of positive and negative selection forces that shaped mitochondrial genomes of Gammarus pisinnus and Gammarus lacustris

We hypothesized that the mitogenome of Gammarus lacustris (GL), native to the Qinghai-Tibet Plateau, might exhibit genetic adaptations to the extreme environmental conditions associated with high altitudes (greater than 3000 m). To test this, we also sequenced the mitogenome of Gammarus pisinnus (GP), whose native range is close to the Tibetan plateau, but at a much lower altitude (200-1500 m). The two mitogenomes exhibited conserved mitochondrial architecture, but low identity between genes (55% atp8 to 76.1% cox1). Standard (homogeneous) phylogenetic models resolved Gammaridae as paraphyletic, but 'heterogeneous' CAT-GTR model as monophyletic. In indirect support of our working hypothesis, GL, GP and Gammarus fossarum exhibit evidence of episodic diversifying selection within the studied Gammaroidea dataset. The mitogenome of GL generally evolves under a strong purifying selection, whereas GP evolves under directional (especially pronounced in atp8) and/or relaxed selection. This is surprising, as GP does not inhabit a unique ecological niche compared to other gammarids. We propose that this rapid evolution of the GP mitogenome may be a reflection of its relatively recent speciation and heightened non-adaptive (putatively metabolic rate-driven) mutational pressures. To test these hypotheses, we urge sequencing mitogenomes of remaining Gammarus species populating the same geographical range as GP.

Hidden cases of tRNA gene duplication and remolding in mitochondrial genomes of amphipods

The evolution of tRNA genes in mitochondrial (mt) genomes is a complex process that includes duplications, degenerations, and transpositions, as well as a specific process of identity change through mutations in the anticodon (tRNA gene remolding or tRNA gene recruitment). Using amphipod-specific tRNA models for annotation, we show that tRNA duplications are more common in the mt genomes of amphipods than what was revealed by previous annotations. Seventeen cases of tRNA gene duplications were detected in the mt genomes of amphipods, and ten of them were tRNA genes that underwent remolding. The additional tRNA gene findings were verified using phylogenetic analysis and genetic distance analysis. The majority of remolded tRNA genes (seven out of ten cases) were found in the mt genomes of endemic amphipod species from Lake Baikal. All additional mt tRNA genes arose independently in the Baikalian amphipods, indicating the unusual plasticity of tRNA gene evolution in these species assemblages. The possible reasons for the unusual abundance of additional tRNA genes in the mt genomes of Baikalian amphipods are discussed. The amphipod-specific tRNA models developed for MiTFi refine existing predictions of tRNA genes in amphipods and reveal additional cases of duplicated tRNA genes overlooked by using less specific Metazoa-wide models. The application of these models for mt tRNA gene prediction will be useful for the correct annotation of mt genomes of amphipods and probably other crustaceans.

Symbiosis, Selection, and Novelty: Freshwater Adaptation in the Unique Sponges of Lake Baikal

Freshwater sponges (Spongillida) are a unique lineage of demosponges that secondarily colonized lakes and rivers and are now found ubiquitously in these ecosystems. They developed specific adaptations to freshwater systems, including the ability to survive extreme thermal ranges, long-lasting dessication, anoxia, and resistance to a variety of pollutants. Although spongillids have colonized all freshwater systems, the family Lubomirskiidae is endemic to Lake Baikal and plays a range of key roles in this ecosystem. Our work compares the genomic content and microbiome of individuals of three species of the Lubomirskiidae, providing hypotheses for how molecular evolution has allowed them to adapt to their unique environments. We have sequenced deep (>92% of the metazoan "Benchmarking Universal Single-Copy Orthologs" [BUSCO] set) transcriptomes from three species of Lubomirskiidae and a draft genome resource for Lubomirskia baikalensis. We note Baikal sponges contain unicellular algal and bacterial symbionts, as well as the dinoflagellate Gyrodinium. We investigated molecular evolution, gene duplication, and novelty in freshwater sponges compared with marine lineages. Sixty one orthogroups have consilient evidence of positive selection. Transporters (e.g., zinc transporter-2), transcription factors (aristaless-related homeobox), and structural proteins (e.g. actin-3), alongside other genes, are under strong evolutionary pressure in freshwater, with duplication driving novelty across the Spongillida, but especially in the Lubomirskiidae. This addition to knowledge of freshwater sponge genetics provides a range of tools for understanding the molecular biology and, in the future, the ecology (e.g., colonization and migration patterns) of these key species.

The mitochondrial genome of UK (non-native) Dikerogammarus haemobaphes (Amphipoda: Gammaridae) informs upon Dikerogammarus evolution, invasions and associated microparasites

The amphipod Dikerogammarus haemobaphes is a high-risk carrier of parasites that impact wildlife in its non-native range. Studies using the mitochondrial genes, Cytochrome Oxidase Sub-Unit 1 (cox1) and small-subunit ribosomal RNA gene (16S), provide some nucleotide detail for understanding the evolution and phylogeography of this species. Despite this, the origins of the invasion remain unknown, as do the origins of its parasites. This study provides the full annotated mitochondrial genome (15,460 bp) of D. haemobaphes, consisting of 2 rRNAs, 24 tRNAs and 14 protein coding genes. Mitochondrial genes from the UK isolate are compared to existing data on NCBI and are used in a concatenated phylogenetic approach and identify D. haemobaphes as an early member of the Gammaridae (Amphipoda). Viral, bacterial, protistan and microsporidian parasites are present across the Gammaridae, including D. haemobaphes, suggesting the ancestor of the Gammaridae harboured related diseases, and that further screening of amphipods is likely to reveal further microparasite diversity. This correlation suggests that other gammarid invaders have the potential to harbour a range of microparasites. The mitochondrial genome of this species will act a resource to facilitate our understanding of geneflow, disease epidemiology and evolutionary history in this invasion-disease model.

Indication of ongoing amphipod speciation in Lake Baikal by genetic structures within endemic species

Background

The ancient Lake Baikal is characterized by an outstanding diversity of endemic faunas with more than 350 amphipod species and subspecies. We determined the genetic diversity within the endemic littoral amphipod species Eulimnogammarus verrucosus, E. cyaneus and E. vittatus and investigated whether within those species genetically separate populations occur across Lake Baikal. Gammarus lacustris from water bodies in the Baikal area was examined for comparison.

Results

Genetic diversities within a species were determined based on fragments of cytochrome c oxidase I (COI) and for E. verrucosus additionally of 18S rDNA. Highly location-specific haplogroups of E. verrucosus and E. vittatus were found at the southern and western shores of Baikal that are separated by the Angara River outflow; E. verrucosus from the eastern shore formed a further, clearly distinct haplotype cluster possibly confined by the Selenga River and Angarskiy Sor deltas. The genetic diversities within these haplogroups were lower than between the different haplogroups. Intraspecific genetic diversities within E. verrucosus and E. vittatus with 13 and 10%, respectively, were similar to interspecies differences indicating the occurrence of cryptic, morphologically highly similar species; for E. verrucosus this was confirmed with 18S rDNA. The haplotypes of E. cyaneus and G. lacustris specimens were with intraspecific genetic distances of 3 and 2%, respectively, more homogeneous indicating no or only recent disruption of gene flow of E. cyaneus across Baikal and recent colonization of water bodies around Baikal by G. lacustris.

Conclusions

Our finding of separation of subgroups of Baikal endemic amphipods to different degrees points to a species-specific ability of dispersal across areas with adverse conditions and to potential geographical dispersal barriers in Lake Baikal.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1470-8) contains supplementary material, which is available to authorized users.

The complete mitochondrial genome of the sand-hopper Trinorchestia longiramus (Amphipoda: Talitridae)

The complete mitochondrial genome of sand-hopper Trinorchestia longiramus was analyzed in this study, which is the first for the genus within the family Talitridae. The mitogenome sequence is 15,401 bp in length containing two ribosomal RNA genes, 22 transfer RNA genes, 13 protein-coding genes, and a control region as found in most amphipods. The gene order showed that T. longiramus has a unique control region location compared to other amphipods. Phylogenetic analysis using the maximum likelihood method positioned T. longiramus within the monophyletic clades of the family Talitridae.

Bursts of amino acid replacements in protein evolution

Evolution can occur both gradually and through alternating episodes of stasis and rapid changes. However, the prevalence and magnitude of fluctuations of the rate of evolution remain obscure. Detecting a rapid burst of changes requires a detailed record of past evolution, so that events that occurred within a short time interval can be identified. Here, we use the phylogenies of the Baikal Lake amphipods and of Catarrhini, which contain very short internal edges which make this task feasible. We detect six salient bursts of evolution of individual proteins during such short time periods, each involving between six and 38 amino acid substitutions. These bursts were extremely unlikely to have occurred neutrally, and were apparently caused by positive selection. On average, in the course of a time interval required for one synonymous substitution per site, a protein undergoes a strong burst of rapid evolution with probability at least approximately 0.01.

Characterization of the mitochondrial genome of an ancient amphipod Halice sp. MT-2017 (Pardaliscidae) from 10,908 m in the Mariana Trench

Small amphipods (Halice sp. MT-2017) with body length <1 cm were collected from the Challenger Deep (~10,920 m below sea level). The divergence time of their lineage was approximately 109 Mya, making this group ancient compared to others under study. The mitochondrial genome of Halice sp. shared the usual gene components of metazoans, comprising 13 protein coding genes (PCGs), 22 transfer RNAs (tRNAs), and 2 ribosomal RNAs (rRNAs). The arrangement of these genes, however, differed greatly from that of other amphipods. Of the 15 genes that were rearranged with respect to the pancrustacean gene pattern, 12 genes (2 PCGs, 2 rRNAs, and 8 tRNAs) were both translocated and strand-reversed. In contrast, the mitochondrial genomes in other amphipods never show so many reordered genes, and in most instances, only tRNAs were involved in strand-reversion-coupled translocation. Other characteristics, including reversed strand nucleotide composition bias, relatively higher composition of non-polar amino acids, and lower evolutionary rate, were also identified. Interestingly, the latter two features were shared with another hadal amphipod, Hirondellea gigas, suggesting their possible associations with the adaptation to deep-sea extreme habitats. Overall, our data provided a useful resource for future studies on the evolutionary and adaptive mechanisms of hadal faunas.

Cryptic species in a well-known habitat: applying taxonomics to the amphipod genus Epimeria (Crustacea, Peracarida)

Taxonomy plays a central role in biological sciences. It provides a communication system for scientists as it aims to enable correct identification of the studied organisms. As a consequence, species descriptions should seek to include as much available information as possible at species level to follow an integrative concept of 'taxonomics'. Here, we describe the cryptic species Epimeria frankei sp. nov. from the North Sea, and also redescribe its sister species, Epimeria cornigera. The morphological information obtained is substantiated by DNA barcodes and complete nuclear 18S rRNA gene sequences. In addition, we provide, for the first time, full mitochondrial genome data as part of a metazoan species description for a holotype, as well as the neotype. This study represents the first successful implementation of the recently proposed concept of taxonomics, using data from high-throughput technologies for integrative taxonomic studies, allowing the highest level of confidence for both biodiversity and ecological research.