MITOCHONDRIAL INTROGRESSION AND INCOMPLETE LINEAGE SORTING THROUGH SPACE AND TIME: PHYLOGENETICS OF CROTAPHYTID LIZARDS

Mito-nuclear discordance across a recent contact zone for California voles

To examine the processes that maintain genetic diversity among closely related taxa, we investigated the dynamics of introgression across a contact zone between two lineages of California voles (Microtus californicus). We tested the prediction that introgression of nuclear loci would be greater than that for mitochondrial loci, assuming ongoing gene flow across the contact zone. We also predicted that genomic markers would show a mosaic pattern of differentiation across this zone, consistent with genomes that are semi-permeable. Using mitochondrial cytochrome b sequences and genome-wide loci developed via ddRAD-seq, we analyzed genetic variation for 10 vole populations distributed along the central California coast; this transect included populations from within the distributions of both parental lineages as well as the putative contact zone. Our analyses revealed that (1) the two lineages examined are relatively young, having diverged ca. 8.5-54 kya, (2) voles from the contact zone in Santa Barbara County did not include F1 or early generation backcrossed individuals, and (3) there appeared to be little to no recurrent gene flow across the contact zone. Introgression patterns for mitochondrial and nuclear markers were not concordant; only mitochondrial markers revealed evidence of introgression, putatively due to historical hybridization. These differences in genetic signatures are intriguing given that the contact zone occurs in a region of continuous vole habitat, with no evidence of past or present physical barriers. Future studies that examine specific isolating mechanisms, such as microhabitat use and mate choice, will facilitate our understanding of how genetic boundaries are maintained in this system.

Synchronous diversification of Sulawesi's iconic artiodactyls driven by recent geological events

The high degree of endemism on Sulawesi has previously been suggested to have vicariant origins, dating back to 40 Ma. Recent studies, however, suggest that much of Sulawesi's fauna assembled over the last 15 Myr. Here, we test the hypothesis that more recent uplift of previously submerged portions of land on Sulawesi promoted diversification and that much of its faunal assemblage is much younger than the island itself. To do so, we combined palaeogeographical reconstructions with genetic and morphometric datasets derived from Sulawesi's three largest mammals: the babirusa, anoa and Sulawesi warty pig. Our results indicate that although these species most likely colonized the area that is now Sulawesi at different times (14 Ma to 2–3 Ma), they experienced an almost synchronous expansion from the central part of the island. Geological reconstructions indicate that this area was above sea level for most of the last 4 Myr, unlike most parts of the island. We conclude that emergence of land on Sulawesi (approx. 1–2 Myr) may have allowed species to expand synchronously. Altogether, our results indicate that the establishment of the highly endemic faunal assemblage on Sulawesi was driven by geological events over the last few million years.

Cenozoic aridization in Central Eurasia shaped diversification of toad-headed agamas (Phrynocephalus; Agamidae, Reptilia)

We hypothesize the phylogenetic relationships of the agamid genus Phrynocephalus to assess how past environmental changes shaped the evolutionary and biogeographic history of these lizards and especially the impact of paleogeography and climatic factors. Phrynocephalus is one of the most diverse and taxonomically confusing lizard genera. As a key element of Palearctic deserts, it serves as a promising model for studies of historical biogeography and formation of arid habitats in Eurasia. We used 51 samples representing 33 of 40 recognized species of Phrynocephalus covering all major areas of the genus. Molecular data included four mtDNA (COI, ND2, ND4, Cytb; 2,703 bp) and four nuDNA protein-coding genes (RAG1, BDNF, AKAP9, NKTR; 4,188 bp). AU-tests were implemented to test for significant differences between mtDNA- and nuDNA-based topologies. A time-calibrated phylogeny was estimated using a Bayesian relaxed molecular clock with nine fossil calibrations. We reconstructed the ancestral area of origin, biogeographic scenarios, body size, and the evolution of habitat preference. Phylogenetic analyses of nuDNA genes recovered a well-resolved and supported topology. Analyses detected significant discordance with the less-supported mtDNA genealogy. The position of Phrynocephalus mystaceus conflicted greatly between the two datasets. MtDNA introgression due to ancient hybridization best explained this result. Monophyletic Phrynocephalus contained three main clades: (I) oviparous species from south-western and Middle Asia; (II) viviparous species of Qinghai–Tibetan Plateau (QTP); and (III) oviparous species of the Caspian Basin, Middle and Central Asia. Phrynocephalus originated in late Oligocene (26.9 Ma) and modern species diversified during the middle Miocene (14.8–13.5 Ma). The reconstruction of ancestral areas indicated that Phrynocephalus originated in Middle East–southern Middle Asia. Body size miniaturization likely occurred early in the history of Phrynocephalus. The common ancestor of Phrynocephalus probably preferred sandy substrates with the inclusion of clay or gravel. The time of Agaminae radiation and origin of Phrynocephalus in the late Oligocene significantly precedes the landbridge between Afro-Arabia and Eurasia in the Early Miocene. Diversification of Phrynocephalus coincides well with the mid-Miocene climatic transition when a rapid cooling of climate drove progressing aridification and the Paratethys salinity crisis. These factors likely triggered the spreading of desert habitats in Central Eurasia, which Phrynocephalus occupied. The origin of the viviparous Tibetan clade has been associated traditionally with uplifting of the QTP; however, further studies are needed to confirm this. Progressing late Miocene aridification, the decrease of the Paratethys Basin, orogenesis, and Plio–Pleistocene climate oscillations likely promoted further diversification within Phrynocephalus. We discuss Phrynocephalus taxonomy in scope of the new analyses.

