Speciation of pelagic zooplankton: Invisible boundaries can drive isolation of oceanic ctenophores

The study of evolution and speciation in non-model systems provides us with an opportunity to expand our understanding of biodiversity in nature. Connectivity studies generally focus on species with obvious boundaries to gene flow, but in open-ocean environments, such boundaries are difficult to identify. Due to the lack of obvious boundaries, speciation and population subdivision in the pelagic environment remain largely unexplained. Comb jellies (Phylum Ctenophora) are mostly planktonic gelatinous invertebrates, many of which are considered to have freely interbreeding distributions worldwide. It is thought that the lobate ctenophore Bolinopsis infundibulum is distributed throughout cooler northern latitudes and B. vitrea warmer. Here, we examined the global population structure for species of Bolinopsis with genetic and morphological data. We found distinct evolutionary patterns within the genus, where B. infundibulum had a broad distribution from northern Pacific to Atlantic waters despite many physical barriers, while other species were geographically segregated despite few barriers. Divergent patterns of speciation within the genus suggest that oceanic currents, sea-level, and geological changes over time can act as either barriers or aids to dispersal in the pelagic environment. Further, we used population genomic data to examine evolution in the open ocean of a distinct lineage of Bolinopsis ctenophores from the North Eastern Pacific. Genetic information and morphological observations validated this as a separate species, Bolinopsis microptera, which was previously described but has recently been called B. infundibulum. We found that populations of B. microptera from California were in cytonuclear discordance, which indicates a secondary contact zone for previously isolated populations. Discordance at this scale is rare, especially in a continuous setting.

How mitonuclear discordance and geographic variation have confounded species boundaries in a widely studied snake

As DNA sequencing technologies and methods for delimiting species with genomic data become more accessible and numerous, researchers have more tools than ever to investigate questions in systematics and phylogeography. However, easy access to sophisticated computational tools is not without its drawbacks. Choosing the right approach for one's question can be challenging when presented with multitudinous options, some of which fail to distinguish between species and intraspecific population structure. Here, we employ a methodology that emphasizes intensive geographic sampling, particularly at contact zones between populations, with a focus on differentiating intraspecific genetic clusters from species in the Pantherophis guttatus complex, a group of North American ratsnakes. Using a mitochondrial marker as well as ddRADseq data, we find evidence of mitonuclear discordance which has contributed to historical confusion about the relationships within this group. Additionally, we identify geographically and genetically structured populations within the species Pantherophis emoryi that are congruent with previously described morphological variation. Importantly, we find that these structured populations within P. emoryi are highly admixed throughout the range of the species and show no evidence of any reproductive isolation. Our data support a revision of the taxonomy of this group, and we recognize two species within the complex and three subspecies within P. emoryi. This study illustrates the importance of thorough sampling of contact zones and consideration of gene flow when delimiting species in widespread complexes containing parapatric lineages.

The conservation value of admixed phenotypes in a critically endangered species complex

In today's environmental crisis, conservationists are increasingly confronted with terminally endangered species whose last few surviving populations may be affected by allelic introgression from closely related species. Yet there is a worrying lack of evidence-based recommendations and solutions for this emerging problem. We analyzed genome-wide DNA markers and plumage variability in a critically endangered insular songbird, the Black-winged Myna (BWM, Acridotheres melanopterus). This species is highly threatened by the illegal wildlife trade, with its wild population numbering in the low hundreds, and its continued survival urgently depending on ex-situ breeding. Its three subspecies occur along a geographic gradient of melanism and are variably interpreted as three species. However, our integrative approach revealed that melanism poorly reflects the pattern of limited genomic differentiation across BWM subspecies. We also uncovered allelic introgression into the most melanistic subspecies, tertius, from the all-black congeneric Javan Myna (A. javanicus), which is native to the same islands. Based on our results, we recommend the establishment of three separate breeding programs to maintain subspecific traits that may confer local adaptation, but with the option of occasional cross-breeding between insurance populations in order to boost genetic diversity and increase overall viability prospects of each breeding program. Our results underscore the importance of evidence-based integrative approaches when determining appropriate conservation units. Given the rapid increase of terminally endangered organisms in need of ex-situ conservation, this study provides an important blueprint for similar programs dealing with phenotypically variable species.

Genome-wide data reveal discordant mitonuclear introgression in the intermediate horseshoe bat (Rhinolophus affinis)

