Divergence with gene flow as facilitated by ecological differences: within-island variation in Darwin's finches

Comparative phylogeography reveals widespread cryptic diversity driven by ecology in Panamanian birds

Widespread species often harbor unrecognized genetic diversity, and investigating the factors associated with such cryptic variation can help us better understand the forces driving diversification. Here, we identify potential cryptic species based on a comprehensive dataset of COI mitochondrial DNA barcodes from 2,333 individual Panamanian birds across 429 species, representing 391 (59%) of the 659 resident landbird species of the country, as well as opportunistically sampled waterbirds. We complement this dataset with additional publicly available mitochondrial loci, such as ND2 and cytochrome b, obtained from whole mitochondrial genomes from 20 taxa. Using barcode identification numbers (BINs), we find putative cryptic species in 19% of landbird species, highlighting hidden diversity in the relatively well-described avifauna of Panama. Whereas some of these mitochondrial divergence events corresponded with recognized geographic features that likely isolated populations, such as the Cordillera Central highlands, the majority (74%) of lowland splits were between eastern and western populations. The timing of these splits are not temporally coincident across taxa, suggesting that historical events, such as the formation of the Isthmus of Panama and Pleistocene climatic cycles, were not the primary drivers of cryptic diversification. Rather, we observed that forest species, understory species, insectivores, and strongly territorial species-all traits associated with lower dispersal ability-were all more likely to have multiple BINs in Panama, suggesting strong ecological associations with cryptic divergence. Additionally, hand-wing index, a proxy for dispersal capability, was significantly lower in species with multiple BINs, indicating that dispersal ability plays an important role in generating diversity in Neotropical birds. Together, these results underscore the need for evolutionary studies of tropical bird communities to consider ecological factors along with geographic explanations, and that even in areas with well-known avifauna, avian diversity may be substantially underestimated.

LAY SUMMARY

- What factors are common among bird species with cryptic diversity in Panama? What role do geography, ecology, phylogeographic history, and other factors play in generating bird diversity?- 19% of widely-sampled bird species form two or more distinct DNA barcode clades, suggesting widespread unrecognized diversity.- Traits associated with reduced dispersal ability, such as use of forest understory, high territoriality, low hand-wing index, and insectivory, were more common in taxa with cryptic diversity. Filogeografía comparada revela amplia diversidad críptica causada por la ecología en las aves de Panamá.

RESUMEN

Especies extendidas frecuentemente tiene diversidad genética no reconocida, y investigando los factores asociados con esta variación críptica puede ayudarnos a entender las fuerzas que impulsan la diversificación. Aquí, identificamos especies crípticas potenciales basadas en un conjunto de datos de códigos de barras de ADN mitocondrial de 2,333 individuos de aves de Panama en 429 especies, representando 391 (59%) de las 659 especies de aves terrestres residentes del país, además de algunas aves acuáticas muestreada de manera oportunista. Adicionalmente, complementamos estos datos con secuencias mitocondriales disponibles públicamente de otros loci, tal como ND2 o citocroma b, obtenidos de los genomas mitocondriales completos de 20 taxones. Utilizando los números de identificación de código de barras (en ingles: BINs), un sistema taxonómico numérico que proporcina una estimación imparcial de la diversidad potencial a nivel de especie, encontramos especies crípticas putativas en 19% de las especies de aves terrestres, lo que destaca la diversidad oculta en la avifauna bien descrita de Panamá. Aunque algunos de estos eventos de divergencia conciden con características geográficas que probablemente aislaron las poblaciones, la mayoría (74%) de la divergencia en las tierras bajas se encuentra entre las poblaciones orientales y occidentales. El tiempo de esta divergencia no coincidió entre los taxones, sugiriendo que eventos históricos tales como la formación del Istmo de Panamá y los ciclos climáticos del pleistoceno, no fueron los principales impulsores de la especiación. En cambio, observamos asociaciones fuertes entre las características ecológicas y la divergencia mitocondriale: las especies del bosque, sotobosque, con una dieta insectívora, y con territorialidad fuerte mostraton múltiple BINs probables. Adicionalmente, el índice mano-ala, que está asociado a la capacidad de dispersión, fue significativamente menor en las especies con BINs multiples, sugiriendo que la capacidad de dispersión tiene un rol importamente en la generación de la diversidad de las aves neotropicales. Estos resultos demonstran la necesidad de que estudios evolutivos de las comunidades de aves tropicales consideren los factores ecológicos en conjunto con las explicaciones geográficos. Palabras clave: biodiversidad tropical, biogeografía, códigos de barras, dispersión, especies crípticas.

Comparative phylogeography reveals widespread cryptic diversity driven by ecology in Panamanian birds

Widespread species often harbor unrecognized genetic diversity, and investigating the factors associated with such cryptic variation can help us better understand the forces driving diversification. Here, we identify potential cryptic species based on a comprehensive dataset of COI mitochondrial DNA barcodes from 2,333 individual Panamanian birds across 429 species, representing 391 (59%) of the 659 resident landbird species of the country, as well as opportunistically sampled waterbirds. We complement this dataset with additional publicly available mitochondrial loci, such as ND2 and cytochrome b, obtained from whole mitochondrial genomes from 20 taxa. Using barcode identification numbers (BINs), we find putative cryptic species in 19% of landbird species, highlighting hidden diversity in the relatively well-described avifauna of Panama. Whereas some of these mitochondrial divergence events corresponded with recognized geographic features that likely isolated populations, such as the Cordillera Central highlands, the majority (74%) of lowland splits were between eastern and western populations. The timing of these splits are not temporally coincident across taxa, suggesting that historical events, such as the formation of the Isthmus of Panama and Pleistocene climatic cycles, were not the primary drivers of cryptic diversification. Rather, we observed that forest species, understory species, insectivores, and strongly territorial species-all traits associated with lower dispersal ability-were all more likely to have multiple BINs in Panama, suggesting strong ecological associations with cryptic divergence. Additionally, hand-wing index, a proxy for dispersal capability, was significantly lower in species with multiple BINs, indicating that dispersal ability plays an important role in generating diversity in Neotropical birds. Together, these results underscore the need for evolutionary studies of tropical bird communities to consider ecological factors along with geographic explanations, and that even in areas with well-known avifauna, avian diversity may be substantially underestimated.

