Testing the molecular and evolutionary causes of a ‘leapfrog’ pattern of geographical variation in coloration

Divergence With Gene Flow and Contrasting Population Size Blur the Species Boundary in Cycas Sect. Asiorientales, as Inferred From Morphology and RAD-Seq Data

The divergence process of incipient species is fascinating but elusive by incomplete lineage sorting or gene flow. Species delimitation is also challenging among those morphologically similar allopatric species, especially when lacking comprehensive data. Cycas sect. Asiorientales, comprised of C. taitungensis and C. revoluta in the Ryukyu Archipelago and Taiwan, diverged recently with continuous gene flow, resulting in a reciprocal paraphyletic relationship. Their previous evolutionary inferences are questioned from few genetic markers, incomplete sampling, and incomprehensive morphological comparison by a long-term taxonomic misconception. By whole range sampling, this study tests the geographic mode of speciation in the two species of Asiorientales by approximate Bayesian computation (ABC) using genome-wide single nucleotide polymorphisms (SNPs). The individual tree was reconstructed to delimit the species and track the gene-flow trajectory. With the comparison of diagnostic morphological traits and genetic data, the allopatric speciation was rejected. Alternatively, continuous but spatially heterogeneous gene flow driven by transoceanic vegetative dispersal and pollen flow with contrasting population sizes blurred their species boundary. On the basis of morphological, genetic, and evolutionary evidence, we synonymized these two Cycas species. This study highlights not only the importance of the Kuroshio Current to species evolution but also the disadvantage of using species with geographically structured genealogies as conservation units.

Divergence, gene flow, and the origin of leapfrog geographic distributions: The history of colour pattern variation in Phyllobates poison-dart frogs

The geographic distribution of phenotypic variation among closely related populations is a valuable source of information about the evolutionary processes that generate and maintain biodiversity. Leapfrog distributions, in which phenotypically similar populations are disjunctly distributed and separated by one or more phenotypically distinct populations, represent geographic replicates for the existence of a phenotype, and are therefore especially informative. These geographic patterns have mostly been studied from phylogenetic perspectives to understand how common ancestry and divergent evolution drive their formation. Other processes, such as gene flow between populations, have not received as much attention. Here, we investigate the roles of divergence and gene flow between populations in the origin and maintenance of a leapfrog distribution in Phyllobates poison frogs. We found evidence for high levels of gene flow between neighbouring populations but not over long distances, indicating that gene flow between populations exhibiting the central phenotype may have a homogenizing effect that maintains their similarity, and that introgression between 'leapfroging' taxa has not played a prominent role as a driver of phenotypic diversity in Phyllobates. Although phylogenetic analyses suggest that the leapfrog distribution was formed through independent evolution of the peripheral (i.e. leapfrogging) populations, the elevated levels of gene flow between geographically close populations poise alternative scenarios, such as the history of phenotypic change becoming decoupled from genome-averaged patterns of divergence, which we cannot rule out. These results highlight the importance of incorporating gene flow between populations into the study of geographic variation in phenotypes, both as a driver of phenotypic diversity and as a confounding factor of phylogeographic inferences.

Spatial predictors of genomic and phenotypic variation differ in a lowland Middle American bird (Icterus gularis)

Spatial patterns of intraspecific variation are shaped by geographical distance among populations, historical changes in gene flow and interactions with local environments. Although these factors are not mutually exclusive and operate on both genomic and phenotypic variation, it is unclear how they affect these two axes of variation. We address this question by exploring the predictors of genomic and phenotypic divergence in Icterus gularis, a broadly distributed Middle American bird that exhibits marked geographical variation in body size across its range. We combined a comprehensive single nucleotide polymorphism and phenotypic data set to test whether genome-wide genetic and phenotypic differentiation are best explained by (i) isolation by distance, (ii) isolation by history or (iii) isolation by environment. We find that the pronounced genetic and phenotypic variation in I. gularis are only partially correlated and differ regarding spatial predictors. Whereas genomic variation is largely explained by historical barriers to gene flow, phenotypic diversity can be best predicted by contemporary environmental heterogeneity. Our genomic analyses reveal strong phylogeographical structure coinciding with the Chivela Pass at the Isthmus of Tehuantepec that was formed during the Pleistocene, when populations were isolated in north-south refugia. In contrast, we found a strong association between body size and environmental variables, such as temperature and precipitation. The relationship between body size and local climate is consistent with a pattern produced by either natural selection or environmental plasticity. Overall, these results provide empirical evidence for why phenotypic and genomic data are often in conflict in taxonomic and phylogeographical studies.

