One, two or more species? Mitonuclear discordance and species delimitation

Congruent patterns of cryptic cladogenesis revealed using RADseq and Sanger sequencing in a velvet worm species complex (Onychophora: Peripatopsidae: Peripatopsis sedgwicki)

In the present study, first generation DNA sequencing (mitochondrial cytochrome c oxidase subunit one, COI) and reduced-representative genomic RADseq data were used to understand the patterns and processes of diversification of the velvet worm, Peripatopsis sedgwicki species complex across its distribution range in South Africa. For the RADseq data, three datasets (two primary and one supplementary) were generated corresponding to 1,259-11,468 SNPs, in order to assess the diversity and phylogeography of the species complex. Tree topologies for the two primary datasets were inferred using maximum likelihood and Bayesian inferences methods. Phylogenetic analyses using the COI datasets retrieved four distinct, well-supported clades within the species complex. Five species delimitation methods applied to the COI data (ASAP, bPTP, bGMYC, STACEY and iBPP) all showed support for the distinction of the Fort Fordyce Nature Reserve specimens. In the main P. sedgwicki species complex, the species delimitation methods revealed a variable number of operational taxonomic units and overestimated the number of putative taxa. Divergence time estimates coupled with the geographic exclusivity of species and phylogeographic results suggest recent cladogenesis during the Plio/Pleistocene. The RADseq data were subjected to a principal components analysis and a discriminant analysis of principal components, under a maximum-likelihood framework. The latter results corroborate the four main clades observed using the COI data, however, applying additional filtering revealed additional diversity. The high overall congruence observed between the RADseq data and COI data suggest that first generation sequence data remain a cheap and effective method for evolutionary studies, although RADseq does provide a far greater resolution of contemporary temporo-spatial patterns.

Population differentiation and dynamics of five pioneer species of Gaultheria from the secondary forests in subtropical China

BACKGROUND

The influence of native secondary succession associated with anthropogenic disturbance on the biodiversity of the forests in subtropical China remains uncertain. In particular, the evolutionary response of small understory shrubs, particularly pioneer species inhabiting continuously disturbed habitats, to topographic heterogeneity and climate change is poorly understood. This study aimed to address this knowledge gap by focusing on the Gaultheria crenulata group, a clade of small pioneer shrubs in subtropical China.

RESULTS

We examined the genetic structure and demographic history of all five species of the G. crenulata group with two maternally inherited chloroplast DNA (cpDNA) fragments and two biparentally inherited low-copy nuclear genes (LCG) over 89 natural populations. We found that the genetic differentiation of this group was influenced by the geomorphological boundary between different regions of China in association with Quaternary climatic events. Despite low overall genetic diversity, we observed an isolation-by-distance (IBD) pattern at a regional scale, rather than isolation-by-environment (IBE), which was attributed to ongoing human disturbance in the region.

CONCLUSION

Our findings suggest that the genetic structure of the G. crenulata group reflects the interplay of geological topography, historical climates, and anthropogenic disturbance during the Pliocene-Pleistocene-Holocene periods in subtropical China. The observed IBD pattern, particularly prominent in western China, highlights the role of limited dispersal and gene flow, possibly influenced by physical barriers or decreased connectivity over geographic distance. Furthermore, the east-to-west trend of gene flow, potentially facilitated by the East Asian monsoon system, underscores the complex interplay of biotic and abiotic factors shaping the genetic dynamics of pioneer species in subtropical China's secondary forests. These findings can be used to assess the impact of environmental changes on the adaptation and persistence of biodiversity in subtropical forest ecosystems.

The efficacy of single mitochondrial genes at reconciling the complete mitogenome phylogeny-a case study on dwarf chameleons

Although genome-scale data generation is becoming more tractable for phylogenetics, there are large quantities of single gene fragment data in public repositories and such data are still being generated. We therefore investigated whether single mitochondrial genes are suitable proxies for phylogenetic reconstruction as compared to the application of full mitogenomes. With near complete taxon sampling for the southern African dwarf chameleons (Bradypodion), we estimated and compared phylogenies for the complete mitogenome with topologies generated from individual mitochondrial genes and various combinations of these genes. Our results show that the topologies produced by single genes (ND2, ND4, ND5, COI, and COIII) were analogous to the complete mitogenome, suggesting that these genes may be reliable markers for generating mitochondrial phylogenies in lieu of generating entire mitogenomes. In contrast, the short fragment of 16S commonly used in herpetological systematics, produced a topology quite dissimilar to the complete mitogenome and its concatenation with ND2 weakened the resolution of ND2. We therefore recommend the avoidance of this 16S fragment in future phylogenetic work.

Species delimitation, discovery and conservation in a tiger beetle species complex despite discordant genetic data

In an age of species declines, delineating and discovering biodiversity is critical for both taxonomic accuracy and conservation. In recent years, there has been a movement away from using exclusively morphological characters to delineate and describe taxa and an increase in the use of molecular markers to describe diversity or through integrative taxonomy, which employs traditional morphological characters, as well as genetic or other data. Tiger beetles are charismatic, of conservation concern, and much work has been done on the morphological delineation of species and subspecies, but few of these taxa have been tested with genetic analyses. In this study, we tested morphologically based taxonomic hypotheses of polymorphic tiger beetles in the Eunota circumpicta (LaFerté-Sénectère, 1841) species complex using multilocus genomic and mtDNA analyses. We find multiple cryptic species within the previous taxonomic concept of Eunota circumpicta, some of which were historically recognized as subspecies. We found that the mtDNA and genomic datasets did not identify the same taxonomic units and that the mtDNA was most at odds with all other genetic and morphological patterns. Overall, we describe new cryptic diversity, which raises important conservation concerns, and provide a working example for testing species and subspecies validity despite discordant data.

The effect of missing data on evolutionary analysis of sequence capture bycatch, with application to an agricultural pest

