Cryptic Species Due to Hybridization: A Combined Approach to Describe a New Species (Carex: Cyperaceae)

Genomic, spatial and morphometric data for discrimination of four species in the Mediterranean Tamus clade of yams (Dioscorea, Dioscoreaceae)

BACKGROUND AND AIMS

Among the numerous pantropical species of the yam genus, Dioscorea, only a small group occurs in the Mediterranean basin, including two narrow Pyrenean endemics (Borderea clade) and two Mediterranean-wide species (D. communis and D. orientalis, Tamus clade). However, several currently unrecognized species and infraspecific taxa have been described in the Tamus clade due to significant morphological variation associated with D. communis. Our overarching aim was to investigate taxon delimitation in the Tamus clade using an integrative approach combining phylogenomic, spatial and morphological data.

METHODS

We analysed 76 herbarium samples using Hyb-Seq genomic capture to sequence 260 low-copy nuclear genes and plastomes, together with morphometric and environmental modelling approaches.

KEY RESULTS

Phylogenomic reconstructions confirmed that the two previously accepted species of the Tamus clade, D. communis and D. orientalis, are monophyletic and form sister clades. Three subclades showing distinctive geographic patterns were identified within D. communis. These subclades were also identifiable from morphometric and climatic data, and introgression patterns were inferred between subclades in the eastern part of the distribution of D. communis.

CONCLUSIONS

We propose a taxonomy that maintains D. orientalis, endemic to the eastern Mediterranean region, and splits D. communis sensu lato into three species: D. edulis, endemic to Macaronesia (Canary Islands and Madeira); D. cretica, endemic to the eastern Mediterranean region; and D. communis sensu stricto, widespread across western and central Europe. Introgression inferred between D. communis s.s. and D. cretica is likely to be explained by their relatively recent speciation at the end of the Miocene, disjunct isolation in eastern and western Mediterranean glacial refugia and a subsequent westward recolonization of D. communis s.s. Our study shows that the use of integrated genomic, spatial and morphological approaches allows a more robust definition of species boundaries and the identification of species that previous systematic studies failed to uncover.

Genetic purity of a rear-edge population of Carex podogyna Franch. et Sav. (Cyperaceae) maintained under interspecific hybridization

Interspecific hybridization is a critical issue in conservation biology because it may drive small populations to extinction through direct or indirect processes. In this study, to develop a conservation strategy for an endangered rear-edge population of Carex podogyna in Ashiu, Kyoto, Japan, we performed a molecular genetic analysis of the wild population and an ex-situ population established from wild seeds. Microsatellite genotypic data revealed a complete loss of genetic diversity in the wild population, suggesting that it has long been prone to genetic drift due to isolation as a small population. In contrast, microsatellite analysis of 13 ex-situ individuals detected multiple alleles that are not harbored in the wild C. podogyna population. Sequence analysis revealed that these individuals are likely natural hybrids between C. podogyna and a co-occurring species, C. curvicollis, although established hybrids have never been found in the natural habitat. Based on our observation of variegated leaves in hybrid individuals, we propose that hybrids have been excluded by natural selection and/or interspecific competition caused by insufficient productivity of photosynthesis, although other genetic and ecological factors may also be influential. Overall, this study indicates that natural mechanisms selectively removing the hybrids have maintained the genetic purity of this rear-edge population of C. podogyna, and also emphasizes the importance of genetic assessment in ex-situ conservation programs.

Hybrid zone of a tree in a Cerrado/Atlantic Forest ecotone as a hotspot of genetic diversity and conservation

The Cerrado, the largest Neotropical savanna, and the Brazilian Atlantic Forest form large ecotonal areas where savanna and forest habitats occupy adjacent patches with closely related species occurring side by side, providing opportunities for hybridization. Here, we investigated the evolutionary divergence between the savanna and forest ecotypes of the widely distributed tree Plathymenia reticulata (n = 233 individuals). Genetic structure analysis of P. reticulata was congruent with the recognition of two ecotypes, whose divergence captured the largest proportion of genetic variance in the data (F CT = 0.222 and F ST = 0.307). The ecotonal areas between the Cerrado and the Atlantic Forest constitute a hybrid zone in which a diversity of hybrid classes was observed, most of them corresponding to second-generation hybrids (F2) or backcrosses. Gene flow occurred mainly toward the forest ecotype. The genetic structure was congruent with isolation by environment, and environmental correlates of divergence were identified. The observed pattern of high genetic divergence between ecotypes may reflect an incipient speciation process in P. reticulata. The low genetic diversity of the P. reticulata forest ecotype indicate that it is threatened in areas with high habitat loss on Atlantic Forest. In addition, the high divergence from the savanna ecotype suggests it should be treated as a different unit of management. The high genetic diversity found in the ecotonal hybrid zone supports the view of ecotones as important areas for the origin and conservation of biodiversity in the Neotropics.

Gene flow creates a mirage of cryptic species in a Southeast Asian spotted stream frog complex

