Bombus terrestris in a mass‐flowering pollinator‐dependent crop: A mutualistic relationship?

Bumblebees (Bombus spp.) rely on an abundant and diverse selection of floral resources to meet their nutritional requirements. In farmed landscapes, mass‐flowering crops can provide an important forage resource for bumblebees, with increased visitation from bumblebees into mass‐flowering crops having an additional benefit to growers who require pollination services. This study explores the mutualistic relationship between Bombus terrestris L. (buff‐tailed bumblebee), a common species in European farmland, and the mass‐flowering crop courgette (Cucurbita pepo L.) to see how effective B. terrestris is at pollinating courgette and in return how courgette may affect B. terrestris colony dynamics. By combining empirical data on nectar and pollen availability with model simulations using the novel bumblebee model Bumble‐BEEHAVE, we were able to quantify and simulate for the first time, the importance of courgette as a mass‐flowering forage resource for bumblebees. Courgette provides vast quantities of nectar to ensure a high visitation rate, which combined with abundant pollen grains, enables B. terrestris to have a high pollination potential. While B. terrestris showed a strong fidelity to courgette flowers for nectar, courgette pollen was not found in any pollen loads from returning foragers. Nonetheless, model simulations showed that early season courgette (nectar) increased the number of hibernating queens, colonies, and adult workers in the modeled landscapes. Synthesis and applications. Courgette has the potential to improve bumblebee population dynamics; however, the lack of evidence of the bees collecting courgette pollen in this study suggests that bees can only benefit from this transient nectar source if alternative floral resources, particularly pollen, are also available to fulfill bees’ nutritional requirements in space and time. Therefore, providing additional forage resources could simultaneously improve pollination services and bumblebee populations.     

Targeted agri‐environment schemes significantly improve the population size of common farmland bumblebee species

Changes in agricultural practice across Europe and North America have been associated with range contractions and local extinction of bumblebees (Bombus spp.). A number of agri‐environment schemes have been implemented to halt and reverse these declines, predominantly revolving around the provision of additional forage plants. Although it has been demonstrated that these schemes can attract substantial numbers of foraging bumblebees, it remains unclear to what extent they actually increase bumblebee populations. We used standardized transect walks and molecular techniques to compare the size of bumblebee populations between Higher Level Stewardship (HLS) farms implementing pollinator‐friendly schemes and Entry Level Stewardship (ELS) control farms. Bumblebee abundance on the transect walks was significantly higher on HLS farms than ELS farms. Molecular analysis suggested maximum foraging ranges of 566 m for Bombus hortorum, 714 m for B. lapidarius, 363 m for B. pascuorum and 799 m for B. terrestris. Substantial differences in maximum foraging range were found within bumblebee species between farm types. Accounting for foraging range differences, B. hortorum (47 vs 13 nests/km²) and B. lapidarius (45 vs 22 nests/km²) were found to nest at significantly greater densities on HLS farms than ELS farms. There were no significant differences between farm type for B. terrestris (88 vs 38 nests/km²) and B. pascuorum (32 vs 39 nests/km²). Across all bumblebee species, HLS management had a significantly positive effect on bumblebee nest density. These results show that targeted agri‐environment schemes that increase the availability of suitable forage can significantly increase the size of wild bumblebee populations.     

Molecular detection of Cotesia vestalis (Hymenoptera: Braconidae) in the beet webworm Loxostege sticticalis L. (Lepidoptera: Crambidae)

Genomic DNA samples from larvae of the beet webworm Loxostege sticticalis collected in the south‐western Russia were used to amplify mitochondrial cytochrome oxidase unit I (COI) gene. In a small proportion of samples, the sequenced product showed considerable heterogeneity due to admixture of a minor sequence. A preliminary BLAST analysis of a 100‐bp‐long fragment of this minor sequence showed its maximal similarity to the COI gene region of Cotesia, a genus of braconid larval endoparasitoids of Lepidoptera. An additional primer was designed to specifically amplify ca 300 bp of the COI gene region from Braconidae. As many as seven of 25 samples were positive by PCR. Sequencing of the amplified products in all these samples showed nucleotide sequence identity to the COI region of Cotesia vestalis (Cotesia plutellae) and the presence of two molecular haplotypes among individual parasitoid samples.     