Sequence capture is a genomic technique that selectively enriches target sequences before high throughput next-generation sequencing, to generate specific sequences of interest. Off-target or 'bycatch' data are often discarded from capture experiments, but can be leveraged to address evolutionary questions under some circumstances. Here, we investigated the effects of missing data on a variety of evolutionary analyses using bycatch from an exon capture experiment on the global pest moth, Helicoverpa armigera. We added > 200 new samples from across Australia in the form of mitogenomes obtained as bycatch from targeted sequence capture, and combined these into an additional larger dataset to total > 1000 mitochondrial cytochrome c oxidase subunit I (COI) sequences across the species' global distribution. Using discriminant analysis of principal components and Bayesian coalescent analyses, we showed that mitogenomes assembled from bycatch with up to 75% missing data were able to return evolutionary inferences consistent with higher coverage datasets and the broader literature surrounding H. armigera. For example, low-coverage sequences broadly supported the delineation of two H. armigera subspecies and also provided new insights into the potential for geographic turnover among these subspecies. However, we also identified key effects of dataset coverage and composition on our results. Thus, low-coverage bycatch data can offer valuable information for population genetic and phylodynamic analyses, but caution is required to ensure the reduced information does not introduce confounding factors, such as sampling biases, that drive inference. We encourage more researchers to consider maximizing the potential of the targeted sequence approach by examining evolutionary questions with their off-target bycatch where possible-especially in cases where no previous mitochondrial data exists-but recommend stratifying data at different genome coverage thresholds to separate sampling effects from genuine genomic signals, and to understand their implications for evolutionary research.

Using genome-wide data to ascertain taxonomic status and assess population genetic structure for Houston toads (Bufo [= Anaxyrus] houstonensis)

The Houston toad (Bufo [= Anaxyrus] houstonensis) is an endangered amphibian with a small geographic range. Land-use changes have primarily driven decline in B. houstonensis with population supplementation predominant among efforts to reduce its current extinction risk. However, there has been historic uncertainty regarding the evolutionary and conservation significance of B. houstonensis. To this end, we used 1170 genome-wide nuclear DNA markers to examine phylogenetic relationships between our focal taxon, representatives of the Nearctic B. americanus group, and B. nebulifer, a sympatric Middle American species. Phylogenetic analyses indicate B. houstonensis is a taxon that is distinct from B. americanus. We corroborated such genetic distinctiveness with an admixture analysis that provided support for recent reproductive isolation between B. americanus and B. houstonensis. However, ABBA-BABA tests for ancient admixture indicated historic gene flow between Nearctic species while no signal of historic gene flow was detected between Nearctic and Middle-American species. We used an admixture analysis to recognize four Management Units (MU) based on observed genetic differentiation within B. houstonensis and recommend captive propagation, population supplementation, and habitat restoration efforts specific to each MU. Our results re-affirm the evolutionary novelty of an endangered relict.

Molecular insights versus morphological traits: rethinking identification of the closely related Angiostrongylus cantonensis and Angiostrongylus malaysiensis

BACKGROUND

The closely related Angiostrongylus cantonensis and Angiostrongylus malaysiensis have been reported to coexist in Thailand and share similar hosts and life cycles. Recently, in an angiostrongyliasis outbreak in Thailand, both A. cantonensis and A. malaysiensis were found in the cerebrospinal fluid of affected patients. Morphological similarities, overlapping distribution, shared hosts and habitats, and the close genetics of the two Angiostrongylus species can complicate accurate species identification. Addressing these challenges, this study aims to evaluate whether a correlation between the morphological and genetic identities of A. cantonensis and A. malaysiensis can improve species identification accuracy.

METHODS

Angiostrongylus spp. specimens from five zoogeographical regions in Thailand were subjected to morphological and molecular identification using the mitochondrial cytochrome b gene and the nuclear internal transcribed spacer 2 region (ITS2). The morphological characters for males and females were then validated using the species identity obtained from the nuclear ITS2 region.

RESULTS

The results revealed that morphological misidentifications between these two closely related species are common due to overlapping morphological characters. Although certain male traits such as body length and width aided species differentiation, female traits were found to be less reliable. Furthermore, hybrid forms (8.2%) were revealed through the ITS2 results, which can further complicate morphological identification. Mito-nuclear discordance was also present in 1.9% of the Angiostrongylus specimens from Thailand, suggesting a complex historical interbreeding between the species.

CONCLUSIONS

Based on our findings, we suggest that nuclear ITS2 is a reliable marker for species identification of A. cantonensis and A. malaysiensis, especially in regions where both species coexist. Additionally, the scope and consequences of hybridization between the two closely related Angiostrongylus species should be further investigated in Thailand.

Cryptic diversity begets challenges and opportunities in biodiversity research

How many species of life are there on Earth? This is a question that we want to know but cannot yet answer. Some scholars speculate that the number of species may reach 2.2 billion when considering cryptic diversity and that each morphology-based insect species may contain an average of 3.1 cryptic species. With nearly two million described species, such high estimates of cryptic diversity would suggest that cryptic species are widespread. The development of molecular species delimitation has led to the discovery of a large number of cryptic species, and cryptic biodiversity has gradually entered our field of vision and attracted more attention. This paper introduces the concept of cryptic species, how they evolve, and methods by which they may be discovered and confirmed, and provides theoretical and methodological guidance for the study of hidden species. A workflow of how to confirm cryptic species is provided. In addition, the importance and reliability of multi-evidence-based integrated taxonomy are reaffirmed as a way to better standardize decision-making processes. Special focus on cryptic diversity and increased funding for taxonomy is needed to ensure that cryptic species in hyperdiverse groups are discoverable and described. An increased focus on cryptic species in the future will naturally arise as more difficult groups are studied, and thereby, we may finally better understand the rules governing the evolution and maintenance of cryptic biodiversity.

Special Issue "Primate Phylogeny and Genetics"

