Hidden diversity of Ctenophora revealed by new mitochondrial COI primers and sequences

The mitochondrial gene cytochrome-c-oxidase subunit 1 (COI) is useful in many taxa for phylogenetics, population genetics, metabarcoding, and rapid species identifications. However, the phylum Ctenophora (comb jellies) has historically been difficult to study due to divergent mitochondrial sequences and the corresponding inability to amplify COI with degenerate and standard COI “barcoding” primers. As a result, there are very few COI sequences available for ctenophores, despite over 200 described species in the phylum. Here, we designed new primers and amplified the COI fragment from members of all major groups of ctenophores, including many undescribed species. Phylogenetic analyses of the resulting COI sequences revealed high diversity within many groups that was not evident from more conserved 18S rDNA sequences, in particular among the Lobata (Ctenophora; Tentaculata; Lobata). The COI phylogenetic results also revealed unexpected community structure within the genus Bolinopsis, suggested new species within the genus Bathocyroe, and supported the ecological and morphological differences of some species such as Lampocteis cruentiventer and similar undescribed lobates (Lampocteis sp. “V” stratified by depth, and “A” differentiated by colour). The newly designed primers reported herein provide important tools to enable researchers to illuminate the diversity of ctenophores worldwide via quick molecular identifications, improve the ability to analyse environmental DNA by improving reference libraries and amplifications, and enable a new breadth of population genetic studies.     

Single nucleotide polymorphism-based species phylogeny of greater fritillary butterflies (Lepidoptera: Nymphalidae: Speyeria) demonstrates widespread mitonuclear discordance

The systematics of Speyeria butterflies has historically been complicated by intraspecific variability that has challenged efforts to delimit species and reconstruct phylogenies. Our study presents a phylogenetic comparison of genomic single nucleotide polymorphisms (SNPs) and mitochondrial COI gene sequences, with comprehensive taxon sampling that includes 15 species and 46 subspecies. Increased sampling of genetic markers and taxa improved the match between genetic clusters, obtained with both phylogenetic and cluster-based analyses, and species previously detected using morphology, as well as showing two species delimitations that may need revision. We also recovered extensive mitonuclear discordance between genomic SNPs and the COI gene, confirming that mitochondrial DNA does not reliably identify several species at broad geographic scales. Resolution of the relationships of Speyeria species demonstrates the importance of sampling variation across the whole genome, and provides an essential foundation for understanding the evolution of this charismatic clade of North American butterflies.     

The subgenus Ashima (Diptera,Drosophildae, Phortica) from China,with DNA barcoding and descriptions of three new species

In the present study, three new species found in Yunnan, Southwest China were described as members of the subgenus Phortica (Ashima): P. (A.) haba An & Chen, sp. nov., P. (A.) montipagana An & Chen, sp. nov. and P. (A.) qingsongi An & Chen, sp. nov. Barcode sequences (partial sequences of the mitochondrial COI gene) were collected from 61 specimens of 16 known and the above-mentioned three new Ashima species. The intra- and interspecific pairwise K-2P (Kimura two-parameter) COI distances were analysed and a phylogenetic tree was constructed based on the barcode sequences. Species delimitation in this subgenus was supported by integrating barcodes with morphological information, in particular for the three new species, considered to be cryptic species. In addition, the diversification of lineages in the subgenus Ashima was proposed to occur in southern China and adjacent areas, suggesting specific adaptation of Ashima species to the high plateau environments.

http://zoobank.org/urn:lsid:zoobank.org:pub:69A81D0E-5993-426D-9B59-A2AC5E52BD84     

Phylogeny and taxonomic reassessment of pikas Ochotona pallasii and O. argentata (Mammalia,Lagomorpha)

