Species delimitation and mitogenome phylogenetics in the subterranean genus Pseudoniphargus (Crustacea: Amphipoda)

Mitochondrial introgression and mito-nuclear discordance obscured the closely related species boundaries in Cletus Stål from China (Heteroptera: Coreidae)

Accurate taxonomy and delimitation are of great importance for pest control strategies and management programs. Here, we focus on Cletus (Insecta: Hemiptera: Coreidae), which includes many crop pests. The species boundaries still conflict and only cytochrome c oxidase subunit I (COI) barcoding has been previously used for molecular studies. We generated new mitochondrial genome and nuclear genome-wide SNPs to explore the species boundaries of 46 Cletus samples from China using multiple species delimitation approaches. All results recovered a monophyly with high support, except for two closely related species in clade I - C. punctiger and C. graminis. Mitochondrial data demonstrated admixture in clade I, while genome-wide SNPs unambiguously identified two separate species, which were confirmed by morphological classification. Inconsistent nuclear and mitochondrial data indicated mito-nuclear discordance. Mitochondrial introgression is the most likely explanation, and more extensive sampling and more comprehensive data are needed to ascertain a pattern. Accurate species delimitation will shed light on species status; thus, an accurate taxonomy is of particular concern, as there is a pressing need to implement precise control of agricultural pests and to perform further research on diversification.

Multilocus phylogeny and species delimitation suggest synonymies of two Lucanus Scopoli, 1763 (Coleoptera, Lucanidae) species names

Phylogenetic relationsships of four nominal Lucanus Scopoli, 1763 species, L.swinhoei Parry, 1874, L.continentalis Zilioli, 1998, L.liuyei Huang & Chen, 2010, and L.wuyishanensis Schenk, 1999, are assessed based on mitochondrial (16S rDNA, COI) and nuclear (28S rDNA, Wingless) genes. The genetic distance is 0.0072 between L.swinhoei and L.continentalis, and 0.0094 between L.wuyishanensis and L.liuyei. Three species-delimitation approaches (ABGD, PTP, and GMYC) consistently showed L.swinhoei + L.continentalis and L.wuyishanensis + L.liuyei as two MOTUs. A new synonymy, L.liuyei = L.wuyishanensis, is proposed. Synonymy of L.swinhoei over L.continentalis is confirmed.

Bermuda’s Walsingham Caves: A Global Hotspot for Anchialine Stygobionts

Bermuda is an Eocene age volcanic island in the western North Atlantic, entirely capped by Pleistocene eolian limestone. The oldest and most highly karstified limestone is a 2 km2 outcrop of the Walsingham Formation containing most of the island’s 150+ caves. Extensive networks of submerged cave passageways, flooded by saltwater, extend under the island. In the early 1980s, cave divers initially discovered an exceptionally rich and diverse anchialine community inhabiting deeper sections of the caves. The fauna inhabiting caves in the Walsingham Tract consists of 78 described species of cave-dwelling invertebrates, including 63 stygobionts and 15 stygophiles. Thus, it represents one of the world’s top hotspots of subterranean biodiversity. Of the anchialine fauna, 65 of the 78 species are endemic to Bermuda, while 66 of the 78 are crustaceans. The majority of the cave species are limited in their distribution to just one or only a few adjacent caves. Due to Bermuda’s high population density, water pollution, construction, limestone quarries, and trash dumping produce severe pressures on cave fauna and groundwater health. Consequently, the IUCN Red List includes 25 of Bermuda’s stygobiont species as critically endangered.

Evidence of cryptic diversity in freshwater Macrobrachium prawns from Indochinese riverine systems revealed by DNA barcode, species delimitation and phylogenetic approaches

The diversity of Indochinese prawns in genus Macrobrachium is enormous due to the habitat diversification and broad tributary networks of two river basins: the Chao Phraya and the Mekong. Despite long-standing interest in SE-Asian decapod diversity, the subregional Macrobrachium fauna is still not yet comprehensively clarified in terms of taxonomic identification or genetic diversification. In this study, integrative taxonomic approaches including morphological examination, DNA barcoding, and molecular species delimitation were used to emphasize the broad scale systematics of Macrobrachium prawns in Indochina. Twenty-seven nominal species were successfully re-verified by traditional and molecular taxonomy. Barcode gap analysis supported broad overlapping of species boundaries. Taxonomic ambiguity of several deposited samples in the public database is related to inter- and intraspecific genetic divergence as indicated by BOLD discordance. Diagnostic nucleotide positions were found in six Macrobrachium species. Eighteen additional putative lineages are herein assigned using the consensus of species delimitation methods. Genetic divergence indicates the possible existence of cryptic species in four morphologically complex and wide-ranging species: M. lanchesteri, M. niphanae, M. sintangense, and some members of the M. pilimanus group. The geographical distribution of some species supports the connections and barriers attributed to paleo-historical events of SE-Asian rivers and land masses. Results of this study show explicitly the importance of freshwater ecosystems in Indochinese subregions, especially for the Mekong River Basin due to its high genetic diversity and species composition found throughout its tributaries.

