Pervasive admixture between eucalypt species has consequences for conservation and assisted migration

Abstract Conservation management often uses information on genetic population structure to assess the importance of local provenancing for ecological restoration and reintroduction programs. For species that do not exhibit complete reproductive isolation, the estimation of population genetic parameters may be influenced by the extent of admixture. Therefore, to avoid perverse outcomes for conservation, genetically informed management strategies must determine whether hybridization between species is relevant, and the extent to which observed population genetic patterns are shaped by interspecific versus intraspecific gene flow. We used genotyping by sequencing to identify over 2,400 informative single nucleotide polymorphisms across 18 populations of Eucalyptus regnans F. Muell., a foundation tree species of montane forests in south‐eastern Australia. We used these data to determine the extent of hybridization with another species, Eucalyptus obliqua L'Hér., and investigate how admixture influences genetic diversity parameters, by estimating metrics of genetic diversity and examining population genetic structure in datasets with and without admixed individuals. We found hybrid individuals at all sites and two highly introgressed populations. Hybrid individuals were not distributed evenly across environmental gradients, with logistic regression identifying hybrids as being associated with temperature. Removal of hybrids resulted in increases in genetic differentiation (F ST), expected heterozygosity, observed heterozygosity and the inbreeding coefficient, and different patterns of isolation by distance. After removal of hybrids and introgressed populations, mountain ash showed very little population genetic structure, with a small effect of isolation by distance, and very low global F ST(0.03). Our study shows that, in plants, decisions around provenancing of individuals for restoration depend on knowledge of whether hybridization is influencing population genetic structure. For species in which most genetic variation is held within populations, there may be little benefit in planning conservation strategies around environmental adaptation of seed sources. The possibility for adaptive introgression may also be relevant when species regularly hybridize.

formed management strategies must determine whether hybridization between species is relevant, and the extent to which observed population genetic patterns are shaped by interspecific versus intraspecific gene flow. We used genotyping by sequencing to identify over 2,400 informative single nucleotide polymorphisms across 18 populations of Eucalyptus regnans F. Muell., a foundation tree species of montane forests in south-eastern Australia. We used these data to determine the extent of hybridization with another species, Eucalyptus obliqua L'Hér., and investigate how admixture influences genetic diversity parameters, by estimating metrics of genetic diversity and examining population genetic structure in datasets with and without admixed individuals. We found hybrid individuals at all sites and two highly introgressed populations. Hybrid individuals were not distributed evenly across environmental gradients, with logistic regression identifying hybrids as being associated with temperature. Removal of hybrids resulted in increases in genetic differentiation (F ST ), expected heterozygosity, observed heterozygosity and the inbreeding coefficient, and different patterns of isolation by distance. After removal of hybrids and introgressed populations, mountain ash showed very little population genetic structure, with a small effect of isolation by distance, and very low global F ST (0.03). Our study shows that, in plants, decisions around provenancing of individuals for restoration depend on knowledge of whether hybridization is influencing population genetic structure. For species in which most genetic variation is held within populations, there may be little benefit in planning conservation strategies around environmental adaptation of seed sources. The possibility for adaptive introgression may also be relevant when species regularly hybridize.

