Convergent evolution of Darwin's finches caused by introgressive hybridization and selection

Abstract Between 1973 and 2003 mean morphological features of the cactus finch, Geospiza scandens, and the medium ground finch, G. fortis, populations on the Galápagos island of Daphne Major were subject to fluctuating directional selection. An increase in bluntness or robustness in the beak of G. scandens after 1990 can only partly be explained by selection. We use 16 microsatellite loci to test predictions of the previously proposed hypothesis that introgressive hybridization contributed to the trend, resulting in genes flowing predominantly from G. fortis to G. scandens. To identify F₁ hybrids and backcrosses we use pedigrees where known, supplemented by the results of assignment tests based on 14 autosomal loci when parents were not known. We analyze changes in morphology and allelic composition in the two populations over a period of 15–20 years. With samples that included F₁ hybrids and backcrosses, the G. scandens population became more similar to the G. fortis population both genetically and morphologically. Gene flow between species was estimated to be three times greater from G. fortis to G. scandens than in the opposite direction, resulting in a 20% reduction in the genetic difference between the species. Nevertheless, removing identified F₁ hybrids and backcrosses from the total sample and reanalyzing the traits did not eliminate the convergence. The two species also converged in beak shape by 22.2% and in body size by 45.5%. A combination of introgressive hybridization and selection jointly provide the best explanation of convergence in morphology and genetic constitution under the changed ecological conditions following a major El Niño event in 1983. The study illustrates how species without postmating barriers to gene exchange can alternate between convergence and divergence when environmental conditions oscillate.     

CULTURAL INHERITANCE OF SONG AND ITS ROLE IN THE EVOLUTION OF DARWIN'S FINCHES

Songs of Darwin's finches were studied on the Galápagos Island of Daphne Major from 1976 to 1995. A single, structurally simple, and unvarying song is sung throughout life by each male of the two common species, Geospiza fortis (medium ground finch) and G. scandens (cactus finch). Songs of the two species differ strongly in quantitative features, and individual variation among males is much broader in G. fortis than in G. scandens. Although there are exceptions, songs of sons strongly resemble the songs of their fathers. They also resemble the songs of their paternal grandfathers, but not their maternal grandfathers, indicating that they are culturally inherited and not genetically inherited. Female G. fortis display a tendency to avoid mating with males that sing the same type of song as their father. They also avoid mating with males that sing heterospecific song, with very rare exceptions. Thus song, an evolving, culturally inherited trait, is an important factor in species recognition and mate choice. It constrains the mating of females to conspecifics, even when there is no genetic penalty to interbreeding, and thus may play a crucial role in species formation by promoting genetic isolation on secondary contact. The barrier is leaky in that occasional errors in song transmission result in misimprinting, which leads to a low incidence of hybridization and introgression. Introgression slows the rate of postzygotic isolation, but can produce individuals in novel genetic and morphological space that can provide the starting point of a new evolutionary trajectory.     

Phenotypic variation and covariation indicate high evolvability of acoustic communication in crickets

Studying the genetic architecture of sexual traits provides insight into the rate and direction at which traits can respond to selection. Traits associated with few loci and limited genetic and phenotypic constraints tend to evolve at high rates typically observed for secondary sexual characters. Here, we examined the genetic architecture of song traits and female song preferences in the field crickets Gryllus rubens and Gryllus texensis. Song and preference data were collected from both species and interspecific F1 and F2 hybrids. We first analysed phenotypic variation to examine interspecific differentiation and trait distributions in parental and hybrid generations. Then, the relative contribution of additive and additive‐dominance variation was estimated. Finally, phenotypic variance–covariance ( P ) matrices were estimated to evaluate the multivariate phenotype available for selection. Song traits and preferences had unimodal trait distributions, and hybrid offspring were intermediate with respect to the parents. We uncovered additive and dominance variation in song traits and preferences. For two song traits, we found evidence for X‐linked inheritance. On the one hand, the observed genetic architecture does not suggest rapid divergence, although sex linkage may have allowed for somewhat higher evolutionary rates. On the other hand, P matrices revealed that multivariate variation in song traits aligned with major dimensions in song preferences, suggesting a strong selection response. We also found strong covariance between the main traits that are sexually selected and traits that are not directly selected by females, providing an explanation for the striking multivariate divergence in male calling songs despite limited divergence in female preferences.     