First complete mitogenomes of Diamesinae, Orthocladiinae, Prodiamesinae, Tanypodinae (Diptera: Chironomidae) and their implication in phylogenetics

BACKGROUND

The mitochondrial genome (mitogenome) has been extensively used for phylogenetic and evolutionary analysis in Diptera, but the study of mitogenome is still scarce in the family Chironomidae.

METHODS

Here, the first complete mitochondrial genomes of four Chironomid species representing Diamesinae, Orthocladiinae, Prodiamesinae and Tanypodinae are presented. Coupled with published mitogenomes of two, a comparative mitochondrial genomic analysis between six subfamilies of Chironomidae was carried out.

RESULTS

Mitogenomes of Chironomidae are conserved in structure, each contains 37 typical genes and a control region, and all genes arrange the same gene order as the ancestral insect mitogenome. Nucleotide composition is highly biased, the control region displayed the highest A + T content. All protein coding genes are under purifying selection, and the ATP8 evolves at the fastest rate. In addition, the phylogenetic analysis covering six subfamilies within Chironomidae was conducted. The monophyly of Chironomidae is strongly supported. However, the topology of six subfamilies based on mitogenomes in this study is inconsistent with previous morphological and molecular studies. This may be due to the high mutation rate of the mitochondrial genetic markers within Chironomidae. Our results indicate that mitogenomes showed poor signals in phylogenetic reconstructions at the subfamily level of Chironomidae.

A multilocus assessment reveals two new synonymies for East Asian Cyclommatus stag beetles (Coleoptera, Lucanidae)

Cyclommatus scutellaris Möllenkamp, 1912, Cyclommatus elsae Kriesche, 1921 and Cyclommatus tamdaoensis Fujita, 2010 are East Asian stag beetle species with long-debated taxonomic relationships due to high intraspecific morphological variability. In this study, we applied multilocus phylogenetic analyses to reassess their relationships. Two mitochondrial genes (16S rDNA, COI) and two nuclear genes (28S rDNA, Wingless) were used to reconstruct the phylogeny through the Bayesian inference (BI) and Maximum Likelihood (ML) methods. Both topologies supported two clades: the clade C. scutellaris was sister to the clade (C. elsae + C. tamdaoensis) with the subclade C. tamdaoensis embedded in the subclade C. elsae. The Kimura 2-parameter (K2P) genetic distance analysis yielded a low mean value (≤0.035) among the three taxa, which was well below the minimum mean value between other Cyclommatus species (≥0.122). We also compared the accuracy and efficiency of two approaches, GMYC and ABGD, in delimitating the three lineages. The result shows that ABGD is a better approach than GMYC. Our molecular data recognizes the three species as different populations of a single species, ranging from Taiwan Island to the continent. Therefore, we propose two new junior synonyms for C. scutellaris: C. tamdaoensis, syn. nov. and C. elsae syn. nov.

HERMES: An improved method to test mitochondrial genome molecular synapomorphies among clades

Mitochondrial chromosomes have diversified among eukaryotes and many different architectures and features are now acknowledged for this genome. Here we present the improved HERMES index, which can measure and quantify the amount of molecular change experienced by mitochondrial genomes. We test the improved approach with ten molecular phylogenetic studies based on complete mitochondrial genomes, representing six bilaterian Phyla. In most cases, HERMES analysis spotted out clades or single species with peculiar molecular synapomorphies, allowing to identify phylogenetic and ecological patterns. The software presented herein handles linear, circular, and multi-chromosome genomes, thus widening the HERMES scope to the complete eukaryotic domain.

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.

Groundwater annelids from Gran Canaria and Fuerteventura (Canary Islands), with the description of two new species of Namanereis (Namanereidinae, Nereididae, Polychaeta)

