Species delimitation under the general lineage concept: an empirical example using wild North American hops (Cannabaceae: Humulus lupulus)

Diversification of the African legless skinks in the subfamily Acontinae (Family Scincidae)

Cladogenic diversification is often explained by referring to climatic oscillations and geomorphic shifts that cause allopatric speciation. In this regard, southern Africa retains a high level of landscape heterogeneity in vegetation, geology, and rainfall patterns. The legless skink subfamily Acontinae occurs broadly across the southern African subcontinent and therefore provides an ideal model group for investigating biogeographic patterns associated with the region. A robust phylogenetic study of the Acontinae with comprehensive coverage and adequate sampling of each taxon has been lacking up until now, resulting in unresolved questions regarding the subfamily's biogeography and evolution. In this study, we used multi-locus genetic markers (three mitochondrial and two nuclear) with comprehensive taxon coverage (all currently recognized Acontinae species) and adequate sampling (multiple specimens for most taxa) of each taxon to infer a phylogeny for the subfamily. The phylogeny retrieved four well-supported clades in Acontias and supported the monophyly of Typhlosaurus. Following the General Lineage Concept (GLC), many long-standing phylogenetic enigmas within Acontias occidentalis and the A. kgalagadi, A. lineatus and A. meleagris species complexes, and within Typhlosaurus were resolved. Our species delimitation analyses suggest the existence of hidden taxa in the A. occidentalis, A. cregoi and A. meleagris species groups, but also suggest that some currently recognized species in the A. lineatus and A. meleagris species groups, and within Typhlosaurus, should be synonymised. We also possibly encountered "ghost introgression" in A. occidentalis. Our inferred species tree revealed a signal of gene flow, which implies possible cross-over in some groups. Fossil evidence calibration dating results showed that the divergence between Typhlosaurus and Acontias was likely influenced by cooling and increasing aridity along the southwest coast in the mid-Oligocene caused by the opening of the Drake Passage. Further cladogenesis observed in Typhlosaurus and Acontias was likely influenced by Miocene cooling, expansion of open habitat, uplifting of the eastern Great Escarpment (GE), and variation in rainfall patterns, together with the effect of the warm Agulhas Current since the early Miocene, the development of the cold Benguela Current since the late Miocene, and their co-effects. The biogeographic pattern of the Acontinae bears close resemblance to that of other herpetofauna (e.g., rain frogs and African vipers) in southern Africa.

Lineage and role in integrative taxonomy of a heterotrophic orchid complex

Lineage-based species definitions applying coalescent approaches to species delimitation have become increasingly popular. Yet, the application of these methods and the recognition of lineage-only definitions have recently been questioned. Species delimitation criteria that explicitly consider both lineages and evidence for ecological role shifts provide an opportunity to incorporate ecologically meaningful data from multiple sources in studies of species boundaries. Here, such criteria were applied to a problematic group of mycoheterotrophic orchids, the Corallorhiza striata complex, analysing genomic, morphological, phenological, reproductive-mode, niche, and fungal host data. A recently developed method for generating genomic polymorphism data-ISSRseq-demonstrates evidence for four distinct lineages, including a previously unidentified lineage in the Coast Ranges and Cascades of California and Oregon, USA. There is divergence in morphology, phenology, reproductive mode, and fungal associates among the four lineages. Integrative analyses, conducted in population assignment and redundancy analysis frameworks, provide evidence of distinct genomic lineages and a similar pattern of divergence in the extended data, albeit with weaker signal. However, none of the extended data sets fully satisfy the condition of a significant role shift, which requires evidence of fixed differences. The four lineages identified in the current study are recognized at the level of variety, short of comprising different species. This study represents the most comprehensive application of lineage + role to date and illustrates the advantages of such an approach.

Differences in Leaf Chemistry and Glandular Trichome Density between Wild Southwestern American Hop (Humulus neomexicanus) and Commercial Hop Cultivars