Most new cryptic species are described using conventional tree- and distance-based species delimitation methods (SDMs), which rely on phylogenetic arrangements and measures of genetic divergence. However, although numerous factors such as population structure and gene flow are known to confound phylogenetic inference and species delimitation, the influence of these processes is not frequently evaluated. Using large numbers of exons, introns, and ultraconserved elements obtained using the FrogCap sequence-capture protocol, we compared conventional SDMs with more robust genomic analyses that assess population structure and gene flow to characterize species boundaries in a Southeast Asian frog complex (Pulchrana picturata). Our results showed that gene flow and introgression can produce phylogenetic patterns and levels of divergence that resemble distinct species (up to 10% divergence in mitochondrial DNA). Hybrid populations were inferred as independent (singleton) clades that were highly divergent from adjacent populations (7%-10%) and unusually similar (<3%) to allopatric populations. Such anomalous patterns are not uncommon in Southeast Asian amphibians, which brings into question whether the high levels of cryptic diversity observed in other amphibian groups reflect distinct cryptic species-or, instead, highly admixed and structured metapopulation lineages. Our results also provide an alternative explanation to the conundrum of divergent (sometimes nonsister) sympatric lineages-a pattern that has been celebrated as indicative of true cryptic speciation. Based on these findings, we recommend that species delimitation of continuously distributed "cryptic" groups should not rely solely on conventional SDMs, but should necessarily examine population structure and gene flow to avoid taxonomic inflation.

The protected tree Dimorphandra wilsonii (Fabaceae) is a population of inter-specific hybrids: recommendations for conservation in the Brazilian Cerrado/Atlantic Forest ecotone

BACKGROUNDS AND AIMS

Dimorphandra wilsonii Rizzini, a critically endangered and protected tree, has a restricted distribution in the ecotone between the Cerrado and the Atlantic Forest in south-eastern Brazil. In this area, it co-occurs with D. mollis Benth., a common tree from the Cerrado, and D. exaltata Schott., a rare tree from the Atlantic Forest. Previous studies of D. wilsonii indicated heterozygosity excess at the individual level. Field observation of some intermediate phenotypes between D. wilsonii and both congeners suggests hybridization of D. wilsonii with D. mollis and/or D. exaltata. Here, we tested the hypothesis that D. wilsonii may have originated from hybridization between D. exaltata and D. mollis. We also performed cytogenetic analysis to examine if the heterozygosity excess could be explained by polyploidy in D. wilsonii.

METHODS

We evaluated the genetic diversity and population structure of D. wilsonii using 11 nuclear simple sequence repeats (SSRs) genotyped in 152 individuals sampled across the taxon's range. We performed comparative genetic analyses using overlapping SSR markers between D. wilsonii and previously published SSR data in D. mollis and D. exaltata to subsequently perform a series of allelic comparisons, multivariate and Bayesian analysis.

KEY RESULTS

Our results suggest that D. wilsonii individuals are most likely to correspond to F1 hybrids between D. exaltata and D. mollis. Cytogenetic analysis indicated that D. wilsonii is diploid with the same chromosome number as D. mollis (2n = 2x = 28).

CONCLUSIONS

Our study raises questions about the taxonomic status and the evolutionary future of D. wilsonii. We suggest that the conservation and management strategy for D. wilsonii should be revised and that it should take into account both parental Dimorphandra species in the ecotone, with special emphasis on the threatened D. exaltata. Finally, this study highlights the value of genetic information for the design of conservation strategies.

Allopatric speciation despite historical gene flow: Divergence and hybridization in Carex furva and C. lucennoiberica (Cyperaceae) inferred from plastid and nuclear RAD-seq data

Gene flow among incipient species can act as a creative or destructive force in the speciation process, generating variation on which natural selection can act while, potentially, undermining population divergence. The flowering plant genus Carex exhibits a rapid and relatively recent radiation with many species limits still unclear. This is the case with the Iberian Peninsula (Spain and Portugal)-endemic C. lucennoiberica, which lay unrecognized within Carex furva until its recent description as a new species. In this study, we test how these species were impacted by interspecific gene flow during speciation. We sampled the full range of distribution of C. furva (15 individuals sampled) and C. lucennoiberica (88 individuals), sequenced two cpDNA regions (atpI-atpH, psbA-trnH) and performed genomic sequencing of 45,100 SNPs using restriction site-associated DNA sequencing (RAD-seq). We utilized a set of partitioned D-statistic tests and demographic analyses to study the degree and direction of introgression. Additionally, we modelled species distributions to reconstruct changes in range distribution during glacial and interglacial periods. Plastid, nuclear and morphological data strongly support divergence between species with subsequent gene flow. Combined with species distribution modelling, these data support a scenario of allopatry leading to species divergence, followed by secondary contact and gene flow due to long-distance dispersal and/or range expansions and contractions in response to Quaternary glacial cycles. We conclude that this is a case of allopatric speciation despite historical secondary contacts, which could have temporally influenced the speciation process, contributing to the knowledge of forces that are driving or counteracting speciation.

Correction: Cryptic Species Due to Hybridization: A Combined Approach to Describe a New Species (Carex: Cyperaceae)

[This corrects the article DOI: 10.1371/journal.pone.0166949.].