Genetic drift does not sufficiently explain patterns of electric signal variation among populations of the mormyrid electric fish Paramormyrops kingsleyae

Communication signals serve crucial survival and reproductive functions. In Gabon, the widely distributed mormyrid fish Paramormyrops kingsleyae emits an electric organ discharge (EOD) signal with a dual role in communication and electrolocation that exhibits remarkable variation: populations of P. kingsleyae have either biphasic or triphasic EODs, a feature that characterizes interspecific signal diversity among the Paramormyrops genus. We quantified variation in EODs of 327 P. kingsleyae from nine populations and compared it to genetic variation estimated from microsatellite loci. We found no correlation between electric signal and genetic distances, suggesting that EOD divergence cannot be explained by drift alone. An alternative hypothesis is that EOD differences are used for mate discrimination, which would require P. kingsleyae be capable of differentiating between divergent EOD waveforms. Using a habituation-dishabituation assay, we found that P. kingsleyae can discriminate between biphasic and triphasic EOD types. Nonetheless, patterns of genetic and electric organ morphology divergence provide evidence for hybridization between these signal types. Although reproductive isolation with respect to signal type is incomplete, our results suggest that EOD variation in P. kingsleyae could be a cue for assortative mating.

Phylogeography indicates incomplete genetic divergence among phenotypically differentiated montane forest populations of Atlapetesalbinucha (Aves, Passerellidae)

The White-naped Brushfinch (Atlapetesalbinucha) comprises up to eight allopatric subspecies mainly identified by the color of the underparts (gray vs. yellow belly). Yellow and gray bellied forms were long considered two different species (A.albinucha and A.gutturalis), but they are presently considered as one polytypic species. Previous studies in the genus Atlapetes have shown that the phylogeny, based on molecular data, is not congruent with characters such as coloration, ecology, or distributional patterns. The phylogeography of A.albinucha was analyzed using two mitochondrial DNA regions from samples including 24 different localities throughout montane areas from eastern Mexico to Colombia. Phylogeographic analyses using Bayesian inference, maximum likelihood and haplotype network revealed incomplete geographic structure. The genetic diversity pattern is congruent with a recent process of expansion, which is also supported by Ecological Niche Models (ENM) constructed for the species and projected into three past scenarios. Overall, the results revealed an incomplete genetic divergence among populations of A.albinucha in spite of the species’ ample range, which contrasts with previous results of phylogeographic patterns in other Neotropical montane forest bird species, suggesting idiosyncratic evolutionary histories for different taxa throughout the region.

Comprehensive molecular phylogeny of barn owls and relatives (Family: Tytonidae), and their six major Pleistocene radiations