Reproductive isolation vs. interbreeding between Gagea lutea (L.) Ker Gawl. and G. pratensis (Pers.) Dumort. (Liliaceae) and their putative hybrids in Mecklenburg‐Western Pomerania (Germany)

Mecklenburg‐Western Pomerania is the main range of two hybridogeneous Gagea taxa, G. megapolitana Henker and G. pomeranica Ruthe, which have the same two parental species, G. lutea (L.) Ker. Gawl. and G. pratensis (Pers.) Dumort. We assessed the degree of reproductive isolation vs. interbreeding between these taxa using data from field observations, crossing experiments and pollen characteristics for nine hybrid, four G. lutea and five G. pratensis populations. Pollen viability was highest in 6x G. lutea and lowest in G. pratensis (most probably 5x), with intermediate figures for the studied hybrids (5x–7x). Despite the assumed anorthoploid states, sexual reproduction (though sometimes very rare) was recorded for all populations in the field and/or in experiments. The crossing experiments revealed that all taxa are also able to hybridize. However, there were differences in the directions (i.e. role as maternal vs. paternal parent) as well as the success of the crossings: the primary hybridization G. lutea x G. pratensis resulted in more seeds than the reverse combination (17.5 vs. 3.3%), but the seed set was highest in backcrosses of the hybrids with G. lutea pollen (41.2%). These differences can be explained by overlap of flowering times, and reduced fertilities due to ploidy levels. The study showed that the taxa of the G. lutea‐pratensis hybrid complex are not yet reproductively isolated but can interbreed and will generate hybrids of higher ranks, forming a hybrid swarm, most probably leading to introgression via backcrosses with G. lutea.     

Inbreeding and interbreeding in Darwin's finches

Studies of inbreeding and interspecific hybridization are generally pursued separately with different metrics. There is a need to integrate them because they have the common goal of seeking an understanding of the genetic and ecological basis of fitness variation in populations. We use mean expected heterozygosity as an axis of variation on which to compare the fitness of inbreeding and hybridizing Darwin's finches (Geospiza scandens and G. fortis) relative to the fitness of matched outbred controls. We find that relative fitness of inbred finches is less than one in the 1991 cohorts of both species. Inbreeding depression is stronger in the species (G. scandens) with the lower genetically effective population size. Relative fitness of hybrids (backcrosses) in the same cohort of G. scandens is greater than one. Evidence of heterosis in G. fortis is mixed. Thus the two interbreeding species displayed somewhat different fitness patterns under the same set of environmental conditions. Hybridization may enhance fitness to different degrees by counteracting the effects of inbreeding depression, by other additive and nonadditive genetic effects, and by producing phenotypes well suited to exploit particular ecological conditions.     

Genetic conservation and management of the California endemic,Torrey pine (Pinus torreyana Parry): Implications of genetic rescue in a genetically depauperate species

Rare species present a challenge under changing environmental conditions as the genetic consequences of rarity may limit species ability to adapt to environmental change. To evaluate the evolutionary potential of a rare species, we assessed variation in traits important to plant fitness using multigenerational common garden experiments. Torrey pine, Pinus torreyana Parry, is one of the rarest pines in the world, restricted to one mainland and one island population. Morphological differentiation between island and mainland populations suggests adaptation to local environments may have contributed to trait variation. The distribution of phenotypic variances within the common garden suggests distinct population‐specific growth trajectories underlay genetic differences, with the island population exhibiting substantially reduced genetic variance for growth relative to the mainland population. Furthermore, F1 hybrids, representing a cross between mainland and island trees, exhibit increased height accumulation and fecundity relative to mainland and island parents. This may indicate genetic rescue via intraspecific hybridization could provide the necessary genetic variation to persist in environments modified as a result of climate change. Long‐term common garden experiments, such as these, provide invaluable resources to assess the distribution of genetic variance that may inform conservation strategies to preserve evolutionary potential of rare species, including genetic rescue.     