New phylogenetic tools and population genetics methods have been developed and vastly advanced over the last decade [...].

Morphological and molecular phylogenetic characterization of Sarcocystis kani sp. nov. and other novel, closely related Sarcocystis spp. infecting small mammals and colubrid snakes in Asia

We investigated the morphology and phylogenetic relationships of novel and previously recognized Sarcocystis spp. infecting small mammals and colubrid snakes in Asia. The nuclear 18S rRNA and mitochondrial cox1 of Sarcocystis sp.1 from mangrove snakes (Boiga dendrophila) in Thailand and Sarcocystis sp.2 from a ricefield rat (Rattus argentiventer) in Sumatra were partially sequenced. Sporocysts of Sarcocystis sp.1 induced development of sarcocysts in experimentally infected rats, which showed a unique ultrastructure that was observed previously by S.P. Kan in rats from Malaysia; therefore, we describe this species as Sarcocystis kani sp. nov. Its integration into the 18S rRNA phylogeny of Sarcocystis spp. cycling between small mammals and colubrid snakes helped clarify relationships among the so-called S. zuoi-complex of molecularly cryptic species: Sarcocystis kani sp. nov., S. sp.2, S. attenuati, S. scandentiborneensis, and S. zuoi were all included in this clade. Tree topology was resolved into dichotomies congruent with the morphological disparities between the taxa. However, cox1 gene sequencing (including newly sequenced S. singaporensis and S. zamani) revealed that Sarcocystis kani, S. attenuati, and S. scandentiborneensis were identical suggesting a recent, common ancestry. To identify other distinctive features, lineage-specific molecular patterns within both genes were examined revealing that all 18S rRNA sequences of the S. zuoi - complex possess a unique, 7-nt long motif in helix 38 of domain V7 that was different in S. clethrionomyelaphis which branched off basally from the complex. Three-dimensional homology modelling of COX1 protein structure identified amino acid substitutions within the barcode area specific for the S. zuoi-complex and substantial divergence in structurally important amino acids between Sarcocystis species of snakes as definitive hosts and other lineages of the Sarcocystidae. We discuss the utility of selected genes for species delimitation of the Sarcocystis spp. under investigation, which probably evolved during recent radiations of their intermediate and definitive hosts.

Mitonuclear discordance results from incomplete lineage sorting, with no detectable evidence for gene flow, in a rapid radiation of Todiramphus kingfishers

Many organisms possess multiple discrete genomes (i.e. nuclear and organellar), which are inherited separately and may have unique and even conflicting evolutionary histories. Phylogenetic reconstructions from these discrete genomes can yield different patterns of relatedness, a phenomenon known as cytonuclear discordance. In many animals, mitonuclear discordance (i.e. discordant evolutionary histories between the nuclear and mitochondrial genomes) has been widely documented, but its causes are often considered idiosyncratic and inscrutable. We show that a case of mitonuclear discordance in Todiramphus kingfishers can be explained by extensive genome-wide incomplete lineage sorting (ILS), likely a result of the explosive diversification history of this genus. For these kingfishers, quartet frequencies reveal that the nuclear genome is dominated by discordant topologies, with none of the internal branches in our consensus nuclear tree recovered in >50% of genome-wide gene trees. Meanwhile, a lack of inter-species shared ancestry, non-significant pairwise tests for gene flow, and little evidence for meaningful migration edges between species, leads to the conclusion that gene flow cannot explain the mitonuclear discordance we observe. This lack of evidence for gene flow combined with evidence for extensive genome-wide gene tree discordance, a hallmark of ILS, leads us to conclude that the mitonuclear discordance we observe likely results from ILS, specifically deep coalescence of the mitochondrial genome. Based on this case study, we hypothesize that similar demographic histories in other 'great speciator' taxa across the Indo-Pacific likely predispose these groups to high levels of ILS and high likelihoods of mitonuclear discordance.

Evaluating the taxonomic status of the large sized Tricolia Risso, 1826 in the Northeast Atlantic and Mediterranean Sea

Despite a long history of taxonomic studies on the genus Tricolia Risso, 1826, there is a shortfall on thorough systematic molecular reviews of the taxon from the NE Atlantic and Mediterranean coasts. Aiming to assess the genetic distinctness among morphospecies and the taxonomic status of currently accepted large sized species in these areas, we conducted a molecular phylogenetic analysis of the genus based on one mitochondrial (cox1) and two nuclear (28S and ITS2) markers. Seven Tricolia species were consistently retrieved in the analyses, including a new genetic lineage in the NE Atlantic designated as Tricolia sp. 1. Molecular analyses revealed that only one species, T. azorica, occurs in the NE Atlantic archipelagos. The sister taxa T. pullus (Mediterranean) and T. picta (NE Atlantic) should be classified as distinct species, instead of subspecies of the T. pullus group (sensu Gofas 1982). Tricolia miniata is also a complex of species in the Mediterranean and future studies across the distribution range are necessary to clarify its status.

Population genomics indicate three different modes of divergence and speciation with gene flow in the green-winged teal duck complex

