Genetic admixture accelerates invasion via provisioning rapid adaptive evolution

Genetic admixture, the intraspecific hybridization among divergent introduced sources, can immediately facilitate colonization via hybrid vigor and profoundly enhance invasion via contributing novel genetic variation to adaption. As hybrid vigor is short‐lived, provisioning adaptation is anticipated to be the dominant and long‐term profit of genetic admixture, but the evidence for this is rare. We employed the 30 years' geographic‐scale invasion of the salt marsh grass, Spartina alterniflora, as an evolutionary experiment and evaluated the consequences of genetic admixture by combining the reciprocal transplant experiment with quantitative and population genetic surveys. Consistent with the documentation, we found that the invasive populations in China had multiple origins from the southern Atlantic coast and the Gulf of Mexico in the US. Interbreeding among these multiple sources generated a “hybrid swarm” that spread throughout the coast of China. In the northern and mid‐latitude China, natural selection greatly enhanced fecundity, plant height and shoot regeneration compared to the native populations. Furthermore, genetic admixture appeared to have broken the negative correlation between plant height and shoot regeneration, which was genetically‐based in the native range, and have facilitated the evolution of super competitive genotypes in the invasive range. In contrast to the evolved northern and mid‐latitude populations, the southern invasive populations showed slight increase of plant height and shoot regeneration compared to the native populations, possibly reflecting the heterotic effect of the intraspecific hybridization. Therefore, our study suggests a critical role of genetic admixture in accelerating the geographic invasion via provisioning rapid adaptive evolution.     

Introgressive hybridization and morphological transgression in the contact zone between two Mediterranean Solea species

Hybrid zones provide natural experiments where new combinations of genotypes and phenotypes are produced. Studying the reshuffling of genotypes and remodeling of phenotypes in these zones is of particular interest to document the building of reproductive isolation and the possible emergence of transgressive phenotypes that can be a source of evolutionary novelties. Here, we specifically investigate the morphological variation patterns associated with introgressive hybridization between two species of sole, Solea senegalensis and Solea aegyptiaca. The relationship between genetic composition at nuclear loci and individual body shape variation was studied in four populations sampled across the hybrid zone located in northern Tunisia. A strong correlation between genetic and phenotypic variation was observed among all individuals but not within populations, including the two most admixed ones. Morphological convergence between parental species was observed close to the contact zone. Nevertheless, the samples taken closest to the hybrid zone also displayed deviant segregation of genotypes and phenotypes, as well as transgressive phenotypes. In these samples, deviant body shape variation could be partly attributed to a reduced condition index, and the distorted genetic composition was most likely due to missing allelic combinations. These results were interpreted as an indication of hybrid breakdown, which likely contributes to postmating reproductive isolation between the two species.     

Detecting local adaptation using the joint sampling of polymorphism data in the parental and derived populations

When a local colonization in a new niche occurs, the new derived population should be subject to different selective pressures from that in the original parental population; consequently it is likely that many loci will be subject to directional selection. In such a quick adaptation event through environmental changes, it is reasonable to consider that selection utilizes genetic variations accumulated in the precolonization phase. This mode of selection from standing variation would play an important role in the evolution of new species. Here, we developed a coalescent-based simulation algorithm to generate patterns of DNA polymorphism in both parental and derived populations. Our simulations demonstrate that selection causes a drastic change in the pattern of polymorphism in the derived population, but not in the parental population. Therefore, for detecting the signature of local adaptation in polymorphism data, it is important to evaluate the data from both parental and derived populations simultaneously.     

Ecological speciation along an elevational gradient in a tropical passerine bird?

Local adaptation of populations along elevational gradients is well known, but conclusive evidence that such divergence has resulted in the origin of distinct species in parapatry remains lacking. We integrated morphological, vocal, genetic and behavioural data to test predictions pertaining to the hypothesis of parapatric ecological speciation associated with elevation in populations of a tropical montane songbird, the Grey‐breasted Wood‐wren (Henicorhina leucophrys: Troglodytidae), from the Sierra Nevada de Santa Marta, Colombia. We confirmed that two distinct populations exist along the elevational gradient. Phylogenetic analyses tentatively indicate that the two populations are not sister taxa, suggesting they did not differentiate from a single ancestor along the gradient, but rather resulted from separate colonization events. The populations showed marked divergence in morphometrics, vocalizations and genetic variation in mitochondrial and nuclear loci, and little to no evidence of hybridization. Individuals of both populations responded more strongly to their own local songs than to songs from another elevation. Although the two forms do not appear to have differentiated locally in parapatry, morphological and vocal divergence along the elevational gradient is consistent with adaptation, suggesting a possible link between adaptive evolution in morphology and songs and the origin of reproductive isolation via a behavioural barrier to gene flow. The adaptive value of phenotypic differences between populations requires additional study.     