The owl family Tytonidae comprises two genera: Phodilus, limited to the forests of central Africa and South-East Asia, and the ubiquitous Tyto. The genus Tyto is majorly represented by the cosmopolitan Common Barn Owl group, with more than 30 subspecies worldwide. Discrete differences in body size and plumage colouration have led to the classification of this family into many species and subspecies, but the taxonomic status and phylogenetic relationships between taxa remain unclear, and in some groups controversial. Although several previous studies attempted to resolve this problem, they have been limited in their taxonomic and geographical coverage, or have relied on restricted molecular evidence and low sample sizes. Based on the most comprehensive sampling to date (16 out of 17 Tyto species, and one out of three Phodilus species), a multi-locus approach using seven mitochondrial and two nuclear markers, and taking advantage of field data and museum collections available worldwide, our main questions in this study were: (1) what are the phylogenetic relationships and classification status of the whole family; (2) when and where did the most important speciation events occur? We confirm that the Common Barn Owl, Tyto alba is divided into three main evolutionary units: the American Barn Owl, T. furcata; the Western Barn Owl, T. alba; and the Eastern Barn Owl, T. javanica, and suggest a Late Miocene (ca. 6 mya) Australasian and African origin of the group. Our results are supported by fossil age information, given that the most recent common ancestor between the Tytonidae genera Phodilus and Tyto was probably from the Oligocene (ca. 28 mya) of Australasia. We finally reveal six major Pleistocene radiations of Tyto, all resulting in wide-range distributions.

On geographic barriers and Pleistocene glaciations: Tracing the diversification of the Russet-crowned Warbler (Myiothlypis coronata) along the Andes

We studied the phylogeography and plumage variation of the Russet-crowned Warbler (Myiothlypis coronata), from Venezuela to Bolivia, with focus on populations from Ecuador and northern Peru. We analyzed sequences of mitochondrial and nuclear genes, geographic distributions, as well as photographs of specimens deposited at museum collections. Phylogenetic analyses identified three major lineages formed by populations from: Venezuela and Colombia (M. c. regulus), Ecuador and northern Peru (M. elata, M. castaneiceps, M. orientalis, M. c. chapmani), and central Peru and Bolivia (M. c. coronata). We found further population structure within M. c. regulus and M. c. coronata, and population structure and complexity of plumage variation within the Ecuador-northern Peru lineage. Time-calibrated trees estimated that most intraspecific variation originated during the Pleistocene; however, this pattern may not be attributed to an increase in diversification rate during that period. We discuss these results in the context of the importance of geographic-ecological barriers in promoting lineage diversification along the Andes and put forward a preliminary taxonomic proposal for major lineages identified in this study.

Drift-driven evolution of electric signals in a Neotropical knifefish

Communication signals are highly diverse traits. This diversity is usually assumed to be shaped by selective forces, whereas the null hypothesis of divergence through drift is often not considered. In Panama, the weakly electric fish Brachyhypopomus occidentalis is widely distributed in multiple independent drainage systems, which provide a natural evolutionary laboratory for the study of genetic and signal divergence in separate populations. We quantified geographic variation in the electric signals of 109 fish from five populations, and compared it to the neutral genetic variation estimated from cytochrome oxidase I (COI) sequences of the same individuals, to test whether drift may be driving divergence of their signals. Signal distances were highly correlated with genetic distances, even after controlling for geographic distances, suggesting that drift alone is sufficient to explain geographic variation in electric signals. Significant differences at smaller geographic scales (within drainages) showed, however, that electric signals may evolve at a faster rate than expected under drift, raising the possibility that additional adaptive forces may be contributing to their evolution. Overall, our data point to stochastic forces as main drivers of signal evolution in this species and extend the role of drift in the evolution of communication systems to fish and electrocommunication.

The dual role of Andean topography in primary divergence: functional and neutral variation among populations of the hummingbird, Metallura tyrianthina

Background

The ridges and valleys of the Andes create physical barriers that limit animal dispersal and cause deterministic local variation in rainfall. This has resulted in physical isolation of animal populations and variation in habitats, each of which has likely contributed to the evolution of high species diversity in the region. However, the relative influences of geographic isolation, ecoclimatic conditions, and their potential interactions remain poorly understood. To address this, we compared patterns of genetic and morphological diversity in Peruvian populations of the hummingbird Metallura tyrianthina.

Results