The processes leading to divergence and speciation can differ broadly among taxa with different life histories. We examine these processes in a small clade of ducks with historically uncertain relationships and species limits. The green-winged teal (Anas crecca) complex is a Holarctic species of dabbling duck currently categorized as three subspecies (Anas crecca crecca, A. c. nimia, and A. c. carolinensis) with a close relative, the yellow-billed teal (Anas flavirostris) from South America. A. c. crecca and A. c. carolinensis are seasonal migrants, while the other taxa are sedentary. We examined divergence and speciation patterns in this group, determining their phylogenetic relationships and the presence and levels of gene flow among lineages using both mitochondrial and genome-wide nuclear DNA obtained from 1,393 ultraconserved element (UCE) loci. Phylogenetic relationships using nuclear DNA among these taxa showed A. c. crecca, A. c. nimia, and A. c. carolinensis clustering together to form one polytomous clade, with A. flavirostris sister to this clade. This relationship can be summarized as (crecca, nimia, carolinensis)(flavirostris). However, whole mitogenomes revealed a different phylogeny: (crecca, nimia)(carolinensis, flavirostris). The best demographic model for key pairwise comparisons supported divergence with gene flow as the probable speciation mechanism in all three contrasts (crecca-nimia, crecca-carolinensis, and carolinensis-flavirostris). Given prior work, gene flow was expected among the Holarctic taxa, but gene flow between North American carolinensis and South American flavirostris (M ∼0.1-0.4 individuals/generation), albeit low, was not expected. Three geographically oriented modes of divergence are likely involved in the diversification of this complex: heteropatric (crecca-nimia), parapatric (crecca-carolinensis), and (mostly) allopatric (carolinensis-flavirostris). Our study shows that ultraconserved elements are a powerful tool for simultaneously studying systematics and population genomics in systems with historically uncertain relationships and species limits.

Unrecognized diversity, genetic structuring, and phylogeography of the genus Triplophysa (Cypriniformes: Nemacheilidae) sheds light on two opposite colonization routes during Quaternary glaciation that occurred in the Qilian Mountains

In recent years, the influence of historical geological and climatic events on the evolution of flora and fauna in the Tibetan Plateau has been a hot research topic. The Qilian Mountain region is one of the most important sources of biodiversity on the Qinghai-Tibet Plateau. Many species existed in the region during the Pleistocene glacial oscillation, and the complex geographical environment provided suitable conditions for the survival of local species. The shrinkage, expansion, and transfer of the distribution range and population size of species have significant effects on genetic diversity and intraspecific differentiation. To reveal the effects of geological uplift and climate oscillation on the evolution of fish populations in the Qilian Mountains, we investigated the genetic structure, phylogenetic relationship, and phylogeographical characteristics of genus Triplophysa species in the Qilian Mountains using the mitochondrial DNA gene (COI), three nuclear genes (RAG1, sRH, and Myh6) and 11 pairs of nuclear microsatellite markers. We collected 11 species of genus Triplophysa living in the Qilian Mountains, among which Triplophysa hsutschouensis and Triplophysa papillosolabiata are widely distributed in the rivers on the northern slope of the Qilian Mountains. There was a high degree of lineage differentiation among species, and the genetic diversity of endemic species was low. The different geographical groups of T. papillosolabiata presented some allogeneic adaptation and differentiation, which was closely related to the changes in the river system. Except for the population expansion event of T. hsutschouensis during the last glacial period of the Qinghai-Tibet Plateau (0.025 MYA), the population sizes of other plateau loach species remained stable without significant population expansion. Starting from the east and west sides of the Qilian Mountains, T. hsutschouensis, and T. papillosolabiata showed two species colonization routes in opposite directions. The geological events of the uplift of the Qinghai-Tibet Plateau and the climatic oscillation of the Quaternary glaciation had a great influence on the genetic structure of the plateau loach in the Qilian Mountains, which promoted the genetic differentiation of the plateau loach and formed some unique new species. The results of this study have important guiding significance for fish habitat protection in the Qilian Mountains.

Cryptic diversity of Crocidura shrews in the savannahs of Eastern and Southern Africa

Crocidura (Eulipotyphla, Soricidae) is the most species-rich genus among mammals, with high cryptic diversity and complicated taxonomy. The hirta-flavescens group of Crocidura represents the most abundant and widespread shrews in savannahs of eastern and southern Africa, making them a suitable phylogeographical model for assessing the role of paleoclimatic changes on current biodiversity in open African habitats. We present the first comprehensive study on the phylogeography, evolutionary history, geographical distribution, systematics, and taxonomy of the group, using the integration of mitochondrial, genome-wide (ddRAD sequencing), morphological and morphometrical data collected from specimens over most of the known geographic distribution. Our genomic data confirmed the monophyly of this group and its sister relationship with the olivieri group of Crocidura. There is a substantial genetic variation within the hirta-flavescens group, with three highly supported clades showing parapatric distribution and which can be distinguished morphologically: C. hirta, distributed in both the Zambezian and Somali-Masai bioregions, C. flavescens, known from South Africa and south-western Zambia, and C. cf. flavescens, which is known to occur only in central and western Tanzania. Morphometric data revealed relatively minor differences between C. hirta and C. cf. flavescens, but they differ in the colouration of the pelage. Diversification of the hirta-flavescens group has most likely happened during phases of grassland expansion and contraction during Plio-Pleistocene climatic cycles. Eastern African Rift system, rivers, and the distinctiveness of Zambezian and Somali-Masai bioregions seem to have also shaped the pattern of their diversity, which is very similar to sympatric rodent species living in open habitats. Finally, we review the group's taxonomy and propose to revalidate C. bloyeti, currently a synonym of C. hirta, including the specimens treated as C. cf. flavescens.

New Color-Patterned Species of Microtendipes Kieffer, 1913 (Diptera: Chironomidae) and a Deep Intraspecific Divergence of Species by DNA Barcodes

The genus Microtendipes Kieffer (Diptera: Chironomidae) has a nearly worldwide distribution, comprising more than 60 species, which are further divided into two species groups based on larval stage. However, species delimitation and identification among the adults of this genus are controversial and uncertain. For instance, previous studies have provided many synonymies based on conspecific color pattern variations in Microtendipes species. Here, we used DNA barcode data to address Microtendipes species delimitation as well as to test whether color pattern variations can be diagnostic characters for interspecific identification. The 151 DNA barcodes used, 51 of which were contributed by our laboratory, represent 21 morphospecies. Species with specific color patterns could be accurately separated based on DNA barcodes. Consequently, the color patterns of adult males could be important diagnostic characters. The average intraspecific and interspecific sequence divergences were 2.8% and 12.5%, respectively, and several species exhibited deep intraspecific divergences higher than 5%. Molecular operational taxonomic units (OTUs) ranged from 21 to 73, based on methods including phylogenetic trees, the assemble species by automatic partitioning method, the Poisson tree process (PTP), and the general mixed Yule-coalescent (GMYC) method. As a result of these analyses, five new species were recognized (M. baishanzuensis sp. nov., M. bimaculatus sp. nov., M. nigrithorax sp. nov., M. robustus sp. nov., and M. wuyiensis sp. nov.).

