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

Lineage Divergence of Dendrolimus punctatus in Southern China Based on Mitochondrial Genome

In southern China, the masson pine caterpillar, Dendrolimus punctatus, has caused serious damage to the Pinus massoniana (Lamb.) pine forests. Here, the whole mitochondrial DNA (mtDNA) was employed to analyze the population evolution of D. punctatus and to understand the process underlying its current phylogenetic pattern. D. punctatus populations within its distribution range in China were categorized into five subgroups: central and eastern China (CEC), southwestern China (SWC), Yibin in Sichuan (SC), Baise in Guangxi (GX), and Luoding in Guangdong (GD), with a high level of haplotype diversity and nucleotide diversity among them. The genetic distances between subgroups are relatively large; however, the genetic distances between populations within the CEC subgroup were relatively small, suggesting that many populations were closely related in this subgroup. The mantel test showed that geographic distance had an important impact on the genetic distance of different geographic populations (r = 0.3633, P < 0.001). The neutrality tests, Bayesian skyline plot, and haplotype network showed that D. punctatus experienced a population expansion around 100,000 years ago. The divergence times of GX/SC, SWC, GD, and CEC were 0.347, 0.236, 0.200, and 0.110 million years ago, respectively. The SWC, CEC, and GD subgroups might have evolved from GX or SC subgroups. The population genetic structure of D. punctatus was closely related to its host tree species, geographic distance among populations, the weak flight capacity, and many eco-environment conditions.

Bayesian Inference of Evolutionary Histories under Time-Dependent Substitution Rates

Many factors complicate the estimation of time scales for phylogenetic histories, requiring increasingly complex evolutionary models and inference procedures. The widespread application of molecular clock dating has led to the insight that evolutionary rate estimates may vary with the time frame of measurement. This is particularly well established for rapidly evolving viruses that can accumulate sequence divergence over years or even months. However, this rapid evolution stands at odds with a relatively high degree of conservation of viruses or endogenous virus elements over much longer time scales. Building on recent insights into time-dependent evolutionary rates, we develop a formal and flexible Bayesian statistical inference approach that accommodates rate variation through time. We evaluate the novel molecular clock model on a foamy virus cospeciation history and a lentivirus evolutionary history and compare the performance to other molecular clock models. For both virus examples, we estimate a similarly strong time-dependent effect that implies rates varying over four orders of magnitude. The application of an analogous codon substitution model does not implicate long-term purifying selection as the cause of this effect. However, selection does appear to affect divergence time estimates for the less deep evolutionary history of the Ebolavirus genus. Finally, we explore the application of our approach on woolly mammoth ancient DNA data, which shows a much weaker, but still important, time-dependent rate effect that has a noticeable impact on node age estimates. Future developments aimed at incorporating more complex evolutionary processes will further add to the broad applicability of our approach.

Phylogeography and population structure of the tsetse fly Glossina pallidipes in Kenya and the Serengeti ecosystem

Glossina pallidipes is the main vector of animal African trypanosomiasis and a potential vector of human African trypanosomiasis in eastern Africa where it poses a large economic burden and public health threat. Vector control efforts have succeeded in reducing infection rates, but recent resurgence in tsetse fly population density raises concerns that vector control programs require improved strategic planning over larger geographic and temporal scales. Detailed knowledge of population structure and dispersal patterns can provide the required information to improve planning. To this end, we investigated the phylogeography and population structure of G. pallidipes over a large spatial scale in Kenya and northern Tanzania using 11 microsatellite loci genotyped in 600 individuals. Our results indicate distinct genetic clusters east and west of the Great Rift Valley, and less distinct clustering of the northwest separate from the southwest (Serengeti ecosystem). Estimates of genetic differentiation and first-generation migration indicated high genetic connectivity within genetic clusters even across large geographic distances of more than 300 km in the east, but only occasional migration among clusters. Patterns of connectivity suggest isolation by distance across genetic breaks but not within genetic clusters, and imply a major role for river basins in facilitating gene flow in G. pallidipes. Effective population size (Ne) estimates and results from Approximate Bayesian Computation further support that there has been recent G. pallidipes population size fluctuations in the Serengeti ecosystem and the northwest during the last century, but also suggest that the full extent of differences in genetic diversity and population dynamics between the east and the west was established over evolutionary time periods (tentatively on the order of millions of years). Findings provide further support that the Serengeti ecosystem and northwestern Kenya represent independent tsetse populations. Additionally, we present evidence that three previously recognized populations (the Mbeere-Meru, Central Kenya and Coastal "fly belts") act as a single population and should be considered as a single unit in vector control.

The impact of Miocene orogeny for the diversification of Caucasian Epeorus (Caucasiron) mayflies (Ephemeroptera: Heptageniidae)

A common hypothesis for the high biodiversity of mountains is the diversification driven by orogeny creating conditions for rapid in situ speciation of resident lineages. The Caucasus is a young mountain system considered as a biodiversity hotspot; however, the origin and evolution of its diversity remain poorly understood. This study focuses on mayflies of the subgenus Caucasiron, one of the most diversified stenotopic mayflies inhabiting various types of streams throughout the Caucasus. Using the time-calibrated phylogeny based on two mitochondrial (COI, 16S) and three nuclear (EF-1α, wg, 28S) gene fragments, we tested the role of Caucasian orogeny in biogeography, diversification patterns, and altitudinal diversification of Caucasiron mayflies. We found that orogeny promoted the lineage diversification of Caucasiron in the Miocene. The highest diversification rate corresponding with the uplift of mountains was followed by a significant slowdown towards the present suggesting minor influence of Pleistocene climatic oscillations on the speciation. The Caucasiron lineages cluster into three principal clades originating in the Upper Miocene. We found a strong support that one of the three clades diversified via allopatric speciation in the Greater Caucasus isolated in the Parathetys Sea. The other two clades originating most likely outside the Greater Caucasus diversified towards high and low altitude, respectively, indicating possible role of climatic factors and/or passive uplift on their differentiation. Current high Caucasiron diversity in the Greater Caucasus is a result of in situ speciation and later immigration from adjacent mountain ranges after the Parathetys Sea retreat. Our phylogeny supported the monophyly of Rhithrogeninae, Epeorus s.l., Caucasiron, and Iron. Epeorus subgenus Ironopsis was found paraphyletic, with its European representatives more closely related to Epeorus s.str. than to Iron. Therefore, we re-arranged taxa treated within Ironopsis to comply with the phylogeny recovered herein.

