Bayesian species delimitation in West African forest geckos (Hemidactylus fasciatus)

Artificial intelligence enables unified analysis of historical and landscape influences on genetic diversity

While genetic variation in any species is potentially shaped by a range of processes, phylogeography and landscape genetics are largely concerned with inferring how environmental conditions and landscape features impact neutral intraspecific diversity. However, even as both disciplines have come to utilize SNP data over the last decades, analytical approaches have remained for the most part focused on either broad-scale inferences of historical processes (phylogeography) or on more localized inferences about environmental and/or landscape features (landscape genetics). Here we demonstrate that an artificial intelligence model-based analytical framework can consider both deeper historical factors and landscape-level processes in an integrated analysis. We implement this framework using data collected from two Brazilian anurans, the Brazilian sibilator frog (Leptodactylus troglodytes) and granular toad (Rhinella granulosa). Our results indicate that historical demographic processes shape most the genetic variation in the sibulator frog, while landscape processes primarily influence variation in the granular toad. The machine learning framework used here allows both historical and landscape processes to be considered equally, rather than requiring researchers to make an a priori decision about which factors are important.

Divergence and reticulation in the Mexican white oaks: ecological and phylogenomic evidence on species limits and phylogenetic networks in the Quercus laeta complex (Fagaceae)

BACKGROUND AND AIMS

Introgressive hybridization poses a challenge to taxonomic and phylogenetic understanding of taxa, particularly when there are high numbers of co-occurring, intercrossable species. The genus Quercus exemplifies this situation. Oaks are highly diverse in sympatry and cross freely, creating syngameons of interfertile species. Although a well-resolved, dated phylogeny is available for the American oak clade, evolutionary relationships within many of the more recently derived clades remain to be defined, particularly for the young and exceptionally diverse Mexican white oak clade. Here, we adopted an approach bridging micro- and macroevolutionary scales to resolve evolutionary relationships in a rapidly diversifying clade endemic to Mexico.

METHODS

Ecological data and sequences of 155 low-copy nuclear genes were used to identify distinct lineages within the Quercus laeta complex. Concatenated and coalescent approaches were used to assess the phylogenetic placement of these lineages relative to the Mexican white oak clade. Phylogenetic network methods were applied to evaluate the timing and genomic significance of recent or historical introgression among lineages.

KEY RESULTS

The Q. laeta complex comprises six well-supported lineages, each restricted geographically and with mostly divergent climatic niches. Species trees corroborated that the different lineages are more closely related to other species of Mexican white oaks than to each other, suggesting that this complex is polyphyletic. Phylogenetic networks estimated events of ancient introgression that involved the ancestors of three present-day Q. laeta lineages.

CONCLUSIONS

The Q. laeta complex is a morphologically and ecologically related group of species rather than a clade. Currently, oak phylogenetics is at a turning point, at which it is necessary to integrate phylogenetics and ecology in broad regional samples to figure out species boundaries. Our study illuminates one of the more complicated of the Mexican white oak groups and lays groundwork for further taxonomic study.

Morphological Species Delimitation in The Western Pond Turtle (Actinemys): Can Machine Learning Methods Aid in Cryptic Species Identification?

As the discovery of cryptic species has increased in frequency, there has been an interest in whether geometric morphometric data can detect fine-scale patterns of variation that can be used to morphologically diagnose such species. We used a combination of geometric morphometric data and an ensemble of five supervised machine learning methods (MLMs) to investigate whether plastron shape can differentiate two putative cryptic turtle species, Actinemys marmorata and Actinemys pallida. Actinemys has been the focus of considerable research due to its biogeographic distribution and conservation status. Despite this work, reliable morphological diagnoses for its two species are still lacking. We validated our approach on two datasets, one consisting of eight morphologically disparate emydid species, the other consisting of two subspecies of Trachemys (T. scripta scripta, T. scripta elegans). The validation tests returned near-perfect classification rates, demonstrating that plastron shape is an effective means for distinguishing taxonomic groups of emydids via MLMs. In contrast, the same methods did not return high classification rates for a set of alternative phylogeographic and morphological binning schemes in Actinemys. All classification hypotheses performed poorly relative to the validation datasets and no single hypothesis was unequivocally supported for Actinemys. Two hypotheses had machine learning performance that was marginally better than our remaining hypotheses. In both cases, those hypotheses favored a two-species split between A. marmorata and A. pallida specimens, lending tentative morphological support to the hypothesis of two Actinemys species. However, the machine learning results also underscore that Actinemys as a whole has lower levels of plastral variation than other turtles within Emydidae, but the reason for this morphological conservatism is unclear.

Molecular and morphometric study of Brazilian populations of Psychodopygus davisi

In this study, we analysed the molecular and morphometric differences of several populations of the putative sand fly vector Psychodopygus davisi (Root, 1934) (Diptera, Psychodidae, Phlebotominae) in Brazil. We amplified the 658 base pair fragments of the DNA barcoding region-cytochrome c oxidase subunit 1 (COI) gene-for 57 specimens of P. davisi and three specimens of Psychodopygus claustrei (Abonnenc, Léger & Fauran, 1979). We merged our data with public sequences of the same species available from GenBank. Then, the combined dataset-87 sequences and 20 localities-was analysed using population structure analysis and different species delimitation approaches. Geometric morphometry of wings was performed for 155 specimens of P. davisi populations from the North, Midwest and Southeast Brazilian regions, analysing the differences in centroid sizes and canonical variates. Molecular analysis indicated high intraspecific genetic distance values for P. davisi (maximum p distance = 5.52%). All algorithms identified P. davisi and P. claustrei as distinct molecular taxonomic units, despite the low interspecific distance (p distance to the nearest neighbour = 4.79%). P. davisi sequences were split into four genetic clusters by population structure analysis and at least five genetic lineages using intermediate scenarios of the species delimitation algorithms. The species validation analysis of BPP strongly supported the five-species model in our dataset. We found high genetic diversity in this taxon, which is in agreement with its wide geographic distribution in Brazil. Furthermore, the wing analysis showed that specimens from the Southeast Region of Brazil are different from those in the North and the Midwest. The evolutionary patterns of P. davisi populations in Brazil suggest the presence of candidate species, which need to be validated in future studies using a more comprehensive approach with both genomic data and morphological characters.

Species Diagnosis and DNA Taxonomy

The use of DNA has helped to improve and speed up species identification and delimitation. However, it also provides new challenges to taxonomists. Incongruence of outcome from various markers and delimitation methods, bias from sampling and skewed species distribution, implemented models, and the choice of methods/priors may mislead results and also may, in conclusion, increase elements of subjectivity in species taxonomy. The lack of direct diagnostic outcome from most contemporary molecular delimitation approaches and the need for a reference to existing and best sampled trait reference systems reveal the need for refining the criteria of species diagnosis and diagnosability in the current framework of nomenclature codes and good practices to avoid nomenclatorial instability, parallel taxonomies, and consequently more and new taxonomic impediment.