Closely related taxa often exhibit mitonuclear discordance attributed to introgression of mitochondrial DNA (mtDNA), yet few studies have considered the underlying causes of mtDNA introgression. Here we test for demographic versus adaptive processes as explanations for mtDNA introgression in three subspecies of the intermediate horseshoe bat (Rhinolophus affinis). We generated sequences of 1692 nuclear genes and 13 mitochondrial protein-coding genes for 48 individuals. Phylogenetic reconstructions based on 320 exon sequences and 2217 single nucleotide polymorphisms (SNPs) both revealed conflicts between the species tree and mtDNA tree. These results, together with geographic patterns of mitonuclear discordance, and shared identical or near-identical mtDNA sequences, suggest extensive introgression of mtDNA between the two parapatric mainland subspecies. Under demographic hypotheses, we would also expect to uncover traces of ncDNA introgression, however, population structure and gene flow analyses revealed little nuclear admixture. Furthermore, we found inconsistent estimates of the timing of population expansion and that of the most recent common ancestor for the clade containing introgressed haplotypes. Without a clear demographic explanation, we also examined whether introgression likely arises from adaptation. We found that five mtDNA genes contained fixed amino acid differences between introgressed and non-introgressed individuals, including putative positive selection found in one codon, although this did not show introgression. While our evidence for rejecting demographic hypotheses is arguably stronger than that for rejecting adaptation, we find no definitive support for either explanation. Future efforts will focus on larger-scale resequencing to decipher the underlying causes of discordant mitonuclear introgression in this system.

Phylogeographic and phenotypic outcomes of brown anole colonization across the Caribbean provide insight into the beginning stages of an adaptive radiation

Some of the most important insights into the ecological and evolutionary processes of diversification and speciation have come from studies of island adaptive radiations, yet relatively little research has examined how these radiations initiate. We suggest that Anolis sagrei is a candidate for understanding the origins of the Caribbean Anolis adaptive radiation and how a colonizing anole species begins to undergo allopatric diversification, phenotypic divergence and, potentially, speciation. We undertook a genomic and morphological analysis of representative populations across the entire native range of A. sagrei, finding that the species originated in the early Pliocene, with the deepest divergence occurring between western and eastern Cuba. Lineages from these two regions subsequently colonized the northern Caribbean. We find that at the broadest scale, populations colonizing areas with fewer closely related competitors tend to evolve larger body size and more lamellae on their toepads. This trend follows expectations for post-colonization divergence from progenitors and convergence in allopatry, whereby populations freed from competition with close relatives evolve towards common morphological and ecological optima. Taken together, our results show a complex history of ancient and recent Cuban diaspora with populations on competitor-poor islands evolving away from their ancestral Cuban populations regardless of their phylogenetic relationships, thus providing insight into the original diversification of colonist anoles at the beginning of the radiation. Our research also supplies an evolutionary framework for the many studies of this increasingly important species in ecological and evolutionary research.

A continent-scale study of the social structure and phylogeography of the bent-wing bat, Miniopterus schreibersii (Mammalia: Chiroptera), using new microsatellite data

Miniopterus schreibersii is a cave-dwelling bat species with a wide distribution in the western Palearctic spanning southern and central Europe, North Africa, Anatolia, the Middle East, and the Caucasus. We investigated the social structure and its effects on the genetic makeup of this species, using 10 nuclear microsatellite markers and a partial fragment of the mitochondrial cytochrome b gene. Samples were examined from the species' entire circum-Mediterranean range. Local structuring that was previously detected among populations of M. schreibersii using mitochondrial markers was not observed for microsatellite markers, indicating male-biased dispersal for the species. Some support was found for postglacial expansions in Europe, with Anatolia potentially acting as the primary refugium during the Last Glacial Maximum (LGM). However, support for this hypothesis is not as strong as that previously detected using mitochondrial DNA markers. This is likely due to the diminishing effect of male-mediated dispersal, replenishing the nuclear diversity faster than the mitochondrial diversity in regions that are relatively far from the glacial refugia.

Mitochondrial Ecophysiology: Assessing the Evolutionary Forces That Shape Mitochondrial Variation

The mitonuclear species concept hypothesizes that incompatibilities between interacting gene products of the nuclear and mitochondrial genomes are a major factor establishing and maintaining species boundaries. However, most of the data available to test this concept come from studies of genetic variation in mitochondrial DNA, and clines in the mitochondrial genome across contact zones can be produced by a variety of forces. Here, we show that using a combination of population genomic analyses of the nuclear and mitochondrial genomes and studies of mitochondrial function can provide insight into the relative roles of neutral processes, adaptive evolution, and mitonuclear incompatibility in establishing and maintaining mitochondrial clines, using Atlantic killifish (Fundulus heteroclitus) as a case study. There is strong evidence for a role of secondary contact following the last glaciation in shaping a steep mitochondrial cline across a contact zone between northern and southern subspecies of killifish, but there is also evidence for a role of adaptive evolution in driving differentiation between the subspecies in a variety of traits from the level of the whole organism to the level of mitochondrial function. In addition, studies are beginning to address the potential for mitonuclear incompatibilities in admixed populations. However, population genomic studies have failed to detect evidence for a strong and pervasive influence of mitonuclear incompatibilities, and we suggest that polygenic selection may be responsible for the complex patterns observed. This case study demonstrates that multiple forces can act together in shaping mitochondrial clines, and illustrates the challenge of disentangling their relative roles.

'Ghost introgression' as a cause of deep mitochondrial divergence in a bird species complex