The female flowers ("cones") of the hop plant (Humulus L.) produce compounds that contribute to the flavor and other properties of beer. Hop leaves and cones produce many of the same compounds, which also confer agronomic traits such as insect and disease resistance. Targeted and untargeted ultraperformance liquid chromatography-quadrupole time-of-flight-mass spectrometry with Waters MS^E technology (UPLC-QTof-MS^E) metabolomics were used to compare leaf phytochemical compositions of greenhouse-grown southwestern American wild Humulus neomexicanus (A. Nelson and Cockerell) Rydb. against a group of commercial hop cultivars consisting of both pure European Humulus lupulus L. and European-North American hybrids. Principal component analysis showed a clear distinction in chemical profiles between the two groups. H. neomexicanus leaves had a significantly higher content of total α acids (p = 4.4 × 10^-9), total bitter acids (p = 2.6 × 10^-6), cohumulone (p = 1.0 × 10^-13), humulone + adhumulone (p = 9.1 × 10^-4), and the prenylflavonoids xanthohumol (p = 0.013) and desmethylxanthohumol (p = 0.029) as well as significantly higher densities of glandular trichomes (p = 1.3 × 10^-6), the biosynthetic site of those compounds. Most flavonol glycosides measured were also significantly more abundant in H. neomexicanus (p = 1.5 × 10^-22 to 0.0027), whereas phenolic acids were consistently, but generally nonsignificantly (p > 0.05), more abundant in the cultivars. The higher bitter acid, prenylflavonoid, and flavonol glycoside content of H. neomexicanus leaves may help to confer more favorable insect and disease-resistance properties.

Reversing extinction trends: new uses of (old) herbarium specimens to accelerate conservation action on threatened species

Although often not collected specifically for the purposes of conservation, herbarium specimens offer sufficient information to reconstruct parameters that are needed to designate a species as 'at-risk' of extinction. While such designations should prompt quick and efficient legal action towards species recovery, such action often lags far behind and is mired in bureaucratic procedure. The increase in online digitization of natural history collections has now led to a surge in the number new studies on the uses of machine learning. These repositories of species occurrences are now equipped with advances that allow for the identification of rare species. The increase in attention devoted to estimating the scope and severity of the threats that lead to the decline of such species will increase our ability to mitigate these threats and reverse the declines, overcoming a current barrier to the recovery of many threatened plant species. Thus far, collected specimens have been used to fill gaps in systematics, range extent, and past genetic diversity. We find that they also offer material with which it is possible to foster species recovery, ecosystem restoration, and de-extinction, and these elements should be used in conjunction with machine learning and citizen science initiatives to mobilize as large a force as possible to counter current extinction trends.

A draft phased assembly of the diploid Cascade hop (Humulus lupulus) genome

Hop (Humulus lupulus L. var Lupulus) is a diploid, dioecious plant with a history of cultivation spanning more than one thousand years. Hop cones are valued for their use in brewing and contain compounds of therapeutic interest including xanthohumol. Efforts to determine how biochemical pathways responsible for desirable traits are regulated have been challenged by the large (2.8 Gb), repetitive, and heterozygous genome of hop. We present a draft haplotype-phased assembly of the Cascade cultivar genome. Our draft assembly and annotation of the Cascade genome is the most extensive representation of the hop genome to date. PacBio long-read sequences from hop were assembled with FALCON and partially phased with FALCON-Unzip. Comparative analysis of haplotype sequences provides insight into selective pressures that have driven evolution in hop. We discovered genes with greater sequence divergence enriched for stress-response, growth, and flowering functions in the draft phased assembly. With improved resolution of long terminal retrotransposons (LTRs) due to long-read sequencing, we found that hop is over 70% repetitive. We identified a homolog of cannabidiolic acid synthase (CBDAS) that is expressed in multiple tissues. The approaches we developed to analyze the draft phased assembly serve to deepen our understanding of the genomic landscape of hop and may have broader applicability to the study of other large, complex genomes.

Species Identification in the Rhododendron vernicosum-R. decorum Species Complex (Ericaceae)

Delimitating species boundaries is the primary aim of biological classification and could be critical for evaluating the evolving process of species and conserving biodiversity. Rhododendron is an iconic group with an extraordinary diversity in southwest China. However, it remains unknown whether the recorded species therein comprise independently evolving lineages or artificially delimitated morphological entities. In this study, we carried out species delimitation of four Rhododendron species in the R. vernicosum-R. decorum species complex based on morphological analyses and population genetic data from nuclear simple sequence repeats (SSR) markers. We randomly selected a total of 105 specimens of different individuals identified as four species across their distributional ranges to examine the statistically distinct phenotypic clusters based on 19 morphological traits. Similarly, we genotyped 55 individuals of four species from 21 populations using 15 SSR markers. The morphological analyses sorted R. decorum and the other three species into two different phenotypic clusters. The genetic clusters were consistent with the morphological clusters. However, we also recovered the third genetic cluster, comprising six R. vernicosum populations and containing the admixed genetic compositions of the other two distinct genetic clusters. This hybrid group was morphologically similar to the typical R. vernicosum (including the samples from its type specimen locality and both R. verruciferum and R. gonggashanense) but with more genetic ancestry from R. decorum. Based on our findings, we identify two distinct species and one putative hybrid group due to introgression in the R. vernicosum-R. decorum species complex. We propose to merge R. verruciferum and R. gonggashanense into R. vernicosum based on genetic compositions and our morphological analyses. The hybrid group inferred from our findings, however, needs further investigations.