Phylogenetic review of the comb-tooth blenny genus Hypleurochilus in the northwest Atlantic and Gulf of Mexico

As some of the smallest vertebrates, yet largest producers of consumed reef biomass, cryptobenthic reef fishes serve a disproportionate role in reef ecosystems and are one of the most poorly understood groups of fish. The blenny genera Hypleurochilus and Parablennius are currently considered paraphyletic and the interrelationships of Parablennius have been the focus of recent phylogenetic studies. However, the interrelationships of Hypleurochilus remain understudied. This genus is transatlantically distributed and comprises 11 species with a convoluted taxonomic history. In this study, relationships for ten Hypleurochilus species are resolved using multi-locus nuclear and mtDNA sequence data, morphological data, and mined COI barcode data.  Mitochondrial and nuclear sequence data from 61 individuals collected from the western Atlantic and northern Gulf of Mexico (N. GoM) delimit seven species into a temperate clade, a tropical clade, and a third distinct lineage. This lineage, herein referred to as H. cf. aequipinnis, may represent a species of Hypleurochilus whose range has expanded into the N. GoM. Inclusion of publicly available COI sequence for an additional three species provides further phylogenetic resolution. H. bananensis forms a new eastern Atlantic clade with H. cf. aequipinnis, providing further evidence for a western Atlantic range expansion. Single marker COI delimitation was unable to elucidate the relationships between H. springeri/H. pseudoaequipinnis and between H. multifilis/H. caudovittatus due to incomplete lineage sorting. Mitochondrial data are also unable to accurately resolve the placement of H. bermudensis. However, a comprehensive approach using multi-locus phylogenetic and species delimitation methods was able to resolve these relationships. While mining publicly available sequence data allowed for the inclusion of an increased number of species in the analysis and a more comprehensive phylogeny, it was not without drawbacks, as a handful of sequences are potentially mis-identified. Overall, we find that the recent divergence of some species within this genus and potential introgression events confound the results of single locus delimitation methods, yet a combination of single and multi-locus analyses has allowed for insights into the biogeography of this genus and uncovered a potential transatlantic range expansion.

Unsupervised machine learning for species delimitation, integrative taxonomy, and biodiversity conservation

Integrative taxonomy, combining data from multiple axes of biologically relevant variation, is a major goal of systematics. Ideally, such taxonomies will derive from similarly integrative species-delimitation analyses. Yet, most current methods rely solely or primarily on molecular data, with other layers often incorporated only in a post hoc qualitative or comparative manner. A major limitation is the difficulty of devising quantitative parametric models linking different datasets in a unified ecological and evolutionary framework. Machine Learning (ML) methods offer flexibility in this arena by easily learning high-dimensional associations between observations (e.g., individual specimens) across a wide array of input features (e.g., genetics, geography, environment, and phenotype) to delimit statistically meaningful clusters. Here, I implement an unsupervised method using Self-Organizing (or "Kohonen") Maps (SOMs) for such purposes. Recent extensions called "SuperSOMs" can integrate multiple layers, each of which exerts independent influence on a two-dimensional output grid via empirically estimated weights. The grid cells are then delimited into K distinct units that can be interpreted as species or other entities. I show empirical examples in salamanders (Desmognathus) and snakes (Storeria) with layers representing alleles, space, climate, and traits. Simulations reveal that the SuperSOM approach can detect K = 1, tends not to over-split, reflects contributions from all layers, and limits large layers (e.g., genetic matrices) from overwhelming other datasets, desirable properties addressing major concerns from previous studies. Finally, I suggest that these and similar methods could integrate conservation-relevant layers such as population trends and human encroachment to delimit management units from an explicitly quantitative framework grounded in the ecology and evolution of species limits and boundaries.

Genomic analysis reveals a cryptic pangolin species

Eight extant species of pangolins are currently recognized. Recent studies found that two mitochondrial haplotypes identified in confiscations in Hong Kong could not be assigned to any known pangolin species, implying the existence of a species. Here, we report that two additional mitochondrial haplotypes identified in independent confiscations from Yunnan align with the putative species haplotypes supporting the existence of this mysterious species/population. To verify the new species scenario we performed a comprehensive analysis of scale characteristics and 138 whole genomes representing all recognized pangolin species and the cryptic new species, 98 of which were generated here. Our morphometric results clearly attributed this cryptic species to Asian pangolins (Manis sp.) and the genomic data provide robust and compelling evidence that it is a pangolin species distinct from those recognized previously, which separated from the Philippine pangolin and Malayan pangolin over 5 Mya. Our study provides a solid genomic basis for its formal recognition as the ninth pangolin species or the fifth Asian one, supporting a new taxonomic classification of pangolins. The effects of glacial climate changes and recent anthropogenic activities driven by illegal trade are inferred to have caused its population decline with the genomic signatures showing low genetic diversity, a high level of inbreeding, and high genetic load. Our finding greatly expands current knowledge of pangolin diversity and evolution and has vital implications for conservation efforts to prevent the extinction of this enigmatic and endangered species from the wild.

Multi-locus phylogeny of the catfish genus Ictalurus Rafinesque, 1820 (Actinopterygii, Siluriformes) and its systematic and evolutionary implications

BACKGROUND

Ictalurus is one of the most representative groups of North American freshwater fishes. Although this group has a well-studied fossil record and has been the subject of several morphological and molecular phylogenetic studies, incomplete taxonomic sampling and insufficient taxonomic studies have produced a rather complex classification, along with intricate patterns of evolutionary history in the genus that are considered unresolved and remain under debate.

RESULTS

Based on four loci and the most comprehensive taxonomic sampling analyzed to date, including currently recognized species, previously synonymized species, undescribed taxa, and poorly studied populations, this study produced a resolved phylogenetic framework that provided plausible species delimitation and an evolutionary time framework for the genus Ictalurus.

CONCLUSIONS

Our phylogenetic hypothesis revealed that Ictalurus comprises at least 13 evolutionary units, partially corroborating the current classification and identifying populations that emerge as putative undescribed taxa. The divergence times of the species indicate that the diversification of Ictalurus dates to the early Oligocene, confirming its status as one of the oldest genera within the family Ictaluridae.

Diversity of Pholcus Spiders (Araneae: Pholcidae) in China's Lüliang Mountains: An Integrated Morphological and Molecular Approach