In the absence of nuclear-genomic differentiation between two populations, deep mitochondrial divergence (DMD) is a form of mito-nuclear discordance. Such instances of DMD are rare and might variably be explained by unusual cases of female-linked selection, by male-biased dispersal, by 'speciation reversal' or by mitochondrial capture through genetic introgression. Here we analyze DMD in an Asian Phylloscopus leaf warbler (Aves: Phylloscopidae) complex. Bioacoustic, morphological and genomic data demonstrate close similarity between the taxa affinis and occisinensis, even though DMD previously led to their classification as two distinct species. Using population genomic and comparative genomic methods on 45 whole genomes, including historical reconstructions of effective population size, genomic peaks of differentiation and genomic linkage, we infer that the form affinis is likely the product of a westward expansion in which it replaced a now-extinct congener that was the donor of its mtDNA and small portions of its nuclear genome. This study provides strong evidence of 'ghost introgression' as the cause of DMD, and we suggest that 'ghost introgression' may be a widely overlooked phenomenon in nature.

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

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

How repeatable is microevolution on islands? Patterns of dispersal and colonization-related plant traits in a phylogeographical context

BACKGROUND AND AIMS

Archipelagos provide a valuable framework for investigating phenotypic evolution under different levels of geographical isolation. Here, we analysed two co-distributed, widespread plant lineages to examine if incipient island differentiation follows parallel patterns of variation in traits related to dispersal and colonization.

METHODS

Twenty-one populations of two anemochorous Canarian endemics, Kleinia neriifolia and Periploca laevigata, were sampled to represent mainland congeners and two contrasting exposures across all the main islands. Leaf size, seed size and dispersability (estimated as diaspore terminal velocity) were characterized in each population. For comparison, dispersability was also measured in four additional anemochorous island species. Plastid DNA data were used to infer genetic structure and to reconstruct the phylogeographical pattern of our focal species.

KEY RESULTS

In both lineages, mainland-island phenotypic divergence probably started within a similar time frame (i.e. Plio-Pleistocene). Island colonization implied parallel increases in leaf size and dispersability, but seed size showed opposite patterns of variation between Kleinia and Periploca species pairs. Furthermore, dispersability in our focal species was low when compared with other island plants, mostly due to large diaspore sizes. At the archipelago scale, island exposure explained a significant variation in leaf size across islands, but not in dispersability or seed size. Combined analyses of genetic and phenotypic data revealed two consistent patterns: (1) extensive within-island but very limited among-island dispersal, and (2) recurrent phenotypic differentiation between older (central) and younger (peripheral) island populations.

CONCLUSIONS

Leaf size follows a more predictable pattern than dispersability, which is affected by stochastic shifts in seed size. Increased dispersability is associated with high population connectivity at the island scale, but does not preclude allopatric divergence among islands. In sum, phenotypic convergent patterns between species suggest a major role of selection, but deviating traits also indicate the potential contribution of random processes, particularly on peripheral islands.

An evaluation of alternative explanations for widespread cytonuclear discordance in annual sunflowers (Helianthus)

Cytonuclear discordance is commonly observed in phylogenetic studies, yet few studies have tested whether these patterns reflect incomplete lineage sorting or organellar introgression. Here, we used whole-chloroplast sequence data in combination with over 1000 nuclear single-nucleotide polymorphisms to clarify the extent of cytonuclear discordance in wild annual sunflowers (Helianthus), and to test alternative explanations for such discordance. Our phylogenetic analyses indicate that cytonuclear discordance is widespread within this group, both in terms of the relationships among species and among individuals within species. Simulations of chloroplast evolution show that incomplete lineage sorting cannot explain these patterns in most cases. Instead, most of the observed discordance is better explained by cytoplasmic introgression. Molecular tests of evolution further indicate that selection may have played a role in driving patterns of plastid variation - although additional experimental work is needed to fully evaluate the importance of selection on organellar variants in different parts of the geographic range. Overall, this study represents one of the most comprehensive tests of the drivers of cytonuclear discordance and highlights the potential for gene flow to lead to extensive organellar introgression in hybridizing taxa.

Testing Hypotheses of Diversification in Panamanian Frogs and Freshwater Fishes Using Hierarchical Approximate Bayesian Computation with Model Averaging

Most Neotropical frog and freshwater fish species sampled to date show phylogeographic breaks along the Pacific coast of the Isthmus of Panama, with lineages in Costa Rica and western Panama isolated from central Panama. We examine temporal patterns of diversification of taxa across this ‘western Panama isthmus’ (WPI) break to test hypotheses about the origin of species geographical distributions and genetic structuring in this region. We tested for synchronous diversification of four codistributed frog taxon-pairs and three fish taxon-pairs sharing the WPI break using hierarchical approximate Bayesian computation with model averaging based on mitochondrial DNA sequences. We also estimated lineage divergence times using full-Bayesian models. Several of our results supported synchronous divergences within the frog and freshwater fish assemblages; however, Bayes factor support was equivocal for or against synchronous or asynchronous diversification. Nevertheless, we infer that frog populations were likely isolated by one or multiple Pliocene–Pleistocene events more recently than predicted by previous models, while fish genetic diversity was structured by Pleistocene events. By integrating our results with external information from geology and elevational sea level modeling, we discuss the implications of our findings for understanding the biogeographical scenario of the diversification of Panamanian frogs and fishes. Consistent with the ‘Bermingham/Martin model’ (Molecular Ecology 1998, 7, 499–517), we conclude that the regional fish assemblage was fractured by processes shaping isthmian landscapes during the Pleistocene glaciations, including drainage basin isolation during lowered sea levels.