Molecular population genetics and phenotypic diversification of two populations of the thermophilic cyanobacterium Mastigocladus laminosus

We investigated the distributions of genetic and phenotypic variation for two Yellowstone National Park populations of the heterocyst-forming cyanobacterium Mastigocladus (Fischerella) laminosus that exhibit dramatic phenotypic differences as a result of environmental differences in nitrogen availability. One population develops heterocysts and fixes nitrogen in situ in response to a deficiency of combined nitrogen in its environment, whereas the other population does neither due to the availability of a preferred nitrogen source. Slowly evolving molecular markers, including the 16S rRNA gene and the downstream internal transcribed spacer, are identical among all laboratory isolates from both populations but belie considerable genetic and phenotypic diversity. The total nucleotide diversity at six nitrogen metabolism loci was roughly three times greater than that observed for the human global population. The two populations are genetically differentiated, although variation in performance on different nitrogen sources among genotypes could not be explained by local adaptation to available nitrogen in the respective environments. Population genetic models suggest that local adaptation is mutation limited but also that the populations are expected to continue to diverge due to low migratory gene flow.     

A multilocus perspective on colonization accompanied by selection and gene flow

The colonization of novel habitats involves complex interactions between founder events, selection, and ongoing migration, and can lead to diverse evolutionary outcomes from local extinction to adaptation to speciation. Although there have been several studies of the demography of colonization of remote habitats, less is known about the demographic consequences of colonization of novel habitats within a continuous species range. Populations of the Eastern Fence Lizard, Sceloporus undulatus, are continuously distributed across two dramatic transitions in substrate color in southern New Mexico and have undergone rapid adaptation following colonization of these novel environments. Blanched forms inhabit the gypsum sand dunes of White Sands and melanic forms are found on the black basalt rocks of the Carrizozo lava flow. Each of these habitats formed within the last 10,000 years, allowing comparison of genetic signatures of population history for two independent colonizations from the same source population. We present evidence on phenotypic variation in lizard color, environmental variation in substrate color, and sequence variation for mitochondrial DNA and 19 independent nuclear loci. To confirm the influence of natural selection and gene flow in this system, we show that phenotypic variation is best explained by environmental variation and that neutral genetic variation is related to distance between populations, not partitioned by habitat. The historical demography of colonization was inferred using an Approximate Bayesian Computation (ABC) framework that incorporates known geological information and allows for ongoing migration with the source population. The inferences differed somewhat between mtDNA and nuclear markers, but overall provided strong evidence of historical size reductions in both white sand and black lava populations at the time of colonization. Populations in both novel habitats appear to have undergone partial but incomplete recovery from the initial bottleneck. Both ABC analyses and measures of mtDNA sequence diversity also suggested that population reductions were more severe in the black lava compared to the white sands habitat. Differences observed between habitats may be explained by differences in colonization time, habitat geometry, and strength or response to natural selection for substrate matching. Finally, effective population size reductions in this system appear to be more dramatic when colonization is accompanied by a change in selection regime. Our analyses are consistent with a demographic cost of adaptation to novel environments and show that it is possible to infer aspects of the historical demography of local adaptation even in the presence of ongoing gene flow.     

Sculpin hybrid zones: natural laboratories for the early stages of speciation

Firmly rooted as we are in the genomic era, it can seem incredible that as recently as 1974, Lewontin declared, 'we know virtually nothing about the genetic changes that occur in species formation'. To the contrary, we now know the genetic architecture of phenotypic differences and reproductive isolation between species for many diverse groups of plants, animals, and fungi. In recent years, detailed genetic analyses have produced a small but growing list of genes that cause reproductive isolation, several of which appear to have diverged by natural selection. Yet, a full accounting of the speciation process requires that we understand the reproductive and ecological properties of natural populations as they begin to diverge genetically, as well as the dynamics of newly evolved barriers to gene flow. One promising approach to this problem is the study of natural hybrid zones, where gene exchange between divergent populations can produce recombinant genotypes in situ . In such individuals, genomic variation might be shaped by introgression at universally adaptive or neutral loci, even as regions associated with local adaptation or reproductive isolation remain divergent. In Nolte et   al . (2009) , the authors take advantage of two independent, recently formed hybrid zones between sculpin species to investigate genome-wide patterns of reproductive isolation. Using a recently developed genomic clines method, the authors identify marker loci that are associated with isolation, and those that show evidence for adaptive introgression. Remarkably, Nolte et   al . (2009) find little similarity between the two hybrid zones in patterns of introgression, a fact that might reflect genetic variation within species or heterogeneous natural selection. In either case, their study system has the potential to provide insight into the early stages of speciation.     