The effects of allospecific mitochondrial genome on the fitness of northern redbelly dace (Chrosomus eos)

Instantaneous mitochondrial introgression events allow the disentangling of the effects of hybridization from those of allospecific mtDNA. Such process frequently occurred in the fish Chrosomus eos, resulting in cybrid individuals composed of a C. eos nuclear genome but with a C. neogaeus mtDNA. This provides a valuable model to address the fundamental question: How well do introgressed individuals perform in their native environment? We infer where de novo production of cybrids occurred to discriminate native environments from those colonized by cybrids in 25 sites from two regions (West-Qc and East-Qc) in Quebec (Canada). We then compared the relative abundance of wild types and cybrids as a measure integrating both fitness and de novo production of cybrids. According to mtDNA variation, 12 introgression events are required to explain the diversity of cybrids. Five cybrid lineages could not be associated with in situ introgression events. This includes one haplotype carried by 93% of the cybrids expected to have colonized West-Qc. These cybrids also displayed a nearly complete allopatric distribution with wild types. We still inferred de novo production of cybrids at seven sites, that accounted for 70% of the cybrids in East-Qc. Wild-type and cybrid individuals coexist in all East-Qc sites while cybrids were less abundant. Allopatry of cybrids restricted to the postglacial expansion suggests the existence of higher fitness for cybrids in specific conditions, allowing for the colonization of different environments and expanding the species' range. However, allospecific mtDNA does not provide a higher fitness to cybrids in their native environment compared to wild types, making the success of an introgressed lineage uncertain.

Genetic and morphological support for possible sympatric origin of fish from subterranean habitats

Two blind Iran cave barbs, Garra typhlops and Garra lorestanensis, exist in sympatry in a single subterranean habitat, raising the hypothesis that they may represent a case of sympatric speciation following a colonization event. Their different mental disc forms have prompted some authors to propose the alternative hypothesis of two separate colonization events. In this study, we analysed a genome-wide panel of 11,257 SNPs genotyped by means of genotyping-by-sequencing combined with mitochondrial cytochrome c oxidase sub-unit I sequence data, field observations and morphological traits to test these two hypotheses. Field data suggest some degree of ecological divergence despite some possible niche overlap such that hybridization is possible. According to both nuclear and mtDNA data, the cave barb species are monophyletic with close phylogenetic relationships with Garra gymnothorax from the Karun-Dez and Karkheh river basins. The historical demography analysis revealed that a model of Isolation-with-Migration (IM) best fitted the data, therefore better supporting a scenario of sympatric origin than that of allopatric isolation followed by secondary contact. Overall, our results offer stronger support to the hypothesis that speciation in the subterranean habitat could have occurred in sympatry following a colonization event from the Karun-Dez-Karkheh basins in the Zagros Mountains of Iran.

Comparative Phylogenomic Assessment of Mitochondrial Introgression among Several Species of Chipmunks (Tamias)

Many species are not completely reproductively isolated, resulting in hybridization and genetic introgression. Organellar genomes, such as those derived from mitochondria (mtDNA) and chloroplasts, introgress frequently in natural systems; however, the forces shaping patterns of introgression are not always clear. Here, we investigate extensive mtDNA introgression in western chipmunks, focusing on species in the Tamias quadrivittatus group from the central and southern Rocky Mountains. Specifically, we investigate the role of selection in driving patterns of introgression. We sequenced 51 mtDNA genomes from six species and combine these sequences with other published genomic data to yield annotated mitochondrial reference genomes for nine species of chipmunks. Genomic characterization was performed using a series of molecular evolutionary and phylogenetic analyses to test protein-coding genes for positive selection. We fit a series of maximum likelihood models using a model-averaging approach, assessed deviations from neutral expectations, and performed additional tests to search for codons under the influence of selection. We found no evidence for positive selection among these genomes, suggesting that selection has not been the driving force of introgression in these species. Thus, extensive mtDNA introgression among several species of chipmunks likely reflects genetic drift of introgressed alleles in historically fluctuating populations.

Geographic patterns of phenotypic diversity in incipient species of North American blister beetles (Coleoptera: Meloidae) are not determined by species niches, but driven by demography along the speciation process

The Epicauta stigmata complex is a group of blister beetles composed of three parapatric or sympatric species that occur in central Mexico to southern USA: E. stigmata, E. uniforma and E. melanochroa. These species are morphologically very similar, and are mainly distinguished by body colour differences. Here we assessed whether phenotypic divergence in coloration patterns define evolutionary units within the complex. We studied the phylogenetic relationships, demographic history and concordances between morphological and ecological traits in the group. The complex apparently had a demographic history of recent population expansion during the last glaciation period 75000 to 9500 years ago. The three species show no reciprocal monophyly, and thus their allospecificity was not confirmed. The current distribution of haplotypes and the genetic divergences in these taxa can be explained by either recent mitochondrial introgression events caused by hybridisation or by incomplete lineage sorting. Colour pattern differences in the complex are not likely a product of local selection acting over a common genetic background. We suggest that phenotypic divergence in colour patterns during an incipient speciation process might be seen as an enhancing factor of cohesion within each of the three evolutionary units.