Repeated divergence of amphibians and reptiles across an elevational gradient in northern Madagascar

How environmental factors shape patterns of biotic diversity in tropical ecosystems is an active field of research, but studies examining the possibility of ecological speciation in terrestrial tropical ecosystems are scarce. We use the isolated rainforest herpetofauna on the Montagne d'Ambre (Amber Mountain) massif in northern Madagascar as a model to explore elevational divergence at the level of populations and communities. Based on intensive sampling and DNA barcoding of amphibians and reptiles along a transect ranging from ca. 470-1470 m above sea level (a.s.l.), we assessed a main peak in species richness at an elevation of ca. 1000 m a.s.l. with 41 species. The proportion of local endemics was highest (about 1/3) at elevations >1100 m a.s.l. Two species of chameleons (Brookesia tuberculata, Calumma linotum) and two species of frogs (Mantidactylus bellyi, M. ambony) studied in depth by newly developed microsatellite markers showed genetic divergence up the slope of the mountain, some quite strong, others very weak, but in each case with genetic breaks between 1100 and 1270 m a.s.l. Genetic clusters were found in transect sections significantly differing in bioclimate and herpetological community composition. A decrease in body size was detected in several species with increasing elevation. The studied rainforest amphibians and reptiles show concordant population genetic differentiation across elevation along with morphological and niche differentiation. Whether this parapatric or microallopatric differentiation will suffice for the completion of speciation is, however, unclear, and available phylogeographic evidence rather suggests that a complex interplay between ecological and allopatric divergence processes is involved in generating the extraordinary species diversity of Madagascar's biota. Our study reveals concordant patterns of diversification among main elevational bands, but suggests that these adaptational processes are only part of the complex of processes leading to species formation, among which geographical isolation is probably also important.

Sexual imprinting leads to speciation in locally adapted populations

Sexual imprinting is widespread in birds and other species but its existence requires explanation. Our results suggest that sexual imprinting leads to speciation in locally-adapted populations if a neutral mating cue-e.g., novel plumage coloration-arises through mutation. Importantly, the mating cue locus is not linked to adaptation loci. Local adaptation is a necessary precursor to speciation and occurs when evolution results in stable genetic polymorphisms with one allele predominating in some areas while others predominate elsewhere. Here we use a deterministic two-niche population genetic model to map the set of migration and selection rates for which polymorphic evolutionary outcomes, i.e., local adaptations, can occur. Approximate equations for the boundaries of the set of polymorphic evolutionary outcomes were derived by Bulmer (American Naturalist, 106, 254, 1972), but our results, obtained by deterministic simulation of the evolutionary process, show that one of Bulmer's equations is inaccurate except when the level of dominance is 0.5, and fails if one of the alleles is dominant. Having an accurate map of the set of migration and selection rates for which polymorphic evolutionary outcomes can occur, we then show using the model of Sibly et al. (Ecology and Evolution, 9, 13506, 2019) that local adaptation in all analyzed cases leads to speciation if a new neutral mating cue arises by mutation. We finish by considering how genome sequencing makes possible testing our model and its predictions.

The terroir of the finch: How spatial and temporal variation shapes phenotypic traits in DARWIN'S finches

The term terroir is used in viticulture to emphasize how the biotic and abiotic characteristics of a local site influence grape physiology and thus the properties of wine. In ecology and evolution, such terroir (i.e., the effect of space or "site") is expected to play an important role in shaping phenotypic traits. Just how important is the pure spatial effect of terroir (e.g., differences between sites that persist across years) in comparison to temporal variation (e.g., differences between years that persist across sites), and the interaction between space and time (e.g., differences between sites change across years)? We answer this question by analyzing beak and body traits of 4388 medium ground finches (Geospiza fortis) collected across 10 years at three locations in Galápagos. Analyses of variance indicated that phenotypic variation was mostly explained by site for beak size (η 2 = 0.42) and body size (η 2 = 0.43), with a smaller contribution for beak shape (η 2 = 0.05) and body shape (η 2 = 0.12), but still higher compared to year and site-by-year effects. As such, the effect of terroir seems to be very strong in Darwin's finches, notwithstanding the oft-emphasized interannual variation. However, these results changed dramatically when we excluded data from Daphne Major, indicating that the strong effect of terroir was mostly driven by that particular population. These phenotypic results were largely paralleled in analyses of environmental variables (rainfall and vegetation indices) expected to shape terroir in this system. These findings affirm the evolutionary importance of terroir, while also revealing its dependence on other factors, such as geographical isolation.

Where did the finch go? Insights from radio telemetry of the medium ground finch (Geospiza fortis)

Movement patterns and habitat selection of animals have important implications for ecology and evolution. Darwin's finches are a classic model system for ecological and evolutionary studies, yet their spatial ecology remains poorly studied. We tagged and radio-tracked five (three females, two males) medium ground finches (Geospiza fortis) to examine the feasibility of telemetry for understanding their movement and habitat use. Based on 143 locations collected during a 3-week period, we analyzed for the first time home-range size and habitat selection patterns of finches at El Garrapatero, an arid coastal ecosystem on Santa Cruz Island (Galápagos). The average 95% home range and 50% core area for G. fortis in the breeding season was 20.54 ha ± 4.04 ha SE and 4.03 ha ± 1.11 ha SE, respectively. For most of the finches, their home range covered a diverse set of habitats. Three finches positively selected the dry-forest habitat, while the other habitats seemed to be either negatively selected or simply neglected by the finches. In addition, we noted a communal roosting behavior in an area close to the ocean, where the vegetation is greener and denser than the more inland dry-forest vegetation. We show that telemetry on Darwin's finches provides valuable data to understand the movement ecology of the species. Based on our results, we propose a series of questions about the ecology and evolution of Darwin's finches that can be addressed using telemetry.