Phylogenetic and variation partitioning analyses showed that geographic isolation rather than climatic dissimilarity explained the greatest proportion of genetic variance. In contrast, bill length variation was explained by climatic seasonality, but not by genetic divergence. We found that mutation-scaled migration rate (m) between persistently humid and semi-humid environments was nearly 20 times higher when the habitats were contiguous (m = 39.9) than when separated by a barrier, the Cordillera de Vilcanota (m = 2.1). Moreover, the population experiencing more gene flow exhibited a lesser degree of bill length divergence despite similar differences in climate.

Conclusions

Geographic isolation is necessary for genetic divergence. Ecological differences, represented here by climate characteristics, are necessary for functional divergence. Gene flow appears to hinder the evolution of functional traits toward local adaptive optima. This suggests that functional diversification requires geographic isolation followed or accompanied by a shift in ecological conditions. Andean topography causes both isolation and climatic variation, underscoring its dual role in biotic diversification.

Electronic supplementary material

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

On the Origin of Pantepui montane biotas: A Perspective Based on the Phylogeny of Aulacorhynchus toucanets

To understand the origin of Pantepui montane biotas, we studied the biogeography of toucanets in the genus Aulacorhynchus. These birds are ideal for analyzing historical relationships among Neotropical montane regions, given their geographic distribution from Mexico south to Bolivia, including northern Venezuela (Cordillera de la Costa), and the Pantepui. Analyses were based on molecular phylogenies using mitochondrial and nuclear DNA sequences. Topology tests were applied to compare alternative hypotheses that may explain the current distribution of Aulacorhynchus toucanets, in the context of previous hypotheses of the origin of Pantepui montane biotas. Biogeographic reconstructions in RASP and Lagrange were used to estimate the ancestral area of the genus, and an analysis in BEAST was used to estimate a time framework for its diversification. A sister relationship between the Pantepui and Andes+Cordillera de la Costa was significantly more likely than topologies indicating other hypothesis for the origin of Pantepui populations. The Andes was inferred as the ancestral area for Aulacorhynchus, and the group has diversified since the late Miocene. The biogeographic patterns found herein, in which the Andes are the source for biotas of other regions, are consistent with those found for flowerpiercers and tanagers, and do not support the hypothesis of the geologically old Pantepui as a source of Neotropical montain diversity. Based on the high potential for cryptic speciation and isolation of Pantepui populations, we consider that phylogenetic studies of additional taxa are important from a conservation perspective.

Delimiting shades of gray: phylogeography of the Northern Fulmar, Fulmarus glacialis

The Northern Fulmar (Fulmarus glacialis) is a common tube-nosed seabird with a disjunct Holarctic range. Taxonomic divisions within the Northern Fulmar have historically been muddled by geographical variation notably including highly polymorphic plumage. Recent molecular analyses (i.e., DNA barcoding) have suggested that genetic divergence between Atlantic and Pacific populations could be on par with those typically observed between species. We employ a multigene phylogenetic analysis to better explore the level of genetic divergence between these populations and to test an old hypothesis on the origin of the modern distribution of color morphs across their range. Additionally, we test whether mutations in the melanocortin-1 receptor gene (MC1R) are associated with dark plumage in the Northern Fulmar. We confirmed that mitochondrial lineages in the Atlantic and Pacific populations are highly divergent, but nuclear markers revealed incomplete lineage sorting. Genetic divergence between these populations is consistent with that observed between many species of Procellariiformes and we recommend elevating these two forms to separate species. We also find that MC1R variation is not associated with color morph but rather is better explained by geographical divergence.

Speciation in an avian complex endemic to the mountains of Middle America (Ergaticus, Aves: Parulidae)