The genomic ancestry, landscape genetics and invasion history of introduced mice in New Zealand

The house mouse (Mus musculus) provides a fascinating system for studying both the genomic basis of reproductive isolation, and the patterns of human-mediated dispersal. New Zealand has a complex history of mouse invasions, and the living descendants of these invaders have genetic ancestry from all three subspecies, although most are primarily descended from M. m. domesticus. We used the GigaMUGA genotyping array (approximately 135 000 loci) to describe the genomic ancestry of 161 mice, sampled from 34 locations from across New Zealand (and one Australian city-Sydney). Of these, two populations, one in the south of the South Island, and one on Chatham Island, showed complete mitochondrial lineage capture, featuring two different lineages of M. m. castaneus mitochondrial DNA but with only M. m. domesticus nuclear ancestry detectable. Mice in the northern and southern parts of the North Island had small traces (approx. 2-3%) of M. m. castaneus nuclear ancestry, and mice in the upper South Island had approximately 7-8% M. m. musculus nuclear ancestry including some Y-chromosomal ancestry-though no detectable M. m. musculus mitochondrial ancestry. This is the most thorough genomic study of introduced populations of house mice yet conducted, and will have relevance to studies of the isolation mechanisms separating subspecies of mice.

Integrating phylogenomic and population genomic patterns in avian lice provides a more complete picture of parasite evolution

Parasite diversity accounts for most of the biodiversity on earth, and is shaped by many processes (e.g., cospeciation, host switching). To identify the effects of the processes that shape parasite diversity, it is ideal to incorporate both deep (phylogenetic) and shallow (population) perspectives. To this end, we developed a novel workflow to obtain phylogenetic and population genetic data from whole genome sequences of body lice parasitizing New World ground-doves. Phylogenies from these data showed consistent, highly resolved species-level relationships for the lice. By comparing the louse and ground-dove phylogenies, we found that over long-term evolutionary scales their phylogenies were largely congruent. Many louse lineages (both species and populations) also demonstrated high host-specificity, suggesting ground-dove divergence is a primary driver of their parasites' diversity. However, the few louse taxa that are generalists are structured according to biogeography at the population level. This suggests dispersal among sympatric hosts has some effect on body louse diversity, but over deeper time scales the parasites eventually sort according to host species. Overall, our results demonstrate that multiple factors explain the patterns of diversity in this group of parasites, and that the effects of these factors can vary over different evolutionary scales. The integrative approach we employed was crucial for uncovering these patterns, and should be broadly applicable to other studies.

The on-again, off-again relationship between mitochondrial genomes and species boundaries

The study of reproductive isolation and species barriers frequently focuses on mitochondrial genomes and has produced two alternative and almost diametrically opposed narratives. On one hand, mtDNA may be at the forefront of speciation events, with co-evolved mitonuclear interactions responsible for some of the earliest genetic incompatibilities arising among isolated populations. On the other hand, there are numerous cases of introgression of mtDNA across species boundaries even when nuclear gene flow is restricted. We argue that these seemingly contradictory patterns can result from a single underlying cause. Specifically, the accumulation of deleterious mutations in mtDNA creates a problem with two alternative evolutionary solutions. In some cases, compensatory or epistatic changes in the nuclear genome may ameliorate the effects of mitochondrial mutations, thereby establishing coadapted mitonuclear genotypes within populations and forming the basis of reproductive incompatibilities between populations. Alternatively, populations with high mitochondrial mutation loads may be rescued by replacement with a more fit, foreign mitochondrial haplotype. Coupled with many nonadaptive mechanisms of introgression that can preferentially affect cytoplasmic genomes, this form of adaptive introgression may contribute to the widespread discordance between mitochondrial and nuclear genealogies. Here, we review recent advances related to mitochondrial introgression and mitonuclear incompatibilities, including the potential for cointrogression of mtDNA and interacting nuclear genes. We also address an emerging controversy over the classic assumption that selection on mitochondrial genomes is inefficient and discuss the mechanisms that lead lineages down alternative evolutionary paths in response to mitochondrial mutation accumulation.

DNA barcoding revises a misidentification on mossy frog: new record and distribution extension of Theloderma corticale Boulenger, 1903 (Amphibia: Anura: Rhacophoridae)

As an endangered animal group, mossy frog (genus Theloderma) has attracted the attention of biologists and wildlife conservationists. Clarifying the taxonomic status and distribution of each species in Theloderma is important to determine the conservation status for each species, establish appropriate conservation strategies and probe the speciation process. Recently, we discovered a medium-sized species of mossy frog of the genus Theloderma in April 2015 during municipal surveys of amphibians in Dayao Mountain of Jinxiu. It was collected from the water-filled tree cavities. However, there remains some uncertainty about the species determination of the mossy frog in the Yinshan station of Dayao Mountain in Guangxi Province, China. Previously, the mossy frog in Guangxi Province was recognized as Th. kwangsiense. In order to clarify the species status of the mossy frog obtained from Guangxi, we sequenced 2414 bp of the 12S and 16S genes in the sample collected from the Dayao Mountain. Combining all the sequence in NCBI, genetic analyses from the data suggest that the sample from the Dayao Mountain is Th. corticale rather than Th. kwangsiense. It is most likely that the most previous studies had wrong species identification. And, it is the first time we use DNA barcoding to prove that the species obtained from Guangxi is a new distribution of Th. corticale.