DNA barcoding commercially important fish species of Turkey

DNA barcoding was used in the identification of 89 commercially important freshwater and marine fish species found in Turkish ichthyofauna. A total of 1765 DNA barcodes using a 654‐bp‐long fragment of the mitochondrial cytochrome c oxidase subunit I gene were generated for 89 commercially important freshwater and marine fish species found in Turkish ichthyofauna. These species belong to 70 genera, 40 families and 19 orders from class Actinopterygii, and all were associated with a distinct DNA barcode. Nine and 12 of the COI barcode clusters represent the first species records submitted to the BOLD and GenBank databases, respectively. All COI barcodes (except sequences of first species records) were matched with reference sequences of expected species, according to morphological identification. Average nucleotide frequencies of the data set were calculated as T = 29.7%, C = 28.2%, A = 23.6% and G = 18.6%. Average pairwise genetic distance among individuals were estimated as 0.32%, 9.62%, 17,90% and 22.40% for conspecific, congeneric, confamilial and within order, respectively. Kimura 2‐parameter genetic distance values were found to increase with taxonomic level. For most of the species analysed in our data set, there is a barcoding gap, and an overlap in the barcoding gap exists for only two genera. Neighbour‐joining trees were drawn based on DNA barcodes and all the specimens clustered in agreement with their taxonomic classification at species level. Results of this study supported DNA barcoding as an efficient molecular tool for a better monitoring, conservation and management of fisheries.     

Identification of Belgian mosquito species (Diptera: Culicidae) by DNA barcoding

Since its introduction in 2003, DNA barcoding has proven to be a promising method for the identification of many taxa, including mosquitoes (Diptera: Culicidae). Many mosquito species are potential vectors of pathogens, and correct identification in all life stages is essential for effective mosquito monitoring and control. To use DNA barcoding for species identification, a reliable and comprehensive reference database of verified DNA sequences is required. Hence, DNA sequence diversity of mosquitoes in Belgium was assessed using a 658 bp fragment of the mitochondrial cytochrome oxidase I (COI) gene, and a reference data set was established. Most species appeared as well‐supported clusters. Intraspecific Kimura 2‐parameter (K2P) distances averaged 0.7%, and the maximum observed K2P distance was 6.2% for Aedes koreicus. A small overlap between intra‐ and interspecific K2P distances for congeneric sequences was observed. Overall, the identification success using best match and the best close match criteria were high, that is above 98%. No clear genetic division was found between the closely related species Aedes annulipes and Aedes cantans, which can be confused using morphological identification only. The members of the Anopheles maculipennis complex, that is Anopheles maculipennis s.s. and An. messeae, were weakly supported as monophyletic taxa. This study showed that DNA barcoding offers a reliable framework for mosquito species identification in Belgium except for some closely related species.     

Microsympatry in cryptic lowland salamanders (Caudata: Plethodontidae: Bolitoglossa subgenus Nanotriton) from north‐western Honduras: implications for taxonomy and regional biogeography

The taxonomy of the diminutive bolitoglossine salamanders of the subgenus Nanotriton has previously been a source of confusion among specialists. At various times, at least three different species of Nanotriton have been reported to occur in Honduras: Bolitoglossa nympha, B. occidentalis and B. rufescens. Phylogenetic analysis of 16S and cyt b from samples from three localities in Honduras (departments of Yoro, Copán and Cortés) confirms that most populations are assignable to B. nympha. Nine samples collected from a single locality on the northern slope of the Sierra de Omoa, between 120 and 190 m in elevation, were found to represent two distinct taxa: B. nympha and B. rufescens. Field examination of these nine specimens found them to be morphologically indistinguishable, including a lack of maxillary teeth in all specimens. As such, in situ identification of the two species in north‐western Honduras is problematic given the occurrence of the two taxa in microsympatry in at least one locality. The discovery of yet another divergent lineage of B. rufescens highlights the need for a taxonomic reassessment within this species complex.     

