Species status and phylogeography of two closely related Coptolabrus species (Coleoptera: Carabidae) in South Korea inferred from mitochondrial and nuclear gene sequences

High molecular diversity in Carabus (Hygrocarabus) variolosus and C. nodulosus

The Carabus subgenus Hygrocarabus contains two taxa: C. variolosus and C. nodulosus, the species or subspecies status of which is handled far from uniform in the literature. Both taxa show a similar morphology, the shape of the tip of the aedeagus provides a reliable morphological marker for identification. We analysed two mitochondrial gene parts (COI-5’ and COI-3’) and a nuclear one (ITS2). High diversity was found showing specific geographical patterns. Introgressive hybridisation was detected but interpreted not as an argument for subspecies status because high genetic distances indicated that it must have taken place in former times. In a laboratory hybridisation experiment, the male did not accept the female of the other taxon, supporting the conclusion that these are separate species. A series of refuges was expected for the period of ice ages. Although only the taxon C. variolosus is listed in Annex II and IV of the EU Habitats Directive, C. nodulosus also falls under this listing, as at the time of including the species into the Annexes in 2004, the two taxa were considered subspecies and hence the listing would include both, independent of later taxonomic revisions.

Phylogeographic patterns of Lygus pratensis (Hemiptera: Miridae): Evidence for weak genetic structure and recent expansion in northwest China

Lygus pratensis (L.) is an important cotton pest in China, especially in the northwest region. Nymphs and adults cause serious quality and yield losses. However, the genetic structure and geographic distribution of L. pratensis is not well known. We analyzed genetic diversity, geographical structure, gene flow, and population dynamics of L. pratensis in northwest China using mitochondrial and nuclear sequence datasets to study phylogeographical patterns and demographic history. L. pratensis (n = 286) were collected at sites across an area spanning 2,180,000 km², including the Xinjiang and Gansu-Ningxia regions. Populations in the two regions could be distinguished based on mitochondrial criteria but the overall genetic structure was weak. The nuclear dataset revealed a lack of diagnostic genetic structure across sample areas. Phylogenetic analysis indicated a lack of population level monophyly that may have been caused by incomplete lineage sorting. The Mantel test showed a significant correlation between genetic and geographic distances among the populations based on the mtDNA data. However the nuclear dataset did not show significant correlation. A high level of gene flow among populations was indicated by migration analysis; human activities may have also facilitated insect movement. The availability of irrigation water and ample cotton hosts makes the Xinjiang region well suited for L. pratensis reproduction. Bayesian skyline plot analysis, star-shaped network, and neutrality tests all indicated that L. pratensis has experienced recent population expansion. Climatic changes and extensive areas occupied by host plants have led to population expansion of L. pratensis. In conclusion, the present distribution and phylogeographic pattern of L. pratensis was influenced by climate, human activities, and availability of plant hosts.

A DNA barcode library for ground beetles (Insecta, Coleoptera, Carabidae) of Germany: The genus Bembidion Latreille, 1802 and allied taxa

As molecular identification method, DNA barcoding based on partial cytochrome c oxidase subunit 1 (COI) sequences has been proven to be a useful tool for species determination in many insect taxa including ground beetles. In this study we tested the effectiveness of DNA barcodes to discriminate species of the ground beetle genus Bembidion and some closely related taxa of Germany. DNA barcodes were obtained from 819 individuals and 78 species, including sequences from previous studies as well as more than 300 new generated DNA barcodes. We found a 1:1 correspondence between BIN and traditionally recognized species for 69 species (89%). Low interspecific distances with maximum pairwise K2P values below 2.2% were found for three species pairs, including two species pairs with haplotype sharing (Bembidion atrocaeruleum/Bembidion varicolor and Bembidion guttula/Bembidion mannerheimii). In contrast to this, deep intraspecific sequence divergences with distinct lineages were revealed for two species (Bembidion geniculatum/Ocys harpaloides). Our study emphasizes the use of DNA barcodes for the identification of the analyzed ground beetles species and represents an important step in building-up a comprehensive barcode library for the Carabidae in Germany and Central Europe as well.

