Speciation and secondary contact in a fossorial island endemic, the São Tomé caecilian

Reproductive isolation in a three-way contact zone

Contact zones between divergent forms within a species provide insight into the role of gene flow in adaptation and speciation. Previous work has focused on contact zones involving only two divergent forms, but in nature, many more than two populations may overlap simultaneously and experience gene flow. Patterns of introgression in wild populations are, therefore, likely much more complicated than is often assumed. We begin to address this gap in current knowledge by investigating patterns of divergence and introgression across a complex natural contact zone. We use phenotypic and genomic data to confirm the existence of a three-way contact zone among divergent freshwater resident, saltwater resident and saltwater migratory three-spined stickleback (Gasterosteus aculeatus) on the island of North Uist, Scottish Western Isles. We find evidence for hybridization, mostly between saltwater resident and saltwater migratory forms. Despite hybridization, genomic analyses reveal pairwise islands of divergence between all forms that are maintained across the contact zone. Genomic cline analyses also provide evidence for selection and/or hybrid incompatibilities in divergent regions. Divergent genomic regions occur across multiple chromosomes and involve many known adaptive loci and several chromosomal inversions. We also identify distinct immune gene expression profiles between forms, but no evidence for transgressive expression in hybrids. Our results suggest that reproductive isolation is maintained in this three-way contact zone, despite some hybridization, and that reduced recombination in chromosomal inversions may play an important role in maintaining this isolation.

Discerning structure versus speciation in phylogeographic analysis of Seepage Salamanders (Desmognathus aeneus) using demography, environment, geography, and phenotype

Numerous mechanisms can drive speciation, including isolation by adaptation, distance, and environment. These forces can promote genetic and phenotypic differentiation of local populations, the formation of phylogeographic lineages, and ultimately, completed speciation. However, conceptually similar mechanisms may also result in stabilizing rather than diversifying selection, leading to lineage integration and the long-term persistence of population structure within genetically cohesive species. Processes that drive the formation and maintenance of geographic genetic diversity while facilitating high rates of migration and limiting phenotypic differentiation may thereby result in population genetic structure that is not accompanied by reproductive isolation. We suggest that this framework can be applied more broadly to address the classic dilemma of "structure" versus "species" when evaluating phylogeographic diversity, unifying population genetics, species delimitation, and the underlying study of speciation. We demonstrate one such instance in the Seepage Salamander (Desmognathus aeneus) from the southeastern United States. Recent studies estimated up to 6.3% mitochondrial divergence and four phylogenomic lineages with broad admixture across geographic hybrid zones, which could potentially represent distinct species supported by our species-delimitation analyses. However, while limited dispersal promotes substantial isolation by distance, microhabitat specificity appears to yield stabilizing selection on a single, uniform, ecologically mediated phenotype. As a result, climatic cycles promote recurrent contact between lineages and repeated instances of high migration through time. Subsequent hybridization is apparently not counteracted by adaptive differentiation limiting introgression, leaving a single unified species with deeply divergent phylogeographic lineages that nonetheless do not appear to represent incipient species.

Habitat-associated Genomic Variation in a Wall Lizard from an Oceanic Island

The lizard Teira dugesii exhibits morphological divergence between beach and inland habitats in the face of gene flow, within the volcanic island of Madeira, Portugal. Here, we analyzed genomic data obtained by genotyping-by-sequencing, which provided 16,378 single nucleotide polymorphisms (SNPs) from 94 individuals sampled from 15 sites across Madeira. Ancient within-island divergence in allopatry appears to have mediated divergence in similar species within other Atlantic islands, but this hypothesis was not supported for T. dugesii. Across all samples, a total of 168 SNPs were classified as statistical outliers using pcadapt and OutFLANK. Redundancy analysis (RDA) revealed that 17 of these outliers were associated with beach/inland habitats. The SNPs were located within 16 sequence tags and 15 of these were homologous with sequences in a 31 Mb region on chromosome 3 of a reference wall lizard genome (the remaining tag could not be associated with any chromosome). We further investigated outliers through contingency analyses of allele frequencies at each of four pairs of adjacent beach-inland sites. The majority of the outliers detected by the RDA were confirmed at two pairs of these matched sites. These analyses also suggested some parallel divergence at different localities. Six other outliers were associated with site elevation, four of which were located on chromosome 5 of the reference genome. Our study lends support to a previous hypothesis that divergent selection between gray shingle beaches and inland regions overcomes gene flow and leads to the observed morphological divergence between populations in these adjacent habitats.