Species‐Specific PCR‐Based Assay for Identification and Detection of Phomopsis (Diaporthe) azadirachtae Causing Die‐Back Disease in Azadirachta indica

Die‐back disease caused by Phomopsis (Diaporthe) azadirachtae is the devastating disease of Azadirachta indica. Accurate identification of P. azadirachtae is always problematic due to morphological plasticity and delayed appearance of conidia. A species‐specific PCR‐based assay was developed for rapid and reliable identification of P. azadirachtae by designing a species‐specific primer‐targeting ITS region of P. azadirachtae isolates. The assay was validated with DNA isolated from different Phomopsis species and other fungal isolates. The PCR assay amplified 313‐bp product from all the isolates of P. azadirachtae and not from any other Phomopsis species or any genera indicating its specificity. The assay successfully detected the pathogen DNA in naturally and artificially infected neem seeds and twigs indicating its applicability in seed quarantine and seed health testing. The sensitivity of the assay was 100 fg when genomic DNA of all isolates was analysed. The PCR‐based assay was 92% effective in comparison with seed plating technique in detecting the pathogen. This is the first report on the development of species‐specific PCR assay for identification and detection of P. azadirachtae. Thus, PCR‐based assay developed is very specific, rapid, confirmatory and sensitive tool for detection of pathogen P. azadirachtae at early stages.     

DNA barcoding of Bradysia (Diptera: Sciaridae) for detection of the immature stages on agricultural crops

We investigated the usefulness of mitochondrial cytochrome c oxidase (COI) DNA barcoding of the genus Bradysia for the detection of immature stages and cryptic species complex. Although the larvae of some species in this genus are agricultural pests, immature stages are rarely identified due to the lack of key morphological characteristics. We constructed partial sequences of the COI gene for 25 species of Bradysia as a first step towards a DNA barcode. Using these data, Bradysia impatiens, B. procera and B. peraffinis were identified from larval specimens collected, respectively, from paprika, ginseng and oak sawdust beds used for cultivating shiitake. Our findings reveal a complex of three species within the B. tilicola group. These species were all identified as important pest B. ocellaris based on the morphology of male genital structures; however, the interspecific genetic divergence of the COI region was significantly greater (16.1–19.4%) than the intraspecific variation in each species. Therefore, B. ocellaris may consist of at least three species. The results demonstrate that COI DNA barcodes are useful for Bradysia species identification.     

Daphnia diversity on the Tibetan Plateau measured by DNA taxonomy

Daphnia on the Tibetan Plateau has been little studied, and information on species diversity and biogeography is lacking. Here, we conducted a 4‐year survey using the barcoding fragment of the mitochondrial COI gene to determine the distribution and diversity of Daphnia species found across the Plateau. Our results show that species richness is higher than previously thought, with total described and provisional species number doubling from 5 to 10. Six of the taxonomic units recovered by DNA taxonomy agreed well with morphology, but DNA barcoding distinguished three clades each for the D. longispina (D. galeata, D. dentifera, and D. longispina) and D. pulex (D. pulex, D. cf. tenebrosa, and D. pulicaria) complexes. The sequence divergence between congeneric species varied within a large range, from 9.25% to 30.71%. The endemic D. tibetana was the most common and widespread species, occurring in 12 hyposaline to mesosaline lakes. The lineage of D. longispina is the first confirmed occurrence in west Tibet.     

First Chemical Analysis and Characterization of the Male Species‐Specific Cephalic Labial‐Gland Secretions of South American Bumblebees