We examined a cranial morphometric data set consisting of 186 specimens from the entire distribution range of Ochotona pallasii sensu lato and O. argentata, as well as 67 complete sequences of the COI gene and 28 sequences of the MGF and PRKCI introns from these and closely allied species. Our results show that the two allopatric morphologically similar taxa composing O. pallasii sensu lato – from Mongolia and adjacent territories and Kazakhstan – are paraphyletic relative to O. argentata. Genetic distances between these three taxa are larger than the intraspecific variation known for the subgenus Pika, in which the species under consideration belong; these distances are even larger than the interspecific differences among closely related species such as O. hyperborea, O. mantchurica and O. hoffmanni. Thus, the three focal taxa are recognized here as distinct species. Inspection of the type specimen of O. pallasii indicated that this specimen was not collected in Kazakhstan, has previously been theorized. The most probable place of the holotype's origin is Russian south‐eastern Altai (Chuyskaya Steppe); whatever its exact origin, it definitively originates from the ‘Mongolian’ taxon. Based on this evidence, the senior synonym for the Kazakh pika is O. opaca Argyropulo, 1930. Thus, we propose to recognize three separate species in the O. pallasii species group: O. pallasii (Mongolia and adjacent territories), O. opaca (eastern Kazakhstan) and O. argentata (Helan Shan Range, China).     

ITS2 secondary structure improves phylogeny estimation in a radiation of blue butterflies of the subgenus Agrodiaetus (Lepidoptera: Lycaenidae: Polyommatus)

Background  

Current molecular phylogenetic studies of Lepidoptera and most other arthropods are predominantly based on mitochondrial genes and a limited number of nuclear genes. The nuclear genes, however, generally do not provide sufficient information for young radiations. ITS2, which has proven to be an excellent nuclear marker for similarly aged radiations in other organisms like fungi and plants, is only rarely used for phylogeny estimation in arthropods, although universal primers exist. This is partly due to difficulties in the alignment of ITS2 sequences in more distant taxa. The present study uses ITS2 secondary structure information to elucidate the phylogeny of a species-rich young radiation of arthropods, the butterfly subgenus Agrodiaetus. One aim is to evaluate the efficiency of ITS2 to resolve the phylogeny of the subgenus in comparison with COI, the most important mitochondrial marker in arthropods. Furthermore, we assess the use of compensatory base changes in ITS2 for the delimitation of species and discuss the prospects of ITS2 as a nuclear marker for barcoding studies.     

Nuclear insertions and heteroplasmy of mitochondrial DNA as two sources of intra‐individual genomic variation in grasshoppers

We examined the level of intra‐individual variation in a region of the mitochondrial genome coding for cytochrome oxydase 1 (COI) in two grasshopper species using a clone‐and‐sequence analysis of hundreds of sequences. In both Locusta migratoria and Chortoicetes terminifera, we found that 60–65% of the clones were unique COI‐like sequences. Among these COI‐like sequences, 70–75% diverged by less than 1% from the real mitochondrial haplotypes, and were likely to represent microheteroplasmic molecules. About 20% of the COI‐like sequences diverged by more than 9% from the mitochondrial haplotypes, and generally included stop codons, suggesting that these sequences were nuclear mitochondrial pseudogenes (NUMTs). Only six sequences, diverging by 2–6% from the mitochondrial haplotypes, were identified as potentially misleading in phylogenetic studies. In addition, we found that five sequences from C. terminifera were associated with mobile elements or repetitive DNA families.     

Analysis of mitochondrial COI sequences of the Diachasmimorpha longicaudata (Hymenoptera: Braconidae) species complex in Thailand

Diachasmimorpha longicaudata is an Opiinae parasitoid used to control tephritid fruit flies, which cause tremendous economic losses of fruits worldwide. In Thailand, D. longicaudata is classified as three sibling species, DLA, DLB and DLBB, based on the morphological and biological species concepts but their genetic variation has not been studied. Therefore, we investigated the genetic differentiation of the mitochondrial COI gene to clarify the ambiguous taxonomy of this species complex. The 603‐bp COI region was sequenced from laboratory‐bred colonies and field‐collected specimens from seven locations representing five geographical regions in Thailand. DLA was associated with the host Bactrocera correcta while DLB and DLBB were associated with Bactrocera dorsalis. The interspecific nucleotide differences of COI sequences among the three groups ranged from 6.70% to 7.62% (Kimura 2‐parameter distance), which adequately separates species complexes within the order Hymenoptera and supports the current sibling species classification. The neighbor joining, maximum likelihood and consensus Bayesian phylogenetic trees constructed from COI sequences revealed that the three sibling species of laboratory and field‐collected D. longicaudata are monophyletic with 100% support. The high genetic variation and molecular phylogeny of the COI sequences were shown to discriminate between the D. longicaudata species examined in this study.     