PLEIOTROPY CAUSES LONG-TERM GENETIC CONSTRAINTS ON LIFE-HISTORY EVOLUTION IN BRASSICA RAPA

Fundamental, long-term genetic trade-offs constrain life-history evolution in wild crucifer populations. I studied patterns of genetic constraint in Brassica rapa by estimating genetic correlations among life-history components by quantitative genetic analyses among ten wild populations, and within four populations. Genetic correlations between age and size at first reproduction were always greater than +0.8 within and among all populations studied. Although quantitative genetics might provide insight about genetic constraints if genetic parameters remain approximately constant, little evidence has been available to determine the constancy of genetic correlations. I found strong and consistent estimates of genetic correlations between life-history components, which were very similar within four natural populations. Population differentiation also showed these same trade-offs, resulting from long-term genetic constraint. For some traits, evolutionary changes among populations were incompatible with a model of genetic drift. Historical patterns of natural selection were inferred from population differentiation, suggesting that correlated response to selection has caused some traits to evolve opposite to the direct forces of natural selection. Comparison with Arabidopsis suggests that these life-history trade-offs are caused by genes that regulate patterns of resource allocation to different components of fitness. Ecological and energetic models may correctly predict these trade-offs because there is little additive genetic variation for rates of resource acquisition, but resource allocation is genetically variable.     

Interspecific hybridization between common and tepary beans: increased hybrid embryo growth,fertility, and efficiency of hybridization through recurrent and congruity backcrossing

Cultivated common bean (Phaseolus vulgaris L.) and tepary bean (Phaseolus acutifolius A. Gray) genotypes possessing desirable agronomic traits were hybridized. The F₁ hybrids were backcrossed twice with the common bean (i.e., recurrent backcrossing). Also, alternate backcrosses with common and tepary beans (i.e., congruity backcrossing) were carried out. Embryo culture was necessary for all initial interspecific crosses, and its requirement was proportionally lower when the common bean was used as the recurrent parent and as the last parent of congruity backcrosses. Modification of the embryo culture technique was necessary to produce congruity hybrids. Effects of both tepary and common bean genotypes on the success rate of hybridization were observed. Tepary accession G 40001 and common bean cultivar ICA Pijao facilitated interspecies hybridization. Growth of hybrid embryos before rescue, recovery of mature hybrid plants, and the vigor and fertility of F₁ hybrids all increased with increased recurrent and congruity backcrosses and intermatings between male-sterile F₁ and selected fertile F₂ plants of the third and fifth congruity backcrosses. Introgression of tepary genes was verified by means of seed protein electrophoretic analysis and morphological markers. The results suggest that congruity backcrossing can help to gradually reduce or overcome P. vulgaris x P. acutifolius hybridization barriers such as genotype incompatibility, early embryo abortion, hybrid sterility, and lower frequencies of hybridization.     

Multivariate analysis of adaptive capacity for upper thermal limits in Drosophila simulans

Thermal tolerance is an important factor influencing the distribution of ectotherms, but our understanding of the ability of species to evolve different thermal limits is limited. Based on univariate measures of adaptive capacity, it has recently been suggested that species may have limited evolutionary potential to extend their upper thermal limits under ramping temperature conditions that better reflect heat stress in nature. To test these findings more broadly, we used a paternal half‐sibling breeding design to estimate the multivariate evolutionary potential for upper thermal limits in Drosophila simulans. We assessed heat tolerance using static (basal and hardened) and ramping assays. Our analyses revealed significant evolutionary potential for all three measures of heat tolerance. Additive genetic variances were significantly different from zero for all three traits. Our G matrix analysis revealed that all three traits would contribute to a response to selection for increased heat tolerance. Significant additive genetic covariances and additive genetic correlations between static basal and hardened heat‐knockdown time, marginally nonsignificant between static basal and ramping heat‐knockdown time, indicate that direct and correlated responses to selection for increased upper thermal limits are possible. Thus, combinations of all three traits will contribute to the evolution of upper thermal limits in response to selection imposed by a warming climate. Reliance on univariate estimates of evolutionary potential may not provide accurate insight into the ability of organisms to evolve upper thermal limits in nature.     