The evolution of signals and reproductive traits involved in the pre‐mating recognition has been in focus of abundant research in several model species, such as bumblebees (genus Bombus). However, the most‐studied bumblebee reproductive trait, the male cephalic labial gland secretions (CLGS), remains unknown among bumblebee species from South America. In this study, the CLGS of five South American bumblebees of the subgenera Thoracobombus (Bombus excellens and B. atratus) and Cullumanobombus (B. rubicundus, B. hortulanus, and B. melaleucus) were investigated, by comparing the chemical compositions of their secretions to those of closely related European species. The results showed an obvious interspecific differentiation in both subgenera. The interspecific differentiation among the species of the Thoracobombus subgenus involved different compounds present at high contents (main compounds), while those of the Cullumanobombus subgenus shared the same main components. This suggests that among the species of the Cullumanobombus subgenus, the differentiation in minor components could lead to species discrimination.     

Deep sympatric mtDNA divergence in the autumnal moth (Epirrita autumnata)

Deep sympatric intraspecific divergence in mtDNA may reflect cryptic species or formerly distinct lineages in the process of remerging. Preliminary results from DNA barcoding of Scandinavian butterflies and moths showed high intraspecific sequence variation in the autumnal moth, Epirrita autumnata. In this study, specimens from different localities in Norway and some samples from Finland and Scotland, with two congeneric species as outgroups, were sequenced with mitochondrial and nuclear markers to resolve the discrepancy found between mtDNA divergence and present species‐level taxonomy. We found five COI sub‐clades within the E. autumnata complex, most of which were sympatric and with little geographic structure. Nuclear markers (ITS2 and Wingless) showed little variation and gave no indications that E. autumnata comprises more than one species. The samples were screened with primers for Wolbachia outer surface gene (wsp) and 12% of the samples tested positive. Two Wolbachia strains were associated with different mtDNA sub‐clades within E. autumnata, which may indicate indirect selection/selective sweeps on haplotypes. Our results demonstrate that deep mtDNA divergences are not synonymous with cryptic speciation and this has important implications for the use of mtDNA in species delimitation, like in DNA barcoding.     

Species‐delimitation and phylogenetic analyses of some cosmopolitan species of Hypnea (Rhodophyta) reveal synonyms and misapplied names to H. cervicornis,including a new species from Brazil

Hypnea has an intricate nomenclatural history due to a wide pantropical distribution and considerable morphological variation. Recent molecular studies have provided further clarification on the systematics of the genus; however, species of uncertain affinities remain due to flawed taxonomic identification. Detailed analyses coupled with literature review indicated a strong relationship among H. aspera, H. cervicornis, H. flexicaulis, and H. tenuis, suggesting a need for further taxonomic studies. Here, we analyzed sequences from two molecular markers (COI‐5P and rbcL) and performed several DNA‐based delimitation methods (mBGD, ABGD, SPN, PTP and GMYC). These molecular approaches were contrasted with morphological and phylogenetic evidence from type specimens and/or topotype collections of related species under a conservative approach. Our results demonstrate that H. aspera and H. flexicaulis represent heterotypic synonyms of H. cervicornis and indicate the existence of a misidentified Hypnea species, widely distributed on the Brazilian coast, described here as a new species: H. brasiliensis. Finally, inconsistencies observed among our results based on six different species delimitation methods evidence the need for adequate sampling and marker choice for different methods.     

A molecular tool to identify Anastatus parasitoids of the brown marmorated stink bug