Identification of tandem repeat families from long-read sequences of Humulus lupulus

Hop (Humulus lupulus L.) is known for its use as a bittering agent in beer and has a rich history of cultivation, beginning in Europe and now spanning the globe. There are five wild varieties worldwide, which may have been introgressed with cultivated varieties. As a dioecious species, its obligate outcrossing, non-Mendelian inheritance, and genomic structural variability have confounded directed breeding efforts. Consequently, understanding the hop genome represents a considerable challenge, requiring additional resources. In order to facilitate investigations into the transmission genetics of hop, we report here a tandem repeat discovery pipeline developed using k-mer filtering and dot plot analysis of PacBio long-read sequences from the hop cultivar Apollo. From this we identified 17 new and distinct tandem repeat sequence families, which represent candidates for FISH probe development. For two of these candidates, HuluTR120 and HuluTR225, we produced oligonucleotide FISH probes from conserved regions of and demonstrated their utility by staining meiotic chromosomes from wild hop, var. neomexicanus to address, for example, questions about hop transmission genetics. Collectively, these tandem repeat sequence families represent new resources suitable for development of additional cytogenomic tools for hop research.

Genetic diversity and biogeographic determinants of population structure in Araucaria angustifolia (Bert.) O. Ktze

Araucaria (Araucaria angustifolia (Bert.) O. Ktze) is a primarily dioecious species threatened with extinction that plays an important social and economic role especially in the southern region of Brazil. The aim of this work is to investigate the diversity and likely determinants of genetic lineages in this species for conservation management. For this, a collection of 30-year-old Araucaria was used. Accessions collected from 12 sites across the species range were analyzed, with ten individuals per site. The SSR genotyping was conducted with 15 loci and the data were analyzed using several complementary approaches. Descriptive statistics among sampling sites were used and diversity was partitioned non-hierarchically to estimate the size and composition of genetic clusters using a Bayesian assignment method. To explore possible biological implications of differences between Niche Models and habitat suitability, a series of statistical procedures were used, and tests were carried out using the software ENM Tools and Maxent. Populations from the southernmost zone showed higher genetic variation and a lower inbreeding coefficient compared to the northernmost zone, which may correlate with their isolation. A positive relation between genetic differentiation and geographic distance was observed. Two genetic groups (southernmost and northernmost zones) were evident. The Niche modelling showed separate ranges for each genetic lineage suggesting that differences in selection pressure may be playing a role in the apparent differentiation and may be adaptive. Finally, an evident correlation was observed between genetic data and habitat suitability. The two distinct groups observed must be considered as independent units for conservation and hybridization in breeding programs.

Climatic niche comparison across a cryptic species complex

According to current molecular evidence, the Chionaspis pinifoliae heterophyllae species complex has been recognized as 10 cryptic species. In this study, we construct potential distribution maps for seven cryptic species based on climatic variables. This was done to assess the main environmental factors that have contributed to the distribution map and test the degree of niche overlap across the seven cryptic species. We used MaxEnt to build the climatic niche models under climatic variables. For these models, the similarities and differences of the niches across the cryptic species were estimated. By comparing the potential distribution model of each cryptic species, our results suggested parapatric, sympatric and allopatry populations for this cryptic species complex. Our results showed high variability in niche overlap, and more often niche conservatism than niche divergence. The current species delimitation of the Chionaspis pinifoliae heterophyllae complex by molecular information and the hypothesis that the niche overlap in the sympatric population is higher than that of the allopatry population were supported based on the findings. This study will provide baseline data and a distribution range to facilitate the further control of these insects and formulate quarantine measures.

3D Molecular Cytology of Hop (Humulus lupulus) Meiotic Chromosomes Reveals Non-disomic Pairing and Segregation, Aneuploidy, and Genomic Structural Variation