Evidence for environmental and ecological selection in a microbe with no geographic limits to gene flow

The ability for organisms to disperse throughout their environment is thought to strongly influence population structure and thus evolution of diversity within species. A decades-long debate surrounds processes that generate and support high microbial diversity, particularly in the ocean. The debate concerns whether diversification occurs primarily through geographic partitioning (where distance limits gene flow) or through environmental selection, and remains unresolved due to lack of empirical data. Here we show that gene flow in a diatom, an ecologically important eukaryotic microbe, is not limited by global-scale geographic distance. Instead, environmental and ecological selection likely play a more significant role than dispersal in generating and maintaining diversity. We detected significantly diverged populations (F_ST > 0.130) and discovered temporal genetic variability at a single site that was on par with spatial genetic variability observed over distances of 15,000 km. Relatedness among populations was decoupled from geographic distance across the global ocean and instead, correlated significantly with water temperature and whole-community chlorophyll a Correlations with temperature point to the importance of environmental selection in structuring populations. Correlations with whole-community chlorophyll a, a proxy for autotrophic biomass, suggest that ecological selection via interactions with other plankton may generate and maintain population genetic structure in marine microbes despite global-scale dispersal. Here, we provide empirical evidence for global gene flow in a marine eukaryotic microbe, suggesting that everything holds the potential to be everywhere, with environmental and ecological selection rather than geography or dispersal dictating the structure and evolution of diversity over space and time.

A combination of long term fragmentation and glacial persistence drove the evolutionary history of the Italian wall lizard Podarcis siculus

BACKGROUND

The current distribution of genetic diversity is the result of a vast array of microevolutionary processes, including short-term demographic and ecological mechanisms and long-term allopatric isolation in response to Quaternary climatic fluctuations. We investigated past processes that drove the population differentiation and spatial genetic distribution of the Italian wall lizard Podarcis siculus by means of sequences of mitochondrial cytb (n = 277 from 115 localities) and nuclear mc1r and β-fibint7genes (n = 262 and n = 91, respectively) from all its distribution range. The pattern emerging from the genetic data was compared with current and past (last glacial maximum) species distribution modeling (SDM).

RESULTS

We identified seven deeply divergent parapatric clades which presumably remained isolated in different refugia scattered mainly throughout the Tyrrhenian coast. Conversely, the Adriatic coast showed only two haplogroups with low genetic variability. These results appear to agree with the SDM prediction at the last glacial maximum (LGM) indicating a narrow area of habitat suitability along the Tyrrhenian coast and much lower suitability along the Adriatic one. However, the considerable land exposure of the Adriatic coastline favored a glacial colonization of the Balkan Peninsula.

CONCLUSIONS

Our population-level historical demography showed a common trend consistent with glacial expansions and regional persistence during the last glacial maximum. This complex genetic signature appears to be inconsistent with the expectation of the expansion-contraction model and post-LGM (re)colonizations from southern refugia. Hence it is one of an increasing number of cases in which these assumptions are not met, indicating that long-term fragmentation and pre-LGM events such as glacial persistence were more prominent in shaping genetic variation in this temperate species.

Molecular phylogeny, biogeography and ecological niche modelling of Cardiocrinum (Liliaceae): insights into the evolutionary history of endemic genera distributed across the Sino-Japanese floristic region

BACKGROUND AND AIMS

The patterns of evolutionary assembly in the Sino-Japanese floristic region (SJFR) remain largely unknown due to a lack of integrative multidimensional studies throughout the region. To address this issue, we elucidated the evolutionary history of Cardiocrinum (Liliaceae), a genus containing four taxa distributed across the SJFR.

METHODS

Fifty-four populations were sampled throughout the geographical range of Cardiocrinum to assess genetic structure, analyse phylogenetic relationships and reconstruct ancestral area based on six chloroplast DNA (cpDNA) fragments and three low copy nuclear genes (LCNG). Ecological niche modelling was used to examine the potential range shifts of Cardiocrinum in response to climatic change.

KEY RESULTS

The molecular data showed high genetic similarity in the cpDNA (98·37 %) and LCNG (94·53 %) sequences. The biogeographical analyses revealed that the ancestor of Cardiocrinum diversified during the late Miocene (approx. 7·32 Mya) in Central China. The ancestor of the C. giganteum lineage dispersed westward to the Himalayas and south-west China with the split between C. giganteum and C. giganteum var. yunnanense occurring around 4·11 Mya consistent with the period of orogeny of the Hengduan Mountains. Some populations of the C. cathayanum lineage dispersed eastward to south Japan via the land bridge approx. 4·97 Mya, providing opportunities for allopatric speciation of C. cordatum The predicted suitable habitats of Cardiocrinum have become smaller and more fragmented since the Last Glacial Maximum.

CONCLUSIONS

Our study provides evidence of a biogeographical pattern of dispersal from Central China to the Himalayas in the west and Japan in the east for genera distributed across the SJFR, and highlights that the orogeny of the Hengduan Mountains and fluctuations of the sea level of the East China Sea played important roles in promoting species divergence.

Genetic homogenization of the nuclear ITS loci across two morphologically distinct gentians in their overlapping distributions in the Qinghai-Tibet Plateau