Illustrated keys and a DNA barcode reference library of the amphibians and terrestrial reptiles (Amphibia, Reptilia) of São Tomé and Príncipe (Gulf of Guinea, West Africa)

The herpetofauna of São Tomé and Príncipe consists of nine species of amphibians, all endemic, and 21 species of terrestrial reptiles, of which 17 are endemic. Our current knowledge regarding its natural history, ecology, and distribution is limited. Here two important tools are provided to support researchers, conservationists, and local authorities in the identification of the country's herpetofauna: an illustrated key to the herpetofauna of the two islands and surroundings islets and a DNA barcode reference library. The keys allow a rapid and unambiguous morphological identification of all occurring species. The DNA barcodes for the entire herpetofauna of the country were produced from 79 specimens, all of which are deposited in museum collections. The barcodes generated are available in online repositories and can be used to provide unambiguous molecular identification of most of the species. Future applications and use of these tools are briefly discussed.

Life on a beach leads to phenotypic divergence despite gene flow for an island lizard

Limited spatial separation within small islands suggests that observed population divergence may occur due to habitat differences without interruption to gene flow but strong evidence of this is scarce. The wall lizard Teira dugesii lives in starkly contrasting shingle beach and inland habitats on the island of Madeira. We used a matched pairs sampling design to examine morphological and genomic divergence between four beach and adjacent (<1 km) inland areas. Beach populations are significantly darker than corresponding inland populations. Geometric morphometric analyses reveal divergence in head morphology: beach lizards have generally wider snouts. Genotyping-by-sequencing allows the rejection of the hypothesis that beach populations form a distinct lineage. Bayesian analyses provide strong support for models that incorporate gene flow, relative to those that do not, replicated at all pairs of matched sites. Madeiran lizards show morphological divergence between habitats in the face of gene flow, revealing how divergence may originate within small islands.

Revealing the cryptic diversity of the widespread and poorly known South American blind snake genus Amerotyphlops (Typhlopidae: Scolecophidia) through integrative taxonomy

Morphological stasis is generally associated with relative constancy in ecological pressures throughout time, producing strong stabilizing selection that retains similar shared morphology. Although climate and vegetation are commonly the main key factors driving diversity and phenotypic diversification in terrestrial vertebrates, fossorial organisms have their morphology mostly defined by their fossorial lifestyle. Among these secretive fossorial organisms, blind snakes of the South American genus Amerotyphlops are considered poorly studied when compared to other taxa. Here, we evaluate the cryptic diversity of Amerotyphlops using phylogenetic and multivariate approaches. We based our phylogenetic analysis on a molecular dataset composed of 12 gene fragments (eight nuclear and four mitochondrial) for 109 species of Typhlopidae. The multivariate analysis was implemented using 36 morphological variables for 377 specimens of Amerotyphlops. Additionally, we contrast our phylogenetic result with the morphological variation found in cranial, external and hemipenial traits. Our phylogenetic results recovered with strong support the following monophyletic groups within Amerotyphlops: (1) a clade formed by A. tasymicris and A. minuisquamus; (2) a clade composed of A. reticulatus; (3) a north-eastern Brazilian clade including A. yonenagae, A. arenensis, A. paucisquamus and A. amoipira; and (4) a clade composed of A. brongersmianus and a complex of cryptic species. Based on these results we describe four new species of Amerotyphlops from north-eastern and south-eastern Brazil, which can be distinguished from the morphologically similar species, A. brongersmianus and A. arenensis.

Phylogeography and evolutionary lineage diversity in the small-eared greater galago, Otolemur garnettii (Primates: Galagidae)

Assessing the true lineage diversity in elusive nocturnal organisms is particularly challenging due to their subtle phenotypic variation in diagnostic traits. The cryptic small-eared greater galago (Otolemur garnettii) offers a great opportunity to test if currently recognized subspecies, suggested by discontinuities in coat colour pattern and geographic barriers, represent distinct evolutionary lineages. To answer this question, we conducted the first population-level phylogeographic study of the species, sampling wild specimens from across almost its entire latitudinal range, including the Zanzibar Archipelago. We applied five species-delimitation algorithms to investigate the genetic diversity and distribution pattern of mitochondrial DNA across the geographic range of three out of four subspecies. Our results suggest that far-northern populations of O. g. lasiotis potentially represent an independently evolving lineage, but populations assigned to O. g. garnettii from Zanzibar Island and of O. g panganiensis from mainland Tanzania do not constitute two independent lineages. A dated phylogeny suggests that this northern clade diverged from all remaining samples approximately 4 Mya. Such old divergence age is in line with the split between many galagid species. This northern lineage could potentially represent an incipient species; however, there is not yet enough evidence to support a new taxonomic status for this unique mitochondrial group.