Spiders of the genus Pholcus were collected for the first time during an expedition to the Lüliang Mountains in Shanxi Province, North China. Phylogenetic analyses of DNA sequence data from COI, H3, wnt, and 28S genes allowed us to group them into nine well-supported clades. We used morphology and four methods of molecular species delimitation, namely Automatic Barcode Gap Discovery (ABGD), the Generalized Mixed Yule Coalescent (GMYC), Bayesian Poisson Tree Processes (bPTP), and Bayesian Phylogenetics and Phylogeography (BPP), to investigate species boundaries. These integrative taxonomic analyses identified the nine clades as nine distinct species, comprising Pholcus luya Peng & Zhang, 2013 and eight other species new to science: Pholcus jiaocheng sp. nov., Pholcus linfen sp. nov., Pholcus lishi sp. nov., Pholcus luliang sp. nov., Pholcus wenshui sp. nov., Pholcus xiangfen sp. nov., Pholcus xuanzhong sp. nov., and Pholcus zhongyang sp. nov. The species occur in geographic proximity and show many morphological similarities. All of them belong to the P. phungiformes species group. The records from the Lüliang Mountains represent the westernmost distribution limit of this species group.

A near-complete species-level phylogeny of uropeltid snakes harnessing historical museum collections as a DNA source

Uropeltidae is a clade of small fossorial snakes (ca. 64 extant species) endemic to peninsular India and Sri Lanka. Uropeltid taxonomy has been confusing, and the status of some species has not been revised for over a century. Attempts to revise uropeltid systematics and undertake evolutionary studies have been hampered by incompletely sampled and incompletely resolved phylogenies. To address this issue, we take advantage of historical museum collections, including type specimens, and apply genome-wide shotgun (GWS) sequencing, along with recent field sampling (using Sanger sequencing) to establish a near-complete multilocus species-level phylogeny (ca. 87% complete at species level). This results in a phylogeny that supports the monophyly of all genera (if Brachyophidium is considered a junior synonym of Teretrurus), and provides a firm platform for future taxonomic revision. Sri Lankan uropeltids are probably monophyletic, indicating a single colonisation event of this island from Indian ancestors. However, the position of Rhinophis goweri (endemic to Eastern Ghats, southern India) is unclear and warrants further investigation, and evidence that it may nest within the Sri Lankan radiation indicates a possible recolonisation event. DNA sequence data and morphology suggest that currently recognised uropeltid species diversity is substantially underestimated. Our study highlights the benefits of integrating museum collections in molecular genetic analyses and their role in understanding the systematics and evolutionary history of understudied organismal groups.

Coalescent-Based Species Delimitation in Herbaceous Bamboos (Bambusoideae, Olyreae) from Eastern Brazil: Implications for Taxonomy and Conservation in a Group with Weak Morphological Divergence Coupled with Low Genetic Diversity

Species delimitation in herbaceous bamboos has been complex and, in some genera, a great part of its diversity has been confirmed only based on genetic information, as is the case of the genus Raddia. It includes nine species, all occurring in Brazil, but only R. portoi predominates in dry forests of the Northeast associated with the Caatinga phytogeographic domain. This species is morphologically close to R. angustifolia, which is known for a single location in the Atlantic Forest in Southern Bahia, and is considered to be threatened by extinction. Besides problems with taxonomic focus, actions for its conservation are complicated because it is not certain if it must be considered an independent species or included in the more widespread R. portoi. In this study, we used coalescent multispecies (MSC) theory approaches combined with genetic structure analyses in an attempt to delimit these two species. Different analyses were congruent and the species delimitation using MSC inferred distinct lineages supporting their recognition as two species. These results solved the taxonomic doubts and also showed the power of these approaches to delimit species as lineages, even in groups with weak morphological divergence and low genetic variability, and also impacting our knowledge for conservation purposes.

Species boundary and phylogeographical pattern provide new insights into the management efforts of Megacopta cribraria (Hemiptera: Plataspidae), a bean bug invading North America

BACKGROUND

Correct identification together with information on distribution range, geographical origin and evolutionary history are the necessary basis for the management and control of invasive species. The bean bug Megacopta cribraria is a crucial agricultural pest of soybean. Recently, M. cribraria has invaded the United States and spread rapidly, causing severe reductions in soybean yields. However, the species boundary and phylogeographical pattern of this invasive bean bug are still unclear.

RESULTS

The results of different species delimitation methods (Automatic Barcode Gap Discovery, Assemble Species by Automatic Partitioning, Bayesian Poisson Tree Processes and Bayesian Phylogenetics and Phylogeography) strongly demonstrated that M. cribraria and Megacopta punctatissima represent the same species. M. punctatissima should not be considered a distinct species but rather a variety of M. cribraria. Phylogenetic analyses revealed three well-supported clades (Southeast Asia [SEA], East Asia continent [EAC] and Japan [JA]) with distinct geographical structures in the M. cribraria-M. punctatissima complex. The SEA clade was at the base of the phylogenetic tree, and the sister relationship between the EAC clade and JA clade was strongly supported. The split between the EAC clade and JA clade occurred at approximately 0.71 Ma, corresponding to the submergence period of the East China Sea land bridge.

CONCLUSION

This study clarified the species boundary between M. cribraria and its closely related species and revealed the phylogeographical pattern and evolutionary history of M. cribraria. The species delimitation and phylogeography results achieved in this study could provide new insights into the monitoring and management of this agricultural pest. © 2022 Society of Chemical Industry.

Genome-wide species delimitation analyses of a silverside fish species complex in central Mexico indicate taxonomic over-splitting

BACKGROUND

Delimiting species across a speciation continuum is a complex task, as the process of species origin is not generally instantaneous. The use of genome-wide data provides unprecedented resolution to address convoluted species delimitation cases, often unraveling cryptic diversity. However, because genome-wide approaches based on the multispecies coalescent model are known to confound population structure with species boundaries, often resulting in taxonomic over-splitting, it has become increasingly evident that species delimitation research must consider multiple lines of evidence. In this study, we used phylogenomic, population genomic, and coalescent-based species delimitation approaches, and examined those in light of morphological and ecological information, to investigate species numbers and boundaries comprising the Chirostoma "humboltianum group" (family Atherinidae). The humboltianum group is a taxonomically controversial species complex where previous morphological and mitochondrial studies produced conflicting species delimitation outcomes. We generated ddRADseq data for 77 individuals representing the nine nominal species in the group, spanning their distribution range in the central Mexican plateau.

RESULTS

Our results conflict with the morphospecies and ecological delimitation hypotheses, identifying four independently evolving lineages organized in three geographically cohesive clades: (i) chapalae and sphyraena groups in Lake Chapala, (ii) estor group in Lakes Pátzcuaro and Zirahuén, and (iii) humboltianum sensu stricto group in Lake Zacapu and Lerma river system.

CONCLUSIONS

Overall, our study provides an atypical example where genome-wide analyses delineate fewer species than previously recognized on the basis of morphology. It also highlights the influence of the geological history of the Chapala-Lerma hydrological system in driving allopatric speciation in the humboltianum group.

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.

