Distinguishing between incomplete lineage sorting and genomic introgressions: complete fixation of allospecific mitochondrial DNA in a sexually reproducing fish (Cobitis; Teleostei), despite clonal reproduction of hybrids

Sperm-dependent asexual species and their role in ecology and evolution

Sexual reproduction is the primary mode of reproduction in eukaryotes, but some organisms have evolved deviations from classical sex and switched to asexuality. These asexual lineages have sometimes been viewed as evolutionary dead ends, but recent research has revealed their importance in many areas of general biology. Our review explores the understudied, yet important mechanisms by which sperm-dependent asexuals that produce non-recombined gametes but rely on their fertilization, can have a significant impact on the evolution of coexisting sexual species and ecosystems. These impacts are concentrated around three major fields. Firstly, sperm-dependent asexuals can potentially impact the gene pool of coexisting sexual species by either restricting their population sizes or by providing bridges for interspecific gene flow whose type and consequences substantially differ from gene flow mechanisms expected under sexual reproduction. Secondly, they may impact on sexuals' diversification rates either directly, by serving as stepping-stones in speciation, or indirectly, by promoting the formation of pre- and postzygotic reproduction barriers among nascent species. Thirdly, they can potentially impact on spatial distribution of species, via direct or indirect (apparent) types of competition and Allee effects. For each such mechanism, we provide empirical examples of how natural sperm-dependent asexuals impact the evolution of their sexual counterparts. In particular, we highlight that these broad effects may last beyond the tenure of the individual asexual lineages causing them, which challenges the traditional perception that asexual lineages are short-lived evolutionary dead ends and minor sideshows. Our review also proposes new research directions to incorporate the aforementioned impacts of sperm-dependent asexuals. These research directions will ultimately enhance our understanding of the evolution of genomes and biological interactions in general.

Quaternary geomorphological and climatic changes associated with the diversification of Iberian freshwater fishes: The case of the genus Cobitis (Cypriniformes, Cobitidae)

We studied the population genetic structure of Cobitis vettonica, an endangered freshwater fish species endemic to the Iberian Peninsula, in order to propose a biogeographic model of the responses of species to the multiple changes that occurred in the Iberian hydrological system during the Quaternary period. We also deciphered the relationship of C. vettonica with its sister species C. paludica, particularly in sympatric areas, and provide genetic information for conservation purposes. To achieve this goal, we analyzed both mitochondrial and nuclear data (the cytochrome b and the nuclear recombination activating 1 genes) and a battery of single‐nucleotide polymorphisms (SNPs) of 248 individuals of C. vettonica or C. paludica from 38 localities, including some sympatric ones, covering the entire distribution area of C. vettonica. We highlight the important role played by the hydrogeomorphological processes and climatic changes that occurred in the Iberian Peninsula during the Quaternary on both the population structure of C. vettonica and its relationship with its sister species C. paludica. Our results support the genetic introgression of populations at the eastern limit of the distribution of C. vettonica. Furthermore, we postulate genetic introgression in sympatric areas. Finally, we propose the establishment or expansion of four Operational Conservation Units (OCUs) for C. vettonica, and highlight the threat faced by its populations due to the low level of genetic diversity detected for some of its populations and genetic introgression with C. paludica, which could eventually displace C. vettonica, resulting in a loss of diversity in this species.

A phylogeny of the genus Limia (Teleostei: Poeciliidae) suggests a single-lake radiation nested in a Caribbean-wide allopatric speciation scenario

OBJECTIVE

The Caribbean is an important global biodiversity hotspot. Adaptive radiations there lead to many speciation events within a limited period and hence are particularly prominent biodiversity generators. A prime example are freshwater fish of the genus Limia, endemic to the Greater Antilles. Within Hispaniola, nine species have been described from a single isolated site, Lake Miragoâne, pointing towards extraordinary sympatric speciation. This study examines the evolutionary history of the Limia species in Lake Miragoâne, relative to their congeners throughout the Caribbean.

RESULTS

For 12 Limia species, we obtained almost complete sequences of the mitochondrial cytochrome b gene, a well-established marker for lower-level taxonomic relationships. We included sequences of six further Limia species from GenBank (total N  = 18 species). Our phylogenies are in concordance with other published phylogenies of Limia. There is strong support that the species found in Lake Miragoâne in Haiti are monophyletic, confirming a recent local radiation. Within Lake Miragoâne, speciation is likely extremely recent, leading to incomplete lineage sorting in the mtDNA. Future studies using multiple unlinked genetic markers are needed to disentangle the relationships within the Lake Miragoâne clade.

Exploring genetic diversity, population structure, and phylogeography in Paracoccidioides species using AFLP markers