Hop (Humulus lupulus L.) is an important crop worldwide, known as the main flavoring ingredient in beer. The diversifying brewing industry demands variation in flavors, superior process properties, and sustainable agronomics, which are the focus of advanced molecular breeding efforts in hops. Hop breeders have been limited in their ability to create strains with desirable traits, however, because of the unusual and unpredictable inheritance patterns and associated non-Mendelian genetic marker segregation. Cytogenetic analysis of meiotic chromosome behavior has also revealed conspicuous and prevalent occurrences of multiple, atypical, non-disomic chromosome complexes, including those involving autosomes in late prophase. To explore the role of meiosis in segregation distortion, we undertook 3D cytogenetic analysis of hop pollen mother cells stained with DAPI and FISH. We used telomere FISH to demonstrate that hop exhibits a normal telomere clustering bouquet. We also identified and characterized a new sub-terminal 180 bp satellite DNA tandem repeat family called HSR0, located proximal to telomeres. Highly variable 5S rDNA FISH patterns within and between plants, together with the detection of anaphase chromosome bridges, reflect extensive departures from normal disomic signal composition and distribution. Subsequent FACS analysis revealed variable DNA content in a cultivated pedigree. Together, these findings implicate multiple phenomena, including aneuploidy, segmental aneuploidy, or chromosome rearrangements, as contributing factors to segregation distortion in hop.

Divergence is not enough: the use of ecological niche models for the validation of taxon boundaries

Delimiting taxon boundaries is crucial for any evolutionary research and conservation regulation. In order to avoid mistaken description of species, the approach of integrative taxonomy recommends considering multidisciplinary lines of evidence, including ecology. Unfortunately, ecological data are often difficult to quantify objectively. Here we test and discuss the potential use of ecological niche models for validating taxon boundaries, using three pairs of closely related plant taxa endemic to the south-western Alps as a case study. We also discuss the application of ecological niche models for species delimitation and the implementation of different approaches. Niche overlap, niche equivalency and niche similarity were assessed both in multidimensional environmental space and in geographic space to look for differences in the niche of three pairs of closely related plant taxa. We detected a high degree of niche differentiation between taxa although this result seems not due to differences in habitat selection. The different statistical tests gave contrasting outcomes between environmental and geographic spaces. According to our results, niche divergence does not seem to support taxon boundaries at species level, but may have had important consequences for local adaptation and in generating phenotypic diversity at intraspecific level. Environmental space analysis should be preferred to geographic space as it provides more clear results. Even if the different analyses widely disagree in their conclusions about taxon boundaries, our study suggests that ecological niche models may help taxonomists to reach a decision.

Non-Mendelian Single-Nucleotide Polymorphism Inheritance and Atypical Meiotic Configurations are Prevalent in Hop

Hop ( L.) breeding programs seek to exploit genetic resources for bitter flavor, aroma, and disease resistance. However, these efforts have been thwarted by segregation distortion including female-biased sex ratios. To better understand the transmission genetics of hop, we genotyped 4512 worldwide accessions of hop, including cultivars, landraces, and over 100 wild accessions using a genotyping-by-sequencing (GBS) approach. From the resulting ∼1.2 million single-nucleotide polymorphisms (SNPs), prequalified GBS markers were validated by inferences in population structures and phylogeny. Analysis of pseudo-testcross (Pt) mapping data from F families revealed mixed patterns of Mendelian and non-Mendelian segregation. Three-dimensional (3D) cytogenetic analysis of late meiotic prophase nuclei from two wild and two cultivated hop revealed conspicuous and prevalent occurrences of multiple, atypical, nondisomic chromosome complexes including autosomes. We used genome-wide association studies (GWAS) and fixation index (F) analysis to demonstrate selection mapping of genetic loci for key traits including sex, bitter acids, and drought tolerance. Among the possible mechanisms underlying the observed segregation distortion from the genomic data analysis, the cytogenetic analysis points to meiotic chromosome behavior as one of the contributing factors. The findings shed light on long-standing questions on the unusual transmission genetics and phenotypic variation in hop, with major implications for breeding, cultivation, and the natural history of .

Species delimitation and phylogeography of the Pectenia species-complex: A misunderstood case of species-pairs in lichenized fungi, where reproduction mode does not delimit lineages

The main goal of this work was to study species boundaries in the genus Pectenia and elucidate the biogeographic history of the four currently accepted species. To accomplish this, we included 92 specimens across the range of Pectenia in Europe and northern Africa. We used three nuclear loci and assessed species circumscription using two Bayesian coalescent-based methods and the Bayes Factor approach. We reviewed the value of reproductive mode and other morphological features as predictors of monophyletic groups. Our results suggest that the production of asexual propagules and sexual structures are not characterizing monophyletic groups. The genus includes two morphologically well-characterized main lineages, where one lineage is composed by two sub-lineages that are with a case of cryptic speciation explained by a biogeographic pattern. We suggest treating the two lineages as two species, which are characterized by lobe structure: Pectenia plumbea and P. atlantica. Both of these species include samples with asexual propagules and apothecia, and thus do not correspond to any of the earlier morphologically defined species. The results of the biogeographic analysis indicate that the Mediterranean basin is the most likely ancestral distribution area of P. plumbea, whereas P. atlantica probably originated in Macaronesia.