Globally, Anastatus species (Hymenoptera: Eupelmidae) are associated with the invasive agricultural pest Halyomorpha halys (Stål) (Hemiptera: Pentatomidae). In Europe, the polyphagous Anastatus bifasciatus (Geoffroy) is the most prevalent native egg parasitoid on H. halys eggs and is currently being tested as a candidate for augmentative biological control. Anastatus bifasciatus frequently displays behavior without oviposition, and induces additional host mortality through oviposition damage and host feeding that is not measured with offspring emergence. This exacerbates accurate assessment of parasitism and host impact, which is crucial for efficacy evaluation as well as for pre‐ and post‐release risk assessment. To address this, a general Anastatus primer set amplifying a 318‐bp fragment within the barcoding region of the cytochrome oxidase I (COI) gene was developed. When challenged with DNA of three Anastatus species —A. bifasciatus, Anastatus japonicus Ashmead, and Anastatus sp.—, five scelionid parasitoid species that might be encountered in the same host environments and 11 pentatomid host species, only Anastatus DNA was successfully amplified. When applied to eggs of the target host, H. halys, and an exemplary non‐target host, Dendrolimus pini L. (Lepidoptera: Lasiocampidae), subjected to host feeding, no Anastatus amplicons were produced. Eggs of the two host species containing A. bifasciatus parasitoid stages, from 1‐h‐old eggs to pupae, and emerged eggs yielded Anastatus fragments. Confirmation of parasitoid presence with dissections and subsequent PCRs with the developed primer pair resulted in 95% success for 1‐h‐old parasitoid eggs. For both host species, field‐exposed sentinel emerged eggs stored dry for 6 months, 100% of the specimens produced Anastatus amplicons. This DNA‐based screening method can be used in combination with conventional methods to better interpret host‐parasitoid and parasitoid‐parasitoid interactions. It will help address ecological questions related to an environmentally friendly approach for the control of H. halys in invaded areas.     

Effective population size as a driver for divergence of an antimicrobial peptide (Hymenoptaecin) in two common European bumblebee species

Social insects are the target of numerous pathogens. This is because the high density of closely‐related individuals frequently interacting with each other enhances the transmission and establishment of pathogens. This high selective pressure results in the rapid evolution of immune genes, which might be counteracted by a reduced effective population size (N_e) lowering the effectiveness of selection. We tested the effect of N_e on the evolutionary rate of an important immune gene for the antimicrobial peptide Hymenoptaecin in two common central European bumblebee species: Bombus terrestris and Bombus lapidarius. Both species are similar in their biology and are expected to be under similar selective pressures because pathogen prevalence does not differ between species. However, previous studies indicated a higher N_e in B. terrestris compared to B. lapidarius. We found high intraspecific variability in the coding sequence but low variability for silent polymorphisms in B. lapidarius. Estimates of long‐ and short‐term N_e were three‐ to four‐fold higher N_e in B. terrestris, although the species did not differ in census population sizes. The difference in N_e might result in less efficient selection and suboptimal adaptation of immune genes (e.g. hymenoptaecin) in B. lapidarius, and thus this species might become less resistant and more tolerant, turning into a superspreader of diseases.     

Detecting an elusive invasive species: a diagnostic PCR to detect Burmese python in Florida waters and an assessment of persistence of environmental DNA

Recent studies have demonstrated that detection of environmental DNA (eDNA) from aquatic vertebrates in water bodies is possible. The Burmese python, Python bivittatus, is a semi‐aquatic, invasive species in Florida where its elusive nature and cryptic coloration make its detection difficult. Our goal was to develop a diagnostic PCR to detect P. bivittatus from water‐borne eDNA, which could assist managers in monitoring this invasive species. First, we used captive P. bivittatus to determine whether reptilian DNA could be isolated and amplified from water samples. We also evaluated the efficacy of two DNA isolation methods and two DNA extraction kits commonly used in eDNA preparation. A fragment of the mitochondrial cytochrome b gene from P. bivittatus was detected in all water samples isolated with the sodium acetate precipitate and the QIAamp DNA Micro Kit. Next, we designed P. bivittatus‐specific primers and assessed the degradation rate of eDNA in water. Our primers did not amplify DNA from closely related species, and we found that P. bivittatus DNA was consistently detectable up to 96 h. Finally, we sampled water from six field sites in south Florida. Samples from five sites, where P. bivittatus has been observed, tested positive for eDNA. The final site was negative and had no prior documented evidence of P. bivittatus. This study shows P. bivittatus eDNA can be isolated from water samples; thus, this method is a new and promising technique for the management of invasive reptiles.     