Recombinant hybrids retain heterozygosity at many loci: new insights into the genomics of reproductive isolation in Populus

The maintenance of species barriers in the face of gene flow is often thought to result from strong selection against intermediate genotypes, thereby preserving genetic differentiation. Most speciation genomic studies thus aim to identify exceptionally divergent loci between populations, but divergence will be affected by many processes other than reproductive isolation (RI) and speciation. Through genomic studies of recombinant hybrids sampled in the wild, genetic variation associated with RI can be observed in situ, because selection against incompatible genotypes will leave detectable patterns of variation in the hybrid genomes. To better understand the mechanisms directly involved in RI, we investigated three natural ‘replicate’ hybrid zones between two divergent Populus species via locus‐specific patterns of ancestry across recombinant hybrid genomes. As expected, genomic patterns in hybrids and their parental species were consistent with the presence of underdominant selection at several genomic regions. Surprisingly, many loci displayed greatly increased between‐species heterozygosity in recombinant hybrids despite striking genetic differentiation between the parental genomes, the opposite of what would be expected with selection against intermediate genotypes. Only a limited, reproducible set of genotypic combinations was present in hybrid genomes across localities. In the absence of clearly delimited ‘hybrid habitats’, our results suggest that complex epistatic interactions within genomes play an important role in advanced stages of RI between these ecologically divergent forest trees. This calls for more genomic studies that test for unusual patterns of genomic ancestry in hybridizing species.     

Rapid appearance of epistasis during adaptive divergence following colonization

Theory predicts that short-term adaptation within populations depends on additive (A) genetic effects, while gene-gene interactions 'epistasis (E)' are important only in long-term evolution. However, few data exist on the genetic architecture of adaptive variation, and the relative importance of A versus non-additive genetic effects continues to be a central controversy of evolutionary biology after more than 70 years of debate. To examine this issue directly, we conducted hybridization experiments between two populations of wild soapberry bugs that have strongly differentiated in 100 or fewer generations following a host plant shift. Contrary to expectation, we found that between-population E and dominance (D) have appeared quickly in the evolution of new phenotypes. Rather than thousands of generations, adaptive gene differences between populations have evolved in tens. Such complex genetic variation could underlie the seemingly extreme rates of evolution that are increasingly reported in many taxa. In the case of the soapberry bug, extraordinary ecological opportunity, rather than mortality, may have created hard selection for genetic variants. Because ultimate division of populations into genetic species depends on epistatic loss of hybrid compatibility, local adaptation based on E may accelerate macro-evolutionary diversification.     

The genomic and ecological context of hybridization affects the probability that symmetrical incompatibilities drive hybrid speciation

Despite examples of homoploid hybrid species, theoretical work describing when, where, and how we expect homoploid hybrid speciation to occur remains relatively rare. Here, I explore the probability of homoploid hybrid speciation due to “symmetrical incompatibilities” under different selective and genetic scenarios. Through simulation, I test how genetic architecture and selection acting on traits that do not themselves generate incompatibilities interact to affect the probability that hybrids evolve symmetrical incompatibilities with their parent species. Unsurprisingly, selection against admixture at “adaptive” loci that are linked to loci that generate incompatibilities tends to reduce the probability of evolving symmetrical incompatibilities. By contrast, selection that favors admixed genotypes at adaptive loci can promote the evolution of symmetrical incompatibilities. The magnitude of these outcomes is affected by the strength of selection, aspects of genetic architecture such as linkage relationships and the linear arrangement of loci along a chromosome, and the amount of hybridization following the formation of a hybrid zone. These results highlight how understanding the nature of selection, aspects of the genetics of traits affecting fitness, and the strength of reproductive isolation between hybridizing taxa can all be used to inform when we expect to observe homoploid hybrid speciation due to symmetrical incompatibilities.     

Temporal variation can facilitate niche evolution in harsh sink environments