Interspecific hybridization and introgression can lead to partial genetic homogenization at certain neutral loci between morphologically distinct species and may obscure the species delimitation using nuclear genes. In this study, we investigated this phenomenon through population genetic survey of two alpine plants (Gentiana siphonantha and G. straminea) in the Qinghai-Tibet Plateau, where the distributions of two species are partly overlapped. We identified two clusters of chloroplast DNA haplotypes which correspond to the two species, and three clusters of ITS ribotypes. In addition to clusters specific to each species, the third ITS cluster, which was most likely derived from hybridization between the other two clusters and subsequent recombination and concerted evolution, was widely shared by two species in their adjacent areas. In contrast to the morphological distinctiveness of the two species, interspecific gene flow possibly led to genetic homogenization at their ITS loci. The new ITS lineage recovered for species in adjacent areas is distinctly different from original lineages found in allopatric areas. These findings may have general implications for our understanding of cryptic changes at some genetic loci caused by interspecific gene flow in the history, and they indicate that species delimitation should be based on a combination of both nuclear and chloroplast DNA sequence variations.

Do genes lie? Mitochondrial capture masks the Red Sea collector urchin's true identity (Echinodermata: Echinoidea: Tripneustes)

Novel COI and bindin sequences of the Red Sea collector echinoid Tripneustes gratilla elatensis are used to show that (1) discordance between mitochondrial and nuclear loci exists in this echinoid genus, (2) Tripneustes gratilla as currently defined possibly comprises a complex of cryptic species, and (3) Red Sea Tripneustes form a genetically distinct clade in the bindin tree, which diverged from other Tripneustes clades at least 2-4million years ago. Morphological reassessment of T. gratilla elatensis shows perfect congruence between identification based on skeletal features and genetic data based on a nuclear marker sequence. Hence the Red Sea Tripneustes subspecies established by Dafni in 1983 is a distinct biological unit. All T. g. elatensis samples analyzed are highly similar to or share mtDNA haplotypes with Philippine T. g. gratilla, as do representatives from other edge-of-range occurrences. This lack of genetic structure in Indo-Pacific Tripneustes is interpreted as a result of wide-spread mitochondrial introgression. New fossil specimens from the Red Sea area confirm the sympatric occurrence of T. g. elatensis and T. g. gratilla in the northern Red Sea during Late Pleistocene, identifying a possible timing for the introgression. In addition, present-day distribution shows a contact zone in the Southern Red Sea (in the Dahlak Archipelago). T. g. elatensis, is yet another example of a Red Sea taxon historically identified as conspecific with its Indo-Pacific relatives, but which turned out to be a morphologically and genetically distinct endemic taxon, suggesting that the level of endemism in the Red Sea may still be underestimated.

Genetic diversity and evolutionary history of the Schizothorax species complex in the Lancang River (upper Mekong)

The genus Schizothorax (Cyprinidae), one of the most diverse genera of ichthyofauna of the Qinghai‐Tibetan Plateau (QTP), is a good candidate for investigating patterns of genetic variation and evolutionary mechanisms. In this study, sequences from the mitochondrial control region, the cytochrome b gene, and two nuclear genes were used to re‐examine the genetic diversity and investigate the evolutionary history of the Schizothorax species complex inhabiting the Lancang River. Three maternal clades were detected in the Schizothorax species complex, but frequent nuclear allele sharing also occurred among the three maternal clades. A discrepancy between topologies of mitochondrial and nuclear loci might result from introgression or/and incomplete lineage sorting. The divergence of the clades of the Schizothorax species complex was closely related to the Late Pliocene and Early Pleistocene orogenesis of the QTP and Southwest Mountains of China. Demographic analyses indicated that the species complex subsequently persisted in situ with stable populations during Pleistocene glacial cycling, which suggested that Pleistocene climate changes did not exert a remarkable influence on the species complex. Our study provides a comprehensive analysis of the genetic diversity and evolutionary history of the Schizothorax species complex in the Lancang River.

Whipworms in humans and pigs: origins and demography

Background

Trichuris suis and T. trichiura are two different whipworm species that infect pigs and humans, respectively. T. suis is found in pigs worldwide while T. trichiura is responsible for nearly 460 million infections in people, mainly in areas of poor sanitation in tropical and subtropical areas. The evolutionary relationship and the historical factors responsible for this worldwide distribution are poorly understood. In this study, we aimed to reconstruct the demographic history of Trichuris in humans and pigs, the evolutionary origin of Trichuris in these hosts and factors responsible for parasite dispersal globally.

Methods

Parts of the mitochondrial nad1 and rrnL genes were sequenced followed by population genetic and phylogenetic analyses. Populations of Trichuris examined were recovered from humans (n = 31), pigs (n = 58) and non-human primates (n = 49) in different countries on different continents, namely Denmark, USA, Uganda, Ecuador, China and St. Kitts (Caribbean). Additional sequences available from GenBank were incorporated into the analyses.

Results

We found no differentiation between human-derived Trichuris in Uganda and the majority of the Trichuris samples from non-human primates suggesting a common African origin of the parasite, which then was transmitted to Asia and further to South America. On the other hand, there was no differentiation between pig-derived Trichuris from Europe and the New World suggesting dispersal relates to human activities by transporting pigs and their parasites through colonisation and trade. Evidence for recent pig transport from China to Ecuador and from Europe to Uganda was also observed from their parasites. In contrast, there was high genetic differentiation between the pig Trichuris in Denmark and China in concordance with the host genetics.