In situ marker‐based assessment of leaf trait evolutionary potential in a marginal European beech population

Evolutionary processes are expected to be crucial for the adaptation of natural populations to environmental changes. In particular, the capacity of rear edge populations to evolve in response to the species limiting conditions remains a major issue that requires to address their evolutionary potential. In situ quantitative genetic studies based on molecular markers offer the possibility to estimate evolutionary potentials manipulating neither the environment nor the individuals on which phenotypes are measured. The goal of this study was to estimate heritability and genetic correlations of a suite of leaf functional traits involved in climate adaptation for a natural population of the tree Fagus sylvatica, growing at the rear edge of the species range. Using two marker‐based quantitative genetics approaches, we obtained consistent and significant estimates of heritability for leaf phenological (phenology of leaf flush), morphological (mass, area, ratio mass/area) and physiological (δ¹³C, nitrogen content) traits. Moreover, we found only one significant positive genetic correlation between leaf area and leaf mass, which likely reflected mechanical constraints. We conclude first that the studied population has considerable genetic diversity for important ecophysiological traits regarding drought adaptation and, second, that genetic correlations are not likely to impose strong genetic constraints to future population evolution. Our results bring important insights into the question of the capacity of rear edge populations to evolve.     

Development of SNP markers identifying European wildcats,domestic cats,and their admixed progeny

Introgression can be an important evolutionary force but it can also lead to species extinction and as such is a crucial issue for species conservation. However, introgression is difficult to detect, morphologically as well as genetically. Hybridization with domestic cats (Felis silvestris catus) is a major concern for the conservation of European wildcats (Felis s. silvestris). The available morphologic and genetic markers for the two Felis subspecies are not sufficient to reliably detect hybrids beyond first generation. Here we present a single nucleotide polymorphism (SNP) based approach that allows the identification of introgressed individuals. Using high‐throughput sequencing of reduced representation libraries we developed a diagnostic marker set containing 48 SNPs (Fst  > 0.8) which allows the identification of wildcats, domestic cats, their hybrids and backcrosses. This allows assessing introgression rate in natural wildcat populations and is key for a better understanding of hybridization processes.     

Hybridization and asymmetric introgression between Cypripedium tibeticum and C. yunnanense in Shangrila County,Yunnan Province,China

Hybridization and introgression are thought to be important for speciation and adaptation in many plants. However, little is known about the hybridization and introgression among Cypripedium species. To investigate the evidence for hybridization and the pattern of introgression between Cypripedium yunnanense and C. tibeticum in Shangrila County, Yunnan Province, China, morphological characters and amplified fragment length polymorphism (AFLP) data for both the species and their putative hybrids were studied. Hand pollination was also performed to verify the crossability of the putative parents. Principal coordinate analysis based on morphological characters and the AFLP data suggested that the putative hybrids were true hybrids of these two Cypripedium species. Analysis with the NewHybrids software indicated that the putative hybrids were F1 generation individuals and backcrosses to C. yunnanese, but no F2 generation was found. Analysis with the Structure software demonstrated asymmetric introgression from C. tibeticum to C. yunnanense. We conclude that natural hybridization and introgression can occur between these two species and that in situ conservation of the parental species is required before fully assessing the evolutionary potential of hybrids.     

Hybridization promotes speciation in Coenonympha butterflies

Hybridization has become a central element in theories of animal evolution during the last decade. New methods in population genomics and statistical model testing now allow the disentangling of the complexity that hybridization brings into key evolutionary processes such as local adaptation, colonization of new environments, species diversification and extinction. We evaluated the consequences of hybridization in a complex of three alpine butterflies in the genus Coenonympha, by combining morphological, genetic and ecological analyses. A series of approximate Bayesian computation procedures based on a large SNP data set strongly suggest that the Darwin's Heath (Coenonympha darwiniana) originated through hybridization between the Pearly Heath (Coenonympha arcania) and the Alpine Heath (Coenonympha gardetta) with different parental contributions. As a result of hybridization, the Darwin's Heath presents an intermediate morphology between the parental species, while its climatic niche seems more similar to the Alpine Heath. Our results also reveal a substantial genetic and morphologic differentiation between the two geographically disjoint Darwin's Heath lineages leading us to propose the splitting of this taxon into two different species.     