Genealogical Position of Japanese Populations of the Globally Distributed Mayfly Cloeon dipterum and Related Species (Ephemeroptera, Baetidae): A Molecular Phylogeographic Analysis

In the present study, we add genetic data of the mayfly Cloeon dipterum collected from the Japanese Islands to the established molecular phylogenetic knowledge in the mitochondrial COI gene of Cloeon mayflies. Cloeon dipterum is a typical cosmopolitan species that includes six intraspecific haplotype groups. The present phylogenetic analysis revealed that haplotypes of the Japanese C. dipterum constitute a seventh group together with a haplotype from Korea. This East Asian group forms a sister group with previously known European and North American haplotype groups (i.e., the clade CT1 to CT3). The present phylogenetic analysis further revealed the occurrence of two described species (C. dipterum and Cloeon ryogokuensis) and possibly three species (Cloeon sp. 1 to 3) in Japan. Consideration is given to the degree of genetic differentiation, divergence time, and differentiation process among these seven genetic groups.

Origins, Divergence, and Contrasting Invasion History of the Sweet Potato Weevil Pests Cylas formicarius (Coleoptera: Brentidae) and Euscepes batatae (Coleoptera: Curculionidae) in the Asia-Pacific

Cylas formicarius F. and Euscepes batatae Waterhouse are the most damaging sweet potato insect pests globally. Both weevils are thought to have invaded the Pacific alongside the movement of sweet potato (Ipomoea batatas (L.) Lam. Convolvulaceae), with C. formicarius having originated in India and E. batatae in Central or South America. Here we compare the genetic relationships between populations of the pests, primarily in the Asia-Pacific, to understand better their contemporary population structure and their historical movement relative to that of sweet potato. Cylas formicarius has divergent mitochondrial lineages that indicate a more complex biogeographic and invasive history than is presently assumed for this insect, suggesting it was widespread across the Asia-Pacific before the arrival of sweet potato. Cylas formicarius must have originally fed on Ipomoea species other than I. batatas but the identity of these species is presently unknown. Cylas formicarius was formerly designated as three species or subspecies and the genetic data presented here suggests that these designations should be reinvestigated. Euscepes batatae has very low genetic diversity which is consistent with its historical association with sweet potato and a recent introduction to the Asia-Pacific from the Americas. The distribution of E. batatae may be narrower than that of C. formicarius in the Asia-Pacific because it has relied relatively more on human-assisted movement. Consequently, E. batatae may become more widespread in the future. Investigating the invasion history of both species will help to understand the probability and nature of future invasions.

Phylogeny, biogeography and systematics of the hyper-diverse blister beetle genus Hycleus (Coleoptera: Meloidae)

Hycleus is a hyper-diverse genus of blister beetles including ~500 species widely distributed in the Old World, currently divided into three "sections" and into 45 "phenetic" species groups according to morphological characters. Recently the monophyly of Hycleus was questioned pointing out its paraphyly with respect to the genera Ceroctis and Paractenodia. In this study, we built a time-calibrated phylogenetic tree based on DNA sequence data from mitochondrial and nuclear genes obtained from 125 species, to understand the phylogenetic relationships among the species of this genus, to infer the biogeographic processes behind their diversification, and to assess their taxonomy and classification. Our results identified four main lineages one of which included the species belonging to Ceroctis and Paractenodia; therefore, both taxa are now referred to Hycleus as new synonyms. The three described sections of Hycleus resulted polyphyletic and are rejected, whereas several species groups represented well supported clades. Hycleus likely originated in Africa during the Early Miocene (~20 Mya), and subsequently spread in Europe and western Asia. Later, in the Late Miocene (~6 Mya) a Saharo-Sindian group branched off from the Palaearctic lineage, whereas the Oriental Region was colonized following a dispersal event through the Arabian Peninsula from the Afrotropical Region (~5 Mya).

Variable Molecular Markers for the Order Mantophasmatodea (Insecta)

The recently discovered insect order Mantophasmatodea currently comprises 19 Southern African species. These mainly occur in allopatry, have high levels of color polymorphism and communicate via species- and gender-specific vibratory signals. High levels of interspecific morphological conservatism mean that cryptic species are likely to be uncovered. These aspects of Mantophasmatodean biology make them an ideal group in which to investigate population divergence due to habitat-specific adaptation, sexual selection, and potentially sensory speciation. Lack of appropriate genetic markers has thus far rendered such studies unfeasible. To address this need, the first microsatellite loci for this order were developed. Fifty polymorphic loci were designed specifically for Karoophasma biedouwense (Austrophasmatidae), out of which 23 were labeled and tested for amplification across the order using 2-3 individuals from 10 species, representing all 4 currently known families. A Bayesian mitochondrially encoded cytochrome c oxidase I (COI) topology was reconstructed and divergence dates within the order were estimated for the first time. Amplification success and levels of polymorphism were compared with genetic divergence and time since divergence. In agreement with studies on vertebrate taxa, both amplification and variability were negatively correlated with distance (temporal and genetic). The high number of informative loci will offer sufficient resolution for both broad level population genetic analysis and individual based pedigree or parentage analyses for most species in Austrophasmatidae, with at least some loci available for the other families. This resource will facilitate research into the evolutionary biology of this understudied but fascinating group.

Persistence of multiple patterns and intraspecific polymorphism in multi-species Müllerian communities of net-winged beetles

BACKGROUND

In contrast to traditional models of purifying selection and a single aposematic signal in Müllerian complexes, some communities of unprofitable prey contain members with multiple aposematic patterns. Processes responsible for diversity in aposematic signaling are poorly understood and large multi-species communities are seldom considered.