The mitochondrial genome of the plant bug Apolygus lucorum (Hemiptera: Miridae): Presently known as the smallest in Heteroptera

The complete mitochondrial (mt) genome of the plant bug, Apolygus lucorum, an important cotton pest, has been sequenced and annotated in this study. The entire circular genome is 14 768 bp in size and represents the smallest in presently known heteropteran mt genomes. The mt genome is encoding for two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, 13 protein coding genes and a control region, and the order, content, codon usage and base organization show similarity to a great extent to the hypothetical ancestral model. All protein coding genes use standard initiation codons ATN. Conventional stop codons TAA and TAG have been assigned to the most protein coding genes; however, COIII, ND4 and ND5 genes show incomplete terminator signal (T). All tRNA genes possess the typical clover leaf structure, but the dihydrouridine arm of tRNA(Ser(AGN)) only forms a simple loop. Secondary structure models of rRNA genes are generally in accordance with the former models, although some differences exist in certain parts. Three intergenic spacers have never been found in sequenced mt genomes of Heteroptera. The phylogenetic study based on protein coding genes is largely congruent with previous phylogenetic work. Both Bayesian inference and maximum likelihood analyses highly support the sister relationship of A. lucorum and Lygus lineolaris, and Miridae presents a sister position to Anthocoridae.

A fuzzy-set-theory-based approach to analyse species membership in DNA barcoding

Reliable assignment of an unknown query sequence to its correct species remains a methodological problem for the growing field of DNA barcoding. While great advances have been achieved recently, species identification from barcodes can still be unreliable if the relevant biodiversity has been insufficiently sampled. We here propose a new notion of species membership for DNA barcoding-fuzzy membership, based on fuzzy set theory-and illustrate its successful application to four real data sets (bats, fishes, butterflies and flies) with more than 5000 random simulations. Two of the data sets comprise especially dense species/population-level samples. In comparison with current DNA barcoding methods, the newly proposed minimum distance (MD) plus fuzzy set approach, and another computationally simple method, 'best close match', outperform two computationally sophisticated Bayesian and BootstrapNJ methods. The new method proposed here has great power in reducing false-positive species identification compared with other methods when conspecifics of the query are absent from the reference database.

Morphometric and molecular analyses of the sand fly species Lutzomyia shannoni (Diptera: Psychodidae: Phlebotominae) collected from seven different geographical areas in the southeastern United States

A morphometric and molecular study of adult male and female Lutzomyia shannoni (Dyar 1929) collected at seven different locations within the southeastern United States was conducted to assess the degree of divergence between the grouped specimens from each location. The collection locations were as follows: Fort Bragg, NC; Fort Campbell, KY; Fort Rucker, AL; Ossabaw Island, GA; Patuxent National Wildlife Research Refuge, MD; Suwannee National Wildlife Refuge, FL; and Baton Rouge, LA. Forty males and forty females from each location were analyzed morphometrically from 54 and 49 character measurements, respectively. In addition, the molecular markers consisting of the partial cytochrome c oxidase subunit I (from 105 sand flies: 15 specimens/collection site) and the partial internal transcribed spacer 2 (from 42 sand flies: six specimens/collection site) were compared. Multivariate analyses indicate that the low degree of variation between the grouped specimens from each collection site prevents the separation of any collection site into an entity that could be interpreted as a distinct population. The molecular analyses were in concordance with the morphometric study as no collection location grouped into a separate population based on the two partial markers. The grouped specimens from each collection site appear to be within the normal variance of the species, indicating a single population in the southeast United States. It is recommended that additional character analyses of L. shannoni based on more molecular markers, behavioral, ecological, and physiological characteristics, be conducted before ruling out the possibility of populations or a cryptic species complex within the southeastern United States.

Revisiting the insect mitochondrial molecular clock: the mid-Aegean trench calibration