Darwin's small and medium ground finches might have taste preferences, but not for human foods

Urbanization is rapidly changing ecological niches. On the inhabited Galapagos Islands, Darwin's finches consume human-introduced foods preferentially; however, it remains unclear why. Here, we presented pastry with flavour profiles typical of human foods (oily, salty and sweet) to small ground finches (Geospiza fuliginosa) and medium ground finches (Geospiza fortis) to test if latent taste preferences might drive the selection of human foods. If human food flavours were consumed more than a neutral or bitter control only at sites with human foods, then we predicted tastes were acquired after urbanization; however, if no site differences were found then this would indicate latent taste preferences. Contrary to both predictions, we found little evidence that human food flavours were preferred compared with control flavours at any site. Instead, finches showed a weak aversion to oily foods, but only at remote (no human foods present) sites. This was further supported by behavioural responses, with beak-wiping occurring more often at remote sites after finches tasted flavours associated with human foods. Our results suggest, therefore, that while Darwin's finches regularly exposed to human foods might have acquired a tolerance to human food flavours, latent taste preferences are unlikely to have played a major role in their dietary response to increased urbanization.

Morphological ghosts of introgression in Darwin's finch populations

Many species of plants, animals, and microorganisms exchange genes well after the point of evolutionary divergence at which taxonomists recognize them as species. Genomes contain signatures of past gene exchange and, in some cases, they reveal a legacy of lineages that no longer exist. But genomic data are not available for many organisms, and particularly problematic for reconstructing and interpreting evolutionary history are communities that have been depleted by extinctions. For these, morphology may substitute for genes, as exemplified by the history of Darwin's finches on the Galápagos islands of Floreana and San Cristóbal. Darwin and companions collected seven specimens of a uniquely large form of Geospiza magnirostris in 1835. The populations became extinct in the next few decades, partly due to destruction of Opuntia cactus by introduced goats, whereas Geospiza fortis has persisted to the present. We used measurements of large samples of G. fortis collected for museums in the period 1891 to 1906 to test for unusually large variances and skewed distributions of beak and body size resulting from introgression. We found strong evidence of hybridization on Floreana but not on San Cristóbal. The skew is in the direction of the absent G. magnirostris We estimate introgression influenced 6% of the frequency distribution that was eroded by selection after G. magnirostris became extinct on these islands. The genetic residuum of an extinct species in an extant one has implications for its future evolution, as well as for a conservation program of reintroductions in extinction-depleted communities.

Intraspecific genomic variation and local adaptation in a young hybrid species

Hybridization increases genetic variation, hence hybrid species may have greater evolutionary potential once their admixed genomes have stabilized and incompatibilities have been purged. Yet, little is known about how such hybrid lineages evolve at the genomic level following their formation, in particular their adaptive potential. Here we investigate how the Italian sparrow (Passer italiae), a homoploid hybrid species, has evolved and locally adapted to its variable environment. Using restriction site-associated DNA sequencing (RAD-seq) on several populations across the Italian peninsula, we evaluate how genomic constraints and novel genetic variation have influenced population divergence and adaptation. We show that population divergence within this hybrid species has evolved in response to climatic variation, suggesting ongoing local adaptation. As found previously in other nonhybrid species, climatic differences appear to increase population differentiation. We also report strong population divergence in a gene known to affect beak morphology. Most of the strongly divergent loci among Italian sparrow populations do not seem to be differentiated between its parent species, the house and Spanish sparrows. Unlike in the hybrid, population divergence within each of the parental taxa has occurred mostly at loci with high allele frequency difference between the parental species, suggesting that novel combinations of parental alleles in the hybrid have not necessarily enhanced its evolutionary potential. Rather, our study suggests that constraints linked to incompatibilities may have restricted the evolution of this admixed genome, both during and after hybrid species formation.

Unexpected diversity and evolutionary lability in root architectural ecomorphs in the rushes of the hyperdiverse Cape flora

Plants use roots to access soil resources, so differences in root traits and their ecological consequences could be a mechanism of species coexistence and niche divergence. Current views of the evolution of root diversity are informed by large-scale evolutionary analyses based on taxonomically coarse sampling and led to the 'root trait phylogenetic conservatism hypothesis'. Here we test this hypothesised conservatism among closely related species, and whether root variation plays an ecological role. We collected root architectural traits for the species-rich Cape rushes (Restionaceae) in the field and from herbaria. We used machine learning to interpolate missing data. Using model-based clustering we classified root syndromes. We modelled the proportion of the syndromes along environmental gradients using assemblages and environmental data of 735 plots. We fitted trait evolutionary models to test for the conservatism hypothesis. We recognised five root syndromes. Responses to environmental gradients are syndrome specific and thus these represent ecomorphs. Trait evolutionary models reveal an evolutionary lability in these ecomorphs. This could present the mechanistic underpinning of the taxonomic radiation of this group which has been linked to repeated habitat shifts. Our results challenge the perspective of strong phylogenetic conservatism and root trait evolution may more generally drive diversification.

Interdependent Phenotypic and Biogeographic Evolution Driven by Biotic Interactions

Biotic interactions are hypothesized to be one of the main processes shaping trait and biogeographic evolution during lineage diversification. Theoretical and empirical evidence suggests that species with similar ecological requirements either spatially exclude each other, by preventing the colonization of competitors or by driving coexisting populations to extinction, or show niche divergence when in sympatry. However, the extent and generality of the effect of interspecific competition in trait and biogeographic evolution has been limited by a dearth of appropriate process-generating models to directly test the effect of biotic interactions. Here, we formulate a phylogenetic parametric model that allows interdependence between trait and biogeographic evolution, thus enabling a direct test of central hypotheses on how biotic interactions shape these evolutionary processes. We adopt a Bayesian data augmentation approach to estimate the joint posterior distribution of trait histories, range histories, and coevolutionary process parameters under this analytically intractable model. Through simulations, we show that our model is capable of distinguishing alternative scenarios of biotic interactions. We apply our model to the radiation of Darwin’s finches—a classic example of adaptive divergence—and find limited support for in situ trait divergence in beak size, but stronger evidence for convergence in traits such as beak shape and tarsus length and for competitive exclusion throughout their evolutionary history. These findings are more consistent with presympatric, rather than postsympatric, niche divergence. Our modeling framework opens new possibilities for testing more complex hypotheses about the processes underlying lineage diversification. More generally, it provides a robust probabilistic methodology to model correlated evolution of continuous and discrete characters. [Bayesian; biotic interactions; competition; data augmentation; historical biogeography; trait evolution.]