Conclusions

We found evidence for an African origin of T. trichiura which were then transmitted with human ancestors to Asia and further to South America. A host shift to pigs may have occurred in Asia from where T. suis seems to have been transmitted globally by a combination of natural host dispersal and anthropogenic factors.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1325-8) contains supplementary material, which is available to authorized users.

Repeated Reticulate Evolution in North American Papilio machaon Group Swallowtail Butterflies

Hybridization between distinct populations or species is increasingly recognized as an important process for generating biodiversity. However, the interaction between hybridization and speciation is complex, and the diverse evolutionary outcomes of hybridization are difficult to differentiate. Here we characterize potential hybridization in a species group of swallowtail butterflies using microsatellites, DNA sequences, and morphology, and assess whether adaptive introgression or homoploid hybrid speciation was the primary process leading to each putative hybrid lineage. Four geographically separated hybrid populations were identified in the Papilio machaon species group. One distinct mitochondrial DNA clade from P. machaon was fixed in three hybrid taxa (P. brevicauda, P. joanae, and P. m. kahli), while one hybrid swarm (P. zelicaon x machaon) exhibited this hybrid mtDNA clade as well as widespread parental mtDNA haplotypes from both parental species. Microsatellite markers and morphology showed variable admixture and intermediacy, ranging from signatures of prolonged differential introgression from the paternal species (P. polyxenes/P. zelicaon) to current gene flow with both parental species. Divergences of the hybrid lineages dated to early- to mid-Pleistocene, suggesting that repeated glaciations and subsequent range shifts of parental species, particularly P. machaon hudsonianus, facilitated initial hybridization. Although each lineage is distinct, P. joanae is the only taxon with sufficient evidence (ecological separation from parental species) to define it as a homoploid hybrid species. The repetition of hybridization in this group provides a valuable foundation for future research on hybridization, and these results emphasize the potential for hybridization to drive speciation in diverse ways.

Mitochondrial introgression via ancient hybridization, and systematics of the Australian endemic pygopodid gecko genus Delma

Of the more than 1500 species of geckos found across six continents, few remain as unfamiliar as the pygopodids - Family Pygopodidae (Gray, 1845). These gekkotans are limited to Australia (44 species) and New Guinea (2 species), but have diverged extensively into the most ecologically diverse limbless radiation save Serpentes. Current phylogenetic understanding of the family has relied almost exclusively on two works, which have produced and synthesized an immense amount of morphological, geographical, and molecular data. However, current interspecific relationships within the largest genus Delma Gray 1831 are based chiefly upon data from two mitochondrial loci (16s, ND2). Here, we reevaluate the interspecific relationships within the genus Delma using two mitochondrial and four nuclear loci (RAG1, MXRA5, MOS, DYNLL1), and identify points of strong conflict between nuclear and mitochondrial genomic data. We address mito-nuclear discordance, and remedy this conflict by recognizing several points of mitochondrial introgression as the result of ancient hybridization events. Owing to the legacy value and intraspecific informativeness, we suggest the continued use of ND2 as a phylogenetic marker. Results identify strong support for species groups, but relationships among these clades, and the placement of several enigmatic taxa remain uncertain. We suggest a more careful review of Delma australis and the 'northwest Australia' clade. Accurately assessing and addressing species richness and relationships within this endemic Australian Gekkotan genus is relevant for understanding patterns of squamate speciation across the region.

Mitochondrial Genome Analyses Suggest Multiple Trichuris Species in Humans, Baboons, and Pigs from Different Geographical Regions

Background

The whipworms Trichuris trichiura and Trichuris suis are two parasitic nematodes of humans and pigs, respectively. Although whipworms in human and non-human primates historically have been referred to as T. trichiura, recent reports suggest that several Trichuris spp. are found in primates.

Methods and Findings

We sequenced and annotated complete mitochondrial genomes of Trichuris recovered from a human in Uganda, an olive baboon in the US, a hamadryas baboon in Denmark, and two pigs from Denmark and Uganda. Comparative analyses using other published mitochondrial genomes of Trichuris recovered from a human and a porcine host in China and from a françois’ leaf-monkey (China) were performed, including phylogenetic analyses and pairwise genetic and amino acid distances. Genetic and protein distances between human Trichuris in Uganda and China were high (~19% and 15%, respectively) suggesting that they represented different species. Trichuris from the olive baboon in US was genetically related to human Trichuris in China, while the other from the hamadryas baboon in Denmark was nearly identical to human Trichuris from Uganda. Baboon-derived Trichuris was genetically distinct from Trichuris from françois’ leaf monkey, suggesting multiple whipworm species circulating among non-human primates. The genetic and protein distances between pig Trichuris from Denmark and other regions were roughly 9% and 6%, respectively, while Chinese and Ugandan whipworms were more closely related.

Conclusion and Significance

Our results indicate that Trichuris species infecting humans and pigs are phylogenetically distinct across geographical regions, which might have important implications for the implementation of suitable and effective control strategies in different regions. Moreover, we provide support for the hypothesis that Trichuris infecting primates represents a complex of cryptic species with some species being able to infect both humans and non-human primates.