Phylogenetic trees in insects are frequently dated by applying a "standard" mitochondrial DNA (mtDNA) clock estimated at 2.3% My(-1), but despite its wide use reliable calibration points have been lacking. Here, we used a well-established biogeographic barrier, the mid-Aegean trench separating the western and eastern Aegean archipelago, to estimate substitution rates in tenebrionid beetles. Cytochrome oxidase I (cox1) for six codistributed genera across 28 islands (444 individuals) on both sides of the mid-Aegean trench revealed 60 independently coalescing entities delimited with a mixed Yule-coalescent model. One representative per entity was used for phylogenetic analysis of mitochondrial (cox1, 16S rRNA) and nuclear (Mp20, 28S rRNA) genes. Six nodes marked geographically congruent east-west splits whose separation was largely contemporaneous and likely to reflect the formation of the mid-Aegean trench at 9-12 Mya. Based on these "known" dates, a divergence rate of 3.54% My(-1) for the cox1 gene (2.69% when combined with the 16S rRNA gene) was obtained under the preferred partitioning scheme and substitution model selected using Bayes factors. An extensive survey suggests that discrepancies in mtDNA substitution rates in the entomological literature can be attributed to the use of different substitution models, the use of different mitochondrial gene regions, mixing of intraspecific with interspecific data, and not accounting for variance in coalescent times or postseparation gene flow. Different treatments of these factors in the literature confound estimates of mtDNA substitution rates in opposing directions and obscure lineage-specific differences in rates when comparing data from various sources.

Nuclear and mitochondrial sequences confirm complex colonization patterns and clear species boundaries for flightless weevils in the Galápagos archipelago

Nuclear sequence data were collected from endemic Galápagos species and an introduced close relative, and contrasted with mitochondrial DNA sequences, continuing investigation into the colonization history and modes of diversification in the weevil genus Galapaganus. The current combined phylogeny together with previously published penalized likelihood age estimates builds a complex picture of the archipelago's colonization history. The present reconstruction relies on submerged platforms to explain the early divergence of the young southern Isabela endemics or the Española or San Cristobal populations. Diversity is later built through inter-island divergence starting on older islands and continuing on two simultaneous tracks towards younger islands. The amount of diversity generated through intra-island processes is skewed towards older islands, suggesting that island age significantly influences diversity. Phylogenetic concordance between nuclear and mitochondrial datasets and well-supported monophyletic species in mitochondrial derived topologies appear to reject the possibility of inter-species hybridization. These clear species boundaries might be related to the tight host associations of adult weevils in discrete ecological zones. If shared hosts facilitate hybridization, then host- or habitat-promoted divergences could prevent it, even in the case of species that share islands, since the altitudinal partitioning of habitats minimizes range overlap.

Mechanical barriers to introgressive hybridization revealed by mitochondrial introgression patterns in Ohomopterus ground beetle assemblages

To reveal the role of diverged body size and genital morphology in reproductive isolation among closely related species, we examined patterns of, and factors limiting, introgressive hybridization between sympatric Ohomopterus ground beetles in central Japan using mitochondrial NADH dehydrogenase subunit 5 (ND5) gene sequences. We sampled 17 local assemblages that consisted of two to five species and estimated levels of interspecific gene flow using the genetic distance, D(A), and maximum-likelihood estimates of gene flow. Sharing of haplotypes or haplotype lineages was detected between six of seven species that occurred in the study areas, indicating mitochondrial introgression. The intensity and direction of mitochondrial gene flow were variable among species pairs. To determine the factors affecting introgression patterns, we tested the relationships between interspecific D(A) and five independent variables: difference in body size, difference in genital size, phylogenetic relatedness (nuclear gene sequence divergence), habitat difference, and species richness of the assemblage. Body and genital size differences contributed significantly to preventing gene flow. Thus, mechanical isolation mechanisms reduce the chance of introgressive hybridization between closely related species. Our results highlight the role of morphological divergence in speciation and assemblage formation processes through mechanical isolation.

Phylogeny of pteromalid parasitic wasps (Hymenoptera: Pteromalidae): Initial evidence from four protein-coding nuclear genes