Genomes From Historic DNA Unveil Massive Hidden Extinction and Terminal Endangerment in a Tropical Asian Songbird Radiation

Quantifying the magnitude of the global extinction crisis is important but remains challenging, as many extinction events pass unnoticed owing to our limited taxonomic knowledge of the world's organisms. The increasing rarity of many taxa renders comprehensive sampling difficult, further compounding the problem. Vertebrate lineages such as birds, which are thought to be taxonomically well understood, are therefore used as indicator groups for mapping and quantifying global extinction. To test whether extinction patterns are adequately gauged in well-studied groups, we implemented ancient-DNA protocols and retrieved whole genomes from the historic DNA of museum specimens in a widely known songbird radiation of shamas (genus Copsychus) that is assumed to be of least conservation concern. We uncovered cryptic diversity and an unexpected degree of hidden extinction and terminal endangerment. Our analyses reveal that >40% of the phylogenetic diversity of this radiation is already either extinct in the wild or nearly so, including the two genomically most distinct members of this group (omissus and nigricauda), which have so far flown under the conservation radar as they have previously been considered subspecies. Comparing the genomes of modern samples with those from roughly a century ago, we also found a significant decrease in genetic diversity and a concomitant increase in homozygosity affecting various taxa, including small-island endemics that are extinct in the wild as well as subspecies that remain widespread across the continental scale. Our application of modern genomic approaches demonstrates elevated levels of allelic and taxonomic diversity loss in a songbird clade that has not been listed as globally threatened, highlighting the importance of ongoing reassessments of extinction incidence even across well-studied animal groups. Key words: extinction, introgression, white-rumped shama, conservation.

The Importance of Contact Zones for Distinguishing Interspecific from Intraspecific Geographic Variation

With limited sampling, geographic variation within a single species can be difficult to distinguish from interspecific variation, confounding our ability to draw accurate species boundaries. We argue that thorough sampling and analysis of contact zones between putative taxa can determine if assortative mating or selection against hybrids exists (supporting the presence of two distinct species), or alternatively if mating is random among genotypes and admixture among adjacent populations is gradual and continuous (supporting geographic variation within a single species). Here, we test two alternative hypotheses for two pairs of named taxa at contact zones within the American milksnake (Lampropeltis triangulum) complex. A prior morphological analysis found areas of gradual intergradation among named taxa, and concluded that the taxa represented geographical races of a single polytypic species. In contrast, a subsequent analysis of gene sequence data, but with limited sampling near the contact zones, hypothesized distinct boundaries between species at the contact zones. At the contact zone between proposed species L. triangulum and L. gentilis, we examined a ∼700 km-wide transect across the states of Kansas and Missouri, with thorough sampling and reduced-representation genomic-level sequencing, to test the two opposing taxonomic hypotheses. Our transect analyses included examinations of population structure, fixed differences, cline-fitting, and an admixture index analysis. These analyses all supported a gradual and continuous geographic cline across a broad intergrade zone between two geographic forms of L. triangulum, thus providing strong support for a single species in this region (and no support for the recognition of L. gentilis as a distinct species). At a second contact zone between proposed species L. triangulum and L. elapsoides (but variously treated as species or subspecies by different researchers) in Kentucky and Tennessee, we re-evaluated morphological data. In this case, the contact zone analysis indicated sympatry and reproductive isolation of the two taxa, and thus strongly supported L. triangulum and L. elapsoides as distinct species. We conclude that detailed studies of contact zones, based on either genetic or morphological data, are essential for distinguishing intraspecific from interspecific variation in the case of widely and continuously distributed taxa.

Genetic variability and the ecology of geographic range: A test of the central-marginal hypothesis in Australian scincid lizards

For many species, both local abundance and regional occupancy are highest near the centre of their geographic distributions. One hypothesis for this pattern is that niche suitability declines with increasing distance from a species geographic centre, such that populations near range margins are characterized by reduced density and increased patchiness. In these smaller edge populations, genetic drift is more powerful, leading to the loss of genetic diversity. This simple verbal model has been formalized as the central-marginal hypothesis, which predicts that core populations should have greater genetic diversity than edge populations. Here, we tested the central-marginal hypothesis using a genomic data set of 25 species-level taxa of Australian scincid lizards in the genera Ctenotus and Lerista. A majority of taxa in our data set showed range-wide patterns of genetic variation consistent with central-marginal hypothesis, and eight of 25 taxa showed significantly greater genetic diversity in the centre of their range. We then explored biological, historical, and methodological factors that might predict which taxa support the central-marginal hypothesis. We found that taxa with the strongest evidence for range expansion were the least likely to follow predictions of the central-marginal hypothesis. The majority of these taxa had range expansions that originated at the range edge, which led to a gradient of decreasing genetic diversity from the range edge to the core, contrary to the central-marginal hypothesis.

Predicting Species Boundaries and Assessing Undescribed Diversity in Pneumocystis, an Obligate Lung Symbiont

Far more biodiversity exists in Fungi than has been described, or could be described in several lifetimes, given current rates of species discovery. Although this problem is widespread taxonomically, our knowledge of animal-associated fungi is especially lacking. Fungi in the genus Pneumocystis are obligate inhabitants of mammal lungs, and they have been detected in a phylogenetically diverse array of species representing many major mammal lineages. The hypothesis that Pneumocystis cospeciate with their mammalian hosts suggests that thousands of Pneumocystis species may exist, potentially equal to the number of mammal species. However, only six species have been described, and the true correspondence of Pneumocystis diversity to host species boundaries is unclear. Here, we use molecular species delimitation to estimate the boundaries of Pneumocystis species sampled from 55 mammal species representing eight orders. Our results suggest that Pneumocystis species often colonize several closely related mammals, especially those in the same genus. Using the newly estimated ratio of fungal to host diversity, we estimate ≈4600 to 6250 Pneumocystis species inhabit the 6495 currently recognized extant mammal species. Additionally, we review the literature and find that only 240 (~3.7%) mammal species have been screened for Pneumocystis, and many detected Pneumocystis lineages are not represented by any genetic data. Although crude, our findings challenge the dominant perspective of strict specificity of Pneumocystis to their mammal hosts and highlight an abundance of undescribed diversity.

Quantitatively defining species boundaries with more efficiency and more biological realism

We introduce a widely applicable species delimitation method based on the multispecies coalescent model that is more efficient and more biologically realistic than existing methods. We extend a threshold-based method to allow the ancestral speciation rate to vary through time as a smooth piecewise function. Furthermore, we introduce the cutting-edge proposal kernels of StarBeast3 to this model, thus enabling rapid species delimitation on large molecular datasets and allowing the use of relaxed molecular clock models. We validate these methods with genomic sequence data and SNP data, and show they are more efficient than existing methods at achieving parameter convergence during Bayesian MCMC. Lastly, we apply these methods to two datasets (Hemidactylus and Galagidae) and find inconsistencies with the published literature. Our methods are powerful for rapid quantitative testing of species boundaries in large multilocus datasets and are implemented as an open source BEAST 2 package called SPEEDEMON.