Host phylogeny, diet, and habitat differentiate the gut microbiomes of Darwin's finches on Santa Cruz Island

Darwin's finches are an iconic example of an adaptive radiation with well-characterized evolutionary history, dietary preferences, and biogeography, offering an unparalleled opportunity to disentangle effects of evolutionary history on host microbiome from other factors like diet and habitat. Here, we characterize the gut microbiome in Darwin's finches, comparing nine species that occupy diverse ecological niches on Santa Cruz island. The finch phylogeny showed moderate congruence with the microbiome, which was comprised mostly of the bacterial phyla Firmicutes, Actinobacteria, and Proteobacteria. Diet, as measured with stable isotope values and foraging observations, also correlated with microbiome differentiation. Additionally, each gut microbial community could easily be classified by the habitat of origin independent of host species. Altogether, these findings are consistent with a model of microbiome assembly in which environmental filtering via diet and habitat are primary determinants of the bacterial taxa present with lesser influence from the evolutionary history between finch species.

Urbanization erodes niche segregation in Darwin's finches

Urbanization is influencing patterns of biological evolution in ways that are only beginning to be explored. One potential effect of urbanization is in modifying ecological resource distributions that underlie niche differences and that thus promote and maintain species diversification. Few studies have assessed such modifications, or their potential evolutionary consequences, in the context of ongoing adaptive radiation. We study this effect in Darwin's finches on the Galápagos Islands, by quantifying feeding preferences and diet niche partitioning across sites with different degrees of urbanization. We found higher finch density in urban sites and that feeding preferences and diets at urban sites skew heavily toward human food items. Furthermore, we show that finches at urban sites appear to be accustomed to the presence of people, compared with birds at sites with few people. In addition, we found that human behavior via the tendency to feed birds at non-urban but tourist sites is likely an important driver of finch preferences for human foods. Site differences in diet and feeding behavior have resulted in larger niche breadth within finch species and wider niche overlap between species at the urban sites. Both factors effectively minimize niche differences that would otherwise facilitate interspecies coexistence. These findings suggest that both human behavior and ongoing urbanization in Galápagos are starting to erode ecological differences that promote and maintain adaptive radiation in Darwin's finches. Smoothing of adaptive landscapes underlying diversification represents a potentially important yet underappreciated consequence of urbanization. Overall, our findings accentuate the fragility of the initial stages of adaptive radiation in Darwin's finches and raise concerns about the fate of the Galápagos ecosystems in the face of increasing urbanization.

Population genomics fits the bill: genetics of adaptive beak variation in Darwin's finches

Darwin's finches are an iconic case of adaptive radiation. The size and shape of their beaks are key adaptive traits related to trophic niche that vary among species and evolve rapidly when the food supply changes. Building on recent studies, a paper in this issue of Molecular Ecology (Chaves et al. ) investigates the genomic basis of beak size variation in sympatric populations of three species of ground finch (Geospiza) by performing a Genome-wide association study using RAD-seq data. The authors find that variation in a small number of markers can explain a substantial proportion of variation in beak size. Some of these markers are in genomic regions that have previously been implicated in beak size variation in Darwin's finches, whereas other markers have not, suggesting both conservation and divergence in the genetic basis of morphological evolution. Overall, the study confirms that loci of large effect are involved in beak size variation, which helps to explain the high heritability and rapid response to selection of this trait. The independent identification of regions containing HMGA2 and DLK1 loci in a GWAS makes them prime targets for functional studies. The study also shows that under the right conditions, RAD-seq can be a viable alternative to genome sequencing for GWAS in wild vertebrate populations.

Evolution within a language: environmental differences contribute to divergence of dialect groups

Background

The processes leading to the diversity of over 7000 present-day languages have been the subject of scholarly interest for centuries. Several factors have been suggested to contribute to the spatial segregation of speaker populations and the subsequent linguistic divergence. However, their formal testing and the quantification of their relative roles is still missing. We focussed here on the early stages of the linguistic divergence process, that is, the divergence of dialects, with a special focus on the ecological settings of the speaker populations. We adopted conceptual and statistical approaches from biological microevolution and parallelled intra-lingual variation with genetic variation within a species. We modelled the roles of geographical distance, differences in environmental and cultural conditions and in administrative history on linguistic divergence at two different levels: between municipal dialects (cf. in biology, between individuals) and between dialect groups (cf. in biology, between populations).

Results

We found that geographical distance and administrative history were important in separating municipal dialects. However, environmental and cultural differences contributed markedly to the divergence of dialect groups. In biology, increase in genetic differences between populations together with environmental differences may suggest genetic differentiation of populations through adaptation to the local environment. However, our interpretation of this result is not that language itself adapts to the environment. Instead, it is based on Homo sapiens being affected by its environment, and its capability to adapt culturally to various environmental conditions. The differences in cultural adaptations arising from environmental heterogeneity could have acted as nonphysical barriers and limited the contacts and communication between groups. As a result, linguistic differentiation may emerge over time in those speaker populations which are, at least partially, separated.

Conclusions

Given that the dialects of isolated speaker populations may eventually evolve into different languages, our result suggests that cultural adaptation to local environment and the associated isolation of speaker populations have contributed to the emergence of the global patterns of linguistic diversity.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1238-6) contains supplementary material, which is available to authorized users.

Evolution of novel mimicry rings facilitated by adaptive introgression in tropical butterflies