Chalcidoidea (approximately 22,000 described species) is the most ecologically diverse superfamily of parasitic Hymenoptera and plays a major role in the biological control of insect pests. However, phylogenetic relationships both within and between chalcidoid families have been poorly understood, particularly for the large family Pteromalidae and relatives. Forty-two taxa, broadly representing Chalcidoidea but concentrated in the 'pteromalid lineage,' were sequenced for 4620 bp of protein-coding sequence from four nuclear genes for which we present new primers. These are: CAD (1719 bp) DDC (708 bp), enolase (1149 bp), and PEPCK (1044 bp). The combined data set was analyzed using parsimony, maximum likelihood, and Bayesian methods. Statistical significance of the apparent non-monophyly of some taxonomic groups on our trees was evaluated using the approximately unbiased test of Shimodaira [Shimodaira, H. 2002. An approximately unbiased test of phylogenetic tree selection. Syst. Biol. 51(3), 492-508]. In accord with previous studies, we find moderate to strong support for monophyly of Chalcidoidea, a sister-group relationship of Mymaridae to the remainder of Chalcidoidea, and a relatively basal placement of Encarsia (Aphelinidae) within the latter. The 'pteromalid lineage' of families is generally recovered as monophyletic, but the hypothesis of monophyly for Pteromalidae, which appear paraphyletic with respect to all other families sampled in that lineage, is decisively rejected (P < 10(-14)). Within Pteromalidae, monophyly was strongly supported for nearly all tribes represented by multiple exemplars, and for two subfamilies. All other multiply-represented subfamilies appeared para- or polyphyletic in our trees, although monophyly was significantly rejected only for Miscogasterinae, Ormocerinae, and Colotrechninae. The limited resolution obtained in the analyses presented here reinforces the idea that reconstruction of pteromalid phylogeny is a difficult problem, possibly due to rapid radiation of many chalcidoid taxa. Initial phylogenetic comparisons of life history traits suggest that the ancestral chalcidoid was small-bodied and parasitized insect eggs.

Phylogeography and introgressive hybridization of the ground beetle Carabus yamato in Japan based on mitochondrial gene sequences

To study the phylogeography of the ground beetle Carabus yamato in Japan, we compared 1,020-bp sequences of the mitochondrial NADH dehydrogenase subunit 5 (ND5) gene from 373 specimens from 37 localities with those of three parapatric species (C. albrechti, C. kimurai, and C. japonicus) that might share mitochondrial lineages with C. yamato through introgressive hybridization. We found 81 haplotypes from C. yamato. Of these, 17 haplotypes were considered to be of an introgressed lineage from C. albrechti, based on the phylogeny and geographic distribution. In addition, one haplotype of C. kimurai was likely an introgressant from C. yamato. Putative introgression events among the four species were restricted to these two directional cases. We analyzed the phylogeography of C. yamato using nested clade phylogeographical analysis and population genetic parameters. The mitochondrial lineages of C. yamato were estimated to have diverged no more than approximately 1.12 million years ago, implying that the estimated historical events occurred after the Early Pleistocene. Carabus yamato was inferred to have experienced a contraction of its distribution range, followed by recent range expansion. Populations in the western and eastern regions, segregated by Ise Bay and the Nobi and Okazaki Plains, diverged in the mitochondrial clades. The northern and most western populations possessed one clade only (except an introgressed lineage), whereas eastern and some southwestern populations possessed several diverged clades, which were considered to be ancestral; these populations may have been associated with refugia during glacial periods.

Nuclear gene sequences resolve species phylogeny and mitochondrial introgression in Leptocarabus beetles showing trans-species polymorphisms

We studied the phylogenetic relationships among Japanese Leptocarabus ground beetles, which show extensive trans-species polymorphisms in mitochondrial gene genealogies. Simultaneous analysis of combined nuclear data with partial sequences from the long-wavelength rhodopsin, wingless, phosphoenolpyruvate carboxykinase, and 28S rRNA genes resolved the relationships among the five species, although separate analyses of these genes provided topologies with low resolution. For both the nuclear gene tree resulting from the combined data from four genes and a mitochondrial cytochrome oxidase subunit I (COI) gene tree, we applied a Bayesian divergence time estimation using a common calibration method to identify mitochondrial introgression events that occurred after speciation. Three mitochondrial lineages shared by two or three species were likely subject to introgression due to interspecific hybridization because the coalescent times for these lineages were much shorter than the corresponding speciation times estimated from nuclear gene sequences. We demonstrated that when species phylogeny is fully resolved with nuclear gene sequence data, comparative analysis of nuclear and mitochondrial gene trees can be used to infer introgressive hybridization events that might cause trans-species polymorphisms in mitochondrial gene trees.