RESULTS

We analyzed the phylogeny and aposematic patterns of closely related Eniclases net-winged beetles in New Guinea using mtDNA and nextRAD data. We suggest three clades of closely related and incompletely reproductively isolated lineages, detail the extent of polymorphism among Eniclases, and categorize their low-contrast aposematic patterns. The warning signal of Eniclases consists of body shape and color, with ambiguous color perception under some circumstances, i.e., when resting on the undersides of leaves. Field observations suggest that perception of the aposematic signal is affected by beetle behavior and environmental conditions. Local communities containing Eniclases consisted of 7-85 metriorrhynchine species assigned to 3-10 colour patterns.

CONCLUSION

As a result, we suggest that under certain light conditions the aposematic colour signal is less apparent than the body shape in net-winged beetle communities. We document variable environmental factors in our study area and highly diverse multi-species communities of other net-winged beetles. Which implies dynamically changing community structure in space and time. Variable environmental conditions and diverse community composition are suggested to be favorable for the persistence of multiple aposematic patterns, imperfect mimics, and intraspecific polymorphism. Further research should identify the relative effect of these factors on purifying selection and the alleles which are responsible for phenotypic differences.

A demonstration of unsupervised machine learning in species delimitation

One major challenge to delimiting species with genetic data is successfully differentiating population structure from species-level divergence, an issue exacerbated in taxa inhabiting naturally fragmented habitats. Many fields of science are now using machine learning, and in evolutionary biology supervised machine learning has recently been used to infer species boundaries. These supervised methods require training data with associated labels. Conversely, unsupervised machine learning (UML) uses inherent data structure and does not require user-specified training labels, potentially providing more objectivity in species delimitation. In the context of integrative taxonomy, we demonstrate the utility of three UML approaches (random forests, variational autoencoders, t-distributed stochastic neighbor embedding) for species delimitation in an arachnid taxon with high population genetic structure (Opiliones, Laniatores, Metanonychus). We find that UML approaches successfully cluster samples according to species-level divergences and not high levels of population structure, while model-based validation methods severely over-split putative species. UML offers intuitive data visualization in two-dimensional space, the ability to accommodate various data types, and has potential in many areas of systematic and evolutionary biology. We argue that machine learning methods are ideally suited for species delimitation and may perform well in many natural systems and across taxa with diverse biological characteristics.

Divergent evolutionary trajectories following speciation in two ectoparasitic honey bee mites

Multispecies host-parasite evolution is common, but how parasites evolve after speciating remains poorly understood. Shared evolutionary history and physiology may propel species along similar evolutionary trajectories whereas pursuing different strategies can reduce competition. We test these scenarios in the economically important association between honey bees and ectoparasitic mites by sequencing the genomes of the sister mite species Varroa destructor and Varroa jacobsoni. These genomes were closely related, with 99.7% sequence identity. Among the 9,628 orthologous genes, 4.8% showed signs of positive selection in at least one species. Divergent selective trajectories were discovered in conserved chemosensory gene families (IGR, SNMP), and Halloween genes (CYP) involved in moulting and reproduction. However, there was little overlap in these gene sets and associated GO terms, indicating different selective regimes operating on each of the parasites. Based on our findings, we suggest that species-specific strategies may be needed to combat evolving parasite communities.

Phylogenetic divergence of island biotas: Molecular dates, extinction, and "relict" lineages

Island formation is a key driver of biological evolution, and several studies have used geological ages of islands to calibrate rates of DNA change. However, many islands are home to "relict" lineages whose divergence apparently pre-dates island age. The geologically dynamic New Zealand (NZ) archipelago sits upon the ancient, largely submerged continent Zealandia, and the origin and age of its distinctive biota have long been contentious. While some researchers have interpreted NZ's biota as equivalent to that of a post-Oligocene island, a recent review of genetic studies identified a sizeable proportion of pre-Oligocene "relict" lineages, concluding that much of the biota survived an incomplete drowning event. Here, we assemble comparable genetic divergence data sets for two recently formed South Pacific archipelagos (Lord Howe; Chatham Islands) and demonstrate similarly substantial proportions of relict lineages. Similar to the NZ biota, our island reviews provide surprisingly little evidence for major genetic divergence "pulses" associated with island emergence. The dominance of Quaternary divergence estimates in all three biotas may highlight the importance of rapid biological turnover and new arrivals in response to recent climatic and/or geological disturbance and change. We provide a schematic model to help account for discrepancies between expected versus observed divergence-date distributions for island biotas, incorporating the effects of both molecular dating error and lineage extinction. We conclude that oceanic islands can represent both evolutionary "cradles" and "museums" and that the presence of apparently archaic island lineages does not preclude dispersal origins.

Host and geography together drive early adaptive radiation of Hawaiian planthoppers

The interactions between insects and their plant host have been implicated in driving diversification of both players. Early arguments highlighted the role of ecological opportunity, with the idea that insects "escape and radiate" on new hosts, with subsequent hypotheses focusing on the interplay between host shifting and host tracking, coupled with isolation and fusion, in generating diversity. Because it is rarely possible to capture the initial stages of diversification, it is particularly difficult to ascertain the relative roles of geographic isolation versus host shifts in initiating the process. The current study examines genetic diversity between populations and hosts within a single species of endemic Hawaiian planthopper, Nesosydne umbratica (Hemiptera, Delphacidae). Given that the species was known as a host generalist occupying unrelated hosts, Clermontia (Campanulaceae) and Pipturus (Urticaceae), we set out to determine the relative importance of geography and host in structuring populations in the early stages of differentiation on the youngest islands of the Hawaiian chain. Results from extensive exon capture data showed that N. umbratica is highly structured, both by geography, with discrete populations on each volcano, and by host plant, with parallel radiations on Clermontia and Pipturus leading to extensive co-occurrence. The marked genetic structure suggests that populations can readily become established on novel hosts provided opportunity; subsequent adaptation allows monopolization of the new host. The results support the role of geographic isolation in structuring populations and with host shifts occurring as discrete events that facilitate subsequent parallel geographic range expansion.

