Cryptic diversity in the pen shell Atrina pectinata (Bivalvia: Pinnidae): high divergence and hybridization revealed by molecular and morphological data

Cryptic species have been increasingly revealed in the marine realm through an analytical approach incorporating multiple lines of evidence (e.g., mtDNA, nuclear genes and morphology). Illustrations of cryptic taxa improve our understanding of species diversity and evolutionary histories within marine animals. The pen shell Atrina pectinata is known to exhibit extensive morphological variations that may harbour cryptic diversity. In this study, we investigated A. pectinata populations along the coast of China and one from Japan to explore possible cryptic diversity and hybridization using a combination of mitochondrial (cytochrome c oxidase subunit I, mtCOI) and nuclear (ribosomal internal transcribed spacer, nrITS) genes as well as morphology. Phylogenetic analyses of mtCOI ‘DNA barcoding gene’ sequences resolved six divergent lineages with intralineage divergences between 0.4% and 0.8%. Interlineage sequence differences ranged from 4.3% to 22.0%, suggesting that six candidate cryptic species are present. The nrITS gene revealed five deep lineages with Kimura 2‐parameter distances of 3.7–30.3%. The five nuclear lineages generally corresponded to mtCOI lineages 1–4 and (5 + 6), suggestive of five distinct evolutionary lineages. Multiple nrITS sequences of significant variance were found within an individual, clearly implying recent hybridization events between/among the evolutionary lineages, which contributed to cytonuclear discordance. Morphologically, five morphotypes matched the five genetic lineages, although the intermediates may well blur the boundaries of different morphotypes. This study demonstrates the importance of combining multiple lines of evidence to explore species cryptic diversity and past evolutionary histories.     

Genetic uniformity and long‐distance clonal dispersal in the invasive androgenetic Corbicula clams

The clam genus Corbicula is an interesting model system to study the evolution of reproductive modes as it includes both sexual and asexual (androgenetic) lineages. While the sexual populations are restricted to the native Asian areas, the androgenetic lineages are widely distributed being also found in America and Europe where they form a major aquatic invasive pest. We investigated the genetic diversity of native and invasive Corbicula populations through a worldwide sampling. The use of mitochondrial and nuclear (microsatellite) markers revealed an extremely low diversity in the invasive populations with only four, undiversified, genetic lineages distributed across Europe and America. On the contrary, in the native populations, both sexual and androgenetic lineages exhibited much higher genetic diversity. Remarkably, the most abundant and widely distributed invasive forms, the so‐called form A and form R found in America and Europe respectively, are fixed for the same single COI (cytochrome c oxydase subunit I) haplotype and same multilocus genotype. This suggests that form R, observed in Europe since the 1980s, derived directly from form A found in America since the 1920s. In addition, this form shares alleles with some Japanese populations, indicating a Japanese origin for this invasive lineage. Finally, our study suggests that few androgenetic Corbicula individuals successfully invaded the non‐native range and then dispersed clonally. This is one striking case of genetic paradox raising the issue of invasive and evolutionary success of genetically undiversified populations.     

Multilocus phylogeny and systematics of Iberian endemic Squalius (Actinopterygii,Leuciscidae)

Inferring the evolutionary history of a group of species can be challenging given the many factors involved. In recent years, the increased availability of sequences of multiple genes per species has spurred the development of new methodologies to analyse multilocus data sets. Two approaches that analyse such data are concatenated supermatrix and coalescent-based species-tree analyses. In this study, we used both of these methods to infer the phylogenetic relationships of Iberian species of the genus Squalius from one mitochondrial and six nuclear genes. We found mitonuclear discordance in the phylogenetic relationships of the group. According to the mitochondrial gene analysis, all species were recovered as monophyletic except S. pyrenaicus; besides, in the concatenated supermatrix analysis of the nuclear markers, this species resolved as polyphyletic with three divergent evolutionary lineages. The coalescent-based nuclear species-tree analysis rendered a well-resolved phylogeny compared with the supermatrix analysis, which was unable to discern between S. carolitertii, S. castellanus and one of the evolutionary lineages of S. pyrenaicus. This result is likely due to the better integration of population uncertainty in the coalescent approach. Furthermore, Bayesian multilocus species delimitation analyses based on a BPP approach strongly supported the distinct nuclear lineages as different species. Nevertheless, the supermatrix analysis was able to obtain well-supported relationships in the divergent lineages with low numbers of individuals. Our study highlights the usefulness of different analytical methodologies to obtain a more complete picture of the evolutionary history of taxa, especially when discordant patterns among genes are found.     

