Cryptic niche conservatism among evolutionary lineages of an invasive lizard

Aim There is increasing evidence that the quality and breadth of ecological niches vary among individuals, populations, evolutionary lineages and therefore also across the range of a species. Sufficient knowledge about niche divergence among clades might thus be crucial for predicting the invasion potential of species. We tested for the first time whether evolutionary lineages of an invasive species vary in their climate niches and invasive potential. Furthermore, we tested whether lineage‐specific models show a better performance than combined models. Location Europe. Methods We used species distribution models (SDMs) based on climatic information at native and invasive ranges to test for intra‐specific niche divergence among mitochondrial DNA (mtDNA) clades of the invasive wall lizard Podarcis muralis. Using DNA barcoding, we assigned 77 invasive populations in Central Europe to eight geographically distinct evolutionary lineages. Niche similarity among lineages was assessed and the predictive power of a combination of clade‐specific SDMs was compared with a combined SDM using the pooled records of all lineages. Results We recorded eight different invasive mtDNA clades in Central Europe. The analysed clades had rather similar realized niches in their native and invasive ranges, whereas inter‐clade niche differentiation was comparatively strong. However, we found only a weak correlation between geographic origin (i.e. mtDNA clade) and invasive occurrences. Clades with narrow realized niches still became successful invaders far outside their native range, most probably due to broader fundamental niches. The combined model using data for all invasive lineages achieved a much better prediction of the invasive potential. Conclusions Our results indicate that the observed niche differentiation among evolutionary lineages is mainly driven by niche realization and not by differences in the fundamental niches. Such cryptic niche conservatism might hamper the success of clade‐specific niche modelling. Cryptic niche conservatism may in general explain the invasion success of species in areas with apparently unsuitable climate.     

Widespread primary,but geographically restricted secondary,human introductions of wall lizards,Podarcis muralis

Establishing the introduction pathways of alien species is a fundamental task in invasion biology. The common wall lizard, Podarcis muralis, has been widely introduced outside of its native range in both Europe and North America, primarily through escaped pets or deliberate release of animals from captive or wild populations. Here, we use Bayesian clustering, approximate Bayesian computation (ABC) methods and network analyses to reconstruct the origin and colonization history of 23 non‐native populations of wall lizards in England. Our analyses show that established populations in southern England originate from at least nine separate sources of animals from native populations in France and Italy. Secondary introductions from previously established non‐native populations were supported for eleven (47%) populations. In contrast to the primary introductions, secondary introductions were highly restricted geographically and appear to have occurred within a limited time frame rather than being increasingly common. Together, these data suggest that extant wall lizard populations in England are the result of isolated accidental and deliberate releases of imported animals since the 1970s, with only local translocation of animals from established non‐native populations. Given that populations introduced as recently as 25 years ago show evidence of having adapted to cool climate, discouraging further translocations may be important to prevent more extensive establishment on the south coast of England.     

Phylogeography and Conservation Genetics of the Common Wall Lizard,Podarcis muralis,on Islands at Its Northern Range

Populations at range limits are often characterized by lower genetic diversity, increased genetic isolation and differentiation relative to populations at the core of geographical ranges. Furthermore, it is increasingly recognized that populations situated at range limits might be the result of human introductions rather than natural dispersal. It is therefore important to document the origin and genetic diversity of marginal populations to establish conservation priorities. In this study, we investigate the phylogeography and genetic structure of peripheral populations of the common European wall lizard, Podarcis muralis, on Jersey (Channel Islands, UK) and in the Chausey archipelago. We sequenced a fragment of the mitochondrial cytochrome b gene in 200 individuals of P. muralis to infer the phylogeography of the island populations using Bayesian approaches. We also genotyped 484 individuals from 21 populations at 10 polymorphic microsatellite loci to evaluate the genetic structure and diversity of island and mainland (Western France) populations. We detected four unique haplotypes in the island populations that formed a sub-clade within the Western France clade. There was a significant reduction in genetic diversity (H_O, H_E and A_R) of the island populations in relation to the mainland. The small fragmented island populations at the northern range margin of the common wall lizard distribution are most likely native, with genetic differentiation reflecting isolation following sea level increase approximately 7000 BP. Genetic diversity is lower on islands than in marginal populations on the mainland, potentially as a result of early founder effects or long-term isolation. The combination of restriction to specific localities and an inability to expand their range into adjacent suitable locations might make the island populations more vulnerable to extinction.     