| INTRODUC TI ON
Substantial biodiversity declines are occurring in many regions of the world due to widespread land clearing, habitat degradation, introduced species and climate change (Evans et al., 2011;Pounds et al., 2006;Woinarski, Burbidge, & Harrison, 2015). Extensive and ongoing land clearing has led to major reductions in forest cover globally (Achard et al., 2014;Reside et al., 2017;Taubert et al., 2018), with synergistic interactions between stressors placing some ecosystems under high threat of rapid collapse or changes in ecosystem state (Brook, Sodhi, & Bradshaw, 2008;Lindenmayer, Hobbs, Likens, Krebs, & Banks, 2011;Lindenmayer & Sato, 2018). With such widespread changes facing ecosystems, it is critical to understand how these stressors interact with the fundamental ecological processes operating within and between foundation species, to adequately manage biodiversity across landscapes.
Using genetic approaches to inform management activities allows conservation efforts to be targeted towards sites of unique genetic composition or adaptive importance, making population genetic studies valuable in many taxa (Ikeda et al., 2017;Maunder, Cowan, Stranc, & Fay, 2001;McCartney-Melstad & Shaffer, 2015;Reynolds et al., 2015). To maximize beneficial outcomes, it is vital that our understanding of the population genetic diversity and structure in target species is as accurate as possible, particularly in the context of a changing environment. For example, understanding patterns of local adaptation across the range of a species is important for developing methods of assisted gene flow to mitigate the impacts of climate change and other threatening processes (Kelly & Phillips, 2016;Supple et al., 2018).
Genetically informed conservation requires a detailed understanding of the spatial distribution of genetic diversity, particularly as it relates to environmental adaptation. Spatial genetic structure and population genetic differentiation are typically considered to be driven by the influences of gene flow, genetic drift and local adaptation (Orsini, Vanoverbeke, Swillen, Mergeay, & Meester, 2013).
However, for species that do not exist in complete reproductive isolation, the estimation of population genetic parameters may be influenced by the extent of hybridization and introgression with closely related species. This could have large implications for the application of genetic data to conservation management of species, for example, by committing resources to putatively distinct populations, when they may actually contain highly admixed individuals.
With the advent of modern DNA genotyping techniques, studies investigating thousands of genetic markers from across the genome are becoming more common (Gaughran et al., 2018;Hand et al., 2015;Harvey, Aleixo, Ribas, & Brumfield, 2017;Hudson, Freeman, Myburg, Potts, & Vaillancourt, 2015), and several studies have investigated patterns of nuclear genetic structure and gene flow across large geographic regions (Hecht, Matala, Hess, & Narum, 2015;Hendricks et al., 2017;Sampson et al., 2018;Shriver et al., 2005). Studies such as these provide critical information for the conservation of populations with unique genetic heritage, identification of areas of adaptive potential for assisted migration and location of source populations or historical refugia (Hecht et al., 2015;Supple et al., 2018). len and seed dispersal, with pollen typically playing a greater role in eucalypts because it tends to disperse further than seeds (Barber, 1965;Petit et al., 2005;Potts & Wiltshire, 1997). Comparisons of the maternally inherited chloroplast and biparentally inherited nuclear DNA have shown that pollen-mediated gene flow can be up to at least 200 times greater than seed-mediated gene flow in some species (Bloomfield, Nevill, Potts, Vaillancourt, & Steane, 2011;Nevill, Bradbury, Williams, Tomlinson, & Krauss, 2014), although at least one study found that gene flow from seed dispersal is practically equivalent to that from pollen dispersal (Jones, Shepherd, Henry, & Delves, 2006).
Given the importance of understanding population genetic structure for conservation, and the knowledge that hybridization in eucalypts is a widespread and common phenomenon (Griffin, Burgess, & Wolf, 1988), we investigated these two aspects in Eucalyptus regnans F. Muell., (mountain ash) one of Australia's most well-known and economically important trees. The existence of hybrids between E. regnans and the frequently co-occurring Eucalyptus obliqua L'Hér. (messmate stringybark) has long been known (Ashton, 1956); however, the extent of hybridization across the range of the species has never been investigated. Similarly, while the chloroplast genetic structure of E. regnans has been studied (Nevill, Bossinger, & Ades, 2010), the structure of the nuclear genome has not. To address these knowledge gaps, our aims were to (a) identify the extent and possible drivers of hybridization across the geographic distribution of E. regnans, (b) describe how identification of admixture using genomic data may influence our understanding of population genetic structure and (c) consider how these factors would influence current management strategies in eucalypts. We address these aims using genotyping by sequencing to obtain large numbers of genomewide genetic markers on individual samples across the natural geographic K E Y W O R D S admixture, eucalypt, Eucalyptus regnans, gene flow, hybrid, mountain ash, population genetics, SNP range of the species. We predict that (a) some individuals and populations will show greater levels of admixture with E. obliqua, (b) levels of admixture will be driven in part by local environmental variables and (c) the inclusion or exclusion of hybrid individuals in population genetic analyses will lead to different strategic outcomes for management. If these predictions are true, there are implications for future studies of population genetic structure and the planning of restoration plantings and assisted gene flow.