Similarity in G matrix structure among natural populations of Arabidopsis lyrata

Understanding the stability of the G matrix in natural populations is fundamental for predicting evolutionary trajectories; yet, the extent of its spatial variation and how this impacts responses to selection remain open questions. With a nested paternal half‐sib crossing design and plants grown in a field experiment, we examined differences in the genetic architecture of flowering time, floral display, and plant size among four Scandinavian populations of Arabidopsis lyrata. Using a multivariate Bayesian framework, we compared the size, shape, and orientation of G matrices and assessed their potential to facilitate or constrain trait evolution. Flowering time, floral display and rosette size varied among populations and significant additive genetic variation within populations indicated potential to evolve in response to selection. Yet, some characters, including flowering start and number of flowers, may not evolve independently because of genetic correlations. Using a multivariate framework, we found few differences in the genetic architecture of traits among populations. G matrices varied mostly in size rather than shape or orientation. Differences in multivariate responses to selection predicted from differences in G were small, suggesting overall matrix similarity and shared constraints to trait evolution among populations.     

Maintenance of strong morphological differentiation despite ongoing natural hybridization between sympatric species of Lomatia (Proteaceae)

Background and Aims

When species cohesion is maintained despite ongoing natural hybridization, many questions are raised about the evolutionary processes operating in the species complex. This study examined the extensive natural hybridization between the Australian native shrubs Lomatia myricoides and L. silaifolia (Proteaceae). These species exhibit striking differences in morphology and ecological preferences, exceeding those found in most studies of hybridization to date.

Methods

Nuclear microsatellite markers (nSSRs), genotyping methods and morphometric analyses were used to uncover patterns of hybridization and the role of gene flow in morphological differentiation between sympatric species.

Key Results

The complexity of hybridization patterns differed markedly between sites, however, signals of introgression were present at all sites. One site provided evidence of a large hybrid swarm and the likely presence of multiple hybrid generations and backcrosses, another site a handful of early generational hybrids and a third site only traces of admixture from a past hybridization event. The presence of cryptic hybrids and a pattern of morphological bimodality amongst hybrids often disguised the extent of underlying genetic admixture.

Conclusions

Distinct parental habitats and phenotypes are expected to form barriers that contribute to the rapid reversion of hybrid populations to their parental character state, due to limited opportunities for hybrid/intermediate advantage. Furthermore, strong genomic filters may facilitate continued gene flow between species without the danger of assimilation. Stochastic fire events facilitate temporal phenological isolation between species and may partly explain the bi-directional and site-specific patterns of hybridization observed. Furthermore, the findings suggest that F₁ hybrids are rare, and backcrosses may occur rapidly following these initial hybridization events.     

Evolution of postzygotic reproductive isolation in galliform birds: analysis of first and second hybrid generations and backcrosses

The process of speciation is a crucial aspect of evolutionary biology. In this study, we analysed the patterns of evolution of postzygotic reproductive isolation in Galliformes using information on hybridization and genetic distance among species. Four main patterns arose: (1) hybrid inviability and sterility in F₁ hybrids increase as species diverge; (2) the presence of geographical overlap does not affect the evolution of postzygotic isolation; (3) the galliforms follow Haldane's rule; (4) hybrid inviability is higher in F₂ than in F₁ hybrids, but does not appear to be increased in the backcrosses. This study contributes to the growing evidence suggesting that the patterns of evolution of postzygotic isolation and the process of speciation are shared among avian groups (and animals in general). In particular, our results support the notion of F₂ hybrid inviability as being key for the maintenance of species genetic integrity when prezygotic isolation barriers are overcome in closely related species, in which postzygotic isolation in the F₁ hybrid might still not be fully developed. To the contrary, hybrids from backcrosses did not show serious inviability problems (at least not more than F₁ hybrids), demonstrating that they could generate gene flow among bird species. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 528–542.     