The desire for variety: Italian wall lizard (Podarcis siculus) populations introduced to the United States via the pet trade are derived from multiple native-range sources

Tests of invasion success often require comparisons between introduced and native populations, but determining the native-range sources for introduced populations can be difficult. Molecular markers can help clarify the geographic extent of native-range sources, helping to identify which populations are appropriate for comparative studies. The Italian Wall Lizard (Podarcis siculus) was introduced multiple times to the United States with extant populations in California, Kansas, New Jersey, and New York. We used phylogeographic analysis of mtDNA sequences (cytb gene) for individuals sampled from these introduced populations and across the native range to identify the number of independent introductions and the location of the source populations. Haplotypes sampled from introduced populations were nested within three geographically distinct, well-supported clades that together encompassed a large portion of the native range. Combining these phylogeographic results with documentation of the introductions revealed putative sources: California individuals are derived from Sicily; Kansas and New York populations are from Tuscany near Florence; and the New Jersey population is likely from the Adriatic coastal region, but a more specific locality is not possible. The pet trade dominates the invasion pathway for P. siculus introductions to the US. The genetically and geographically diverse sampling of its native range may be driven by the desire for phenotypic variety in the pet trade, a hypothesis that needs future testing.     

Adaptive responses to cool climate promotes persistence of a non-native lizard

Successful establishment and range expansion of non-native species often require rapid accommodation of novel environments. Here, we use common-garden experiments to demonstrate parallel adaptive evolutionary response to a cool climate in populations of wall lizards (Podarcis muralis) introduced from southern Europe into England. Low soil temperatures in the introduced range delay hatching, which generates directional selection for a shorter incubation period. Non-native lizards from two separate lineages have responded to this selection by retaining their embryos for longer before oviposition—hence reducing the time needed to complete embryogenesis in the nest—and by an increased developmental rate at low temperatures. This divergence mirrors local adaptation across latitudes and altitudes within widely distributed species and suggests that evolutionary responses to climate can be very rapid. When extrapolated to soil temperatures encountered in nests within the introduced range, embryo retention and faster developmental rate result in one to several weeks earlier emergence compared with the ancestral state. We show that this difference translates into substantial survival benefits for offspring. This should promote short- and long-term persistence of non-native populations, and ultimately enable expansion into areas that would be unattainable with incubation duration representative of the native range.     

Molecular markers reconstruct the invasion history of variable leaf watermilfoil (<Emphasis Type="Italic">Myriophyllum heterophyllum</Emphasis>) and distinguish it from closely related species

Genetic variation is increasingly recognized as an important factor influencing the establishment and spread of introduced species, and depends on both the introduction history and partitioning of genetic variation within and among potential source populations. We examine patterns of genetic variation in native and introduced populations of variable leaf watermilfoil, Myriophyllum heterophyllum, using chloroplast (trnL-F) and ribosomal (ITS) DNA sequences, as well as amplified fragment length polymorphisms (AFLPs). We identify a strong phylogeographic break distinguishing populations located on the Atlantic Coastal Plain (ACP) versus other (“Continental”) portions of the native range. Within these distinct biogeographic regions, we also find genetic variation to be strongly partitioned among populations as analysis of molecular variance (AMOVA) partitioned 91 and 75% of cpDNA and ITS diversity among populations, respectively. We demonstrate that the introduced ranges of variable leaf watermilfoil (northeastern and western US) result from multiple independent introductions from a variety of source populations, including lineages from both the ACP and Continental portions of the native range. In addition, we used our molecular markers to demonstrate that variable leaf watermilfoil is genetically distinct from three closely-related species that it is morphologically similar to. In particular, we demonstrate that M. heterophyllum is clearly distinct from a morphologically similar native species in the western US, M. hippuroides—whose distinctiveness from M. heterophyllum has been questioned—and therefore confirm the introduction of M. heterophyllum in the western US. Furthermore, we provide the first evidence for hybridization between these two species. Finally, our molecular markers identify previously unrecognized genetic variation in these four species, and therefore demonstrate the need for further taxonomic investigation.     

Effects of exotic species on evolutionary diversification