Deeply divergent sympatric mitochondrial lineages of the earthworm Lumbricus rubellus are not reproductively isolated

Background

The accurate delimitation of species is essential to numerous areas of biological research. An unbiased assessment of the diversity, including the cryptic diversity, is of particular importance for the below ground fauna, a major component of global biodiversity. On the British Isles, the epigeic earthworm Lumbricus rubellus, which is a sentinel species in soil ecotoxicology, consists of two cryptic taxa that are differentiated in both the nuclear and the mitochondrial (mtDNA) genomes. Recently, several deeply divergent mtDNA lineages were detected in mainland Europe, but whether these earthworms also constitute cryptic species remains unclear. This information is important from an evolutionary perspective, but it is also essential for the interpretation and the design of ecotoxicological projects. In this study, we used genome-wide RADseq data to assess the reproductive isolation of the divergent mitochondrial lineages of L. rubellus that occur in sympatry in multiple localities in Central Europe.

Results

We identified five divergent (up to 16 % net p-distance) mitochondrial lineages of L. rubellus in sympatry. Because the clustering of the RADseq data was according to the population of origin and not the mtDNA lineage, reproductive isolation among the mtDNA lineages was not likely. Although each population contained multiple mtDNA lineages, subdivisions within the populations were not observed for the nuclear genome. The lack of fixed differences and sharing of the overwhelming majority of nuclear polymorphisms between localities, indicated that the populations did not constitute allopatric species. The nucleotide diversity within the populations was high, 0.7–0.8 %.

Conclusions

The deeply divergent mtDNA sympatric lineages of L. rubellus in Central Europe were not reproductively isolated groups. The earthworm L. rubellus, which is represented by several mtDNA lineages in continental Europe, apparently is a single highly polymorphic species rather than a complex of several cryptic species. This study demonstrated the critical importance of the use of multilocus nuclear data for the unbiased assessment of cryptic diversity and for the delimitation of species in soil invertebrates.

Electronic supplementary material

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

Limitations of climatic data for inferring species boundaries: insights from speckled rattlesnakes

Phenotypes, DNA, and measures of ecological differences are widely used in species delimitation. Although rarely defined in such studies, ecological divergence is almost always approximated using multivariate climatic data associated with sets of specimens (i.e., the "climatic niche"); the justification for this approach is that species-specific climatic envelopes act as surrogates for physiological tolerances. Using identical statistical procedures, we evaluated the usefulness and validity of the climate-as-proxy assumption by comparing performance of genetic (nDNA SNPs and mitochondrial DNA), phenotypic, and climatic data for objective species delimitation in the speckled rattlesnake (Crotalus mitchellii) complex. Ordination and clustering patterns were largely congruent among intrinsic (heritable) traits (nDNA, mtDNA, phenotype), and discordance is explained by biological processes (e.g., ontogeny, hybridization). In contrast, climatic data did not produce biologically meaningful clusters that were congruent with any intrinsic dataset, but rather corresponded to regional differences in atmospheric circulation and climate, indicating an absence of inherent taxonomic signal in these data. Surrogating climate for physiological tolerances adds artificial weight to evidence of species boundaries, as these data are irrelevant for that purpose. Based on the evidence from congruent clustering of intrinsic datasets, we recommend that three subspecies of C. mitchellii be recognized as species: C. angelensis, C. mitchellii, and C. Pyrrhus.

Delineating ecological boundaries of Hanuman langur species complex in peninsular India using MaxEnt modeling approach

Hanuman langur is one of the widely distributed and extensively studied non-human diurnal primates in India. Until recently it was believed to be a single species - Semnopithecus entellus. Recent molecular and morphological studies suggest that the Hanuman langurs consists of at least three species S. entellus, S. hypoleucos and S. priam. Furthermore, morphological studies suggested that both S. hypoleucos and S. priam have at least three subspecies in each. We explored the use of ecological niche modeling (ENM) to confirm the validity of these seven taxa and an additional taxon S. johnii belonging to the same genus. MaxEnt modeling tool was used with 19 bioclimatic, 12 vegetation and 6 hydrological environmental layers. We reduced total environmental variables to 14 layers after testing for collinearity and an independent test for model prediction was done using ENMTools. A total of 196 non-overlapping data points from primary and secondary sources were used as inputs for ENM. Results showed eight distinct ecological boundaries, corroborating the eight taxa mentioned above thereby confirming validity of these eight taxa. The study, for the first time provided ecological variables that determined the ecological requirements and distribution of members of the Hanuman langur species complex in the Indian peninsula.