A wide hybrid zone mediated by precipitation contributed to confused geographical structure of Scutiger boulengeri

Confused geographical structure of a population and mitonuclear discordance are shaped by a combination of rapid changes in population demographics and shifts in ecology. In this study, we generated a time-calibrated phylogeny of Scutiger boulengeri, an endemic Xizang alpine toad occurring in mountain streams on the Qinghai-Xizang (Tibet) Plateau (QTP). Based on three mitochondrial DNA (mtDNA) genes, eight clades were assigned to three deeply divergent lineages. Analysis of nuclear DNA (nuDNA) genes revealed three distinct clusters without geographic structure, indicating significantly high rates of gene flow. Coalescent theory framework analysis (approximate Bayesian computation model DIYABC and Migrate-N) suggested that divergence of the main intraspecific clusters was the result of hybridization after secondary contact in the Holocene around 0.59 million years ago (Ma). The ratio of mtDNA F _ST (fixation index) to nuDNA F _ST was 2.3, thus failing to show male-biased dispersal. Geographic cline analysis showed that a wide hybrid zone was initially established in southwestern China, without significant reproductive isolation but with strong introgression in S. boulengeri, suggesting high hybrid fitness. Furthermore, mtDNA genes exhibited isolation by distance (IBD) while nuDNA genes exhibited significant isolation by environment (IBE). Results suggested that mitonuclear discordance may have initially been caused by geographic isolation, followed by precipitation-mediated hybridization, producing a wide hybrid zone and geographic structure confusion of nuDNA genes in S. boulengeri. This study indicated that complicated historical processes may have led to specific genetic patterns, with a specific climate factor facilitating gene flow in the system.

Co-introduction of Dolicirroplectanum lacustre, a monogenean gill parasite of the invasive Nile perch Lates niloticus: intraspecific diversification and mitonuclear discordance in native versus introduced areas

The Nile perch (Lates niloticus) is a notorious invasive species. The introductions of Nile perch into several lakes and rivers in the Lake Victoria region led to the impoverishment of trophic food webs, particularly well documented in Lake Victoria. Additionally, its parasites were co-introduced, including Dolicirroplectanum lacustre (Monogenea, Diplectanidae). Dolicirroplectanum lacustre is the single monogenean gill parasite of latid fishes (Lates spp.) inhabiting several major African freshwater systems. We examined the intra-specific diversification of D. lacustre from Lates niloticus in Lake Albert, Uganda (native range) and Lake Victoria (introduced range) by assessing morphological and genetic differentiation, and microhabitat preference. We expected reduced morphological and genetic diversity for D. lacustre in Lake Victoria compared with Lake Albert, as a result of the historical introductions. We found that D. lacustre displayed high morphological variability within and between African freshwaters, with two morphotypes identified, as in former studies. The single shared morphotype between Lake Albert and Lake Victoria displayed similar levels of haplotype and nucleotide diversity between the lakes. Mitonuclear discordance within the morphotypes of D. lacustre indicates an incomplete reproductive barrier between the morphotypes. The diversification in the mitochondrial gene portion is directly linked with the morphotypes, while the nuclear gene portions indicate conspecificity. Based on our results, we reported reduced genetic and morphological diversity, potentially being a result of a founder effect in Lake Victoria.

Extremely divergent COI sequences within an amphipod species complex: A possible role for endosymbionts?

Some heritable endosymbionts can affect host mtDNA evolution in various ways. Amphipods host diverse endosymbionts, but whether their mtDNA has been influenced by these endosymbionts has yet to be considered. Here, we investigated the role of endosymbionts (microsporidians and Rickettsia) in explaining highly divergent COI sequences in Paracalliope fluviatilis species complex, the most common freshwater amphipods in New Zealand. We first contrasted phylogeographic patterns using COI, ITS, and 28S sequences. While molecular species delimitation methods based on 28S sequences supported 3-4 potential species (N, C, SA, and SB) among freshwater lineages, COI sequences supported 17-27 putative species reflecting high inter-population divergence. The deep divergence between NC and S lineages (~20%; 28S) and the substitution saturation on the 3rd codon position of COI detected even within one lineage (SA) indicate a very high level of morphological stasis. Interestingly, individuals infected and uninfected by Rickettsia comprised divergent COI lineages in one of four populations tested, suggesting a potential influence of endosymbionts in mtDNA patterns. We propose several plausible explanations for divergent COI lineages, although they would need further testing with multiple lines of evidence. Lastly, due to common morphological stasis and the presence of endosymbionts, phylogeographic patterns of amphipods based on mtDNA should be interpreted with caution.

Integrative species delimitation reveals fine-scale allopatric speciation in a good-flying insect: a case study on

Alpha taxonomy is fundamental for many biological fields. Delineation of the species boundary, however, can be challenging in a species complex, where different species share a similar morphology and diagnostic characters may not be available. In this context, integrative approaches that incorporate molecular and morphological data sets, and account for speciation history can be helpful to alpha taxonomy. Different approaches to species delimitation based on different assumptions are complementary and by integrating the results from multiple approaches we can generate a more reliable and objective taxonomic decision. In this study, we applied three molecular approaches to species delimitation and inferred the demographic history based on an isolation with migration model to test a morphologically based taxonomic hypothesis for the Cylindera pseudocylindriformis complex. We discuss the association between genetic divergence and microhabitat specialisation, and further corroborate that C. subtilis sp. nov. is a valid new species by integrating the results from model-based species delimitation and the genealogical divergence index. We argue that genetic endemism can occur at a small geographic scale, even in a winged insect like tiger beetles. Our results also indicated that there may still be undocumented species diversity of Taiwanese Cylindera remaining to be discovered. ZooBank LSID: urn:lsid:zoobank.org:pub:9DEC1432-365C-4872-8D06-73B95F30624F