Linked selection, differential introgression and recombination rate variation promote heterogeneous divergence in a pair of yellow croakers

Understanding the mechanisms underlying heterogeneous genomic divergence is of particular interest in evolutionary biology. Highly differentiated genomic regions, known as genomic islands, often evolve between diverging lineages. These genomic islands may be related to selection promoting adaptation or reproductive isolation. Based on whole genome assembly and genome-wide RAD sequencing in a pair of yellow croakers (genus: Larimichthys), we investigated the evolutionary processes shaping genomic landscapes of divergence. Demographic modelling indicated that the two species diverged following a secondary contact scenario, where differential introgression and linked selection were suggested to be involved in heterogeneous genomic divergence. We identified reduced recombination rate in genomic islands and a relatively good conservation of both genetic diversity and recombination landscapes between species, which highlight the roles of linked selection and recombination rate variation in promoting heterogeneous divergence in the common ancestral lineage of the two species. In addition, we found a positive correlation between differentiation (F_ST ) and absolute sequence divergence (D_xy ), and elevated D_xy in genomic islands, indicating that the genomic landscape of divergence was not shaped by linked selection alone. Restricted gene flow in highly differentiated regions has probably remodelled the landscape of heterogeneous genomic divergence. This study highlights that highly differentiated genomic regions can also arise from a combination of linked selection and differential gene flow in interaction with varying recombination rates.

Phylogenomics, introgression, and demographic history of South American true toads (Rhinella)

The effects of genetic introgression on species boundaries and how they affect species' integrity and persistence over evolutionary time have received increased attention. The increasing availability of genomic data has revealed contrasting patterns of gene flow across genomic regions, which impose challenges to inferences of evolutionary relationships and of patterns of genetic admixture across lineages. By characterizing patterns of variation across thousands of genomic loci in a widespread complex of true toads (Rhinella), we assess the true extent of genetic introgression across species thought to hybridize to extreme degrees based on natural history observations and multi-locus analyses. Comprehensive geographic sampling of five large-ranged Neotropical taxa revealed multiple distinct evolutionary lineages that span large geographic areas and, at times, distinct biomes. The inferred major clades and genetic clusters largely correspond to currently recognized taxa; however, we also found evidence of cryptic diversity within taxa. While previous phylogenetic studies revealed extensive mito-nuclear discordance, our genetic clustering analyses uncovered several admixed individuals within major genetic groups. Accordingly, historical demographic analyses supported that the evolutionary history of these toads involved cross-taxon gene flow both at ancient and recent times. Lastly, ABBA-BABA tests revealed widespread allele sharing across species boundaries, a pattern that can be confidently attributed to genetic introgression as opposed to incomplete lineage sorting. These results confirm previous assertions that the evolutionary history of Rhinella was characterized by various levels of hybridization even across environmentally heterogeneous regions, posing exciting questions about what factors prevent complete fusion of diverging yet highly interdependent evolutionary lineages.

Saving collections: taxonomic revision of the herpetological collection of the Instituto de Investigação Científica Tropical, Lisbon (Portugal) with a protocol to rescue abandoned collections

The herpetological collections of the Instituto de Investigação Científica Tropical, Lisbon, are amongst the most important collections from the former Portuguese territories in Africa and Asia. The collection comprises more than 5000 preserved specimens, including type specimens of nine taxa, Trachylepis adamastor, Trachypelis thomensis, Panaspis thomensis, Naja peroescobari, Dalophia angolensis, Hemidactylus nzingae, Boaedon fradei, Platysaurus maculatus maculatus, and Platysaurus maculatus lineicauda. The collection was abandoned in the early years of 2000s and was at risk of being lost. In this paper the entire collection is reviewed, a catalogue provided of the extant specimens, and a brief account of the history of herpetological research at IICT given. Details are also provided on the recovery of the collection and a protocol to rescue abandoned collections.