Understanding the genetic basis of phenotypic variation and the mechanisms involved in the evolution of adaptive novelty, especially in adaptive radiations, is a major goal in evolutionary biology. Here, we used whole-genome sequence data to investigate the origin of the yellow hindwing bar in the Heliconius cydno radiation. We found modular variation associated with hindwing phenotype in two narrow noncoding regions upstream and downstream of the cortex gene, which was recently identified as a pigmentation pattern controller in multiple species of Heliconius. Genetic variation at each of these modules suggests an independent control of the dorsal and ventral hindwing patterning, with the upstream module associated with the ventral phenotype and the downstream module with the dorsal one. Furthermore, we detected introgression between H. cydno and its closely related species Heliconius melpomene in these modules, likely allowing both species to participate in novel mimicry rings. In sum, our findings support the role of regulatory modularity coupled with adaptive introgression as an elegant mechanism by which novel phenotypic combinations can evolve and fuel an adaptive radiation.

Epigenetic variation between urban and rural populations of Darwin's finches

Background

The molecular basis of evolutionary change is assumed to be genetic variation. However, growing evidence suggests that epigenetic mechanisms, such as DNA methylation, may also be involved in rapid adaptation to new environments. An important first step in evaluating this hypothesis is to test for the presence of epigenetic variation between natural populations living under different environmental conditions.

Results

In the current study we explored variation between populations of Darwin’s finches, which comprise one of the best-studied examples of adaptive radiation. We tested for morphological, genetic, and epigenetic differences between adjacent “urban” and “rural” populations of each of two species of ground finches, Geospiza fortis and G. fuliginosa, on Santa Cruz Island in the Galápagos. Using data collected from more than 1000 birds, we found significant morphological differences between populations of G. fortis, but not G. fuliginosa. We did not find large size copy number variation (CNV) genetic differences between populations of either species. However, other genetic variants were not investigated. In contrast, we did find dramatic epigenetic differences between the urban and rural populations of both species, based on DNA methylation analysis. We explored genomic features and gene associations of the differentially DNA methylated regions (DMR), as well as their possible functional significance.

Conclusions

In summary, our study documents local population epigenetic variation within each of two species of Darwin’s finches.

Electronic supplementary material

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

Absence of population structure across elevational gradients despite large phenotypic variation in mountain chickadees (Poecile gambeli)