THE AUTOMATION AND EVALUATION OF NESTED CLADE PHYLOGEOGRAPHIC ANALYSIS

Nested clade phylogeographic analysis (NCPA) is a popular method for reconstructing the demographic history of spatially distributed populations from genetic data. Although some parts of the analysis are automated, there is no unique and widely followed algorithm for doing this in its entirety, beginning with the data, and ending with the inferences drawn from the data. This article describes a method that automates NCPA, thereby providing a framework for replicating analyses in an objective way. To do so, a number of decisions need to be made so that the automated implementation is representative of previous analyses. We review how the NCPA procedure has evolved since its inception and conclude that there is scope for some variability in the manual application of NCPA. We apply the automated software to three published datasets previously analyzed manually and replicate many details of the manual analyses, suggesting that the current algorithm is representative of how a typical user will perform NCPA. We simulate a large number of replicate datasets for geographically distributed, but entirely random-mating, populations. These are then analyzed using the automated NCPA algorithm. Results indicate that NCPA tends to give a high frequency of false positives. In our simulations we observe that 14% of the clades give a conclusive inference that a demographic event has occurred, and that 75% of the datasets have at least one clade that gives such an inference. This is mainly due to the generation of multiple statistics per clade, of which only one is required to be significant to apply the inference key. We survey the inferences that have been made in recent publications and show that the most commonly inferred processes (restricted gene flow with isolation by distance and contiguous range expansion) are those that are commonly inferred in our simulations. However, published datasets typically yield a richer set of inferences with NCPA than obtained in our random-mating simulations, and further testing of NCPA with models of structured populations is necessary to examine its accuracy.

Evolution and diversification of the forest and hypogean ground-beetle genus Trechus in the Canary Islands

The beetle genus Trechus (Carabidae) is represented in the Macaronesian Islands by 43 endemic species. The Canary Islands have 16 endemic species, with two adapted to hypogean life. Phylogenetic relationships among 177 individuals of 38 Canarian, Madeiran, Azorean and continental Trechus species were examined using mitochondrial DNA and nuclear internal transcribed spacer 2 (ITS2) sequence data. Results show two main lineages in the Canaries: one comprising two sister groups with species from the laurel forest of La Gomera and Tenerife, and the other containing the single species from Gran Canaria and a species complex in the four western islands including two troglobites. Calibrations were applied to a linearized tree using a relaxed molecular clock method to estimate the major evolutionary divergence times of the Canarian Trechus species. Although the species assemblage in this archipelago is relatively ancient (7-8 million years), much of the species diversity is recent. Transition to the hypogean environment is more consistent with the "adaptive shift" rather than with the competing "climatic relict" hypothesis.

Comparative phylogeography of three Leptocarabus ground beetle species in South Korea, based on the mitochondrial COI and nuclear 28S rRNA genes

We analyzed the intraspecific gene genealogies of three Leptocarabus ground beetle species (L. seishinensis, L. semiopacus, L. koreanus) in South Korea using sequence data from the mitochondrial cytochrome oxidase subunit I (COI) and nuclear 28S rRNA (28S) genes, and compared phylogeographical patterns among the species. The COI data detected significant genetic differentiation among local populations of all three species, whereas the 28S data showed genetic differentiation only for L. seishinensis. The clearest differentiation of L. seishinensis among local populations was between the northern and southern regions in the COI clades, whereas the 28S clade, which likely indicates relatively ancient events, revealed a range expansion across the northern and southern regions. Leptocarabus semiopacus had the most shallow differentiation of the COI haplotypes, and some clades occurred across the northern and southern regions. In L. koreanus, four diverged COI clades occurred in different regions, with partial overlaps. We discuss the difference in phylogeographical patterns among these Leptocarabus species, as well as between these and other groups of carabid beetles in South Korea.