Introducing SPEEDEMON, a package for BEAST 2 that better defines species boundaries based on molecular data demonstrated on gecko and loris datasets.

Cryptic species in a colorful genus: Integrative taxonomy of the bush robins (Aves, Muscicapidae, Tarsiger) suggests two overlooked species

Several cryptic avian species have been validated by recent integrative taxonomic efforts in the Sino-Himalayan mountains, indicating that avian diversity in this global biodiversity hotspot may be underestimated. In the present study, we investigated species limits in the genus Tarsiger, the bush robins, a group of montane forest specialists with high species richness in the Sino-Himalayan region. Based on comprehensive sampling of all 11 subspecies of the six currently recognized species, we applied an integrative taxonomic approach by combining multilocus, acoustic, plumage and morphometric analyses. Our results reveal that the isolated north-central Chinese populations of Tarsiger cyanurus, described as the subspecies albocoeruleus but usually considered invalid, is distinctive in genetics and vocalisation, but only marginally differentiated in morphology. We also found the Taiwan endemic T. indicus formosanus to be distinctive in genetics, song and morphology from T. i. indicus and T. i. yunnanensis of the Sino-Himalayan mountains. Moreover, Bayesian species delimitation using BPP suggests that both albocoeruleus and formosanus merit full species status. We propose their treatment as 'Qilian Bluetail' T. albocoeruleus and 'Taiwan Bush Robin' T. formosanus, respectively.

Multilocus Data Analysis Reveal the Diversity of Cryptic Species in the Tillandsia ionantha (Bromeliaceae: Tillansiodeae) Complex

Independent evolutionary lineages or species that lack phenotypic variation as an operative criterion for their delimitation are known as cryptic species. However, these have been delimited using other data sources and analysis. The aims of this study are: (1) to evaluate the divergence of the populations of the T. ionantha complex; and (2) to delimit the species using multilocus data, phylogenetic analysis and the coalescent model. Phylogenetic analyses, genetic diversity and population structure, and isolation by distance analysis were performed. A multispecies coalescent analysis to delimit the species was conducted. Phylogenetic analysis showed that T. ionantha is polyphyletic composed of eight evolutionary lineages. Haplotype distribution and genetic differentiation analysis detected strong population structure and high values of genetic differentiation among populations. The positive correlation between genetic differences with geographic distance indicate that the populations are evolving under the model of isolation by distance. The coalescent multispecies analysis performed with starBEAST supports the recognition of eight lineages as different species. Only three out of the eight species have morphological characters good enough to recognize them as different species, while five of them are cryptic species. Tillandsia scaposa and T. vanhyningii are corroborated as independent lineages, and T. ionantha var. stricta changed status to the species level.

Genome and cuticular hydrocarbon-based species delimitation shed light on potential drivers of speciation in a Neotropical ant species complex

Geographic separation that leads to the evolution of reproductive isolation between populations generally is considered the most common form of speciation. However, speciation may also occur in the absence of geographic barriers due to phenotypic and genotypic factors such as chemical cue divergence, mating signal divergence, and mitonuclear conflict. Here, we performed an integrative study based on two genome-wide techniques (3RAD and ultraconserved elements) coupled with cuticular hydrocarbon (CHC) and mitochondrial (mt) DNA sequence data, to assess the species limits within the Ectatomma ruidum species complex, a widespread and conspicuous group of Neotropical ants for which heteroplasmy (i.e., presence of multiple mtDNA variants in an individual) has been recently discovered in some populations from southeast Mexico. Our analyses indicate the existence of at least five distinct species in this complex: two widely distributed across the Neotropics, and three that are restricted to southeast Mexico and that apparently have high levels of heteroplasmy. We found that species boundaries in the complex did not coincide with geographic barriers. We therefore consider possible roles of alternative drivers that may have promoted the observed patterns of speciation, including mitonuclear incompatibility, CHC differentiation, and colony structure. Our study highlights the importance of simultaneously assessing different sources of evidence to disentangle the species limits of taxa with complicated evolutionary histories.

No link between population isolation and speciation rate in squamate reptiles

Rates of species formation vary widely across the tree of life and contribute to massive disparities in species richness among clades. This variation can emerge from differences in metapopulation-level processes that affect the rates at which lineages diverge, persist, and evolve reproductive barriers and ecological differentiation. For example, populations that evolve reproductive barriers quickly should form new species at faster rates than populations that acquire reproductive barriers more slowly. This expectation implicitly links microevolutionary processes (the evolution of populations) and macroevolutionary patterns (the profound disparity in speciation rate across taxa). Here, leveraging extensive field sampling from the Neotropical Cerrado biome in a biogeographically controlled natural experiment, we test the role of an important microevolutionary process-the propensity for population isolation-as a control on speciation rate in lizards and snakes. By quantifying population genomic structure across a set of codistributed taxa with extensive and phylogenetically independent variation in speciation rate, we show that broad-scale patterns of species formation are decoupled from demographic and genetic processes that promote the formation of population isolates. Population isolation is likely a critical stage of speciation for many taxa, but our results suggest that interspecific variability in the propensity for isolation has little influence on speciation rates. These results suggest that other stages of speciation-including the rate at which reproductive barriers evolve and the extent to which newly formed populations persist-are likely to play a larger role than population isolation in controlling speciation rate variation in squamates.

Ecological divergence and synchronous Pleistocene diversification in the widespread South American butter frog complex

Phylogeographic studies primarily focus on the major role of landscape topography in driving lineage diversification. However, populational phylogeographic breaks may also occur as a result of either niche conservatism or divergence, in the absence of geographic barriers to gene flow. Furthermore, these two factors are not mutually exclusive and can act in concert, making it challenging to evaluate their relative importance on explaining genetic variation in nature. Herein, we use sequences of two mitochondrial and four nuclear genes to investigate the timing and diversification patterns of species pertaining to the Leptodactylus latrans complex, which harbors four morphologically cryptic species with broad distributions across environmental gradients in eastern South America. The origin of this species complex dates back to the late Miocene (ca. 5.5 Mya), but most diversification events occurred synchronically during the late Pleistocene likely as the result of ecological divergence driven by Quaternary climatic oscillations. Further, significant patterns of environmental niche divergences among species in the L. latrans complex imply that ecological isolation is the primary mode of genetic diversification, mostly because phylogenetic breaks are associated with environmental transitions rather than topographic barriers at both species and populational scale. We provided new insights about diversification patterns and processes within a species complex of broadly and continuously distributed group of frogs along South America.