Paracoccidioidomycosis (PCM) is a life-threatening systemic fungal infection acquired after inhalation of Paracoccidioides propagules from the environment. The main agents include members of the P. brasiliensis complex (phylogenetically-defined species S1, PS2, PS3, and PS4) and P. lutzii. DNA-sequencing of protein-coding loci (e.g., GP43, ARF, and TUB1) is the reference method for recognizing Paracoccidioides species due to a lack of robust phenotypic markers. Thus, developing new molecular markers that are informative and cost-effective is key to providing quality information to explore genetic diversity within Paracoccidioides. We report using new amplified fragment length polymorphism (AFLP) markers and mating-type analysis for genotyping Paracoccidioides species. The bioinformatic analysis generated 144 in silico AFLP profiles, highlighting two discriminatory primer pairs combinations (#1 EcoRI-AC/MseI-CT and #2 EcoRI-AT/MseI-CT). The combinations #1 and #2 were used in vitro to genotype 165 Paracoccidioides isolates recovered from across a vast area of South America. Considering the overall scored AFLP markers in vitro (67-87 fragments), the values of polymorphism information content (PIC = 0.3345-0.3456), marker index (MI = 0.0018), effective multiplex ratio (E = 44.6788-60.3818), resolving power (Rp = 22.3152-34.3152), discriminating power (D = 0.5183-0.5553), expected heterozygosity (H = 0.4247-0.4443), and mean heterozygosity (H avp  = 0.00002-0.00004), demonstrated the utility of AFLP markers to speciate Paracoccidioides and to dissect both deep and fine-scale genetic structures. Analysis of molecular variance (AMOVA) revealed that the total genetic variance (65-66 %) was due to variability among P. brasiliensis complex and P. lutzii (PhiPT = 0.651-0.658, P < 0.0001), supporting a highly structured population. Heterothallism was the exclusive mating strategy, and the distributions of MAT1-1 or MAT1-2 idiomorphs were not significantly skewed (1:1 ratio) for P. brasiliensis s. str. (χ2 = 1.025; P = 0.3113), P. venezuelensis (χ2 = 0.692; P = 0.4054), and P. lutzii (χ2 = 0.027; P = 0.8694), supporting random mating within each species. In contrast, skewed distributions were found for P. americana (χ2 = 8.909; P = 0.0028) and P. restrepiensis (χ2 = 4.571; P = 0.0325) with a preponderance of MAT1-1. Geographical distributions confirmed that P. americana, P. restrepiensis, and P. lutzii are more widespread than previously thought. P. brasiliensis s. str. is by far the most widely occurring lineage in Latin America countries, occurring in all regions of Brazil. Our new DNA fingerprint assay proved to be rapid, reproducible, and highly discriminatory, to give insights into the taxonomy, ecology, and epidemiology of Paracoccidioides species, guiding disease-control strategies to mitigate PCM.

Molecular phylogeny for the Neotropical freshwater stingrays (Myliobatiformes: Potamotrygoninae) reveals limitations of traditional taxonomy

The subfamily Potamotrygoninae, the only extant clade of elasmobranchs exclusive to freshwater environments, encompasses four genera and 38 species distributed across almost every major South American river basin. Despite their importance in the ornamental fish trade, the taxonomy and evolutionary relationships within potamotrygonines have not yet been resolved. Here, we present a comprehensive molecular phylogeny for the Neotropical freshwater stingrays, based on extensive species and population sampling (35 species and &gt; 350 individuals from drainages across South America). Our phylogeny corroborates the monophyly of the genera Paratrygon and Heliotrygon and the monophyly of the Potamotrygon + Plesiotrygon clade. Within the Potamotrygon + Plesiotrygon clade, we identify a core Potamotrygon clade characterized by short branches, low nodal support and incongruence with current species-level taxonomy. In the core Potamotrygon clade, specimens of widespread species, such as Potamotrygon motoro and Potamotrygon orbignyi, do not form monophyletic lineages; instead, specimens from these species are often closely related to those of other species from the same river basins. These patterns could be caused by inaccurate taxonomy, hybridization, incomplete lineage sorting and rapid diversification. We discuss the conservation of Neotropical freshwater stingrays from a phylogenetic perspective and suggest ways to prioritize potamotrygonid conservation efforts with respect to endemism and evolutionary distinctiveness.

Insight into the Genetic Population Structure of Wild Red Foxes in Poland Reveals Low Risk of Genetic Introgression from Escaped Farm Red Foxes

In this study we assessed the level of genetic introgression between red foxes bred on fur farms in Poland and the native wild population. We also evaluated the impact of a geographic barrier and isolation by distance on gene flow between two isolated subpopulations of the native red fox and their genetic differentiation. Nuclear and mitochondrial DNA was collected from a total of 308 individuals (200 farm and 108 wild red foxes) to study non-native allele flow from farm into wild red fox populations. Genetic structure analyses performed using 24 autosomal microsatellites showed two genetic clusters as being the most probable number of distinct populations. No strong admixture signals between farm and wild red foxes were detected, and significant genetic differentiation was identified between the two groups. This was also apparent from the mtDNA analysis. None of the concatenated haplotypes detected in farm foxes was found in wild animals. The consequence of this was that the haplotype network displayed two genetically distinct groups: farm foxes were completely separated from native ones. Neither the River Vistula nor isolation by distance had a significant impact on gene flow between the separated wild red fox subpopulations. The results of our research indicate a low probability of genetic introgression between farm and native red foxes, and no threat to the genetic integrity of this species.

No evidence for accumulation of deleterious mutations and fitness degradation in clonal fish hybrids: Abandoning sex without regrets

Despite its inherent costs, sexual reproduction is ubiquitous in nature, and the mechanisms to protect it from a competitive displacement by asexuality remain unclear. Popular mutation-based explanations, like the Muller's ratchet and the Kondrashov's hatchet, assume that purifying selection may not halt the accumulation of deleterious mutations in the nonrecombining genomes, ultimately leading to their degeneration. However, empirical evidence is scarce and it remains particularly unclear whether mutational degradation proceeds fast enough to ensure the decay of clonal organisms and to prevent them from outcompeting their sexual counterparts. To test this hypothesis, we jointly analysed the exome sequences and the fitness-related phenotypic traits of the sexually reproducing fish species and their clonal hybrids, whose evolutionary ages ranged from F1 generations to 300 ky. As expected, mutations tended to accumulate in the clonal genomes in a time-dependent manner. However, contrary to the predictions, we found no trend towards increased nonsynonymity of mutations acquired by clones, nor higher radicality of their amino acid substitutions. Moreover, there was no evidence for fitness degeneration in the old clones compared with that in the younger ones. In summary, although an efficacy of purifying selection may still be reduced in the asexual genomes, our data indicate that its efficiency is not drastically decreased. Even the oldest investigated clone was found to be too young to suffer fitness consequences from a mutation accumulation. This suggests that mechanisms other than mutation accumulation may be needed to explain the competitive advantage of sex in the short term.

Can Mitogenomes of the Northern Wheatear (Oenanthe oenanthe) Reconstruct Its Phylogeography and Reveal the Origin of Migrant Birds?

The Northern Wheatear (Oenanthe oenanthe, including the nominate and the two subspecies O. o. leucorhoa and O. o. libanotica) and the Seebohm’s Wheatear (Oenanthe seebohmi) are today regarded as two distinct species. Before, all four taxa were regarded as four subspecies of the Northern Wheatear. Their classification has exclusively been based on ecological and morphological traits, while their molecular characterization is still missing. With this study, we used next-generation sequencing to assemble 117 complete mitochondrial genomes covering O. o. oenanthe, O. o. leucorhoa and O. seebohmi. We compared the resolution power of each individual mitochondrial marker and concatenated marker sets to reconstruct the phylogeny and estimate speciation times of three taxa. Moreover, we tried to identify the origin of migratory wheatears caught on Helgoland (Germany) and on Crete (Greece). Mitogenome analysis revealed two different ancient lineages that separated around 400,000 years ago. Both lineages consisted of a mix of subspecies and species. The phylogenetic trees, as well as haplotype networks are incongruent with the present morphology-based classification. Mitogenome could not distinguish these presumed species. The genetic panmixia among present populations and taxa might be the consequence of mitochondrial introgression between ancient wheatear populations.

Cytogenetic Characterization of Seven Novel satDNA Markers in Two Species of Spined Loaches (Cobitis) and Their Clonal Hybrids

Interspecific hybridization is a powerful evolutionary force. However, the investigation of hybrids requires the application of methodologies that provide efficient and indubitable identification of both parental subgenomes in hybrid individuals. Repetitive DNA, and especially the satellite DNA sequences (satDNA), can rapidly diverge even between closely related species, hence providing a useful tool for cytogenetic investigations of hybrids. Recent progress in whole-genome sequencing (WGS) offers unprecedented possibilities for the development of new tools for species determination, including identification of species-specific satDNA markers. In this study, we focused on spined loaches (Cobitis, Teleostei), a group of fishes with frequent interspecific hybridization. Using the WGS of one species, C. elongatoides, we identified seven satDNA markers, which were mapped by fluorescence in situ hybridization on mitotic and lampbrush chromosomes of C. elongatoides, C. taenia and their triploid hybrids (C. elongatoides × 2C. taenia). Two of these markers were chromosome-specific in both species, one had centromeric localization in multiple chromosomes and four had variable patterns between tested species. Our study provided a novel set of cytogenetic markers for Cobitis species and demonstrated that NGS-based development of satDNA cytogenetic markers may provide a very efficient and easy tool for the investigation of hybrid genomes, cell ploidy, and karyotype evolution.

A molecular assessment of species diversity in Tympanopleura and Ageneiosus catfishes (Auchenipteridae: Siluriformes)

In order to test the congruence of genetic data to the morphologically defined Neotropical catfish genera Tympanopleura and Ageneiosus and explore species diversity, we generated 17 DNA barcodes from five of six species of Tympanopleura and 12 of 13 species of Ageneiosus. To discriminate limits between species, an automatic barcode gap discovery (ABGD), a generalised mixed yule-coalescent model (GYMC) and fixed distance thresholds Kimura two-parameter (K2P; 3%) were used to discriminate putative species limits from the DNA barcodes. The ABGD, GMYC and K2P methods agreed by each generating 13 clusters: six in Tympanopleura (five nominal plus one undescribed species) and seven in Ageneiosus. These clusters corresponded broadly to the described species, except in the case of the Ageneiosus ucayalensis group (A. akamai, A. dentatus, A. intrusus, A. ucayalensis, A. uranophthalmus and A. vittatus). Haplotype sharing and low divergences may have prevented molecular methods from distinguishing these species. We hypothesise that this is the result of a recent radiation of a sympatric species group distributed throughout the Amazon Basin. One putative new species of Tympanopleura was also supported by the molecular data. These results taken together highlight the utility of molecular methods such as DNA barcoding in understanding patterns of diversification across large geographic areas and in recognising overlooked diversity.

OSF-Builder: A New Tool for Constructing and Representing Evolutionary Histories Involving Introgression

Introgression is an evolutionary process which provides an important source of innovation for evolution. Although various methods have been used to detect introgression, very few methods are currently available for constructing evolutionary histories involving introgression. In this article, we propose a new method for constructing such evolutionary histories whose starting point is a species forest (consisting of a collection of lineage trees, usually arising as a collection of clades or monophyletic groups in a species tree), and a gene tree for a specific allele of interest, or allele tree for short. Our method is based on representing introgression in terms of a certain "overlay" of the allele tree over the lineage trees, called an overlaid species forest (OSF). OSFs are similar to phylogenetic networks although a key difference is that they typically have multiple roots because each monophyletic group in the species tree has a different point of origin. Employing a new model for introgression, we derive an efficient algorithm for building OSFs called OSF-Builder that is guaranteed to return an optimal OSF in the sense that the number of potential introgression events is minimized. As well as using simulations to assess the performance of OSF-Builder, we illustrate its use on a butterfly data set in which introgression has been previously inferred. The OSF-Builder software is available for download from https://www.uea.ac.uk/computing/software/OSF-Builder.

Bidirectional mitochondrial introgression between Korean cobitid fish mediated by hybridogenetic hybrids

Genomic introgression through interspecific hybridization has been observed in some species of the freshwater fish family Cobitidae. Within this family, a Cobitis hankugensis-Iksookimia longicorpa diploid-triploid hybrid species complex on the Korean peninsula is unique in displaying hybridogenesis, a unisexual reproduction mode that allows hybrids to mediate the transfer of mitochondrial DNA (but not nuclear DNA) between the two parent species. However, populations of the parental species in the wild have never been examined for the potential effect of introgression on their genomes. To address the genetic consequences of unisexual hybridization on the parental species, we examined genetic structure of the two parental species, C. hankugensis and I. longicorpa, in three independent natural habitats where they coexist with their hybrid complex using DNA sequence data of one mitochondrial gene and three nuclear genes. We found that mitochondrial introgression between the two species was extensive in all the examined localities, while there was no evidence of nuclear introgression across the species boundary. This result indicates that the hybridogenetic individuals mediate mitochondrial introgression from one species to the other, producing mito-nuclear mosaic genomes such as C. hankugensis nuclear genomes associated with I. longicorpa mitochondrial DNA and the reverse. The direction and degree of introgression varied among the three localities, but the underlying mechanisms for this observation proved elusive. Introgression might depend on which species serves as the predominant sperm or ovum donor or the environmental conditions of the localities. The present study suggests that introgressive hybridization between pure C. hankugensis and I. longicorpa species is highly likely where the two species co-occur with hybridogenetic individuals, but the consequence of introgression could be variable due to the history and environmental characteristics of particular populations across the parental species' ranges.

Diversity of spined loaches from Asia Minor in a phylogenetic context (Teleostei: Cobitidae)

Accurate determination of species diversity in areas of high endemicity, particularly those lacking comprehensive systematic knowledge, represents a challenge for both taxonomists and conservationists. This need is particularly evident in areas greatly affected by anthropogenic disturbances such as the Eastern Mediterranean and its freshwater environments. To improve our knowledge of Eastern Mediterranean freshwater fishes, we phylogenetically studied Western Palearctic Cobitis species, focusing on those found in Turkey. Overall, our results provide a robust framework to assess the number of species of Cobitis. Phylogenetic reconstructions based on mitochondrial (cyt b) and nuclear (RAG1) sequences show seven major clades (Clades 1-7) grouping all Western Palearctic Cobitis species, except C. melanoleuca. In general, each major clade comprises Cobitis species that inhabit geographically close areas and have similar secondary sexual characters. Multiple divergent lineages were identified in our analyses, some of which were highly divergent such as the ones inhabiting Turkish freshwaters. Moreover, in some analyses, several of the identified lineages were incongruent with a priori defined species. Furthermore, our analyses identified eight potentially new candidate species, six that had been suggested in previous studies and two that are reported here for the first time. Our results reveal Turkey as the area with the greatest diversity of spined loaches in the Mediterranean.

Modeling Hybridization Under the Network Multispecies Coalescent

Simultaneously modeling hybridization and the multispecies coalescent is becoming increasingly common, and inference of species networks in this context is now implemented in several software packages. This article addresses some of the conceptual issues and decisions to be made in this modeling, including whether or not to use branch lengths and issues with model identifiability. This article is based on a talk given at a Spotlight Session at Evolution 2017 meeting in Portland, Oregon. This session included several talks about modeling hybridization and gene flow in the presence of incomplete lineage sorting. Other talks given at this meeting are also included in this special issue of Systematic Biology.

Plasmodium vivax-like genome sequences shed new insights into Plasmodium vivax biology and evolution

Although Plasmodium vivax is responsible for the majority of malaria infections outside Africa, little is known about its evolution and pathway to humans. Its closest genetic relative, P. vivax-like, was discovered in African great apes and is hypothesized to have given rise to P. vivax in humans. To unravel the evolutionary history and adaptation of P. vivax to different host environments, we generated using long- and short-read sequence technologies 2 new P. vivax-like reference genomes and 9 additional P. vivax-like genotypes. Analyses show that the genomes of P. vivax and P. vivax-like are highly similar and colinear within the core regions. Phylogenetic analyses clearly show that P. vivax-like parasites form a genetically distinct clade from P. vivax. Concerning the relative divergence dating, we show that the evolution of P. vivax in humans did not occur at the same time as the other agents of human malaria, thus suggesting that the transfer of Plasmodium parasites to humans happened several times independently over the history of the Homo genus. We further identify several key genes that exhibit signatures of positive selection exclusively in the human P. vivax parasites. Two of these genes have been identified to also be under positive selection in the other main human malaria agent, P. falciparum, thus suggesting their key role in the evolution of the ability of these parasites to infect humans or their anthropophilic vectors. Finally, we demonstrate that some gene families important for red blood cell (RBC) invasion (a key step of the life cycle of these parasites) have undergone lineage-specific evolution in the human parasite (e.g., reticulocyte-binding proteins [RBPs]).

Genome-specific histories of divergence and introgression between an allopolyploid unisexual salamander lineage and two ancestral sexual species

Quantifying introgression between sexual species and polyploid lineages traditionally thought to be asexual is an important step in understanding what drives the longevity of putatively asexual groups. Here, we capitalize on three recent innovations-ultraconserved element (UCE) sequencing, bioinformatic techniques for identifying genome-specific variation in polyploids, and model-based methods for evaluating historical gene flow-to measure the extent and tempo of introgression over the evolutionary history of an allopolyploid lineage of all-female salamanders and two ancestral sexual species. Our analyses support a scenario in which the genomes sampled in unisexual salamanders last shared a common ancestor with genomes in their parental species ∼3.4 million years ago, followed by a period of divergence between homologous genomes. Recently, secondary introgression has occurred at different times with each sexual species during the last 500,000 years. Sustained introgression of sexual genomes into the unisexual lineage is the defining characteristic of their reproductive mode, but this study provides the first evidence that unisexual genomes have undergone long periods of divergence without introgression. Unlike other sperm-dependent taxa in which introgression is rare, the alternating periods of divergence and introgression between unisexual salamanders and their sexual relatives could explain why these salamanders are among the oldest described unisexual animals.

Why species delimitation matters for fungal ecology: Colletotrichum diversity on wild and cultivated cashew in Brazil

Anthracnose is one of the most important plant diseases globally, occurring on a wide range of cultivated and wild host species. This study aimed to identify the Colletotrichum species associated with cashew anthracnose in Brazil, determine their phylogenetic relationships and geographical distribution, and provide some insight into the factors that may be influencing community composition. Colletotrichum isolates collected from symptomatic leaves, stems, inflorescences, and fruit of cultivated and wild cashew, across four Brazilian biomes, were identified as Colletotrichum chrysophilum, Colletotrichum fragariae, Colletotrichum fructicola, Colletotrichum gloeosporioides sensu stricto, Colletotrichum queenslandicum, Colletotrichum siamense and Colletotrichum tropicale. Colletotrichum siamense was the most dominant species. The greatest species richness was associated with cultivated cashew; leaves harbored more species than the other organs; the Atlantic Forest encompassed more species than the other biomes; and Pernambuco was the most species-rich location. However, accounting for the relative abundance of Colletotrichum species and differences in sample size across strata, the interpretation of which community is most diverse depends on how species are delimited. The present study provides valuable information about the Colletotrichum/cashew pathosystem, sheds light on the causal agents identification,and highlights the impact that species delimitation can have on ecological studies of fungi.

Molecular systematics reveals multiple lineages and cryptic speciation in the freshwater crayfish Parastacus brasiliensis (von Martens, 1869) (Crustacea : Decapoda : Parastacidae)

The characterisation of intraspecific genetic diversity in representatives of the South American crayfish genus Parastacus Huxley, 1879 is here carried out for the first time by comparing populations of Parastacus brasiliensis (von Martens, 1869) as currently defined. Phylogenetic reconstructions based on mitochondrial and nuclear markers indicate the existence of multiple lineages, of which only one can be considered as P. brasiliensis sensu stricto. In addition, there are seven other lineages, one of which is the subspecies Parastacus brasiliensis promatensis Fontoura &amp; Conter, 2008, which is here elevated to species level. We thereby increase to 14 the number of recognised species within Parastacus in South America. Genetic distances among P. brasiliensis sensu stricto and the lineages ‘A’, ‘B’ and ‘C’ increase with geographical distances, suggesting isolation by distance as an important driver of diversification, and eventually speciation, in these burrowing crayfishes. Parastacus brasiliensis occurs mainly in the Guaíba Lake basin and studied populations show limited connectivity and gene flow, probably due to habitat fragmentation. On the basis of genetic and distribution data, the conservation status of P. brasiliensis sensu stricto is confirmed as Near Threatened (NT). These findings lead us to encourage the establishment of preservation areas for isolated populations. The importance of preserving the newly recognised distinct gene pools in order to maintain overall genetic diversity is emphasised.

Population genomics and geographical parthenogenesis in Japanese harvestmen (Opiliones, Sclerosomatidae, Leiobunum)

Naturally occurring population variation in reproductive mode presents an opportunity for researchers to test hypotheses regarding the evolution of sex. Asexual reproduction frequently assumes a geographical pattern, in which parthenogenesis-dominated populations are more broadly dispersed than their sexual conspecifics. We evaluate the geographical distribution of genomic signatures associated with parthenogenesis using nuclear and mitochondrial DNA sequence data from two Japanese harvestman sister taxa, Leiobunum manubriatum and Leiobunum globosum. Asexual reproduction is putatively facultative in these species, and female-biased localities are common in habitat margins. Past karyotypic and current cytometric work indicates L. globosum is entirely tetraploid, while L. manubriatum may be either diploid or tetraploid. We estimated species phylogeny, genetic differentiation, diversity, and mitonuclear discordance in females collected across the species range in order to identify range expansion toward marginal habitat, potential for hybrid origin, and persistence of asexual lineages. Our results point to northward expansion of a tetraploid ancestor of L. manubriatum and L. globosum, coupled with support for greater male gene flow in southern L. manubriatum localities. Specimens from localities in the Tohoku and Hokkaido regions were indistinct, particularly those of L. globosum, potentially due to little mitochondrial differentiation or haplotypic variation. Although L. manubriatum overlaps with L. globosum across its entire range, L. globosum was reconstructed as monophyletic with strong support using mtDNA, and marginal support with nuclear loci. Ultimately, we find evidence for continued sexual reproduction in both species and describe opportunities to clarify the rate and mechanism of parthenogenesis.

Interspecies introgressive hybridization in spiny frogs Quasipaa (Family Dicroglossidae) revealed by analyses on multiple mitochondrial and nuclear genes

Introgression may lead to discordant patterns of variation among loci and traits. For example, previous phylogeographic studies on the genus Quasipaa detected signs of genetic introgression from genetically and morphologically divergent Quasipaa shini or Quasipaa spinosa. In this study, we used mitochondrial and nuclear DNA sequence data to verify the widespread introgressive hybridization in the closely related species of the genus Quasipaa, evaluate the level of genetic diversity, and reveal the formation mechanism of introgressive hybridization. In Longsheng, Guangxi Province, signs of asymmetrical nuclear introgression were detected between Quasipaa boulengeri and Q. shini. Unidirectional mitochondrial introgression was revealed from Q. spinosa to Q. shini. By contrast, bidirectional mitochondrial gene introgression was detected between Q. spinosa and Q. shini in Lushan, Jiangxi Province. Our study also detected ancient hybridizations between a female Q. spinosa and a male Q. jiulongensis in Zhejiang Province. Analyses on mitochondrial and nuclear genes verified three candidate cryptic species in Q. spinosa, and a cryptic species may also exist in Q. boulengeri. However, no evidence of introgressive hybridization was found between Q. spinosa and Q. boulengeri. Quasipaa exilispinosa from all the sampling localities appeared to be deeply divergent from other communities. Our results suggest widespread introgressive hybridization in closely related species of Quasipaa and provide a fundamental basis for illumination of the forming mechanism of introgressive hybridization, classification of species, and biodiversity assessment in Quasipaa.

Hybrid asexuality as a primary postzygotic barrier between nascent species: On the interconnection between asexuality, hybridization and speciation

Although sexual reproduction is ubiquitous throughout nature, the molecular machinery behind it has been repeatedly disrupted during evolution, leading to the emergence of asexual lineages in all eukaryotic phyla. Despite intensive research, little is known about what causes the switch from sexual reproduction to asexuality. Interspecific hybridization is one of the candidate explanations, but the reasons for the apparent association between hybridization and asexuality remain unclear. In this study, we combined cross-breeding experiments with population genetic and phylogenomic approaches to reveal the history of speciation and asexuality evolution in European spined loaches (Cobitis). Contemporary species readily hybridize in hybrid zones, but produce infertile males and fertile but clonally reproducing females that cannot mediate introgressions. However, our analysis of exome data indicates that intensive gene flow between species has occurred in the past. Crossings among species with various genetic distances showed that, while distantly related species produced asexual females and sterile males, closely related species produce sexually reproducing hybrids of both sexes. Our results suggest that hybridization leads to sexual hybrids at the initial stages of speciation, but as the species diverge further, the gradual accumulation of reproductive incompatibilities between species could distort their gametogenesis towards asexuality. Interestingly, comparative analysis of published data revealed that hybrid asexuality generally evolves at lower genetic divergences than hybrid sterility or inviability. Given that hybrid asexuality effectively restricts gene flow, it may establish a primary reproductive barrier earlier during diversification than other "classical" forms of postzygotic incompatibilities. Hybrid asexuality may thus indirectly contribute to the speciation process.

Characterization of hybridization within a secondary contact region of the inshore fish, Bostrychus sinensis, in the East China Sea

The northwest Pacific marginal seas are a primary center of phylogeographic and evolutionary research, because of their dynamic geographic history of falling and rising sea levels during the glaciations and interglaciations of the last one million years. Here we present new molecular and morphological data for geographic samples of the four-eyed sleeper (Bostrychus sinensis), which reinforce the evidence for secondary contact and hybridization between two phylogeographic lineages in the East China Sea. Specifically, we find that the secondary contact region is characterized by a low frequency of hybridization, where mitochondrial DNA introgression is relatively common, whereas F₁ hybrids are correspondingly scarce. Furthermore, the adult standard lengths of the two phylogeographic lineages vary geographically in a manner that is consistent with reproductive character displacement. Collectively, the molecular and morphological data document that sleeper hybridization conforms to the classic "tension zone" model, where alleles are lost via reduced hybrid viability and/or positive assortative mating but are then replenished by dispersal from south of the secondary contact region. They also indicate that the two phylogeographic lineages are at an incipient stage of the speciation process. These results and conclusions for the four-eyed sleeper are presented as a case study for future research on the vicariance, secondary contact, and hybridization of marine groups in the northwest Pacific marginal seas.

Origins of the Greenland shark (Somniosus microcephalus): Impacts of ice-olation and introgression

Herein, we use genetic data from 277 sleeper sharks to perform coalescent‐based modeling to test the hypothesis of early Quaternary emergence of the Greenland shark (Somniosus microcephalus) from ancestral sleeper sharks in the Canadian Arctic‐Subarctic region. Our results show that morphologically cryptic somniosids S. microcephalus and Somniosus pacificus can be genetically distinguished using combined mitochondrial and nuclear DNA markers. Our data confirm the presence of genetically admixed individuals in the Canadian Arctic and sub‐Arctic, and temperate Eastern Atlantic regions, suggesting introgressive hybridization upon secondary contact following the initial species divergence. Conservative substitution rates fitted to an Isolation with Migration (IM) model indicate a likely species divergence time of 2.34 Ma, using the mitochondrial sequence DNA, which in conjunction with the geographic distribution of admixtures and Pacific signatures likely indicates speciation associated with processes other than the closing of the Isthmus of Panama. This time span coincides with further planetary cooling in the early Quaternary period followed by the onset of oscillating glacial‐interglacial cycles. We propose that the initial S. microcephalus–S. pacificus split, and subsequent hybridization events, were likely associated with the onset of Pleistocene glacial oscillations, whereby fluctuating sea levels constrained connectivity among Arctic oceanic basins, Arctic marginal seas, and the North Atlantic Ocean. Our data demonstrates support for the evolutionary consequences of oscillatory vicariance via transient oceanic isolation with subsequent secondary contact associated with fluctuating sea levels throughout the Quaternary period—which may serve as a model for the origins of Arctic marine fauna on a broad taxonomic scale.

Genome-scale data reveal that endemic Poecilia populations from small sulphidic springs display no evidence of inbreeding

Populations with limited ranges can be highly vulnerable to changes in their environment and are, thus, of high conservation concern. Populations that experience human-induced range reductions are often highly inbred and lack genetic diversity, but it is unknown whether this is also the case for populations with naturally small ranges. The fishes Poecilia sulphuraria (listed as critically endangered) and Poecilia thermalis, which are endemic to small hydrogen sulphide-rich springs in southern Mexico, are examples of such populations with inherently small habitats. We used geometric morphometrics and population genetics to quantify phenotypic and genetic variation within and among two populations of P. sulphuraria and one population of P. thermalis. Principal component analyses revealed phenotypic and genetic differences among the populations. Evidence for inbreeding was low compared to populations that have undergone habitat reduction. The genetic data were also used to infer the demographic history of these populations to obtain estimates for effective population sizes and migration rates. Effective population sizes were large given the small habitats of these populations. Our results imply that these three endemic extremophile populations should each be considered separately for conservation purposes. Additionally, this study suggests that populations in naturally small habitats may have lower rates of inbreeding and higher genetic diversity than expected, and therefore may be better equipped to handle environmental perturbations than anticipated. We caution, however, that the inferred lack of inbreeding and the large effective population sizes could potentially be a result of colonization by genetically diverse ancestors.

New insight into hybridization and unidirectional introgression between Ammodytes japonicus and Ammodytes heian (Trachiniformes, Ammodytidae)

Based on northern (NOL) and southern (SOL) mitochondrial lineages, recently, it proposed the new species Ammodytes heian and revived the species name Ammodytes japonicus to describe sand lances from the northwestern Pacific Ocean. This study used molecular methods to investigate genetic relationships between the two sand lance species in Korea and Japan. In total, 154 specimens were collected from four locations in Korea (Baengnyeongdo in the Yellow Sea, Tongyeong in the Korean Strait, and Jumunjin and Gijang in the East Sea), and 50 specimens were collected from a single location in Japan (Wakkanai in the Okhotsk Sea). Mitochondrial DNA analysis demonstrated that the individuals from Baengnyeongdo and Tongyeong all belonged to the SOL, whereas those from Gijang, Jumunjin, and Wakkanai included individuals from both the NOL and SOL (over 75% NOL). Population structure analyses were performed on the same individuals using seven microsatellite DNA markers. The population structure analysis based on 201 specimens identified two clusters (named as northern group and southern group), with the admixture proportion (q) of < 0.1 for the northern group in the Backyeongdo and Tongyeong sand lances and < 0.1 for the southern group in the Wakkanai sand lances. The high heterogeneity indicated that the former was probably A. japonicus and the latter probably A. heian. However, the admixture proportion in the Jumunjin and Gijang sand lances was 0.71–0.75 for the southern group, indicating that hybridization and unidirectional introgression from SOL to NOL occurs in southwestern margin of the East Sea. Our findings illustrate the speciation process based on different patterns of gene flow between Korean and Japanese sand lance, which is strongly influenced by both the paleo-climatic change and the contemporary local oceanic current pattern.

Terrestrial mountain islands and Pleistocene climate fluctuations as motors for speciation: A case study on the genus Pseudovelia (Hemiptera: Veliidae)

This study investigated the influences of geographic isolation and climate fluctuation on the genetic diversity, speciation, and biogeography of the genus Pseudovelia (Hemiptera: Veliidae) in subtropical China and tropic Indo-China Peninsula. Species nucleotide and haplotype diversities decreased with reduction in species distribution limits. The gene tree was congruent with the taxonomy of monophyly, except for four species, P. contorta, P. extensa, P. tibialis tibialis, and P. vittiformis. The conflicts between the genes and species tree could be due to long-term isolation and incomplete lineage sorting. Diversification analysis showed that the diversification rate (0.08 sp/My shifted to 0.5 sp/My) changed at 2.1 Ma, which occurred in the early Pleistocene period. Ancestral area reconstruction suggested that subtropical species possibly evolved from the tropics region (i.e., Indo-China Peninsula). Results implied that narrow endemics harbored relatively low genetic diversity because of small effective population and genetic drift. Radiation of subtropical Pseudovelia species was rapidly promoted by Pleistocene climate fluctuations and geographic isolation. The acute rising of the Hengduan Mountain with the entire uplift of the Qinghai-Tibet Plateau induced the initial differentiation of Pseudovelia species. These results highlighted the importance of geographical isolation and climate changes in promoting speciation in mountain habitat islands.

Marked hybridization and introgression in Ophrys sect. Pseudophrys in the western Iberian Peninsula

PREMISE OF THE STUDY

Orchids in the genus Ophrys represent extraordinary cases of tight coevolution between plants and their pollinators, and as a result, they present opportunities for studying hybridization, or a lack thereof, during speciation. However, few studies assess the real effect of hybridization in diversification. The three most representative species of section Pseudophrys in the western Iberian Peninsula-O. dyris, O. fusca, and O. lutea-were chosen to study evolutionary relationships and examine speciation.

METHODS

Using eight specific nuclear microsatellite loci, 357 individuals from 28 locations were studied; 142 of these samples were also studied with four plastid microsatellite loci. Data were analyzed using Bayesian cluster analysis, a median-joint network, and multivariate analysis.

KEY RESULTS

Many O. dyris and O. fusca specimens had three or four alleles and were therefore treated as tetraploid. Ophrys dyris is poorly genetically separated from O. fusca, and pure populations are rare. Ophrys fusca and O. lutea are distinct, but hybrids/introgressed individuals were detected in most of the populations and supported by plastid haplotypes. Ophrys fusca is subdivided into three well-delimited genetic lineages with a strict geographic correspondence confirmed by plastid haplotypes.

CONCLUSIONS

Because postzygotic barriers are weak, leakage in this highly specialized orchid-pollinator system contributes to hybridization and introgression. These leakages may have occurred during periods of past climate change, promoting homogenization and the potential for generations of new biodiversity via production of novel genotypes/phenotypes interacting with pollinators.

Differential introgression suggests candidate beneficial and barrier loci between two parapatric subspecies of Pearson's horseshoe bat Rhinolophus pearsoni

Observations that rates of introgression between taxa can vary across loci are increasingly common. Here, we test for differential locus-wise introgression in 2 parapatric subspecies of Pearson′s horseshoe bat (Rhinolophus pearsoni chinensis and R. p. pearsoni). To efficiently identify putative speciation genes and/or beneficial genes in our current system, we used a candidate gene approach by including loci from X chromosome that are suggested to be more likely involved in reproductive isolation in other organisms and loci underlying hearing that have been suggested to spread across the hybrid zone in another congeneric species. Phylogenetic and coalescent analyses were performed at 2 X-linked, 4 hearing genes, as well as 2 other autosomal loci individually. Likelihood ratio tests could not reject the model of zero gene flow at 2 X-linked and 2 autosomal genes. In contrast, gene flow was supported at 3 of 4 hearing genes. While this introgression could be adaptive, we cannot rule out stochastic processes. Our results highlight the utility of the candidate gene approach in searching for speciation genes and/or beneficial genes across the species boundary in natural populations.

Asexual Reproduction Does Not Apparently Increase the Rate of Chromosomal Evolution: Karyotype Stability in Diploid and Triploid Clonal Hybrid Fish (Cobitis, Cypriniformes, Teleostei)

Interspecific hybridization, polyploidization and transitions from sexuality to asexuality considerably affect organismal genomes. Especially the last mentioned process has been assumed to play a significant role in the initiation of chromosomal rearrangements, causing increased rates of karyotype evolution. We used cytogenetic analysis and molecular dating of cladogenetic events to compare the rate of changes of chromosome morphology and karyotype in asexually and sexually reproducing counterparts in European spined loach fish (Cobitis). We studied metaphases of three sexually reproducing species and their diploid and polyploid hybrid clones of different age of origin. The material includes artificial F1 hybrid strains, representatives of lineage originated in Holocene epoch, and also individuals of an oldest known age to date (roughly 0.37 MYA). Thereafter we applied GISH technique as a marker to differentiate parental chromosomal sets in hybrids. Although the sexual species accumulated remarkable chromosomal rearrangements after their speciation, we observed no differences in chromosome numbers and/or morphology among karyotypes of asexual hybrids. These hybrids possess chromosome sets originating from respective parental species with no cytogenetically detectable recombinations, suggesting their integrity even in a long term. The switch to asexual reproduction thus did not provoke any significant acceleration of the rate of chromosomal evolution in Cobitis. Asexual animals described in other case studies reproduce ameiotically, while Cobitis hybrids described here produce eggs likely through modified meiosis. Therefore, our findings indicate that the effect of asexuality on the rate of chromosomal change may be context-dependent rather than universal and related to particular type of asexual reproduction.

Untangling reticulate evolutionary relationships among New World and Hawaiian mints (Stachydeae, Lamiaceae)

The phenomenon of polyploidy and hybridization usually results in novel genetic combinations, leading to complex, reticulate evolution and incongruence among gene trees, which in turn may show different phylogenetic histories than the inherent species tree. The largest tribe within the subfamily Lamioideae (Lamiaceae), Stachydeae, which includes the globally distributed Stachys, and one of the largest Hawaiian angiosperm radiations, the endemic mints, is a widespread and taxonomically challenging lineage displaying a wide spectrum of morphological and chromosomal diversity. Previous molecular phylogenetic studies have showed that while the Hawaiian mints group with Mexican-South American Stachys based on chloroplast DNA sequence data, nuclear ribosomal DNA (nrDNA) sequences suggest that they are most closely related to temperate North American Stachys. Here, we have utilized five independently inherited, low-copy nuclear loci, and a variety of phylogenetic methods, including multi-locus coalescence-based tree reconstructions, to provide insight into the complex origins and evolutionary relationships between the New World Stachys and the Hawaiian mints. Our results demonstrate incongruence between individual gene trees, grouping the Hawaiian mints with both temperate North American and Meso-South American Stachys clades. However, our multi-locus coalescence tree is concurrent with previous nrDNA results placing them within the temperate North American Stachys clade. Our results point toward a possible allopolyploid hybrid origin of the Hawaiian mints arising from temperate North American and Meso-South American ancestors, as well as a reticulate origin for South American Stachys. As such, our study is another significant step toward further understanding the putative parentage and the potential influence of hybridization and incomplete lineage sorting in giving rise to this insular plant lineage, which following colonization underwent rapid morphological and ecological diversification.

Phylogenetics and phylogeography of a long-legged harvestman (Arachnida : Opiliones) in the Brazilian Atlantic Rain Forest reveals poor dispersal, low diversity and extensive mitochondrial introgression

We used DNA sequence data to test the morphology-based taxonomy and examine the biogeography of the Brazilian Atlantic Rain Forest genus Promitobates. Most species are well differentiated morphologically, and a previous morphological phylogeny recovered the genus as monophyletic. However, some of these species have overlapping geographical distributions and considerable intraspecific variation, perhaps representing a species complex. Mitochondrial (12S rRNA and COI) and nuclear (ITS2 and 28S) genes were sequenced from 132 specimens collected from 27 localities. The results are consistent with significant mitochondrial introgression among the species P. ornatus, P. hatschbachi, P. lager, P. bellus and P. intermedius (the ‘P. ornatus species complex’), with one specimen identified as a hybrid between P. nigripes and this complex. A phylogeographic study of the complex was conducted using mitochondrial haplotypes. This revealed remarkably poor dispersal among populations, with only one case of a shared haplotype, and very low genetic diversity. The phylogeny showed a clear break between populations on either side of a narrow region of forest, suggesting an important historical event separated these lineages in the genus. The analyses also pointed to population breaks that date back several millions of years or extremely small effective population sizes, depending on the mutation rate.