Multilocus genealogies reveal multiple cryptic species and biogeographical complexity in the California turret spider Antrodiaetus riversi (Mygalomorphae, Antrodiaetidae)

Antrodiaetus riversi (Araneae, Antrodiaetidae) is a dispersal-limited, habitat specialized mygalomorph spider species endemic to mesic woodlands of northern and central California. This species occupies a disjunct distribution, with populations in the Sierra Nevada and Coast Ranges, separated by the inhospitable Central Valley. Previous studies of morphological and allozyme variation have suggested that these populations may constitute cryptic species. We investigated the phylogeography of A. riversi using both nuclear and mitochondrial DNA sequences, collected for a comprehensive population sample. These data reveal the presence of at least five species in the A. riversi complex - these species are deeply diverged, and genealogically exclusive in both nuclear and mitochondrial genomes. Each of these species is characterized by extreme population subdivision and deep phylogeographical structuring, consistent with minimal gene flow across the dissected Californian landscape. Three species are restricted to the Coast Ranges, one to high altitudes of the central Sierran Nevada, and one species is found in both ranges. These species have allopatric distributions, although species parapatry is hypothesized to occur in several areas. Species diversification appears to have pulsed in the Late Miocene/Early Pliocene, a timing consistent with biogeographical reconstructions for many Californian taxa, and a time of turbulent geological activity in the region.

Ribosomal RNA as molecular barcodes: a simple correlation analysis without sequence alignment

MOTIVATION

We explored the feasibility of using unaligned rRNA gene sequences as DNA barcodes, based on correlation analysis of composition vectors (CVs) derived from nucleotide strings. We tested this method with seven rRNA (including 12, 16, 18, 26 and 28S) datasets from a wide variety of organisms (from archaea to tetrapods) at taxonomic levels ranging from class to species.

RESULT

Our results indicate that grouping of taxa based on CV analysis is always in good agreement with the phylogenetic trees generated by traditional approaches, although in some cases the relationships among the higher systemic groups may differ. The effectiveness of our analysis might be related to the length and divergence among sequences in a dataset. Nevertheless, the correct grouping of sequences and accurate assignment of unknown taxa make our analysis a reliable and convenient approach in analyzing unaligned sequence datasets of various rRNAs for barcoding purposes.

AVAILABILITY

The newly designed software (CVTree 1.0) is publicly available at the Composition Vector Tree (CVTree) web server http://cvtree.cbi.pku.edu.cn.

Utility of Nuclear Allele Networks for the Analysis of Closely Related Species in the Genus Carabus, Subgenus Ohomopterus

Nuclear DNA sequence data for diploid organisms are potentially a rich source of phylogenetic information for disentangling the evolutionary relationships of closely related organisms, but present special phylogenetic problems owing to difficulties arising from heterozygosity and recombination. We analyzed allelic relationships for two nuclear gene regions (phosphoenolpyruvate carboxykinase and elongation factor-1a), along with a mitochondrial gene region (NADH dehydrogenase subunit 5), for an assemblage of closely related species of carabid beetles (Carabus subgenus Ohomopterus). We used a network approach to examine whether the nuclear gene sequences provide substantial phylogenetic information on species relationships and evolutionary history. The mitochondrial gene genealogy strongly contradicted the morphological species boundary as a result of introgression of heterospecific mitochondria. Two nuclear gene regions showed high allelic diversity within species, and this diversity was partially attributable to recombination between various alleles and high variability in the intron region. Shared nuclear alleles among species were rare and were considered to represent shared ancestral polymorphism. Despite the presence of recombination, nuclear allelic networks recovered species monophyly more often and presented genetic differentiation patterns (low to high) among species more clearly. Overall, nuclear gene networks provide clear evidence for separate biological species and information on the phylogenetic relationships among closely related carabid beetles.