PERCHED AT THE MITO‐NUCLEAR CROSSROADS: DIVERGENT MITOCHONDRIAL LINEAGES CORRELATE WITH ENVIRONMENT IN THE FACE OF ONGOING NUCLEAR GENE FLOW IN AN AUSTRALIAN BIRD

Relationships among multilocus genetic variation, geography, and environment can reveal how evolutionary processes affect genomes. We examined the evolution of an Australian bird, the eastern yellow robin Eopsaltria australis, using mitochondrial (mtDNA) and nuclear (nDNA) genetic markers, and bioclimatic variables. In southeastern Australia, two divergent mtDNA lineages occur east and west of the Great Dividing Range, perpendicular to latitudinal nDNA structure. We evaluated alternative scenarios to explain this striking discordance in landscape genetic patterning. Stochastic mtDNA lineage sorting can be rejected because the mtDNA lineages are essentially distinct geographically for > 1500 km. Vicariance is unlikely: the Great Dividing Range is neither a current barrier nor was it at the Last Glacial Maximum according to species distribution modeling; nuclear gene flow inferred from coalescent analysis affirms this. Female philopatry contradicts known female‐biased dispersal. Contrasting mtDNA and nDNA demographies indicate their evolutionary histories are decoupled. Distance‐based redundancy analysis, in which environmental temperatures explain mtDNA variance above that explained by geographic position and isolation‐by‐distance, favors a nonneutral explanation for mitochondrial phylogeographic patterning. Thus, observed mito‐nuclear discordance accords with environmental selection on a female‐linked trait, such as mtDNA, mtDNA–nDNA interactions or genes on W‐chromosome, driving mitochondrial divergence in the presence of nuclear gene flow.     

Genetic and chaetal variation in Nais worms (Annelida,Clitellata, Naididae)

The genus Nais is a group of oligochaetous clitellates, common in eutrophic freshwater habitats. About 30 species are described. Species identification is based primarily on chaetal characters, which are often subtle, inconsistent, and even overlapping between nominal species. We investigated the correlation between genetic variation and chaetal morphology in this genus. Eighty‐one individuals from Europe, North America, and China were included in the study. Seventy‐five of these were preserved as vouchers. They were scrutinized with regard to chaetal morphology, and ten different morphotypes were identified. Three molecular markers, two mitochondrial (the COI gene and 16S rDNA) and one nuclear (the ITS region), were used to establish the genetic lineages in the material. Genetic variation was found to be largely congruent with chaetal character patterns. However, at least nine separately evolving lineages (all supported by mitochondrial as well as nuclear data) correspond to at most six nominal species. Four morphotypes/lineages are recognized as Nais barbata, Nais christinae, Nais elinguis, and Nais stolci, respectively, whereas five, or possibly more, lineages represent a morphological continuum covering the variation of the Nais communis/variabilis complex. Thus, cryptic speciation is revealed. Our results indicate that a taxonomic revision of the genus will be needed in the future.     

Nuclear genes (but not mitochondrial DNA barcodes) reveal real species: Evidence from the Brenthis fritillary butterflies (Lepidoptera,Nymphalidae)

Discordance between entities revealed by nuclear versus mitochondrial genes is a common phenomenon in evolutionary and taxonomic studies. However, little attention has been paid to analysis of how such discordant entities correspond to traditional species detected through investigation of their morphology, ecology, and distribution. Here, we used one mitochondrial (COI, DNA barcode fragment) and four nuclear (CAD, Ca‐ATPase, arginine kinase, wg) genes to analyze the genetic structure of the taxonomically well‐studied butterfly genus Brenthis (Lepidoptera, Nymphalidae). Analysis of COI revealed multiple diverged allopatric and sympatric mitochondrial lineages within the known Brenthis species hinting at possible presence of unrecognized cryptic species. However, these multiple‐species hypotheses were not supported by further studies of nuclear genes and phenotypic traits. The discovered mitochondrial lineages did not correspond to the clusters revealed by nuclear genes. Simultaneously, we found a complete congruence between (a) traditional species boundaries, (b) clusters recognized by nuclear genes, and (c) clusters identified via cladistic analysis of phenotypic traits (genitalia and wing pattern characters, ecological preferences, and chromosome numbers). We conclude that in case of the genus Brenthis, nuclear genes rather than mtDNA barcodes reveal real species boundaries. Additionally, we suggest to support each DNA barcode‐based taxonomic conclusion by analysis of phased alleles of nuclear genes, avoiding widely used practice of nuclear and mitochondrial genes concatenation without any examination of interaction of these different types of data.     