Montane habitats are characterized by predictably rapid heterogeneity along elevational gradients and are useful for investigating the consequences of environmental heterogeneity for local adaptation and population genetic structure. Food-caching mountain chickadees inhabit a continuous elevation gradient in the Sierra Nevada, and birds living at harsher, high elevations have better spatial memory ability and exhibit differences in male song structure and female mate preference compared to birds inhabiting milder, low elevations. While high elevation birds breed, on average, two weeks later than low elevation birds, the extent of gene flow between elevations is unknown. Despite phenotypic variation and indirect evidence for local adaptation, population genetic analyses based on 18 073 single nucleotide polymorphisms across three transects of high and low elevation populations provided no evidence for genetic differentiation. Analyses based on individual genotypes revealed no patterns of clustering, pairwise estimates of genetic differentiation (F_ST, Nei's D) were very low, and AMOVA revealed no evidence for genetic variation structured by transect or by low and high elevation sites within transects. In addition, we found no consistent evidence for strong parallel allele frequency divergence between low and high elevation sites within the three transects. Large elevation-related phenotypic variation may be maintained by strong selection despite gene flow and future work should focus on the mechanisms underlying such variation.

Ongoing Speciation and Gene Flow between Taxonomically Challenging Trochulus Species Complex (Gastropoda: Hygromiidae)

Geographical isolation, selection and genetic drift can cause the geographical diversification of populations and lead to speciation. Land snail species in the genus Trochulus show overlaps in geographical ranges as well as in morphology, but genetic data do not always support the species-level taxonomy based on morphological characters. Such a group offers an excellent opportunity to explore the processes involved. We have addressed the problem by determining the status of the restricted endemic T. graminicola within the larger context of Trochulus taxonomy. We used an integrated approach based on morphological features, ecological preferences and two molecular markers: mitochondrial COI sequences and microsatellites. Comparison of these results demonstrated: (i) conchological distinction of T. striolatus and T. sericeus; (ii) anatomical, ecological and genetic differentiation of T. graminicola and (iii) concordance between morphological characters and mtDNA markers in T. striolatus. Moreover, our data showed an intricate evolutionary history within the genus Trochulus, which can be best explained by: (i) recent or ongoing gene flow between taxa or (ii) their large ancestral polymorphism. Both of these hypotheses suggest that diversification within this group of snails has occurred relatively recently. The mismatches between species defined on morphology and on molecular genetics indicate the complexity of the processes involved in the diversification of this genus.

Patterns of gene flow and selection across multiple species of Acrocephalus warblers: footprints of parallel selection on the Z chromosome

Background

Understanding the mechanisms and selective forces leading to adaptive radiations and origin of biodiversity is a major goal of evolutionary biology. Acrocephalus warblers are small passerines that underwent an adaptive radiation in the last approximately 10 million years that gave rise to 37 extant species, many of which still hybridize in nature. Acrocephalus warblers have served as model organisms for a wide variety of ecological and behavioral studies, yet our knowledge of mechanisms and selective forces driving their radiation is limited. Here we studied patterns of interspecific gene flow and selection across three European Acrocephalus warblers to get a first insight into mechanisms of radiation of this avian group.

Results

We analyzed nucleotide variation at eight nuclear loci in three hybridizing Acrocephalus species with overlapping breeding ranges in Europe. Using an isolation-with-migration model for multiple populations, we found evidence for unidirectional gene flow from A. scirpaceus to A. palustris and from A. palustris to A. dumetorum. Gene flow was higher between genetically more closely related A. scirpaceus and A. palustris than between ecologically more similar A. palustris and A. dumetorum, suggesting that gradual accumulation of intrinsic barriers rather than divergent ecological selection are more efficient in restricting interspecific gene flow in Acrocephalus warblers. Although levels of genetic differentiation between different species pairs were in general not correlated, we found signatures of apparently independent instances of positive selection at the same two Z-linked loci in multiple species.

Conclusions

Our study brings the first evidence that gene flow occurred during Acrocephalus radiation and not only between sister species. Interspecific gene flow could thus be an important source of genetic variation in individual Acrocephalus species and could have accelerated adaptive evolution and speciation rate in this avian group by creating novel genetic combinations and new phenotypes. Independent instances of positive selection at the same loci in multiple species indicate an interesting possibility that the same loci might have contributed to reproductive isolation in several speciation events.

Electronic supplementary material

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

Nuclear introgression without mitochondrial introgression in two turtle species exhibiting sex-specific trophic differentiation

Despite the presence of reproductive barriers between species, interspecific gene introgression has been documented in a range of natural systems. Comparing patterns of genetic introgression in biparental versus matrilineal markers can potentially reveal sex-specific barriers to interspecific gene flow. Hybridization has been documented in the freshwater turtles Graptemys geographica and G. pseudogeographica, whose ranges are largely sympatric. Morphological differentiation between the species is restricted to females, with female G. geographica possessing large heads and jaws compared to the narrow heads of G. pseudogeographica females. If hybrid females are morphologically intermediate, they may be less successful at exploiting parental feeding niches, thereby limiting the introgression of maternally inherited, but not biparental, molecular markers. We paired sequence data with stable isotope analysis and examined sex-specific genetic introgression and trophic differentiation in sympatric populations of G. geographica and G. pseudogeographica. We observed introgression from G. pseudogeographica into G. geographica at three nuclear loci, but not at the mitochondrial locus. Analysis of ∂(15)N and ∂(13)C was consistent with species differences in trophic positioning in females, but not males. These results suggest that ecological divergence in females may reduce the opportunity for gene flow in this system.

Population genomics of local adaptation versus speciation in coral reef fishes (Hypoplectrus spp, Serranidae)

Are the population genomic patterns underlying local adaptation and the early stages of speciation similar? Addressing this question requires a system in which (i) local adaptation and the early stages of speciation can be clearly identified and distinguished, (ii) the amount of genetic divergence driven by the two processes is similar, and (iii) comparisons can be repeated both taxonomically (for local adaptation) and geographically (for speciation). Here, we report just such a situation in the hamlets (Hypoplectrus spp), brightly colored reef fishes from the wider Caribbean. Close to 100,000 SNPs genotyped in 126 individuals from three sympatric species sampled in three repeated populations provide genome-wide levels of divergence that are comparable among allopatric populations (F st estimate = 0.0042) and sympatric species (F st estimate = 0.0038). Population genetic, clustering, and phylogenetic analyses reveal very similar patterns for local adaptation and speciation, with a large fraction of the genome undifferentiated (F st estimate ≈ 0), a very small proportion of F st outlier loci (0.05-0.07%), and remarkably few repeated outliers (1-3). Nevertheless, different loci appear to be involved in the two processes in Hypoplectrus, with only 7% of the most differentiated SNPs and outliers shared between populations and species comparisons. In particular, a tropomyosin (Tpm4) and a previously identified hox (HoxCa) locus emerge as candidate loci (repeated outliers) for local adaptation and speciation, respectively. We conclude that marine populations may be locally adapted notwithstanding shallow levels of genetic divergence, and that from a population genomic perspective, this process does not appear to differ fundamentally from the early stages of speciation.

Adult survival selection in relation to multilocus heterozygosity and body size in a tropical bird species, the Zenaida dove, Zenaida aurita

Both phenotypic and genetic quality can influence the survival of individuals through time, although their relative influences are rarely addressed simultaneously. Here we used capture-mark-recapture modelling to assess the influence of both multilocus heterozygosity (MLH) and body size on apparent adult survival in a tropical bird species, the Zenaida dove, Zenaida aurita, using a sample of 391 individuals genotyped at 11 microsatellites, while controlling for the effects of sex. No effect of body size on either adult survival or capture rate was found. In the best model, survival was a logit linear function of MLH, whereas detection probability was a sex-dependent logit linear function of the logarithm of field effort, increasing with time and affected by a random individual effect. Using a Bayesian approach, we found that MLH explained 1.14% of the total deviance, as expected from theory and previous studies of heterozygosity-fitness correlations, with no evidence for local effects. However, results from capture-mark-recapture modelling indicated that expected longevity varied from 4.8 years in the least heterozygous individuals (MLH = 0.37) to 10.6 years in the most heterozygous ones (MLH = 1), thus suggesting that MLH had potentially a substantial effect on survival. We discuss our results in relation to current hypotheses about the origin of heterozygosity-fitness correlations.

Environment, but not genetic divergence, influences geographic variation in colour morph frequencies in a lizard

Background

Identifying the causes of intraspecific phenotypic variation is essential for understanding evolutionary processes that maintain diversity and promote speciation. In polymorphic species, the relative frequencies of discrete morphs often vary geographically; yet the drivers of spatial variation in morph frequencies are seldom known. Here, we test the relative importance of gene flow and natural selection to identify the causes of geographic variation in colour morph frequencies in the Australian tawny dragon lizard, Ctenophorus decresii.

Results

Populations of C. decresii are polymorphic for male throat coloration and all populations surveyed shared the same four morphs but differed in the relative frequencies of morphs. Despite genetic structure among populations, there was no relationship between genetic similarity or geographic proximity and similarity in morph frequencies. However, we detected remarkably strong associations between morph frequencies and two environmental variables (mean annual aridity index and vegetation cover), which together explained approximately 45 % of the total variance in morph frequencies.

Conclusions

Spatial variation in selection appears to play an important role in shaping morph frequency patterns in C. decresii. Selection associated with differences in local environmental conditions, combined with relatively low levels of gene flow, is expected to favour population divergence in morph composition, but may be counteracted by negative frequency-dependent selection favouring rare morphs.

Electronic supplementary material

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

Landscape genetic connectivity in a riparian foundation tree is jointly driven by climatic gradients and river networks

Fremont cottonwood (Populus fremonti) is a foundation riparian tree species that drives community structure and ecosystem processes in southwestern U.S. ecosystems. Despite its ecological importance, little is known about the ecological and environmental processes that shape its genetic diversity, structure, and landscape connectivity. Here, we combined molecular analyses of 82 populations including 1312 individual trees dispersed over the species' geographical distribution. We reduced the data set to 40 populations and 743 individuals to eliminate admixture with a sibling species, and used multivariate restricted optimization and reciprocal causal modeling to evaluate the effects of river network connectivity and climatic gradients on gene flow. Our results confirmed the following: First, gene flow of Fremont cottonwood is jointly controlled by the connectivity of the river network and gradients of seasonal precipitation. Second, gene flow is facilitated by mid-sized to large rivers, and is resisted by small streams and terrestrial uplands, with resistance to gene flow decreasing with river size. Third, genetic differentiation increases with cumulative differences in winter and spring precipitation. Our results suggest that ongoing fragmentation of riparian habitats will lead to a loss of landscape-level genetic connectivity, leading to increased inbreeding and the concomitant loss of genetic diversity in a foundation species. These genetic effects will cascade to a much larger community of organisms, some of which are threatened and endangered.

Darwin's finches and their diet niches: the sympatric coexistence of imperfect generalists

Adaptive radiation can be strongly influenced by interspecific competition for resources, which can lead to diverse outcomes ranging from competitive exclusion to character displacement. In each case, sympatric species are expected to evolve into distinct ecological niches, such as different food types, yet this expectation is not always met when such species are examined in nature. The most common hypotheses to account for the coexistence of species with substantial diet overlap rest on temporal variation in niches (often diets). Yet spatial variation in niche overlap might also be important, pointing to the need for spatiotemporal analyses of diet and diet overlap between closely related species persisting in sympatry. We here perform such an analysis by characterizing the diets of, and diet overlap among, four sympatric Darwin's ground finch species at three sites and over 5 years on a single Galápagos island (Santa Cruz). We find that the different species have broadly similar and overlapping diets - they are to some extent generalists and opportunists - yet we also find that each species retains some 'private' resources for which their morphologies are best suited. Importantly, use of these private resources increased considerably, and diet overlap decreased accordingly, when the availability of preferred shared foods, such as arthropods, was reduced during drought conditions. Spatial variation in food resources was also important. These results together suggest that the ground finches are 'imperfect generalists' that use overlapping resources under benign conditions (in space or time), but then retreat to resources for which they are best adapted during periods of food limitation. These conditions likely promote local and regional coexistence.

Geographic variation in animal colour polymorphisms and its role in speciation

Polymorphic species, in which multiple variants coexist within a population, are often used as model systems in evolutionary biology. Recent research has been dominated by the hypothesis that polymorphism can be a precursor to speciation. To date, the majority of research regarding polymorphism and speciation has focused on whether polymorphism is maintained within a population or whether morphs within populations may diverge to form separate species (sympatric speciation); however, the geographical context of speciation in polymorphic systems is likely to be both diverse and complex. In this review, we draw attention to the geographic variation in morph composition and frequencies that characterises many, if not most polymorphic species. Recent theoretical and empirical developments suggest that such variation in the number, type and frequency of morphs present among populations can increase the probability of speciation. Thus, the geographical context of a polymorphism requires a greater research focus. Here, we review the prevalence, causes and evolutionary consequences of geographic variation in polymorphism in colour-polymorphic animal species. The prevalence and nature of geographic variation in polymorphism suggests that polymorphism may be a precursor to and facilitate speciation more commonly than appreciated previously. We argue that a better understanding of the processes generating geographic variation in polymorphism is vital to understanding how polymorphism can promote speciation.

Regional differences in seasonal timing of rainfall discriminate between genetically distinct East African giraffe taxa

Masai (Giraffa tippelskirchi), Reticulated (G. reticulata) and Rothschild's (G. camelopardalis) giraffe lineages in East Africa are morphologically and genetically distinct, yet in Kenya their ranges abut. This raises the question of how divergence is maintained among populations of a large mammal capable of long-distance travel, and which readily hybridize in zoos. Here we test four hypotheses concerning the maintenance of the phylogeographic boundaries among the three taxa: 1) isolation-by-distance; 2) physical barriers to dispersal; 3) general habitat differences resulting in habitat segregation; or 4) regional differences in the seasonal timing of rainfall, and resultant timing of browse availability. We used satellite remotely sensed and climate data to characterize the environment at the locations of genotyped giraffes. Canonical variate analysis, random forest algorithms, and generalized dissimilarity modelling were employed in a landscape genetics framework to identify the predictor variables that best explained giraffes' genetic divergence. We found that regional differences in the timing of precipitation, and resulting green-up associated with the abundance of browse, effectively discriminate between taxa. Local habitat conditions, topographic and human-induced barriers, and geographic distance did not aid in discriminating among lineages. Our results suggest that selection associated with regional timing of events in the annual climatic cycle may help maintain genetic and phenotypic divergence in giraffes. We discuss potential mechanisms of maintaining divergence, and suggest that synchronization of reproduction with seasonal rainfall cycles that are geographically distinct may contribute to reproductive isolation. Coordination of weaning with green-up cycles could minimize the costs of lactation and predation on the young. Our findings are consistent with theory and empirical results demonstrating the efficacy of seasonal or phenologically dictated selection pressures in contributing to the reproductive isolation of parapatric populations.

Multilocus analysis of divergence and introgression in sympatric and allopatric sibling species of the Lutzomyia longipalpis complex in Brazil

BACKGROUND

Lutzomyia longipalpis, the main vector of visceral leishmaniasis in Latin America, is a complex of sibling species. In Brazil, a number of very closely related sibling species have been revealed by the analyses of copulation songs, sex pheromones and molecular markers. However, the level of divergence and gene flow between the sibling species remains unclear. Brazilian populations of this vector can be divided in two main groups: one producing Burst-type songs and the Cembrene-1 pheromone and a second more diverse group producing various Pulse song subtypes and different pheromones.

METHODOLOGY/PRINCIPAL FINDINGS

We analyzed 21 nuclear loci in two pairs of Brazilian populations: two sympatric populations from the Sobral locality (1S and 2S) in northeastern Brazil and two allopatric populations from the Lapinha and Pancas localities in southeastern Brazil. Pancas and Sobral 2S are populations of the Burst/Cembrene-1 species while Lapinha and Sobral 1S are two putative incipient species producing the same pheromone and similar Pulse song subtypes. The multilocus analysis strongly suggests the occurrence of gene flow during the divergence between the sibling species, with different levels of introgression between loci. Moreover, this differential introgression is asymmetrical, with estimated gene flow being higher in the direction of the Burst/Cembrene-1 species.

CONCLUSIONS/SIGNIFICANCE

The results indicate that introgressive hybridization has been a crucial phenomenon in shaping the genome of the L. longipalpis complex. This has possible epidemiological implications and is particularly interesting considering the potential for increased introgression caused by man-made environmental changes and the current trend of leishmaniasis urbanization in Brazil.

Extremely reduced dispersal and gene flow in an island bird

The Réunion grey white-eye, Zosterops borbonicus, a passerine bird endemic to Réunion Island in the Mascarene archipelago, represents an extreme case of microgeographical plumage colour variation in birds, with four distinct colour forms occupying different parts of this small island (2512 km(2)). To understand whether such population differentiation may reflect low levels of dispersal and gene flow at a very small spatial scale, we examined population structure and gene flow by analysing variation at 11 microsatellite loci among four geographically close localities (<26 km apart) sampled within the distribution range of one of the colour forms, the brown-headed brown form. Our results revealed levels of genetic differentiation that are exceptionally high for birds at such a small spatial scale. This strong population structure appears to reflect low levels of historical and contemporary gene flow among populations, unless very close geographically (<10 km). Thus, we suggest that the Réunion grey white-eye shows an extremely reduced propensity to disperse, which is likely to be related to behavioural processes.

Local adaptation within a hybrid species

Ecological divergence among populations may be strongly influenced by their genetic background. For instance, genetic admixture through introgressive hybridization or hybrid speciation is likely to affect the genetic variation and evolvability of phenotypic traits. We studied geographic variation in two beak dimensions and three other phenotypic traits of the Italian sparrow (Passer italiae), a young hybrid species formed through interbreeding between house sparrows (P. domesticus) and Spanish sparrows (P. hispaniolensis). We found that beak morphology was strongly influenced by precipitation regimes and that it appeared to be the target of divergent selection within Italian sparrows. Interestingly, however, the degree of parental genetic contribution in the hybrid species had no effect on phenotypic beak variation. Moreover, beak height divergence may mediate genetic differentiation between populations, consistent with isolation-by-adaptation within this hybrid species. The study illustrates how hybrid species may be relatively unconstrained by their admixed genetic background, allowing them to adapt rapidly to environmental variation.

Widespread evidence for incipient ecological speciation: a meta-analysis of isolation-by-ecology

Ecologically mediated selection has increasingly become recognised as an important driver of speciation. The correlation between neutral genetic differentiation and environmental or phenotypic divergence among populations, to which we collectively refer to as isolation-by-ecology (IBE), is an indicator of ecological speciation. In a meta-analysis framework, we determined the strength and commonality of IBE in nature. On the basis of 106 studies, we calculated a mean effect size of IBE with and without controlling for spatial autocorrelation among populations. Effect sizes were 0.34 (95% CI 0.24-0.42) and 0.26 (95% CI 0.13-0.37), respectively, indicating that an average of 5% of the neutral genetic differentiation among populations was explained purely by ecological contrast. Importantly, spatial autocorrelation reduced IBE correlations for environmental variables, but not for phenotypes. Through simulation, we showed how the influence of isolation-by-distance and spatial autocorrelation of ecological variables can result in false positives or underestimated correlations if not accounted for in the IBE model. Collectively, this meta-analysis showed that ecologically induced genetic divergence is pervasive across time-scales and taxa, and largely independent of the choice of molecular marker. We discuss the importance of these results in the context of adaptation and ecological speciation and suggest future research avenues.

Evidence for asymmetrical divergence-gene flow of nuclear loci, but not mitochondrial loci, between seabird sister species: blue-footed (Sula nebouxii) and Peruvian (S. variegata) boobies

Understanding the process of speciation requires understanding how gene flow influences divergence. Recent analyses indicate that divergence can take place despite gene flow and that the sex chromosomes can exhibit different levels of gene flow than autosomes and mitochondrial DNA. Using an eight marker dataset including autosomal, z-linked, and mitochondrial loci we tested the hypothesis that blue-footed (Sula nebouxii) and Peruvian (S. variegata) boobies diverged from their common ancestor with gene flow, paying specific attention to the differences in gene flow estimates from nuclear and mitochondrial markers. We found no gene flow at mitochondrial markers, but found evidence from the combined autosomal and z-linked dataset that blue-footed and Peruvian boobies experienced asymmetrical gene flow during or after their initial divergence, predominantly from Peruvian boobies into blue-footed boobies. This gene exchange may have occurred either sporadically between periods of allopatry, or regularly throughout the divergence process. Our results add to growing evidence that diverging species can remain distinct but exchange genes.

The use of the MegaBACE for sequencing and genotype analysis

Despite the advent of next generation sequencing techniques, which provide access to an enormous amount of genomic information in a relatively short time, the conventional Sanger sequencing and microsatellite genotyping analyses present a straightforward method to answer clearly defined questions in population genetics, phylogeography, or forensics. The MegaBACE is a platform that provides both applications with equally reliable performance. In this overview, protocols for the classical techniques of Sanger sequencing and microsatellite genotyping are described. This chapter aims to supply the user of the MegaBACE with methodological tools and some "insider" knowledge of this highly sensitive apparatus.