Exotic species invasions create almost ideal conditions for promoting evolutionary diversification: establishment of allopatric populations in new environmental conditions; altered ecological opportunities for native species; and new opportunities for hybridization between previously allopatric taxa. Here, we review recent studies of the evolutionary consequences of species invasions, revealing abundant and widespread examples of exotic species promoting evolutionary diversification via increased genetic differentiation among populations of both exotic and native species and the creation of new hybrid lineages. Our review indicates that, although the well-documented reductions to biodiversity caused by exotic species might outweigh the increases resulting from diversification, a complete understanding of the net effects of exotic species on biodiversity in the long term will require consideration of both.     

Genetic diversity and historical biogeography of the Maltese wall lizard,Podarcis filfolensis (Squamata: Lacertidae)

Podarcis filfolensis is an endemic lizard from the Maltese archipelago. There is evidence of human-mediated decline and even extirpation of some insular populations of this species. However, information about the intraspecific genetic diversity and phylogeographic patterns of this species is limited. Here we analyze genetic markers from a multi-locus dataset (mtDNA, 2,533 bp; nuclear c-mos gene, 353 bp; 11 microsatellites) for individuals from extant populations of P. filfolensis. Despite generally low genetic variability, two main mitochondrial groupings were clearly identified. In general, individuals from the main island of Malta were genetically distinct from those from Gozo, Comino, Cominotto and Small Blue Lagoon Rock, and also from Linosa and Lampione individuals. Three genetic clusters were detected based on microsatellite data: one was found at higher frequency on Malta, while the other two included samples from the remaining islands, showing some concordance with the mtDNA pattern. A time-calibrated Bayesian tree for the principal mitochondrial lineages indicated strong statistical support for two P. filfolensis lineages that originated in the Pleistocene (105.4–869 Ka). We show that these lineages largely meet the criteria for recognition as evolutionary significant units despite some recent admixture (possibly due to recent translocations between islands). Human disturbance, low genetic variability, evidence of bottlenecks and extirpation on one island indicate that a thorough review of the current conservation status of P. filfolensis would be timely.     

Climate and habitat configuration limit range expansion and patterns of dispersal in a non‐native lizard

Invasive species are one of the main causes of biodiversity loss worldwide. As introduced, populations increase in abundance and geographical range, so does the potential for negative impacts on native communities. As such, there is a need to better understand the processes driving range expansion as species become established in recipient landscapes. Through an investigation into capacity for population growth and range expansion of introduced populations of a non‐native lizard (Podarcis muralis), we aimed to demonstrate how multi‐scale factors influence spatial spread, population growth, and invasion potential in introduced species. We collated location records of P. muralis presence in England, UK through data collected from field surveys and a citizen science campaign. We used these data as input for presence‐background models to predict areas of climate suitability at a national‐scale (5 km resolution), and fine‐scale habitat suitability at the local scale (2 m resolution). We then integrated local models into an individual‐based modeling platform to simulate population dynamics and forecast range expansion for 10 populations in heterogeneous landscapes. National‐scale models indicated climate suitability has restricted the species to the southern parts of the UK, primarily by a latitudinal cline in overwintering conditions. Patterns of population growth and range expansion were related to differences in local landscape configuration and heterogeneity. Growth curves suggest populations could be in the early stages of exponential growth. However, annual rates of range expansion are predicted to be low (5–16 m). We conclude that extensive nationwide range expansion through secondary introduction is likely to be restricted by currently unsuitable climate beyond southern regions of the UK. However, exponential growth of local populations in habitats providing transport pathways is likely to increase opportunities for regional expansion. The broad habitat niche of P. muralis, coupled with configuration of habitat patches in the landscape, allows populations to increase locally with minimal dispersal.     

Is Gene Flow Promoting the Reversal of Pleistocene Divergence in the Mountain Chickadee (Poecile gambeli)?