Molecular systematics of threadfin breams and relatives (Teleostei,Nemipteridae)

Threadfin breams and relatives of the family Nemipteridae comprise 69 currently recognized species in five genera. They are found in the tropical and subtropical Indo‐West Pacific and most are commercially important. Using recently developed molecule‐based approaches exploiting DNA sequence variation among species/specimens, this study reconstructed a comprehensive phylogeny of the Nemipteridae, examined the validity of species and explored the cryptic diversity of the family, and tested previous phylogenetic hypotheses. A combined data set (105 taxa from 41 morphospecies) with newly determined sequences from two nuclear genes (RAG1 and RH) and one mitochondrial gene (COI), and a data set with only COI gene sequences (329 newly obtained plus 328 from public databases from a total of 53 morphospecies) were used in the phylogenetic analysis. The latter was further used for species delimitation analyses with two different tools to explore species diversity. Our phylogenetic results showed that all the currently recognized genera were monophyletic. The monotypic genus Scaevius is the sister group of Pentapodus and they together are sister to Nemipterus. These three genera combined to form the sister group of the clade comprising Parascolopsis and Scolopsis. The validity of most of the examined species was confirmed except in some cases. The combined evidence from the results of different analyses revealed a gap in our existing knowledge of species diversity in the Nemipteridae. We found several currently recognized species contain multiple separately evolving metapopulation lineages within species; some lineages should be considered as new species for further assignment. Finally, some problematic sequences deposited in public databases (probably due to misidentification) were also revised in this study to improve the accuracy for prospective DNA barcoding work on nemipterid fishes.     

First report on Wolbachia endosymbiosis in freshwater Aphelocheirus aestivalis (Heteroptera: Aphelocheiridae) and its potential impact on genetic diversity of host

Alphaproteobacteria Wolbachia have been described as endosymbionts of approximately half of all aquatic insect species. These bacteria might affect not only reproduction but also the genetic diversity of its hosts. In the present study we identified Wolbachia endosymbiosis in freshwater true bug Aphelocheirus aestivalis F., 1794 (Heteroptera: Aphelocheiridae). Despite the fact that A. aestivalis is widely distributed in Europe, it occurs rather locally, often in isolated populations. Taking into account that Wolbachia, close relationships and past demographic phenomena could affect the genetic diversity of its host, we analyzed mitochondrial (COI and 16S) and nuclear (internal transcribed spacer 2) markers determined for A. aestivalis individuals collected from five populations. Moreover, we compared obtained COI sequences with those deposited in GenBank. Analyses revealed low genetic differentiation among samples tested, as well as low variation among determined COI sequences and those downloaded from the database. Although Wolbachia infection could correlate with decreasing mitochondrial diversity of its host, we suggest that low genetic variation observed in tested A. aestivalis samples (at both mitochondrial and nuclear levels) is a result of populations’ close relationships, past demographic phenomena or is characteristic for this species. Detailed analysis of the wsp gene fragment revealed two distinct strains of Wolbachia infecting A. aestivalis. Both of them belong to supergroup A, also found in other arthropods.     

Discordance between morphological and molecular species boundaries among Caribbean species of the reef sponge Callyspongia

Sponges are among the most species‐rich and ecologically important taxa on coral reefs, yet documenting their diversity is difficult due to the simplicity and plasticity of their morphological characters. Genetic attempts to identify species are hampered by the slow rate of mitochondrial sequence evolution characteristic of sponges and some other basal metazoans. Here we determine species boundaries of the Caribbean coral reef sponge genus Callyspongia using a multilocus, model‐based approach. Based on sequence data from one mitochondrial (COI), one ribosomal (28S), and two single‐copy nuclear protein‐coding genes, we found evolutionarily distinct lineages were not concordant with current species designations in Callyspongia. While C. fallax, C. tenerrima, and C. plicifera were reciprocally monophyletic, four taxa with different morphologies (C. armigera, C. longissima, C. eschrichtii, and C. vaginalis) formed a monophyletic group and genetic distances among these taxa overlapped distances within them. A model‐based method of species delimitation supported collapsing these four into a single evolutionary lineage. Variation in spicule size among these four taxa was partitioned geographically, not by current species designations, indicating that in Callyspongia, these key taxonomic characters are poor indicators of genetic differentiation. Taken together, our results suggest a complex relationship between morphology and species boundaries in sponges.     