The Canary Islands are the richest volcanic region in the world in subterranean adapted fauna, followed by the Hawaiian Islands and the Undara Cave in Australia. Most of the subterranean adapted aquatic fauna from the Canary Islands is restricted to the anchialine environments in La Corona lava tube in Lanzarote, while the oligohaline stygobiont fauna, usually found in groundwater or interstitial freshwaters, is scarcer and represented by a few species of amphipods, copepods, and a single polychaete annelid recorded from Fuerteventura and doubtfully identified as Namanereis hummelincki (non Augener, 1933) (Hartmann-Schröder, 1988). Two new species of polychaete annelids belonging to the subfamily Namanereidinae are described from Gran Canaria and Fuerteventura. Both species live in groundwater, are depigmented and eyeless, and have bifid jaws. Although they are seemingly more related to each other than to other members of the bifid-jaw group, Namanereis canariarumsp. nov. can be diagnosed by its relatively longer tentacular cirri and parapodial dorsal cirri, as well as the presence of pseudospiniger chaetae. In contrast, Namanereis llanetensissp. nov. has shorter cirri and usually lacks pseudospiniger chaetae. Namanereis canariarumsp. nov. and Namanereis llanetensissp. nov. increase to 20 the total number of currently described species within this enigmatic genus. More than half of those species are adapted to live in groundwaters.

Mitochondrial genomes of twelve species of hyperdiverse Trigonopterus weevils

Mitochondrial genomes of twelve species of Trigonopterus weevils are presented, ten of them complete. We describe their gene order and molecular features and test their potential for reconstructing the phylogeny of this hyperdiverse genus comprising > 1,000 species. The complete mitochondrial genomes examined herein ranged from 16,501 bp to 21,007 bp in length, with an average AT content of 64.2% to 69.7%. Composition frequencies and skews were generally lower across species for atp6, cox1-3, and cob genes, while atp8 and genes coded on the minus strand showed much higher divergence at both nucleotide and amino acid levels. Most variation within genes was found at the codon level with high variation at third codon sites across species, and with lesser degree at the coding strand level. Two large non-coding regions were found, CR1 (between rrnS and trnI genes) and CR2 (between trnI and trnQ), but both with large variability in length; this peculiar structure of the non-coding region may be a derived character of Curculionoidea. The nad1 and cob genes exhibited an unusually high interspecific length variation of up to 24 bp near the 3' end. This pattern was probably caused by a single evolutionary event since both genes are only separated by trnS2 and length variation is extremely rare in mitochondrial protein coding genes. We inferred phylogenetic trees using protein coding gene sequences implementing both maximum likelihood and Bayesian approaches, each for both nucleotide and amino acid sequences. While some clades could be retrieved from all reconstructions with high confidence, there were also a number of differences and relatively low support for some basal nodes. The best partition scheme of the 13 protein coding sequences obtained by IQTREE suggested that phylogenetic signal is more accurate by splitting sequence variation at the codon site level as well as coding strand, rather than at the gene level. This result corroborated the different patterns found in Trigonopterus regarding to A+T frequencies and AT and GC skews that also greatly diverge at the codon site and coding strand levels.

Mitogenome phylogenetics in the genus Palaemon (Crustacea: Decapoda) sheds light on species crypticism in the rockpool shrimp P. elegans

The genus Palaemon comprises worldwide marine and freshwater shrimps and prawns, and some of them are ecologically or commercially important species. Palaemon is not currently a monophyletic group, so phylogenetics and systematics are constantly changing. Species crypticism has been pointed out in several Palaemon species, being the clearest evidence in the European rockpool shrimp P. elegans. Here we sequenced and described seven European Palaemon mitochondrial genomes. The mitochondrial protein-coding genes were used, along with those of three other Palaemon species, to perform mitogenome phylogenetic analyses to clarify the evolutionary relationships within the genus, and particularly to shed light on the cryptic species found within P. elegans. The Messinian Salinity Crisis (5.3-5.9 Ma, late Miocene) was proposed to be the origin of this cryptic species and it was used as aged constraint for calibration analysis. We provide the largest and the first time-calibrated mitogenome phylogeny of the genus Palaemon and mitogenome substitution rate was estimated (1.59% per million years) in Decapoda for the first time. Our results highlighted the need for future systematics changes in Palaemon and crypticism in P. elegans was confirmed. Mitochondrial genome and cox1 (1.41%) substitution rate estimates matched those published elsewhere, arguing that the Messinian Salinity Crisis was a plausible event driving the split between P. elegans and its cryptic species. Molecular dating suggested that Pleistocene glaciations were likely involved in the differentiation between the Atlantic and Mediterranean populations of P. elegans. On the contrary, the divergence between the Atlantic and Mediterranean populations of the common littoral shrimp P. serratus was greater and dated to be much older (4.5-12.3 Ma, Plio-Miocene), so we considered that they could represent two separated species. Therefore, species crypticism in the genus Palaemon seems to be a common phenomenon.

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.