Phylogenetic analysis of the Neotropical Albitarsis Complex based on mitogenome data

BACKGROUND

Some of the most important malaria vectors in South America belong to the Albitarsis Complex (Culicidae; Anophelinae; Anopheles). Understanding the origin, nature, and geographical distribution of species diversity in this important complex has important implications for vector incrimination, control, and management, and for modelling future responses to climate change, deforestation, and human population expansion. This study attempts to further explore species diversity and evolutionary history in the Albitarsis Complex by undertaking a characterization and phylogenetic analysis of the mitogenome of all 10 putative taxa in the Albitarsis Complex.

METHODS

Mitogenome assembly and annotation allowed for feature comparison among Albitarsis Complex and Anopheles species. Selection analysis was conducted across all 13 protein-coding genes. Maximum likelihood and Bayesian inference methods were used to construct gene and species trees, respectively. Bayesian methods were also used to jointly estimate species delimitation and species trees.

RESULTS

Gene composition and order were conserved across species within the complex. Unique signatures of positive selection were detected in two species-Anopheles janconnae and An. albitarsis G-which may have played a role in the recent and rapid diversification of the complex. The COI gene phylogeny does not fully recover the mitogenome phylogeny, and a multispecies coalescent-based phylogeny shows that considerable uncertainty exists through much of the mitogenome species tree. The origin of divergence in the complex dates to the Pliocene/Pleistocene boundary, and divergence within the distinct northern South American clade is estimated at approximately 1 million years ago. Neither the phylogenetic trees nor the delimitation approach rejected the 10-species hypothesis, although the analyses could not exclude the possibility that four putative species with scant a priori support (An. albitarsis G, An. albitarsis H, An. albitarsis I, and An. albitarsis J), represent population-level, rather than species-level, splits.

CONCLUSION

The lack of resolution in much of the species tree and the limitations of the delimitation analysis warrant future studies on the complex using genome-wide data and the inclusion of additional specimens, particularly from two putative species, An. albitarsis I and An. albitarsis J.

Incipient speciation and its impact on taxonomic decision: a case study using a sky island sister-species pair of stag beetles (Lucanidae: Lucanus)

Species delimitation can be difficult when the divergence between focal taxa is in the incipient stage of speciation, because conflicting results are expected among different data sets, and the species limits can differ depending on the species concept applied. We studied speciation history and investigated the impact on taxonomic decision-making when using different types of data in a Taiwanese endemic sister-species pair of stag beetles, Lucanus miwai and Lucanus yulaoensis, from sky island habitats. We showed that the two geographical taxa can be diagnosed by male mandibular shape. We found two mitochondrial co1 lineages with pairwise sequence divergence > 3%; however, L. miwai might not be monophyletic. The result of our multispecies coalescent-based species delimitation using five nuclear loci supported the evolutionary independence of the two sister species, but the calculated values of the genealogical divergence index (gdi) corresponded to the ambiguous zone of species delimitation. We also showed that post-divergence gene flow is unlikely. Our study demonstrates challenges in the delineation of incipient species, but shows the importance of understanding the speciation history and adopting integrative approaches to reconcile seemingly conflicting results before making evolutionarily relevant taxonomic decisions.

Molecular species delimitation of shrub frogs of the genus Pseudophilautus (Anura, Rhacophoridae)

Sri Lanka is an amphibian hotspot of global significance. Its anuran fauna is dominated by the shrub frogs of the genus Pseudophilautus. Except for one small clade of four species in Peninsular India, these cool-wet adapted frogs, numbering some 59 extant species, are distributed mainly across the montane and lowland rain forests of the island. With species described primarily by morphological means, the diversification has never yet been subjected to a molecular species delimitation analysis, a procedure now routinely applied in taxonomy. Here we test the species boundaries of Pseudophilautus in the context of the phylogenetic species concept (PSC). We use all the putative species for which credible molecular data are available (nDNA-Rag-1; mt-DNA- 12S rRNA, 16S rRNA) to build a well resolved phylogeny, which is subjected to species delimitation analyses. The ABGD, bPTP, mPTP and bGMYC species delimitation methods applied to the 16S rRNA frog barcoding gene (for all species), 12S rRNA and Rag-1 nDNA grouped P. procax and P. abundus; P. hallidayi and P. fergusonianus; P. reticulatus and P. pappilosus; P. pleurotaenia and P. hoipolloi; P. hoffmani and P. asankai; P. silvaticus and P. limbus; P. dilmah and P. hankeni; P. fulvus and P. silus.. Surprisingly, all analyses recovered 14 unidentified potential new species as well. The geophylogeny affirms a distribution across the island's aseasonal 'wet zone' and its three principal hill ranges, suggestive of allopatric speciation playing a dominant role, especially between mountain masses. Among the species that are merged by the delimitation analyses, a pattern leading towards a model of parapatric speciation emerges-ongoing speciation in the presence of gene flow. This delimitation analysis reinforces the species hypotheses, paving the way to a reasonable understanding of Sri Lankan Pseudophilautus, enabling both deeper analyses and conservation efforts of this remarkable diversification. http://zoobank.org/urn:lsid:zoobank.org:pub:DA869B6B-870A-4ED3-BF5D-5AA3F69DDD27.

Type specimen sequencing, multilocus analyses, and species delimitation methods recognize the cosmopolitan Corallina berteroi and establish the northern Japanese C. yendoi sp. nov. (Corallinaceae, Rhodophyta)

A partial rbcL sequence of the lectotype specimen of Corallina berteroi shows that it is the earliest available name for C. ferreyrae. Multilocus species delimitation analyses (ABGD, SPN, GMYC, bPTP, and BPP) using independent or concatenated COI, psbA, and rbcL sequences recognized one, two, or three species in this complex, but only with weak support for each species hypothesis. Conservatively, we recognize a single worldwide species in this complex of what appears to be multiple, evolving populations. Included in this species, besides C. ferreyrae, are C. caespitosa, the morphologically distinct C. melobesioides, and, based on a partial rbcL sequence of the holotype specimen, C. pinnatifolia. Corallina berteroi, not C. officinalis, is the cosmopolitan temperate species found thus far in the NE Atlantic, Mediterranean Sea, warm temperate NW Atlantic and NE Pacific, cold temperate SW Atlantic (Falkland Islands), cold and warm temperate SE Pacific, NW Pacific and southern Australia. Also proposed is C. yendoi sp. nov. from Hokkaido, Japan, which was recognized as distinct by 10 of the 13 species discrimination analyses, including the multilocus BPP.

Solving the Coral Species Delimitation Conundrum

