The impact of distance and a shifting temperature gradient on genetic connectivity across a heterogeneous landscape
© Rossetto et al; licensee BioMed Central Ltd. 2011
Received: 16 February 2011
Accepted: 18 May 2011
Published: 18 May 2011
Inter-population distance and differences in breeding times are barriers to reproduction that can contribute to genotypic differentiation between populations. Temporal changes in environmental conditions and local selective processes can further contribute to the establishment of reproductive barriers. Telopea speciosissima (Proteaceae) is an excellent subject for studying the effect of geographic, edaphic and phenological heterogeneity on genotypic differentiation because previous studies show that these factors are correlated with morphological variation. Molecular, morphological and environmental datasets were combined to characterise the relative influence of these factors on inter-population differentiation, and Bayesian analyses were used to investigate current levels of admixture between differentiated genomes.
A landscape genetic approach involving molecular and morphological analyses identified three endpoints of differentiated population groups: coastal, upland and southern. The southern populations, isolated from the other populations by an edaphic barrier, show low migration and no evidence of admixture with other populations. Amongst the northern populations, coastal and upland populations are connected along a skewed altitudinal gradient by genetically intermediate populations. The strong association between temperature and flowering time in Telopea speciosissima was shown to maintain a temporally unstable reproductive barrier between coastal and upland populations.
Substrate-mediated allopatry appears to be responsible for long-term genetic isolation of the southern populations. However, the temperature-dependent reproductive barrier between upland and coastal populations bears the genetic signature of temporal adjustments. The extreme climatic events of the last glacial maximum are likely to have caused more complete allochronic isolation between upland and coastal populations, as well as exerting increased selective pressure upon local genomes. However, at intermediate altitudes, current climatic conditions allow for the incorporation of alleles from previously distinct genomes, generating new, intermediate genomic assemblages and possibly increasing overall adaptive potential.
The identification of reproductive boundaries between genetically differentiated populations can provide useful cues to the factors influencing population-level connectivity and micro-evolutionary processes. Geographic and/or habitat isolation are among the most important factors responsible for establishing the levels of reproductive segregation that lead to measurable genetic differentiation . As dispersal of individuals and gametes is generally more likely among geographically close demes, landscape discontinuities and distributional gaps result in increased neutral divergence . Allochronic separation (divergence in breeding times) can also result in non-random mating and contribute to increased genetic variance across the landscape . A range of environmental factors (particularly temperature and photoperiod) can affect flowering phenology and consequently cause temporal reproductive isolation .
Local selective processes can further contribute to the development of pre- or post-zygotic barriers to reproduction. In such circumstances gene flow is restricted by assortative mating caused by the reduced competitiveness of migrants and/or admixed individuals within differentiated habitats [5, 6]. Often, more than one factor affects reproductive isolation between populations and to further complicate matters the relative importance of historical processes needs to be taken into account. As a result, the combination of historical and geographical analyses is increasingly used to differentiate between past and present landscape-level connectivity .
Temporal changes in local environmental conditions, such as those experienced during the climatic cycles of the Quaternary, can cause localised bottlenecks and extinctions of populations that were previously part of a continuous distributional range . Such periodical contractions can lead to temporal allopatry, drift and vicariant differentiation. As environmental conditions improve landscape-level connectivity can be re-established. These contraction/expansion cycles are influenced by the species' ecology and the landscape features that characterise their habitat, and the genetic structure measured across current-day distributions often reflects the interactions between historical biology and regional geography [9, 10].
Unravelling the sequence and tempo of the events establishing reproductive barriers and affecting genetic structure can provide interesting insights into micro- and macro-evolutionary processes, as well as support the development of improved conservation strategies that are mindful of evolutionary potential. This is particularly relevant in a climate change context where predicted climatic shifts could have a significant impact on inter-population connectivity and micro-evolutionary potential.
Telopea speciosissima R.Br. (the Waratah, Proteaceae) is a particularly good model for studying the relative influence of climatic and environmental conditions on between-population connectivity: it is distributed along altitudinal and latitudinal gradients; it has a limited flowering season and marked differences in flowering times between populations ; and its distribution is interrupted by edaphic barriers. Furthermore, an analysis of morphometric variation in T. speciosissima sensu lato  showed that what had been considered to be a single, geographically widespread species actually comprised two allopatric, edaphically differentiated taxa: T. speciosissima on sandstones of the Sydney Basin and T. aspera on the northern granites of the Gibraltar Range (New South Wales, Australia). The major axis of variation from ordinations separated these two species as distinct clusters but secondary ordination axes revealed considerable residual variation, suggesting that differentiation was also present among southern T. speciosissima populations.
In this study we aim to detect the genetic signatures that identify the geographic, edaphic and climatic factors leading to differentiation between Waratah populations. In particular, we ask the following questions: is there congruence between morphological and genetic differentiation; how do environmental variables (geographic, edaphic, climatic) differentially impact on genetic structure; and can we use species-wide analyses to identify directionality of gene flow among differentiated population groups, and investigate the temporal strength of reproductive barriers?
Individual plants have long generation times as primary regeneration is through pyrogenesis from large lignotubers. Waratahs flower in a short season, generally lasting six weeks in spring, early in low-lying warmer coastal areas and later in higher cooler areas. The highly conspicuous red inflorescences are primarily bird pollinated . Flowers are protandrous, pollen shedding onto the pollen presenter (modified style tip) several days before receptivity, and have been shown to be exclusively out-crossing [14, 15]. Seed production peaks several years after fire regeneration . Seeds have no dormancy mechanisms and germinate readily given adequate moisture conditions .
Sampling, DNA extractions and PCRs
We aimed to obtain a suitable account of the genetic diversity across the entire distribution of the species, as well as obtain a sample that would overlap with that of previous morphological studies . Leaf material was collected and total genomic DNA was extracted using DNeasy® 96 plant kits (QIAGEN®, Hilden, Germany). Seven simple sequence repeat (SSR) loci specifically developed for T. speciosissima  were used in this study. Six loci (TS03bgt, TS04bgt, TS12bgt, TS18bgt, TS23bgt, TS27bgt) were dimers (CT repeats) while TS16bgt included a dimer (CT) and a tetramer (CTCA). PCR conditions followed those reported by Porter et al. . To test genotyping accuracy PCRs were repeated for over 20% of the reactions. Less than 5% of those repeats identified errors that needed confirmation through replication of PCR and genotyping.
Allelic distribution was measured, and to avoid potential bias caused by uneven sampling , a standardized estimate of allelic richness independent of sample size  was calculated using the program FSTAT 2.9.3 . Measures of expected (He) and observed (Ho) heterozygosities, and the inbreeding coefficient (FIS; ) were also calculated using FSTAT 2.9.3, with the probability of FIS being greater than zero determined after 10,000 permutations, and 95% confidence interval determined after 15,000 bootstraps. Hardy-Weinberg and linkage equilibriums were assessed using the exact test with the program GENEPOP 3.2a , and significance levels were determined after 1,000 batches of 5,000 iterations each. Sequential Bonferroni corrections were applied to significance tests .
Shannon's information index of population subdivision, SHUA, was used to obtain an estimate of between-population dispersal based on an average effective population size of 100 individuals (N100m). SHUA provides a robust estimation of genetic exchange over a wider range of dispersal values than FST, because it is more sensitive to rare alleles, it can incorporate unequal numbers of individuals per population and shows no dependence between measures at different hierarchical levels [24, 25].
We used the Bayesian clustering method described by Pritchard et al.  and implemented in STUCTURE 2.3.3 to identify (in the absence of preliminary information on population boundaries) genetically differentiated groups of individuals. After a preliminary test aimed at finding a suitable range for K and the optimal burn-in period, we tested K from 1 to 18 through five independent runs. All runs were based on 7.5 × 105 Markov chain Monte Carlo iterations after a burn-in period of 3 × 105 iterations without prior information on the locality of origin of the individuals sampled. The admixture frequency model was run under the assumption of correlated allele frequencies to improve clustering of closely related populations . The optimal number of clusters was verified using the ΔK statistical approach suggested by Evanno et al. .
Analysis of molecular variance (AMOVA; ) was used to quantify variance components and the significance of the genetic subdivisions identified by the Bayesian test, as well as a range of other relevant groupings. Finally, principal coordinate analysis (PCoA) was used to produce a comparative graphical representation of genotypic similarity based on genetic distances among populations. Finally, to reveal the possible influence of isolation by distance (IBD) in explaining between-population genetic differentiation, Mantel tests (999 permutations) were performed between matrices of linearised pairwise FST values (FST/(1- FST) and
ln-adjusted pairwise geographical distances. AMOVA, PCoA, Mantel and SHUA were run using GenAlEx6.3 .
We used the Bayesian model implemented in NewHybrids  to detect genetic intermediate individuals resulting from the interbreeding between genetically distinct populations, and distinguish from similarity caused by incomplete lineage sorting. The original tests and simulations by Anderson & Thompson  show that admixture can be detected without the need for diagnostic alleles. Although a high number of informative loci produces much better posterior probabilities (PP) for assigning to hybrid categories, when FST >0.2 (as is the case in this study) a smaller number of loci is sufficient for detecting admixture events .
Our objective was not to identify and quantify specific admixture categories but rather to ascertain the mixing of differentiated upland and coastal genomes. To do this we conducted a series of preliminary analyses that culminated in a run using seven populations: the three with highest PP of assignment to the upland group, the three with highest PP of assignment to the coastal group, and the four at intermediate altitude. Jeffreys priors were used with a burn-in of 3.5 × 105 sweeps followed by 7.5 × 105 sweeps. Posterior probability of assignment as pure, F1s, F2s, and backcrosses were initially measured and than combined to obtain proportions of admixture vs. pure.
Morphology and phenology
We selected morphometric variables that had high loadings (>|0.95|) on the second axis of the PCoA of Crisp & Weston  and scored them for of the available herbarium specimens collected from localities included in our genetic sampling. All variables were measurements or counts of leaf attributes and all were represented by the mean of the measurements or counts from five leaves. We used the same method of ordination used to ordinate our genetic data (PCoA, but using squared Euclidean distances instead of genetic distance). One variable, leaf width, was used to represent size variation. This was also transformed to produce a size-independent shape variable by taking the natural logarithm of the ratio of leaf length to leaf width. Two other variables were similarly transformed.
The date of first flowering of individual plants from six populations across the altitudinal range of T. speciosissima (70, 180, 200, 410, 760, 980 m) was recorded during previous site visits (n= 256). Flowering dates were converted to day of the year and correlated to altitude of collection site. Data on flowering time of mixed seedlings grown in three plantations at different altitudes in the Sydney region (280, 580, 780 m) were extracted from Dupree & Goodwin  and similarly compared (n = 731). Finally, flowering date of plants collected from wild sites and grown in a common garden environment (at Mount Annan Botanic Garden, Sydney) was recorded and compared with altitude of provenance (180, 400, 500, 760 m; n= 85).
The climatic and environmental data presented in Figures 1 and 2 were created using ArcView V3.3 and ArcGIS V9.2. Temperature and rainfall data are from WorldClim Version 1.4 (release 3) average for 1960-1990. Data are 30 arc-second resolution grid (approximately 0.7 km2 - equivalent to 930 m vertical × 750-775 m horizontal resolution).
Genetic diversity in T. speciosissima
Population genetic statistics for T. speciosissima
33 13 S
151 12 E
33 24 S
150 13 E
33 27 S
150 39 E
Bells Line of Road
33 30 S
150 16 E
33 31 S
150 37 E
33 32 S
151 17 E
33 37 S
151 17 E
33 39 S
151 11 E
33 04 S
151 20 E
33 46 S
150 23 E
34 27 S
150 27 E
34 38 S
150 39 E
Waterfall Flat (Royal NP)
34 09 S
151 00 E
Bottle Forest (Royal NP)
34 05 S
151 00 E
Curra Moors (Royal NP)
34 09 S
151 05 E
35 22 S
150 28 E
35 04 S
150 25 E
35 08 S
150 41 E
Genetic structure in T. speciosissima
Genetic structure and gene flow in T. speciosissima
Average pair-wise values
T. speciosissima mean
Upland vs. rest
Coastal vs. rest
Southern vs. rest
Southern vs. rest
Upland and Coastal vs. all others
Upland vs. Coastal vs. Southern
Coastal vs. Intermediate
Upland vs. Coastal vs. Southern vs. all others
Upland vs. Intermediate
Upland vs. Southern vs. all others
Connectivity and admixture
A recent study by Roberts et al.  showed that relaxing assignment threshold to pure from 95% PP to 80% PP made little difference to their findings. As a result, here we define as pure individuals, those with >80% PP of being assigned to a particular population group. In the three selected coastal populations 89% of individuals were assigned to pure coastal. This was in contrast with the three selected upland populations where 51% of individuals were assigned to pure upland. Interestingly, the Newnes population showed no evidence of admixture in contrast to the geographically close Bells Line of Road population in which 35% of individuals had a PP of >80% of being admixed. The situation for the intermediate populations was more complex: a single individual at Mountain Lagoon was identified as pure upland, while 52% of individuals were assigned to pure coastal. At the opposite end of the spectrum all individuals were assigned to pure coastal at Carrington Falls.
Morphology and phenology
PCoA of the morphometric data produced a similar ordination to the PCoA of genetic data in that specimens from upland, southern and coastal groups of populations formed clusters that occupied the same relative positions to the corresponding clusters in the genetic analysis, albeit with some overlap as would be expected from characters that are both phenotypically labile and ontogenetically variable (Figure 3b; ).
First flowering in five wild locations occurred when the 10 day mean maximum was 20°C and the minimum was 9°C (Australian Bureau of Meteorology), regardless of altitude. The timing of the temperature shift necessary for flowering varied across the altitudinal gradient, and corresponded to the genetic differentiation observed between coastal and upland populations (Figures 1 and 2).
The analysis of flowering time from 256 flowering events from seven natural locations distributed along the altitudinal gradient showed a very strong association between altitude and flowering time (R2 = 0.929; P < 0.0001). For instance, average flowering time at Evans Lookout (980 m) was on Julian day 312, while at the coastal site at Patonga (180 m) it was on day 261. There was a weaker but significant relationship between flowering time and altitude in three cultivated populations sourced from mixed genotypes (R2 = 0.471; P < 0.0001), and data from the common garden environment suggested that altitude of original collection (sourced from three different altitudes) had no influence on flowering within a controlled environment (R2 = 0.011; P = 0.174).
Genetic diversity and structure
Levels of genetic diversity were generally consistent across all sampled Waratah populations, irrespective of their size and location. Heterozygosity is likely to be maintained in small populations through preferential outcrossing (T. speciosissima is highly receptive to cross pollen; [37, 38]), and measured values were similar to those recorded for other geographically localised, bird pollinated species of Proteaceae (eg. Banksia hookeriana; ).
Although high gene flow maintains equilibrium at a local scale, between-population connectivity appears to be less consistent in T. speciosissima. PCoA and Bayesian assignment identified three major population groups (southern, upland, coastal) supported by AMOVA and between-group FST values. Morphological analysis discovered similar patterns of geographic variation, confirming strong differentiation in the south and identifying a previously unreported upland morphotype.
Understanding the processes leading to genetic and morphological differentiation requires the identification of the barriers that cause partial or complete reproductive isolation . Edaphic conditions and geographic distance are likely to represent important isolating mechanisms for the Waratah. Shale-derived soils of the Cumberland Plain form an island of unsuitable habitat at the heart of the overall distributional range of the species (Figure 2d), and the southern populations are located at the edge of the sandstones to which the species is restricted. Furthermore, the southern group grows on Permian sediments (rather than on Triassic rocks) and is isolated by a barrier of volcanic-derived soils. The combination of geographic distance and marginal habitat could therefore explain the considerable genetic and morphological differentiation detected in the south.
Although the terminally winged seeds of T. speciosissima are not easily dispersed and tend to germinate in proximity to the maternal plant , it is unlikely that limited dispersal is the only mechanism restricting gene flow. The Waratah is bird and mammal pollinated  and our data show that gene flow across large distances is possible (within the coastal group for example), thus suggesting that factors other than dispersal impact on between-population connectivity.
Altitude, temperature, phenology and local genomes
Our study detected considerable genetic structure along the altitudinal gradient and a close association between altitude, temperature and flowering time in T. speciosissima. Flowering in the Waratah is triggered by a shift in average maximum temperature to around 20°C, and this shift occurs gradually along the altitudinal gradient (Figures 1, 2). Coastal populations flower early (20°C average attained in September) and upland populations flower late (20°C average attained in November). The temperature gradient associated with topography impacts on phenology and forms a temporal isolating barrier between coastal and upland populations.
If temporal reproductive isolation is the predominant factor separating coastal and upland genomes we would expect it to be incomplete across transitional zones where flowering overlap is likely. The Waratah populations that are intermediate in altitude and flowering time are also genetically intermediate (Figures 3, 4), suggesting that they operate as a contact zone between the upland and coastal genomes (as also supported by the estimated number of migrants; Table 2).
In contrast, there was very low migration and no evidence of admixture between the southern and any of the other populations, suggesting that the temporal barriers to gene flow along the altitudinal gradient are easier to overcome than those constrained by geographic distance and edaphic factors. While geographic isolation can be difficult to overcome without suitable means of dispersal , allochronic isolation can be quickly bypassed following changes in climatic conditions that enable phenological overlap.
Current vs. historical isolation
Admixture can take place without the need for complete sympatry, particularly when pollen exchange is mediated by highly mobile vectors . As a result, unless the Waratah populations at intermediate altitude have established relatively recently a more even allelic spread should be expected along the East/West distributional axis, especially given the level of migration between the intermediate and the other two population groups. Two main factors could lead to the strong genetic structure observed: long-term fluctuations in climatic conditions that affect the strength of phenological isolation, and adaptation to local conditions that lead to habitat isolation.
If a strong association between flowering time and local conditions exists, a change in these conditions will directly impact on reproductive isolation . Our common-garden phenological data emphasize the predominance of a genetically controlled response to shifting temperature conditions along the altitudinal gradient. As a result, the temperature differential that separates the upland from the coastal genomes represents an adjustable reproductive barrier. In current climatic conditions, geographic proximity and partial phenological overlap allow for genetic exchange to occur and result in the observed admixture front.
However, during glacial maxima temperature gradients were more intense, and in some upland areas temperatures were 10°C lower than present . Hesse et al.  suggested that during the LGM the Newnes Plateau (where some of the upland Waratah populations are located) experienced extreme climatic conditions. The activity of aeolian sand dunes between 30 and 15Kya suggests that during the LGM (and previous glacial peaks) the area was dominated by grasslands . As a result, surviving T. speciosissima populations would have experienced significantly cooler and drier conditions than those experienced by the populations living at lower altitudes, and allochronic isolation would have been stronger or complete. In such climatically distinct areas, assortative mating would have been a powerful factor in developing and maintaining the genomic differences that were measured in this study.
Exposure to differentiated environmental conditions could also have established selective reproductive barriers and reduced the viability of migrant and/or admixed genotypes. The fact that the Newnes population is not admixed (unlike other upland populations) suggests that vertebrate-mediated pollen dispersal is not uniform across what is fairly continuous habitat. Such pattern suggests that habitat isolation might also contribute to the genetic differentiation observed between upland and coastal populations, a hypothesis we are currently testing via cross-pollination, reciprocal transplant and controlled-environment experiments.
Our study on a long-lived woody shrub distributed along altitudinal and latitudinal gradients, provides evidence for the contrasting temporal stability of different reproductive barriers. While the association between geographic and genetic isolation appears to be stable, the temperature-dependent reproductive barrier between upland and coastal population bears the genetic signature of temporal adjustments. The detection of localised admixture events and the availability of morphological, phenological and climatic (current and historical) information helped us identify a transitional zone established as a result of the post-glacial alleviation of extreme differences in climatic conditions.
These temporal changes resulted in the merger of reproductively isolated and differentiated genomes, a process that can have deleterious (sterility, homogenisation, outbreeding depression) or beneficial (the incorporation of novel diversity favouring adaptation to changing conditions) evolutionary consequences . The long-term consequences of this genomic reshuffling in the Waratah are yet to be clarified but provide a valuable framework for in-depth adaptation studies. Understanding how changes in connectivity driven by climatic fluctuations affect natural populations has important implications for the management of biodiversity in view of predicted environmental change, and will be the focus of long-term in-situ experimentation and genomic studies on the Waratah.
This project was funded by the Hermon Slade Foundation (HSF 07-1). The authors would like to thank M. Baum, D. Benson, M. Crisp and P. Rymer for useful discussions.
- Coyne JA, Orr HA: Speciation. 2004, Sunderland, Massachusetts: Sinauer AssociatesGoogle Scholar
- Templeton AR: Population genetics and microevolutionary theory. Book Population genetics and microevolutionary theory. 2006, City: Wiley Online Library, (Editor ed.^eds.)View ArticleGoogle Scholar
- Devaux C, Lande R: Incipient allochronic speciation due to non-selective assortative mating by flowering time, mutation and genetic drift. P Roy Soc B. 2008, 275: 2723-2732. 10.1098/rspb.2008.0882.View ArticleGoogle Scholar
- Elzinga JA, Atlan A, Biere A, Gigord L, Weis AE, Bernasconi G: Time after time: flowering phenology and biotic interactions. Trends Ecol Evol. 2007, 22: 432-439. 10.1016/j.tree.2007.05.006.View ArticlePubMedGoogle Scholar
- Nosil P, Vines TH, Funk DJ: Reproductive isolation caused by natural selection against immigrants from divergent habitats. Evolution. 2005, 59: 705-719.PubMedGoogle Scholar
- Lowry DB, Modliszewski JL, Wright KM, Wu CA, Willis JH: The strength and genetic basis of reproductive isolating barriers in flowering plants. Philos T Roy Soc B. 2008, 363: 3009-3021. 10.1098/rstb.2008.0064.View ArticleGoogle Scholar
- Sork VL, Waits L: Contributions of landscape genetics-approaches, insights, and future potential. Mol Ecol. 2010, 19: 3489-3495. 10.1111/j.1365-294X.2010.04786.x.View ArticlePubMedGoogle Scholar
- Willis KJ, Niklas KJ: The role of Quaternary environmental change in plant macroevolution: the exception or the rule?. Phil Trans B. 2004, 359: 159-172. 10.1098/rstb.2003.1387.View ArticleGoogle Scholar
- Hewitt G: The genetic legacy of the Quaternary ice ages. Nature. 2000, 405: 907-913. 10.1038/35016000.View ArticlePubMedGoogle Scholar
- Petit RJ, Aguinagalde I, de Beaulieu JL, Bittkau C, Brewer S, Cheddadi R, Ennos R, Fineschi S, Grivet D, Lascoux M: Glacial Refugia: hotspots but not melting pots of genetic diversity. Science. 2003, 300 (5625): 1563-1565. 10.1126/science.1083264.View ArticlePubMedGoogle Scholar
- Offord CA: Floral- and leaf- character variation in waratah (Telopea speciosissima, family Proteaceae), with reference to conservation and horticultural improvement. Aust J Bot. 2007, 55: 684-691. 10.1071/BT06211.View ArticleGoogle Scholar
- Crisp MD, Weston PH: Geographic and ontogenetic variation in morphology of Australian Waratahs (Telopea: Proteaceae). Syst Biol. 1993, 42: 49-76.View ArticleGoogle Scholar
- Pyke GH: Pollination biology of Telopea speciosissima. Waratahs-their biology, cultivation and conservation. Edited by: Armstrong JA. 1987, Canberra, Australia: Australian National Botanic Gardens Occasional Publication, 9: 35-43.Google Scholar
- Goldingay RL, Whelan RJ: The influence of pollinators on fruit positioning in the Australian shrub Telopea speciosissima (Proteaceae). Oikos. 1993, 68: 501-509. 10.2307/3544918.View ArticleGoogle Scholar
- Offord CA: An examination of the reproductive biology of Telopea speciosissima (Proteaceae) with emphasis on the nature of protandry and the role of self-pollination in fruit set. Int J Plant Sci. 2004, 165: 73-83. 10.1086/380745.View ArticleGoogle Scholar
- Denham AJ, Auld TD: Flowering, seed dispersal, seed predation and seedling recruitment in two pyrogenic flowering resprouters. Aust J Bot. 2002, 50: 545-557. 10.1071/BT02009.View ArticleGoogle Scholar
- Porter C, Rymer PD, Rossetto M: Isolation and characterization of microsatellite markers for the waratah, Telopea speciosissima (Proteaceae). Mol Ecol N. 2006, 6: 446-448. 10.1111/j.1471-8286.2006.01266.x.View ArticleGoogle Scholar
- Leberg PL: Estimating allelic richness: effects of sample size and bottlenecks. Mol Ecol. 2002, 11: 2445-2449.View ArticlePubMedGoogle Scholar
- El Mousadik A, Petit RJ: High level of genetic differentiation for allelic richness among populations of the argan tree [Argania spinosa (L.) Skeels] endemic to Morocco. Theor Appl Genet. 1996, 92: 832-839. 10.1007/BF00221895.View ArticlePubMedGoogle Scholar
- Goudet J: FSTAT (version 1.2): a computer program to calculate F-statistics. J Hered. 1995, 86: 485-486.Google Scholar
- Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.View ArticleGoogle Scholar
- Raymond M, Rousset F: GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Hered. 1995, 86: 248-Google Scholar
- Rice WR: Analyzing tables of statistical tests. Evolution. 1989, 43: 223-225. 10.2307/2409177.View ArticleGoogle Scholar
- Sherwin WB, Jabot F, Rush R, Rossetto M: Measurement of biological information with applications from genes to landscapes. Mol Ecol. 2006, 15: 2857-2869. 10.1111/j.1365-294X.2006.02992.x.View ArticlePubMedGoogle Scholar
- Rossetto M, Kooyman R, Sherwin W, Jones R: Dispersal limitations, rather than bottlenecks or habitat specificity, can restrict the distribution of rare and endemic rainforest trees. Am J Bot. 2008, 95: 321-329. 10.3732/ajb.95.3.321.View ArticlePubMedGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMedPubMed CentralGoogle Scholar
- Falush D, Stephens M, Pritchard JK: Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003, 164: 1567-1587.PubMedPubMed CentralGoogle Scholar
- Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.View ArticlePubMedGoogle Scholar
- Excoffier L, Smouse PE, Quattro JM: Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992, 131: 479-PubMedPubMed CentralGoogle Scholar
- Peakall ROD, Smouse PE: GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol N. 2006, 6: 288-295. 10.1111/j.1471-8286.2005.01155.x.View ArticleGoogle Scholar
- Anderson EC, Thompson EA: A model-based method for identifying species hybrids using multilocus genetic data. Genetics. 2002, 160: 1217-1229.PubMedPubMed CentralGoogle Scholar
- Vähä JP, Primmer CR: Efficiency of model-based Bayesian methods for detecting hybrid individuals under different hybridization scenarios and with different numbers of loci. Mol Ecol. 2006, 15: 63-72.View ArticlePubMedGoogle Scholar
- Dupree SA, Goodwin PB: Effect of temperature, daylength and growth regulators on flowering of Protea, Telopea and Leucospermum. Acta Hort. 1990, 264: 79-86.View ArticleGoogle Scholar
- Thompson S, Lamothe M, Meirmans PG, Perinet P, Isabel N: Repeated unidirectional introgression towards Populus balsamifera in contact zones of exotic and native poplars. Mol Ecol. 2010, 19: 132-145.View ArticlePubMedGoogle Scholar
- Quintela M, Thulin CG, Höglund J: Detecting hybridization between willow grouse (Lagopus lagopus) and rock ptarmigan (L. muta) in Central Sweden through Bayesian admixture analyses and mtDNA screening. Cons Gen. 2010, 11: 557-569. 10.1007/s10592-009-0040-9.View ArticleGoogle Scholar
- Roberts DG, Gray CA, West RJ, Ayre DJ: Marine genetic swamping: hybrids replace an obligately estuarine fish. Mol Ecol. 2010, 19: 508-520. 10.1111/j.1365-294X.2009.04501.x.View ArticlePubMedGoogle Scholar
- Whelan RJ, Goldingay RL: Factors affecting fruit-set in Telopea speciosissima (Proteaceae): the importance of pollen limitation. J Ecol. 1989, 1123-1134.Google Scholar
- Goldingay RL: Further assessment of pollen limitation in the waratah (Telopea speciosissima). Aust J Bot. 2000, 48: 209-214. 10.1071/BT97110.View ArticleGoogle Scholar
- He T, Lamont BB: High microsatellite genetic diversity fails to predict greater population resistance to extreme drought. Cons Genet. 2010, 11: 1445-1451. 10.1007/s10592-009-9971-4.View ArticleGoogle Scholar
- Rossetto M, Crayn D, Ford A, Mellick R, Sommerville K: The influence of environment and life-history traits on the distribution of genes and individuals: a comparative study of 11 rainforest trees. Mol Ecol. 2009, 18: 1422-1438. 10.1111/j.1365-294X.2009.04111.x.View ArticlePubMedGoogle Scholar
- Hendry AP, Day T: Population structure attributable to reproductive time: isolation by time and adaptation by time. Mol Ecol. 2005, 14: 901-916. 10.1111/j.1365-294X.2005.02480.x.View ArticlePubMedGoogle Scholar
- Barrows TT, Stone JO, Fifield LK, Cresswell RG: Late Pleistocene glaciation of the Kosciuszko Massif, Snowy Mountains, Australia. Quaternary Res. 2001, 55: 179-189. 10.1006/qres.2001.2216.View ArticleGoogle Scholar
- Hesse PP, Humphreys GS, Selkirk PM, Adamson DA, Gore DB, Nobes DC, Price DM, Schwenninger JL, Smith B, Tulau M: Late Quaternary aeolian dunes on the presently humid Blue Mountains, Eastern Australia. Quatern Int. 2003, 108: 13-32. 10.1016/S1040-6182(02)00191-X.View ArticleGoogle Scholar
- Arnold ML, Martin NH: Hybrid fitness across time and habitats. Trends Ecol Evol. 2010, 25: 530-536. 10.1016/j.tree.2010.06.005.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.