Selection and environmental adaptation along a path to speciation in the Tibetan frog Nanorana parkeri

Tibetan frogs, Nanorana parkeri, are differentiated genetically but not morphologically along geographical and elevational gradients in a challenging environment, presenting a unique opportunity to investigate processes leading to speciation. Analyses of whole genomes of 63 frogs reveal population structuring and historical demography, characterized by highly restricted gene flow in a narrow geographic zone lying between matrilines West (W) and East (E). A population found only along a single tributary of the Yalu Zangbu River has the mitogenome only of E, whereas nuclear genes of W comprise 89-95% of the nuclear genome. Selection accounts for 579 broadly scattered, highly divergent regions (HDRs) of the genome, which involve 365 genes. These genes fall into 51 gene ontology (GO) functional classes, 14 of which are likely to be important in driving reproductive isolation. GO enrichment analyses of E reveal many overrepresented functional categories associated with adaptation to high elevations, including blood circulation, response to hypoxia, and UV radiation. Four genes, including DNAJC8 in the brain, TNNC1 and ADORA1 in the heart, and LAMB3 in the lung, differ in levels of expression between low- and high-elevation populations. High-altitude adaptation plays an important role in maintaining and driving continuing divergence and reproductive isolation. Use of total genomes enabled recognition of selection and adaptation in and between populations, as well as documentation of evolution along a stepped cline toward speciation.

Sex chromosomes and speciation in birds and other ZW systems

Theory and empirical patterns suggest a disproportionate role for sex chromosomes in evolution and speciation. Focusing on ZW sex determination (females ZW, males ZZ; the system in birds, many snakes, and lepidopterans), I review how evolutionary dynamics are expected to differ between the Z, W and the autosomes, discuss how these differences may lead to a greater role of the sex chromosomes in speciation and use data from birds to compare relative evolutionary rates of sex chromosomes and autosomes. Neutral mutations, partially or completely recessive beneficial mutations, and deleterious mutations under many conditions are expected to accumulate faster on the Z than on autosomes. Sexually antagonistic polymorphisms are expected to arise on the Z, raising the possibility of the spread of preference alleles. The faster accumulation of many types of mutations and the potential for complex evolutionary dynamics of sexually antagonistic traits and preferences contribute to a role for the Z chromosome in speciation. A quantitative comparison among a wide variety of bird species shows that the Z tends to have less within-population diversity and greater between-species differentiation than the autosomes, likely due to both adaptive evolution and a greater rate of fixation of deleterious alleles. The W chromosome also shows strong potential to be involved in speciation, in part because of its co-inheritance with the mitochondrial genome. While theory and empirical evidence suggest a disproportionate role for sex chromosomes in speciation, the importance of sex chromosomes is moderated by their small size compared to the whole genome.

Our love-hate relationship with DNA barcodes, the Y2K problem, and the search for next generation barcodes

DNA barcodes are very useful for species identification especially when identification by traditional morphological characters is difficult. However, the short mitochondrial and chloroplast barcodes currently in use often fail to distinguish between closely related species, are prone to lateral transfer, and provide inadequate phylogenetic resolution, particularly at deeper nodes. The deficiencies of short barcode identifiers are similar to the deficiencies of the short year identifiers that caused the Y2K problem in computer science. The resolution of the Y2K problem was to increase the size of the year identifiers. The performance of conventional mitochondrial COI barcodes for phylogenetics was compared with the performance of complete mitochondrial genomes and nuclear ribosomal RNA repeats obtained by genome skimming for a set of caddisfly taxa (Insect Order Trichoptera). The analysis focused on Trichoptera Family Hydropsychidae, the net-spinning caddisflies, which demonstrates many of the frustrating limitations of current barcodes. To conduct phylogenetic comparisons, complete mitochondrial genomes (15 kb each) and nuclear ribosomal repeats (9 kb each) from six caddisfly species were sequenced, assembled, and are reported for the first time. These sequences were analyzed in comparison with eight previously published trichopteran mitochondrial genomes and two triochopteran rRNA repeats, plus outgroup sequences from sister clade Lepidoptera (butterflies and moths). COI trees were not well-resolved, had low bootstrap support, and differed in topology from prior phylogenetic analyses of the Trichoptera. Phylogenetic trees based on mitochondrial genomes or rRNA repeats were well-resolved with high bootstrap support and were largely congruent with each other. Because they are easily sequenced by genome skimming, provide robust phylogenetic resolution at various phylogenetic depths, can better distinguish between closely related species, and (in the case of mitochondrial genomes), are backwards compatible with existing mitochondrial barcodes, it is proposed that mitochondrial genomes and rRNA repeats be used as next generation DNA barcodes.