Distinguishing coral species is not only crucial for physiological, ecological and evolutionary studies, but also to enable effective management of threatened reef ecosystems. However, traditional hypotheses that delineate coral species based on morphological traits from the coral skeleton are frequently at odds with tree-based molecular approaches. Additionally, a dearth of species-level molecular markers has made species delimitation particularly challenging in species-rich coral genera, leading to the widespread assumption that inter-specific hybridization might be responsible for this apparent conundrum. Here, we used three lines of evidence - morphology, breeding trials and molecular approaches - to identify species boundaries in a group of ecologically important tabular Acropora corals. In contrast to previous studies, our morphological analyses yielded groups that were congruent with experimental crosses as well as with coalescent-based and allele sharing-based multilocus approaches to species delimitation. Our results suggest that species of the genus Acropora are reproductively isolated and independently evolving units that can be distinguished morphologically. These findings not only pave the way for a taxonomic revision of coral species, but also outline an approach that can provide a solid basis to address species delimitation and provide conservation support to a wide variety of keystone organisms.

Sympatric and independently evolving lineages in the Thoropa miliaris - T. taophora species complex (Anura: Cycloramphidae)

Species delimitation can be challenging and affected by subjectivity. Sibling lineages that occur in sympatry constitute good candidates for species delimitation regardless of the adopted species concept. The Thoropa miliaris + T. taophora species complex exhibits high genetic diversity distributed in several lineages that occur sympatrically in the southeastern Atlantic Forest of Brazil. We used 414 loci obtained by anchored hybrid enrichment to characterize genetic variation in the Thoropa miliaris species group (T. saxatilis, T megatympanum, T. miliaris, and T. taophora), combining assignment analyses with traditional and coalescent phylogeny reconstruction. We also investigated evolutionary independence in co-occurring lineages by estimating gene flow, and validated lineages under the multispecies coalescent. We recovered most previously described lineages as unique populations in assignment analyses; exceptions include two lineages within T. miliaris that are further substructured, and the merging of all T. taophora lineages. We found very low probabilities of gene flow between sympatric lineages, suggesting independent evolution. Species tree inferences and species delimitation yielded resolved relationships and indicate that all lineages constitute putative species that diverged during the Pliocene and Pleistocene, later than previously estimated.

Enlightening the black and white: species delimitation and UNITE species hypothesis testing in the Russula albonigra species complex

Russula albonigra is considered a well-known species, morphologically delimited by the context of the basidiomata blackening without intermediate reddening, and the menthol-cooling taste of the lamellae. It is supposed to have a broad ecological range and a large distribution area. A thorough molecular analysis based on four nuclear markers (ITS, LSU, RPB2 and TEF1-α) shows this traditional concept of R. albonigra s. lat. represents a species complex consisting of at least five European, three North American, and one Chinese species. Morphological study shows traditional characters used to delimit R. albonigra are not always reliable. Therefore, a new delimitation of the R. albonigra complex is proposed and a key to the described European species of R. subgen. Compactae is presented. A lectotype and an epitype are designated for R. albonigra and three new European species are described: R. ambusta, R. nigrifacta, and R. ustulata. Different thresholds of UNITE species hypotheses were tested against the taxonomic data. The distance threshold of 0.5% gives a perfect match to the phylogenetically defined species within the R. albonigra complex. Publicly available sequence data can contribute to species delimitation and increase our knowledge on ecology and distribution, but the pitfalls are short and low quality sequences.

Contrasting genetic signal of recolonization after rainforest fragmentation in African trees with different dispersal abilities

Although today the forest cover is continuous in Central Africa, this may have not always been the case, as the scarce fossil record in this region suggests that arid conditions might have significantly reduced tree density during the ice ages. Our aim was to investigate whether the dry ice age periods left a genetic signature on tree species that can be used to infer the date of the past fragmentation of the rainforest. We sequenced reduced representation libraries of 182 samples representing five widespread legume trees and seven outgroups. Phylogenetic analyses identified an early divergent lineage for all species in West Africa (Upper Guinea) and two clades in Central Africa: Lower Guinea-North and Lower Guinea-South. As the structure separating the Northern and Southern clades-congruent across species-cannot be explained by geographic barriers, we tested other hypotheses with demographic model testing using δαδι. The best estimates indicate that the two clades split between the Upper Pliocene and the Pleistocene, a date compatible with forest fragmentation driven by ice age climatic oscillations. Furthermore, we found remarkably older split dates for the shade-tolerant tree species with nonassisted seed dispersal than for light-demanding species with long-distance wind-dispersed seeds. Different recolonization abilities after recurrent cycles of forest fragmentation seem to explain why species with long-distance dispersal show more recent genetic admixture between the two clades than species with limited seed dispersal. Despite their old history, our results depict the African rainforests as a dynamic biome where tree species have expanded relatively recently after the last glaciation.

The biogeography of Dromiciops in southern South America: Middle Miocene transgressions, speciation and associations with Nothofagus