Male and female secondary sexual traits show different patterns of quantitative genetic and environmental variation in the horned beetle Onthophagus sagittarius

The expression of secondary sexual traits in females has often been attributed to a correlated response to selection on male traits. In rare cases, females have secondary sexual traits that are not homologous structures to secondary sexual traits in males and are thus less likely to have evolved in females because of correlated selection. In this study, we used the dung beetle Onthophagus sagittarius, a species with sex‐specific horns, to examine the environmental and quantitative genetic control of horn expression in males and females. Offspring subjected to different brood mass manipulations (dung addition/removal) were found to differ significantly in body size. Brood mass manipulation also had a significant effect on the length of male horns; however, female horn length was found to be relatively impervious to the treatment, showing stronger patterns of additive genetic variance than males. We found no correlations between horn expression in males and females. We therefore conclude that the horns of O. sagittarius females are unlikely to result from genetic correlations between males and females. Rather, our data suggest that they may be under independent genetic control.     

Molecular and morphological evidence of natural interspecific hybridization between the uncommon Eucalyptus aggregata and the widespread E. rubida and E. viminalis

Human activities can promote increased hybridization in the genus Eucalyptus with potentially detrimental consequences for the persistence of rare species. However, many hybrid combinations have not been investigated with combined use of genetic markers and morphology. We assessed the efficiency of the STRUCTURE program and morphological intermediacy for identifying hybrids between the uncommon tree, Eucalyptus aggregata, which putatively hybridizes with the common congeners, E. rubida and E. viminalis in south-eastern Australia. We sampled 1,005 seedlings across 27 populations, all seedlings were genotyped at 6 allozyme loci and scored for 22 stem and leaf characters. Both marker sets confirmed that E. aggregata is hybridizing with both E. rubida and E. viminalis. Allozymes revealed hybrids from E. aggregata trees in 88% of populations and hybrids comprised 7.3% of all seedlings. Both genetics and morphology indicated that ~50% were likely to be F₁ hybrids, and both simulations and morphological characteristics indicated that the remainder were mostly backcrosses. Morphological analysis correctly distinguished 71% of F₁ hybrids from parentals and was least accurate when dealing with potential backcrosses (50% success). Hence, techniques using genetic data (no prior information) and the assessment of appropriate admixture thresholds through simulations provided the most accurate estimates of hybrid frequency. In this study, potential introgression and the high frequency of hybrids in small populations (~30%), suggests that hybridization should be considered in the management and conservation of E. aggregata.     

MUTATIONAL CONTRIBUTIONS TO GENETIC VARIANCE-COVARIANCE MATRICES: AN EXPERIMENTAL APPROACH USING INDUCED MUTATIONS IN ARABIDOPSIS THALIANA

Genetic potential for evolutionary change and covariational constraints are typically summarized as the genetic variance-covariance matrix G , and there is currently debate over the extent to which G remains effectively constant during the course of adaptive evolution. However, G provides only a temporally restricted view of constraints that ignores possible biases in how new mutations affect multivariate phenotypes. We used chemical mutagenesis to study the effect of mutations as summarized by the mutational covariance matrix, M , in Arabidopsis thaliana. By introducing mutations into three isogenic strains of A. thaliana, we were able to quantify M directly as the genetic variance-covariance matrix of mutagenized lines. Induced mutations generally did not alter the means of the six morphology and life-history traits we measured, but they did affect the levels of available genetic variation and the covariances among traits. However, these effects were not consistent among the three isogenic lines; that is, there were significant differences among the lines in both the number of mutations produced by ethyl-methane-sulfonate treatment and the M matrices they induced. The evolutionary implications of the dependence of M on the number of mutations, the particular genetic background, and the mutagenic sampling of loci in the genome are discussed in light of commonly applied models of multivariate evolution and the potential for the genetic architecture itself to change in ways that facilitate the coordinated evolution of complex phenotypes.