Linking large-scale genetic structure of three Argynnini butterfly species to geography and environment

Understanding which factors and processes are associated with genetic differentiation within and among species remains a major goal in evolutionary biology. To explore differences and similarities in genetic structure and its association with geographical and climatic factors in sympatric sister species, we conducted a large‐scale (>32° latitude and >36° longitude) comparative phylogeographical study on three Argynnini butterfly species (Speyeria aglaja, Fabriciana adippe and F. niobe) that have similar life histories, but differ in ecological generalism and dispersal abilities. Analyses of nuclear (ddRAD‐sequencing derived SNP markers) and mitochondrial (COI sequences) data revealed differences between species in genetic structure and how genetic differentiation was associated with climatic factors (temperature, solar radiation, precipitation, wind speed). Geographical proximity accounted for much of the variation in nuclear and mitochondrial structure and evolutionary relationships in F. adippe and F. niobe, but only explained the pattern observed in the nuclear data in S. aglaja, for which mitonuclear discordance was documented. In all species, Iberian and Balkan individuals formed genetic clusters, suggesting isolation in glacial refugia and limited postglacial expansion. Solar radiation and precipitation were associated with the genetic structure on a regional scale in all species, but the specific combinations of environmental and geographical factors linked to variation within species were unique, pointing to species‐specific responses to common environments. Our findings show that the species share similar colonization histories, and that the same ecological factors, such as niche breadth and dispersal capacity, covary with genetic differentiation within these species to some extent, thereby highlighting the importance of comparative phylogeographical studies in sympatric sister species.

The dynamics and the timeline of speciation in the gall-forming aphid Geoica spp. within and among Pistacia host tree species

Trees of the genus Pistacia serve as obligate hosts for gall-forming aphids (Hemiptera, Aphididae, Fordini). Each aphid species induces a characteristic gall on a single Pistacia host species. The genus Geoica (Fordini) induce similar spherical closed galls on the lower side of the leaflet's midvein, on different Pistacia species. Two species of Pistacia trees that harbor Geoica galls grow naturally in Israel: P. palaestina and P. atlantica. We analyzed the phylogeny and the genetic structure of the Geoica species complex in Israel, and assessed the genetic differentiation and the level of host plant specificity of the aphids between P. atlantica and P. palaestina. We found that the splitting of the genus between P. atlantica and P. palaestina is estimated to have occurred 24-25 Ma (the Oligocene/Miocene boundary). Five different haplotypes suggesting five different species have been further speciating among Geoica spp., galling on P. atlantica, and an additional three species, on P. palaestina.

Mitonuclear discordance and patterns of reproductive isolation in a complex of simultaneously hermaphroditic species, the Allolobophora chlorotica case study

Historical events of population fragmentation, expansion and admixture over geological time may result in complex patterns of reproductive isolation and may explain why, for some taxa, the study of mitochondrial (mt) and nuclear (nu) genetic data results in discordant evolutionary patterns. Complex patterns of taxonomic diversity were recently revealed in earthworms for which distribution is largely the result of paleogeographical events. Here, we investigated reproductive isolation patterns in a complex of cryptic species of earthworms in which discordant patterns between mt and nu genetic lineages were previously revealed, the Allolobophora chlorotica aggregate. Using four nu microsatellite markers and a fragment of the cytochrome c oxidase subunit I mt gene, we carried out a parentage analysis to investigate the mating patterns (i) between individuals belonging to two divergent mt lineages that cannot be distinguished with nu markers and (ii) between individuals belonging to lineages that are differentiated both at the mt and nu levels. Amongst the 157 field-collected individuals, 66 adults were used in cross-breeding experiments to form 22 trios based on their assignment to a mt lineage, and 453 obtained juveniles were genotyped. We showed that adults that mated with both their potential mates in the trio produced significantly more juveniles. In crosses between lineages that diverged exclusively at the mt level, a sex-specific pattern of reproduction characteristic to each lineage was observed, suggesting a possible conflict of interest concerning the use of male/female function between mating partners. In crosses between lineages that diverged both at the mt and nu levels, a high production of cocoons was counterbalanced by a low hatching rate, suggesting a post-zygotic reproductive isolation. Different degrees of reproductive isolation, from differential sex allocation to post-zygotic isolation, were thus revealed. Lineages appear to be at different stages in the speciation process, which likely explain the observed opposite patterns of mitonuclear congruence.

Species delimitation and mitonuclear discordance within a species complex of biting midges

The inability to distinguish between species can be a serious problem in groups responsible for pathogen transmission. Culicoides biting midges transmit many pathogenic agents infecting wildlife and livestock. In North America, the C. variipennis species complex contains three currently recognized species, only one of which is a known vector, but limited species-specific characters have hindered vector surveillance. Here, genomic data were used to investigate population structure and genetic differentiation within this species complex. Single nucleotide polymorphism data were generated for 206 individuals originating from 17 locations throughout the United States and Canada. Clustering analyses suggest the occurrence of two additional cryptic species within this complex. All five species were significantly differentiated in both sympatry and allopatry. Evidence of hybridization was detected in three different species pairings indicating incomplete reproductive isolation. Additionally, COI sequences were used to identify the hybrid parentage of these individuals, which illuminated discordance between the divergence of the mitochondrial and nuclear datasets.

Evolutionary history of a desert perennial Arnebia szechenyi (Boraginaceae): Intraspecific divergence, regional expansion and asymmetric gene flow