Trichuris trichiura and Trichuris suis are whipworms found in humans and pigs, respectively, causing morbidity in humans and being associated with production losses in pigs. Although Trichuris from non-human primates is attributed to T. trichiura, hence considered the same species as the one infecting humans, several recent reports question this assumption. Morphologically similar parasites that have a wide global distribution and/or those capable of infecting several host species may comprise several ‘hidden’ species. In this study, we sequenced, annotated, and compared the mitochondrial genomes (including published genomes) of Trichuris obtained from different hosts in different geographical regions, including humans (Uganda and China), pigs (China, Uganda, and Denmark) and two types of non-human primates (baboons and françois’ leaf monkey). We found high genetic distinctiveness between human Trichuris from China and Uganda. Likewise, pig Trichuris from Denmark and other regions also showed considerable, although lower, genetic diversity. This suggests that both pig- and human-derived Trichuris may represent different species with potential differences in endemicity, which may have important implications for implementing effective control strategies. Our data also suggests that Trichuris infecting primates comprises several species and may be transmitted from non-human primates to humans.

DNA barcoding and phylogenetic relationships in Anatidae

Mitochondrial cytochrome c oxidase subunit I (COI) has been used as a powerful marker in a variety of phylogenetic studies. According to studies of bird species, the 694-bp sequence of the mitochondrial gene encoding COI is extremely useful for species identification and phylogeny. In the present study, we analyzed the COI barcodes of 79 species from 26 genera belonging to the Anatidae family. Sixty-six species (83.54%) of the species were identified correctly from their DNA barcodes. The remaining 13 species shared barcodes sequences with closely related species. Kimura two-parameter (K2P) distances were calculated between barcodes. The average genetic distance between species was 41 times higher compared to the average genetic distance within species. Neighbor-joining method was used to construct a phylogenetic tree, which grouped all of the genera into three divergent clades. Dendrocygna and Nomonyx + Oxyura were identified as early offshoots of the Anatidae. All the remaining taxa fell into two clades that correspond to the two subfamilies Anserinae and Anatiane. Based on our results, DNA barcoding is an effective molecular tool for Anatidae species identification and phylogenetic inference.

Taxonomic challenges in freshwater fishes: a mismatch between morphology and DNA barcoding in fish of the north-eastern part of the Congo basin

This study evaluates the utility of DNA barcoding to traditional morphology-based species identifications for the fish fauna of the north-eastern Congo basin. We compared DNA sequences (COI) of 821 samples from 206 morphologically identified species. Best match, best close match and all species barcoding analyses resulted in a rather low identification success of 87.5%, 84.5% and 64.1%, respectively. The ratio 'nearest-neighbour distance/maximum intraspecific divergence' was lower than 1 for 26.1% of the samples, indicating possible taxonomic problems. In ten genera, belonging to six families, the number of species inferred from mtDNA data exceeded the number of species identified using morphological features; and in four cases indications of possible synonymy were detected. Finally, the DNA barcodes confirmed previously known identification problems within certain genera of the Clariidae, Cyprinidae and Mormyridae. Our results underscore the large number of taxonomic problems lingering in the taxonomy of the fish fauna of the Congo basin and illustrate why DNA barcodes will contribute to future efforts to compile a reliable taxonomic inventory of the Congo basin fish fauna. Therefore, the obtained barcodes were deposited in the reference barcode library of the Barcode of Life Initiative.

Phylogeography of Pogonomyrmex barbatus and P. rugosus harvester ants with genetic and environmental caste determination

We present a phylogeographic study of at least six reproductively isolated lineages of new world harvester ants within the Pogonomyrmex barbatus and P. rugosus species group. The genetic and geographic relationships within this clade are complex: Four of the identified lineages show genetic caste determination (GCD) and are divided into two pairs. Each pair has evolved under a mutualistic system that necessitates sympatry. These paired lineages are dependent upon one another because their GCD requires interlineage matings for the production of F1 hybrid workers, and intralineage matings are required to produce queens. This GCD system maintains genetic isolation among these interdependent lineages, while simultaneously requiring co-expansion and emigration as their distributions have changed over time. It has also been demonstrated that three of these four GCD lineages have undergone historical hybridization, but the narrower sampling range of previous studies has left questions on the hybrid parentage, breadth, and age of these groups. Thus, reconstructing the phylogenetic and geographic history of this group allows us to evaluate past insights and hypotheses and to plan future inquiries in a more complete historical biogeographic context. Using mitochondrial DNA sequences sampled across most of the morphospecies' ranges in the U.S.A. and Mexico, we conducted a detailed phylogeographic study. Remarkably, our results indicate that one of the GCD lineage pairs has experienced a dramatic range expansion, despite the genetic load and fitness costs of the GCD system. Our analyses also reveal a complex pattern of vicariance and dispersal in Pogonomyrmex harvester ants that is largely concordant with models of late Miocene, Pliocene, and Pleistocene range shifts among various arid-adapted taxa in North America.

Genetic evaluation of the evolutionary distinctness of a federally endangered butterfly, Lange's Metalmark

Background

The Mormon Metalmark (Apodemia mormo) species complex occurs as isolated and phenotypically variable colonies in dryland areas across western North America. Lange’s Metalmark, A. m. langei, one of the 17 subspecies taxonomically recognized in the complex, is federally listed under the U.S. Endangered Species Act of 1973. Metalmark taxa have traditionally been described based on phenotypic and ecological characteristics, and it is unknown how well this nomenclature reflects their genetic and evolutionary distinctiveness. Genetic variation in six microsatellite loci and mitochondrial cytochrome oxidase subunit I sequence was used to assess the population structure of the A. mormo species complex across 69 localities, and to evaluate A. m. langei’s qualifications as an Evolutionarily Significant Unit.