The Pleistocene glacial cycles left a genetic legacy on taxa throughout the world; however, the persistence of genetic lineages that diverged during these cycles is dependent upon levels of gene flow and introgression. The consequences of secondary contact among taxa may reveal new insights into the history of the Pleistocene’s genetic legacy. Here, we use phylogeographic methods, using 20 nuclear loci from regional populations, to infer the consequences of secondary contact following divergence in the Mountain Chickadee (Poecile gambeli). Analysis of nuclear data identified two geographically-structured genetic groups, largely concordant with results from a previous mitochondrial DNA (mtDNA) study. Additionally, the estimated multilocus divergence times indicate a Pleistocene divergence, and are highly concordant with mtDNA. The previous mtDNA study showed a paucity of sympatry between clades, while nuclear patterns of gene flow show highly varied patterns between populations. The observed pattern of gene flow, from coalescent-based analyses, indicates southern populations in both clades exhibit little gene flow within or between clades, while northern populations are experiencing higher gene flow within and between clades. If this pattern were to persist, it is possible the historical legacy of Pleistocene divergence may be preserved in the southern populations only, and the northern populations would become a genetically diverse hybrid species.     

Molecular phylogeography and population structure of a mid-elevation montane frog Leptobrachium ailaonicum in a fragmented habitat of southwest China

Leptobrachium ailaonicum is a vulnerable anuran restricted to a patchy distribution associated with small mountain streams surrounded by forested slopes at mid-elevations (approximately 2000–2600 m) in the subtropical Mount Wuliang and Mount Ailao ranges in southwest China (Yunnan Province) and northern Vietnam. Given high habitat specificity and lack of suitable habitat in lower elevations between these ranges, we hypothesized limited gene flow between populations throughout its range. We used two mitochondrial genes to construct a phylogeographic pattern within this species in order to test our hypothesis. We also examined whether this phylogeographic pattern is a response to past geological events and/or climatic oscillations. A total of 1989 base pairs were obtained from 81 individuals of nine populations yielding 51 unique haplotypes. Both Bayesian and maximum parsimony phylogenetic analyses revealed four deeply divergent and reciprocally monophyletic mtDNA lineages that approximately correspond to four geographical regions separated by deep river valleys. These results suggest a long history of allopatric separation by vicariance. The distinct geographic distributions of four major clades and the estimated divergence time suggest spatial and temporal separations that coincide with climatic and paleogeographic changes following the orogeny and uplift of Mount Ailao during the late Miocene to mid Pliocene in southwest China. At the southern distribution, the presence of two sympatric yet differentiated clades in two areas are interpreted as a result of secondary contact between previously allopatric populations during cooler Pleistocene glacial cycles. Analysis of molecular variance indicates that most of the observed genetic variation occurs among the four regions implying long-term interruption of maternal gene flow, suggesting that L. ailaonicum may represent more than one distinct species and should at least be separated into four management units corresponding to these four geographic lineages for conservation.     

Molecular survey of morphological subspecies reveals new mitochondrial lineages in Podarcis muralis (Squamata: Lacertidae) from the Tuscan Archipelago (Italy)

Recent analyses of molecular markers have significantly revised the traditional taxonomy of Podarcis species (Squamata: Lacertidae), leading to critically reconsider the taxonomic value of several subspecies described only on morphological bases. In fact, lizards often exhibit high morphological plasticity both at the intra‐specific and the intra‐population level, especially on islands, where phenotypic divergences are mainly due to local adaptation, rather than to evolutionary differentiation. The Common wall lizard Podarcis muralis exhibits high morphological variability in biometry, pholidosis values and colour pattern. Molecular analyses have confirmed the key role played by the Italian Peninsula as a multi‐glacial refuge for P. muralis, pointing out the lack of congruence between mitochondrial lineages and the four peninsular subspecies currently recognized. Here, we analyse a portion of the protein‐encoding cytochrome b gene in the seven subspecies described for the Tuscan Archipelago (Italy), in order to test whether the mitochondrial haplotypes match the morphologically based taxonomy proposed for Common wall lizard. We also compare our haplotypes with all the others from the Italian Peninsula to investigate the presence of unique genetic lineages in insular populations. Our results do not agree completely with the subspecific division based on morphology. In particular, the phylogenetic analyses show that at least four subspecies are characterized by very similar haplotypes and fall into the same monophyletic clade, whereas the other three subspecies are closer to peninsular populations from central Italy. From these results, we conclude that at least some subspecies could be better regarded as simple eco‐phenotypes; in addition, we provide an explanation for the distinctiveness of exclusive lineages found in the archipelago, which constituted a refuge for this species during last glacial periods.     

Chemical communication,sexual selection,and introgression in wall lizards