The complex interactions of historical, geological and climatic events on plant evolution have been an important research focus for many years. However, the role of desert formation and expansion in shaping the genetic structures and demographic histories of plants occurring in arid areas has not been well explored. In the present study, we investigated the phylogeography of Arnebia szechenyi, a desert herb showing a near-circular distribution surrounding the Tengger Desert in Northwest China. We measured genetic diversity of populations using three maternally inherited chloroplast DNA (cpDNA) fragments and seven bi-paternally inherited nuclear DNA (nDNA) loci that were sequenced from individuals collected from 16 natural populations across its range and modelled current and historical potential habitats of the species. Our data indicated a considerably high level of genetic variation within A. szechenyi and noteworthy asymmetry in historical migration from the east to the west. Moreover, two nuclear genetic groups of populations were revealed, corresponding to the two geographic regions separated by the Tengger Desert. However, analysis of cpDNA data did not show significant geographic structure. The most plausible explanation for the discrepancy between our findings based on cpDNA and nDNA data is that A. szechenyi populations experienced long periods of geographic isolation followed by range expansion, which would have promoted generalized recombination of the nuclear genome. Our findings further highlight the important role that the Tengger Desert, together with the Helan Mountains, has played in the evolution of desert plants and the preservation of biodiversity in arid Northwest China.

Demographic inferences and climatic niche modelling shed light on the evolutionary history of the emblematic cold-adapted Apollo butterfly at regional scale

Cold-adapted species escape climate warming by latitudinal and/or altitudinal range shifts, and currently occur in Southern Europe in isolated mountain ranges within "sky islands". Here, we studied the genetic structure of the Apollo butterfly in five such sky islands (above 1,000 m) in France, and infer its demographic history since the last interglacial, using single nucleotide polymorphisms (ddRADseq SNPs). The Auvergne and Alps populations show strong genetic differentiation but not alpine massifs, although separated by deep valleys. Combining three complementary demographic inference methods and species distribution models (SDMs) we show that the LIG period was highly unfavourable for Apollo that probably survived in small population in the highest summits of Auvergne. The population shifted downslope and expanded eastward between LIG and LGM throughout the large climatically suitable Rhône valley between the glaciated summits of Auvergne and Alps. The Auvergne and Alps populations started diverging before the LGM but remained largely connected till the mid-Holocene. Population decline in Auvergne was more gradual but started before (~7 kya vs. 800 ya), and was much stronger with current population size ten times lower than in the Alps. In the Alps, the low genetic structure and limited evidence for isolation by distance suggest a nonequilibrium metapopulation functioning. The core Apollo population experienced cycles of contraction-expansion with climate fluctuations with largely interconnected populations overtime according to a "metapopulation-pulsar" functioning. This study demonstrates the power of combining demographic inferences and SDMs to determine past and future evolutionary trajectories of an endangered species at a regional scale.

Are members of the Anopheles fluviatilis complex conspecific?

Anopheles fluviatilis sensu lato, a primary malaria vector in India, has been identified to be comprised of four cryptic species, provisionally designated as species S, T, U and V. However, Kumar et al. (Mol Ecol Resour, 2013;13:354-61) considered all of the then known three members of this species complex (S, T and U) conspecific. The specific status of species S and T was refuted based on the lack of sufficient barcode gap in mitochondrial-CO1 and the perceived presence of heterozygotes in populations as detected through one of the two species-specific PCR assays employed for the cryptic species identification. The existence of species U was refuted claiming that earlier investigations have already refuted their existence. Here we discuss problems associated with the CO1-based barcode approach for delimitation of cryptic species, the perceived heterozygosity between species S and T based on a species-specific PCR assay, and interpretation of published reports. We demonstrated that fixed differences do exist in the ITS2-rDNA sequence of species S and T with no evidence of heterozygotes in sympatric populations and, that the observed heterozygosity by Kumar et al. in the ITS2-based species diagnostic PCR is due to the high mispriming tendency of the T-specific primer with species S. We infer that mitochondrial DNA-based 'barcoding gap', an arbitrary threshold recommended for species delimitation, alone, is inadequate to delimit the members of An. fluviatilis complex.

Evidence for multiple colonisations and Wolbachia infections shaping the genetic structure of the widespread butterfly Polyommatus icarus in the British Isles

The paradigm of isolation in southern refugia during glacial periods followed by expansions during interglacials, producing limited genetic differentiation in northern areas, dominates European phylogeography. However, the existence of complex structured populations in formerly glaciated areas, and islands connected to mainland areas during glacial maxima, call for alternative explanations. We reconstructed the mtDNA phylogeography of the widespread Polyommatus Icarus butterfly with an emphasis on the formerly glaciated and connected British Isles. We found distinct geographical structuring of CO1 haplogroups, with an ancient lineage restricted to the marginal European areas, including Northern Scotland and Outer Hebrides. Population genomic analyses, using ddRADSeq genomic markers, also reveal substantial genetic structuring within Britain. However, there is negligble mito-nuclear concordance consistent with independent demographic histories of mitochondrial versus nuclear DNA. While mtDNA-Wolbachia associations in northern Britain could account for the geographic structuring of mtDNA across most of the British Isles, for nuclear DNA markers (derived from ddRADseq data) butterflies from France cluster between northern and southern British populations - an observation consistent with a scenario of multiple recolonisation. Taken together our results suggest that contemporary mtDNA structuring in the British Isles (and potentially elsewhere in Europe) largely results from Wolbachia infections, however, nuclear genomic structuring suggests a history of at least two distinct colonisations. This two-stage colonisation scenario has previously been put forth to explain genetic diversity and structuring in other British flora and fauna. Additionally, we also present preliminary evidence for potential Wolbachia-induced feminization in the Outer Hebrides.

Evolutionary history of Sundaland shrews (Eulipotyphla: Soricidae: Crocidura) with a focus on Borneo