Phylogenetic analysis,taxonomic revision,and dental ontogeny of the Cretaceous Zhelestidae (Mammalia: Eutheria)

Iphisa elegans Gray, 1851 is a ground‐dwelling lizard widespread over Amazonia that displays a broadly conserved external morphology over its range. This wide geographical distribution and conservation of body form contrasts with the expected poor dispersal ability of the species, the tumultuous past of Amazonia, and the previously documented prevalence of cryptic species in widespread terrestrial organisms in this region. Here we investigate this homogeneity by examining hemipenial morphology and conducting phylogenetic analyses of mitochondrial (CYTB) and nuclear (C‐MOS) DNA sequence data from 49 individuals sampled across Amazonia. We detected remarkable variation in hemipenial morphology within this species, with multiple cases of sympatric occurrence of distinct hemipenial morphotypes. Phylogenetic analyses revealed highly divergent lineages corroborating the patterns suggested by the hemipenial morphotypes, including co‐occurrence of different lineages. The degrees of genetic and morphological distinctness, as well as instances of sympatry among mtDNA lineages/morphotypes without nuDNA allele sharing, suggest that I. elegans is a complex of cryptic species. An extensive and integrative taxonomic revision of the I. elegans complex throughout its wide geographical range is needed. © 2012 The Linnean Society of London, Zoological Journal of the Linnean Society, 2012, 166 , 361–376.     

Phylogeography and the origins of range disjunctions in a north temperate fish,the pygmy whitefish (Prosopium coulterii), inferred from mitochondrial and nuclear DNA sequence analysis

Aim To investigate the degree of phylogeographical divergence within pygmy whitefish (Prosopium coulterii) and to test hypotheses concerning the origin of disjunct populations within North America. Location North America from western Alaska to Lake Superior. Methods Mitochondrial (ATPase subunit VI) and nuclear (ITS‐1, ITS‐2) DNA sequence variation was assessed across the species’ North American range to test for the existence of distinct phylogeographical groupings of pygmy whitefish associated with known glacial refugia. Coalescent simulations of the mitochondrial DNA (mtDNA) data were used to test hypotheses of population structure. Results This species is composed of two monophyletic mitochondrial clades across its North American range. The two mtDNA clades differed by an average 3.3% nucleotide sequence divergence. These clades were also distinguished by ITS‐2, but the relationships among lineages were not resolved by the ITS‐1 analysis. Coalescent analyses rejected the null hypothesis of the current disjunct distributions being a result of fragmentation of a single widespread ancestral lineage across a variety of effective population sizes and divergence times. Main conclusions The current range disjunctions of pygmy whitefish in North America probably resulted from isolation, genetic divergence, and selective dispersal from at least two major Pleistocene glacial refugia: Beringia and Cascadia. More recent isolation and dispersal from an upper Mississippi refugium is suggested by relationships within one of the clades and by distributional evidence from co‐distributed species. The Beringian and Cascadian refugia have played major roles in the zoogeography of Nearctic temperate aquatics, but the roles of smaller refugia appear more variable among other species.     

Untangling the evolutionary history of a highly polymorphic species: introgressive hybridization and high genetic structure in the desert cichlid fish Herichtys minckleyi

Understanding the origin of biodiversity requires knowledge on the evolutionary processes that drive divergence and speciation, as well as on the processes constraining it. Intraspecific polymorphisms can provide insight into the mechanisms that generate and maintain phenotypic, behavioural and life history diversification, and can help us understand not only the processes that lead to speciation but also the processes that prevent local fixation of morphs. The ‘desert cichlid’ Herichtys minckleyi is a highly polymorphic species endemic to a biodiversity hotspot in northern Mexico, the Cuatro Ciénegas valley. This species is polymorphic in body shape and trophic apparatus, and eco‐morphotypes coexist in small spring‐fed lagoons across the valley. We investigated the genetic structure of these polymorphisms and their phylogeographic history by analysing the entire control region of the mitochondrial DNA and 10 nuclear microsatellite markers in several populations from different sites and morphs. We found two very divergent mitochondrial lineages that most likely predate the closing of the valley and are not associated with morphotypes or sites. One of these lineages is also found in the sister species Herichthys cyanoguttatus. Data from neutral microsatellite markers suggest that most lagoons or drainages constitute their own genetic cluster with sympatric eco‐morphotypes forming panmictic populations. Alternative mechanisms such as phenotypic plasticity and a few loci controlled traits provide possible explanations for the sympatric coexistence of discrete nonoverlapping eco‐morphotypes with apparent lack of barriers to gene flow within multiple lagoons and drainages.     

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

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