Divergence in communication systems should influence the likelihood that individuals from different lineages interbreed, and consequently shape the direction and rate of hybridization. Here, we studied the role of chemical communication in hybridization, and its contribution to asymmetric and sexually selected introgression between two lineages of the common wall lizard (Podarcis muralis). Males of the two lineages differed in the chemical composition of their femoral secretions. Chemical profiles provided information regarding male secondary sexual characters, but the associations were variable and inconsistent between lineages. In experimental contact zones, chemical composition was weakly associated with male reproductive success, and did not predict the likelihood of hybridization. Consistent with these results, introgression of chemical profiles in a natural hybrid zone resembled that of neutral nuclear genetic markers overall, but one compound in particular (tocopherol methyl ether) matched closely the introgression of visual sexual characters. These results imply that associations among male chemical profiles, sexual characters, and reproductive success largely reflect transient and environmentally driven effects, and that genetic divergence in chemical composition is largely neutral. We therefore suggest that femoral secretions in wall lizards primarily provide information about residency and individual identity rather than function as sexual signals.     

The shaping of mitochondrial DNA phylogeographic patterns of the brown hare (Lepus europaeus) under the combined influence of Late Pleistocene climatic fluctuations and anthropogenic translocations

The phylogeographic structure of the brown hare (Lepus europaeus) was studied by analysing mtDNA control region sequences of 98 individuals from continental and insular Greece, Bulgaria, Cyprus and northern Israel, together with 44 published sequences from Italy and central Europe. We found two distinct clades separated by an average nucleotide divergence of 6.6%, which may correspond to a Balkan and to an Asia Minor refugium. The estimated time of separation of the two clades was dated back to 105,000- 490,000 years ago. These two clades coexist in the area of northeastern Greece and Bulgaria, most likely as a result of a post-glacial northward expansion. Within the southern Balkan refugium, network analyses showed geographical structuring, which supports the hypothesis of several isolated Late Pleistocene populations. The central European and Italian populations appear to have originated from a non-detected northern Balkan population that was genetically closely related to some northern Greek populations, as a result of postglacial expansion, translocations or a combination of both. Moreover, several cases of ancient and recent translocations by humans were detected, especially for some island populations, while the eastern Aegean islands off the Asia Minor coast were most likely colonized naturally through Late Pleistocene land bridge connection. The genetic analysis presented here provides a framework for designing proper conservation and management guidelines for this species.     

Male behaviour drives assortative reproduction during the initial stage of secondary contact

Phenotypic divergence in allopatry can facilitate speciation by reducing the likelihood that individuals of different lineages hybridize during secondary contact. However, few studies have established the causes of reproductive isolation in the crucial early stages of secondary contact. Here, we establish behavioural causes of assortative reproduction between two phenotypically divergent lineages of the European wall lizard (Podarcis muralis), which have recently come into secondary contact. Parentage was highly assortative in experimental contact zones. However, despite pronounced divergence in male phenotypes, including chemical and visual sexual signals, there was no evidence that females discriminated between males of the two lineages in staged interactions or under naturalistic free‐ranging conditions. Instead, assortative reproduction was driven by male mate preferences and, to a lesser extent, male–male competition. The effects were more pronounced when the habitat structure promoted high lizard densities. These results emphasize that assortative reproduction can occur in the absence of female choice and that male behaviour may play an important role in limiting hybridization during the initial stages of secondary contact.     

Phylogeography and pleistocene evolution in the North American black bear

To determine the extent of phylogeographic structuring in North American black bear (Ursus americanus) populations, we examined mitochondrial DNA sequences (n = 118) and restriction fragment length polymorphism profiles (n = 258) in individuals from 16 localities. Among the bears examined, 19 lineages falling into two highly divergent clades were identified. The clades differ at 5.0% of nucleotide positions, a distance consistent with an origin 1.8 MYA, and have different but overlapping geographical distributions. Areas of clade cooccurrence show that eastern and western populations are currently mixing, but regional differences in lineage distribution suggest that mixing has begun only recently. The long-term population history of black bears appears to be characterized predominantly by long-term regional isolation followed by recent contact and hybridization. Congruence between the pattern of diversity observed in black bears and patterns of forest refuge formation during the Pleistocene supports earlier speculation that Pleistocene forest fragmentations underlie a common pattern in the phylogeography of North American forest taxa.      

Spatial mixing of mitochondrial lineages and greater genetic diversity in some invasive populations of the American mink (<Emphasis Type="Italic">Neovison vison</Emphasis>) compared to native populations