Loss and gain of sexual reproduction in the same stick insect

The outcome of competition between different reproductive strategies within a single species can be used to infer selective advantage of the winning strategy. Where multiple populations have independently lost or gained sexual reproduction it is possible to investigate whether the advantage is contingent on local conditions. In the New Zealand stick insect Clitarchus hookeri, three populations are distinguished by recent change in reproductive strategy and we determine their likely origins. One parthenogenetic population has established in the United Kingdom and we provide evidence that sexual reproduction has been lost in this population. We identify the sexual population from which the parthenogenetic population was derived, but show that the UK females have a post‐mating barrier to fertilisation. We also demonstrate that two sexual populations have recently arisen in New Zealand within the natural range of the mtDNA lineage that otherwise characterizes parthenogenesis in this species. We infer independent origins of males at these two locations using microsatellite genotypes. In one population, a mixture of local and nonlocal alleles suggested males were the result of invasion. Males in another population were most probably the result of loss of an X chromosome that produced a male phenotype in situ. Two successful switches in reproductive strategy suggest local competitive advantage for outcrossing over parthenogenetic reproduction. Clitarchus hookeri provides remarkable evidence of repeated and rapid changes in reproductive strategy, with competitive outcomes dependent on local conditions.

Untangling the Hypogeococcus pungens species complex (Hemiptera: Pseudococcidae) for Argentina, Australia, and Puerto Rico based on host plant associations and genetic evidence

Hypogeococcus pungens, a mealybug native of southern South America, is devastating native cacti in Puerto Rico and threatening cactus diversity in the Caribbean, and potentially in Central and North America. The taxonomic status of H. pungens is controversial since it has been reported feeding not only on Cactaceae but also on other plant families throughout its distribution range. However, in Australia, where the species had been exported from Argentina to control weedy American cacti, it was never found on host plants other than Cactaceae. These conflicting pieces of evidence not only cast doubt on the species identity that invaded Puerto Rico, but also have a negative impact on the search for natural enemies to be used in biological control programs against this pest. Here we present reproductive incompatibility and phylogenetic evidences that give support to the hypothesis that H. pungens is a species complex in which divergence appears to be driven by the host plants. The nuclear EF1α and 18S and the mitochondrial COI genes were used as markers to evaluate the phylogenetic relationships among H. pungens populations collected in Argentina, Australia and Puerto Rico feeding on Cactaceae and/or Amaranthaceae. Additionally, we conducted reciprocal crosses between mealybugs from both hosts. Species delimitation analysis revealed two well-supported putative species within H. pungens, one including mealybugs feeding on Amaranthaceae (H. pungens sensu stricto), and a new undescribed species using Cactaceae as hosts. Additionally, we found asymmetric reproductive incompatibility between these putative species suggesting recent reproductive isolation. The Bayesian species delimitation also suggested that the Australian mealybug population may derive from another undescribed species. Overall, the patterns of genetic differentiation may be interpreted as the result of recent speciation events prompted by host plant shifts. Finally, the finding of a single haplotype in the Puerto Rico population suggests only one invasive event. We still need to identify the geographical origin of the pest in order to enable the use of biological control to reduce the threat to cacti diversity in the Caribbean.

Diatraea saccharalis history of colonization in the Americas. The case for human-mediated dispersal

The sugarcane borer moth, Diatraea saccharalis, is one of the most important pests of sugarcane and maize crops in the Western Hemisphere. The pest is widespread throughout South and Central America, the Caribbean region and the southern United States. One of the most intriguing features of D. saccharalis population dynamics is the high rate of range expansion reported in recent years. To shed light on the history of colonization of D. saccharalis, we investigated the genetic structure and diversity in American populations using single nucleotide polymorphism (SNPs) markers throughout the genome and sequences of the mitochondrial gene cytochrome oxidase (COI). Our primary goal was to propose possible dispersal routes from the putative center of origin that can explain the spatial pattern of genetic diversity. Our findings showed a clear correspondence between genetic structure and the geographical distributions of this pest insect on the American continents. The clustering analyses indicated three distinct groups: one composed of Brazilian populations, a second group composed of populations from El Salvador, Mexico, Texas and Louisiana and a third group composed of the Florida population. The predicted time of divergence predates the agriculture expansion period, but the pattern of distribution of haplotype diversity suggests that human-mediated movement was most likely the factor responsible for the widespread distribution in the Americas. The study of the early history of D. saccharalis promotes a better understanding of range expansion, the history of invasion, and demographic patterns of pest populations in the Americas.

Mitochondrial Genomics Reveals Shared Phylogeographic Patterns and Demographic History among Three Periodical Cicada Species Groups

The mass application of whole mitogenome (MG) sequencing has great potential for resolving complex phylogeographic patterns that cannot be resolved by partial mitogenomic sequences or nuclear markers. North American periodical cicadas (Magicicada) are well known for their periodical mass emergence at 17- and 13-year intervals in the north and south, respectively. Magicicada comprises three species groups, each containing one 17-year species and one or two 13-year species. Within each life cycle, single-aged cohorts, called broods, of periodical cicadas emerge in different years, and most broods contain members of all three species groups. There are 12 and three extant broods of 17- and 13-year cicadas, respectively. The phylogeographic relationships among the populations and broods within the species groups have not been clearly resolved. We analyzed 125 whole MG sequences from all broods and seven species within three species groups to ascertain the divergence history of the geographic and allochronic populations and their life cycles. Our mitogenomic phylogeny analysis clearly revealed that each of the three species groups had largely similar phylogeographic subdivisions (east, middle, and west) and demographic histories (rapid population expansion after the last glacial period). The mitogenomic phylogeny also partly resolved the brood diversification process, which could be explained by hypothetical temporary life cycle shifts, and showed that none of the 13- and 17-year species within the species groups was monophyletic, possibly due to gene flow between them. Our findings clearly reveal phylogeographic structures in the three Magicicada species groups, demonstrating the advantage of whole MG sequence data in phylogeographic studies.