Divergent clades or cryptic species? Mito-nuclear discordance in a Daphnia species complex

BACKGROUND

Genetically divergent cryptic species are frequently detected by molecular methods. These discoveries are often a byproduct of molecular barcoding studies in which fragments of a selected marker are used for species identification. Highly divergent mitochondrial lineages and putative cryptic species are even detected in intensively studied animal taxa, such as the crustacean genus Daphnia. Recently, eleven such lineages, exhibiting genetic distances comparable to levels observed among well-defined species, were recorded in the D. longispina species complex, a group that contains several key taxa of freshwater ecosystems. We tested if three of those lineages represent indeed distinct species, by analyzing patterns of variation of ten nuclear microsatellite markers in six populations.

RESULTS

We observed a discordant pattern between mitochondrial and nuclear DNA, as all individuals carrying one of the divergent mitochondrial lineages grouped at the nuclear level with widespread, well-recognized species coexisting at the same localities (Daphnia galeata, D. longispina, and D. cucullata).

CONCLUSIONS

A likely explanation for this pattern is the introgression of the mitochondrial genome of undescribed taxa into the common species, either in the distant past or after long-distance dispersal. The occurrence of highly divergent but rare mtDNA lineages in the gene pool of widespread species would suggest that hybridization and introgression in the D. longispina species complex is frequent even across large phylogenetic distances, and that discoveries of such distinct clades must be interpreted with caution. However, maintenance of ancient polymorphisms through selection is another plausible alternative that may cause the observed discordance and cannot be entirely excluded.

Digging up the roots of an insular hotspot of genetic diversity: decoupled mito-nuclear histories in the evolution of the Corsican-Sardinian endemic lizard Podarcis tiliguerta

BACKGROUND

Mediterranean islands host a disproportionately high level of biodiversity and endemisms. Growing phylogeographic evidence on island endemics has unveiled unexpectedly complex patterns of intra-island diversification, which originated at diverse spatial and temporal scales. We investigated multilocus genetic variation of the Corsican-Sardinian endemic lizard Podarcis tiliguerta with the aim of shedding more light on the evolutionary processes underlying the origin of Mediterranean island biodiversity.

RESULTS

We analysed DNA sequences of mitochondrial (12S and nd4) and nuclear (acm4 and mc1r) gene fragments in 174 individuals of P. tiliguerta from 81 localities across the full range of the species in a geographic and genealogical framework. We found surprisingly high genetic diversity both at mitochondrial and nuclear loci. Seventeen reciprocally monophyletic allopatric mitochondrial haplogroups were sharply divided into four main mitochondrial lineages (two in Corsica and two in Sardinia) of Miocene origin. In contrast, shallow divergence and shared diversity within and between islands was observed at the nuclear loci. We evaluated alternative biogeographic and evolutionary scenarios to explain such profound discordance in spatial and phylogenetic patterning between mitochondrial and nuclear genomes. While neutral models provided unparsimonious explanations for the observed pattern, the hypothesis of environmental selection driving mitochondrial divergence in the presence of nuclear gene flow is favoured.

CONCLUSIONS

Our study on the genetic variation of P. tiliguerta reveals surprising levels of diversity underlining a complex phylogeographic pattern with a striking example of mito-nuclear discordance. These findings have profound implications, not only for the taxonomy and conservation of P. tiliguerta. Growing evidence on deep mitochondrial breaks in absence of geographic barriers and of climatic factors associated to genetic variation of Corsican-Sardinian endemics warrants additional investigation on the potential role of environmental selection driving the evolution of diversity hotspots within Mediterranean islands.

Islands within an island: Population genetic structure of the endemic Sardinian newt, Euproctus platycephalus

The identification of historic and contemporary barriers to dispersal is central to the conservation of endangered amphibians, but may be hindered by their complex life history and elusive nature. The complementary information generated by mitochondrial (mtDNA) and microsatellite markers generates a valuable tool in elucidating population structure and the impact of habitat fragmentation. We applied this approach to the study of an endangered montane newt, Euproctus platycephalus. Endemic to the Mediterranean island of Sardinia, it is threatened by anthropogenic activity, disease, and climate change. We have demonstrated a clear hierarchy of structure across genetically divergent and spatially distinct subpopulations. Divergence between three main mountain regions dominated genetic partitioning with both markers. Mitochondrial phylogeography revealed a deep division dating to ca. 1 million years ago (Mya), isolating the northern region, and further differentiation between the central and southern regions ca. 0.5 Mya, suggesting an association with Pleistocene severe glacial oscillations. Our findings are consistent with a model of southward range expansion during glacial periods, with postglacial range retraction to montane habitat and subsequent genetic isolation. Microsatellite markers revealed further strong population structure, demonstrating significant divergence within the central region, and partial differentiation within the south. The northern population showed reduced genetic diversity. Discordance between mitochondrial and microsatellite markers at this scale indicated a further complexity of population structure, in keeping with male-biased dispersal and female philopatry. Our study underscores the need to elucidate cryptic population structure in the ecology and conservation strategies for endangered island-restricted amphibians, especially in the context of disease and climate change.