Molecular phylogeny and diversity of penaeid shrimps (Crustacea: Decapoda) from South-East Asian waters

The decapod family Penaeidae comprises most of the economically important marine shrimp species. Its members are widespread throughout the world, with its highest species diversity centred in the Indo-West Pacific region. Despite this importance, their taxonomy, classification and phylogeny are not yet settled due in part to incongruence among hypotheses proposed from molecular versus morphological studies. In this study, using a thorough taxonomic sampling of especially the South-East Asian species, we aim to (a) utilize a reconstructed phylogeny to test the monophyly of the Penaeidae and its currently recognized genera and (b) explore its species diversity in South-East Asian waters. To infer the phylogeny, a combined gene data set (including 109 ingroup and six outgroup taxa) of mitochondrial genes, COI and 16S rRNA, and two nuclear genes, NaK and PEPCK, was utilized. To explore its diversity, another data set that included 371 COI gene sequences (231 newly generated and 140 retrieved from public sources) was compiled and subsequently analysed with two different tools (ABGD and bPTP) for species delimitation. Other than supporting the non-monophyly of the Penaeidae with the Sicyoniidae nested within the penaeid tribe Trachypenaeini, the genera Penaeus, Mierspenaeopsis and Parapenaeopsis were also revealed to be polyphyletic. Our species delimitation analysis inferred that 94 putative species actually existed within the 71 morphospecies reviewed, indicating an underestimated biodiversity in this family and the potential presence of new species within the following morphospecies: Kishinouyepenaeopsis cornuta, K. incisa, Mierspenaeopsis sculptilis, M hardwicki, Parapenaeopsis coromandelica and Penaeus monodon.     

Molecular Phylogeny of the Chipmunk Genus Tamias Based on the Mitochondrial Cytochrome Oxidase Subunit II Gene

Complete sequences of the mitochondrial cytochrome oxidase subunit II gene were used to construct a phylogeny for 21 of the 25 currently recognized chipmunk species. Phylogenetic analyses indicate that T. striatus (subgenus Tamias, eastern United States) and T. sibiricus (subgenus Eutamias, Asia) are distantly related to the other species (subgenus Neotamias), which constitute a western North American radiation. We discuss and compare our molecular phylogeny to previous taxonomies and present a suggested classification of the species groups for the subgenus Neotamias.     

Molecular phylogeny of the land snail genus Alopia (Gastropoda: Clausiliidae) reveals multiple inversions of chirality

Whereas the vast majority of gastropods possess dextral shell and body organization, members of the Clausiliidae family are almost exclusively sinistral. Within this group a unique feature of the alpine genus Alopia is the comparable representation of sinistral and dextral taxa, and the existence of enantiomorph taxon pairs that appear to differ only in their chirality. We carried out a molecular phylogenetic study, using mitochondrial cytochrome c oxidase subunit I (COI) gene sequences, in order to find out whether chiral inversions are more frequent in this genus than in other genera of land snails. Our results revealed multiple independent inversions in the evolutionary history of Alopia and a close genetic relationship between members of the enantiomorph pairs. The inferred COI phylogeny also provided valuable clues for the taxonomic division and zoogeographical evaluation of Alopia species. The high number of inverse forms indicates unstable fixation of the coiling direction. This deficiency and the availability of enantiomorph pairs may make Alopia species attractive experimental models for genetic studies aimed at elucidating the molecular basis of chiral stability. © 2013 The Linnean Society of London     

Developing a new molecular marker for aphid species identification: Evaluation of eleven candidate genes with species-level sampling