Niche divergence versus neutral processes: combined environmental and genetic analyses identify contrasting patterns of differentiation in recently diverged pine species

Background and Aims

Solving relationships of recently diverged taxa, poses a challenge due to shared polymorphism and weak reproductive barriers. Multiple lines of evidence are needed to identify independently evolving lineages. This is especially true of long-lived species with large effective population sizes, and slow rates of lineage sorting. North American pines are an interesting group to test this multiple approach. Our aim is to combine cytoplasmic genetic markers with environmental information to clarify species boundaries and relationships of the species complex of Pinus flexilis, Pinus ayacahuite, and Pinus strobiformis.

Methods

Mitochondrial and chloroplast sequences were combined with previously obtained microsatellite data and contrasted with environmental information to reconstruct phylogenetic relationships of the species complex. Ecological niche models were compared to test if ecological divergence is significant among species.

Key Results and Conclusion

Separately, both genetic and ecological evidence support a clear differentiation of all three species but with different topology, but also reveal an ancestral contact zone between P. strobiformis and P. ayacahuite. The marked ecological differentiation of P. flexilis suggests that ecological speciation has occurred in this lineage, but this is not reflected in neutral markers. The inclusion of environmental traits in phylogenetic reconstruction improved the resolution of internal branches. We suggest that combining environmental and genetic information would be useful for species delimitation and phylogenetic studies in other recently diverged species complexes.

New metrics for comparison of taxonomies reveal striking discrepancies among species delimitation methods in Madascincus lizards

Delimiting and describing species is fundamental to numerous biological disciplines such as evolution, macroecology, and conservation. Delimiting species as independent evolutionary lineages may and often does yield different outcomes depending on the species criteria applied, but methods should be chosen that minimize the inference of objectively erroneous species limits. Several protocols exploit single-gene or multi-gene coalescence statistics, assignment tests or other rationales related to nuclear DNA (nDNA) allele sharing to automatically delimit species. We apply seven different species delimitation protocols to a taxonomically confusing group of Malagasy lizards (Madascincus), and compare the resulting taxonomies with two newly developed metrics: the Taxonomic index of congruence C tax which quantifies the congruence between two taxonomies, and the Relative taxonomic resolving power index R tax which quantifies the potential of an approach to capture a high number of species boundaries. The protocols differed in the total number of species proposed, between 9 and 34, and were also highly incongruent in placing species boundaries. The Generalized Mixed Yule-Coalescent approach captured the highest number of potential species boundaries but many of these were clearly contradicted by extensive nDNA admixture between sympatric mitochondrial DNA (mtDNA) haplotype lineages. Delimiting species as phenotypically diagnosable mtDNA clades failed to detect two cryptic species that are unambiguous due to a lack of nDNA gene flow despite sympatry. We also consider the high number of species boundaries and their placement by multi-gene Bayesian species delimitation as poorly reliable whereas the Bayesian assignment test approach provided a species delimitation highly congruent with integrative taxonomic practice. The present study illustrates the trade-off in taxonomy between reliability (favored by conservative approaches) and resolving power (favored by inflationist approaches). Quantifying excessive splitting is more difficult than quantifying excessive lumping, suggesting a priority for conservative taxonomies in which errors are more liable to be detected and corrected by subsequent studies.

Spatiotemporal reconstruction of the Aquilegia rapid radiation through next-generation sequencing of rapidly evolving cpDNA regions

Aquilegia is a well-known model system in the field of evolutionary biology, but obtaining a resolved and well-supported phylogenetic reconstruction for the genus has been hindered by its recent and rapid diversification. Here, we applied 454 next-generation sequencing to PCR amplicons of 21 of the most rapidly evolving regions of the plastome to generate c. 24 kb of sequences from each of 84 individuals from throughout the genus. The resulting phylogeny has well-supported resolution of the main lineages of the genus, although recent diversification such as in the European taxa remains unresolved. By producing a chronogram of the whole Ranunculaceae family based on published data, we inferred calibration points for dating the Aquilegia radiation. The genus originated in the upper Miocene c. 6.9 million yr ago (Ma) in Eastern Asia, and diversification occurred c. 4.8 Ma with the split of two main clades, one colonizing North America, and the other Western Eurasia through the mountains of Central Asia. This was followed by a back-to-Asia migration, originating from the European stock using a North Asian route. These results provide the first backbone phylogeny and spatiotemporal reconstruction of the Aquilegia radiation, and constitute a robust framework to address the adaptative nature of speciation within the group.