Divergence maintained by climatic selection despite recurrent gene flow: a case study of Castanopsis carlesii (Fagaceae)

Local adaptation to different environments has the potential to maintain divergence between populations despite recurrent gene flow and is an important driver for generating biological diversity. In this study, we investigate the role of adaptation in the maintenance of two parapatric varieties of a forest tree. We used sequence variation of chloroplastic DNA and restriction site-associated DNA to investigate the genetic structure of two varieties of Castanopsis carlesii in subtropical China and relate it to climatic variation. We used niche reconstruction methods to investigate niche differentiation between the two varieties and to estimate the past distribution of this species. A deep divergence was observed between the two varieties, but evidence of introgression and genetic admixture was detected in two phenotypically and geographically intermediate populations. Niche reconstruction suggests that the distribution of the two varieties was disjunct during periods of global cooling and that the two varieties occupy significantly different niches. The genetic structure was mainly driven by environmental factors, and 13 outlier loci under divergent selection were correlated with climatic variation. These results suggest that the two varieties evolved in allopatry and came back into secondary contact after the last glacial maximum and that they are an evolutionary example of divergence maintained by climatic selection despite recurrent gene flow.

A comprehensive phylogeography of the Hyles euphorbiae complex (Lepidoptera: Sphingidae) indicates a 'glacial refuge belt'

We test the morphology based hypothesis that the Western Palaearctic spurge hawkmoths represent two species, the Eurasian H. euphorbiae and Afro-Macaronesian H. tithymali. It has been suggested that these species merged into several hybrid swarm populations, although a mitochondrial phylogeography revealed substructure with local differentiation. We analysed a three-gene mt-dataset (889 individuals) and 12 microsatellite loci (892 individuals). Microsatellite analyses revealed an overall weak differentiation and corroborated the superordinate division into two clusters. The data indicate that the populations studied belong to only one species according to the biological species concept, refuting the opening hypothesis. A future taxonomic revision appears necessary to reflect the division into two subgroups. Ancestral mitochondrial polymorphisms are retained in H. euphorbiae, indicating gene flow within a broad 'glacial refuge belt' and ongoing postglacial gene flow. Diverse patterns of extensive mito-nuclear discordance in the Mediterranean and the Middle East presumably evolved by more recent processes. This discordance indicates introgression of H. tithymali-related mitochondrial haplogroups, accompanied (to a lesser degree) by nuclear alleles, into Italian and Aegean H. euphorbiae populations as recently as the late Holocene. The complex mosaic of divergence and reintegration is assumed to have been influenced by locally differing environmental barriers to gene flow.

High genetic diversity in the Culex pipiens complex from a West Nile Virus epidemic area in Southern Europe

BACKGROUND

The Culex pipiens complex includes the most widespread mosquito species in the world. Cx. pipiens is the primary vector of the West Nile Virus (WNV) in Europe and North America. Cases of WNV have been recorded in Italy since 1998. In particular, wet areas along the Po River are considered some of the most WNV affected areas in Italy. Here, we analyzed the genetic structure of ten Cx. pipiens populations collected in the last part of the Po River including the Delta area.

METHODS

We assessed the genetic variability of two mitochondrial markers, cytochrome oxidase 1 (COI) and 2 (COII), for a total of 1200 bp, and one nuclear marker, a fragment of acetylcholinesterase-2 (ace-2), 502 bp long. The effect of the landscape features was evaluated comparing haplotype and nucleotide diversity with the landscape composition.

RESULTS

The analysis showed a high genetic diversity in both COI and COII gene fragments mainly shared by the populations in the Delta area. The COI-COII network showed that the set of haplotypes found was grouped into three main supported lineages with the higher genetic variability gathered in two of the three lineages. By contrast, ace-2 fragment did not show the same differentiation, displaying alleles grouped in a single clade. Finally, a positive correlation between mitochondrial diversity and natural wetland areas was found.

CONCLUSIONS

The high mitochondrial genetic diversity found in Cx. pipiens populations from the Po River Delta contrasts with the low variability of inland populations. The different patterns of genetic diversity found comparing mitochondrial and nuclear markers could be explained by factors such as differences in effective population size between markers, sex biased dispersal or lower fitness of dispersing females. Moreover, the correlation between genetic diversity and wetland areas is consistent with ecosystem stability and lack of insecticide pressure characteristic of this habitat. The mtDNA polymorphism found in the Po River Delta is even more interesting due to possible linkages between the mitochondrial lineages and different biting behaviors of the mosquitoes influencing their vector ability of arboviral infections.

Determining threshold values for barcoding fungi: lessons from Cortinarius (Basidiomycota), a highly diverse and widespread ectomycorrhizal genus