| Study area and species
Eucalyptus regnans grows in wet forests of the south-east Australian states of Victoria and Tasmania. It is the tallest angiosperm in the world, with reliable records of individuals exceeding 100 m (Beale, 2007;Hardy, 1918Hardy, , 1935. It is also a serotinous obligate seeder, requiring high-intensity fires to open the understorey, create fertile ash beds and stimulate the mass release of seeds from the forest canopy (Ashton, 1981b;Ashton & Chinner, 1999). Without fire, trees are typically unable to produce offspring that survive to maturity, primarily due to predation of seeds by ants (Ashton, 1979;O'Dowd & Gill, 1984), low availability of light (Gilbert, 1959), browsing of seedlings by herbivores and fungal infection of seedlings (Ashton & Macauley, 1972).
Eucalyptus regnans is patchily distributed through a 700 km by 500 km area, growing only where climatic conditions are suitable (Cochrane, 1969). It reaches its highest elevations (>1,100 m ASL) in the northernmost part of its range, on the Errinundra Plateau, and grows near to sea level in some southern parts of its Tasmanian distribution. As the island of Tasmania has been separated from the Australian mainland by more than 200 km for over ten thousand years (Duncan, Worth, Jordan, Jones, & Vaillancourt, 2016;Lambeck, Rouby, Purcell, Sun, & Sambridge, 2014), it is assumed that there has been very little or no gene flow between E. regnans stands in these two regions for at least that length of time.
In a number of locations throughout Victoria and Tasmania, trees displaying intermediate characteristics between E. regnans and other species have been recorded. These specimens have been identified mostly as hybrid individuals between E. regnans and E. obliqua, and, less commonly, E. regnans and E. macrorhyncha (Ashton, 1958(Ashton, , 1981aAshton & Sandiford, 1988). At least two individuals have also been found that appear to be tri-hybrids-the result of a E. regnansx obliqua hybrid individual mating with a E. macrorhyncha (Yorke & Ashton, 1982). As red stringybark does not occur naturally in Tasmania, E. regnans × macrorhyncha hybrids do not occur there.

| Sample collection
We collected 387 E. regnans tissue samples from across its geographic distribution (Figure 1). At each of 16 sites, we walked a transect collecting tissue from trees spaced at least 20 m apart, until we had sampled 20 trees. We targeted trees with diameters at breast height of more than 60 cm, to avoid sampling younger trees that were propagated after the practice of reseeding logged coupes using seed of nonlocal provenance became common practice (Flint & Fagg, 2007). As the combined effects of logging and wildfires have caused a reduction in the size and frequency of old undisturbed patches of trees (Lindenmayer, Blanchard, Blair, McBurney, & Banks, 2016), a linear transect of fixed length was sometimes impossible.
At seven of the sites, we collected a second sample from the 20th tree, to serve as a technical replicate from the field. We were also able to incorporate an extra 42 E. regnans samples collected during fieldwork for other studies into some analyses, taken from various locations (Supporting Information Table S1). Twenty-one E. obliqua samples, taken from the Cathedral Range region in Victoria, were also sequenced to allow us to determine the extent of hybridization between the two species. All samples were putatively identified as E. regnans or E. obliqua using purported diagnostic morphological characters (Brooker & Kleinig, 2006).
Tissue collected was mostly cambium, obtained by cutting through the rough and smooth bark using a machete, and slicing off a 10 × 5 × 0.2 cm strip. A small number of samples were leaf tissue, collected by climbing trees using standard arborist techniques and picking two fresh, growing leaves. All samples were air-dried in individual brown paper bags and then stored at 4°C prior to DNA extraction.

| Sample preparation and sequencing
Approximately 600 mg of each tissue sample was chilled to −65°C and homogenized using an automated tissue grinding machine (Labman Max Planck Cryogenic Grinder Dispenser, Labman Automation & Custom Robotics), before storage at −18°C to await DNA extraction.
Samples were ordered randomly, and whole genomic DNA was

| Demultiplexing and initial filters
Of the 408 samples (387 E. regnans and 21 E. obliqua), sequencing resulted in nearly 1.49 billion read pairs. We demultiplexed reads using exact matches and combinatorial index mode with Axe (Murray & Borevitz, 2018) and were unable to assign 7% of read pairs to a sample. We then used BBDuk to remove adapters and quality-trim (Phred score Q = 30) reads at both ends, and NextGenMap (Sedlazeck, Rescheneder, & von Haeseler, 2013) to align reads to the E. grandis v2.0 reference genome (Bartholome et al., 2015;Myburg et al., 2014). We used SAMtools (Li et al., 2009) to convert the dataset into sample-specific Binary Alignment/Map (BAM) files and sort reads.
To create a sample-by-SNP matrix, we used the "ANGSD" soft- Two separate filtering strategies were conducted on this dataset using the statistical software package R (R Core Team, 2017). The first of these developed a set of SNPs for investigating the extent of admixture with E. obliqua. The second method was used to investigate whether hybridization influences population genetic structure and isolation by distance across the geographic distribution of E. regnans. Each of these filters has been discussed in the methods of the relevant analysis. formed a PCA using the indpca function of the "hierfstat" (Goudet, 2005) package in R. The first two principal components for each sample were plotted using the "ggplot2" (Wickham, 2009) package, with samples coloured by the mean level of admixture with E. obliqua, as determined above.

| Extent of admixture
Next, we calculated mean admixture proportions for every population and mapped this across the landscape using the "ggplot2" and "scatterpie" (Yu, 2018) packages.

| Influence of admixture on population structure
To determine the level of influence that unrealized hybridization can have on population structure, we first filtered SNPs using call rate (≥0.4), minor allele frequency (MAF ≥ 0.01) and observed heterozygosity (≤0.5). The likelihood that each SNP does not deviate from Hardy-Weinberg equilibrium (HWE) was checked using the HWChisqStats function of the "Hardy-Weinberg" (Graffelman, 2015) package. Any SNP out of HWE in more than three populations (where n ≥ 15) was removed from further analysis. In addition, the snpgdsLDpruning function in the "SNPRelate" (Zheng et al., 2012) package was used to prune out SNPs using a linkage disequilibrium threshold of 0.5 and a sliding window of 5,000 bp.
Filtering was done on two groups: (a) all samples including hybrids (but excluding the reference E. obliqua samples) and (b) pure E. regnans individuals (i.e., those with <10% admixture with E. obliqua). In the admixture-inclusive group, 2,474 SNPs were retained.
For the pure group, 2,481 SNPs were retained. Any samples with more than 50% missing data were removed from each dataset, leaving 323 samples and 228 samples in the admixture-inclusive and admixture-free groups, respectively. Visual inspection of hierarchical clustering dendrograms, created using the hclust function on a Euclidean genetic distance matrix, showed that all biological replicates were closely paired in both datasets (Supporting Information Figure S2), indicating they were reliable and contamination was not likely to be a factor.
Standard population genetic parameters were calculated for all loci in both groups, using the basic.stats function in the "hierfstat" (Goudet, 2005)  were created using the cmdscale function, with populations plotted using the "ggplot2" package. To investigate isolation by distance, we calculated pairwise F ST between populations using the "Nei87" method of the genet.dist function in the "hierfstat" package, and regressed these distances against pairwise population geographic distance, calculated using the earth.dist function in the "fossil" (Vavrek, 2011) package. Genetic differentiation and geographic distance were transformed to F ST /(1 − F ST ) and log(geographic distance), respectively, to allow for linear interpretation (Rousset, 1997).

| Spatial structure and population genetics in Eucalyptus regnans
Spatial structuring of genotypes was investigated using two methods; firstly, a Mantel test (Mantel, 1967) comparing pairwise genetic distance with the natural logarithm of pairwise geographic distance was performed using the mantel.rtest function in the R package "ade4" (Dray & Dufour, 2007). We then used the spatial structure analysis function in GenAlEx v6.503 (Peakall & Smouse, 2006 to determine the maximum geographic distance at which genotypes show significant spatial autocorrelation. This was interpreted as the maximum distance at which the lower confidence interval of the spatial autocorrelation coefficient r was greater than zero on the y-axis.
To investigate population genetics within E. regnans, we calcu- We then calculated the pairwise genetic distance between all individuals in these populations using the "Dch" method of the genet.
dist function. This distance matrix was then read into GenAlEx to conduct an analysis of molecular variance (AMOVA) with 999 permutations to determine the amount of genetic variation explained within and among sites.

| Extent of admixture
Despite NGSadmix using a different approach (likelihoods instead of called genotypes) and much larger number of loci, there was a very strong correlation of ancestry coefficients computed using that method and both STRUCTURE (r = 0.94, p < 0.001) and SNPRelate (0.95, p < 0.001), with bar plots showing near-identical patterns between methods (Supporting Information Figure   S1). There was an even stronger correlation between the re- had coefficients between 0.6 and 0.9, and 16 had coefficients of greater than 0.9.
Principal components analysis showed that highly admixed samples (e.g., those from Wilsons Promontory) were able to be easily identified due to their separation from the bulk of the

| Environmental association with admixture
Model selection showed that hybrid individuals were not randomly distributed across all environmental variables. Of the 11 top models (ΔAICc < 2), MINTEMP and MAXTEMP were identified in all 11, suggesting that these two variables had the strongest effects on the probability of hybrid occurrence. Sites with a high MAXTEMP (hot summers) and sites with a low MINTEMP (cold winters) had lower probabilities of hybrid occurrence, whereas sites with a high MINTEMP and sites with a low MAXTEMP had higher probabilities of hybrid occurrence (Figure 4). Variables that occurred in fewer top models included NORTH, RAD, PHOS, RAIN, AGE and TWI; however, none of these variables showed any clear trends with hybrid occurrence (Supporting Information Figure S3). EAST was not identified in any top models.

| Influence of admixture on population structure
The removal of hybrid individuals from the dataset modified the values of the genetic parameters investigated ( Figure 5)

| Spatial structure and population genetics in Eucalyptus regnans
After removing all hybrids and discounting sites that had <10 sam- the Errinundra Plateau site as being relatively distinct from the rest of the sites.

| D ISCUSS I ON
Using a genomewide SNP dataset, we investigated the extent of hybridization between two widespread forest trees and the way in which admixture influences estimation of genetic parameters and interpretation of population genetic structure. We found that hybrids were not distributed evenly across geographic or environmental space, with some populations showing more admixture than others, and a strong association of hybrid occurrence with summer temperatures. Once hybrids were removed, we found very little evidence of population genetic differentiation or local isolation of alleles, with high levels of gene flow, slow generation times and a lack of historical geographic isolation probably responsible for these findings.
This study has implications for our understanding of the pro- and rapid adaptation to novel conditions. It is in this context that we discuss our findings and explore their implications for studies of population genetics and genetically informed conservation.

| Hybridization
Had we not considered the possibility of hybrid individuals occurring in our dataset, we would have included introgressed sites such as Wilsons Promontory and Bruny Island in the analyses, as well as 41 hybrid samples scattered through our other sites. Many of the hybrid individuals showed no obvious morphological differences (e.g., in bud and fruit shape or rough bark height and thickness) to pure specimens, with morphological approaches now generally regarded as a poor way of identifying hybrids (Field et al., 2009;McKinnon, Smith, & Potts, 2010;Rhymer & Simberloff, 1996;Schwabe, Neale, & McGlaughlin, 2015). Despite this, the Wilsons Promontory samples, which are morphologically more similar to E. regnans but genetically much closer to E. obliqua, had a different growth habit to typical E. regnans individuals, appearing stunted in form. This was initially thought to be the result of phenotypic plasticity rather than genetic architecture, as the trees do not resemble typical E. obliqua either. In our study, removing clear outliers using PCA or a similar method would have retained more than half of the hybrid samples. Eucalypts are typically preferentially outcrossing, open-pollinated, and often found in sympatry with multiple congeneric species, which may partly explain why more than half of all species form natural hybrid combinations (Griffin et al., 1988;Potts, Barbour, Hingston, & Vaillancourt, 2003). In addition, many of these combinations can occur at relatively high frequencies within populations (Field et al., 2009;McKinnon et al., 2010) However, despite the knowledge that hybridization is not uncommon, it is often not considered in population genetic studies of eucalypts. We also note the possibility that for some species there may be multiple hybrid combinations with other species, which may vary regionally (Griffin et al., 1988). In this study, we only considered admixture with a single species; however, gene flow with red stringybark may regularly occur in particular sites within Victoria (Ashton & Sandiford, 1988).
In the case of E. regnans, hybridization with E. obliqua appears to be a more pervasive phenomenon than previously realized, with all sampled sites containing at least one hybrid individual, and two sites where more than half of the samples were hybrids.
The Wilsons Promontory individuals contained very high levels of E. obliqua ancestry, despite their greater morphological similarity to E. regnans. To understand why variation in the level of admixture between populations occurs, examination of the factors controlling gene flow between eucalypt species is necessary. Previous research has identified three predominant drivers of hybridization in eucalypts, including the extent of geographic isolation, the degree of overlap in flowering times, and the level of phylogenetic divergence between species (Barbour, Potts, Vaillancourt, & Tibbits, 2006;Butcher, McDonald, & Bell, 2009;Field et al., 2009;McKinnon et al., 2010;Potts et al., 2003), although strong abiotic (e.g., climatic or geological) gradients may also be explanatory in some cases (Pryor, 1976). When stands of one species are geographically isolated from conspecifics, there is a greater opportunity for pollen from another species to successfully pollinate flowers in the isolated stand, termed "pollen swamping" (Ellstrand & Elam, 1993;Field et al., 2009). Additionally, environmental conditions can influence the degree of overlap in flowering times, with temperature previously shown to be an important driver of flowering phenology in eucalypts (Hudson, Kim, & Keatley, 2010;Law, Mackowski, Schoer, & Tweedie, 2000;Rawal, Kasel, Keatley, & Nitschke, 2015). As E. obliqua is a predominantly summer-flow-

| Patterns of genetic structure
The low levels of genetic population structure that we observed are not atypical for eucalypts, with a number of studies (Broadhurst, Mellick, Knerr, Li, & Supple, 2018;Dillon et al., 2015;Gauli, Steane, Vaillancourt, & Potts, 2014;Sampson et al., 2018;Supple et al., 2018) finding that geographic structure does not contribute greatly to population differentiation. Similarly, the F ST values that we observed between sites are low but comparable to those found in other eucalypts (Sampson et al., 2018;Supple et al., 2018;Yeoh, Bell, Foley, Wallis, & Moran, 2012), and low geographic structuring of genetic diversity is not unusual in widespread forest trees, such as Pinus taeda L. (Eckert et al., 2010) and Quercus lobataNée (Grivet, Sork, Westfall, & Davis, 2008;Sork et al., 2010).
As gene flow in trees is often effected predominantly through pollen dispersal (Sork, 2016), there is often a clear difference in population structure between the nuclear genome (which is inherited biparentally) and genetic components that are inherited maternally, for example, the chloroplast in angiosperms (Sampson et al., 2018). In eucalypts, seed dispersal is highly restricted, with individual trees distributing seeds tens of metres, and stand edges typically only capable of moving about 1-2 m/year (Booth, 2017).
While it is difficult to ascertain the upper limit of the dispersal curve, pollen appears to regularly disperse hundreds of metres to kilometres (Bloomfield et al., 2011;Broadhurst, 2013;Byrne, Elliott, Yates, & Coates, 2008;Sampson et al., 2018). This explains why, when investigating chloroplast structure for phylogeographic purposes, Nevill et al. (2010) found a highly structured genetic pattern in E. regnans, whereas our (nuclear-based) results show very little population structuring. Chloroplast DNA is highly conserved, with genetic structure typically reflecting historical patterns of dispersal and colonization. Further, only a very small proportion of chloroplast variation is typically contained within coding regions (Young, Lanzatella, Sarath, & Tobias, 2011), suggesting that there is unlikely to be substantial levels of local selection acting on the chloroplast.    (Kuparinen, Savolainen, & Schurr, 2010). Previous studies have shown that when selection is very high, high levels of dispersal can maximize local adaptation, and when selection pressure is low (but present), local adaptation is highest under moderate levels of dispersal (Banks, Davies, & Cary, 2017;Forester, Jones, Joost, Landguth, & Lasky, 2016). This is due to migration providing the genetic diversity for selection to act upon, but also potentially overriding the effects of selective processes when selection is weak.

| Implications for seed provenancing
Seed used for native vegetation restoration activities has historically been collected from small local geographic areas (Broadhurst et al., 2008), because of the perceived risk of introducing genotypes that are not adapted to local conditions (Hamilton, 2001). Over the past two decades, this "local is best" approach has been criticized for a number of reasons, including the failure to consider changing environmental/climatic conditions and not incorporating enough genetic variability (Broadhurst et al., 2008;Choi, 2007 (Fitzpatrick, Gerberich, Kronenberger, Angeloni, & Funk, 2015), but also that the potential benefits of gene flow are large and often outweigh the risk of negative impacts associated with anthropogenic disturbances or novel selective pressures (Fisher, Garner, & Walker, 2009;Harrisson et al., 2016). In many cases, the preservation of genetic uniqueness and taxonomic integrity is no longer considered scientifically justifiable (Ralls et al., 2018). Further, as rare new variants are unlikely to be the cause of beneficial adaptation to local conditions (Alberto et al., 2013;Savolainen, Lascoux, & Merilä, 2013), adaptive alleles are still likely to be present in seed collected from distant localities, simply at different frequencies. Reductions in the frequency of some adaptive alleles are typically of minor importance and can be naturally corrected over a small number of generations (Fitzpatrick et al., 2015;Ralls et al., 2018). Thus, we suggest that sourcing seed from trees across a range of environments (local or distant) will ensure adaptive potential for restoration into a changing and challenging future.
There are other considerations for seed sourcing that also need to be considered. For example, climate-adjusted provenancing has been suggested as a way of assisting plant populations with adaptation to climate change (Prober et al., 2015). To achieve this, seed sourcing for restoration would have to consider the future climate at a site and predict the best source populations to match that climate. However, there is one important consideration for this approach that has not been adequately considered in the recent literature: the issue of hybridization and introgression. The possibility that some populations contain a higher degree of admixture with another species cannot be ignored, because interspecific gene flow is one way in which species may cope with rapidly changing abiotic environments (De La Torre, Wang et al., 2014). By selecting seed from populations that have "adapted" to different climatic conditions, managers may actually be selecting seed that is simply admixed with another species.
For example, our data indicate that the occurrence of E. regnans × obliqua hybrids is correlated with temperatures of the hottest and coolest months, so collection of seed based on temperature may inadvertently result in the collection of hybrid genotypes, or even near-pure E. obliqua, which is likely not the desired outcome of management actions. It should also be noted that while our sampling avoided young trees, contemporary seed crops may show different levels of admixture in some populations. For example, those in landscapes where logging and fire have modified the extent and age structure of E. regnans compared to E. obliqua, or where changes in temperature have increased overlap in flowering times. To be certain of the hybrid status at a proposed seed collection locality, genotyping of seedlings would be necessary.

| CON CLUS ION
Our population genomic analysis of E. regnans found widespread admixture of varying levels with a congener, suggesting regular hybridization throughout the range of the species. As many genera of plants are known to form natural hybrids, it is critical that admixture and its role in the adaptive process is considered appropriately in population genetic studies, as introgressed populations may skew genetic analyses and affect management decisions. The combination of widespread hybridization and high levels of gene flow in E. regnans, with similar results having been found for a number of other eucalypt species, suggests that introgressive adaptation through porous genomes may be a common way for this taxon to adapt to rapid environmental change in climate and fire regimes. Selection on hybrids expressing traits harvested from sympatric congeners may allow for rapid adaptive change to new conditions. Furthermore, as the occurrence of hybrid individuals was not distributed evenly across geographic or climatic space, the use of climatic variables to select genotypes for assisted migration may not be the most appropriate way to manage eucalypts for conservation purposes and requires more detailed consideration. Taylor and Oliver King for their assistance and company during fieldwork, and Wade Blanchard for statistical advice. We also appreciate the helpful and constructive comments provided by three anonymous reviewers.

CO N FLI C T O F I NTE R E S T
None declared.

DATA A RCH I V I N G
Sequence data for all individuals are deposited to the NCBI