- Open Access
The high Andes, gene flow and a stable hybrid zone shape the genetic structure of a wide-ranging South American parrot
Frontiers in Zoology volume 8, Article number: 16 (2011)
While the gene flow in some organisms is strongly affected by physical barriers and geographical distance, other highly mobile species are able to overcome such constraints. In southern South America, the Andes (here up to 6,900 m) may constitute a formidable barrier to dispersal. In addition, this region was affected by cycles of intercalating arid/moist periods during the Upper/Late Pleistocene and Holocene. These factors may have been crucial in driving the phylogeographic structure of the vertebrate fauna of the region. Here we test these hypotheses in the burrowing parrot Cyanoliseus patagonus (Aves, Psittaciformes) across its wide distributional range in Chile and Argentina.
Our data show a Chilean origin for this species, with a single migration event across the Andes during the Upper/Late Pleistocene, which gave rise to all extant Argentinean mitochondrial lineages. Analyses suggest a complex population structure for burrowing parrots in Argentina, which includes a hybrid zone that has remained stable for several thousand years. Within this zone, introgression by expanding haplotypes has resulted in the evolution of an intermediate phenotype. Multivariate regressions show that present day climatic variables have a strong influence on the distribution of genetic heterogeneity, accounting for almost half of the variation in the data.
Here we show how huge barriers like the Andes and the regional environmental conditions imposed constraints on the ability of a parrot species to colonise new habitats, affecting the way in which populations diverged and thus, genetic structure. When contact between divergent populations was re-established, a stable hybrid zone was formed, functioning as a channel for genetic exchange between populations.
Current molecular genetic methods allow the understanding of the genetic structure underlying different populations of a species with previously unforeseen resolution [e.g. [1–3]]. As a result, it is possible to undertake fundamental investigations in ecology and evolution, like the study of the influence of past and current environmental conditions, together with ecological barriers, in shaping the population structure of wild animal species. These studies provide a unique opportunity to understand how species have evolved and how they are organised across landscapes . The constraints that heterogeneous landscapes (e.g. barriers, resource distribution) and environmental conditions (e.g. climate) impose on the ability of animals to colonise new habitats have genetic implications affecting the structure, dynamics and persistence of populations [e.g. [5–9]]. Thus, significant genetic structuring can be expected among populations where gene flow is restricted [e.g. [5, 10]].
However, when contact between divergent populations is re-established, hybrid zones can form with resultant consequences for the evolutionary trajectories of the interbreeding populations [e.g. ]. Two scenarios could potentially occur: if hybrid fitness is high, introgression will be widespread and hybridising populations may become panmictic over time, replacing the original populations ; but if interbreeding is limited in geographical range, hybridising populations may experience genetic exchange without panmixia [e.g. ]. Under such circumstances hybrid zones may persist over time and function as channels for genetic exchange between the populations, increasing overall levels of genetic and phenotypic diversity [e.g. [14–17]].
Although Pleistocene climate conditions played an important role in initiating major phylogeographical structuring in today's fauna , very little is known about their effects on southern South American terrestrial vertebrates. Phylogeographic studies in this region are scarce and have been conducted on rodents [19–21], lizards [22–24], amphibians [9, 25], and a bird species . Several of these studies suggested that the phylogeographic patterns observed, such as past fragmentation, range expansion, and secondary contact, could in part be understood in light of Pleistocene climate conditions. There is evidence of several cycles of more arid conditions intercalated with moist periods during the Pleistocene and Holocene of Southern South America [e.g. [27–30]], which were related to glaciation events , influencing the vegetation distribution  and the fauna depending on it.
One species that may allow hypotheses about gene flow both across a heterogeneous distribution, and into hybrid zones to be tested is the burrowing parrot (Cyanoliseus patagonus) (Aves, Psittaciformes). This species is distributed in a particularly heterogeneous arid to semi-arid landscape, across an extensive ~1,000,000 km2 range in Chile and Argentina. A previous study  suggested that precipitation and temperature restrict burrowing parrot distribution in Argentina. According to this study , burrowing parrots are restricted to an area with median annual precipitation up to 600 mm, and annual average temperatures of no less than 8°C. However, this topic merits further research, as the study  was based on outdated distributional data and on the plain interpretation of maps, without detailed statistical analyses. The most dominant feature of this region is the high Andes, a mountain range that attains an altitude of up to 6,900 m and appears to separate burrowing parrot populations in Chile from those in Argentina. The predominant ecosystem on the Chilean side of this region is the 'Matorral', where vegetation is adapted to the generally dry conditions of a Mediterranean climate zone [34–36]. On the Argentinean side the semi-desert scrubland known as the 'Monte' is predominant, and this occurs from Patagonia to the North-west of Argentina [ and references therein]. Even though there are no extrinsic barriers to parrot dispersal in Argentina, more than 2,300 km separate the southernmost and northernmost burrowing parrot populations there. In addition, burrowing parrots breed in sandstone, limestone or earthen cliffs or "barrancas" (gorges or ravines), where they excavate nest burrows and form colonies [e.g. ]. These cliffs are heterogeneously located in the driest parts of the burrowing parrots range, being commonly found along permanent or temporary rivers, lakeshores, and the seacoast. Given the dryness of this environment, burrowing parrot colonies are never far from freshwater on which they are completely dependent, as they need to drink several times per day [37, 38] (Figure S1). These specific requirements for nest sites, which are spread over thousands of square kilometres, and water, together with the colossal barrier of the Andes, may favour the isolation of burrowing parrot breeding sites and a complex population structure driven by genetic drift.
Due to the heterogeneity of habitats within this species' range, four burrowing parrot sub-species have been proposed, three of which are found in Argentina: C. p. patagonus in Patagonia, C. p. andinus in the Cuyo region to the west and north-west, and C. p. conlara ranging in the San Luis region between the former two (hereafter patagonus, andinus, conlara) [34, 39, 40]. The sub-species C. p. bloxami (hereafter bloxami) is found on the Andean foothills of Central Chile [34, 41]. Three of the sub-species, namely andinus, patagonus and bloxami, are clearly morphologically distinct (size and plumage coloration) , while some authors [34, 40, 42] considered conlara a hybrid, owing to its intermediate geographic location and phenotype, between patagonus and andinus. Little is known about the genetic structure of burrowing parrots and how this corresponds to the morphological sub-species described above. A previous study  attempted to address this using seven microsatellite markers and suggested moderate differentiation between bloxami and all other subspecies, but differentiation within Argentinean samples was not detected. Analysis of a larger sample using a uniparental marker such a mitochondrial DNA may increase the resolution of genetic structure in this species.
Burrowing parrots are currently threatened by intense collection of birds for the pet trade , unjustified persecution as a crop pest [41, 44] and strong habitat loss and degradation, particularly in the Monte ecosystem . The latter could strongly reduce connectivity among the populations, enhancing isolation. As key species in the Monte, any negative impact on burrowing parrots could potentially affect other species since their abandoned and semi-collapsed nests provide breeding space to many other cavity nesters (such as insects, reptiles, birds and small mammals) .
Given the marked phenotypic differences between the populations on both sides of the Andes we hypothesise that gene flow across this high mountain range, the largest barrier in the region, must be severely restricted. We tested this hypothesis using three mtDNA loci, in a large scale sampling effort covering almost the entire species range. We also aimed to uncover the underlying population structure of this species, determine their geographic origins and suggest possible routes of colonisation. We also used these data to determine if an andinus-patagonus hybrid zone exists. Lastly, considering the conservation value of this key species, its potentially restricted distribution with respect to climatic factors, and the unchecked degradation of their preferred habitats, we aim to ascertain the extent to which ecological and climatic factors influence their population structure.
Recently, various aspects of the breeding biology of this species have been investigated extensively, providing the necessary framework for this study [e.g. [37, 38, 47–57]]. Fieldwork was carried out from November to December 2007 (Argentina), February 2008 (northern Chile) and from October to November 2008 (Argentina). Thirty-three colonies and eleven roosting places of the four previously proposed sub-species (bloxami, andinus, conlara and patagonus) were visited and naturally moulted feathers were non-invasively collected (Figure 1). Since burrowing parrots moult their primary feathers at the beginning of their breeding season, from November onwards [38, 58], collection is usually straightforward, as feathers tend to accumulate at the bottom of the cliffs with colonies. Taxonomic assignment of colonies was conducted following previous studies [34, 39]. The southernmost populations of bloxami in Chile were not accessible.
DNA was extracted from the feather quill, for a subset of 3 - 29 individuals from each of 44 populations using the DNAeasy blood and tissue kit (Qiagen, Germany). Mitochondrial genes cytochrome b (cytb), cytochrome oxidase subunit I (COI) and ATPase subunits 6/8 were amplified via PCR using the following mitochondrial-specific primers: L15424 5'- ATCCCATTCCACCAT ACTACTC, H15767 5'- ATGAAGGGATGTTCTACTACTGGTTG-3'( cytb, 586bp) , COI-F 5'- CCTGCAGGAGGAGGAGAYCC-3', COI-A 5'- AGTATAAGCGTCTGGGTAGTC-3' (COI, 455bp) , and CO2GQL 5'- GGACAATGCTCAGAAATCTGCGG-3', CO3HMH 5'- CATGGGCTGGGGTCRACTATGTG-3' (ATP6/8, 818bp) . PCRs were conducted in 20µl reaction volumes containing 100 ng DNA template, 10 mM of each primer, 10 mM dNTPs (Roth, Karlsruhe), 3.125 mM MgCl, 5 U Taq Polymerase (Qiagen Taq Polymerase Core Kit) in a 1× PCR reaction buffer. PCR commenced with denaturation at 94°C for 3 minutes, followed by 35 cycles of denaturation at 95°C for 45 s, annealing at 55°C (CO1, ATP6/8) or 52°C (cytb) for 45 s and extension at 72°C for 45 s. A final extension step (10 min at 72°C) concluded the PCR. Products were checked on a 1% agarose gel and purified of primers and excess dNTPs using either exonuclease-shrimp alkaline phosphatase (Fermentas Life Sciences) or the MinElu-te PCR Purification Kit (Qiagen), following the manufacturers specifications. PCR products were then sequenced in both directions by either Qiagen (Qiagen GmbH, Sample & Assay Technologies, Hilden, Germany) or by using Big Dye chemistry (Applied Biosystems) and run on an AB 3130xl genetic analyser (Applied Biosystems). Sequences were checked, edited and aligned in CLC DNA Workbench 5.5 (CLC bio). A total of 150 individuals from 41 of the attempted 44 locations were suitable for downstream analyses (Table 1, 2, 3 and 4; see also Additional file 1 Table S1).
Genetic variation was quantified as the number of haplotypes, nucleotide and haplotype diversity per population, and was determined using DnaSP v5 . The degree of population structuring was ascertained using a Bayesian population assignment model that assumes no-admixture and that loci are linked BAPS 5 . Unlike other clustering algorithms, this BAPS module allows the assumption of linkage between loci, thereby enabling population assignments using multilocus mtDNA data. The simulation was run ten times for Kmax values of five, ten and 20 potential populations.
Phylogenetic structuring among haplotypes was investigated by maximum likelihood (ML) in Treefinder [64, 65], which first determined the best-fit substitution model for each gene partition. ML analyses was carried out assuming the HKY+G (cytb), HKY+I (CO1) and J2+G (CO2/3) models, each optimising rate, frequency and heterogeneity parameters directly from the data. Models including a rate-heterogeneity parameter (+G) assumed five gamma categories. The significance of nodal bi-partitions was determined by 1,000 bootstrap replicates, from which an 80% majority rule consensus tree was constructed. We also constructed a median-joining network using Network 18.104.22.168  for a graphical representation of the unrooted relationships and frequency of haplotypes.
Demographic parameters based on the coalescent can be useful in testing hypotheses of population history. In particular, statistics based on the distribution of pair-wise differences (mismatch distribution)  between individuals in a population as well as the detection of selection among selectively neutral loci may be used as signatures of past population expansion events. We used Arlequin 3.5  to calculate the mismatch distributions and Fu's Fs for all populations identified by the clustering and phylogenetic methods above.
We roughly dated each node using a fossil calibration technique that incorporates rate smoothing for lineages with unequal rates of evolution. By constraining the tree topology to that of the 80% majority rule consensus tree calculated previously, we generated a range of node height estimates by running a bootstrap analyses 1,000 times in Treefinder, using the original starting parameters. The variance in node height was then taken as the spread of this variable in our data set. The fossil record for C. patagonus is extremely scarce, with three known fossils dated to 126 kyr [69–72]. A much earlier fossil that can also be attributed to the genus Cyanoliseus dates back to 750 kya [70–73]. However, since the nearest outgroup taxon for which molecular data were available belonged to the genus Diopsittaca, it was not possible to utilise the older fossil date in our calibration. We therefore assumed a minimum or latest date of 126 kyr for the coalescence of all C. patagonus lineages. We then performed the calibration analysis in Treefinder and, due to the availability of a single calibration date, we used local rate minimum deformation rate-smoothing, to account for the possibility of differing rates of lineage evolution within the phylogeny. 95% confidence intervals were generated from the spread of node heights.
We investigated the possibility of differing mutation rates among the three gene partitions in order to independently date each node using available mutation rates. However, since reliable mutation rates are currently only available for the avian cytochrome b gene , we partitioned the sequence data into cyt b and CO1+CO2/3 and performed identical analyses in BEAST 1.5.0 , using the parameter-rich GTR+G+I model (with 5 gamma categories) to check if cyt b mutation rates may be applied to other mitochondrial data. Plotting the resulting cyt b node heights against those of CO1+CO2/3 returned a correlation co-efficient (R2) value of 0.5858 (p < 0.001, data not shown), suggesting significantly differing rates of evolution between the two gene partitions. As a rate-smoothing option is unavailable in BEAST, we proceeded with only the cytochrome b data, using both the average avian mutation rate (µ) of 2.1%/million years (myr) and a Psittaciform-specific rate of 3.4%/myr .
Given the possibility of a relatively recent divergence and the apparent reliance of burrowing parrots on their habitat and associated climate [33, 37], we also investigated the extent to which taxonomic, ecological and climatic factors influenced the distribution of genetic heterogeneity in the data. We fitted a linear model to our data using DISTLM 5 , thereby testing the influence of 32 variables (4 taxonomic, 7 ecological and 21 climatic; Table 5 and Table S1) using multiple matrix regressions. The climatic variables were selected taking into account the possible restriction of burrowing parrots distribution to certain temperature and precipitation ranges as suggested in a previous study . The 21 temperature- and precipitation-related climatic variables were obtained, for each of the burrowing parrot colonies, with DIVA-GIS 22.214.171.124. [http://www.diva-gis.org], which are based on the WorldClim database [http://www.worldclim.org], version 1.3, at 2.5 minutes resolution . We first assessed the marginal genetic variation explained by each set of variables separately. Since some of these variables are likely to co-vary with geography (latitude and longitude of sample locations), we controlled for this influence by also reporting the conditional residual variation remaining after the influence of geography was subtracted. Still controlling for geography, we then tested all variables together using a forward selection approach , which sequentially determines the variables that explained the majority of the marginal variation in the data. The significance of all DISTLM regressions was determined by 9,999 permutations.
To explore the influence of the geographic landscape in more detail, we used a spatial clustering algorithm, implemented in Tess 2.3.1 . This method performs Bayesian clustering given a set of input sampling locations. The program constructed Dirichlet cells around each sampling location to produce a Voronoi tessalation. Multilocus DNA sequences were then used to statistically infer population structure based on this Voronoi neighbourhood system using a hidden Markov random field prior. Parameter space was explored by Markov chain simulations. We ran the program ten times again with a Kmax of five, ten and 20, assuming no-admixture at 50,000 sweeps each, with a burn-in of 10,000 sweeps. We also included 23 dummy points, well outside the extant distribution C. patagonus, in order to restrict the simulation to within the species range. The resulting hard-clustering image represents the geographical clustering of individuals in landscape space, given their sampling origins and their multilocus genotypes. To allow the optimal tessalation to be viewed in geographic space, the hard clustering image with the lowest deviance information criterion (DIC) score was imported into ArcGIS 9.3, and superimposed onto a map of the region.
Complete sequences for all three mitochondrial fragments were obtained for 150 individuals from 41 locations (Tables 1, 2, 3 & 4, and Additional file 1 Table S1), representing 7 ecoregions and encompassing the entire species range. However, since the quality of DNA isolated from shed feathers was low, this number was less than half of the 327 individuals from which DNA was extracted. Genetic variation was also low, with only 81 polymorphic sites from a concatenated sequence of 1,859 bp, resulting in 51 unique haplotypes (Tables 1, 2, 3 & 4, and Additional file 2 Table S2). This was reflected in a species-wide haplotype diversity of 0.943 and nucleotide diversity of 0.00530 (Table 4), and given varying locality samples sizes these values ranged from 0.500 - 1.000 and from 0 - 0.0360 respectively (Tables 1, 2, & 3). Considering taxonomic designations, conlara was found to have the highest nucleotide diversity. Non-conlara sampling localities, contained similar levels of diversity, with the exception of eight localities with a nucleotide diversity of less than 0.001. This is reduced to only four localities, if we exclude those with fewer than 3 individuals.
Despite low diversity, Bayesian clustering consistently structured the entire sample into four population clusters regardless of the Kmax prior, and the distributions of these are depicted in Figure 1. One of the four populations corresponded exactly to the bloxami phenotype, and was found exclusively on the Chilean side of the species range (inset, Figure 1). Within Argentina, only two members of the northern andinus phenotype did not fall within an Andinus population cluster, and the patagonus phenotype was divided into two populations, hereafter called Patagonus1 and Patagonus2. There was no support for the existence of the intermediate conlara, as this phenotype clustered either within the Andinus or Patagonus2 populations, with a single conlara clustered within Patagonus1.
Maximum likelihood analysis recovered a fully resolved haplotype phylogeny (Figure 2A) again showing four population groups, but with Bloxami as basal and distinct from all other populations. Within Argentina, the northern Andinus forms a sister relationship with the entire Patagonus population, with the two most ancestral haplotypes of the latter sampled in the Cuyo region (locations 3 to 9, and 16, Figure 1), among phenotypically andinus populations. In Patagonia, the patagonus subspecies is divided into the genetically distinct populations Patagonus1 and Patagonus2, yet without distinguishing phenotypic characteristics. A median-joining haplotype network (Figure 2B) confirmed the large divergence separating populations on either side of the Andes, but also demonstrated that all three Argentinean populations contain star-shaped haplotype clusters in which several less-frequent haplotypes are very closely related to a single common haplotype. Both phylogeny and network showed that individuals identified as conlara or "undetermined" (where phenotype could not be assigned to andinus or patagonus with certainty) belonged to either the Andinus (11/33 individuals; 5/14 haplotypes) or Patagonus (Patagonus1: 1/55 individuals, 1/15 haplotypes; Patagonus2: 18/51 individuals, 6/16 haplotypes) populations (Figure 2 A/B), confirming a hybrid origin for this phenotype. The six undetermined individuals (5 haplotypes) that were sampled within the conlara range from a population near the town of San Luis and from the Quinto river were most closely related to confirmed conlara haplotypes on the median-joining network (Figure 2B), showing that both populations are admixed and therefore belong within conlara.
Mismatch distributions of pair-wise nucleotide differences (Figure 2C) showed largely unimodal distributions among Argentinean populations, but was multimodal within the Chilean Bloxami. Bloxami was also the only population where the hypothesis of a sudden population expansion was rejected (p(SSD) = 0.0290). Furthermore, negative and significant Fu's Fs values, also indicating population expansion, were recovered for Andinus, Patagonus1 and Patagonus2, with the slightly negative value for Bloxami not significant. Assuming a fossil calibration of at least 126 kyr for the coalescence of all C. patagonus lineages, we estimated dates for all population nodes to between 9 and 55 kya (Additional file 3 Table S3). Applying the average avian mutation rate (µ = 2.1%/myr) to the cytochrome b data returned higher nodal divergence estimates with much larger 95% confidence intervals. However, the psittaciform-specific mutation rate (µ = 3.4%/myr) yielded estimates very similar to those of the fossil calibration method. Owing to the larger variance of BEAST estimates, the 95% interval of almost all fossil calibration dates was contained entirely within the psittaciform-rate interval (Additional file 3 Table S3). The only exception being node E, although both fossil and µ estimates did overlap. Indeed, if a divergence time (T) of 126 kyr for all burrowing parrots was substituted into the equation µ = node height/2T, a mutation rate of 3.5%/myr would result, implying that the cytochrome b gene evolves more rapidly among psittaciform species. The 95% confidence intervals of the fossil calibrated population nodes reflect the stochasticity in the maximum likelihood estimation and the lack of internal calibration nodes. Their accuracy, therefore, relies heavily on that of the 126 kyr fossil estimate, for which there exist three independently dated fossils [69–72]. In contrast, the much larger variance of BEAST estimates may reflect the limitations of using only a third of the total available data. We therefore continue with the more accurate fossil calibrated estimates (Figure 2A), and we view these as the latest range of dates by which each lineage splitting event could have occurred.
Influences on genetic variation
Since extant burrowing parrot populations appear to have evolved relatively recently, one might expect to find genetic structure to be associated to present-day ecological and environmental variables. The amount of marginal variation explained by predictor variables was highest for the suite of 21 climatic variables (85%) followed by phenotype (22%) and lastly ecoregions (0.81%). This pattern held true even when the variation attributable to geography was removed in conditional tests (Table 5). When each climatic variable was tested individually, temperature seasonality, isothermality, temperature annual range, mean temperature of warmest quarter and maximum temperature of warmest month could all account for more than 15% of the conditional variation in the data set. When autocorrelations between these variables were taken into account in sequential tests, only temperature seasonality (28%) and precipitation seasonality (11%) were found to explain more than 10% of the genetic variation (Table 5).
The 44% of genetic variation explained by geography was partitioned into the four already-defined population clusters by TESS software. The optimal tessellation cluster was obtained consistently for all runs and for all starting Kmax values and this partitioned Chilean from Argentinean populations, as have previous analyses (Figure 2). Spatial structuring within Argentina describes genetically Andinus haplotypes in the Cuyo region, to the north of the species range, whereas in Patagonia, Patagonus1 inhabits the entire Atlantic coast and a thin wedge between southern limit of Andinus and the north of the Río Colorado valley (Figure 3). The Patagonus2 population occupies a non-contiguous area along the entire course of the Río Colorado and except for a population in which Patagonus1 is particularly dominant, also the entire course of the Río Negro until just before it enters the Atlantic Ocean. Some Patagonus2 haplotypes also range north into the San Luis region, which is the southern-eastern end of the Andinus distribution, where hybridization with Andinus results in the conlara phenotype.
Origins and gene flow across the Andes
Although we were not able to access populations to the south of the Bloxami range in Chile, it is nevertheless apparent that very limited gene flow across the Andes has rendered the Bloxami populations both genetically and phenotypically distinct from all other burrowing parrots. This further corroborates a previous study using microsatellite loci that showed a clear separation between Chilean and Argentinean populations, but no structure within Argentina . By contrast, cluster analyses of 3 mitochondrial fragments suggest a complex population structure for burrowing parrots in Argentina. Since the Bloxami split is at the root of the species phylogeny, and that this is the only population to not bear any of the signatures of a recent population expansion (Figure 2B/C), we suggest a Chilean origin for this species, with a single migration event across the Andes that gave rise to all extant Argentinean mitochondrial lineages (Figure 3). Although the sampled Chilean lineages coalesce to only 12 - 34 kya, their southern range was not sampled and the Bloxami mismatch distribution (Figure 2C) suggests a population with no significant changes in effective population size over time. Furthermore, the network of Bloxami lineages (black circles, Figure 2B) describes a population with several haplotypes in low frequency, none of which are related to each other by fewer than two mutational steps, implying that intermediate and other more distant haplotypes, that could have increased the Bloxami coalescence time, were either lost or not sampled.
In contrast, the histories of all populations in Argentina are characterised by recent (between 9 - 55 kya) expansion (Figure 2 B/C). The migration event across the Andes therefore dates roughly to between 55 - 126 kya, to a time when the southern and central reaches of that mountain range were heavily glaciated [e.g. [28, 31, 32, 80]], making a crossing over the lower (< 2,500 m) passes in these areas impossible for want of water and shelter until at least 14 kya . Instead, we suggest a crossing over one of the few intermediate (2,500 - 4,000 m) mountain passes to the north of the present-day species distribution (Figure 3). Although rare, recent evidence of exactly such a crossing from Chile into the Mendoza region of Argentina by a Peruvian pelican (Pelecanus thagus), where all passes are higher than 2,500 m, shows that bird migration across the high Andes is biologically possible . The resulting founder effect of an Andean crossing may also help to explain the contrasting size and plumage differences that exist between bloxami and andinus.
Gene and phenotype divergence in Argentina
The initial divergence of Andinus-Patagonus mitochondrial lineages during the late Pleistocene 34 - 55 kya represents a clear asynchrony between genetic and phenotypic divergence times within this species. The two most ancestral haplotypes assigned to the Patagonus1 population (Figure 2A/B) are phenotypically andinus individuals, sampled among andinus locations in the Cuyo region, at least 500 km from the andinus-patagonus hybrid zone in San Luis. Since older Patagonus haplotypes were not detected, despite a sample of 106 patagonus individuals from across their entire distribution, we propose that phenotypic divergence between Andinus and Patagonus lagged genetic divergence, by 3 - 41 kyr, which is the minimum-maximum time difference between the initial Andinus- Patagonus and the Patagonus1-Patagonus2 divergences. Further evidence that phenotypic differences take longer to evolve in this species is the genetic diversification of Patagonus1 from Patagonus2 within Patagonia, which occurred no later than 9 kya, yet the two populations still remain phenotypically indistinguishable from each other. This lag in phenotypic divergence may result from greater lineage sorting among mitochondrial genomes, and whether it holds true for slower evolving unlinked regions such as nuclear introns remains to be investigated.
Colonisation and diversification
Summarising the results and available information, we propose a model for the colonisation and diversification events to have shaped the distribution of genetic variation in burrowing parrots (Figure 3). The presence of ancestral Patagonus haplotypes within the range of andinus implies that genetic diversification of these two populations occurred in the Cuyo region, with phenotypic diversification occurring later (after 21 kya, see Figure 2A) as Patagonus expanded south-east (see model, Figure 3). Since burrowing parrots are completely dependent on water and suitable strata in which to form nests [e.g. ], we predicted the population structure to be a product of isolation and hence local genetic drift. Our results suggest, rather, a low degree of isolation since most sampling localities are similarly diverse, indicating gene flow between them. This is most likely owing to suitable habitats being continuously distributed along river courses and because burrowing parrots are able to make one to four foraging trips to a distance of up to 66 km from their colonies in a single day . Hence, we conclude that the south-eastern expansion of genetically Patagonus individuals probably followed the courses of the Bermejo-Desaguadero-Atuel, Colorado and Negro river systems, which were already in existence by that time [28, 82]. By 14 - 31 kya, Patagonus was established along river courses, and along the Atlantic coast of Patagonia (Figure 3). The end of the last glacial maximum (LGM) ~14 kya saw the establishment of wetter, more suitable habitats in north-western Patagonia [e.g. [27–30, 80]] and possibly promoted further diversification into Patagonus2 between 9 - 21 kya. Spatial clustering shows that the Holocene expansion of Patagonus2 lineages has rendered this population dominant in much of north-western Patagonia, between the Andes and the Atlantic (Figure 3), especially along and to the south of the Río Negro, where only few populations contain Patagonus1 haplotypes in high frequency.
Our results convincingly demonstrate the existence of a hybrid zone in the San Luis region, where Patagonus (Patagonus2 in all but one case) haplotypes expanding into the south-eastern range of Andinus resulted in the evolution of an intermediate phenotype: conlara, and this is also the most genetically diverse of the four burrowing parrot taxonomic groupings. The dynamics of hybrid zones are of great interest, especially their potential to give rise to new species or populations [e.g. ]. While our results warrant a thorough study of the Andinus-Patagonus hybrid zone using nuclear DNA markers, mitochondrial data reveal several processes of interest. We found no significant difference between the proportion of conlara individuals belonging to Andinus and Patagonus2 (χ2 = 0.003, df = 1, P = 0.960) and Andinus and all Patagonus (χ2 = 2.7, df = 1, P = 0.102) suggesting, firstly that there is no sex bias in dispersing Patagonus individuals, and secondly that there was no bias by resident Andinus against the choice of invading Patagonus individuals as potential mates. From our phylogeny (Figure 2A), we also conclude that the observed level of introgression resulted from at least four, possibly independent, migration events by Patagonus individuals. Within the Patagonus2 clade, four conlara-undetermined haplotypes comprise an entire sub-clade, suggesting that they resulted from a single successful migration event, and given the depth of this sub-clade, Patagonus2 introgression may date back to at least the latter half of the Holocene (inference, Figure 2A). Similarly, the most ancestral of Andinus haplotypes was carried by a conlara individual and dates back to the coalescence of the Andinus clade (14 - 34 kya). Together with the morphological conformity of the conlara phenotype within its range , this result strongly suggests a hybrid zone that has remained stable for several thousand years. The inability of seven polymorphic microsatellite loci to discriminate among Argentinean populations  also suggests a hybrid zone that has remained stable long enough for homoplasy to mask phylogeographic signal in hypervariable markers. Since phenotypic variation across the andinus-conlara-patagonus continuum, appears to be discretely, rather than continuously partitioned, implying that gene flow out of the hybrid zone is lower that gene flow into it, the possibility of the conlara phenotype being selectively advantageous in the San Luis region cannot be ruled out.
The effect of climate
Our data show that up to 48% of the genetic variation can be accounted for by present-day climatic variables, considerably greater than the variation accounted for phenotype and ecoregions. This is even more surprising since ecoregions would normally be expected to co-vary with climate. This lack of variation attributable to ecoregion might be explained by the fact that the two major constraints on burrowing parrots habitat suitability, namely availability of water and cliffs, are not included in ecoregion definitions. Even the significant correlation with present-day climate appears counterintuitive since climatic conditions in the past are more likely to have had an influence over population structure. However, if we consider that the majority of the haplotypes within Patagonus1/2 and Andinus most likely evolved after the LGM 14 kya or more recently (Figure 2A), and that local climatic conditions have remained relatively similar in the last 8 kyr [e.g. [28–30, 80]], this result is less surprising. Conservation implications of this are particularly important in the present-day reality of climate change.
The lack of gene flow makes the high Andes an important barrier to migration in burrowing parrots, and possibly in other bird species [e.g. ]. This also renders the isolated Bloxami population genetically and phenotypically distinct. This evolutionary significance is important from a conservation and management perspective. Burrowing parrots are listed as ''threatened'' species in the vertebrate red list of Chile and as such are legally protected . This is because only 5,000 - 6,000 bloxami are distributed in the IV and VII regions of the country [41, 85]. The size and uniqueness of this population means that further reductions should be avoided.
The situation differs considerably in Argentina, where burrowing parrots are officially considered an abundant agricultural pest (National Law of Defence of Agricultural Production 6704/63), despite agricultural damage occurring only locally [40, 86], with very little actual crop damage recorded [37, 44]. Owing to its persecution as a pest species [ and references therein], several colonies have been destroyed or severely reduced in size, including the extirpation of the largest known colony of some 50,000 nests . Collection of burrowing parrots for the pet trade was also encouraged  and population reductions continued to reach levels considerable enough for the regional government of the Río Negro province to ban all hunting and trade (resolutions 23-DF-2004, 24-DF-2004, Dirección de Fauna de la Provincia de Río Negro, Argentina), thereby extending legal protection to all but seven Patagonian colonies. This protection effectively includes the bulk of the Patagonus population, approximately 40,000 nests, where genetic diversity is partitioned into two genetically distinct, yet phenotypically indistinguishable populations, which are impossible to manage separately. It should be noted, however, that 37,000 of these nests are located in a single colony - El Cóndor (sampling location 39, Figure 1) , which is located in an area undergoing habitat degradation that is estimated to be ten times higher (3.7% annually)  than the world average of 0.4% . The continued existence of the burrowing parrot in Patagonia is therefore uncertain.
A negative side effect of the recent protection of Patagonian populations is the noticeable increase in commercial value of Andinus populations in the Cuyo region. All available data from this study (Tables 1, 2 & 3) and the literature , together with unpublished data of 1 colony for which we did not obtain sequences, La Manga stream, La Rioja, 290 nests) show that the total Andinus population numbers no more than 2,000 nests. In contrast to Patagonus, where individuals are genetically but not phenotypically distinguishable, Andinus populations are distinctive both genetically (Figures 2 and 3) and phenotypically  from Patagonus, comprising an evolutionary significant unit that appears to be kept isolated by the Andes to the West and a stable hybrid zone to the South-east. We suggest a complete stop of trade in the Cuyo region and the development of conservation measures, particularly of the cliffs with colonies, which are crucial for the survival this population.
Most alterations in the environment, like climate change, are potential sources of new or intensified directional selection on traits important for the fitness of the species living in it [e.g. ]. Evolutionary responses take place on a time scale comparable to that of changes in climate, but the degree of adaptation will depend on the interplay of natural selection with processes such as gene flow, genetic drift, mutation and demography [e.g. ]. The tight link between genetic variation and climatic variables here reported, in light of the present-day reality of climate change, could lead to important conservation implications in burrowing parrots. Climate change is likely to impose selection pressures on traits important for fitness [e.g. ], thus affecting the chances of persistence of this species. Furthermore, climate change could differently affect the four population clusters detected throughout the species range and, in addition, the populations could vary in the rate of adaptation. Consequently the outcome of climate change, on top of other environmental constraints (e.g. presence of water and Monte vegetation, occurrence of cliffs), could be particularly important for some of the populations, in particular the currently endangered Bloxami and Andinus populations. Under added selection pressure, such as that imposed by ongoing climate change, populations can respond in roughly three ways: 1) by shifting in abundance and distribution, 2) by going extinct, or 3) by evolving . Even if predicting which one, or even which combinations of them, will occur is difficult [92, 93], some likely scenarios can be expected for the burrowing parrot population clusters. Shifts in the distribution could be possible in the case of Bloxami, distributed in Chile. Due to intensive poaching, several cliffs along the historical distribution of the species (Figure 1), which traditional contained colonies, have been found to be empty at present [41, 85] (JFM pers. observ.). Provided that water and natural food items are still available in those places, and that the current small size of Bloxami is not further reduced, these colonies could be reoccupied by burrowing parrots. The situation appears to be quite different for the size-reduced Andinus population in the Cuyo region of Argentina (Additional file 4 Figure S1). In this area, only few suitable breeding places (high cliffs close to water and Monte vegetation) are available, leaving this population only two possibilities in case of strong or too sudden climate change: going extinct or evolving. For Patagonus1 and Patagonus2, two genetically distinct, yet phenotypically indistinguishable populations, the situation appears to differ again. In a few occasions, individuals from several colonies of these populations were found to breed in nests dug in human-build structures like shafts of mines, abandoned adobe buildings and wells for collecting water for cattle . Additionally, some suitable breeding places are available in the southernmost as well in the easternmost areas of the historical distribution of the species. Provided the future occurrence of a more benign climate in the South and a release of human-induced pressure in the East, some of those areas could be used by burrowing parrots. However, the currently rapid habitat degradation in the region inhabited by Patagonus1 and Patagonus2 (see above) makes this possibility very uncertain, as in fragmented landscapes, rapid climate change has the potential to overwhelm the capacity for adaptation of the populations and dramatically alter their genetic composition . Altogether, how closely adaptation can be expected to accommodate climate change and the additional pressure of habitat loss and fragmentation, remains a question for further research.
Faria PJ, Guedes NMR, Yamashita C, Martuscelli P, Miyaki CY: Genetic variation and population structure of the endangered Hyacinth Macaw (Anodorhynchus hyacinthinus): implications for conservation. Biodiversity and Conservation. 2008, 17: 765-779. 10.1007/s10531-007-9312-1.
Joseph L, Dolman G, Donnellan S, Saint KM, Berg ML, Bennett ATD: Where and when does a ring start and end? Testing the ring-species hypothesis in a species complex of Australian parrots. Proc R Soc Lond B. 2008, 275: 2431-2440. 10.1098/rspb.2008.0765.
Ribas CC, Miyaki CY, Cracraft J: Phylogenetic relationships, diversification and biogeography in Neotropical Brotogeris parakeets. J Biogeogr. 2009, 36: 1712-1729. 10.1111/j.1365-2699.2009.02131.x.
Holderegger R, Wagner HH: Landscape genetics. BioScience. 2008, 58: 199-207. 10.1641/B580306.
Hanski I: Metapopulation dynamics. Nature. 1998, 396: 41-49. 10.1038/23876.
Gibbs JP: Demography versus habitat fragmentation as determinants of genetic variation in wild populations. Biol Conserv. 2001, 100: 15-20. 10.1016/S0006-3207(00)00203-2.
Holzhauer SIJ, Ekschmitt K, Sander A-C, Dauber J, Wolters V: Effect of historic landscape change on the genetic structure of the bush-cricket Metrioptera roeseli. Landsc Ecol. 2006, 21: 891-899. 10.1007/s10980-005-0438-9.
Keyghobadi N: The genetic implications of habitat fragmentation for animals. Can J Zool. 2007, 85: 1049-1064. 10.1139/Z07-095.
Koscinski D, Yates AG, Handford P, Lougheed SC: Effects of landscape and history on diversification of a montane, stream-breeding amphibian. J Biogeogr. 2009, 36: 255-265. 10.1111/j.1365-2699.2008.02009.x.
Frankham R: Do island populations have less genetic variation than mainland populations?. Heredity. 1997, 78: 311-327. 10.1038/hdy.1997.46.
Grant PR, Grant BR: Hybridization and bird species. Science. 1992, 256: 193-197. 10.1126/science.256.5054.193.
Bell DA: Genetic differentiation, geographic variation and hybridization in gulls of the Larus glaucescens-occidentalis complex. Condor. 1996, 98: 527-546. 10.2307/1369566.
Parsons J, Olson SL, Braun MJ: Unidirectional spread of secondary sexual plumage traits across an avian hybrid zone. Science. 1993, 260: 1643-1646. 10.1126/science.260.5114.1643.
Grant PR, Grant BR: Phenotypic and genetic effects of hybridization in Darwin's Finches. Evolution. 1994, 48: 297-316. 10.2307/2410094.
Rhymer JM, Williams MJ, Braun MJ: Mitochondrial analysis of gene flow between New Zealand mallards (Anas platyrhynchos) and grey ducks (A. superciliosa). Auk. 1994, 111: 970-978.
Pierotti R, Annett CA: Hybridization and male parental investment in birds. Condor. 1993, 95: 670-679. 10.2307/1369610.
Dowling TE, Secor CL: The role of hybridization and introgression in the diversification of animals. Annu Rev Ecol Syst. 1997, 28: 593-619. 10.1146/annurev.ecolsys.28.1.593.
Avise JC, Walker D: Pleistocene phylogeographic effects on avian populations and the speciation process. Proc R Soc Lond B. 1998, 265: 457-463. 10.1098/rspb.1998.0317.
Kim I, Phillips CJ, Monjeau JA, Birney EC, Noack K, Pumo DE, Sikes RS, Dole JA: Habitat islands, genetic diversity, and gene flow in a Patagonian rodent. Molecular Ecology. 1998, 7: 667-678. 10.1046/j.1365-294x.1998.00369.x.
Martinez JJ, Gonzalez-Ittig RE, Theiler GR, Ojeda R, Lanzone C, Ojeda A, Gardenal CN: Patterns of speciation in two sibling species of Graomys (Rodentia, Cricetidae) based on mtDNA sequences. J Zool Syst Evol Res. 2010, 48: 159-166. 10.1111/j.1439-0469.2009.00539.x.
Ojeda AA: Phylogeography and genetic variation in the South American rodent Tympanoctomys barrerae (Rodentia: Octodontidae). J Mammalogy. 2010, 91: 302-313. 10.1644/09-MAMM-A-177.1.
Morando M, Avila LJ, Baker J, Sites JW: Phylogeny and phylogeography of the Liolaemus darwinii complex (Squamata: Liolaemidae): Evidence for introgression and incomplete lineage sorting. Evolution. 2004, 58: 842-861.
Avila LJ, Morando M, Sites JW: Congeneric phylogeography: hypothesizing species limits and evolutionary processes in Patagonian lizards of the Liolaemus boulengeri group (Squamata: Liolaemini). Biol J Linn Soc. 2006, 89: 241-275. 10.1111/j.1095-8312.2006.00666.x.
Yoke MM, Morando M, Avila LJ, Sites JW: Phylogeography and genetic structure in the Cnemidophorus longicauda complex (Squamata, Teiidae). Herpetologica. 2006, 62: 420-434. 10.1655/0018-0831(2006)62[424:PAGSIT]2.0.CO;2.
Koscinski D, Handford P, Tubaro PL, Sharp S, Lougheed SC: Pleistocene climatic cycling and diversification of the Andean treefrog, Hypsiboas andinus. Mol Ecol. 2008, 17: 2012-2025. 10.1111/j.1365-294X.2008.03733.x.
Alonso Roldan V, Rossi Fraire H, Navarro JL, Gardenal CM, Martella MB: Genetic structure of Greater Rhea (Rhea americana) populations in two regions with different land-uses in central Argentina. Emu. 2009, 109: 214-221. 10.1071/MU08073.
Iriondo MH, Garcia NO: Climatic variations in the Argentine Plains during last 18.000 years. Palaeogeo, Palaeoclim, Palaeoecol. 1993, 101: 209-220. 10.1016/0031-0182(93)90013-9.
Iriondo M: La Pampa. Climas cuaternarios en America del Sur. Edited by: Argollo J, Mourguiart P. 1995, La Paz, Bolivia: Institut Français de Recherche Scientifique pour le Développement en Coopération, 283-306.
Tonni EP, Cione AL: Los mamíferos como indicadores de cambios climáticos en el Cuatemario de la región pampeana de la Argentina. Climas cuaternarios en America del Sur. Edited by: Argollo J, Mourguiart P. 1995, La Paz, Bolivia: Institut Français de Recherche Scientifique pour le Développement en Coopération, 319-326.
Cione AL, Tonni EP, Soibelzon L: The broken zig-zag: Late Cenozoic large mammal and tortoise extinction in South America. Rev Mus Argentino Cienc Nat, n s. 2003, 5: 1-19.
Clapperton CM: The glaciation of the Andes. Quat Scie Rev. 1983, 2: 83-155. 10.1016/0277-3791(83)90005-7.
Vuilleumier BS: Pleistocene changes in the flora and fauna of South America. Science. 1971, 173: 771-780. 10.1126/science.173.3999.771.
Tambussi C, Acosta Hospitaleche C, Horlent N: La avifauna del Cuaternario de Argentina: inferencias paleoambientales a partir del registro de los Psittacidae. Mon Soc Hist Nat Balears. 2007, 14: 69-80.
Darrieu CA: Las razas geográficas de Cyanoliseus patagonus (Aves: Psittacidae). Neotropica. 1980, 26: 207-216.
Badano EI, Cavieres LA, Molina-Montenegro MA, Quiroz CL: Slope aspect influences plant association patterns in the Mediterranean matorral of central Chile. J Arid Environm. 2005, 62: 93-108. 10.1016/j.jaridenv.2004.10.012.
López RP, Larrea Alcázar D, Macía MJ: The arid and dry plant formations of South America and their floristic connections: new data, new interpretation?. Darwiniana. 2006, 44: 18-31.
Masello JF, Pagnossin ML, Sommer C, Quillfeldt P: Population size, provisioning frequency, flock size and foraging range at the largest known colony of Psittaciformes: the Burrowing Parrots of the north-eastern Patagonian coastal cliffs. Emu. 2006, 106: 69-79. 10.1071/MU04047.
Masello JF, Lubjuhn T, Quillfeldt P: Is the structural and psittacofulvin-based colouration of wild Burrowing Parrots Cyanoliseus patagonus condition dependent?. J Avian Biol. 2008, 39: 653-662. 10.1111/j.1600-048X.2008.04417.x.
Nores M, Yzurieta D: Especiación en las Sierras Pampeanas de Córdoba y San Luis (Argentina), con descripción de siete nuevas subespecies de aves. Hornero. 1983, 88-102. N° Extraordinario
Bucher EH, Rinaldi S: Distribución y situación actual del loro barranquero (Cyanoliseus patagonus) en la Argentina. Vida Silv Neotrop. 1986, 1: 55-61.
Rojas Martínez ME: Estudio de la interacción entre las poblaciones de loro tricahue Cyanoliseus patagonus bloxami, y la actividad agrícola en las comunas de Vicuña y Monte Patria, Región de Coquimbo, Chile. 2008, Santiago de Chile: Servicio Agrícola y Ganadero, Ministerio de Agricultura, Gobierno de Chile
Bó NA: Notas preliminares sobre la avifauna del nordeste de San Luis. Hornero. 1965, 10: 251-268.
Klauke N, Masello JF, Quillfeldt P, Segelbacher G: Isolation of tetranucleotide microsatellite loci in the burrowing parrot (Cyanoliseus patagonus). Journal of Ornithology. 2009, 150: 921-924. 10.1007/s10336-009-0423-1.
Failla M, Seijas VA, Quillfeldt P, Masello JF: Potencial impacto del loro barranquero (Cyanoliseus patagonus): evaluación de percepción de daño en Patagonia Nordeste, Argentina. Gestión Ambiental. 2008, 16: 27-40.
Pezzola A, Winschel C, Sánchez R: Estudio multitemporal de la degradación del monte nativo en el partido de Patagones - Buenos Aires. 2004, INTA-Ascasubi, Ascasubi, Buenos Aires, Argentina
Masello JF, Marchesan M, Quillfeldt P: Zehn Jahre Forschung in der größten Papageienkolonie der Welt - Teil 1. Die Felsensittiche im Nordosten Patagoniens. Papageien. 2008, H.12/2008: 426-429.
Masello JF, Pagnossin GA, Palleiro GE, Quillfeldt P: Use of miniature security cameras to record behaviour of burrow-nesting birds. Vogelwarte. 2001, 41: 150-154.
Masello JF, Sramkova A, Quillfeldt P, Epplen JT, Lubjuhn T: Genetic monogamy in burrowing parrots Cyanoliseus patagonus?. J Avian Biol. 2002, 33: 99-103. 10.1034/j.1600-048X.2002.330116.x.
Masello JF, Pagnossin ML, Lubjuhn T, Quillfeldt P: Ornamental non-carotenoid red feathers of wild Burrowing Parrots. Ecol Res. 2004, 19: 421-432. 10.1111/j.1440-1703.2004.00653.x.
Masello JF, Choconi RG, Sehgal RMN, Tell LA, Quillfeldt P: Blood and intestinal parasites in wild Psittaciformes: a case study of Burrowing Parrots (Cyanoliseus patagonus). Orn Neotropical. 2006, 17: 515-529.
Masello JF, Choconi RG, Helmer M, Kremberg T, Lubjuhn T, Quillfeldt P: Do leucocytes reflect condition in nestling burrowing parrots (Cyanoliseus patagonus) in the wild?. Comp Biochem Physiol A. 2009, 152: 176-181. 10.1016/j.cbpa.2008.09.018.
Masello JF, Lubjuhn T, Quillfeldt P: Hidden dichromatism in Burrowing Parrots Cyanoliseus patagonus as revealed by spectrometric colour analysis. Hornero. 2009, 24: 47-55.
Masello JF, Quillfeldt P: Chick growth and breeding success of the Burrowing Parrot. Condor. 2002, 104: 574-586. 10.1650/0010-5422(2002)104[0574:CGABSO]2.0.CO;2.
Masello JF, Quillfeldt P: Body size, body condition and ornamental feathers of Burrowing Parrots: Variation between years and sexes, assortative mating and influences on breeding success. Emu. 2003, 103: 149-161. 10.1071/MU02036.
Masello JF, Quillfeldt P: Are haematological parameters related to body condition, ornamentation and breeding success in wild burrowing parrots Cyanoliseus patagonus?. J Avian Biol. 2004, 35: 445-454. 10.1111/j.0908-8857.2004.03278.x.
Masello JF, Quillfeldt P: Consequences of La Niña phase of ENSO for the survival and growth of nestling Burrowing Parrots on the Atlantic coast of South America. Emu. 2004, 104: 337-346. 10.1071/MU03039.
Masello JF, Quillfeldt P: Klimawandel und Brutverhalten: erfolgreich brüten in wechselhafter Umwelt? Eine Fallstudie am Felsensittich. Vogelwarte. 2008, 46: 302-303.
Bucher EH, Bertin MA, Santamaria AB: Reproduction and molt in the burrowing parrot. Wilson Bull. 1987, 99: 107-109.
Astuti D, Azuma N, Suzuki H, Higashi S: Phylogenetic relationships within parrots (Psittacidae) inferred from mitochondrial cytochrome-b gene sequences. Zool Sci. 2006, 23: 191-198. 10.2108/zsj.23.191.
Palumbi SR: Nucleic acids II: The polymerase chain reaction. Molecular Systematics. 1996, Sunderland, Massachusetts: Sinauer Associates, Inc, 2
Eberhard JR, Bermingham E: Phylogeny and biogeography of the Amazona ochrocephala (Aves: Psittacidae) complex. Auk. 2004, 121: 318-332. 10.1642/0004-8038(2004)121[0318:PABOTA]2.0.CO;2.
Librado P, Rozas J: DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25: 1451-1452. 10.1093/bioinformatics/btp187.
Corander J, Tang J: Bayesian analysis of population structure based on linked molecular information. Mathematical Biosciences. 2007, 205: 19-31. 10.1016/j.mbs.2006.09.015.
Jobb G, von Haeseler A, Strimmer K: TREEFINDER: a powerful graphical analysis environment for molecular phylogenetics. BMC Evolutionary Biology. 2004, 4: e18-10.1186/1471-2148-4-18.
Jobb G: TREEFINDER. (Version of October 2008). 2008. Munich. Distributed by the author, [http://www.treefinder.de]
Bandelt HJ, Forster P, Röhl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48.
Rogers AR, Harpending H: Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol. 1992, 9: 552-569.
Excoffier L, Lischer HEL: Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Res. 2010, 10: 564-567. 10.1111/j.1755-0998.2010.02847.x.
Tonni EP, Laza JH: Las aves de la Fauna Local Paso de Otero (Pleistoceno Tardio) en la Provincia de Buenos Aires. Su significación ecológica, climática y zoogeográfica. Ameghiniana. 1980, 17: 313-322.
Tambussi CP, Noriega JI, Tonni EP: Late Cenozoic birds of Buenos Aires Province (Argentina): an attempt to document quantitative faunal changes. Palaeogeo, Palaeoclim, Palaeoecol. 1993, 101: 117-129. 10.1016/0031-0182(93)90155-C.
Tambussi C, Acosta Hospitaleche C, Rinderknecht A, Ubilla M: Parrots (Aves, Psittaciformes) in the Pleistocene of Uruguay. Ameghiniana. 2009, 46: 431-435.
Acosta Hospitaleche CI, Tambussi CP: Cyanoliseus patagonopsis nov. sp. (Aves, Psittaciformes) del Pleistoceno de Punta Hermengo, provincia de Buenos Aires. Cyanoliseus patagonopsis nov. sp. (Aves, Psittaciformes) of the Pleistocene of Punta Hermengo, Buenos Aires Province. Ameghiniana. 2006, 43: 249-253.
Tonni EP: Cyanoliseus ensenadensis (Cattoi) nov. comb. (Aves, Psittaciformes) del Pleistoceno Medio de la provincia de Buenos Aires. Ameghiniana. 1972, 2: 145-148.
Weir JT, Schluter D: Calibrating the avian molecular clock. Molecular Ecology. 2008, 17: 2321-2328. 10.1111/j.1365-294X.2008.03742.x.
Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. Bmc Evolutionary Biology. 2007, 7: 214-10.1186/1471-2148-7-214.
Anderson MJ: DISTLM v.5: a FORTRAN computer program to calculate a distance-based multivariate analysis for a linear model. 2004, (v.5). Auckland, New Zealand, Department of Statistics, University of Auckland, [http://www.stat.auckland.ac.nz/~mja/Programs.htm#Mine]
Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A: Very high resolution interpolated climate surfaces for global land areas. Int J Climatol. 2005, 25: 1965-1978. 10.1002/joc.1276.
Anderson MJ: DISTLM forward: a FORTRAN computer program to calculate a distance-based multivariate analysis for a linear model using forward selection. 2003, Auckland, New Zealand, Department of Statistics, University of Auckland, [http://www.stat.auckland.ac.nz/~mja/Programs.htm#Mine]
François O, Ancelet S, Guillot G: Bayesian clustering using hidden Markov random fields in spatial population genetics. Genetics. 2006, 174: 805-816. 10.1534/genetics.106.059923.
McCulloch RD, Bentley MJ, Purves RS, Hulton NRJ, Sugden DE, Clapperton CM: Climatic inferences from glacial and palaeoecological evidence at the last glacial termination, southern South America. J Quat Scie. 2000, 15: 409-417. 10.1002/1099-1417(200005)15:4<409::AID-JQS539>3.0.CO;2-#.
Elías A: Nuevo registro del pelícano pardo (Pelecanus thagus) en la Argentina. Nuestras Aves. 2008, 53: 48-
Suriano JM, Ferro G, Dalponte M: Convenio Dirección de Minería de Río Negro-Segemar: Geología y Recursos Minerales de la Hoja 3966-IV, Choele Choel y del Sector Rionegrino de las Hojas 3963-III, Colonia Juliá; 4163-I-III, General Conesa y 4163-II-IV, Viedma. 1999, Viedma
Bruno G, Torres R: Distribución, status, caracterización del hábitat, y diagnóstico de uso, del loro barranquero (Cyanoliseus patagonus) en el centro de Argentina. 2006, Córdoba, Proyecto Elé
Glade A, (ed): Libro rojo de los vertebrados terrestres de Chile. 1993, Santiago de Chile: CONAF, 2
Naranjo Navia J: Análisis de las variaciones poblacionales en la especie Cyanoliseus patagonus bloxami, (loro tricahue) en la VII Región de Chile. 2005, Dinámica de Poblaciones, Universidad de Chile, Santiago de Chile, Report, 20 pp
Bucher EH: Neotropical parrots as agricultural pests. New world parrots in crisis: solutions from conservation biology. Edited by: Beissinger SR, Snyder NFR. 1992, New York and London: Smithsonian Institution Press, 201-219.
Voitzuk L: El loro barranquero, plaga agrícola del sud de la provincia de Buenos Aires. Bol Fitosanitario. 1975, 48: 31-34.
Moschione FN, Banchs RA: Proyecto Calas. Una experiencia de manejo adaptativo para el aprovechamiento sustentable de psitácidos y como estrategia de conservación de sus habitats en la Argentina. Manejo de Fauna Silvestre en la Argentina. Programas de uso sustentable. Edited by: Bolkovic ML, Ramadori D. 2006, Buenos Aires: Dirección de Fauna Silvestre, 27-37.
Balmford A, Green RE, Jenkins M: Measuring the changing state of nature. Trends Ecol Evol. 2003, 18: 326-330. 10.1016/S0169-5347(03)00067-3.
Moschione F, González MA: Distribución del Loro Barranquero Cyanoliseus patagonus andinus en la Provincia de Salta. Libro de resúmenes. XI Reunión Argentina de Ornitología. 2005, Aves Argentinas, Buenos Aires, 10-
Gienapp P, Teplitsky C, Alho JS, Mills JA, Merilä J: Climate change and evolution: disentangling environmental and genetic responses. Molecular Ecology. 2008, 17: 167-178. 10.1111/j.1365-294X.2007.03413.x.
Davis MB, Shaw RG, Etterson JR: Evolutionary responses to changing climate. Ecology. 2005, 86: 1704-1714. 10.1890/03-0788.
Holt RD: The Microevolutionary Consequences of Climate Change. Trends in Ecology & Evolution. 1990, 5: 311-315. 10.1016/0169-5347(90)90088-U.
Masello JF, Quillfeldt P: La colonia de loros barranqueros en la costa rionegrina de El Cóndor. Un patrimonio mundial. Las mesetas patagónicas que caen al mar: la costa rionegrina. Edited by: Masera RF, Lew J, Serra Peirano G. 2005, Viedma: Ministerio de Familia, Gobierno de Río Negro, Argentina, 349-371.
Jump AS, Peñuelas J: Running to stand still: adaptation and the response of plants to rapid climate change. Ecology Letters. 2005, 8: 1010-1020. 10.1111/j.1461-0248.2005.00796.x.
We wish to thank the following persons for helping with sample collection: Orlando Amaya, Brent Barrett, Julio Bufelli, Mirta Carbajal, Gabriel Cortéz, Francesca Cunninghame, Miguel Escobar, Esteban Fernández, Víctor Fratto, Pablo Giovine, Patricia González, Anna Grasse, Guillermo Luna Jorquera, Mariano Linares, Fabián Llanos, Walter Marcial, Gabriela Murga Velasco, Enrique Narvaes, Juanjo Navarro, Sandra Rivera, Oscar Saá, Sergio Sánchez, Lorenzo Sympson, Adalberto Taux, Ricardo Torres, and Luis Valenzuela Riquelme. We also would like to thank the following curators for their help: Friederike Woog and Iris Heynen (Staatliches Museum für Naturkunde Stuttgart), Jean-Louis Martelli (Musée des Confluences, Centre de Conservation et d'Etude des Collections, Lyon), Sylke Frahnert (Museum für Naturkunde Berlin), Marie-Dominique Wandhammer and Dominique Nitka (Musée zoologique de Strasbourg), Alexander Haas and Jakob Hallermann (Zoologisches Museum, Universität Hamburg), Bernd Nicolai (Museum Heineanum, Halberstadt) and Dietrich Heidecke (Zoologisches Sammlungen Halle, Martin-Luther-Universität Halle Wittenberg). This project was partially supported by the World Parrot Trust, the German Ornithologists' Society (Deutsche Ornithologen-Gesellschaft), the Liz Claiborne Art Ortenberg Foundation, the Wildlife Conservation Society, and the Deutsche Forschungsgemeinschaft (QU148-1, Germany). We wish to thank Ann Michels for helping with the application for the necessary permits in Chile. We especially thank Jamie Gilardi, Graham Harris and Bill Conway for crucial support. In Argentina, the present study was carried out under permission of the Dirección de Fauna de la Provincia de Río Negro, Argentina (Exp. no. 143089-DF-98), and the Dirección de Fauna Silvestre, Secretaría de Ambiente de la Nación, Argentina (CITES permits).
The authors declare that they have no competing interests.
JFM and PQ conceived and designed the study. YM, GS and MS planned the genetic analyses. JFM, PQ, MF and MC carried out the extensive fieldwork. GKM, NK and JFM generated the molecular data. JFM, NK and YM participated in bioinformatic analyses. YM, JFM and PQ were responsible for data analysis and drafted the manuscript. GS and MS reviewed the final draft of the manuscript. All authors read and approved the final manuscript.