We examine the impact of temporal variation on adaptive evolution in "sink" environments, where a species encounters conditions outside its niche. Sink populations persist because of recurrent immigration from sources. Prior studies have highlighted the importance of demographic constraints on adaptive evolution in sinks and revealed that adaptation is less likely in harsher sinks. We examine two complementary models of population and evolutionary dynamics in sinks: a continuous-state quantitative-genetics model and an individual-based model. In the former, genetic variance is fixed; in the latter, genetic variance varies because of mutation, drift, and sampling. In both models, a population in a constant harsh sink environment can exist in alternative states: local maladaptation (phenotype comparable to immigrants from the source) or adaptation (phenotype near the local optimum). Temporal variation permits transitions between these states. We show that moderate amounts of temporal variation can facilitate adaptive evolution in sinks, permitting niche evolution, particularly for slow or autocorrelated variation. Such patterns of temporal variation may particularly pertain to sinks caused by biotic interactions (e.g., predation). Our results are relevant to the evolutionary dynamics of species' ranges, the fate of exotic invasive species, and the evolutionary emergence of infectious diseases into novel hosts.     

Genetic basis of adaptation in Arabidopsis thaliana: local adaptation at the seed dormancy QTL DOG1

Local adaptation provides an opportunity to study the genetic basis of adaptation and investigate the allelic architecture of adaptive genes. We study delay of germination 1 (DOG1), a gene controlling natural variation in seed dormancy in Arabidopsis thaliana and investigate evolution of dormancy in 41 populations distributed in four regions separated by natural barriers. Using F(ST) and Q(ST) comparisons, we compare variation at DOG1 with neutral markers and quantitative variation in seed dormancy. Patterns of genetic differentiation among populations suggest that the gene DOG1 contributes to local adaptation. Although Q(ST) for seed dormancy is not different from F(ST) for neutral markers, a correlation with variation in summer precipitation supports that seed dormancy is adaptive. We characterize dormancy variation in several F(2) -populations and show that a series of functionally distinct alleles segregate at the DOG1 locus. Theoretical models have shown that the number and effect of alleles segregatin at quantitative trait loci (QTL) have important consequences for adaptation. Our results provide support to models postulating a large number of alleles at quantitative trait loci involved in adaptation.     

Parasite‐mediated selection in a natural metapopulation of Daphnia magna

Parasite‐mediated selection varying across time and space in metapopulations is expected to result in host local adaptation and the maintenance of genetic diversity in disease‐related traits. However, nonadaptive processes like migration and extinction‐(re)colonization dynamics might interfere with adaptive evolution. Understanding how adaptive and nonadaptive processes interact to shape genetic variability in life‐history and disease‐related traits can provide important insights into their evolution in subdivided populations. Here we investigate signatures of spatially fluctuating, parasite‐mediated selection in a natural metapopulation of Daphnia magna. Host genotypes from infected and uninfected populations were genotyped at microsatellite markers, and phenotyped for life‐history and disease traits in common garden experiments. Combining phenotypic and genotypic data a Q_ST–F_ST‐like analysis was conducted to test for signatures of parasite mediated selection. We observed high variation within and among populations for phenotypic traits, but neither an indication of host local adaptation nor a cost of resistance. Infected populations have a higher gene diversity (Hs) than uninfected populations and Hs is strongly positively correlated with fitness. These results suggest a strong parasite effect on reducing population level inbreeding. We discuss how stochastic processes related to frequent extinction‐(re)colonization dynamics as well as host and parasite migration impede the evolution of resistance in the infected populations. We suggest that the genetic and phenotypic patterns of variation are a product of dynamic changes in the host gene pool caused by the interaction of colonization bottlenecks, inbreeding, immigration, hybrid vigor, rare host genotype advantage and parasitism. Our study highlights the effect of the parasite in ameliorating the negative fitness consequences caused by the high drift load in this metapopulation.     

Intraspecific hybridization and the recovery of fitness in the native legume Chamaecrista fasciculata