The mitochondrial cytochrome c oxidase subunit I (COI) gene has been utilized as a molecular marker for aphid species identification. However, this gene has sometimes resulted in misidentification because of low interspecific genetic divergences between some species pairs. In this study, to propose new molecular markers for the family Aphididae, we first screened 2289 sequences of 11 genes (COI, COII, CytB, ATP6, lrRNA, srRNA, ITS1, ITS2, EF1a, 18S, and 28S) collected from the GenBank. Among the 11 genes, ATP6 gene revealed the largest genetic divergence among congeneric species with the smallest divergence among conspecific individuals; in contrast, species pairs with low genetic divergences (< 1%) were not observed. Secondly, for statistically testing the usefulness of ATP6 gene in species identification, we analyzed genetic distances between all of the combinations of 32 individuals of 20 species for both COI and ATP6 genes. The ATP6 gene showed lower intraspecific (on average 0.08%) and higher interspecific (on average 8.28%) genetic distances than the COI gene (on average 0.19% and 6.24%, respectively) for the same pairs of individuals. This study corroborates the usefulness of the ATP6 gene as a new molecular marker that could improve the misidentification problems that are inherent with the COI gene.     

DNA sequences of the mitochondrial COI gene have low levels of divergence among deep-sea octocorals (Cnidaria: Anthozoa)

We are analyzing genetic diversity in deep-seamount octocorals with the ultimate goal of studying the effect of retention and dispersal of larvae on genetic population structure. Here we report on the sequence diversity of the mitochondrial cytochrome oxidase I (COI) gene among 11 species. Uncorrected pairwise sequence divergences ranged from 0.4–10.3% for comparisons among species spanning the intrageneric to interordinal levels. Relative to other invertebrates, these divergences are very low, suggesting that COI may not be useful as a genetic marker for studying dispersal among deep-sea octocoral populations. Possible explanations for the reduced rates of divergence observed include a lower rate of evolution for octocoral mitochondrial genomes and the presence of a gene, mtMSH, which may code for a mitochondrial DNA mismatch-repair system. We report the finding of mtMSH in three deep-sea octocorals (Acanthogorgia sp., Corallium ducale, and Paramuricea sp.), which brings the total published observations of this gene to six species, all in the subclass Octocorallia.     

Dissimilarity of Species and Forms of Planktonic Neocalanus Copepods Using Mitochondrial COI, 12S,Nuclear ITS,and 28S Gene Sequences

Background

A total of six Neocalanus species inhabit the oceans of the world. Of these, three species plus form variants (N. cristatus, N. plumchrus, N. flemingeri large form, and N. flemingeri small form), which constitute a monophyletic group among Neocalanus copepods, occur in the Northwestern Pacific off Japan. In the present study, we have tried to discriminate the three species plus form variants of Neocalanus copepods based on sequences of four DNA marker regions.

Methodology/Principal Findings

Discrimination was performed based on the DNA sequence information from four genetic markers, including the mitochondrial COI, 12S, nuclear ITS, and 28S gene regions. Sequence dissimilarity was compared using both distance- and character-based approaches. As a result, all three species were confirmed to be distinct based on the four genetic marker regions. On the contrary, distinction of the form variants was only confirmed based on DNA sequence of the mitochondrial COI gene region.

Conclusions/Significance

Although discrimination was not successful for the form variants based on the mitochondrial 12S, nuclear ITS, and 28S genes, diagnostic nucleotide sequence characters were observed in their mitochondrial COI gene sequences. Therefore, these form variants are considered to be an important unit of evolution below the species level, and constitute a part of the Neocalanus biodiversity.     

Assessing the diversity and distribution of Cephalothrix species (Nemertea: Palaeonemertea) in European waters by comparing different species delimitation methods

