# Within-range translocations and their consequences in European larch

Stefanie Wagner, Sascha Liepelt, Sophie Gerber, Rémy J. Petit

► To cite this version:

Stefanie Wagner, Sascha Liepelt, Sophie Gerber, Rémy J. Petit. Within-range translocations and their consequences in European larch. PLoS ONE, 2015, 10 (5), 17 p. 10.1371/journal.pone.0127516. [https://doi.org/10.1371/journal.pone.0127516](https://doi.org/10.1371/journal.pone.0127516)

HAL Id: hal-02641244

https://hal.inrae.fr/hal-02641244v1

Submitted on 28 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
RESEARCH ARTICLE

# Within-Range Translocations and Their Consequences in European Larch

Stefanie Wagner¹,², Sascha Liepelt³, Sophie Gerber¹,², Rémy J. Petit¹,²*

1 INRA, UMR1202 BIOGECO, F-33610, Cestas, France, 2 Univ. Bordeaux, BIOGECO, UMR 1202, F-33400, talence, France, 3 University of Marburg, Faculty of Biology, Conservation Biology, D-35032, Marburg, Germany

* petit@pierroton.inra.fr

<figure><img src="image_2.png" /></figure>

CrossMark
click for updates

## Abstract

In contrast to biological invasions, translocations of individuals within a species range are understudied, due to difficulties in systematically detecting them. This results in limited knowledge about the corresponding processes and uncertainties regarding the status of ex- tant populations. European larch, a forest tree whose fragmented native distribution is re- stricted to the Alps and to other Central European mountains, has been massively planted for at least 300 years. Here we focus on the genetic characterization of translocations hav- ing taken place within its native range. Microsatellite variation at 13 nuclear loci and se- quence data of two mitochondrial DNA fragments were analyzed on the basis of a comprehensive range-wide population sample. Two complementary methods (GENECLASS and STRUCTURE) were used to infer translocation events based on nuclear data whereas mi- tochondrial data were used for validation of these inferences. Using GENECLASS, we found translocation events in a majority of populations. Additional cases of translocation and many instances of admixture were identified using STRUCTURE, thanks to the clear-cut ances- tral genetic structure detected in this species. In particular, a strong divide between Alpine and Central European populations, also apparent at mitochondrial markers, helped uncover details on translocation events and related processes. Translocations and associated ad- mixture events were found to be heterogeneously distributed across the species range, with a particularly high frequency in Central Europe. Furthermore, translocations frequently in- volved multiple geographic sources, some of which were over-represented. Our study illus- trates the importance of range-wide investigations for tracing translocations back to their origins and for revealing some of their consequences. It provides some first clues for devel- oping suitable conservation and management strategies.

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Wagner S, Liepelt S, Gerber S, Petit RJ (2015) Within-Range Translocations and Their Consequences in European Larch. PLoS ONE 10(5): e0127516. doi:10.1371/journal.pone.0127516

**Academic Editor:** João Pinto, Instituto de Higiene e Medicina Tropical, PORTUGAL

**Received:** September 30, 2014

**Accepted:** April 16, 2015

**Published:** May 22, 2015

**Copyright:** © 2015 Wagner et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** Sample locations and microsatellite data: DRYAD doi: 10.5061/dryad.h25hj.

**Funding:** This project was funded by the German Research Foundation (DFG) and by the EU-funded project Trees4Future.

**Competing Interests:** The authors have declared that no competing interests exist.

## Introduction

Translocations are defined as intentional or accidental anthropogenic movements of living or- ganisms from one place to another ([1], S1 Text). Translocations taking place within a species range have long been performed for agricultural and horticultural purposes and are used
increasingly for restocking, reintroduction and assisted migration [2–6]. Like biological invasions, within-range translocations are raising questions on the existence of single versus multiple introduction events, the source of non-native material, and the speed of adaptation to local conditions. However, admixture between native and non-native material is a specificity of within-range translocations with unique evolutionary implications [2, 4, 5]. In principle, translocations can have positive consequences (increased genetic diversity, performance and density of local populations) or detrimental ones (competition, hybridization, genetic swamping and outbreeding depression) [2, 4, 7, 8].

Differentiating native from non-native material can provide critical information for the conservation of genetic resources [2, 9, 10]. Unfortunately, historical information on source and release areas of transferred material is often insufficient or lacking and non-native individuals cannot in most cases be differentiated from native ones by morphological features, especially when admixture between native and non-native material has already taken place [11]. Hence, to date, in spite of their widespread occurrence and increasing frequency across the planet, within-range translocations of undomesticated species have remained understudied [7].

Among the few existing studies, some have focused on admixture between native and non-native stocks of fishes for conservation purposes by relying on Bayesian approaches [12–16]. These studies have demonstrated that restocking can have contrasting consequences, with hybridization rates depending not only on stocking intensity but also on environmental factors. Forest trees have also been the object of a few investigations that focussed on tracing non-native material introduced into native populations [17–19]. These studies, relying on maternally inherited organelle markers, have demonstrate that where translocations had taken place, genetic diversity had increased and genetic differentiation and spatial genetic structure had decreased [18]. They also illustrate that it is not possible to identify a priori a native reference sample completely exempt of translocations, as forest history cannot in most cases be reconstructed in detail at the individual stand level [5, 10, 20]. To date, however, no study has explored translocations from a range-wide perspective using a systematic approach.

European forest trees represent particularly good models to investigate within-range translocations. Historical records indicate that forest reproductive material (typically seed or seedlings) has been traded across trees' native ranges for centuries. However, the overall impact is not clear, as native material might or not be better adapted than non-native one. Moreover, the onset of massive plantations is often quite recent and admixture proceeds on a rather long time scale due to the long generation times of trees [5, 9, 10, 20–22]. Among trees, European larch (*Larix decidua* Mill.) is a good candidate to study within-range translocations due to its restricted pollen dispersal ability that should result in high population structure [23, 24]. Its current fragmented native distribution range was reached in the early Holocene [25, 26], which must also have favoured the build-up of a distinct genetic structure that should help identify translocations. *Larix decidua* is found at high altitudes in the Alps and in other Central European mountains (Carpathians, Sudetes) as well as in the Polish lowlands (collectively referred to as Central Europe; see S1 Fig). In Central Europe, population sizes are rather small due to a particularly pronounced fragmentation in this part of the range. According to historical documents, plantations of larch in Europe have been performed since the sixteenth century and have culminated in the nineteenth century [27]. They have reached such a popularity that forecasters used to speak of a "Lärchenmanie" (German) or a "manie du mélèze" (French) [27–31]. Reports indicate that seed sources from the Tyrol (Austrian and Italian Central Alps) and the Sudetes (mountains at the border between Germany, Czech Republic and Poland) have been used preferentially [29]. A consequence of these widespread translocations was the concomitant expansion of larch canker that contributed to the decline of larch plantations during the twentieth century. The outbreak of this disease has been attributed to the use of maladapted
canker-sensitive Alpine material in afforestation (e.g. [32]). In contrast, material from Central Europe, especially from the Sudetes Mountains, is known to be particularly canker-resistant.

The present study focuses on the following questions: Can translocations be accurately identified in *L. decidua*, to the point that we can safely reconstruct the ancient genetic structure and identify native material? Are translocations more frequently detected in some parts of the range than in others? Have some seed sources been overrepresented among the transferred material? Have multiple geographic sources been used in plantations? Have translocations been followed by admixture? To answer these questions, we rely on molecular markers and population genetic methods (GENECLASS and STRUCTURE [33–35]) to systematically detect translocation and admixture events. We establish a reference of native genotypes, compare the ancient genetic structure with the extant genetic structure, and study the intensity, distribution and sources of translocations.

## Material and Methods

### Plant material and DNA isolation

Ascertaining the extent of the native range of *Larix decidua*, in particular in Central Europe, is difficult due to its long plantation history. With permission of the Nordwestdeutsche Forstliche Versuchsanstalt (NW-FVA), we relied on a comprehensive collection of populations sampled in 1957–1958 within the native range to establish an international provenance trial (IUFRO, [27, 32, 36]). Although collecting material *ex situ* rather than *in situ* can add error sources such as seed lot contamination, this particular collection seemed very valuable for our study as it is based on a precise compilation of historical information about the species native range and on large seed samples that should be representative of the original stands [27, 29, 32, 36, 37]. Moreover, some populations sampled at the time might have become eradicated or intermixed with exotic plantations during the more than 50 years that have elapsed since the establishment of this reference collection. In 2010, we collected 32 populations from this trial, including three populations originating from the Swietokrzyskie Mountains and from south of Warsaw in Poland that correspond to the *polonica* variety [38] listed as vulnerable according to IUCN standards (S1 Table, S1 Fig). To improve the coverage of the species native range, we collected eight other populations *in situ*. For comparison, we sampled five populations from beyond the native range as well as one *L. sibirica* population from a Swedish provenance trial. In each population, phloem or needles were sampled from 24 individuals. Sampling technique and DNA isolation was as described before [39].

### Microsatellite genotyping

We genotyped 24 individuals per population at 13 highly variable microsatellite loci combined in two multiplexes [39]. Genotype scoring and allele binning based on raw allele sizes had been described before [40]. No correction for null alleles was performed as they had been shown to be rare [39].

### Mitochondrial DNA sequencing

We sequenced two mtDNA regions (UBC460 and *atpA*) from eight individuals per population. Primer sequences for UBC460 are provided in [41] while primer sequences for *atpA* (*L3atp-2-1* 5'- GCGGCTGCCTATAGATACGA-3' and *L3atpfl/1* 5'- GCTACCGAGGCAGATATGGA-3') are based on published sequences [42]. PCR reactions were performed in a volume of 25 µl with 1× PCR buffer, 0.24 mM of each dNTP, 0.16 mg/ml of BSA, 0.2 µM of each primer, 1 U of Taq polymerase (DreamTaq Green, Fermentas) and 20 ng of genomic DNA. PCR programs started
with an initial denaturation step at 94°C for 4 min followed by 30 (UBC460) or 28 (atpA) cycles of denaturation at 94°C for 1 min, annealing at 65°C (UBC460) or 57°C (atpA) for 45 s and elongation at 72°C for 2 min 50 s (UBC460) or 1 min 50 s (atpA). The cycles were followed by a final elongation step at 72°C for 10 min. PCR products were sent to LGC Genomics (Berlin, Germany) for clean-up and sequencing. Sequence editing and alignment were performed using CodonCode Aligner Version 4.0.3. Haplotypes were detected after concatenating both sequenced regions and their relative frequencies were calculated using GenAlEx 6.41 [43]. For reconstructing minimum spanning networks, we focused on the variable sites, counting each indel as a single event, regardless of its size. The atpA region contained a minisatellite with two different repeat motifs and a variable number of repeats. This variation was coded into single nucleotide positions to compute a minimum spanning network using TCS 1.21 [44] with a fixed 12-step connection limit.

## Recent migrant detection

The software GENECLASS [33] has been designed to identify recent migrants using nuclear multi-locus genotypes. It derives the likelihood of each individual to belong to the population in which it has been sampled. We used it to identify translocation events based on microsatellite genotypes, thereby assuming that non-native individuals occur at low frequency within native populations and that admixture has remained limited. We used simulations and a resampling method [45] to generate 10,000 individuals based on allele frequency distribution for each population. Type I error to falsely identify an individual as a migrant was set at 0.01. A stepwise approach was used to improve detection rates in populations harbouring multiple migrants. The procedure was stopped after three successive runs to keep type I error low.

## Nuclear cluster analysis

Bayesian cluster analysis relied on STRUCTURE 2.1 [34, 35]. The number of clusters (K) was set from one to nine, with 10 runs for each K value. We used a burn-in of 200,000 and 1,000,000 iterations and explored both the no-admixture and the admixture model. To determine the most likely value of K, we computed the posterior probabilities and the second order rate of change L"(K) with STRUCTURE HARVESTER [46, 47]. We studied the genetic distances between clusters by deriving a neighbour joining tree from the pairwise distance matrix of the STRUCTURE output using the *R* package APE [48]. In addition, we ran a Principle Coordinate Analysis [43] to identify major genetic groups independently.

## Detection of translocation and admixture events

In each population, we consider as resulting from translocation events all purebred individuals (i.e. individuals assigned to one ancestral group at a given threshold score based on the STRUCTURE results) that do not belong to the most frequent ancestral group. By so doing, we assume that (i) prior to translocations a strong ancestral genetic structure existed, with no or limited admixture before the onset of translocations by humans, that (ii) native individuals of each population originate from a single ancestral group, and that (iii) native individuals still represent the majority of individuals. We used simulations to study risks associated with that procedure.

As a starting point to work on genetic exchange, we decided to focus on only two genetic groups, those corresponding to the first axis of the Principal Coordinate Analysis, to accurately detect translocation and admixture events. To confirm our expectation that these two cluster groups are sufficiently differentiated for this purpose, we simulated crosses within the corresponding cluster groups using HYBRIDLAB [49]. The simulated individuals were analysed by
STRUCTURE and then assigned back to the two cluster groups to evaluate assignment errors. Two different assignment thresholds were explored (0.875 and 0.9375, for details see S2 Fig). In both cases, a selection of the original genotypes was used as learning samples (USEPOPINFO option) and the simulated genotypes were included as supplementary genotypes, i.e. they were not part of the analysis.

## Cytonuclear association

We studied cytonuclear association to investigate if the two cluster groups have only recently become in contact, which would point to human-induced translocations. The rationale is that, if two populations that initially differ at both organelle and nuclear genomes come into contact, individuals exhibiting organelle genomes of one population and nuclear genomes typical for the other population will be generated within a few generations through crossbreeding, assuming selective neutrality. Conversely, if not all possible cytonuclear combinations are detected in a given population, very recent human-induced mixing can be inferred.

## Reconstruction of the ancestral genetic structure

We combined data from GENECLASS and STRUCTURE to identify all individuals that are not fully native (recent migrants identified by GENECLASS and non-native purebred and admixed individuals identified by STRUCTURE). We then compared summary statistics before and after excluding these individuals, focusing on allelic richness $N_A$, gene diversity $H_S$, inbreeding coefficient $F_{IS}$ and fixation index $F_{ST}$ measured using FSTAT version 2.9.3.2 [50] with 15,000 permutations. For mtDNA data, we compared the diversity index $h_S$ and the coefficient of genetic differentiation $G_{ST}$ [51] using PERMUT (http://www.pierroton.inra.fr/genetics/labo/Software/).

## Patterns of translocation

To study if the two regions (Alps and Central Europe) had experienced asymmetric exchanges of genetic material, we compared the proportion of translocation events identified with STRUCTURE in each region using a z-test. Then, to investigate the origin of the material used for translocation, we compared the genetic composition of non-native individuals (recent migrants detected by GENECLASS and non-native purebreds detected by STRUCTURE) and native individuals (all remaining purebreds). In this analysis, to avoid including false positives, admixed individuals identified by STRUCTURE were not included unless they had been identified by GENECLASS as being recent migrants. Finally, to study if translocations had involved multiple geographic sources, we counted the number of genetic clusters identified by STRUCTURE in populations with versus without evidence for translocation. An individual was considered to belong to a given nuclear cluster if its ancestry was mainly found in this cluster (assignment value >0.5). In populations without evidence for translocation, we considered all clusters except one (the native one); in populations with evidence for translocation, we considered all clusters except two, the native and the most frequent non-native one. We used a Wilcoxon rank-sum test based on the test statistic W to compare the mean number of extra-clusters between these two groups of populations [52].

## Patterns of admixture

First, we compared the proportion of admixture events identified with STRUCTURE in the Alpine versus the Central European region using a z-test. Second, we tested if translocations eventually result in admixture using two different approaches. In a first approach, we predicted that the proportion of admixed individuals should be higher in populations in which translocation
events have been detected than in the remaining populations, as contact between populations necessarily precedes admixture. To test this expectation, populations were grouped in two categories: populations with evidence for translocation and populations without such evidence, as explained above. Mean admixture level was then compared between these categories (after excluding individuals that had provided evidence for translocation), using a one-sided *t*-test. In a second approach, we relied on mitochondrial data to study if translocations eventually result in admixture. As individuals of non-native origin should typically be present at low frequency and will then be more likely to mate with local native individuals, whereas local native individuals should mostly mate with each other, the proportion of individuals with admixed genomes should be higher among individuals exhibiting non-native mitochondrial haplotypes than among individuals possessing native haplotypes. These two proportions were measured in each region (Alps and Central Europe) and compared using a *z*-test.

## Results

### Microsatellite genotyping and mitochondrial DNA sequencing

A sample of 1026 *Larix decidua* individuals from 45 populations were successfully genotyped at 13 microsatellite loci. We detected on average 22 alleles per locus. Sequencing of the two mtDNA fragments in 363 individuals from 43 populations led to the detection of 22 mtDNA haplotypes. The observed molecular variation was a combination of nucleotide substitutions, insertions/deletions as well as copy number variation of the minisatellite (S2 Table). To illustrate the relationships between haplotypes we constructed a minimum spanning network (Fig 1A). The eight *Larix sibirica* samples are characterized by a single private haplotype (H23) that has a central position in the network. The *L.decidua* haplotypes form two groups (i.e. lineages), each composed of two frequent haplotypes and some less frequent ones. In lineage 1 (H16, H17, H18 and H22), haplotypes H16 and H18 were found in 58% of the individuals. In lineage 2 (H8, H9, H10, H19, H20, H21 and 12 additional rare haplotypes differing in copy number of one minisatellite motif), haplotypes H9 and H10 occurred in 28% of the individuals (Fig 1B).

### Recent migrant detection

Three nested runs of *GENECLASS* led to the detection of 79 (9%) recent migrants in 33 of the 40 native populations (S3 Fig). Many populations were characterized by one migrant (10 populations), two migrants (9 populations) or three migrants (8 populations). In addition, there were three populations with four migrants and another three with five migrants. In seven populations no migrant was detected.

### Bayesian cluster analyses

Cluster analyses of the 1026 *Larix decidua* genotypes revealed support for the existence of three ancestral groups (*K* = 3) using the model without admixture (S4 Fig). The model with admixture identified a similar solution at *K* = 3 but detected another solution at *K* = 7 supported by the posterior probabilities and the second order rate of change (S5 Fig, Fig 2A and 2B). This latter solution identifies clusters that are hierarchically nested within the three larger clusters. We selected the solution at *K* = 7, as it provides increased resolution when reconstructing the ancestral genetic structure across the species range. To study translocations, we grouped the seven clusters into two particularly strongly differentiated groups: group 1 made of clusters 1–4 (corresponding to clusters 1+2 in the solution at *K* = 3, these two clusters being the most closely related) and group 2 made of clusters 5–7 (corresponding to cluster 3 in the solution at *K* = 3).
<figure><img src="image_4.png" /><figcaption>**Fig 1. Results of mitochondrial analyses.** Minimum spanning network (A), distribution of combined mitochondrial haplotypes before (B) and after (C) translocation removal. In (A) circles represent haplotypes colored by their codes and scaled to their frequencies. Haplotypes caused by minisatellite variation are in purple. White circles correspond to predicted haplotypes that were not observed. Branches correspond to single mutations regardless of their length. In (B) and (C) circles represent the haplotype composition of the 40 presumably native populations and three populations sampled outside the range (n°72, 73 and 79) (~8 individuals/population).</figcaption></figure>

doi:10.1371/journal.pone.0127516.g001
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Nuclear clusters detected by STRUCTURE. Neighbor joining tree (A) and distribution of the seven SSR clusters before (B) and after (C) removing presumably non-native genotypes. In (A) each rectangle represents a cluster. In (B) and (C) pie charts represent cluster composition of the 40 presumably native populations and 5 populations sampled outside the range (n°72, 73, 78, 79 and 80) (~24 individuals/population). Individuals with *q* values > 0.5 are colored, the remaining ones are in white.</figcaption></figure>

doi:10.1371/journal.pone.0127516.g002

These two groups were also distinguished on the first axis of the principle coordinate analysis, which explains 38% of the variation, compared to 23% for the second axis (S6 Fig). Simulations
<figure><img src="image_6.png" /><figcaption>**Fig 3.** Non-native purebreds and admixture events detected with STRUCTURE. Populations sampled within the native range and five populations sampled outside the native range (n°72, 73, 78, 79 and 80).</figcaption></figure>

doi:10.1371/journal.pone.0127516.g003

confirmed that assignment power is sufficient to differentiate between these two cluster groups. The two thresholds used to assign genotypes into purebred and admixed category led to similar inferences, so we relied on the more relaxed threshold (0.875) (S3 Table,). When simulating crosses between individuals belonging to the same cluster group, the resulting genotypes were never falsely assigned to the other cluster group. Moreover, the risk to falsely assign them to the admixed category was quite low: the rate of false positives was 13.3% of the simulated individuals in group 1 and 12.3% in group 2; even lower rates were observed under a more realistic scenario (S3 Table). These results confirm that assignment power is sufficient to differentiate between individuals originating from the two cluster groups (clusters 1–4: group 1; and 5–7: group 2). Applying the assignment threshold of ≥0.875 for each of the two cluster groups on the 40 populations, we detected 18 (2%) purebred individuals in populations otherwise dominated by the other cluster group (Fig 3). These 18 individuals, distributed in six Alpine and four Central European populations, were considered to represent recent translocation events. In five populations, one individual considered to result from translocation was detected, in four populations two such individuals were detected and in one population (population 51) as many as five individuals were detected. We also found 141 admixed individuals (16%) in 32 of the 40 populations (Fig 3).

## Cytonuclear analysis

We checked the association between nuclear assignments and mitochondrial lineages for all studied individuals sampled within the range that had been studied with both markers. We found 15 out of 40 populations harbouring both mitochondrial lineages. These 15 populations included the ten populations within which we had also detected both types of nuclear purebreds. Yet, at the individual level, all of the 166 nuclear genotypes assigned to group 1
<table><caption>Table 1. Counts of mito-nuclear genotype combinations found in 40 populations.</caption><thead><tr><th></th><th></th><th>nuclear group 1 purebred</th><th>nuclear group 2 purebred</th><th>admixed</th><th>total</th></tr></thead><tbody><tr><td>All (Alps +</td><td>mt lineage 1</td><td>166</td><td>1</td><td>37</td><td>204</td></tr><tr><td>Central Europe)</td><td>mt lineage 2</td><td>0</td><td>96</td><td>36</td><td>132</td></tr><tr><td>Alps</td><td>mt lineage 1</td><td>160</td><td>0</td><td>24</td><td>184</td></tr><tr><td></td><td>mt lineage 2</td><td>0</td><td>7</td><td>7</td><td>14</td></tr><tr><td>Central Europe</td><td>mt lineage 1</td><td>6</td><td>1</td><td>13</td><td>20</td></tr><tr><td></td><td>mt lineage 2</td><td>0</td><td>89</td><td>29</td><td>118</td></tr></tbody></table>

doi:10.1371/journal.pone.0127516.t001

purebreds carried mtDNA haplotypes from lineage 1, and all but one of the 97 group 2 purebred individual carried mtDNA haplotypes from lineage 2, indicating nearly total cytonuclear association (Table 1). This suggests that the contact between these divergent genetic groups is recent, thereby supporting our conclusion that individuals assigned to the locally uncommon nuclear group or mitochondrial lineage are not native.

## Reconstruction of the ancestral genetic structure

By combining results from both *GENECLASS* and *STRUCTURE*, we detected 208 (23%) presumably non-native individuals distributed in 38 of the 40 populations. There were only 30 individuals identified by both methods (six purebred and 24 admixed individuals). In particular, using *GENECLASS*, we identified 49 migrants not detected with *STRUCTURE*. As expected, *GENECLASS* was less efficient at detecting migrants in populations characterized by high rates of admixture (S7 Fig). Geographic structure obtained after removing these 208 individuals is shown in Figs 1C and 2C. For nuclear data, global summary statistics before and after removing these individuals differed significantly for all measures, with the trimmed dataset having lower within population diversity, increased differentiation among populations and being closer to Hardy-Weinberg equilibrium than the original dataset ($N_A = 6.2$ vs. 6.6, $H_S = 0.74$ vs. 0.76, $F_{ST} = 0.11$ vs. 0.08, $F_{IS} = 0.02$ vs. 0.05; $p < 0.001$ in all cases). Note however that a few populations still had comparatively high heterozygote deficit after removing individuals considered to be non-native ($F_{IS} > 0.08$ for populations 3, 53, 66 and 84, S4 Table,). Genetic diversity at mtDNA markers was also lower and genetic structure higher in the trimmed dataset ($h_S = 0.25$, $G_{ST} = 0.67$) than in the original one ($h_S = 0.37$, $G_{ST} = 0.52$). In fact, after trimming, most populations became dominated by a single genotypic cluster. Two large regions of the range, the Alps and Central Europe, turned out to be composed of individuals characterized by reciprocally exclusive nuclear clusters and mitochondrial lineages: all populations of the Alpine region were dominated by individuals with nuclear genotypes assigned to cluster group 1 and characterized by mitochondrial lineage 1, and all populations from Central Europe were dominated by individuals with nuclear genotypes assigned to cluster group 2 and characterized by mitochondrial lineage 2, making it possible to compare in detail the translocations across these two regions.

## Patterns of translocation

We found evidence for translocation in larch populations collected both *in situ* and *ex situ* in provenance tests, indicating that the collection mode did not strongly influence results. The proportion of individuals considered to originate from translocation events was not significantly different in the Alpine and in the Central European region (2% versus 3%, $z = -1.4$, $p = 0.17$). Considering results from both *GENECLASS* and *STRUCTURE*, we found that one of the four Alpine clusters (cluster 2) was over-represented in non-native material from Central Europe (50%; S8 Fig). This cluster is made up of individuals from populations of the Austrian and Italian
Central Alps, including from the Tyrol (S1 Fig). Genotypes with this cluster were also overrepresented when examining translocations within the Alpine region (S8 Fig). In contrast, none of the three Central European clusters was over-represented when considering translocations from Central Europe to the Alps or within Central Europe. To check if each population typically harbours translocations from a single source or from several ones, we focussed on the Alpine region where more populations had been sampled and where the geographic structure is stronger. We found a mean of 1.8 non-native additional clusters in the six populations with evidence for translocation, compared to only 0.8 in the 18 populations without evidence for translocation. This difference is significant (*W* > 148, *p* < 0.001), suggesting that when translocations occur in a given population, they often involve material from more than one source.

## Patterns of admixture

The proportion of admixed individuals was significantly lower in Alpine populations than in Central European ones (10% versus 23%, *z* = -1.5, *p* < 10⁻⁵). Furthermore, mean admixture values were significantly higher in populations with evidence for translocation than in populations without evidence for translocation (Alps: *q* = 0.08 versus 0.05, *p* < 10⁻⁴, Central Europe: *q* = 0.19 versus 0.09, *p* < 10⁻⁴, Fig 4). Finally, admixture levels were higher in individuals harbouring the non-native mitochondrial lineage: in the Alps, 50% (7/14) of the individuals with mtDNA haplotypes from lineage 2 were admixed, compared to only 13% (24/184) for individuals with mtDNA haplotypes from lineage 1 (*z* = -3.7, *p* < 10⁻³). Similarly, in populations from Central Europe, 65% (13/20) of the individuals with mtDNA haplotypes from lineage 1 were admixed compared to only 25% (29/118) for individuals with mtDNA haplotypes from lineage 2 (*z* = -3.6, *p* < 10⁻³).

## Genetic composition of populations sampled outside the species range

All five populations were located in Central Europe. Yet 59% of all individuals sampled were made of pure Alpine material, i.e. a much higher proportion than for the 16 reference populations from Central Europe (3%, Fig 3). There were also 15 admixed individuals in these five populations. In the three populations where mitochondrial data was available (populations 72, 73 and 79), strict cytonuclear association was observed, as in the native range, pointing to limited introgression and hence recent plantations.

<figure><img src="image_7.png" /><figcaption>**Fig 4.** Admixture in Central European populations. Populations with evidence for recent translocation (A) and without evidence (B). Bars represent individual *q*-values (group 2 in white, group 1 in grey).</figcaption></figure>

doi:10.1371/journal.pone.0127516.g004
Discussion

Selecting suitable methods to detect translocations and investigate their consequences is chal-
lenging, as it implies differentiating conspecific individuals that originate from different parts
of the range in a context where admixture has already taken place. Yet, this detection task is of
fundamental importance to determine the status of extant populations [2, 9, 12, 14, 53, 54]. In
this study, we took advantage of several suitable features of *Larix decidua*, including the persis-
tence of native genotypes in all parts of the range and a strong geographic structure at both nu-
clear and mitochondrial markers.

To systematically detect translocations, we first relied on a classical assignment method to exclude reference populations as the origin of individuals on the basis of multilocus genotype data (GENECLASS, [33]). The advantage of this approach is that there is no need to investigate the entire genetic structure of the species. However, in populations where high rates of admixture had been detected, GENECLASS was less effective at detecting introduced material, showing that this method also has limitations, even if some admixed individuals were identified as migrants. While recent migrants can in principle occur by natural means, they were conservative-ly considered to represent non-native material.

The other method used to identify translocations relies on the extant range-wide genetic
structure identified using STRUCTURE [35]. We showed using simulations (HYBRIDLAB, [49]) that
assignment power was sufficient to detect translocations. These translocations could also be
validated with mitochondrial data: individual genotypes identified with nuclear DNA data as
being non-native were also consistently shown to be non-native using mitochondrial data. Yet
only a few generations of crossbreeding would have sufficed to reshuffle mitochondrial lineages
and nuclear clusters through hybridization and backcrosses, suggesting that genetic exchanges
between the two ancestral genetic groups have been limited. This point supports our assertion
that the corresponding individuals found in minority in a population consist indeed mostly of
recent long-distance translocation events. The limitation of this second approach is that it de-
pends on our ability to detect a strong ancestral genetic structure. In our study, we conserva-
tively considered only two ancestral groups of populations, to reduce assignment errors,
thereby trading power for accuracy. Hence, only a subset of the translocation events could be
identified, thus underestimating the frequency of translocations and biasing the results towards
the identification of long distance translocation events. This should be kept in mind when in-
terpreting the results.

The two approaches used to detect translocation and admixture events turned out to be complementary as they identified largely non-overlapping sets of presumably non-native individuals. Once all these individuals had been removed, a clear-cut geographic pattern emerged. In particular, all populations from the Alps turned out to be genetically different from Central European populations. Hence, we were not compelled to resort to post hoc arguments to dismiss populations from our reference set. In fact, some populations from Central Europe had a very high proportion of non-native genotypes, and the pre-plantation genetic structure that we inferred is probably still not completely devoid of exotic material, as shown by the persistence of relatively high heterozygote deficit in some populations.

We found that translocations had preferentially involved one of the four Alpine genetic clusters. In the populations sampled outside of the species range, genotypes assigned to this cluster were also particularly frequent. This finding is in agreement with historical information indicating that *L. decidua* seeds from the Tyrol had been intensively used in plantations [29]. In contrast, we could not confirm a preferential use for plantations of seeds from the Sudetes, as documented in the literature [26].
Our results also show that translocations had typically involved multiple geographic sources. This result is reminiscent of earlier studies in oaks based on chloroplast DNA, where planted stands were shown to have much greater diversity than native ones, likely as a consequence of the use of mixed seed lots during plantations [55]. Another possible explanation for this observation is that there had been several episodes of plantations at each site, involving different seed sources.

We found strong support for a link between translocation and admixture, as expected if admixture events had taken place subsequently to translocations, once planted trees have become mature. First, we found an increased proportion of admixed individuals in populations where one or more non-native purebred genotypes had been detected, pointing to the recent onset of admixture between native and non-native trees. Second, individuals with non-native mitochondrial haplotypes also had elevated rates of admixture. Although some admixed individuals are likely false positives, the previous findings as well as the high proportion and broad distribution of admixture suggest that in many places at least one generation has elapsed since the first trees were introduced, ruling out seed lot contamination in plantations as the main source of translocation events. The particularly high admixture rate detected in Central Europe makes sense as the distribution of *L. decidua* in Central Europe is more fragmented than in the Alps and populations there are of rather small size, making them more susceptible to genetic swamping [4, 5]. Translocations and the associated admixture process uncovered in this study are probably a result of the intensive periods of larch planting during the eighteenth and nineteenth century [27–31]. Hence, our results shed new light on the “Lärchenmanie” by uncovering its irreversible long-term consequences through a modification of the genetic structure of the species, including in areas where threatened gene pools have been recognized (*cf*. *L. decidua* variety *polonica*). Finally, the results for populations sampled beyond the native distribution range show that we are now in the position to start investigate their make-up and subsequent history.

In conclusion, our study illustrates the importance of investigating translocations across species range to efficiently and accurately trace non-native material of different geographic origins and to uncover native material as well as species ancient history. In the future, such approaches will provide essential information for conservation of genetic resources. For instance, in *L. decidua*, we could now in principle identify canker-sensitive Alpine material introduced across the range and native Polish populations threatened by introgression, using greater sample sizes. While further progress should be achieved through the development of new demo-genetic methods that simultaneously identify introductions, admixture events and their timing, existing population genetic tools can already be used more widely, thereby revealing the frequency and prevalence of translocations across the globe as well as their consequences.

**Supporting Information**

**S1 Fig. Current distribution range of *Larix decidua* and studied populations.**
(DOCX)

**S2 Fig. Simulation study.**
(DOCX)

**S3 Fig. Distribution of migrants detected by GENECLASS.**
(DOCX)

**S4 Fig. Individual bar plots of STRUCTURE model with and without admixture.**
(XLSX)
S5 Fig. Determining the number of populations with STRUCTURE using Evanno's method.
(DOCX)

S6 Fig. Results of Principle Coordinate Analyses.
(DOCX)

S7 Fig. Comparison GENECLASS and STRUCTURE results.
(DOCX)

S8 Fig. Genetic composition of presumed natural and non-native material.
(XLSX)

S1 Table. Sampling information of the studied *Larix decidua* populations.
(DOCX)

S2 Table. Nucleotide and minisatellite variation of 23 mitochondrial haplotypes.
(DOCX)

S3 Table. False positives detection by simulations.
(DOCX)

S4 Table. FIS before and after removing individuals considered to be non-native.
(DOCX)

S1 Text. Definitions.
(DOCX)

**Acknowledgments**

This project was funded by the German Research Foundation (DFG) and by the EU-funded project Trees4Future. Genotyping was performed in the Genome-Transcriptome facility of the Functional Genomic Centre of Bordeaux, with support from the National Research Agency (ANR-10-EQPX-16 XYLOFOREST). We thank Franck Salin, Sarah Monllor and Christophe Boury for support and assistance during genotyping and Vladimir Semerikov, Christina Huneck and Christina Mengel for their help in mitochondrial marker design. We thank Eric Welk for providing distributional data on *Larix decidua*. For assistance in the various fields of the project we thank François Ehrenmann, Raphaël Leblois, Mathieu Gauthier, Vanina Guérin, Patrick Léger and Erwan Guichoux, Luc Pâques, Jutta Mäsgen, Laurent Bouffier, Frédéric Lagane, Christoph Sperisen, Yuri Hayda, Ihor Neyko, Guy Bour, Czesław Koziol, Philippe Rozenberg, Leopoldo Sanchez, Frédéric Millier, Jean-Paul Charpentier, Owe Martinsson, Lars Karlman, Hermann Spellmann, the Ukrainian Research Institute for Mountain forestry, the Polish Osusznica Forest Institute, and the Polish Maskulińskie Forest Institute. For helpful comments on the manuscript we thank Arndt Hampe, Guillaume de Lafontaine and Stéphanie Mariette.

**Author Contributions**

Conceived and designed the experiments: SW RJP SL. Performed the experiments: SW SL. Analyzed the data: SW RJP SL. Contributed reagents/materials/analysis tools: SW RJP SL. Wrote the paper: SW RJP SL SG.

**References**

1. IUCN (1987) IUCN position statement on the translocation of living organisms: introductions, re-introduction, and re-stocking. IUCN, Gland, Switzerland. Available: http://intranet.iucn.org/webfiles/doc/
SSC/SSCwebsite/Policy_statements/IUCN_Position_Statement_on_Translocation_of_Living_Organisms.pdf.

2. Allendorf FW, Leary RF, Spruell P, Wenburg JK (2001) The problems with hybrids: setting conservation guidelines. *Trends in Ecology & Evolution* 16: 613–622.

3. Armstrong DP, Seddon PJ (2008) Directions in reintroduction biology. *Trends in Ecology & Evolution* 23: 20–25.

4. Hufford KM, Mazer SJ (2003) Plant ecotypes: genetic differentiation in the age of ecological restoration. *Trends in Ecology & Evolution* 18: 147–155.

5. Lefèvre F (2004) Human impacts on forest genetic resources in the temperate zone: an updated review. *Forest Ecology and Management* 197: 257–271.

6. Seddon PJ (2010) From reintroduction to assisted colonization: moving along the conservation translocation spectrum. *Restoration Ecology* 18: 796–802.

7. Largiadèr CR (2007) Hybridization and introgression between native and alien species. In: Nentwig W., editor editors. *Biol Invasions*. Berlin, Heidelberg: Springer pp. 275–292. PMID: 18001943

8. Saltonstall K (2011) Remnant native *Phragmites australis* maintains genetic diversity despite multiple threats. *Conservation Genetics* 12: 1027–1033.

9. Koskela J, Lefèvre F, Schueler S, Kraigher H, Olrik DC, et al. (2013) Translating conservation genetics into management: Pan-European minimum requirements for dynamic conservation units of forest tree genetic diversity. *Biological Conservation* 157: 39–49.

10. Bradshaw RH (2004) Past anthropogenic influence on European forests and some possible genetic consequences. *Forest Ecology and Management* 197: 203–212.

11. Geburek T (1997) Isozymes and DNA markers in gene conservation of forest trees. *Biodiversity & Conservation* 6: 1639–1654.

12. Gum B, Gross R, Kuehn R (2006) Discriminating the impact of recent human mediated stock transfer from historical gene flow on genetic structure of European grayling *Thymallus thymallus* L. *Journal of Fish Biology* 69: 115–135.

13. Marie A D, Bernatchez L, Garant D (2012) Environmental factors correlate with hybridization in stocked brook charr (*Salvelinus fontinalis*). *Canadian Journal of Fisheries and Aquatic Sciences* 69: 884–893.

14. Neville HM, Dunham JB (2011) Patterns of hybridization of nonnative cutthroat trout and hatchery rainbow trout with native redband trout in the Boise River, Idaho. *North American Journal of Fisheries Management* 31: 1163–1176.

15. Scribner K, Page K, Bartron M (2000) Hybridization in freshwater fishes: a review of case studies and cytonuclear methods of biological inference. *Reviews in Fish Biology and Fisheries* 10: 293–323.

16. Winkler KA, Pamminger-Lahnsteiner B, Wanzenböck J, Weiss S (2011) Hybridization and restricted gene flow between native and introduced stocks of Alpine whitefish (*Coregonus* sp.) across multiple environments. *Molecular Ecology* 20: 456–472. doi: 10.1111/j.1365-294X.2010.04961.x PMID: 21199024

17. Deguilloux MF, Pemonge MH, Bertel L, Kremer A, Petit RJ (2003) Checking the geographical origin of oak wood: molecular and statistical tools. *Molecular Ecology* 12: 1629–1636. PMID: 12755890

18. König AO, Ziegenhagen B, van Dam BC, Csaikl UM, Coart E, et al. (2002) Chloroplast DNA variation of oaks in western Central Europe and genetic consequences of human influences. *Forest Ecology and Management* 156: 147–166.

19. Lowe A, Munro R, Samuel S, Cottrell J (2004) The utility and limitations of chloroplast DNA analysis for identifying native British oak stands and for guiding replanting strategy. *Forestry* 77: 335–347.

20. Lefèvre F (2002) Le peuplier noir: une ressource génétique à l'interface entre habitats naturels d'intérêt communautaire et sylviculture intensive. *Dossiers de l'Environnement de l'INRA* 21: 125–132. doi: 10.1016/j.annder.2008.07.055 PMID: 19061651

21. König A (2005) Provenance research: evaluating the spatial pattern of genetic variation. In: Geburek T. and Turok J., editors. *Conservation and Management of Forest Genetic Resources in Europe*. Zvolen, Slovakia: Arbora Publishers. pp. 275–328.

22. Petit RJ, Hampe A (2006) Some evolutionary consequences of being a tree. *Annual Review of Ecology, Evolution, and Systematics* 37: 187–214.

23. Sjögren P, Van der Knaap WO, Huusko A and van Leeuwen JFN (2008) Pollen productivity, dispersal, and correction factors for the major tree taxa in the Swiss Alps based on pollen trap results. *Review of Palaeobotany and Palynology* 152: 200–210.

24. Sjögren P, van der Knaap WO, Kaplan JO, van Leeuwen JFN, Ammann B (2008) A pilote stuy on pollen representation of mountain valley vegetation in the central Alps. *Review of Palaeobotany and Palynology* 149: 208–218.
25. Huntley B, Birks HJB (1983) An atlas of past and present pollen maps for Europe: 0–13000 years ago. Cambridge: Cambridge University Press. 667 p.

26. Lang G (1994) Quartäre Vegetationsgeschichte Europas: Methoden und Ergebnisse. Jena: Gustav Fischer Verlag. 462 p.

27. Rubner K (1953) Die pflanzengeographischen Grundlagen des Waldbaus. Radebeul, Berlin: Neumann. 583 p.

28. Fourchy P (1952) Etudes sur l'Ecologie et la Sylviculture du Mélèze (Larix europaea D. C.). Écologie du mélèze particulièrement dans les Alpes Françaises. Nancy: Ecole Nationale des Eaux et Forêts. pp. 137.

29. Münch E (1936) Das Lärchensterben. Forstwissenschaftliches Zentralblatt 58: 641–671.

30. Schober R (1949) Die Lärche—Eine ertragskundlich-biologische Untersuchung. Hannover: M. & H. Schaper Verlag. 285 p. doi: 10.1007/s10147-012-0476-7 PMID: 23011101

31. Tschermak L (1935) Die natürliche Verbreitung der Lärche in den Ostalpen. Wien: Julius Springer Verlag. 361 p.

32. Schober R (1977) Vom II. internationalen Lärchenprovenienzversuch: Ein Beitrag zur Lärchenherkunftsfrage. Frankfurt am Main: J.D. Sauerländer's Verlag. 385 p.

33. Piry S, Alapetite A, Cornuet JM, Paetkau D, Baudouin L, et al. (2004) GENECLASS2: A software for genetic assignment and first-generation migrant detection. Journal of Heredity 95: 536–539. PMID: 15475402

34. Falush D, Stephens M, Pritchard JK (2003) Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics 164: 1567–1587. PMID: 12930761

35. Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–959. PMID: 10835412

36. Schober R (1985) Neue Ergebnisse des II. Internationalen Lärchenprovenienzversuches von 1958/1959 nach Aufnahmen von Teilversuchen in 11 europäischen Ländern und den U.S.A. Schriften aus der Forstlichen Fakultät der Universität Göttingen und der Niedersächsischen Forstlichen Versuchsanstalt. Frankfurt am Main: J.D. Sauerländer's Verlag pp. 1–164.

37. Pluess AR (2011) Pursuing glacier retreat: genetic structure of a rapidly expanding *Larix decidua* population. Molecular Ecology 20: 473–485. doi: 10.1111/j.1365-294X.2010.04972.x PMID: 21199030

38. Lewandowski A (1997) Genetic relationships between European and Siberian larch, *Larix* spp. (Pinaeae), studied by allozymes. Is the Polish larch a hybrid between these two species? Plant Systematics and Evolution 204: 65–73.

39. Wagner S, Gerber S, Petit RJ (2012) Two highly informative dinucleotide SSR multiplexes for the conifer *Larix decidua* (European larch). Molecular Ecology Resources 12: 717–725. doi: 10.1111/j.1755-0998.2012.03139.x PMID: 22487486

40. Guichoux E, Lagache L, Wagner S, Chaumeil P, Léger P, et al. (2011) Current trends in microsatellite genotyping. Molecular Ecology Resources 11: 591–611. doi: 10.1111/j.1755-0998.2011.03014.x PMID: 21565126

41. Semerikov VL, Vendramin GG, Sebastiani F, Lascoux M (2006) RAPD-derived, PCR-based mitochon-drial markers from *Larix* species and their usefulness in phylogeny. Conservation Genetics 7: 621–625.

42. Polezhaeva MA, Lascoux M, Semerikov VL (2010) Cytoplasmic DNA variation and biogeography of *Larix* Mill. in Northeast Asia. Molecular Ecology 19: 1239–1252. doi: 10.1111/j.1365-294X.2010.04552.x PMID: 20163546

43. Peakall ROD, Smouse PE (2006) GENEALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Molecular Ecology Notes 6: 288–295.

44. Clement M, Posada D, Crandall KA (2000) TCS: a computer program to estimate gene genealogies. Molecular Ecology 9: 1657–1659. PMID: 11050560

45. Paetkau D, Slade R, Burden M, Estoup A (2004) Genetic assignment methods for the direct, real-time estimation of migration rate: a simulation-based exploration of accuracy and power. Molecular Ecology 13: 55–65. PMID: 14653788

46. Earl D, vonHoldt B (2012) STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conservation Genetics Resources 4: 359–361.

47. Evanno G, Regnaut S, Goudet J (2005) Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Molecular Ecology 14: 2611–2620. PMID: 15969739

48. Paradis E, Claude J, Strimmer K (2004) APE: analyses of phylogenetics and evolution in R language. Bioinformatics 20: 289–290. PMID: 14734327
49. Nielsen EE, Bach LA, Kotlicki P (2006) Hybridlab (version 1.0): a program for generating simulated hybridations from population samples. Molecular Ecology Notes 6: 971-973.

50. Goudet J (2001) FSTAT, a program to estimate and test gene diversities and fixation indices (version 2.9.3). Available: http://www.unil.ch/izea/softwares/fstat.html. Updated from Goudet (1995).

51. Pons O, Petit R (1996) Measuring and testing genetic differentiation with ordered versus unordered alleles. Genetics 144: 1237-1245. PMID: 8913764

52. Dodge Y (2008) The Concise Encyclopedia of Statistics New York: Springer. 622 p.

53. Hufford MB, Lubinsky P, Pyhjäjärvi T, Devengenzo MT, Ellstrand NC, et al. (2013) The genomic signature of crop-wild introgression in maize. PLOS Genetics 9: e1003477. doi: 10.1371/journal.pgen.1003477 PMID: 23671421

54. Rosenthal DM, Ramakrishnan AP, Cruzan MB (2008) Evidence for multiple sources of invasion and intraspecific hybridization in *Brachypodium sylvaticum* (Hudson) Beauv. in North America. Molecular Ecology 17: 4657-4669. doi: 10.1111/j.1365-294X.2008.03844.x PMID: 18627455

55. König AO, Ziegenhagen B, van Dam BC, Csaikl UM, Coart E, et al. (2002) Chloroplast DNA variation of oaks in western Central Europe and genetic consequences of human influences. Forest Ecology and Management 156: 147-166.