Genetic incompatibilities and low offspring fitness are characteristic outcomes of hybridization between species. Yet, the creative potential of recombination following hybridization continues to be debated. Here we quantify the outcome of hybridization and recombination between adaptively divergent populations of the North American legume Chamaecrista fasciculata in a large-scale field experiment. Previously, hybrids between these populations demonstrated hybrid breakdown, suggesting the expression of adaptive epistatic interactions underlying population genetic differentiation. However, the outcome of hybridization ultimately rests on the performance of even later generation recombinants. In experiments that compared the performance of recombinant F6 and F2 generations with nonrecombinant F1 and parental genotypes, we observed that increasing recombination had contrasting effects on different life-history components. Lifetime fitness, defined as the product of survivorship and reproduction, showed a strong recovery of fitness in the F6. The overall gain in fitness with increased recombination suggests that hybridization and recombination may provide the necessary genetic variation for adaptive evolution within species. We discuss the mechanisms that may account for the gain in fitness with recombination, and explore the implications for hybrid speciation and phenotypic evolution.     

Patterns of genetic diversity and candidate genes for ecological divergence in a homoploid hybrid sunflower, Helianthus anomalus

Natural hybridization accompanied by a shift in niche preference by hybrid genotypes can lead to hybrid speciation. Natural selection may cause the fixation of advantageous alleles in the ecologically diverged hybrids, and the loci experiencing selection should exhibit a reduction in allelic diversity relative to neutral loci. Here, we analyzed patterns of genetic diversity at 59 microsatellite loci associated with expressed sequence tags (ESTs) in a homoploid hybrid sunflower species, Helianthus anomalus. We used two indices, ln RV and ln RH, to compare variation and heterozygosity (respectively) at each locus between the hybrid species and its two parental species, H. annuus and H. petiolaris. Mean values of ln RV and ln RH were significantly lower than zero, which implies that H. anomalus experienced a population bottleneck during its recent evolutionary history. After correcting for the apparent bottleneck, we found six loci with a significant reduction in variation or with heterozygosity in the hybrid species, compared to one or both of the parental species. These loci should be viewed as a ranked list of candidate loci, pending further sequencing and functional analyses. Sequence data were generated for two of the candidate loci, but population genetics tests failed to detect deviations from neutral evolution at either locus. Nonetheless, a greater than eight-fold excess of nonsynonymous substitutions was found near a putative N-myristoylation motif at the second locus (HT998), and likelihood-based models indicated that the protein has been under selection in H. anomalus in the past and, perhaps, in one or both parental species. Finally, our data suggest that selective sweeps may have united populations of H. anomalus isolated by a mountain range, indicating that even low gene-flow species may be held together by the spread of advantageous alleles.     

A powerful regression-based method for admixture mapping of isolation across the genome of hybrids

We propose a novel method that uses natural admixture between divergent lineages (hybridization) to investigate the genetic architecture of reproductive isolation and adaptive introgression. Our method employs multinomial regression to estimate genomic clines and to quantify introgression for individual loci relative to the genomic background (clines in genotype frequency along a genomic admixture gradient). Loci with patterns of introgression that deviate significantly from null expectations based on the remainder of the genome are potentially subject to selection and thus of interest to understanding adaptation and the evolution of reproductive isolation. Using simulations, we show that different forms of selection modify these genomic clines in predictable ways and that our method has good power to detect moderate to strong selection for multiple forms of selection. Using individual-based simulations, we demonstrate that our method generally has a low false positive rate, except when genetic drift is particularly pronounced (e.g. low population size, low migration rates from parental populations, and substantial time since initial admixture). Additional individual-based simulations reveal that moderate selection against heterozygotes can be detected as much as 50 c m away from the focal locus directly experiencing selection, but is not detected at unlinked loci. Finally, we apply our analytical method to previously published data sets from a mouse ( Mus musculus and M. domesticus ) and two sunflower ( Helianthus petiolaris and H. annuus ) hybrid zones. This method should be applicable to numerous species that are currently the focus of research in evolution and ecology and should help bring about new insights regarding the processes underlying the origin and maintenance of biological diversity.     

Adaptive Introgression Facilitates Adaptation to High Latitudes in European Aspen (Populus tremula L.)

Understanding local adaptation has become a key research area given the ongoing climate challenge and the concomitant requirement to conserve genetic resources. Perennial plants, such as forest trees, are good models to study local adaptation given their wide geographic distribution, largely outcrossing mating systems, and demographic histories. We evaluated signatures of local adaptation in European aspen (Populus tremula) across Europe by means of whole-genome resequencing of a collection of 411 individual trees. We dissected admixture patterns between aspen lineages and observed a strong genomic mosaicism in Scandinavian trees, evidencing different colonization trajectories into the peninsula from Russia, Central and Western Europe. As a consequence of the secondary contacts between populations after the last glacial maximum, we detected an adaptive introgression event in a genome region of ∼500 kb in chromosome 10, harboring a large-effect locus that has previously been shown to contribute to adaptation to the short growing seasons characteristic of Northern Scandinavia. Demographic simulations and ancestry inference suggest an Eastern origin—probably Russian—of the adaptive Nordic allele which nowadays is present in a homozygous state at the north of Scandinavia. The strength of introgression and positive selection signatures in this region is a unique feature in the genome. Furthermore, we detected signals of balancing selection, shared across regional populations, that highlight the importance of standing variation as a primary source of alleles that facilitate local adaptation. Our results, therefore, emphasize the importance of migration–selection balance underlying the genetic architecture of key adaptive quantitative traits.     