Delimiting genetic units in Neotropical toads under incomplete lineage sorting and hybridization

Background

Delimiting genetic units is useful to enhance taxonomic discovery and is often the first step toward understanding evolutionary mechanisms generating diversification. The six species within the Rhinella crucifer group of toads were defined under morphological criteria alone. Previous data suggest limited correspondence of these species to mitochondrial lineages, and morphological intergradation at transitions between forms suggests hybridization. Here we extensively sampled populations throughout the geographic distribution of the group and analyzed mitochondrial and nuclear sequence data to delimit genetic units using tree–based and allele frequency–based approaches.

Results

These approaches yielded complementary results, with allele frequency-based methods performing unexpectedly well given the limited number of loci examined. Both mitochondrial and nuclear markers supported a genetic structure of five units within the group, with three of the inferred units distributed within its main range, while two other units occur in separate isolates. The inferred units are mostly discordant with currently described forms: unequivocal association exists for only two of the six species in the group. Genetic evidence for hybridization exists for two pairs of units, with clear cyto–nuclear allele mixing observed in one case.

Conclusions

Our results confirmed that current taxonomy does not represent evolutionary units in the Rhinella crucifer group. Correspondence between genetically distinguishable units and the currently recognized species is only possible for Rhinella henseli and R. inopina. The recognition of other species relies on the reassessment of the geographic range of R. crucifer, the examination of the type series of R. ornata for hybrids, and on the use of additional markers to verify the genetic distinctiveness of R. abei. We state that R. pombali should not remain a valid species since its description appears to be based on hybrids, and that the name R. pombali should be considered a synonym of both R. crucifer and R. ornata. The fifth inferred but undescribed genetic unit may represent a new species. Our results underscore the potential of the R. crucifer species group to contribute to a better understanding of diversification processes and hybridization patterns in the Neotropics, and provide the basis for future evolutionary and taxonomic studies.

Identifying multiple origins of polyploid taxa: a multilocus study of the hybrid cloak fern (Astrolepis integerrima; Pteridaceae)

PREMISE OF THE STUDY

Molecular studies have shown that multiple origins of polyploid taxa are the rule rather than the exception. To understand the distribution and ecology of polyploid species and the evolutionary significance of polyploidy in general, it is important to delineate these independently derived lineages as accurately as possible. Although gene flow among polyploid lineages and backcrossing to their diploid parents often confound this process, such post origin gene flow is very infrequent in asexual polyploids. In this study, we estimate the number of independent origins of the apomictic allopolyploid fern Astrolepis integerrima, a morphologically heterogeneous species most common in the southwestern United States and Mexico, with outlying populations in the southeastern United States and the Caribbean.

METHODS

Plastid DNA sequence and AFLP data were obtained from 33 A. integerrima individuals. Phylogenetic analysis of the sequence data and multidimensional clustering of the AFLP data were used to identify independently derived lineages.

KEY RESULTS

Analysis of the two datasets identified 10 genetic groups within the 33 analyzed samples. These groups suggest a minimum of 10 origins of A. integerrima in the northern portion of its range, with both putative parents functioning as maternal donors, both supplying unreduced gametes, and both contributing a significant portion of their genetic diversity to the hybrids.

CONCLUSIONS

Our results highlight the extreme cryptic genetic diversity and systematic complexity that can underlie a single polyploid taxon.

Genetic diversity of Anopheles triannulatus s.l. (Diptera: Culicidae) from northwestern and southeastern Colombia

Anopheles triannulatus s.l. is a species complex, however in Colombia its taxonomic status is unclear. This study was conducted to understand the level of genetic differentiation or population structure of specimens of An. triannulatus s.l. from northwestern and southeastern Colombia. Cytochrome oxidase subunit I (COI) and internal transcribed spacer (ITS2) sequence analyses suggested high genetic differentiation between the NW and SE populations. A TCS network and Bayesian inference analysis based on 814 bp of COI showed two main groups: group I included samples from the NW and group II samples from the SE. Two main ITS2-polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) patterns were found. Pattern I is present in both the NW and SE, and pattern II is found in the SE specimens. To further elucidate the taxonomic status of An. triannulatus s.l. in Colombia and how these COI lineages are related to the Triannulatus Complex species, the evaluation of immature stages, male genitalia, and additional mitochondrial and nuclear markers will be needed.