Different distance-based threshold selection approaches were used to assess and compare use of the internal transcribed spacer (ITS) region to distinguish among 901 Cortinarius species represented by >3000 collections. Sources of error associated with genetic markers and selection approaches were explored and evaluated using MOTUs from genus and lineage based-alignments. Our study indicates that 1%-2% more species can be distinguished by using the full-length ITS barcode as compared to either the ITS1 or ITS2 regions alone. Optimal threshold values for different picking approaches and genetic marker lengths inferred from a subset of species containing major lineages ranged from 97.0% to 99.5% sequence similarity using clustering optimization and UNITE SH, and from 1% to 2% sequence dissimilarity with CROP. Errors for the optimal cutoff ranged from 0% to 70%, and these can be reduced to a maximum of 22% when excluding species lacking a barcode gap. A threshold value of 99% is suitable for distinguishing species in the majority of lineages in the genus using the entire ITS region but only 90% of the species could be identified using just the ITS1 or ITS2 region. Prior identification of species, lacking barcode gaps and their subsequent separate analyses, maximized the accuracy of threshold approaches.

A tale of two haplotype groups: Evaluating the New World Junonia ring species hypothesis using the distribution of divergent COI haplotypes

The New World Junonia butterflies are a possible ring species with a circum-Caribbean distribution. Previous reports suggest a steady transition between North and South American forms in Mesoamerica, but in Cuba the forms were thought to co-exist without interbreeding representing the overlapping ends of the ring. Three criteria establish the existence of a ring species: a ring-shaped geographic distribution, gene flow among intervening forms, and genetic isolation in the region of range overlap. We evaluated mitochondrial cytochrome oxidase I haplotypes in Junonia from 9 species in the Western Hemisphere to test the Junonia ring species hypothesis. Junonia species are generally not monophyletic with respect to COI haplotypes, which are shared across species. However, two major COI haplotype groups exist. Group A predominates in South America, and Group B predominates in North and Central America. Therefore, COI haplotypes can be used to assess the degree of genetic influence a population receives from each continent. Junonia shows a ring-shaped distribution around the Caribbean, and evidence is consistent with gene flow among forms of Junonia, including those from Mesoamerica. However, we detected no discontinuity in gene flow in Cuba or elsewhere in the Caribbean consistent with genetic isolation in the region of overlap. Though sampling is still very limited in the critical region, the only remaining possiblity for a circum-Caribbean discontinuity in gene flow is at the Isthmus of Panama, where there may be a transition from 98% Group B haplotypes in Costa Rica to 85-100% Group A haplotypes in South America.

Three-dimensional post-glacial expansion and diversification of an exploited oceanic fish

Vertical divergence in marine organisms is being increasingly documented, yet much remains to be carried out to understand the role of depth in the context of phylogeographic reconstruction and the identification of management units. An ideal study system to address this issue is the beaked redfish, Sebastes mentella – one of four species of ‘redfish’ occurring in the North Atlantic – which is known for a widely distributed ‘shallow‐pelagic’ oceanic type inhabiting waters between 250 and 550 m, and a more localized ‘deep‐pelagic’ population dwelling between 550 and 800 m, in the oceanic habitat of the Irminger Sea. Here, we investigate the extent of population structure in relation to both depth and geographic spread of oceanic beaked redfish throughout most of its distribution range. By sequencing the mitochondrial control region of 261 redfish collected over a decadal interval, and combining 160 rhodopsin coding nuclear sequences and previously genotyped microsatellite data, we map the existence of two strongly divergent evolutionary lineages with significantly different distribution patterns and historical demography, and whose genetic variance is mostly explained by depth. Combined genetic data, analysed via independent approaches, are consistent with a Late Pleistocene lineage split, where segregation by depth probably resulted from the interplay of climatic and oceanographic processes with life history and behavioural traits. The ongoing process of diversification in North Atlantic S. mentella may serve as an ‘hourglass’ to understand speciation and adaptive radiation in Sebastes and in other marine taxa distributed across a depth gradient.

On the adaptive value of cytoplasmic genomes in plants

Is DNA variation maintained in organelle genomes selectively neutral? The answer to this question has important implications for many aspects of ecology and evolution. While traditionally the answer has been 'yes', recent studies in animals have shown that, on the contrary, mitochondrial DNA polymorphism is frequently adaptive. In plants, however, the neutrality assumption has not been strongly challenged. Here, we begin with a critical evaluation of arguments in favour of this long-held view. We then discuss the latest empirical evidence for the opposing prediction that sequence variation in plant cytoplasmic genomes is frequently adaptive. While outstanding research progress is being made towards understanding this fundamental topic, we highlight the need for studies that combine information ranging from field experiments to physiology to molecular evolutionary biology. Such an interdisciplinary approach provides a means for determining the frequency, drivers and evolutionary significance of adaptive organelle DNA variation.

Population differentiation and species formation in the deep sea: the potential role of environmental gradients and depth