The genetic characteristics of introduced populations have a relevant impact on their ability to establish and spread. The American mink (Neovison vison), native to North America, is an important invasive species in the Iberian Peninsula. Here, we used mitochondrial DNA sequences data to investigate the genetic diversity and phylogeographic structure of invasive versus native populations of this species. We also evaluated whether genetic diversity in invasive populations could be explained by the genetic characteristics of the native sources from which they derived. Phylogenetic analysis revealed two major lineages in the native range, which indicated a clear separation between western and eastern populations. On the contrary, we found no evidence of genetic structure in the invasive range. This was probably the result of the diverse origins of the released specimens and the rapid expansion and encounters of the introduced populations. We detected spatial mixing of both North American lineages in several sampling localities of the north central area of the Iberian Peninsula, giving rise to high levels of genetic diversity in some areas compared to North American populations. This could potentially lead to higher fitness of these individuals and thus increase the population viability and invasiveness of this species. These results point to the need to better study the populations in which lineages mix and, if necessary, intensify control efforts in them.     

Rapid genetic assimilation of native wall lizard populations (Podarcis muralis) through extensive hybridization with introduced lineages

The Common Wall Lizard (Podarcis muralis) has established more than 150 non-native populations in Central Europe, stemming from eight geographically distinct evolutionary lineages. While the majority of these introduced populations are found outside the native range, some of these populations also exist at the northern range margin in southwestern Germany. To (i) infer the level of hybridization in contact zones of alien and native lineages; and (ii) compare the genetic diversity among purebred introduced, native and hybrid populations, we used a combination of maternally inherited markers (mtDNA: cytb) and Mendelian markers (microsatellites). Our results suggest a rapid genetic assimilation of native populations by strong introgression from introduced lineages. Discordant patterns of mtDNA and nDNA variation within hybrid populations may be explained by directed mate choice of females towards males of alien lineages. In contrast to previous studies, we found a nonlinear relationship between genetic diversity and admixture level. The genetic diversity of hybrid populations was substantially higher than in introduced and native populations belonging to a single lineage, but rapidly reaching a plateau of high genetic diversity at an admixture level of two. However, even introduced populations with low founder sizes and from one source population retained moderate levels of genetic diversity and no evidence for a genetic bottleneck was found. The extent of introgression and the dominance of alien haplotypes in mixed populations indicate that introductions of non-native lineages represent a serious threat to the genetic integrity of native populations due to the rapid creation of hybrid swarms.     

Comparative Analyses of Plastid Sequences between Native and Introduced Populations of Aquatic Weeds Elodea canadensis and E. nuttallii

Non-indigenous species (NIS) are species living outside their historic or native range. Invasive NIS often cause severe environmental impacts, and may have large economical and social consequences. Elodea (Hydrocharitaceae) is a New World genus with at least five submerged aquatic angiosperm species living in fresh water environments. Our aim was to survey the geographical distribution of cpDNA haplotypes within the native and introduced ranges of invasive aquatic weeds Elodea canadensis and E. nuttallii and to reconstruct the spreading histories of these invasive species. In order to reveal informative chloroplast (cp) genome regions for phylogeographic analyses, we compared the plastid sequences of native and introduced individuals of E. canadensis. In total, we found 235 variable sites (186 SNPs, 47 indels and two inversions) between the two plastid sequences consisting of 112,193 bp and developed primers flanking the most variable genomic areas. These 29 primer pairs were used to compare the level and pattern of intraspecific variation within E. canadensis to interspecific variation between E. canadensis and E. nuttallii. Nine potentially informative primer pairs were used to analyze the phylogeographic structure of both Elodea species, based on 70 E. canadensis and 25 E. nuttallii individuals covering native and introduced distributions. On the whole, the level of variation between the two Elodea species was 53% higher than that within E. canadensis. In our phylogeographic analysis, only a single haplotype was found in the introduced range in both species. These haplotypes H1 (E. canadensis) and A (E. nuttallii) were also widespread in the native range, covering the majority of native populations analyzed. Therefore, we were not able to identify either the geographic origin of the introduced populations or test the hypothesis of single versus multiple introductions. The divergence between E. canadensis haplotypes was surprisingly high, and future research may clarify mechanisms that structure native E. canadensis populations.