The implementation of the phylogeographic approach for the study of biodiversity is critical in poorly sampled regions like the montane systems of Middle America, as complex evolutionary histories often result in the presence of independent lineages not properly considered by traditional taxonomy. Herein we sequenced 2370 bp of mtDNA (ND2, cyt b and ATPase) from 81 individuals of Ergaticus, a complex of birds endemic to the montane forests of Middle America. Although current taxonomy recognizes two species, the results reveal considerable genetic structure with the presence of four mtDNA lineages. Two of these lineages within Ergaticus ruber evidence the need of a revaluation of the species limits for this taxon. The general phylogeographic pattern can be explained as a consequence of relative isolation of the populations in different mountain ranges separated by low elevation barriers. Most population groups did not show signals of demographic expansion with the exception of the one corresponding to clade 1. The divergence time estimates point to the Pleistocene as an important time period for the diversification of this complex.

Evolution of arid zone birds in Australia: leapfrog distribution patterns and mesic-arid connections in quail-thrush (Cinclosoma, Cinclosomatidae)

The quail-thrush, Cinclosoma, include between five and seven species distributed broadly across arid and semi-arid inland Australia, mesic forests of south-eastern Australia and New Guinea. It has been suggested that the arid zone species of quail-thrush arose from forest ancestors as Australia changed from a warm wet climate to a cooler drier climate since the late-Miocene. We generated multilocus (mitochondrial ND2 and eight nuclear loci) gene and species trees with complete taxon sampling of Cinclosoma to investigate evolutionary relationships and species status of some taxa. Topologies reconstructed in congruent, highly-resolved gene trees and species trees that supported the recognition of seven species. Ancestral state reconstruction and divergence time estimates suggest that arid-adapted taxa radiated in parallel with a drying climate and changing habitat. A 'leapfrog' distribution in phenotypes of arid zone taxa was likely a result of ancestral retention of inconspicuous (or camouflaged) plumage patterns. A specimen-based report from 1968 of hybridization between non-sister taxa Cinclosoma castanotum and Cinclosoma marginatum was verified using molecular analysis on specimens collected at the same locality 40 years later. We discuss the implications of hybridization to the evolution of this species group.

Selection and geographic isolation influence hummingbird speciation: genetic, acoustic and morphological divergence in the wedge-tailed sabrewing (Campylopterus curvipennis)

BACKGROUND

Mesoamerica is one of the most threatened biodiversity hotspots in the world, yet we are far from understanding the geologic history and the processes driving population divergence and speciation for most endemic taxa. In species with highly differentiated populations selective and/or neutral factors can induce rapid changes to traits involved in mate choice, promoting reproductive isolation between allopatric populations that can eventually lead to speciation. We present the results of genetic differentiation, and explore drift and selection effects in promoting acoustic and morphological divergence among populations of Campylopterus curvipennis, a lekking hummingbird with an extraordinary vocal variability across Mesoamerica.

RESULTS

Analyses of two mitochondrial genes and ten microsatellite loci genotyped for 160 individuals revealed the presence of three lineages with no contemporary gene flow: C. c. curvipennis, C. c. excellens, and C. c. pampa disjunctly distributed in the Sierra Madre Oriental, the Tuxtlas region and the Yucatan Peninsula, respectively. Sequence mtDNA and microsatellite data were congruent with two diversification events: an old vicariance event at the Isthmus of Tehuantepec (c. 1.4 Ma), and a more recent Pleistocene split, isolating populations in the Tuxtlas region. Hummingbirds of the excellens group were larger, and those of the pampa group had shorter bills, and lineages that have been isolated the longest shared fewer syllables and differed in spectral and temporal traits of a shared syllable. Coalescent simulations showed that fixation of song types has occurred faster than expected under neutrality but the null hypothesis that morphological divergence resulted from drift was not rejected.

CONCLUSIONS

Our phylogeographic analyses uncovered the presence of three Mesoamerican wedge-tailed sabrewing lineages, which diverged at different time scales. These results highlight the importance of the Isthmus of Tehuantepec and more recent Pleistocene climatic events in driving isolation and population divergence. Coalescent analyses of the evolution of phenotypic traits suggest that selection is driving song evolution in wedge-tailed sabrewings but drift could not be rejected as a possibility for morphological divergence.