The current distribution of the flora and fauna of southern South America is the result of drastic geological events that occurred during the last 20 million years, including marine transgressions, glaciations and active vulcanism. All these have been associated with fragmentation, isolation and subsequent expansion of the biota, south of 35°S, such as the temperate rainforest. This forest is mostly dominated by Nothofagus trees and is the habitat of the relict marsupial monito del monte, genus Dromiciops, sole survivor of the order Microbiotheria. Preliminary analyses using mtDNA proposed the existence of three main Dromiciops lineages, distributed latitudinally, whose divergence was initially attributed to recent Pleistocene glaciations. Using fossil-calibrated dating on nuclear and mitochondrial genes, here we reevaluate this hypothesis and report an older (Miocene) biogeographic history for the genus. We performed phylogenetic reconstructions using sequences from two mitochondrial DNA and four nuclear DNA genes in 159 specimens from 31 sites across Chile and Argentina. Our phylogenetic analysis resolved three main clades with discrete geographic distributions. The oldest and most differentiated clade corresponds to that of the northern distribution (35.2°S to 39.3°S), which should be considered a distinct species (D. bozinovici, sensu D'Elía et al. 2016). According to our estimations, this species shared a common ancestor with D. gliroides (southern clades) about ~13 million years ago. Divergence time estimates for the southern clades (39.6°S to 42.0°S) ranged from 9.57 to 6.5 Mya. A strong genetic structure was also detected within and between clades. Demographic analyses suggest population size stability for the northern clade (D. bozinovici), and recent demographic expansions for the central and southern clades. All together, our results suggest that the diversification of Dromiciops were initiated by the Middle Miocene transgression (MMT), the massive marine flooding that covered several lowlands of the western face of Los Andes between 37 and 48°S. The MMT resulted from an increase in global sea levels at the Miocene climatic optimum, which shaped the biogeographic origin of several species, including Nothofagus forests, the habitat of Dromiciops.

Three new species of Trichoderma in the Harzianum and Longibrachiatum lineages from Peruvian cacao crop soils based on an integrative approach

The hyperdiverse genus Trichoderma is one of most useful groups of microbes for a number of human activities, and their accurate identification is crucial. The structural simplicity and lack of distinctive phenotypic variation in this group enable the use of DNA-based species delimitation methods in combination with phylogenies (and morphology when feasible) to establish well-supported boundaries among species. Our study employed a multilocus phylogeny and four DNA-based methods (automated barcode gap discovery [ABGD], statistical parsimony [SPN], generalized mixed Yule coalescent [GMYC], and Bayesian phylogenetics and phylogeography [BPP]) for four molecular markers (acl1, act, rpb2, and tef1) to delimit species of two lineages of Trichoderma. Although incongruence among these methods was observed in our analyses, the genetic distance (ABGD) and coalescence (BPP) methods and the multilocus phylogeny strongly supported and confirmed recognition of 108 and 39 different species in the Harzianum and Longibrachiatum lineages, including three new species associated with cacao farms in northern Peru, namely, T.awajun, sp. nov., T. jaklitschii, sp. nov., and T. peruvianum, sp. nov. Morphological distinctions between the new species and their close relatives are primarily related to growth rates, colony appearance, and size of phialides and conidia. This study confirmed that an integrative approach (DNA-based methods, multilocus phylogeny, and phenotype) is more likely to reliably verify supported species boundaries in Trichoderma.

WP, Sterkhova VV, Kusamba C, Rödel M, Penner J, Peterson AT, Brown RM. Rivers, not refugia, drove diversification in arboreal, sub-Saharan African snakes

The relative roles of rivers versus refugia in shaping the high levels of species diversity in tropical rainforests have been widely debated for decades. Only recently has it become possible to take an integrative approach to test predictions derived from these hypotheses using genomic sequencing and paleo-species distribution modeling. Herein, we tested the predictions of the classic river, refuge, and river-refuge hypotheses on diversification in the arboreal sub-Saharan African snake genus Toxicodryas. We used dated phylogeographic inferences, population clustering analyses, demographic model selection, and paleo-distribution modeling to conduct a phylogenomic and historical demographic analysis of this genus. Our results revealed significant population genetic structure within both Toxicodryas species, corresponding geographically to river barriers and divergence times from the mid-Miocene to Pliocene. Our demographic analyses supported the interpretation that rivers are indications of strong barriers to gene flow among populations since their divergence. Additionally, we found no support for a major contraction of suitable habitat during the last glacial maximum, allowing us to reject both the refuge and river-refuge hypotheses in favor of the river-barrier hypothesis. Based on conservative interpretations of our species delimitation analyses with the Sanger and ddRAD data sets, two new cryptic species are identified from east-central Africa. This study highlights the complexity of diversification dynamics in the African tropics and the advantages of integrative approaches to studying speciation in tropical regions.

Molecular Phylogeny of Cryptonanus (Didelphidae: Thylamyini): Evidence for a recent and complex diversification in South American open biomes

Systematic revisions of South American marsupials have contributed to our knowledge about genus and species diversity in the last decades, including studies of the most recently described genus Cryptonanus (Didelphidae), currently comprising four recognized species. Herein we provide the first phylogeny for these mouse opossums based on comprehensive sampling, including representatives from all nominal taxa, encompassing most of the geographic distribution of the genus while also extending its known range. The taxonomic status of Cryptonanus species was explored by analyses of multiple mitochondrial and nuclear DNA markers to assess phylogenetic relationships and to provide divergence time estimates, species delimitations and biogeographical hypotheses. Cryptonanus monophyly remained highly supported despite the inclusion of abundant new data from more than a hundred specimens, comprising 10 independent evolutionary lineages. Species-complexes within valid nominal taxa reveal higher species richness in the genus. Based on divergence estimates from a dated phylogeny, we suggest that Cryptonanus diversified along the Quaternary, with speciation events occurring well into the Pleistocene. The best supported biogeographical hypothesis endorses speciation by vicariance and subset speciation across open formations in shaping the evolutionary history of this didelphid genus, strongly associated with dry tropical landscapes of South America.

Incorporating the speciation process into species delimitation

The "multispecies" coalescent (MSC) model that underlies many genomic species-delimitation approaches is problematic because it does not distinguish between genetic structure associated with species versus that of populations within species. Consequently, as both the genomic and spatial resolution of data increases, a proliferation of artifactual species results as within-species population lineages, detected due to restrictions in gene flow, are identified as distinct species. The toll of this extends beyond systematic studies, getting magnified across the many disciplines that rely upon an accurate framework of identified species. Here we present the first of a new class of approaches that addresses this issue by incorporating an extended speciation process for species delimitation. We model the formation of population lineages and their subsequent development into independent species as separate processes and provide for a way to incorporate current understanding of the species boundaries in the system through specification of species identities of a subset of population lineages. As a result, species boundaries and within-species lineages boundaries can be discriminated across the entire system, and species identities can be assigned to the remaining lineages of unknown affinities with quantified probabilities. In addition to the identification of species units in nature, the primary goal of species delimitation, the incorporation of a speciation model also allows us insights into the links between population and species-level processes. By explicitly accounting for restrictions in gene flow not only between, but also within, species, we also address the limits of genetic data for delimiting species. Specifically, while genetic data alone is not sufficient for accurate delimitation, when considered in conjunction with other information we are able to not only learn about species boundaries, but also about the tempo of the speciation process itself.

Genetic, morphological and acoustic differentiation of African trident bats (Rhinonycteridae: Triaenops)

Rhinonycteridae (trident bats) are a small Palaeotropical family of insectivorous bats allied to Hipposideridae. Their taxonomy has been in a state of flux. Here, we use mitochondrial and nuclear sequences to evaluate species relationships, confirming the monophyly of both Triaenops and Paratriaenops. Although most Triaenops afer specimens are recovered as a group, mitochondrial analyses strongly support some Kenyan individuals as members of Triaenops persicus. Analyses of four nuclear introns (ACOX2, COPS7A, RODGI and STAT5A) strongly support the mitochondrial topology. Morphometric analysis of the skull, external morphology and echolocation calls confirm that the Triaenops from the Rift Valley in Kenya (Nakuru, Baringo and Pokot counties) are distinct from typical T. afer in coastal (Kilifi and Kwale counties) or interior (Laikipia and Makueni counties) colonies. We interpret these analyses to indicate that two species of Triaenops occur in East Africa: T. afer in coastal regions along the Indian Ocean and in the highlands of central Kenya and Ethiopia, and T. persicus in the Rift Valley of Kenya. Although they appear widely disjunct from Middle Eastern populations, Kenyan T. persicus might be more widely distributed in the Rift Valley; they are somewhat differentiated from Middle Eastern populations in terms of both cranial morphology and vocalizations.