The cephalic labial gland secretions of two socially parasitic bumblebees Bombus hyperboreus (Alpinobombus) and Bombus inexspectatus (Thoracobombus) question their inquiline strategy

Social parasitic Hymenopterans have evolved morphological, chemical, and behavioral adaptations to overcome the sophisticated recognition and defense systems of their social host to invade host nests and exploit their worker force. In bumblebees, social parasitism appeared in at least 3 subgenera independently: in the subgenus Psithyrus consisting entirely of parasitic species, in the subgenus Alpinobombus with Bombus hyperboreus, and in the subgenus Thoracobombus with B. inexspectatus. Cuckoo bumblebee males utilize species‐specific cephalic labial gland secretions for mating purposes that can impact their inquiline strategy. We performed cephalic labial gland secretions in B. hyperboreus, B. inexspectatus and their hosts. Males of both parasitic species exhibited high species specific levels of cephalic gland secretions, including different main compounds. Our results showed no chemical mimicry in the cephalic gland secretions between inquilines and their host and we did not identify the repellent compounds already known in other cuckoo bumblebees.     

Comparative analyses of plastid and AFLP data suggest different colonization history and asymmetric hybridization between Betula pubescens and B. nana

Birches (Betula spp.) hybridize readily, confounding genetic signatures of refugial isolation and postglacial migration. We aimed to distinguish hybridization from range‐shift processes in the two widespread and cold‐adapted species B. nana and B. pubescens, previously shown to share a similarly east–west‐structured variation in plastid DNA (pDNA). We sampled the two species throughout their ranges and included reference samples of five other Betula species and putative hybrids. We analysed 901 individual plants using mainly nuclear high‐resolution markers (amplified fragment length polymorphisms; AFLPs); a subset of 64 plants was also sequenced for two pDNA regions. Whereas the pDNA variation as expected was largely shared between B. nana and B. pubescens, the two species were distinctly differentiated at AFLP loci. In B. nana, both the AFLP and pDNA results corroborated the former pDNA‐based hypothesis that it expanded from at least two major refugia in Eurasia, one south of and one east of the North European ice sheets. In contrast, B. pubescens showed a striking lack of geographic structuring of its AFLP variation. We identified a weak but significant increase in nuclear (AFLP) gene flow from B. nana into B. pubescens with increasing latitude, suggesting hybridization has been most frequent at the postglacial expansion front of B. pubescens and that hybrids mainly backcrossed to B. pubescens. Incongruence between pDNA and AFLP variation in B. pubescens can be explained by efficient expansion from a single large refugium combined with leading‐edge hybridization and plastid capture from B. nana during colonization of new territory already occupied by this more cold‐tolerant species.     

DNA barcoding of Corydalis,the most taxonomically complicated genus of Papaveraceae

The genus Corydalis is recognized as one of the most taxonomically challenging plant taxa. It is mainly distributed in the Himalaya–Hengduan Mountains, a global biodiversity hotspot. To date, no effective solution for species discrimination and taxonomic assignment in Corydalis has been developed. In this study, five nuclear and chloroplast DNA regions, ITS, ITS2, matK, rbcL, and psbA‐trnH, were preliminarily assessed based on their ability to discriminate Corydalis to eliminate inefficient regions, and the three regions showing good performance (ITS, ITS2 and matK) were then evaluated in 131 samples representing 28 species of 11 sections of four subgenera in Corydalis using three analytical methods (NJ, ML, MP tree; K2P‐distance and BLAST). The results showed that the various approaches exhibit different species identification power and that BLAST shows the best performance among the tested approaches. A comparison of different barcodes indicated that among the single barcodes, ITS (65.2%) exhibited the highest identification success rate and that the combination of ITS + matK (69.6%) provided the highest species resolution among all single barcodes and their combinations. Three Pharmacopoeia‐recorded medicinal plants and their materia medica were identified successfully based on the ITS and ITS2 regions. In the phylogenetic analysis, the sections Thalictrifoliae, Sophorocapnos, Racemosae, Aulacostigma, and Corydalis formed well‐supported separate lineages. We thus hypothesize that the five sections should be classified as an independent subgenus and that the genus should be divided into three subgenera. In this study, DNA barcoding provided relatively high species discrimination power, indicating that it can be used for species discrimination in this taxonomically complicated genus and as a potential tool for the authentication of materia medica belonging to Corydalis.     