AFLPs and Mitochondrial Haplotypes Reveal Local Adaptation to Extreme Thermal Environments in a Freshwater Gastropod

The way environmental variation shapes neutral and adaptive genetic variation in natural populations is a key issue in evolutionary biology. Genome scans allow the identification of the genetic basis of local adaptation without previous knowledge of genetic variation or traits under selection. Candidate loci for divergent adaptation are expected to show higher F_ST than neutral loci influenced solely by random genetic drift, migration and mutation. The comparison of spatial patterns of neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection among populations living in contrasting environments. Using the gastropod Radix balthica as a system, we analyzed 376 AFLP markers and 25 mtDNA COI haplotypes for candidate loci and associations with local adaptation among contrasting thermal environments in Lake Mývatn, a volcanic lake in northern Iceland. We found that 2% of the analysed AFLP markers were under directional selection and 12% of the mitochondrial haplotypes correlated with differing thermal habitats. The genetic networks were concordant for AFLP markers and mitochondrial haplotypes, depicting distinct topologies at neutral and candidate loci. Neutral topologies were characterized by intense gene flow revealed by dense nets with edges connecting contrasting thermal habitats, whereas the connections at candidate loci were mostly restricted to populations within each thermal habitat and the number of edges decreased with temperature. Our results suggest microgeographic adaptation within Lake Mývatn and highlight the utility of genome scans in detecting adaptive divergence.     

Hybrid genome evolution by transposition

Species hybridization is reviewed focusing on its role as a source of evolutionary novelties. Contrary to the view that hybrids are lineages devoid of evolutionary value, a number of case studies are given that show how hybrids are responsible for reticulate evolution that may lead to the origin of new species. Hybrid evolution is mediated by extensive genome repatterning followed by rapid stabilization and fixation of highly adapted genotypes. Some well-documented cases demonstrate that bursts of transposition follow hybridization and may contribute to the genetic instability observed after hybridization. The mechanism that triggers transposition in hybrids is largely unknown, but coupling of hybrid transposition and demethylation has been observed in mammals and plants. A natural scenario is proposed in which marginal small hybrid populations undergo transposition mediated genome reorganizations accompanied by exogenous and endogenous selection that, in concert with drift, lead to rapid fixation of high fitness hybrid genotypes. These genotypes may represent parental introgressed species or be entirely new species.     

Epigenetic patterns newly established after interspecific hybridization in natural populations of Solanum

Interspecific hybridization is known for triggering genetic and epigenetic changes, such as modifications on DNA methylation patterns and impact on phenotypic plasticity and ecological adaptation. Wild potatoes (Solanum, section Petota) are adapted to multiple habitats along the Andes, and natural hybridizations have proven to be a common feature among species of this group. Solanum × rechei, a recently formed hybrid that grows sympatrically with the parental species S. kurtzianum and S. microdontum, represents an ideal model for studying the ecologically and evolutionary importance of hybridization in generating of epigenetic variability. Genetic and epigenetic variability and their correlation with morphological variation were investigated in wild and ex situ conserved populations of these three wild potato species using amplified fragment length polymorphism (AFLP) and methylation‐sensitive amplified polymorphism (MSAP) techniques. We observed that novel methylation patterns doubled the number of novel genetic patterns in the hybrid and that the morphological variability measured on 30 characters had a higher correlation with the epigenetic than with the genetic variability. Statistical comparison of methylation levels suggested that the interspecific hybridization induces genome demethylation in the hybrids. A Bayesian analysis of the genetic data reveled the hybrid nature of S. × rechei, with genotypes displaying high levels of admixture with the parental species, while the epigenetic information assigned S. × rechei to its own cluster with low admixture. These findings suggested that after the hybridization event, a novel epigenetic pattern was rapidly established, which might influence the phenotypic plasticity and adaptation of the hybrid to new environments.