A mirage of cryptic species: Genomics uncover striking mitonuclear discordance in the butterfly Thymelicus sylvestris

Mitochondrial DNA (mtDNA) sequencing has led to an unprecedented rise in the identification of cryptic species. However, it is widely acknowledged that nuclear DNA (nuDNA) sequence data are also necessary to properly define species boundaries. Next generation sequencing techniques provide a wealth of nuclear genomic data, which can be used to ascertain both the evolutionary history and taxonomic status of putative cryptic species. Here, we focus on the intriguing case of the butterfly Thymelicus sylvestris (Lepidoptera: Hesperiidae). We identified six deeply diverged mitochondrial lineages; three distributed all across Europe and found in sympatry, suggesting a potential case of cryptic species. We then sequenced these six lineages using double-digest restriction-site associated DNA sequencing (ddRADseq). Nuclear genomic loci contradicted mtDNA patterns and genotypes generally clustered according to geography, i.e., a pattern expected under the assumption of postglacial recolonization from different refugia. Further analyses indicated that this strong mtDNA/nuDNA discrepancy cannot be explained by incomplete lineage sorting, sex-biased asymmetries, NUMTs, natural selection, introgression or Wolbachia-mediated genetic sweeps. We suggest that this mitonuclear discordance was caused by long periods of geographic isolation followed by range expansions, homogenizing the nuclear but not the mitochondrial genome. These results highlight T. sylvestris as a potential case of multiple despeciation and/or lineage fusion events. We finally argue, since mtDNA and nuDNA do not necessarily follow the same mechanisms of evolution, their respective evolutionary history reflects complementary aspects of past demographic and biogeographic events.

Herbivore range expansion triggers adaptation in a subsequently-associated third trophic level species and shared microbial symbionts

Invasive species may change the life history strategies, distribution, genetic configuration and trophic interactions of native species. The diamondback moth, Plutella xylostella L., is an invasive herbivore attacking cultivated and wild brassica plants worldwide. Here we present phylogeographic analyses of P. xylostella and one of its major parasitoids, Cotesia vestalis, using mitochondrial markers, revealing the genetic diversity and evolutionary history of these two species. We find evidence that C. vestalis originated in Southwest China, then adapted to P. xylostella as a new host by ecological sorting as P. xylostella expanded its geographic range into this region. Associated with the expansion of P. xylostella, Wolbachia symbionts were introduced into local populations of the parasitoid through horizontal transfer from its newly associated host. Insights into the evolutionary history and phylogeographic system of the herbivore and its parasitoid provide an important basis for better understanding the impacts of biological invasion on genetic configuration of local species.

Evolutionary persistence of insect bunyavirus infection despite host acquisition and expression of the viral nucleoprotein gene

How insects combat RNA virus infection is a subject of intensive research owing to its importance in insect health, virus evolution, and disease transmission. In recent years, a pair of potentially linked phenomena have come to light as a result of this work-first, the pervasive production of viral DNA from exogenous nonretroviral RNA in infected individuals, and second, the widespread distribution of nonretroviral integrated RNA virus sequences (NIRVs) in the genomes of diverse eukaryotes. The evolutionary consequences of NIRVs for viruses are unclear and the field would benefit from studies of natural virus infections co-occurring with recent integrations, an exceedingly rare circumstance in the literature. Here, we provide evidence that a novel insect-infecting phasmavirus (Order Bunyavirales) has been persisting in a phantom midge host, Chaoborus americanus, for millions of years. Interestingly, the infection persists despite the host's acquisition (during the Pliocene), fixation, and expression of the viral nucleoprotein gene. We show that virus prevalence and geographic distribution are high and broad, comparable to the host-specific infections reported in other phantom midges. Short-read mapping analyses identified a lower abundance of the nucleoprotein-encoding genome segment in this virus relative to related viruses. Finally, the novel virus has facilitated the first substitution rate estimation for insect-infecting phasmaviruses. Over a period of approximately 16 million years, we find rates of (0.6 - 1.6) × 10^-7 substitutions per site per year in protein coding genes, extraordinarily low for negative-sense RNA viruses, but consistent with the few estimates produced over comparable evolutionary timescales.

Phylogeography and species distribution modelling of Cryptocephalusbarii (Coleoptera: Chrysomelidae): is this alpine endemic species close to extinction?

The alternation of glacial and interglacial cycles of the Quaternary period contributed in shaping the current species distribution. Cold-adapted organisms experienced range expansion and contraction in response to the temperature decrease and increase, respectively. In this study, a fragment of the mitochondrial marker COI was used to investigate the phylogeography of Cryptocephalusbarii, a cold-adapted alpine leaf beetle species endemic of Orobie Alps, northern Italy. The relationships among populations, their divergence time, and the most probable migration model were estimated and are discussed in light of the Pleistocene climate oscillations. Through a species distribution modelling analysis, the current habitat suitability was assessed and the distribution in a future global warming scenario predicted. The main divergence events that led to the actual population structure took place from ~750,000 to ~150,000 years ago, almost following the pattern of the climate oscillations that led to the increase of the connections between the populations during cold periods and the isolation on massifs in warm periods. The most supported migration model suggests that the species survived to past adverse climatic conditions within refugia inside and at the limit of the actual range. The species distribution modelling analysis showed that C.barii is extremely sensitive to air temperature variations, thus the increase of temperature caused by global warming will reduce the suitable areas within the species range, leading to its possible extinction in the next 50 years. Cryptocephalusbarii is a representative case of how cold adapted and limited distributed species have been and could be affected by climate change, that highlights the implementation of conservation actions.

Restriction-site associated DNA markers provide new insights into the evolutionary history of the bark beetle genus Dendroctonus