Results

We discovered substantial genetic divergence within the species complex, especially across the Continental Divide, with population genetic structure corresponding more closely with geographic proximity and local isolation than with taxonomic divisions originally based on wing color and pattern characters. Lange’s Metalmark was as genetically divergent as several other locally isolated populations in California, and even the unique phenotype that warranted subspecific and conservation status is reminiscent of the morphological variation found in some other populations.

Conclusions

This study is the first genetic treatment of the A. mormo complex across western North America and potentially provides a foundation for reassessing the taxonomy of the group. Furthermore, these results illustrate the utility of molecular markers to aid in demarcation of biological units below the species level. From a conservation point of view, Apodemia mormo langei’s diagnostic taxonomic characteristics may, by themselves, not support its evolutionary significance, which has implications for its formal listing as an Endangered Species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0354-9) contains supplementary material, which is available to authorized users.

Interbreeding among deeply divergent mitochondrial lineages in the American cockroach (Periplaneta americana)

DNA barcoding promises to be a useful tool to identify pest species assuming adequate representation of genetic variants in a reference library. Here we examined mitochondrial DNA barcodes in a global urban pest, the American cockroach (Periplaneta americana). Our sampling effort generated 284 cockroach specimens, most from New York City, plus 15 additional U.S. states and six other countries, enabling the first large-scale survey of P. americana barcode variation. Periplaneta americana barcode sequences (n = 247, including 24 GenBank records) formed a monophyletic lineage separate from other Periplaneta species. We found three distinct P. americana haplogroups with relatively small differences within (≤0.6%) and larger differences among groups (2.4%-4.7%). This could be interpreted as indicative of multiple cryptic species. However, nuclear DNA sequences (n = 77 specimens) revealed extensive gene flow among mitochondrial haplogroups, confirming a single species. This unusual genetic pattern likely reflects multiple introductions from genetically divergent source populations, followed by interbreeding in the invasive range. Our findings highlight the need for comprehensive reference databases in DNA barcoding studies, especially when dealing with invasive populations that might be derived from multiple genetically distinct source populations.

Montane and coastal species diversification in the economically important Mexican grasshopper genus Sphenarium (Orthoptera: Pyrgomorphidae)

The genus Sphenarium (Pyrgomorphidae) is a small group of grasshoppers endemic to México and Guatemala that are economically and culturally important both as a food source and as agricultural pests. However, its taxonomy has been largely neglected mainly due to its conserved interspecific external morphology and the considerable intraspecific variation in colour pattern of some taxa. Here we examined morphological as well as mitochondrial and nuclear DNA sequence data to assess the species boundaries and evolutionary history in Sphenarium. Our morphological identification and DNA sequence-based species delimitation, carried out with three different approaches (DNA barcoding, general mixed Yule-coalescent model, Bayesian species delimitation), all recovered a higher number of putative species of Sphenarium than previously recognised. We unambiguously delimit seven species, and between five and ten additional species depending on the data/method analysed. Phylogenetic relationships within the genus strongly support two main clades, one exclusively montane, the other coastal. Divergence time estimates suggest late Miocene to Pliocene ages for the origin and most of the early diversification events in the genus, which were probably influenced by the formation of the Trans-Mexican Volcanic Belt. A series of Pleistocene events could have led to the current species diversification in both montane and coastal regions. This study not only reveals an overlooked species richness for the most popular edible insect in Mexico, but also highlights the influence of the dynamic geological and climatic history of the region in shaping its current diversity.

Hybridization and mitochondrial genome introgression between Rana chensinensis and R. kukunoris

Mitochondrial genome (mito-genome) introgression among metazoans is commonplace, and several biological processes may promote such introgression. We examined two proposed processes for the mito-genome introgression between Rana chensinensis and R. kukunoris: natural hybridization and sex-biased dispersal. We sampled 477 individuals from 28 sites in the potential hybrid zone in the western Tsinling Mountains. Mitochondrial gene (cyt-b) trees were used to examine the introgression events. Microsatellite DNA loci, cyt-b and morphological data were used to identify hybrids and to examine the extent of natural hybridization. We detected rampant bidirectional introgressions, both ancient and recent, between the two species. Furthermore, we found a wide hybrid zone, and frequent and asymmetric hybridization. The hybrid zone cline analysis revealed a clear mitochondrial-nuclear discordance; while most nuclear markers displayed similar and steep clines, cyt-b had a displaced cline centre and a more gradual and wider cline. We also detected strong and asymmetric historical maternal gene flow across the hybrid zone. This widespread hybridization and detected low mito-nuclear conflicts may, at least partially, explain the high frequency of introgression. Lastly, microsatellite data and population genetic methods were used to assess sex-biased dispersal. A weak pattern of female-biased dispersal was detected in both species, suggesting it may not play an important role in the observed introgression. Our data are consistent with the hybridization hypothesis, but support for the sex-biased dispersal hypothesis is weak. We further suggest that selective advantages of the R. kukunoris-type mito-genome in thermal adaptation may also contribute to the introgression between the two species.