Large-scale species delimitation method for hyperdiverse groups

Accelerating the description of biodiversity is a major challenge as extinction rates increase. Integrative taxonomy combining molecular, morphological, ecological and geographical data is seen as the best route to reliably identify species. Classic molluscan taxonomic methodology proposes primary species hypotheses (PSHs) based on shell morphology. However, in hyperdiverse groups, such as the molluscan family Turridae, where most of the species remain unknown and for which homoplasy and plasticity of morphological characters is common, shell-based PSHs can be arduous. A four-pronged approach was employed to generate robust species hypotheses of a 1000 specimen South-West Pacific Turridae data set in which: (i) analysis of COI DNA Barcode gene is coupled with (ii) species delimitation tools GMYC (General Mixed Yule Coalescence Method) and ABGD (Automatic Barcode Gap Discovery) to propose PSHs that are then (iii) visualized using Klee diagrams and (iv) evaluated with additional evidence, such as nuclear gene rRNA 28S, morphological characters, geographical and bathymetrical distribution to determine conclusive secondary species hypotheses (SSHs). The integrative taxonomy approach applied identified 87 Turridae species, more than doubling the amount previously known in the Gemmula genus. In contrast to a predominantly shell-based morphological approach, which over the last 30 years proposed only 13 new species names for the Turridae genus Gemmula, the integrative approach described here identified 27 novel species hypotheses not linked to available species names in the literature. The formalized strategy applied here outlines an effective and reproducible protocol for large-scale species delimitation of hyperdiverse groups.

Genetic population structure in the Antarctic benthos: insights from the widespread amphipod, Orchomenella franklini

Currently there is very limited understanding of genetic population structure in the Antarctic benthos. We conducted one of the first studies of microsatellite variation in an Antarctic benthic invertebrate, using the ubiquitous amphipod Orchomenella franklini (Walker, 1903). Seven microsatellite loci were used to assess genetic structure on three spatial scales: sites (100 s of metres), locations (1-10 kilometres) and regions (1000 s of kilometres) sampled in East Antarctica at Casey and Davis stations. Considerable genetic diversity was revealed, which varied between the two regions and also between polluted and unpolluted sites. Genetic differentiation among all populations was highly significant (F(ST) = 0.086, R(ST) = 0.139, p<0.001) consistent with the brooding mode of development in O. franklini. Hierarchical AMOVA revealed that the majority of the genetic subdivision occurred across the largest geographical scale, with N(e)m≈1 suggesting insufficient gene flow to prevent independent evolution of the two regions, i.e., Casey and Davis are effectively isolated. Isolation by distance was detected at smaller scales and indicates that gene flow in O. franklini occurs primarily through stepping-stone dispersal. Three of the microsatellite loci showed signs of selection, providing evidence that localised adaptation may occur within the Antarctic benthos. These results provide insights into processes of speciation in Antarctic brooders, and will help inform the design of spatial management initiatives recently endorsed for the Antarctic benthos.

Evaluating multiple criteria for species delimitation: an empirical example using Hawaiian palms (Arecaceae: Pritchardia)

BACKGROUND

Robust species delimitations are fundamental for conservation, evolutionary, and systematic studies, but they can be difficult to estimate, particularly in rapid and recent radiations. The consensus that species concepts aim to identify evolutionarily distinct lineages is clear, but the criteria used to distinguish evolutionary lineages differ based on the perceived importance of the various characteristics of evolving populations. We examined three different species-delimitation criteria (monophyly, absence of genetic intermediates, and diagnosability) to determine whether currently recognized species of Hawaiian Pritchardia are distinct lineages.

RESULTS

Data from plastid and nuclear genes, microsatellite loci, and morphological characters resulted in various levels of lineage subdivision that were likely caused by differing evolutionary rates between data sources. Additionally, taxonomic entities may be confounded because of the effects of incomplete lineage sorting and/or gene flow. A coalescent species tree was largely congruent with the simultaneous analysis, consistent with the idea that incomplete lineage sorting did not mislead our results. Furthermore, gene flow among populations of sympatric lineages likely explains the admixture and lack of resolution between those groups.

CONCLUSIONS

Delimiting Hawaiian Pritchardia species remains difficult but the ability to understand the influence of the evolutionary processes of incomplete lineage sorting and hybridization allow for mechanisms driving species diversity to be inferred. These processes likely extend to speciation in other Hawaiian angiosperm groups and the biota in general and must be explicitly accounted for in species delimitation.