The bark beetle genus Dendroctonus contains some of the most economically important pests of conifers worldwide. Despite many attempts, there is no agreement today on the phylogenetic relationships within the genus, which limits our understanding of its evolutionary history. Here, using restriction-site associated DNA (RAD) markers from 70 specimens representing 17 species (85% of the known diversity) we inferred the phylogeny of the genus, its time of origin and biogeographic history, as well as the evolution of key ecological traits (host plants, larval behavior and adults' attack strategies). For all combinations of tested parameters (from 6444 to 23,570 RAD tags analyzed), the same, fully resolved topology was inferred. Our analyses suggest that the most recent common ancestor (mrca) of all extant Dendroctonus species was widely distributed across eastern Palearctic and western Nearctic during the early Miocene, from where species dispersed to other Holarctic regions. A first main inter-continental vicariance event occurred during early Miocene isolating the ancestors of D. armandi in the Palearctic, which was followed by the radiation of the main Dendroctonus lineages in North America. During the Late Miocene, the ancestor of the 'rufipennis' species group colonized north-east Palearctic regions from western North America, which was followed by a second main inter-continental vicariance event isolating Pleistocene populations in Asia (D. micans) and western North America (D. murrayanae and D. punctatus). The present study supports previous hypotheses explaining intercontinental range disjunctions across the Northern Hemisphere by the fragmentation of a continuous distribution due to climatic cooling, host range fragmentation and geological changes during the late Cenozoic. The reconstruction of ancestral ecological traits indicates that the mrca bored individual galleries and mass attacked the boles of pines. The gregarious feeding behavior of the larvae as well as the individual attack of the base of trees have apparently independently evolved twice in North America (in the 'rufipennis' and the 'valens' species groups), which suggests a higher adaptive potential than previously thought and may be of interest for plant protection and biodiversity conservation in a rapidly changing world.

Molecular phylogeny reveals the gradual evolutionary transition to soft-bodiedness in click-beetles and identifies sub-Saharan Africa as a cradle of diversity for Drilini (Coleoptera: Elateridae)

Drilini are soft-bodied predatory click-beetles (Elateridae: Agrypninae) with incompletely metamorphosed females. Due to divergent morphology, their classification has been contentious. We present the first densely sampled molecular phylogeny of Drilini based on nuclear and mitochondrial markers. Altogether, 44 species, representing all genera, were analysed using maximum likelihood and the Bayesian approach. Molecular analyses recovered five major clades that were also well supported by morphology. Afrotropical lineages mark deep splits. Most Palearctic species belong to a terminal clade. A few species of the predominantly Afrotropical Selasia are distributed from Arabia up to the Himalayas and Thailand. The origin of Drilini is dated to the Late Eocene (~35.5 Mya) and rapid radiation is identified from the Eocene/Oligocene to the Middle Miocene. We describe the gradual transformation of male morphological traits, e.g. the level of sclerotization, structure of mouthparts, loss of the thoracic interlocking mechanism, shortened elytra and expanded larviform abdomen. Five new genera, Austroselasia, Habeshaselasia, Illubaboria, Malacodrilus and Mashaselasia are proposed. Latoselasia, previously a subgenus of Wittmerselasia, is elevated to the genus rank. Five new species, Habeshaselasia nekemtensis, H. illubaborensis, Illubaboria bicolor, Malacodrilus hajeki and Mashaselasia aethiopica are described.

Molecular phylogeny of Entomobrya (Collembola: Entomobryidae) from China: Color pattern groups and multiple origins

Highly diversified colorations among springtails (Collembola) have been widely used for species diagnosis, but their phylogenetic significance is poorly known. We addressed this issue in the largest Entomobryinae genus Entomobrya, which possesses variable color patterns among species. The relationships within the genus and to other genera have also rarely been studied. Based on material mainly from China, we have conducted a multilocus phylogeny and topology tests with likelihood and Bayesian algorithms, and accordingly demonstrated the non-monophyly of Chinese Entomobrya. The division of five clades, including Entomobrya and several related genera, coincided well with five types of colorations, respectively. Further analyses of divergence time and historical biogeography revealed that Chinese Entomobrya originated mainly from Palearctic (northern and western) China in the Paleocene and Eocene. This study highlights the great phylogenetic values as well as taxonomic uses of coloration in Chinese Entomobrya. Multiple phylogenetic and biogeographic origins of Entomobrya imply its complicated relationships with both scaled and unscaled genera of Entomobryinae.

Genetic basis of species-specific genitalia reveals role in species diversification

The diversity of genital morphology among closely related animals with internal fertilization is well known, but the genetic backgrounds are unclear. Here, we show that, in Carabus (Ohomopterus) beetles showing correlated evolution of male and female genital parts, only a few major quantitative trait loci (QTLs) determine differences in genital dimensions between sister species, and sequence divergence is pronounced in the genomic regions containing genital QTLs. The major QTLs for male and female genital dimensions reside in different locations within the same linkage group, implying that coevolution between the sexes is only loosely constrained and can respond to sexually antagonistic selection. The same genomic regions containing the major QTLs show elevated divergence between three pairs of parapatric species with marked differences in genital parts. Our study demonstrates that species diversification can follow coevolution of genitalia between the sexes, even without tight linkage of loci affecting male and female genital dimensions.

Phylogeographic structures of the host insects of Ophiocordyceps sinensis

A fungus-insect complex, known as DongChong XiaCao, is formed from the infection of the hepialid larvae by the fungus Ophiocordyceps sinensis, which is endemic to the Qinghai-Tibetan Plateau (QTP). Due to previously limited sample collection size, the data about the diversity and structure of the host insect was insufficient and lacked details. The purpose of this study was aimed to discuss the diversity and phylogeography of the host insects of O. sinensis with a large-scale sampling. The mitochondrial cytochrome oxidase I gene (cox1) was sequenced and analyzed among 710 samples representing 88 geographic locations. 205 haplotypes of cox1 were identified from all the 710 samples and 4 phylogenetic clades with 12 subclades were identified. Instead of the single latitude-based divergence suggested previously, three distribution patterns were deduced to correspond to the phylogeographic structures, including but not limited to the co-existence of a wide and specific local phylogeographic distribution structures. Two separate genetic diversity and differentiation centers, namely the northwestern Yunnan and the southeastern Tibet were identified. Dating analyses from three calibrations supported that the divergence of the 4 clades occurred in the Oligocene-Miocene period (30.54-13.66 million years ago) (Ma), which were connected with the second and third geological movements of the QTP (17-25, 8-13 Ma). Our results provide a more detailed understanding of the divergence and distribution patterns of the host insects of O. sinensis.

Towards a global DNA barcode reference library for quarantine identifications of lepidopteran stemborers, with an emphasis on sugarcane pests

Lepidopteran stemborers are among the most damaging agricultural pests worldwide, able to reduce crop yields by up to 40%. Sugarcane is the world’s most prolific crop, and several stemborer species from the families Noctuidae, Tortricidae, Crambidae and Pyralidae attack sugarcane. Australia is currently free of the most damaging stemborers, but biosecurity efforts are hampered by the difficulty in morphologically distinguishing stemborer species. Here we assess the utility of DNA barcoding in identifying stemborer pest species. We review the current state of the COI barcode sequence library for sugarcane stemborers, assembling a dataset of 1297 sequences from 64 species. Sequences were from specimens collected and identified in this study, downloaded from BOLD or requested from other authors. We performed species delimitation analyses to assess species diversity and the effectiveness of barcoding in this group. Seven species exhibited <0.03 K2P interspecific diversity, indicating that diagnostic barcoding will work well in most of the studied taxa. We identified 24 instances of identification errors in the online database, which has hampered unambiguous stemborer identification using barcodes. Instances of very high within-species diversity indicate that nuclear markers (e.g. 18S, 28S) and additional morphological data (genitalia dissection of all lineages) are needed to confirm species boundaries.

Phylogeographic structure and ecological niche modelling reveal signals of isolation and postglacial colonisation in the European stag beetle

Lucanus cervus (L.), the stag beetle, is a saproxylic beetle species distributed widely across Europe. Throughout its distribution the species has exhibited pronounced declines and is widely considered threatened. Conservation efforts may be hindered by the lack of population genetic data and understanding of the spatial scale of population connectivity. To address this knowledge gap this research details the first broad scale phylogeographic study of L. cervus based on mitochondrial DNA (mtDNA) sequencing and microsatellite analysis of samples collected from 121 localities across Europe. Genetic data were complemented by palaeo-distribution models of spatial occupancy during the Last Glacial Maximum to strengthen inferences of refugial areas. A salient feature of the mtDNA was the identification of two lineages. Lineage I was widespread across Europe while lineage II was confined to Greece. Microsatellites supported the differentiation of the Greek samples and alongside palaeo-distribution models indicated this area was a glacial refuge. The genetic endemism of the Greek samples, and demographic results compatible with no signatures of spatial expansion likely reflects restricted dispersal into and out of the area. Lineage I exhibited a shallow star like phylogeny compatible with rapid population expansion across Europe. Demographic analysis indicated such expansions occurred after the Last Glacial Maximum. Nuclear diversity and hindcast species distribution models indicated a central Italian refuge for lineage I. Palaeo-distribution modelling results also suggested a western refuge in northern Iberia and south-west France. In conclusion the results provide evidence of glacial divergence in stag beetle while also suggesting high, at least on evolutionary timescales, gene flow across most of Europe. The data also provide a neutral genetic framework against which patterns of phenotypic variation may be assessed.

The Malacca Strait separates distinct faunas of poorly-flying Cautires net-winged beetles

We investigated the spatial and temporal patterns of Cautires diversification on the Malay Peninsula and Sumatra to understand if the narrow and frequently dry Malacca Strait separates different faunas. Moreover, we analyzed the origin of Cautires in Malayan and Sumatran mountains. We sampled 18 localities and present the mtDNA-based phylogeny of 76 species represented by 388 individuals. The phylogenetic tree was dated using mtDNA evolution rates and the ancestral ranges were estimated using the maximum likelihood approach. The phylogeny identified multiple lineages on the Malay Peninsula since the Upper Eocene (35 million years ago, mya) and a delayed evolution of diversity in Sumatra since the Upper Oligocene (26 mya). A limited number of colonization events across the Malacca Strait was identified up to the Pliocene and more intensive faunal exchange since the Pleistocene. The early colonization events were commonly followed by in situ diversification. As a result, the Malacca Strait now separates two faunas with a high species-level turnover. The montane fauna diversified in a limited space and seldom took part in colonization events across the Strait. Besides isolation by open sea or a savannah corridor, mimetic patterns could decrease the colonization capacity of Cautires. The Malay fauna is phylogenetically more diverse and has a higher value if conservation priorities should be defined.

Recalibrating the molecular clock for Arctic marine invertebrates based on DNA barcodes

Divergence times for species assemblages of Arctic marine invertebrates have often been estimated using a standard rate (1.4%/MY) of molecular evolution calibrated using a single sister pair of tropical crustaceans. Because rates of molecular evolution vary among taxa and environments, it is essential to obtain clock calibrations from northern lineages. The recurrent opening and closure of the Bering Strait provide an exceptional opportunity for clock calibration. Here, we apply the iterative calibration approach to investigate patterns of molecular divergence among lineages of northern marine molluscs and arthropods using publicly available sequences of the cytochrome c oxidase subunit I (COI) gene and compare these results with previous estimates of trans-Bering divergences for echinoderms and polychaetes. The wide range of Kimura two-parameter (K2P) divergences among 73 trans-Bering sister pairs (0.12%–16.89%) supports multiple pulses of migration through the Strait. Overall, the results indicate a rate of K2P divergence of 3.2%/MY in molluscs, 5%–5.2%/MY in arthropods, and 3.5%–4.7%/MY in polychaetes. While these rates are considerably higher than the often-adopted 1.4%/MY rate, they are similar to calibrations (3%–5%/MY) in several other studies of marine invertebrates. This upward revision in rates means there is a need both to reevaluate the evolutionary history of marine lineages and to reexamine the impact of prior climatic changes upon the diversification of marine life.

Molecular phylogeny and diversification timing of the Nemouridae family (Insecta, Plecoptera) in the Japanese Archipelago

The generation of the high species diversity of insects in Japan was profoundly influenced by the formation of the Japanese Archipelago. We explored the species diversification and biogeographical history of the Nemouridae Billberg, 1820 family in the Japanese Archipelago using mitochondrial DNA and nuclear DNA markers. We collected 49 species among four genera: Indonemoura Baumann, 1975; Protonemura Kempny, 1898; Amphinemura, Ris 1902 and Nemoura Latreille, 1796 in Japan, China, South Korea and North America. We estimated their divergence times-based on three molecular clock node calibrations-using Bayesian phylogeography approaches. Our results suggested that Japanese Archipelago formation events resulted in diversification events in the middle of the Cretaceous (<120 Ma), speciation in the Paleogene (<50 Ma) and intra-species diversification segregated into eastern and western Japan of the Fossa Magna region at late Neogene (20 Ma). The Indonemoura samples were genetically separated into two clades-that of Mainland China and that of Japan. The Japanese clade clustered with the Nemouridae species from North America, suggesting the possibility of a colonisation event prior to the formation of the Japanese Archipelago. We believe that our results enhanced the understanding both of the origin of the species and of local species distribution in the Japanese Archipelago.

Genetic divergence, population differentiation and phylogeography of the cicada Subpsaltria yangi based on molecular and acoustic data: an example of the early stage of speciation?

Background

Geographical isolation combined with historical climatic fluctuations have been identified as two major factors that contribute to the formation of new species. On the other hand, biotic factors such as competition and predation are also able to drive the evolution and diversification of organisms. To determine whether geographical barriers contributed to population divergence or speciation in the rare endemic cicada Subpsaltria yangi the population differentiation, genetic structure and phylogeography of the species were investigated in the Loess Plateau and adjacent areas of northwestern China by analysing mitochondrial and nuclear DNA and comparing the calling song structure of 161 male individuals.

Results

The results reveal a low level of genetic differentiation and relatively simple phylogeographic structure for this species, but two independent clades corresponding to geographically isolated populations were recognised. Genetic and geographical distances were significantly correlated among lineages. Results of divergence-time estimation are consistent with a scenario of isolation due to glacial refugia and interglacial climate oscillation in northwestern China. Significant genetic divergence was found between the population occurring in the Helan Mountains and other populations, and recent population expansion has occurred in the Helan Mountains and/or adjacent areas. This population is also significantly different in calling song structure from other populations.

Conclusions

Geographical barriers (i.e., the deserts and semi-deserts surrounding the Helan Mountains), possibly coupled with related ecological differences, may have driven population divergence and allopatric speciation. This provides a possible example of incipient speciation in Cicadidae, improves understanding of population differentiation, acoustic signal diversification and phylogeographic relationships of this rare cicada species of conservation concern, and informs future studies on population differentiation, speciation and phylogeography of other insects with a high degree of endemism in the Helan Mountains and adjacent areas.

Electronic supplementary material

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

Game of clones: Is Wolbachia inducing speciation in a weevil with a mixed reproductive mode?

Parthenogenesis is widely distributed in Metazoa but it is especially frequent in weevils (Coleoptera, Curculionidae) with one fifth of all known cases. Previous studies have shown that in the tribe Naupactini parthenogenetic reproduction most likely originated with an infection of the endoparasitic bacterium Wolbachia pipientis. In particular, Pantomorus postfasciatus possess a mixed reproductive mode: some populations have males while in others they are absent, and females produce clones by thelytoky. To better understand this scenario, we studied the population structure and infection status in 64 individuals of P. postfasciatus from Argentina and Brazil. We sequenced two mitochondrial (COI and COII) and one nuclear (ITS-1) fragments and obtained two very divergent haplogroups, one corresponding to the sexual populations uninfected with Wolbachia, and another conforming a monophyletic parthenogenetic (or presumptively parthenogenetic) and infected clade. Each of these haplogroups was identified as an independently evolutionary unit by all species delimitation analyses accomplished: multilocus *BEAST and BP&P, and single locus GMYC and K/θ rule. Additionally, present evidence suggests that Wolbachia infection occurred at least twice in all-female populations of P. postfasciatus with two different bacterial strains. Speciation mediated by Wolbachia is a recently described phenomenon and the case of P. postfasciatus is the first known case in a diplo-diploid insect. A model that describes how thelytoky-inducing phenotypes of Wolbachia could generate new lineages is discussed.

Assessing Genomic Admixture between Cryptic Plutella Moth Species following Secondary Contact

Cryptic species are genetically distinct taxa without obvious variation in morphology and are occasionally discovered using molecular or sequence data sets of populations previously thought to be a single species. The world-wide Brassica pest, Plutella xylostella (diamondback moth), has been a problematic insect in Australia since 1882, yet a morphologically cryptic species with apparent endemism (P. australiana) was only recognized in 2013. Plutella xylostella and P. australiana are able to hybridize under laboratory conditions, and it was unknown whether introgression of adaptive traits could occur in the field to improve fitness and potentially increase pressure on agriculture. Phylogenetic reconstruction of 29 nuclear genomes confirmed P. xylostella and P. australiana are divergent, and molecular dating with 13 mitochondrial genes estimated a common Plutella ancestor 1.96 ± 0.175 Ma. Sympatric Australian populations and allopatric Hawaiian P. xylostella populations were used to test whether neutral or adaptive introgression had occurred between the two Australian species. We used three approaches to test for genomic admixture in empirical and simulated data sets including 1) the f3 statistic at the level of the population, 2) pairwise comparisons of Nei's absolute genetic divergence (dXY) between populations, and 3) changes in phylogenetic branch lengths between individuals across 50-kb genomic windows. These complementary approaches all supported reproductive isolation of the Plutella species in Australia, despite their ability to hybridize. Finally, we highlight the most divergent genomic regions between the two cryptic Plutella species and find they contain genes involved with processes including digestion, detoxification, and DNA binding.