Unravelling evolutionary lineages among South African velvet worms (Onychophora: Peripatopsis) provides evidence for widespread cryptic speciation

The endemic South African velvet worm genus Peripatopsis currently contains eight recognized species described from variable morphological characters and the current taxonomy is unsatisfactory. In an attempt to investigate evolutionary relationships within Peripatopsis, we collected 137 individuals from 34 sample localities for six of the eight species. Sequence data derived from two partial mitochondrial (mt)DNA gene loci (COI and 12S rRNA), as well as partial sequence data from the ribosomal nuclear 18S rDNA locus in combination with gross morphological characters and scanning electron microscopy (SEM), was used to examine evolutionary relationships. Phylogenetic relationships were investigated using minimum evolution (ME) and Bayesian inferences (BI). Additionally, we also undertook a maximum likelihood (ML) analyses on the combined DNA sequence data set. The combined DNA evidence topologies derived from the ME, BI, and ML was highly congruent and was characterized by the presence of multiple lineages within recognized taxa. Peripatopsis clavigera, Peripatopsis moseleyi, and Peripatopsis sedgwicki each comprised two evolutionary lineages; Peripatopsis capensis comprised three; and Peripatopsis balfouri comprised six operational taxonomic units respectively. Genealogical exclusivity at both mtDNA and nuclear DNA among the geographically coherent groups coupled with pronounced sequence divergence suggested a two‐fold increase in the number of species within Peripatopsis. Previously used gross morphological characters (such as the number of leg pairs and colour) were either highly variable within operational taxonomic units, or were invariant, suggesting that alternative morphological characters are necessary for species discrimination. SEM results revealed potentially useful diagnostic characters that can discriminate between at least discriminate some of the newly‐identified lineages. © 2009 The Linnean Society of London, Biological Journal of the Linnean Society 2009, 97 , 200–216.     

Clustering of fungal community internal transcribed spacer sequence data obscures taxonomic diversity

Next‐generation DNA sequencing has enabled a rapid expansion in the size of molecular fungal ecology studies employing the nuclear internal transcribed spacer (ITS) region. Many sequence‐processing pipelines and protocols require sequence clustering to generate operational taxonomic units (OTUs) based on sequence similarity as a step to reduce total data quantity and complexity prior to taxonomic assignment. However, the consequences of ITS sequence clustering in regard to sample taxonomic coverage have not been carefully examined. Here we demonstrate that typically used clustering thresholds for fungal ITS sequences result in statistically significant losses in taxonomic coverage. Analyses using environmentally derived fungal sequences indicated an average of 3.1% of species went undetected (P < 0.05) if the sequences were denoised and clustered at a 97% threshold prior to taxonomic assignment. Additionally, an in silico analysis using a reference fungal ITS database suggested that approximately 25% of species went undetected if the sequences were clustered prior to taxonomic assignment. Finally, analysis of sequences derived from pure‐cultured fungal isolates of known identity indicated sequence denoising and clustering were not critical in improving identification accuracy.     

Divergence, hybridization, and recombination in the mitochondrial genome of the human pathogenic yeast Cryptococcus gattii

The inheritance of mitochondrial genes and genomes are uniparental in most sexual eukaryotes. This pattern of inheritance makes mitochondrial genomes in natural populations effectively clonal. Here, we examined the mitochondrial population genetics of the emerging human pathogenic fungus Cryptococcus gattii . The DNA sequences for five mitochondrial DNA fragments were obtained from each of 50 isolates belonging to two evolutionary divergent lineages, VGI and VGII. Our analyses revealed a greater sequence diversity within VGI than that within VGII, consistent with observations of the nuclear genes. The combined analyses of all five gene fragments indicated significant divergence between VGI and VGII. However, the five individual genealogies showed different relationships among the isolates, consistent with recent hybridization and mitochondrial gene transfer between the two lineages. Population genetic analyses of the multilocus data identified evidence for predominantly clonal mitochondrial population structures within both lineages. Interestingly, there were clear signatures of recombination among mitochondrial genes within the VGII lineage. Our analyses suggest historical mitochondrial genome divergence within C. gattii , but there is evidence for recent hybridization and recombination in the mitochondrial genome of this important human yeast pathogen.     

