Past climate change on Sky Islands drives novelty in a core developmental gene network and its phenotype
© Favéet al. 2015
Received: 13 July 2015
Accepted: 6 August 2015
Published: 4 September 2015
A fundamental and enduring problem in evolutionary biology is to understand how populations differentiate in the wild, yet little is known about what role organismal development plays in this process. Organismal development integrates environmental inputs with the action of gene regulatory networks to generate the phenotype. Core developmental gene networks have been highly conserved for millions of years across all animals, and therefore, organismal development may bias variation available for selection to work on. Biased variation may facilitate repeatable phenotypic responses when exposed to similar environmental inputs and ecological changes. To gain a more complete understanding of population differentiation in the wild, we integrated evolutionary developmental biology with population genetics, morphology, paleoecology and ecology. This integration was made possible by studying how populations of the ant species Monomorium emersoni respond to climatic and ecological changes across five ‘Sky Islands’ in Arizona, which are mountain ranges separated by vast ‘seas’ of desert. Sky Islands represent a replicated natural experiment allowing us to determine how repeatable is the response of M. emersoni populations to climate and ecological changes at the phenotypic, developmental, and gene network levels.
We show that a core developmental gene network and its phenotype has kept pace with ecological and climate change on each Sky Island over the last ∼90,000 years before present (BP). This response has produced two types of evolutionary change within an ant species: one type is unpredictable and contingent on the pattern of isolation of Sky lsland populations by climate warming, resulting in slight changes in gene expression, organ growth, and morphology. The other type is predictable and deterministic, resulting in the repeated evolution of a novel wingless queen phenotype and its underlying gene network in response to habitat changes induced by climate warming.
Our findings reveal dynamics of developmental gene network evolution in wild populations. This holds important implications: (1) for understanding how phenotypic novelty is generated in the wild; (2) for providing a possible bridge between micro- and macroevolution; and (3) for understanding how development mediates the response of organisms to past, and potentially, future climate change.
How populations differentiate within species is a fundamental problem in evolutionary biology that is key for uncovering the processes that generate biological diversity, including microevolution, speciation, and the emergence of phenotypic novelty [1–6]. Traditionally, biologists have used the tools of population and quantitative genetics, two fields central to the modern evolutionary synthesis, to study population differentiation [1, 7–9]. These fields have significantly advanced our knowledge of variation and change in the frequency of alleles and phenotypes (quantitative and discrete) within and between populations . They have also uncovered signatures of natural selection and genetic drift and are being used to identify loci responsible for adaptive and non-adaptive phenotypes driving the evolution of populations [8, 11–14]. However, a largely unexplored dimension of this problem is what role, if any, does organismal development play in the process of population differentiation in the wild [3, 15, 16].
Organismal development integrates the action of gene regulatory networks with environmental inputs to generate the phenotype [3, 17–19], and therefore, may play a key role in facilitating phenotypic differentiation of populations and species exposed to ecological changes [3, 17, 19]. Major advances in the field of evolutionary developmental biology, such as the discovery of Hox genes, have revealed an unexpected degree of evolutionary conservation of developmental regulatory genes across the animal kingdom [7, 20]. These regulatory genes, which are transcription factors and signalling molecules, are assembled in hierarchically organized networks [18, 20]. A primary function of developmental networks is body plan formation [18, 20]. For example, the expression, structure, function, and regulation of the developmental genes Sonic hedgehog and patched have been conserved for hundreds of millions of years across vertebrates and invertebrates . Both Sonic hedgehog and its downstream target patched regulate anterior/posterior patterning in the developing limb of both fruit flies and chickens . The high degree of conservation and structure of developmental networks suggests that these networks may bias the variation that selection can act upon. Such bias may facilitate repeatable phenotypic responses when populations are exposed to similar environmental inputs or encounter similar ecological changes [1, 3, 22, 23]. Here we use an integrative approach that combines multiple levels of organization (gene networks, development, and phenotype) as well as multiple fields (population genetics, paleoecology and ecology, morphology, and evolutionary developmental biology) to gain a more complete understanding of population differentiation in the wild.
The paleoecological record of the Sky Islands is one of the most complete records available [24, 25]. Numerous studies have reconstructed the major ecological and climatic changes that have occurred in the region by using the stratification of pollen grains and middens of pack rats in the fossil record. These studies show that during the Pleistocene glacial periods ∼90,000–20,000 years BP when atmospheric temperature around the region was substantially cooler, there was a largely continuous forest landscape connecting Sky Island mountain ranges [24, 25]. The region then began warming up between ∼20,000–10,000 years BP, and atmospheric temperatures increased substantially, resulting in the formation of vast deserts that isolated the Sky Islands [24, 25]. Species that occurred at lower elevations in the intervening valleys during the Pleistocene glacial periods, such as Ponderosa pine trees and Douglas fir, are currently only found at high-elevations. These past climatic changes in the region have generally caused shifts in elevation, expansions, and contractions of populations in numerous species [28–32]. Each Sky Island is now considered to be one of the steepest terrestrial ecological gradients in North America, from low-elevation deserts, through mid-elevation oak-juniper woodlands, to high elevation coniferous forests . Several ecological factors vary steeply along these altitudinal clines, in particular temperature, precipitation and habitat fragmentation [33–35]. Environmental changes that would normally be observed across a latitudinal gradient stretching thousands of kilometres can be observed across a few hundred meters of elevational change in the Sky Islands [26, 27, 36]. This poses a unique adaptive challenge for organisms on Sky Islands that have distributional ranges that encompass more than one altitudinal zone [26, 27, 36]. These steep gradients therefore behave as powerful sensors of climatic and environmental changes because noticeable changes in communities occur rapidly and over short geographical distances [33–35].
Ants are an ecologically dominant group  that have long been used as bioindicators of environmental change [38, 39] and climate warming [40, 41]. In most ant species, the colony is made up of a single winged queen that performs most of the reproduction in the colony and her wingless workers that perform almost all the other tasks . Throughout most of the year, the queen produces eggs that develop into wingless workers, but in response to specific environmental conditions associated with spring time, such as the gradual increase in daylight and temperature, they begin producing virigin queens and males that use their fully functional wings to disperse and participate in mating flights away from the mother colony. Immediately after these mating flights, males die while queens tear off their wings to found a new nest underground . This life history strategy is known as ‘independent colony foundation’ [42–45]. In some ant species, however, an alternative life history strategy called ‘dependent colony foundation’ has evolved [42–45]. In these species, there are multiple queens in a single ant colony in which some or all of the queens are wingless. Virgin queens that are wingless mate on the ground and disperse on foot aided and accompanied by wingless workers from the mother colony. This strategy dramatically decreases mortality of wingless queens, but at the expense of long distance dispersal. Wingless queens have evolved multiple times independently within ants [43–45], suggesting that wingless queens represent an adaptive life-history strategy. It takes less energy to produce wingless queens, and selection may favor the evolution of wingless queens in particular ecological conditions where aerial dispersal is risky [29, 46, 47].
Colonies of the ant species M. emersoni are distributed throughout the Sky Islands, where we discovered the presence of two alternative queen phenotypes, winged and wingless (Fig. 1 c and d) . A single M. emersoni colony can be headed by one or multiple queens, ranging from 1 to more than 70, and colonies can be composed either exclusively by winged queens, exclusively by wingless queens, or a mixture of winged and wingless queens . Current evidence suggests that determination of alternative queen phenotypes in a colony occurs during development and can be influenced by both environmental and genetic factors [43, 49, 50]. On each Sky Island, M. emersoni occurs sparsely at lower elevations, where its abundance is very low and is associated with creek or woody patches near mountain ranges . At mid-elevations of each Sky Island, from ∼1500–2300 meters, M. emersoni occurs at high densities in oak-juniper woodlands, whereas at the highest elevations of each Sky Island, from 2300 m and above, M. emersoni occurs in patches of high density in coniferous forests dominated by Pinyon pine and Douglas fir .
Studying how population differentiation occurs in M. emersoni along several naturally replicated ecological gradients on five Sky Islands in Arizona provides a unique opportunity to ask whether organismal development facilitates repeatable, and thus, predictable responses to similar ecological and climate changes on each Sky Island. Furthermore, the interaction between ecology and organismal development in the context of Sky Islands is key for understanding population dynamics of gene networks and its impacts on phenotypic evolution on microevolutionary timescales, as well as for understanding how past changes in developmental gene networks may affect the future course of evolution in response to climate change.
A combination of landscape fragmentation and elevation best explain the distribution of wingless queens along ecological gradients. Based on previous studies, we generated 12 models (M01 to M12) representing alternative hypotheses to explain the variation in frequency of wingless queens along the ecological gradient of each Sky Islands. The 12 models were ranked using the method of  implemented in R (Additional file 1: Table S2) and model probabilities (Prob) are reported, with m06 scoring as the highest ranking model
Underlying biological hypotheses
Elevation influences a wide range of ecological factors and may account for their joint effects.
Atmospheric pressure is lower at higher altitude and efficient flight may require costly adaptations .
Insects will experience a decrease in annual thermal budget for growth and development. Critical temperature thresholds for growth, development and activity will be exceeded less frequently .
Temperature seasonality, the variability in temperature over the year, is indicative a less stable environment. Flightless insects are known to be found in higher frequencies in more stable habitats .
Flying insects will have a tendency to be blown away from patchy habitats, thus only wingless forms remain [46, 47]. The probability of finding a suitable habitat is decreased when habitats are patchy, thus dispersal risk is increased .
The need for dispersal in high productivity habitats is reduced because the conditions for survival are met and the relative risk of long distance dispersal becomes too high 
Frag + Elev
Both fragmentation and elevation influence the distribution of wingless queens.
Frag + Elev + Frag x Elev
The effect of fragmentation depends on elevation: fragmentation may only be affecting the distribution of wingless queens at certain altitudes.
Frag + Temp
Both fragmentation and temperature may influence the distribution of wingless queens.
Frag + Temp + Frag x Temp
The effect of fragmentation depends on temperature: fragmentation may only be affecting the distribution of wingless queens at certaine temperature. Habitat fragmentation and temperature are know to be potentially interrelated .
Frag + Prod
Both fragmentation and productivity can influence the distribution of wingless queens.
Frag + Prod + Frag x Prod
Habitat fragmentation can affect productivity. Fragmentation may affect wingless queen distribution in habitats with a certain productivity level only. Habitat fragmentation and biomass growth are know to be potentially interrelated .
Elev + Temp + Prod + Frag + Seas
Structure-rich model, including the effects of elevation, temperature, productivity, fragmentation and seasonality. Interactions were omitted.
Our data show that emergence of the wingless queen phenotype occurred repeatedly by parallel evolution in response to similar ecological conditions arising on each of the Sky Islands during climatic warming. First, modern gene flow between mountain ranges is absent or low; the degree of genetic differentiation (F ST ) is high between collection sites on different Sky Islands, but is lower within a Sky Island (Additional file 1: Table S3 and S5). Second, phylogenetic and redundancy analyses based respectively on multilocus AFLP data and mtDNA sequences show that M. emersoni colonized each Sky Island independently (Fig. 3 b, Additional file 1: Figure S3 and S4). We find that populations on each Sky Island contain unique haplotype variants, while no shared haplotypes were recovered among Sky Islands, indicating that unique haplotypes originated within each Sky Island (Additional file 1: Figure S3). Finally, our results show that geography dictates the patterns of genetic differentiation at neutral loci: winged and wingless queens from the same Sky Island share a more recent ancestor than queens from different Sky Islands (Additional file 1: Figure S3 and S4), and genetic variation among queens is best explained by geography (collection site) rather than habitat or phenotype (Fig. 3 b). All together, this supports an incipient process of parallel differentiation of the populations and suggests that M. emersoni wingless queens emerged independently on each Sky Island when each population encountered similar but newly established ecological conditions, namely, fragmented and high elevation habitats.
The presence of wing loss in queens in Monomorium minimum-group species. Summary of known queen phenotypes for Monomorium minimum-group species that are native to the United States and Mexico; species are listed alphabetically
Winged and wingless c,d
Winged and wingless a,b
Winged and wingless a
Winged and wingless a,c
M. sp. AZ-01
M. sp. AZ-03
Winged and wingless d
M. sp. cf. ergatogyna- Baja
We then determined whether the expression of these genes has been interrupted in wingless queens, and found that their expression is altered relative to winged queens (Figs. 5, 6 and 7 and Additional file 1: Figure S6). In wingless queens from different Sky Islands, we found that the 4 genes involved in wing tissue patterning (ubx, en, cut, wg) show similar aberrant patterns of expression (Fig. 5 e to t and Additional file 1: Figure S6C to H). These shared alterations reveal that wingless queens, regardless of their geographic origin and extent of vestigial wing disc growth (Fig. 4 a), interrupt wing patterning at least in part through similar changes in gene expression. Therefore, although organ growth retains the signature of historical events caused by past demographic changes, alteration of the underlying expression of wing patterning genes is recurrent among populations on each Sky Island.
In contrast to genes involved in wing patterning, we observed that the expression of mef2 - a key transcription factor involved in wing muscle development  - differs among populations and shows a signature of the demographic split between the southern and the northern Sky Islands. Mef2, which is expressed in the nuclei of myoblasts and appears as a dotted pattern in wing discs, shows expression differences between wingless queens from northern and southern Sky Islands (Fig. 6 f and j versus N and R, and Additional file 1: Figure S7). In the northern Sky Islands, Mef2 expression is restricted to the outer edge of the vestigial disc in wingless queens, whereas in the southern Sky Islands, it is expressed more extensively and covers the central region of vestigial discs. Thus, like patterns of wing disc growth (Fig. 4 h), Mef2 expression patterns precisely recapitulate the history of contacts between Sky Islands. Furthermore, these differences in the pattern of Mef2 expression in vestigial wing discs are likely to have evolved neutrally because wingless queens lack flight muscles regardless of their origin (Fig. 6 g, h, k, l, o, p, s and t). This suggests that another mechanism, such as the interruption of other genes in the network or apoptosis of myoblasts , prevents mef2 expressing cells from further differentiating into flight muscles.
A paleoecological reconstruction for the parallel evolution of a developmental network and its phenotype in response to past climate change on Sky Islands
We have provided evidence for the following scenario, which shows that in response to both direct and indirect effects of past climate change on Sky Islands, a developmental gene network and its corresponding phenotype can respond to changing ecological conditions on short time scales (Fig. 8): Paleoecological reconstructions of the Arizona Sky Islands [24, 25] show that during the last glaciation between about 90,000 to 20,000 years BP, a continuous forest habitat across the landscape prevailed, except for high elevations which were inhospitable alpine habitats. The developmental potential to produce wingless queens in M. emersoni was present during this period. During or immediately following the vicariant event ∼80,000 years BP causing the demographic split between northern and southern Sky Islands, these M. emersoni populations accumulated unexpressed differences in the genes that control wing growth and flight muscle development. Subsequently, during climatic warming ∼10,000 years BP, the landscape became fragmented and new habitats arose at high elevations that were suitable for colonization by M. emersoni. As populations encountered these new habitats on each of the Sky Islands, the wingless phenotype evolved independently and became more frequent through selection on the phenotypic expression of the developmental potential for producing wingless queens. The parallel evolution of wingless queens on each Sky Island occurred through recurrent interruptions of similar wing-patterning genes in the network (ubx, en, cut, wg), and at the same time, induced the release of the unexpressed differences in vestigial wing growth and mef-2 expression that accumulated in populations during or after the demographic split between the northern and southern Sky Islands ∼80,000 years BP. Climate warming during the last ∼10,000 years BP formed vast deserts between Sky Islands thereby isolating the five Sky Island populations from one another and forced M. emersoni populations to adapt physiologically to the warm conditions at the lower elevations of their range. This may have led to the rapid evolution of the vestigial wing hinge and its underlying exd expression on each Sky Island.
Novelty as a mosaic of contingency and determinism
We show that developmental systems can facilitate two types of evolutionary change that contribute to the differentiation of populations within a species: one that is governed by chance events, is unpredictable, and has produced small differences in phenotype and genotype: wing hinge morphology, wing disc growth and expression of exd and mef-2. In contrast, the other type is repeatable and has produced large differences in phenotype and genotype: the wingless queen. The evolution of wingless queens represents an evolutionary novelty - it is the gain of a developmental switch, the complete loss of wings and wing muscles, multiple alterations in the expression of the core developmental network underlying wing development, and the gain of a new life history strategy (dependent colony foundation). These findings offer a new perspective for the ongoing discussion of novelty generation in phenotypic evolution . On the one hand, it has been argued that adaptive solutions through natural selection are deterministic and that evolution is largely predictable , while on the other, it has been argued that historical contingency, which is the unique chain of prior events that determine which of several evolution paths a lineage will follow, renders evolution unpredictable . Several studies, from a phylogenetically diverse array of organisms, show that the parallel evolution of phenotypic traits often occurs through recurrent alterations in the same genes [12, 23, 71–74]. Yet, at the same time, these and other studies report slight differences across multiple levels of biological organization: phenotype, development, and genetics [75–77]. In vitro studies on microbial populations suggest that these similarities and differences underlying the parallel evolution of traits are the consequence of a complex interplay between deterministic processes and historical events [78–80]. We show, in natural populations within a species, that the emergence of novelty at multiple biological levels can be generated by the combination of ancestral developmental potentials facilitating parallel evolution in response to similar ecological conditions and chance events that create slight variations on these repeatable evolutionary novelties.
A mosaic of contingent and deterministic processes at multiple biological levels may be a general property of the evolution of novelty. Most, if not all, organisms experience similar population dynamics over evolutionary time and parallel evolution within and between species is a nearly universal pattern throughout the tree of life [69, 81]. Many independently evolved traits, such as, pigmentation spots on Drosophila wings [53, 82], pharyngeal jaws of cichlid fish , hind limbs of anoles lizards , coat color of wild mice , and coloration on butterfly wings , exhibit patterns consistent with our findings in that they show slight variations in phenotype and genotype. Our results suggest that these slight variations are generated by microevolutionary forces and historical contingency acting between populations, together causing the gradual accumulation and evolution of unexpressed changes in regulatory pathways. These unexpressed differences become expressed during the process of parallel evolution leading to novel variations on independently evolved traits. This process driving parallel phenotypic novelty within species may translate into parallel phenotypic novelty between species, and as in the case pharyngeal jaws of cichlid fish , slight variations accumulated in one of the lineages may potentially fuel phenotypic diversification.
Implications for bridging Micro- and Macroevolution
One of the biggest challenges facing evolutionary theory today is to explain how microevolution translates into macroevolution. Evolutionary biologists have addressed this problem from the perspective of paleontology [6, 87], population and quantitative genetics , and evolutionary developmental biology [16, 89, 90], producing a range of different perspectives on the translation between micro- and macroevolution. One perspective, which took hold during the neo-Darwinian synthesis in the 1940’s, proposes no distinction micro- and macroevolution, where macroevolution is simply the product of a specific kind of microevolution (mutations of small effect, individual level selection and drift) extrapolated over long temporal scales . This perspective, however, is not well supported in the face of recent advances in quantitative genetics [8, 91], experimental evolution , developmental plasticity  and paleontology . These recent advances show that: (1) dynamics of evolutionary change observed between individuals within populations cannot be extrapolated to explain the dynamics of change observed between species within clades [87, 93]; and (2) alleles of large effect [8, 91] and large discrete phenotypic variants [3, 54, 94] can be fixed in populations. Indeed, most modern textbooks of evolution now distinguish micro- from macroevolution and use the species boundary as the dividing line between them, most often defining microevolution as ‘evolution below the species level’ - small changes, adaptive or neutral, that arise within populations, while macroevolution is most often defined as ‘evolution above the species level’ - large changes that arise within populations leading to the origin of species or higher taxonomic groups [2, 95–98]. Although all changes must necessarily arise within populations, macroevolutionary changes are thought to arise during or after populations have undergone speciation, and therefore, characterize differences observed between closely-related species or higher taxonomic groups [3, 93].
Our findings, however, suggest that the species boundary may not adequately separate micro- and macroevolutionary changes. We observed macroevolution-like change (emergence of wingless queens) between populations within the same species (M. emersoni), which is comparable to the variation observed between species in the genus Monomorium (Table 2). Our observations are consistent with an idea known as ‘intraspecific macroevolution’, which proposes that micro- and macroevolution are distinct, but that speciation is not the primary cause of large phenotypic changes between species leading to macroevolution . It proposes that both small and large phenotypic changes occur within the same species, and that the existence of alternative phenotypes within lineages facilitates macroevolution by reducing the negative fitness consequences associated with the appearance and fixation of large phenotypic variants .
Finally, it has been proposed that changes in hierarchically organized developmental networks can simultaneously produce the large- and small-scale changes associated with macro- and microevolution [99, 100]. Genes in the network involved in core developmental processes (‘kernels’ sensu Davidson, ), such as body plan formation, are thought to produce large-scale variants when changed, be highly stable during evolution, and are associated with macroevolution. In contrast, genes in the network involved in terminal differentiation processes, like pigment formation, are expected to produce small-scale variants, be more labile during evolution, and are associated with microevolution [99, 100]. Our results are consistent but show that the core developmental network responsible for wing formation facilitates the evolution of both types of changes within a species in response to ecological change. Based on all of our findings, we propose that the parallel phenotypic novelty within species may translate into parallel phenotypic novelty between species providing a possible bridge between micro- and macroevolution.
The role of organismal development in past and future climate change
Understanding how developmental systems respond to past climate change is critical to facing the challenges of global warming over the next few decades. Integrating the responses of developmental systems with environmental changes will be necessary for building accurate projection models to predict species distributions under climatic change scenarios. Predicting species distributions accurately can affect decision making about conservation biology strategies, agriculture and public health issues. Current projection models integrate demographic processes such as migration and population size changes, but few integrate the possible response of developmental systems to environmental changes, which as we show, may profoundly affect the life history and dispersion abilities of a species. Ultimately, such changes directly or indirectly affect species interactions and ecosystem processes . Sky Island mountain forests have been used as models and indicators of climate change. They have responded to past  and present  global warming by contracting towards the higher elevations , and climatic and ecological models predict that warming and desertification in Arizona will continue over the next 50 years [33, 34, 102]. Several species found in the mountainous areas of the American Southwest already show signs of changes in phenology in response to recent global warming [103, 104]. These dramatic changes in environmental conditions will be integrated by the developmental systems of different organisms, and as we show, can generate both predictable and unpredictable phenotypic responses. Generally integrating the predictable part of the response of developmental systems into climate change studies should hold the promise of making more accurate predictions for organismal responses to climate change.
Since the modern synthesis, our view of how population differentiation occurs in the wild is largely based on the advances made in population and quantitative genetics, especially with the advent and integration of genomics tools. This view is one of changes in allele and phenotypic frequencies in response to adaptive and non-adaptive forces. By integrating organismal development and paleoecology on Sky Islands, we enrich this view by showing that the emergence of novelty within populations occurs through a combination contingent and deterministic processes. This combination results in the parallel evolution of traits with slight variations at multiple levels of biological organization (phenotype, growth, developmental networks). We propose that this parallel phenotypic novelty within the species may represent the springboards for translation into macroevolutionary changes between species. Studying these evolutionary processes on the Sky Islands, with their rich fauna provides replicated ecological gradients, are a good model system for understanding how development influences an organism’s response to past and future climate change.
Materials and methods
Animal collection and culturing
We discovered the M. emersoni wingless phenotype in the Chiricahuas mountain range in 1999. We found that this wingless phenotype exists on four other Sky Island mountains in Southeastern Arizona, USA and we collected colonies from 14 sites (Additional file 1: Table S3). We preserved colonies in 95 % ethanol or kept them alive in plastic boxes with glass test tubes filled with water constrained by a cotton plug. We fed colonies with a mix of crickets, mealworms and Bhatkar-Whitcomb diet . The colonies were maintained at 27 ° C and 70 % humidity with a day/night cycle of 12 hours. We induced production of winged and wingless queens by exposing them to a cycle of temperature range: the temperature of the growth chamber was reduced by 1 ° C for 17 days until is reached 10 ° C. We kept the chamber at 10 ° C for 2 weeks, and then increased the temperature gradually by 1 ° C per day for 17 days until the original conditions were restored.
We performed all statistical analyses using R, unless otherwise stated.
Inferring modern and historical population demography
Sequencing and genotyping
We extracted DNA from 318 M.emersoni queens, each from a different colony, using thoracic tissue from preserved individuals in ethanol. We amplified by PCR two mitochondrial genes (cytochrome oxidase 1 (CO1) and 2(CO2)), and sequenced them in both directions on an ABI3730 (for primers, see ) for a total of 1350 base pairs (bp) for 70 individuals. We also amplified and sequenced 1020 bp from a nuclear gene, hymenoptaecin, including one intron. For each Sky Island, we isolated alleles from hymenoptaecin to create a pool of DNA for each Sky Island by mixing DNA from about 20 queens from different colonies. This method has been shown to be accurate to estimate allele frequencies of a population [107, 108]. We amplified all alleles in the DNA pool simultaneously by PCR and cloned them into a PGEM-T vector. We picked 15 bacterial clones per mountain and we sequenced them with M13 primers in both directions. This allowed us to obtain the exact haplotype phase for every nuclear allele sampled and sequenced. All sequences were assembled and aligned with Geneious Pro V.5.4.3 and checked by eye. GenBank accession numbers are KT356282-KT356545.
We also genotyped all individuals for Amplified Fragment Length Polymorphisms (AFLP) markers, as described in . We used 4 primer combinations to generate 157 unambiguous loci: M CTA and E ACC , M CTA and E ACG , M CAG and E ACC , and M CAG and E ACG . We ran the PCR products on a LiCor DNA analyser 4300 and the same researcher scored them blindly using the accompanying Saga software. To estimate the error rate in the procedure, about 10 % of the samples were amplified and scored twice, and we recovered a similar error rate to what has been reported in other studies .
Estimation of modern gene flow
We estimated F ST values with the Bayesian method developed by  and implemented in AFLP-surv . We removed from the analysis sites where we collected less than 15 individual queens from different colonies. P-values for the significance of each F ST were estimated after applying Boneferroni corrections for multiple comparisons. Pairwise F ST estimates for all 14 sites are reported in Additional file 1: Table S5. Stars represent parameters that could not be estimated with confidence. To test for the presence of recent gene flow, we performed a Mantel test  to calculate the correlation between the three-dimensional geographic distance and the genetic differentiation (F ST estimates) for all pairwise comparisons among sites. We expected a significant correlation to occur if restricted but non-zero gene flow existed among Sky Islands, or if Sky Islands diverged recently from a continuous ancestral population without recent gene flow. On the other hand, we expected a non-significant correlation for the two following demographic scenarios: a complete panmixia among most sites, or a highly reduced gene flow for a long time among most sites. We tested the Mantel coefficient for significance using 10,000 random permutations on the matrices.
Estimation of historical gene flow and times of divergence
We used the Bayesian search strategy implemented in MIGRATE-n 3.2.15  to estimate the effective population sizes, the average historical migration rates, the historical migration rates and the coalescence times for each population. We also evaluated 2 models of population subdivision by comparing their marginal likelihoods using the modified thermodynamic integration (TH) implemented in MIGRATE-n . The first model considered all sampling sites as part of one panmictic unit while the second model considered each mountain range as independent units (limited or absence of gene flow). Each MIGRATE-n run consisted of 10 replicates of 4 MCMC heated short chains (heating scheme: 1.0, 1.2, 3.0, 6.0) and one long chain of 5,000,000 sampled trees where 500,000 burn-in trees were discarded to ensure proper sampling of the parameter space. All parameter estimates showed narrow posterior probability distributions, indicating high confidence in the estimates. All comparisons involving the Chiricahuas as the sending population show clear outlier values and were thus not considered. Θ values from the haploid and maternally-transmitted mitochondrial loci were corrected for comparison with a diploid locus to obtain an estimate of 4Ne μ. We used the direct estimate of mtDNA mutation rate calculated by  and the immune gene mutation rate estimate in Drosophila from  to calculate migration rates and times of divergence. MIGRATE-n also implements the Skyline plot method for detecting changes in migration rates over time, thus allowing the reconstruction of historical migrant exchanges between mountain ranges. The analysis was run separately for each pair of geographically adjacent mountain ranges (Chiricahuas-Huachucas, Huachucas-Catalinas, Catalinas-Pinals, Pinals-Pinaleños and Pinaleños-Chiricahuas) with the same run parameters as described above.
Testing for competing evolutionary scenarios of population splits
We used Approximate Bayesian Computation implemented in the program DIYABC v.126.96.36.199 beta  to test competing hypotheses of lineage divergence. Instead of estimating the likelihood of a model from an MCMC approach, DIYABC uses approximate Bayesian computation , where similarity of summary statistics between observed and simulated data sets are compared. Data can be simulated for several scenarios, where population sizes, times of divergence, population size change and admixture can be specifically defined.
In order to focus on testing the most probable scenarios, we built up our hypotheses based on the results of the Skyline plot analysis from MIGRATE-n (Additional file 1: Table S4 and Additional file 1: Figure S1). We compared three divergence scenarios; (A) a simultaneous isolation of the five mountain ranges from a common ancestor, (B) a sequential divergence of Sky Islands based on the results of the Skyline plot analysis using as priors the estimated times where no more gene flow is observed between two Sky Islands, and (C) a third scenario based on scenario B, but with the Huachucas and the Catalinas interchanged on the topology, thereby affecting the relationships within the Southern and Northern Sky Islands (Additional file 1: Figure S2).
The prior probability uniform distribution on the oldest divergence event (t4) ranged from 10,000 to 200,000 generations in the past and the prior probability uniform distribution on the effective population (Ne) size ranged from 1000 to 200,000 (coherent with MIGRATE-n estimates of Ne). The sequence of subsequent splits was forced to be kept in this specific order when simulating the data and the time interval prior probability distribution for each of these respective splits were t4 =1 to 50,000, t3 =10,000 to 70,000, t2 =25,000 to 100,000 and t1 =10,000 to 200,000 generations in the past. The prior on Ne was kept constant through time and was sampled from a uniform distribution ranging from 1000–100,000. Most importantly, three conditions were imposed for the sampling of parameters for the population splitting to occur in the specified order: t1 ≥ t2, t2 ≥ t3 and t3 ≥ t4. 300,000 simulations were performed and the direct estimate method was used to estimate the posterior probabilities of each scenario.
Parallel evolution of M.emersoni populations
Inferring mitochondrial haplotype network
To get insight about whether the wingless phenotype evolved repeatedly in parallel or has a single origin, we first infered genealogical realtionships among queen haplotypes across the Sky Islands populations. A mitochondrial DNA network was constructed using the R package ape. Divergent haplotypes (12 and 20 steps away from the central haplotype were discovered in the Chiricahuas but left out of the analysis for simplicity of the representation of the network.
A neighbor-joining tree was constructed using AFLP data with the phylogenetic package Phylip . The tree was constructed using sequentially the executables Genedist, Neighbor, Seqboot, Consense and Drawtree. Bootstrap values over 70 are reported on the tree.
We first calculated the pairwise mtDNA genetic distance between each pairs of queens and performed a multidimensional scaling to obtain a set of synthetic variables that best represented the pairwise distances between records. To quantify habitat type in a discrete manner, we performed a principal component analysis using all environmental variables on all collection sites, assigning a habitat type (1 to 5) to sites with high similarity. We built three explanatory matrices: collection site, phenotype (winged or wingless) and environment (habitat type). We partitioned the variation in genetic distance with respect to the three explanatory matrices using a redundancy analysis ordination.
Adaptation to climatic changes
Tests for temperature tolerance
Heat resistance of worker ants was measured using a knockdown assay. About 20 individuals workers were collected from inside the nest of each colony, placed into 8 mL glass vials, and incubated at 50 ° C for one hour into a hybridization oven. Resistance was scored for each ant as the time taken to be knock down and mortality after treatment was recorded for each sample. Cold resistance was evaluated as recovery from a chill coma. About 20 individuals workers were collected from inside the nest of each colony, placed into small petri dishes, and immersed for 16 hours into a box filled with ice and water. This treatment led to rapid immobilization of all individuals. The temperature was kept constant around 4 ° C. Petri dishes were then exposed to ambient temperature and recovery was scored for each worker as the time taken for muscular coordination to come back.
Environmental factors the winged and wingless phenotypes may be adapting to along the ecological gradients
Explanatory variables for queen phenotype distribution
We calculated the percentage of wingless queens present at each site. This value represents the proportion of wingless queens that were successful in founding or persisting in a colony, and thus, have survived the dispersal/founding selection period. For each site, we measured the latitude and elevation in the field, and we obtained climatic variables from WorldClim Atlas and the vegetational growth average and vegetational growth peak were obtained from the National Atlas (Additional file 1: Table S1). The Worldclim Atlas and the National Atlas data have a 1 km2 resolution, which is, from our observations, of the same magnitude as the area of our average sampling sites.
We performed an information-theoretic approach to evaluate the performance of 12 models explaining the frequency distribution of winged and wingless queens along the ecological gradients. This approach, widely used in ecology, does not reject or accept models, but rather ranks a priori models and provides posterior probabilities that allow for the interpretation of the performance of each model. Using the findings of previous studies in the literature, we constructed 10 biologically plausible models that could explain which ecological factor(s) have the most influence on distribution of wingless queens along the ecological gradients. We used the method and R code of  to fit and rank models: we fitted models using the R function lm, and we calculated the sample-size corrected Akaike Information Criterion (AICc) for each model. We controlled for population differentiation using the score of the first component from a principal component analysis on the mtDNA pairwise genetic distances among all individual queens. We ranked the models based on AICc and calculated their posterior and conditional probabilities for further interpretation.
Interruptions of the gene network controlling wing development
We collected queen larvae at late last larval instar, just prior to the prepupal stage, following the criteria of . We fixed larvae for 2h in a 4 % formaldehyde solution in PEM buffer , and dissected them under a Zeiss Discovery V12 Stereomicroscope to expose their thoracic imaginal discs by removing obstructive tissues.
Wing imaginal disc morphology and cell division
We calculated surface area of the forewing imaginal disc and of a leg imaginal disc (in μm2) using AxioVision software (Carl Zeiss Canada Ltd., Toronto, Ontario, Canada). For all samples that were undamaged we counted the number of Phosphohistone expressing cells in the forewing imaginal disc from a z-stack image. To test for a difference in static allometry of leg and wing disc among populations, we performed a Standardized Major Axis bivariate line-fitting analysis (SAM) and tested for slope differences  (Additional file 1: Figure S6A). We also performed an analysis of variance on the number of cells undergoing mitosis per wing disc unit area for each population (Fig. 4).
Expression patterns of genes of the wing patterningnetwork
Immunohistochemistry: we performed antibody stainings for Engrailed (En), Ultrabithorax (Ubx), Cut (Cut) and Extradenticle (Exd) proteins following the protocol of . Antibody stainings for Myosin enhancing factor-2 (Mef2) and Phosphohistone-3 (PH3) were conducted similarly, except that the tissues were treated with PAT 1 % instead of PTW 1 % to increase the antibody penetration in the tissue. We systematically included winged queen larvae in each staining reaction containing wingless queen larvae as positive controls and we over developed the staining reaction for larvae of both winged and wingless queens to confirm that the absence of expression was real and not an artifact. We imaged the stained tissues with a Zeiss microscope using the AxioVision software. Expression of genes in the network that controls wing development is highly conserved across winged castes of ants [63, 65]. Nonetheless, we confirm that winged queens from several collection sites do not show any expression differences (data not shown).
Whole mount in-situ hybridization: We isolated a fragment of the gene wingless by PCR using the primers of  on cDNA synthesized from embryonic and larval RNA. We cloned and sequenced this fragment to confirm its identity. Genbank accession number is KT361189. We followed the whole mount in-situ protocol from .
3D X-ray microtomography
In order to describe in more detail the differences in external and internal anatomy of the thorax, we made 3D X-ray microtomography images (microCT scans) of winged queens and 11 wingless queens of different geographic origins. Ants were stained with elemental iodine (1 % w/v in absolute ethanol) to enhance the X-ray contrast of the soft tissues . The scans were performed on an Xradia MicroCT 3D imaging system with voxel sizes between 1.8 and 4.2 μm. We examined sagittal, coronal, and transverse sections for the presence of direct and indirect flight muscles. We applied a surface rendering algorithm using OsiriX to visualize the 3D reconstruction of the external anatomy of the thorax.
The population sets of DNA sequences generated in this study are available on GenBank (accession numbers KT356282-KT356545). Other data used in this study are available upon request to the authors.
We thank J. Losos, G. Wray, R. Barrett, T. Sanger, Y. Idaghdour, M.B. Dijkstra, A.P. Hendry, L. Nilson, J-P. Lessard, A. Gonzalez, J. Marcus, K Chaudhary, and Abouheif Lab members for comments. We also thank E. Furlong for the anti-mef2 antibody and R. Keller for technical advice. This work was supported by grants from Fonds de la Recherche du Québec - Nature et Technologies, Canada Research Chairs program, and National Sciences and Engineering Research Council (NSERC) to E.A., NSERC postgraduate scholarship to M-J.F.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Schluter D, Nagel L. Parallel speciation by natural selection. American Naturalist. 1995; 146:292–301.View ArticleGoogle Scholar
- Futuyma DJ. Evolutionary Biology: Sinauer; 1998.Google Scholar
- Eberhard MJW. Developmental Plasticity and Evolution. New-York: Oxford University Press; 2003.Google Scholar
- Darwin C. On the Origins of Species by Means of Natural Selection. London: Murray; 1859.Google Scholar
- Mayr E. Systematics and the Origin of Species, from the Viewpoint of a Zoologist: Harvard University Press; 1942.Google Scholar
- Gould SJ, Eldredge N. Punctuated equilibria: an alternative to phyletic gradualism In: Schopf TJM, editor. Models in Paleobiology. San Francisco: Freeman Cooper: 1972. p. 82–115.Google Scholar
- Raff RA. The Shape of Life: Genes, Development, and the Evolution of Animal Form. Chicago: University of Chicago Press; 1996.Google Scholar
- Orr HA. The population genetics of adaptation: the distribution of factors fixed during adaptive evolution. Evolution. 1998;:935–49.Google Scholar
- Hartl DL, Clark AG, Clark AG. Principles of Population Genetics, vol. 116. Sunderland: Sinauer associates; 1997.Google Scholar
- Nielsen R. Molecular signatures of natural selection. Annu. Rev. Genet. 2005; 39:197–218.PubMedView ArticleGoogle Scholar
- Jones FC, Chan YF, Schmutz J, Grimwood J, Brady SD, Southwick AM, et al. A genome-wide snp genotyping array reveals patterns of global and repeated species-pair divergence in sticklebacks. Current Biology. 2011; 22(1).Google Scholar
- Chan Y, Marks M, Jones F, Villarreal G, Shapiro M, Brady S, et al. Adaptive evolution of pelvic reduction in sticklebacks by recurrent deletion of a pitx1 enhancer. Science. 2010; 327.5963:302–5.View ArticleGoogle Scholar
- Arnegard ME, McGee MD, Matthews B, Marchinko KB, Conte GL, Kabir S, et al. Genetics of ecological divergence during speciation. Nature. 2014.Google Scholar
- Soria-Carrasco V, Gompert Z, Comeault AA, Farkas TE, Parchman TL, Johnston JS, et al. Stick insect genomes reveal natural selection’s role in parallel speciation. Science. 2014; 344(6185):738–42.PubMedView ArticleGoogle Scholar
- Pfennig DW, Wund MA, Snell-Rood EC, Cruickshank T, Schlichting CD, Moczek AP. Phenotypic plasticity’s impacts on diversification and speciation. Trends in Ecology & Evolution. 2010; 25(8):459–67.View ArticleGoogle Scholar
- Nunes MD, Arif S, Schlötterer C, McGregor AP. A perspective on micro-evo-devo: progress and potential. Genetics. 2013; 195(3):625–34.PubMed CentralPubMedView ArticleGoogle Scholar
- Gilbert SF. Developmental Biology. Sunderland: Sinauer Associates; 2010.Google Scholar
- Davidson EH. The Regulatory Genome: Gene Regulatory Networks in Development and Evolution. Burlington: Academic Press; 2010.Google Scholar
- Gilbert SF, Epel D. Ecological Developmental Biology. Sunderland: Sinauer Associates; 2009.Google Scholar
- Carroll SB, Grenier JK, Weatherbee SD. From DNA to Diversity: Molecular Genetics and the Evolution of Animal Design. Cambridge: Blackwell; 2001.Google Scholar
- Marigo V, Scott MP, Johnson RL, Goodrich LV, Tabin CJ. Conservation in hedgehog signaling: induction of a chicken patched homolog by sonic hedgehog in the developing limb. Development. 1996; 122(4):1225–33.PubMedGoogle Scholar
- Stern DL. Evolution, Development, & the Predictable Genome: Roberts and Company Publishers; 2011.Google Scholar
- Stern D, Orgogozo V. Is genetic evolution predictable?Science. 2009; 323(5915):746–51.PubMed CentralPubMedView ArticleGoogle Scholar
- Lomolino M, Brown J, Davis R. Islands biogeography of montane forest mammals in the american southwest. Ecology. 1989;:180–94.Google Scholar
- Thompson RS, Anderson KH. Biomes of western North America at 18,000, 6000 and 0 C-14 yr bp reconstructed from pollen and packrat midden data. Journal of Biogeography. 2000; 27(3):555–84.View ArticleGoogle Scholar
- Warshall P. The madrean sky island archipelago: a planetary overview In: DeBano L, Ffolliott P, Ortega-Rubio A, Gottfried G, Hamre R, Edminster C, editors. RM-GTR-264 : Biodiversity and Management of the Madrean Archipelago: the Sky Islands of Southwestern United States and Northwestern Mexico. Fort Collins: US Department of Agriculture, Forest Service: 1994. p. 6–18.Google Scholar
- Gottfried GJ, Gebow BS, Eskew LG, Edminster CB. Connecting Mountain Islands and Desert Seas: Biodiversity and Management of the Madrean Archipelago II. Proc. RMRS-P-36. Fort Collins, CO: U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station: 1–6.Google Scholar
- Ober K, Matthews B, Ferrieri A, Kuhn S. The evolution and age of populations of Scaphinotus petersi Roeschke on Arizona Sky Islands (coleoptera, carabidae, cychrini). ZooKeys. 2011; 147:183.PubMedView ArticleGoogle Scholar
- Heinze J, Rueppell O. The frequency of multi-queen colonies increases with altitude in a nearctic ant. Ecological Entomology. 2014;39:4.Google Scholar
- Masta S. Phylogeography of the jumping spider Habronattus pugillis (Araneae : Salticidae): Recent vicariance of sky island populations?Evolution. 2000; 54(5):1699–711.PubMedView ArticleGoogle Scholar
- Smith CI, Farrell BD. Range expansions in the flightless longhorn cactus beetles, Moneilema gigas and Moneilema armatum, in response to pleistocene climate changes. Molecular Ecology. 2005; 14(4):1025–44.PubMedView ArticleGoogle Scholar
- Perez-Alquicira J, Molina-Freaner F, Pinero D, Weller S, Martinez-Meyer E, Rozas J, et al. The role of historical factors and natural selection in the evolution of breeding systems of Oxalis alpina in the sonoran desert sky islands. Journal of Evolutionary Biology. 2010; 23(10):2163–175.PubMedView ArticleGoogle Scholar
- Kupfer JA, John A, Cairns DM. The suitability of montane ecotones as indicators of global climatic change. Progress in Physical Geography. 1996; 20(3):253–72.View ArticleGoogle Scholar
- Kupfer JA, Balmat J, Smith JL. Shifts in the potential distribution of sky island plant communities in response to climate change. USDA Forest Service Proceedings. 2005;:485–90.Google Scholar
- Brusca RC, Wiens JF, Meyer WM, Eble J, Franklin K, Overpeck JT, et al. Dramatic response to climate change in the southwest: Robert Whittaker’s 1963 Arizona mountain plant transect revisited. Ecology and Evolution. 2013; 3(10):3307–319.PubMed CentralPubMedGoogle Scholar
- Garfin G, Jardine A, Merideth R, Black M, LeRoy S. Assessment of Climate Change in the Southwest United States: Island Press; 2013.Google Scholar
- Hölldobler B, Wilson EO. The Ants. Cambridge: The Belknap Press of Harvard University Press; 1990.View ArticleGoogle Scholar
- Andersen AN, Majer JD. Ants show the way down under: invertebrates as bioindicators in land management. Frontiers in Ecology and the Environment. 2004; 2(6):291–8.View ArticleGoogle Scholar
- Agosti D. Ants: Standard Methods for Measuring and Monitoring Biodiversity. Biological diversity handbook series: Smithsonian Institution Press; 2000.Google Scholar
- Diamond SE, Nichols LM, McCoy N, Hirsch C, Pelini SL, Sanders NJ, et al. A physiological trait-based approach to predicting the responses of species to experimental climate warming. Ecology. 2012; 93(11):2305–312.PubMedView ArticleGoogle Scholar
- Pelini SL, Diamond SE, MacLean H, Ellison AM, Gotelli NJ, Sanders NJ, et al. Common garden experiments reveal uncommon responses across temperatures, locations, and species of ants. Ecology and Evolution. 2012; 2(12):3009–015.PubMed CentralPubMedView ArticleGoogle Scholar
- Hölldobler B, Wilson EO. The number of queens: an important trait in ant evolution. Naturwissenschaften. 1977; 64(1):8–15.View ArticleGoogle Scholar
- Heinze J, Keller L. Alternative reproductive strategies: a queen perspective in ants. Trends in Ecology & Evolution. 2000; 15(12):508–12.View ArticleGoogle Scholar
- Peeters C. Convergent evolution of wingless queens reproductive across all subfamilies of ants, and sporadic loss of winged queens (hymenoptera: Formicidae). Myrmecological News. 2012; 16:75–91.Google Scholar
- Peeters C, Ito F. Colony dispersal and the evolution of queen morphology in social hymenoptera. Annual Review of Entomology. 2001; 46(1):601–30.PubMedView ArticleGoogle Scholar
- Heinze J, Buschinger A. Queen polymorphism in Leptothorax Species-A - its genetic and ecological background (hymenoptera, formicidae). Insectes Sociaux. 1989; 36(2):139–55.View ArticleGoogle Scholar
- Roff DA. The evolution of flightlessness in insects. Ecological Monographs. 1990;:389–421.Google Scholar
- Favé M-J. The role of past climate change in driving novelties in sky island ant populations. Montreal: PhD thesis, McGill University; 2013.Google Scholar
- Cremer S, Heinze J. Stress grows wings: Environmental induction of winged dispersal males in cardiocondyla ants. Current Biology. 2003; 13(3):219–223.PubMedView ArticleGoogle Scholar
- Mitchell HK, Petersen NS. Developmental abnormalities in drosophila induced by heat shock. Developmental genetics. 1982; 3(2):91–102.View ArticleGoogle Scholar
- Whittingham MJ, Stephens PA, Bradbury RB, Freckleton RP. Why do we still use stepwise modelling in ecology and behaviour?Journal of animal ecology. 2006; 75(5):1182–9.PubMedView ArticleGoogle Scholar
- Whiting MF, Bradler S, Maxwell T. Loss and recovery of wings in stick insects. Nature. 2003; 421(6920):264–7.PubMedView ArticleGoogle Scholar
- Prud’Homme B, Gompel N, Rokas A, Kassner VA, Williams TM, Yeh SD, et al. Repeated morphological evolution through cis-regulatory changes in a pleiotropic gene. Nature. 2006; 440(7087):1050–3.PubMedView ArticleGoogle Scholar
- Rajakumar R, Mauro D, Dijkstra M, Huang M, Wheeler D, Hiou-Tim F et al. Ancestral developmental potential facilitates parallel evolution in ants. Science. 2012; 335(6064):79–82.PubMedView ArticleGoogle Scholar
- Wund MA, Valena S, Wood S, Baker JA. Ancestral plasticity and allometry in threespine stickleback reveal phenotypes associated with derived, freshwater ecotypes. Biological Journal of the Linnean Society. 2012; 105(3):573–83.PubMed CentralPubMedView ArticleGoogle Scholar
- Losos JB. Convergence, adaptation, and constraint. Evolution. 2011; 65(7):1827–40.PubMedView ArticleGoogle Scholar
- Abouheif E, Favé M-J, Ibarrarán-Viniegra AS, Lesoway MP, Rafiqi AM, Rajakumar R. Eco-evo-devo: The time has come. In: Ecological Genomics. Springer: 2014. p. 107–25.Google Scholar
- Molet M, Wheeler DE, Peeters C. Evolution of novel mosaic castes in ants: Modularity, phenotypic plasticity, and colonial buffering. The American Naturalist. 2012; 180(3):328–41.PubMedView ArticleGoogle Scholar
- (Cohen SM, Bate M, Arias AM, editors.)1993. The Development of Drosophila Melanogaster vol. 1. Cold Spring Harbor Laboratory Pr.Google Scholar
- Weatherbee SD, Halder G, Kim J, Hudson A, Carroll S. Ultrabithorax regulates genes at several levels of the wing-patterning hierarchy to shape the development of the Drosophila haltere. Genes & Development. 1998; 12(10):1474–82.View ArticleGoogle Scholar
- Cripps RM, Black BL, Zhao B, Lien CL, Schulz RA, Olson EN. The myogenic regulatory gene mef2 is a direct target for transcriptional activation by twist during Drosophila myogenesis. Genes & Development. 1998; 12(3):422–34.View ArticleGoogle Scholar
- Held Jr LI. Imaginal Discs: the Genetic and Cellular Logic of Pattern Formation, vol. 39. New-York: Cambridge University Press; 2005.Google Scholar
- Abouheif E, Wray GA. Evolution of the gene network underlying wing polyphenism in ants. Science. 2002; 297(5579):249–52.PubMedView ArticleGoogle Scholar
- Tomoyasu Y, Arakane Y, Kramer KJ, Denell RE. Repeated co-options of exoskeleton formation during wing-to-elytron evolution in beetles. Current Biology. 2009; 19(24):2057–065.PubMedView ArticleGoogle Scholar
- Shbailat S, Khila A, Abouheif E. Correlations between spatiotemporal changes in gene expression and apoptosis underlie wing polyphenism in the ant Pheidole morrisi. Evolution and Development. 2010; 12(6):579–90.View ArticleGoogle Scholar
- Carroll SB, Gates J, Keys DN, Paddock SW, Panganiban G, Selegue JE, et al. Pattern formation and eyespot determination in butterfly wings. Science (New York, NY). 1994; 265(5168):109.View ArticleGoogle Scholar
- Shbailat SJ, Abouheif E. The wing-patterning network in the wingless castes of myrmicine and formicine ant species is a mix of evolutionarily labile and non-labile genes. Journal of Experimental Zoology Part B: Molecular and Developmental Evolution. 2013; 320(2):1–10.View ArticleGoogle Scholar
- Peterson T, Müller GB. What is evolutionary novelty? process versus character based definitions. Journal of Experimental Zoology. Part B: Molecular and Developmental Evolution. 2013; 320(6):345–350.Google Scholar
- Conway Morris S. Life’s Solution: Inevitable Humans in a Lonely Universe. New-York: Cambridge Univ Press; 2003.View ArticleGoogle Scholar
- Gould SJ. Wonderful Life: the Burgess Shale and the Nature of History: Random House; 2000.Google Scholar
- Colosimo PF, Hosemann KE, Balabhadra S, Villarreal G, Dickson M, Grimwood J, et al. Widespread parallel evolution in sticklebacks by repeated fixation of ectodysplasin alleles. Science. 2005; 307.5717:1928–33.View ArticleGoogle Scholar
- Jeong S, Rokas A, Carroll SB. Regulation of body pigmentation by the abdominal-b hox protein and its gain and loss in drosophila evolution. Cell. 2006; 125(7):1387–99.PubMedView ArticleGoogle Scholar
- Protas ME, Hersey C, Kochanek D, Zhou Y, Wilkens H, Jeffery WR, et al. Genetic analysis of cavefish reveals molecular convergence in the evolution of albinism. Nature genetics. 2005; 38(1):107–11.PubMedView ArticleGoogle Scholar
- Sucena E, Delon I, Jones I, Payre F, Stern DL. Regulatory evolution of shavenbaby/ovo underlies multiple cases of morphological parallelism. Nature. 2003; 424(6951):935–8.PubMedView ArticleGoogle Scholar
- Tenaillon O, Rodríguez-Verdugo A, Gaut RL, McDonald P, Bennett AF, Long AD, et al. The molecular diversity of adaptive convergence. Science. 2012; 335(6067):457–61.PubMedView ArticleGoogle Scholar
- Green DA, Extavour CG. Convergent evolution of a reproductive trait through distinct developmental mechanisms in drosophila. Developmental biology. 2012; 372(1):120–30.PubMedView ArticleGoogle Scholar
- Manceau M, Domingues VS, Linnen CR, Rosenblum EB, Hoekstra HE. Convergence in pigmentation at multiple levels: mutations, genes and function. Philosophical Transactions of the Royal Society B: Biological Sciences. 2010; 365(1552):2439–450.View ArticleGoogle Scholar
- Travisano M, Mongold JA, Bennett AF, Lenski RE. Experimental tests of the roles of adaptation, chance, and history in evolution. Science. 1995; 267(5194):87–90.PubMedView ArticleGoogle Scholar
- MacLean RC, Bell G. Divergent evolution during an experimental adaptive radiation. Proceedings of the Royal Society of London Series B-Biological Sciences. 2003; 270(1524):1645–50.View ArticleGoogle Scholar
- Blount ZD, Borland CZ, Lenski RE. Historical contingency and the evolution of a key innovation in an experimental population of escherichia coli. Proceedings of the National Academy of Sciences. 2008; 105(23):7899.View ArticleGoogle Scholar
- McGhee GR. Convergent Evolution: Limited Forms Most Beautiful: MIT Press; 2011.Google Scholar
- Wittkopp PJ, Stewart EE, Arnold LL, Neidert AH, Haerum BK, Thompson EM, et al. Intraspecific polymorphism to interspecific divergence: genetics of pigmentation in drosophila. Science. 2009; 326(5952):540–4.PubMedView ArticleGoogle Scholar
- Albertson RC, Streelman JT, Kocher TD, Yelick PC. Integration and evolution of the cichlid mandible: the molecular basis of alternate feeding strategies. Proceedings of the National Academy of Sciences of the United States of America. 2005; 102(45):16287.PubMed CentralPubMedView ArticleGoogle Scholar
- Kolbe J, Leal M, Schoener T, Spiller D, Losos J. Founder effects persist despite adaptive differentiation: A field experiment with lizards. Science. 2012; 335(6072):1086–9.PubMedView ArticleGoogle Scholar
- Steiner CC, Römpler H, Boettger LM, Schöneberg T, Hoekstra HE. The genetic basis of phenotypic convergence in beach mice: similar pigment patterns but different genes. Molecular Biology and Evolution. 2009; 26(1):35–45.PubMedView ArticleGoogle Scholar
- Brakefield PM, Gates J, Keys D, Kesbeke F, Wijngaarden PJ, Monteiro A, et al. Development, plasticity and evolution of butterfly eyespot patterns. Nature. 1996; 384(6606):236–42.PubMedView ArticleGoogle Scholar
- Jablonski D. Scale and hierarchy in macroevolution. Palaeontology. 2007; 50(1):87–109.View ArticleGoogle Scholar
- Charlesworth B, Lande R, Slatkin M. A neo-darwinian commentary on macroevolution. Evolution. 1982;:474–98.Google Scholar
- Abouheif E. Parallelism as the pattern and process of mesoevolution. Evolution and Development. 2008; 10(1):3–5.PubMedView ArticleGoogle Scholar
- Arthur W. Micro-, macro-, and megaevolution In: Hall BK, Olsen WM, editors. Keywords and Concepts in Evolutional Developmental Biology. Harvard University Press: 2007. p. 496.Google Scholar
- Leroi AM. The scale independence of evolution. Evolution & development. 2000; 2(2):67–77.View ArticleGoogle Scholar
- Meyer JR, Dobias DT, Weitz JS, Barrick JE, Quick RT, Lenski RE. Repeatability and contingency in the evolution of a key innovation in phage lambda. Science. 2012; 335(6067):428–32.PubMed CentralPubMedView ArticleGoogle Scholar
- Erwin DH. Macroevolution is more than repeated rounds of microevolution. Evolution & Development. 2000; 2(2):78–84.View ArticleGoogle Scholar
- Gibson G, Hogness DS. Effect of polymorphism in the drosophila regulatory gene ultrabithorax on homeotic stability. Science. 1996; 271(5246):200–3.PubMedView ArticleGoogle Scholar
- Zimmer C, Emlen DJ. Evolution: Making Sense of Life: Roberts; 2013.Google Scholar
- Hall B. Strickberger’s Evolution: Jones & Bartlett Learning; 2008.Google Scholar
- Ridley M. Evolution, 3rd Edition: Blackwell Publishing; 2003.Google Scholar
- Stern DL. Evolution, Development, & the Predictable Genome. Greenwood Village: Roberts and Company Publishers; 2011.Google Scholar
- Erwin DH, Davidson EH. The evolution of hierarchical gene regulatory networks. Nature Reviews Genetics. 2009; 10(2):141–8.PubMedView ArticleGoogle Scholar
- Davidson EH, Erwin DH. Evolutionary innovation and stability in animal gene networks. Journal of Experimental Zoology Part B: Molecular and Developmental Evolution. 2010; 314(3):182–6.Google Scholar
- Parmesan C. Ecological and evolutionary responses to recent climate change. Annual Review of Ecology, Evolution, and Systematics. 2006;:637–69.Google Scholar
- Seager R, Ting M, Held I, Kushnir Y, Lu J, Vecchi G, et al. Model projections of an imminent transition to a more arid climate in southwestern North America. Science. 2007; 316(5828):1181–4.PubMedView ArticleGoogle Scholar
- Brown JL, Li SH, Bhagabati N. Long-term trend toward earlier breeding in an american bird: A response to global warming?. Proceedings of the National Academy of Sciences. 1999; 96(10):5565–569.View ArticleGoogle Scholar
- Inouye DW, Barr B, Armitage KB, Inouye BD. Climate change is affecting altitudinal migrants and hibernating species. Proceedings of the National Academy of Sciences. 2000; 97(4):1630–33.View ArticleGoogle Scholar
- Bhatkar A, Whitcomb WH. Artificial diet for rearing various species of ants. Fla. Entomol. 1970; 53:229–32.View ArticleGoogle Scholar
- Moreau C, Bell C, Vila R, Archibald S, Pierce N. Phylogeny of the ants: Diversification in the age of angiosperms. Science. 2006; 312.5770:101–104.View ArticleGoogle Scholar
- Zhu Y, Bergland AO, González J, Petrov DA. Empirical validation of pooled whole genome population re-sequencing in Drosophila melanogaster. PloS one. 2012; 7(7):41901.View ArticleGoogle Scholar
- Skalski GT, Couch CR, Garber AF, Weir BS, Sullivan CV. Evaluation of dna pooling for the estimation of microsatellite allele frequencies: a case study using striped bass (Morone saxatilis). Genetics. 2006; 173(2):863–75.PubMed CentralPubMedView ArticleGoogle Scholar
- Vos P, Hogers R, Bleeker M, Reijans M, Vandelee T, Hornes M, et al. Aflp - a new technique for dna-fingerprinting. Nucleic Acids Research. 1995; 23(21):4407–414.PubMed CentralPubMedView ArticleGoogle Scholar
- Bonin A, Ehrich D, Manel S. Statistical analysis of amplified fragment length polymorphism data: a toolbox for molecular ecologists and evolutionists. Molecular Ecology. 2007; 16(18):3737–758.PubMedView ArticleGoogle Scholar
- Zhivotovsky LA. Estimating population structure in diploids with multilocus dominant dna markers. Molecular Ecology. 1999; 8(6):907–13.PubMedView ArticleGoogle Scholar
- Vekemans X, Beauwens T, Lemaire M, Roldán-Ruiz I. Data from amplified fragment length polymorphism (aflp) markers show indication of size homoplasy and of a relationship between degree of homoplasy and fragment size. Molecular ecology. 2002; 11(1):139–51.PubMedView ArticleGoogle Scholar
- Mantel N. Detection of disease clustering and a generalized regression approach. Cancer Research. 1967; 27:209–20.PubMedGoogle Scholar
- Beerli P. Comparison of bayesian and maximum-likelihood inference of population genetic parameters. Bioinformatics. 2006; 22.3:341–345.View ArticleGoogle Scholar
- Beerli P, Palczewski M. Unified framework to evaluate panmixia and migration direction among multiple sampling locations. Genetics. 2010; 185(1):313–463.PubMed CentralPubMedView ArticleGoogle Scholar
- Haag-Liautard C, Coffey N, Houle D, Lynch M, Charlesworth B, Keightley PD. Direct estimation of the mitochondrial dna mutation rate in Drosophila melanogaster. PLoS biology. 2008; 6(8):204.View ArticleGoogle Scholar
- Obbard DJ, Welch JJ, Kim KW, Jiggins FM. Quantifying adaptive evolution in the drosophila immune system. PLoS genetics. 2009; 5(10):1000698.View ArticleGoogle Scholar
- Cornuet JM, Ravigne V, Estoup A. Inference on population history and model checking using dna sequence and microsatellite data with the software diyabc (v1.0). Bmc Bioinformatics. 2010; 11(401).Google Scholar
- Beaumont MA, Zhang WY, Balding DJ. Approximate bayesian computation in population genetics. Genetics. 2002; 162(4):2025–035.PubMed CentralPubMedGoogle Scholar
- Felsenstein J. Phylip-phylogeny inference package (version 3.2). Distributed by the author: Department of Genome Sciences, University of Washington, Seattle; 2004.Google Scholar
- Correa C, Hendry AP. Invasive salmonids and lake order interact in the decline of puye grande Galaxias platei in western patagonia lakes. Ecological Applications. 2012; 22(3):828–42.PubMedView ArticleGoogle Scholar
- Wheeler DE, Nijhout HF. Imaginal wing discs in larvae of the soldier caste of Pheidole bicarinata vinelandica forel (hymenoptera: Formicidae). International Journal of Insect Morphology and Embryology. 1981; 10(2):131–9.View ArticleGoogle Scholar
- Patel N. Imaging neuronal subsets and other cell-types in whole-mount drosophila embryos and larvae using antibody probes. In: Methods in Cell Biology. San Diego: Academic Press Inc.: 1994. p. 445–87.Google Scholar
- Warton DI, Wright IJ, Falster DS, Westoby M. Bivariate line-fitting methods for allometry. Biological Reviews. 2006; 81(2):259–91.PubMedView ArticleGoogle Scholar
- Tautz D, Pfeifle C. A non-radioactive insitu hybridization method for the localization of specific rnas in drosophila embryos reveals translational control of the segmentation gene hunchback. Chromosoma. 1989; 98(2):81–5.PubMedView ArticleGoogle Scholar
- Metscher BD. Microct for comparative morphology: simple staining methods allow high-contrast 3d imaging of diverse non-mineralized animal tissues. BMC physiology. 2009; 9(1):11.PubMed CentralPubMedView ArticleGoogle Scholar
- Mani MS, Giddings LE. Ecology of Highlands: Dr. W. Junk.; 1980.Google Scholar
- Hodkinson ID. Terrestrial insects along elevation gradients: species and community responses to altitude. Biological Reviews. 2005; 80(3):489–513.PubMedView ArticleGoogle Scholar
- Harrison RG. Dispersal polymorphisms in insects. Annual Review of Ecology and Systematics. 1980;:95–118.Google Scholar
- Peeters C, Molet M. Colonial reproduction and life histories In: Lach L, Parr C, Abbott K, editors. Ant Ecology. New-York: Oxford University Press: 2010. p. 159–76.Google Scholar
- Murcia C. Edge effects in fragmented forests: implications for conservation. Trends in Ecology & Evolution. 1995; 10(2):58–62.View ArticleGoogle Scholar