Soft‐bodied marine taxa, like ribbon worms (Nemertea), often lack clear diagnostic morphological characters impeding traditional species delimitation. Therefore, recent studies concentrated on molecular genetic methods to solve taxonomic issues. Different delimitation methods were employed to explore species boundaries and the presence of cryptic species. However, the performance of the different delimitation methods needs to be tested. A particularly promising nemertean genus in this regard is the palaeonemertean genus Cephalothrix that is commonly found in European waters. In order to gain information on the number and distribution of European cephalotrichids and to test different tree‐based and non‐tree‐based delimitation methods, we analyzed a dataset comprising the barcoding region of the mitochondrial cytochrome c oxidase subunit I (COI) of 215 European Cephalothrix specimens, of which 78 were collected for this study. Our results show the presence of 12–13 European lineages of which several can be assigned to known European species. Analyzing a second dataset comprising 74 additional sequences from the Pacific and the Atlantic Oceans helped identify some of the unassigned European specimens. One resulting clade seems to represent a non‐native introduced Cephalothrix species, while another has never been recorded from Europe before. In our analysis, especially the tree‐based methods and the phylogenetic analysis proved to be a useful tool when delimiting species. It remains unclear whether the different identified clades result from cryptic speciation or from a high genetic variability of the COI gene.     

Comparison of Chironomus usenicus and Chironomus curabilis with species of the group plumosus (Diptera) inferred from the mitochondrial DNA Gene COI and by the polytene chromosomes banding pattern

The sequence of a 595-bp fragment of the mitochondrial COI gene was determined for the species Chironomus usenicus and Chironomus curabilis of the genus Chironomus. Phylogenetic reconstructions based on the analysis of the COI gene sequence coincide on the whole with cytogenetic data, permitting Ch. usenicus and Ch. curabilis to be regarded as members of the group plumosus. Chironomus usenicus and Ch. plumosus have identical COI gene sequences. Two hypotheses explaining this identity are considered: inheritance of mtDNA from one of the parental species during hybridogenesis and horizontal transfer of mitochondrial genes.     

Exploring phylogenetic informativeness and nuclear copies of mitochondrial DNA (numts) in three commonly used mitochondrial genes: mitochondrial phylogeny of peppermint,cleaner, and semi‐terrestrial shrimps (Caridea: Lysmata,Exhippolysmata, and Merguia)

Of paramount importance to studies that profit from molecular trees is the accuracy and robustness of the reconstructed phylogenies. Causes of systematic error that can mislead phylogenetic methods include nuclear copies of mitochondrial DNA (numts) and low phylogenetic informativeness (PI). Herein, numts and PI were explored in three mitochondrial genes commonly used for phylogenetic reconstruction: 16S, 12S, and cytochrome c oxidase I (COI). Shrimps from the genera Lysmata, Exhippolysmata, and Merguia were used as a model system. The existence of: (1) multiple bands on gels of COI and 12S polymerase chain reaction (PCR) products from various species; (2) double peaks, background noise, and ambiguity in sequence chromatograms of COI and 12S PCR products that produced a single clear band in other species; and (3) indels, stop codons, and considerable composition bias in COI‐like cloned sequences of one problematic species (Lysmata seticaudata), was interpreted as evidence of pervasive non‐functional nuclear copies of mitochondrial DNA (numts) of the targeted COI (and probably 12S) mtDNA fragment. The information content of the three mtDNA markers studied was investigated using PI profiling, spectral analysis, and neighbour‐nets. Marker‐specific PI profiles suggested that the COI marker has the highest information content and greatest power for resolving both shallow and deep nodes in trees depicting the phylogenetic relationship among the species studied. Nonetheless, spectral analysis of splits and neighbour‐nets suggested that the 16S and 12S markers were equally or even more powerful than the COI marker for resolving nodes at all phylogenetic levels. Altogether, these analyses suggest that all three mtDNA markers are equally useful for resolving phylogenetic relationships in the shrimps studied, and that PI profiling is not necessarily useful to estimate overall gene utility. A ‘total‐evidence’ phylogenetic analysis that included 34 species and used a concatenated data set of 1403 characters (from reliable 16S, 12S and COI sequences), demonstrated that the genus Lysmata is paraphyletic, and that the monophyletic clade comprising species of Lysmata and Exhippolysmata can be divided into four well‐supported subclades (Neotropical, Cleaner, Cosmopolitan, and Morphovariable). © 2013 The Linnean Society of London