The hyperdiverse shrew genus Crocidura is one of few small mammal genera distributed across Sundaland and all of its boundaries. This represents a rare opportunity to study the geological history of this region through the evolutionary history of these shrews. We generate a phylogeny of all recognized species of Sundaland Crocidura and show that most speciation events took place during the Pleistocene, prior to the inundation of the Sunda Shelf around 400 000 years ago. We find east–west differentiation within two separate lineages on Borneo, and that the current taxonomy of its two endemic species does not reflect evolutionary history, but ecophenotypic variation of plastic traits related to elevation. Sulawesi shrews are monophyletic, with a single notable exception: the black-footed shrew (C. nigripes). We show that the black-footed shrew diverged from its relatives on Borneo recently, suggesting a human-assisted breach of Wallace’s line. Overall, the number of Crocidura species, especially on Borneo, probably remains an underestimate.

Integrative analysis supports a new species of the Oecomys catherinae complex (Rodentia, Cricetidae) from Amazonia

The sigmodontine rodent Oecomys catherinae has been recognized as a species complex that includes five different lineages distributed throughout Brazil south of the Amazon river, named in literature according to their geographic position (central, eastern, northern, western, and westernmost). The western lineage also has been referred as Oecomys aff. catherinae in the literature, because it may represent a distinct species. Herein, we used molecular and morphological data to investigate the possible species distinctiveness within the O. catherinae complex. From molecular analyses, we recovered the western lineage as a monophyletic group with high genetic divergence from the other lineages of the complex. This lineage can be differentiated from other lineages of the O. catherinae complex and other congeners through a combination of mensural and morphological characters that includes: smaller size; posterior margin of nasals usually surpassing the maxillary–frontal suture; supraorbital crest moderately developed, similar to the temporal crest; parietal slightly expanded laterally; subsquamosal fenestra present; alisphenoid strut absent; and anterior cingulum of M1 present. The western lineage occurs in the southern Amazonia, in sympatry with the westernmost lineage and other species of Oecomys. Herein we describe this lineage as a new species, increasing to 19 the number of species within Oecomys.

Genetic and morphological variation of Oxymycterus (Rodentia: Sigmodontinae) in the Brazilian Atlantic Forest

We present a new assessment of the genetic and morphological variation within Oxymycterus quaestor Thomas, 1903, which currently includes the junior synonyms O. judex Thomas, 1909 and O. misionalis Sanborn, 1931. We integrate distinct lines of evidence, including variation of mitochondrial (Cytochrome b [Cytb]) and nuclear (intron 7 of beta fibrinogen gene [Fgb]) sequences, and the assessment of skull quantitative traits based on geometric morphometrics, throughout the Atlantic Forest of Southeastern-Southern Brazil, Argentina, and Paraguay. Phylogenetic relationships based on Cytb indicate that O. quaestor is structured in four well-supported clades (lineages A–D), one of them (lineage C) including topotypes of a previously associated nominal form (O. judex). However, these Cytb lineages exhibit lower levels of differentiation based on the Fgb locus, and are not recovered in the genealogies of this nuclear marker, representing a case of mitonuclear discordance. The Cytb lineages also broadly overlapped in the morphospace both in skull shape and size, which sustain the current wider concept of O. quaestor as one single young species (0.947 Myr) that is recently expanding, and ultimately branching out, in the Atlantic Forest.

Species complex diversification by host plant use in an herbivorous insect: The source of Puerto Rican cactus mealybug pest and implications for biological control

Cryptic taxa have often been observed in the form of host‐associated species that diverged as the result of adaptation to alternate host plants. Untangling cryptic diversity in species complexes that encompass invasive species is a mandatory task for pest management. Moreover, investigating the evolutionary history of a species complex may help to understand the drivers of their diversification. The mealybug Hypogeococcus pungens was believed to be a polyphagous species from South America and has been reported as a pest devastating native cacti in Puerto Rico, also threatening cactus diversity in the Caribbean and North America. There is neither certainty about the identity of the pest nor the source population from South America. Recent studies pointed to substantial genetic differentiation among local populations, suggesting that H. pungens is a species complex. In this study, we used a combination of genome‐wide SNPs and mtDNA variation to investigate species diversity within H. pungens sensu lato to establish host plant ranges of each one of the putative members of the complex, to evaluate whether the pattern of host plant association drove diversification in the species complex, and to determine the source population of the Puerto Rican cactus pest. Our results suggested that H. pungens comprises at least five different species, each one strongly associated with specific host plants. We also established that the Puerto Rican cactus pest derives from southeastern Brazilian mealybugs. This is an important achievement because it will help to design reliable strategies for biological control using natural enemies of the pest from its native range.

A thirteen-million-year divergence between two lineages of Indonesian coelacanths

Coelacanth fishes of the genus Latimeria are the only surviving representatives of a basal lineage of vertebrates that originated more than 400 million years ago. Yet, much remains to be unveiled about the diversity and evolutionary history of these ‘living fossils’ using new molecular data, including the possibility of ‘cryptic’ species or unknown lineages. Here, we report the discovery of a new specimen in eastern Indonesia allegedly belonging to the species L. menadoensis. Although this specimen was found about 750 km from the known geographical distribution of the species, we found that the molecular divergence between this specimen and others of L. menadoensis was great: 1.8% compared to 0.04% among individuals of L. chalumnae, the other living species of coelacanth. Molecular dating analyses suggested a divergence date of ca. 13 million years ago between the two populations of Indonesian coelacanths. We elaborate a biogeographical scenario to explain the observed genetic divergence of Indonesian coelacanth populations based on oceanic currents and the tectonic history of the region over Miocene to recent. We hypothesize that several populations of coelacanths are likely to live further east of the present capture location, with potentially a new species that remains to be described. Based on this, we call for an international effort to take appropriate measures to protect these fascinating but vulnerable vertebrates which represent among the longest branches on the Tree of Life.