Ecological speciation probably plays a more prominent role in diversification than previously thought, particularly in marine ecosystems where dispersal potential is great and where few obvious barriers to gene flow exist. This may be especially true in the deep sea where allopatric speciation seems insufficient to account for the rich and largely endemic fauna. Ecologically driven population differentiation and speciation are likely to be most prevalent along environmental gradients, such as those attending changes in depth. We quantified patterns of genetic variation along a depth gradient (1600-3800m) in the western North Atlantic for a protobranch bivalve (Nuculaatacellana) to test for population divergence. Multilocus analyses indicated a sharp discontinuity across a narrow depth range, with extremely low gene flow inferred between shallow and deep populations for thousands of generations. Phylogeographical discordance occurred between nuclear and mitochondrial loci as might be expected during the early stages of species formation. Because the geographic distance between divergent populations is small and no obvious dispersal barriers exist in this region, we suggest the divergence might reflect ecologically driven selection mediated by environmental correlates of the depth gradient. As inferred for numerous shallow-water species, environmental gradients that parallel changes in depth may play a key role in the genesis and adaptive radiation of the deep-water fauna.

Evolutionary branching in complex landscapes

Divergent adaptation to different environments can promote speciation, and it is thus important to consider spatial structure in models of speciation. Earlier theoretical work, however, has been limited to particularly simple types of spatial structure (linear environmental gradients and spatially discrete metapopulations), leaving unaddressed the effects of more realistic patterns of landscape heterogeneity, such as nonlinear gradients and spatially continuous patchiness. To elucidate the consequences of such complex landscapes, we adapt an established spatially explicit individual-based model of evolutionary branching. We show that branching is most probable at intermediate levels of various types of heterogeneity and that different types of heterogeneity have, to some extent, additive effects in promoting branching. In contrast to such additivity, we find a novel refugium effect in which refugia in hostile environments provide opportunities for colonization, thus increasing the probability of branching in patchy landscapes. Effects of patchiness depend on the scale of patches relative to dispersal. Providing a needed connection to empirical research on biodiversity and conservation policy, we introduce empirically accessible spatial environmental metrics that quantitatively predict a landscape's branching propensity.

Gene flow at major transitional areas in sea bass (Dicentrarchus labrax) and the possible emergence of a hybrid swarm

The population genetic structure of sea bass (Dicentrarchus labrax) along a transect from the Atlantic Ocean (AO) to the Eastern Mediterranean (EM) Sea differs from that of most other marine taxa in this area. Three populations (AO, Western Mediterranean [WM], EM) are recognized today, which were originally two allopatric populations. How two ancestral genetic units have evolved into three distinct units has not been addressed yet. Therefore, to investigate mechanisms that lead to the emergence of the central WM population, its current status, and its connectivity with the two parental populations, we applied 20 nuclear loci that were either gene associated or gene independent. Results confirmed the existence of three distinct gene pools, with higher differentiation at two transitional areas, the Almeria-Oran Front (AOF) and of the Siculo-Tunisian Strait (STS), than within any population. Significant linkage disequilibrium and heterozygote excess indicated that the STS is probably another tension zone, as already described for the AOF. Neutrality tests fail to reveal marker loci that could be driven by selection within or among metapopulations, except for locus DLA0068. Collectively, results support that the central WM population arose by trapping two tensions zones at distinct geographic locations of limited connectivity. Population assignment further revealed that WM individuals were more introgressed than individuals from the other two metapopulations. This suggests that this population might result from hybrid swarming, and was or is still seeded by genes received through the filter of each tension zone.

Deep sympatric mitochondrial divergence without reproductive isolation in the common redstart Phoenicurus phoenicurus

Mitochondrial DNA usually shows low sequence variation within and high sequence divergence among species, which makes it a useful marker for phylogenetic inference and DNA barcoding. A previous study on the common redstart (Phoenicurus phoenicurus) revealed two very different mtDNA haplogroups (5% K2P distance). This divergence is comparable to that among many sister species; however, both haplogroups coexist and interbreed in Europe today. Herein, we describe the phylogeographic pattern of these lineages and test hypotheses for how such high diversity in mtDNA has evolved. We found no evidence for mitochondrial pseudogenes confirming that both haplotypes are of mitochondrial origin. When testing for possible reproductive barriers, we found no evidence for lineage-specific assortative mating and no difference in sperm morphology, indicating that they are not examples of cryptic species, nor likely to reflect the early stages of speciation. A gene tree based on a short fragment of cytochrome c oxidase subunit 1 from the common redstart and 10 other Phoenicurus species, showed no introgression from any of the extant congenerics. However, introgression from an extinct congeneric cannot be excluded. Sequences from two nuclear introns did not show a similar differentiation into two distinct groups. Mismatch distributions indicated that the lineages have undergone similar demographic changes. Taken together, these results confirm that deeply divergent mitochondrial lineages can coexist in biological species. Sympatric mtDNA divergences are relatively rare in birds, but the fact that they occur argues against the use of threshold mtDNA divergences in species delineation.