HOST‐SPECIALIST LINEAGES DOMINATE THE ADAPTIVE RADIATION OF REEF CORAL ENDOSYMBIONTS

Bursts in species diversification are well documented among animals and plants, yet few studies have assessed recent adaptive radiations of eukaryotic microbes. Consequently, we examined the radiation of the most ecologically dominant group of endosymbiotic dinoflagellates found in reef‐building corals, Symbiodinium Clade C, using nuclear ribosomal (ITS2), chloroplast (psbA^ncr), and multilocus microsatellite genotyping. Through a hierarchical analysis of high‐resolution genetic data, we assessed whether ecologically distinct Symbiodinium, differentiated by seemingly equivocal rDNA sequence differences, are independent species lineages. We also considered the role of host specificity in Symbiodinium speciation and the correspondence between endosymbiont diversification and Caribbean paleo‐history. According to phylogenetic, biological, and ecological species concepts, Symbiodinium Clade C comprises many distinct species. Although regional factors contributed to population‐genetic structuring of these lineages, Symbiodinium diversification was mainly driven by host specialization. By combining patterns of the endosymbiont's host specificity, water depth distribution, and phylogeography with paleo‐historical signals of climate change, we inferred that present‐day species diversity on Atlantic coral reefs stemmed mostly from a post‐Miocene adaptive radiation. Host‐generalist progenitors spread, specialized, and diversified during the ensuing epochs of prolonged global cooling and change in reef‐faunal assemblages. Our evolutionary reconstruction thus suggests that Symbiodinium undergoes “boom and bust” phases in diversification and extinction during major climate shifts.     

Two deep evolutionary lineages in the circumtropical glasseye Heteropriacanthus cruentatus (Teleostei,Priacanthidae) with admixture in the south‐western Indian Ocean

A phylogeographic study of the circumtropical glasseye Heteropriacanthus cruentatus was conducted. Molecular analyses indicate two mitochondrial cytochrome c oxidase subunit I (coI) lineages that are 10·4% divergent: one in the western Atlantic (Caribbean) and another that was detected across the Indo‐Pacific. A fixed single nucleotide polymorphism (SNP) was detected at a nuclear locus (S7 ribosomal protein) and is consistent with this finding. There is evidence of recent dispersal from the Atlantic to the Indian Ocean with individuals of mixed lineages detected in South Africa and the Mozambique Channel. Using coalescent analyses of the mitochondrial dataset, time of divergence between lineages was estimated to be c. 15·3 million years. The deep divergence between these two lineages indicates distinct evolutionary units, however, due to the lack of morphological differences and evidence of hybridization between lineages, taxonomic revision is not suggested at this time.     

Evolutionary history of tall fescue morphotypes inferred from molecular phylogenetics of the <Emphasis Type="Italic">Lolium</Emphasis>-<Emphasis Type="Italic">Festuca</Emphasis>species complex

Background  

The agriculturally important pasture grass tall fescue (Festuca arundinacea Schreb. syn. Lolium arundinaceum (Schreb.) Darbysh.) is an outbreeding allohexaploid, that may be more accurately described as a species complex consisting of three major (Continental, Mediterranean and rhizomatous) morphotypes. Observation of hybrid infertility in some crossing combinations between morphotypes suggests the possibility of independent origins from different diploid progenitors. This study aims to clarify the evolutionary relationships between each tall fescue morphotype through phylogenetic analysis using two low-copy nuclear genes (encoding plastid acetyl-CoA carboxylase [Acc1] and centroradialis [CEN]), the nuclear ribosomal DNA internal transcribed spacer (rDNA ITS) and the chloroplast DNA (cpDNA) genome-located matK gene. Other taxa within the closely related Lolium-Festuca species complex were also included in the study, to increase understanding of evolutionary processes in a taxonomic group characterised by multiple inter-specific hybridisation events.     