Evaluation of 11 single‐locus and seven multilocus DNA barcodes in Lamium L. (Lamiaceae)

The aim of this work was to evaluate the suitability of selected DNA regions in the barcoding of plants, based on the species belonging to the genus Lamium (Lamiaceae). For this purpose, nine chloroplast barcodes, that is, accD, matK, rbcL, rpoA, rpoB, rpoC1, rpoC2, trnH‐psbA, trnL‐trnF, as well as ITS nuclear region, and intron of mitochondrial nad5 gene were tested. Among the single‐locus barcodes, most effective in the identification of Lamium species was the trnH‐psbA spacer and matK gene. The high level of variability and resolving power was also observed in the case of rpoA and rpoC2 genes. Despite the high interspecies variability of ITS region, it turned out to be inapplicable in Lamium identification. An important disadvantage of ITS as a barcode is a limitation of its use in polyploid plants, samples contaminated with fungal material or samples with partially degraded DNA. We have also evaluated five‐two‐locus and two‐three‐locus barcode regions created from a combination of most effective single loci. The best‐performing barcode combinations were matK + trnH‐psbA and matK + rpoA. Both of them had equally high discriminative power to identify Lamium species.     

Vicariance of Pyrocoelia fireflies (Coleoptera: Lampyridae) in the Ryukyu islands,Japan

We performed molecular phylogenetic analyses based on the mitochondrial COI gene (687 bp) and the nuclear 28S rRNA gene (715 bp) and reconstructed phylogenetic trees of the Pyrocoelia fireflies in the Ryukyu Islands and eastern Asia. Age calibration was done using a robust geological constraint: the Okinawa trough and associated straits began to rift at 1.55 Ma, isolating the Ryukyu Islands from the Chinese continent, Japanese islands, Taiwan island and some of the islands from each other. We suggest that the physical isolation of these islands began to generate the allopatric speciation within these islands, so the timing of this isolation was assigned to an appropriate node. The topology is completely concordant among phylogenetic trees reconstructed using MEGA (maximum‐likelihood), raxmlGUI (maximum‐likelihood) and BEAST (Bayesian inference; including combined analysis of COI and 28S rRNA genes). Two lineages are recognized, related to their emergence time; spring to summer, and autumn. In each lineage, vicariance is inferred to have begun at 1.55 Ma from our phylogenetic and geological analyses. In lineage 1, P. oshimana (Amami), P. matsumurai (Okinawa), P. discicollis (W. Japan), P. fumosa (E. Japan) and P. abdominalis (Yaeyama) were differentiated. In lineage 2, P. rufa (Tsushima and Korea), P. miyako (Miyako‐jima), P. atripennis (Ishigaki‐jima) and P. praetexta (Taiwan and HongKong) were differentiated. Pyrocoelia analis (Taiwan and China) emerges throughout the year except for winter, and constitutes another lineage. We suggest that Pyrocoelia fireflies differentiated at 2 Ma to generate these three lineages. The base substitution rate for the COI gene is estimated as 4.48% Myr^–1 and that for the 28S rRNA gene is 0.394%  Myr^–1, and these rates were used in a combined BEAty analysis in BEAST.