Cytonuclear discordance and the species status of Myotis myotis and Myotis blythii (Chiroptera)

We used both highly variable mitochondrial and nuclear loci to investigate the large mouse‐eared bat species complex in the Western Palaearctic to clarify their systematic position. Although mitochondrial lineages show no species segregation and some haplotypes are shared between Myotis myotis and Myotis blythii sensu lato, Bayesian clustering methods based on multilocus genotypes indicate highly concordant nuclear and morphological species assignment. These multilocus, nuclear analyses detected only a single putative F₁ hybrid in the extensive areas of sympatry sampled, thus confirming the biological species status of M. myotis and M. blythii s.l. We propose that the strong cytonuclear discordance in these species complex results from a combination of prior spatial isolation of the two species in different glacial refugia, followed by a succession of mitochondrial introgression events that occurred during the eastward and westward expansions of M. myotis and of M. blythii, respectively. The nuclear markers further indicate the presence of a notable genetic discontinuity within M. myotis that broadly separates populations into an eastern and a western component with an overlap zone in the Balkans. This eastern and western discontinuity is also apparent in the mitochondrial lineages with the D haplogroup largely confined to samples found in Thrace and Asia Minor. None of these genetic discontinuities correspond to the distribution of the two commonly recognized M. myotis subspecies (myotis and macrocephalicus). We also show that distinct morphological subspecies within M. blythii (oxygnathus, omari, risorius and lesviacus) in Europe and the near‐East do not correlate with significant evolutionarily units, whether identified by mitochondrial or nuclear data and thus only represent local morphological variants with little taxonomic relevance.     

Molecular characterization of isolates of Beauveria bassiana obtained from overwintering and summer populations of Sunn Pest (Eurygaster integriceps)

Aim:  To examine whether isolates of the entomopathogenic fungus Beauveria bassiana are more closely associated to their summer hosts compared with overwintering hosts, with recently developed molecular tools based on mitochondrial regions. Methods and Results:  Primers for the traditional ITS1‐5·8S‐ITS2 region and two mitochondrial intergenic regions, namely, nad3‐atp9 and atp6‐rns, were used. All amplified products were sequenced, aligned and Neighbour‐Joining (NJ), parsimony and Bayesian phylogenetic inference analyses were performed. The isolates examined were grouped with very good support into three distinct groups, two of them showed geographical correlation, but no clear association to their host. Conclusions:  The mitochondrial intergenic regions used were more informative than the nuclear ITS1‐5·8S‐ITS2 sequences. The sequence variability observed, that allowed the phylogenetic placement of the isolates into distinct groups, depended on the geographical origin of the isolates and can be exploited for designing group‐specific and isolate‐specific primers for their genetic fingerprinting. No clear associations with summer Sunn Pest populations were observed. Significance and Impact of the Study:  Studies on the genetic variability of biocontrol agents like B. bassiana are indispensable for the development of molecular tools for their future monitoring.     

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.     

Multiple glacial refugia and contemporary dispersal shape the genetic structure of an endemic amphibian from the Pyrenees

Historical factors (colonization scenarios, demographic oscillations) and contemporary processes (population connectivity, current population size) largely contribute to shaping species’ present‐day genetic diversity and structure. In this study, we use a combination of mitochondrial and nuclear DNA markers to understand the role of Quaternary climatic oscillations and present‐day gene flow dynamics in determining the genetic diversity and structure of the newt Calotriton asper (Al. Dugès, 1852), endemic to the Pyrenees. Mitochondrial DNA did not show a clear phylogeographic pattern and presented low levels of variation. In contrast, microsatellites revealed five major genetic lineages with admixture patterns at their boundaries. Approximate Bayesian computation analyses and linear models indicated that the five lineages likely underwent separate evolutionary histories and can be tracked back to distinct glacial refugia. Lineage differentiation started around the Last Glacial Maximum at three focal areas (western, central and eastern Pyrenees) and extended through the end of the Last Glacial Period in the central Pyrenees, where it led to the formation of two more lineages. Our data revealed no evidence of recent dispersal between lineages, whereas borders likely represent zones of secondary contact following expansion from multiple refugia. Finally, we did not find genetic evidence of sex‐biased dispersal. This work highlights the importance of integrating past evolutionary processes and present‐day gene flow and dispersal dynamics, together with multilocus approaches, to gain insights into what shaped the current genetic attributes of amphibians living in montane habitats.