Reconstructing patterns of reticulate evolution in plants

Until recently, rigorously reconstructing the many hybrid speciation events in plants has not been practical because of the limited number of molecular markers available for plant phylogenetic reconstruction and the lack of good, biologically based methods for inferring reticulation (network) events. This situation should change rapidly with the development of multiple nuclear markers for phylogenetic reconstruction and new methods for reconstructing reticulate evolution. These developments will necessitate a much greater incorporation of population genetics into phylogenetic reconstruction than has been common. Population genetic events such as gene duplication coupled with lineage sorting and meiotic and sexual recombination have always had the potential to affect phylogenetic inference. For tree reconstruction, these problems are usually minimized by using uniparental markers and nuclear markers that undergo rapid concerted evolution. Because reconstruction of reticulate speciation events will require nuclear markers that lack these characteristics, effects of population genetics on phylogenetic inference will need to be addressed directly. Current models and methods that allow hybrid speciation to be detected and reconstructed are discussed, with a focus on how lineage sorting and meiotic and sexual recombination affect network reconstruction. Approaches that would allow inference of phylogenetic networks in their presence are suggested.     

A tale of two genomes: contrasting patterns of phylogeographic structure in a widely distributed bat

One of the most widely distributed bats in the New World, the big brown bat (Eptesicus fuscus) exhibits well-documented geographic variation in morphology and life history traits, suggesting the potential for significant phylogeographic structure as well as adaptive differentiation among populations. In a pattern broadly consistent with morphologically defined subspecies, we found deeply divergent mitochondrial lineages restricted to different geographic regions. In contrast, sequence data from two nuclear loci suggest a general lack of regional genetic structure except for peripheral populations in the Caribbean and Mexico/South America. Coalescent analyses suggest that the striking difference in population structure between genomes cannot be attributed solely to different rates of lineage sorting, but is likely due to male-mediated gene flow homogenizing nuclear genetic diversity across most of the continental range. Despite this ongoing gene flow, selection has apparently been effective in producing and maintaining adaptive differentiation among populations, while strong female site fidelity, maintained over the course of millions of years, has produced remarkably deep divergence among geographically isolated matrilines. Our results highlight the importance of evaluating multiple genetic markers for a more complete understanding of population structure and history.     

Inferring phylogeny despite incomplete lineage sorting

It is now well known that incomplete lineage sorting can cause serious difficulties for phylogenetic inference, but little attention has been paid to methods that attempt to overcome these difficulties by explicitly considering the processes that produce them. Here we explore approaches to phylogenetic inference designed to consider retention and sorting of ancestral polymorphism. We examine how the reconstructability of a species (or population) phylogeny is affected by (a) the number of loci used to estimate the phylogeny and (b) the number of individuals sampled per species. Even in difficult cases with considerable incomplete lineage sorting (times between divergences less than 1 N(e) generations), we found the reconstructed species trees matched the "true" species trees in at least three out of five partitions, as long as a reasonable number of individuals per species were sampled. We also studied the tradeoff between sampling more loci versus more individuals. Although increasing the number of loci gives more accurate trees for a given sampling effort with deeper species trees (e.g., total depth of 10 N(e) generations), sampling more individuals often gives better results than sampling more loci with shallower species trees (e.g., depth = 1 N(e)). Taken together, these results demonstrate that gene sequences retain enough signal to achieve an accurate estimate of phylogeny despite widespread incomplete lineage sorting. Continued improvement in our methods to reconstruct phylogeny near the species level will require a shift to a compound model that considers not only nucleotide or character state substitutions, but also the population genetics processes of lineage sorting. [Coalescence; divergence; population; speciation.].     

Parentage and sibship inference from markers in polyploids

Many plants and some animal species are polyploids. Nondisomically inherited markers (e.g. microsatellites) in such species cannot be analysed directly by standard population genetics methods developed for diploid species. One solution is to transform the polyploid codominant genotypes to pseudodiploid‐dominant genotypes, which can then be analysed by standard methods for various purposes such as spatial genetic structure, individual relatedness and relationship. Although this data transformation approach has been used repeatedly in the literature, no systematic study has been conducted to investigate how efficient it is, how much marker information is lost and thus how much analysis accuracy is reduced. More specifically, it is unknown whether or not the transformed data can be used to infer parentage and sibship jointly, and how different sampling schemes (number and polymorphism of markers, number of individuals) and ploidy level affect the inference accuracy. This study analyses both simulated and empirical data to examine the effects of polyploid levels, actual pedigree structures and marker number and polymorphism on the accuracy of joint parentage and sibship assignments in polyploid species. We show that sibship, parentage and selfing rates in polyploids can be inferred accurately from a typical set of microsatellite loci. We also show that inferences can be substantially improved by allowing for a small genotyping error rate to accommodate the distortion in assumed Mendelian inheritance of the converted markers when large sibship groups are involved. The results are discussed in the context of polyploid data analysis in molecular ecology.     

Characterizing genic and nongenic molecular markers: comparison of microsatellites and SNPs

The implications of transitioning to single nucleotide polymorphism (SNPs) from microsatellite markers (MSs) have been investigated in a number of population genetics studies, but the effect of genomic location on the amount of information each type of marker reveals has not been explored in detail. We developed novel SNP markers flanking 1 kb regions of 13 genic (within gene or <1 kb away from gene) and 13 nongenic (>10 kb from annotated gene) MSs in the threespine stickleback genome to obtain comparable data for both types of markers. We analysed patterns of genetic diversity and divergence on various geographic scales after converting the SNP loci within each genomic region into haplotypes. Marker type (SNP haplotype or MS) and location (genic or nongenic) significantly affected most estimates of population diversity and divergence. Between‐lineage divergence was significantly higher in SNP haplotypes (genic and nongenic), however, within‐lineage divergence was similar between marker types. Most divergence and diversity measures were uncorrelated between markers, except for population differentiation which was correlated between MSs and SNP haplotypes (both genic and nongenic). Broad‐scale population structure and assignment were similarly resolved by both marker types, however, only the MSs were able to delimit fine‐scale population structuring, particularly when genic and nongenic markers were combined. These results demonstrate that estimates of genetic variability and differentiation among populations can be strongly influenced by marker type, their genomic location in relation to genes and by the interaction of these two factors. This highlights the importance of having an awareness of the inherent strengths and limitations associated with different molecular tools to select the most appropriate methods for accurately addressing various ecological and evolutionary questions.     

Incomplete lineage sorting impacts the inference of macroevolutionary regimes from molecular phylogenies when concatenation is employed: An analysis based on Cetacea

Interest in methods that estimate speciation and extinction rates from molecular phylogenies has increased over the last decade. The application of such methods requires reliable estimates of tree topology and node ages, which are frequently obtained using standard phylogenetic inference combining concatenated loci and molecular dating. However, this practice disregards population‐level processes that generate gene tree/species tree discordance. We evaluated the impact of employing concatenation and coalescent‐based phylogeny inference in recovering the correct macroevolutionary regime using simulated data based on the well‐established diversification rate shift of delphinids in Cetacea. We found that under scenarios of strong incomplete lineage sorting, macroevolutionary analysis of phylogenies inferred by concatenating loci failed to recover the delphinid diversification shift, while the coalescent‐based tree consistently retrieved the correct rate regime. We suggest that ignoring microevolutionary processes reduces the power of methods that estimate macroevolutionary regimes from molecular data.     

Discordance between phylogenetics and coalescent-based divergence modelling: exploring phylogeographic patterns of speciation in the Carex macrocephala species complex

We fit a molecular data set, consisting of the rpL16 cpDNA marker and eight microsatellite loci, to the isolation-with-migration model as implemented in IM a to test a well-supported phylogenetic hypothesis of relationships within the Carex macrocephala species complex (Cyperaceae). The phylogenetic hypothesis suggests C. macrocephala from North America is reciprocally monophyletic and is sister to a reciprocally monophyletic clade of C. kobomugi . The North American C. macrocephala and C. kobomugi clade form a sister clade with a lineage of Asian C. macrocephala , thereby forming a paraphyletic C. macrocephala species. Not only does the phylogenetic hypothesis suggest C. macrocephala is paraphyletic, but it also suggests that the two lineages which share a partially overlapping distribution, Asian C. macrocephala and C. kobomugi , are not the most closely related. To test these relationships, we used coalescent-based population genetic models to infer divergence time for each lineage pair within the species complex. The coalescent-based models account for the stochastic forces which drive population divergence, and can account for the lineage sorting that occurs prior to lineage divergence. A drawback to phylogenetic-based phylogeographical analyses is that they do not account for stochastic lineage sorting that occurs between gene divergence and lineage divergence. By comparing the relative divergence time of the three main lineages within this group, Asian C. macrocephala , North American C. macrocephala , and C. kobomugi , we concluded that the phylogenetic hypothesis is incorrect, and the divergence between these lineages occurred during the Late Pleistocene epoch.     

Distinguishing adaptive from nonadaptive genetic differentiation: comparison of Q(ST) and F(ST) at two spatial scales

Genetic differentiation in 20 hierarchically sampled populations of wild barley was analyzed with quantitative traits, allozymes and Random Amplified Polymorphic DNAs (RAPDs), and compared for three marker types at two hierarchical levels. Regional subdivision for both molecular markers was much lower than for quantitative traits. For both allozymes and RAPDs, most loci exhibited minor or no regional differentiation, and the relatively high overall estimates of the latter were due to several loci with exceptionally high regional differentiation. The allozyme- and RAPD-specific patterns of differentiation were concordant in general with one another, but not with quantitative trait differentiation. Divergent selection on quantitative traits inferred from very high regional Q(ST) was in full agreement with our previous results obtained from a test of local adaptation and multilevel selection analysis. In contrast, most variation in allozyme and RAPD variation was neutral, although several allozyme loci and RAPD markers were exceptional in their levels of regional differentiation. However, it is not possible to answer the question whether these exceptional loci are directly involved in the response to selection pressure or merely linked to the selected loci. The fact that Q(ST) and F(ST) did not differ at the population scale, that is, within regions, but differed at the regional scale, for which local adaptation has been previously shown, implies that comparison of the level of subdivision in quantitative traits, as compared with molecular markers, is indicative of adaptive population differentiation only when sampling is carried out at the appropriate scale.     

Mapping the genomic architecture of adaptive traits with interspecific introgressive origin: a coalescent-based approach

Recent studies of eukaryotes including human and Neandertal, mice, and butterflies have highlighted the major role that interspecific introgression has played in adaptive trait evolution. A common question arises in each case: what is the genomic architecture of the introgressed traits? One common approach that can be used to address this question is association mapping, which looks for genotypic markers that have significant statistical association with a trait. It is well understood that sample relatedness can be a confounding factor in association mapping studies if not properly accounted for. Introgression and other evolutionary processes (e.g., incomplete lineage sorting) typically introduce variation among local genealogies, which can also differ from global sample structure measured across all genomic loci. In contrast, state-of-the-art association mapping methods assume fixed sample relatedness across the genome, which can lead to spurious inference. We therefore propose a new association mapping method called Coal-Map, which uses coalescent-based models to capture local genealogical variation alongside global sample structure. Using simulated and empirical data reflecting a range of evolutionary scenarios, we compare the performance of Coal-Map against EIGENSTRAT, a leading association mapping method in terms of its popularity, power, and type I error control. Our empirical data makes use of hundreds of mouse genomes for which adaptive interspecific introgression has recently been described. We found that Coal-Map’s performance is comparable or better than EIGENSTRAT in terms of statistical power and false positive rate. Coal-Map’s performance advantage was greatest on model conditions that most closely resembled empirically observed scenarios of adaptive introgression. These conditions had: (1) causal SNPs contained in one or a few introgressed genomic loci and (2) varying rates of gene flow – from high rates to very low rates where incomplete lineage sorting dominated as a primary cause of local genealogical variation.

     

Bayesian Inference of Reticulate Phylogenies under the Multispecies Network Coalescent

The multispecies coalescent (MSC) is a statistical framework that models how gene genealogies grow within the branches of a species tree. The field of computational phylogenetics has witnessed an explosion in the development of methods for species tree inference under MSC, owing mainly to the accumulating evidence of incomplete lineage sorting in phylogenomic analyses. However, the evolutionary history of a set of genomes, or species, could be reticulate due to the occurrence of evolutionary processes such as hybridization or horizontal gene transfer. We report on a novel method for Bayesian inference of genome and species phylogenies under the multispecies network coalescent (MSNC). This framework models gene evolution within the branches of a phylogenetic network, thus incorporating reticulate evolutionary processes, such as hybridization, in addition to incomplete lineage sorting. As phylogenetic networks with different numbers of reticulation events correspond to points of different dimensions in the space of models, we devise a reversible-jump Markov chain Monte Carlo (RJMCMC) technique for sampling the posterior distribution of phylogenetic networks under MSNC. We implemented the methods in the publicly available, open-source software package PhyloNet and studied their performance on simulated and biological data. The work extends the reach of Bayesian inference to phylogenetic networks and enables new evolutionary analyses that account for reticulation.     

Analysis of population genetic structure with RAPD markers

Recent advances in the application of the polymerase chain reaction make it possible to score individuals at a large number of loci. The RAPD (random amplified polymorphic DNA) method is one such technique that has attracted widespread interest. The analysis of population structure with RAPD data is hampered by the lack of complete genotypic information resulting from dominance, since this enhances the sampling variance associated with single loci as well as induces bias in parameter estimation. We present estimators for several population-genetic parameters (gene and genotype frequencies, within- and between-population heterozygosities, degree of inbreeding and population subdivision, and degree of individual relatedness) along with expressions for their sampling variances. Although completely unbiased estimators do not appear to be possible with RAPDs, several steps are suggested that will insure that the bias in parameter estimates is negligible. To achieve the same degree of statistical power, on the order of 2 to 10 times more individuals need to be sampled per locus when dominant markers are relied upon, as compared to codominant (RFLP, isozyme) markers. Moreover, to avoid bias in parameter estimation, the marker alleles for most of these loci should be in relatively low frequency. Due to the need for pruning loci with low-frequency null alleles, more loci also need to be sampled with RAPDs than with more conventional markers, and some problems of bias cannot be completely eliminated.     

Comparative analysis of AFLPs and SSRs efficiency in resolving population genetic structure of Mediterranean Solea vulgaris

The performance of different molecular markers in the assessment of population structure was tested using samples of Solea vulgaris collected in the Mediterranean within and outside the hypothetical dispersal ability of the species. A total of 172 individuals belonging to four population samples were analysed using 15 microsatellites [simple sequence repeats (SSRs)] and 153 amplified fragment length polymorphisms (AFLPs). Considering the global qualitative patterns, we found a correlation between SSRs and AFLPs in detecting genetic differentiation among samples. However, on a small geographical scale, AFLPs were able to discriminate individuals from neighbouring populations whereas SSRs were not, and the percentage of individuals correctly assigned to their population of origin was higher with AFLPs than with SSRs. The high number of loci analysed with the AFLP technique could increase the probability to include outlier loci in the analysis; however, the neutrality test performed on our data set did not show evidence of selection acting on the S. vulgaris samples. Even if the choice of the molecular marker depends mainly on the biological question to be addressed, the higher power of discrimination and the comparative technical ease of obtaining data from AFLPs with respect to SSRs suggest the use of AFLPs for many population genetics studies.     

Does gene flow destroy phylogenetic signal? The performance of three methods for estimating species phylogenies in the presence of gene flow

Incomplete lineage sorting has been documented across a diverse set of taxa ranging from song birds to conifers. Such patterns are expected theoretically for species characterized by certain life history characteristics (e.g. long generation times) and those influenced by certain historical demographic events (e.g. recent divergences). A number of methods to estimate the underlying species phylogeny from a set of gene trees have been proposed and shown to be effective when incomplete lineage sorting has occurred. The further effects of gene flow on those methods, however, remain to be investigated. Here, we focus on the performance of three methods of species tree inference, ESP-COAL, minimizing deep coalescence (MDC), and concatenation, when incomplete lineage sorting and gene flow jointly confound the relationship between gene and species trees. Performance was investigated using Monte Carlo coalescent simulations under four models (n-island, stepping stone, parapatric, and allopatric) and three magnitudes of gene flow (N_em = 0.01, 0.10, 1.00). Although results varied by the model and magnitude of gene flow, methods incorporating aspects of the coalescent process (ESP-COAL and MDC) performed well, with probabilities of identifying the correct species tree topology typically increasing to greater than 0.75 when five more loci are sampled. The only exceptions to that pattern included gene flow at moderate to high magnitudes under the n-island and stepping stone models. Concatenation performs poorly relative to the other methods. We extend these results to a discussion of the importance of species and population phylogenies to the fields of molecular systematics and phylogeography using an empirical example from Rhododendron.     

Distinguishing noise from signal in patterns of genomic divergence in a highly polymorphic avian radiation

Recently diverged taxa provide the opportunity to search for the genetic basis of the phenotypes that distinguish them. Genomic scans aim to identify loci that are diverged with respect to an otherwise weakly differentiated genetic background. These loci are candidates for being past targets of selection because they behave differently from the rest of the genome that has either not yet differentiated or that may cross species barriers through introgressive hybridization. Here we use a reduced‐representation genomic approach to explore divergence among six species of southern capuchino seedeaters, a group of recently radiated sympatric passerine birds in the genus Sporophila. For the first time in these taxa, we discovered a small proportion of markers that appeared differentiated among species. However, when assessing the significance of these signatures of divergence, we found that similar patterns can also be recovered from random grouping of individuals representing different species. A detailed demographic inference indicates that genetic differences among Sporophila species could be the consequence of neutral processes, which include a very large ancestral effective population size that accentuates the effects of incomplete lineage sorting. As these neutral phenomena can generate genomic scan patterns that mimic those of markers involved in speciation and phenotypic differentiation, they highlight the need for caution when ascertaining and interpreting differentiated markers between species, especially when large numbers of markers are surveyed. Our study provides new insights into the demography of the southern capuchino radiation and proposes controls to distinguish signal from noise in similar genomic scans.     

A Bayesian approach for constructing genetic maps when markers are miscoded

The advent of molecular markers has created opportunities for a better understanding of quantitative inheritance and for developing novel strategies for genetic improvement of agricultural species, using information on quantitative trait loci (QTL). A QTL analysis relies on accurate genetic marker maps. At present, most statistical methods used for map construction ignore the fact that molecular data may be read with error. Often, however, there is ambiguity about some marker genotypes. A Bayesian MCMC approach for inferences about a genetic marker map when random miscoding of genotypes occurs is presented, and simulated and real data sets are analyzed. The results suggest that unless there is strong reason to believe that genotypes are ascertained without error, the proposed approach provides more reliable inference on the genetic map.     

Variation for neutral markers is correlated with variation for quantitative traits in the plant pathogenic fungus Mycosphaerella graminicola

We compared genetic variation and population differentiation at RFLP marker loci with seven quantitative characters including fungicide resistance, temperature sensitivity, pycnidial size, pycnidial density, colony size, percentage of leaves covered by pycnidia (PLACP) and percentage of leaves covered by lesions (PLACL) in Mycosphaerella graminicola populations sampled from four regions. Wide variation in population differentiation was found across the quantitative traits assayed. Fungicide resistance, temperature sensitivity, and PLACP displayed a significantly higher Q(ST) than G(ST), consistent with selection for local adaptation, while pycnidial size, pycnidial density and colony size displayed a lower or significantly lower Q(ST) than G(ST), consistent with constraining selection. There was not a statistical difference between Q(ST) and G(ST) in PLACL. We also found a positive and significant correlation between genetic variation in molecular marker loci and quantitative traits at the multitrait scale, suggesting that estimates of overall genetic variation for quantitative traits in M. graminicola could be derived from analysis of the molecular genetic markers.     

Evaluating the Influence of the Microsatellite Marker Set on the Genetic Structure Inferred in Pyrus communis L.

Fingerprinting information can be used to elucidate in a robust manner the genetic structure of germplasm collections, allowing a more rational and fine assessment of genetic resources. Bayesian model-based approaches are nowadays majorly preferred to infer genetic structure, but it is still largely unresolved how marker sets should be built in order to obtain a robust inference. The objective was to evaluate, in Pyrus germplasm collections, the influence of the SSR marker set size on the genetic structure inferred, also evaluating the influence of the criterion used to select those markers. Inferences were performed considering an increasing number of SSR markers that ranged from just two up to 25, incorporated one at a time into the analysis. The influence of the number of SSR markers used was evaluated comparing the number of populations and the strength of the signal detected, and also the similarity of the genotype assignments to populations between analyses. In order to test if those results were influenced by the criterion used to select the SSRs, several choosing scenarios based on the discrimination power or the fixation index values of the SSRs were tested. Our results indicate that population structure could be inferred accurately once a certain SSR number threshold was reached, which depended on the underlying structure within the genotypes, but the method used to select the markers included on each set appeared not to be very relevant. The minimum number of SSRs required to provide robust structure inferences and adequate measurements of the differentiation, even when low differentiation levels exist within populations, was proved similar to that of the complete list of recommended markers for fingerprinting. When a SSR set size similar to the minimum marker sets recommended for fingerprinting it is used, only major divisions or moderate (F _ST>0.05) differentiation of the germplasm are detected.     

Phylogenetic, morphological, and chemotaxonomic incongruence in the North American endemic genus Echinacea

The study of recently formed species is important because it can help us to better understand organismal divergence and the speciation process. However, these species often present difficult challenges in the field of molecular phylogenetics because the processes that drive molecular divergence can lag behind phenotypic divergence. In the current study we show that species of the recently diverged North American endemic genus of purple coneflower, Echinacea, have low levels of molecular divergence. Data from three nuclear loci and two plastid loci provide neither resolved topologies nor congruent hypotheses about species-level relationships. This lack of phylogenetic resolution is likely due to the combined effects of incomplete lineage sorting, hybridization, and backcrossing following secondary contact. The poor resolution provided by molecular markers contrasts previous studies that found well-resolved and taxonomically supported relationships from metabolic and morphological data. These results suggest that phenotypic canalization, resulting in identifiable morphological species, has occurred rapidly within Echinacea. Conversely, molecular signals have been distorted by gene flow and incomplete lineage sorting. Here we explore the impact of natural history on the genetic organization and phylogenetic relationships of Echinacea.     

Contrasting patterns of genetic structure in two species of the coral trout Plectropomus (Serranidae) from east and west Australia: introgressive hybridisation or ancestral polymorphisms

Inter-specific genetic relationships among regional populations of two species of grouper (Plectropomus maculatus and Plectropomus leopardus) were examined using mitochondrial and nuclear markers. mtDNA revealed contrasting regional inter-specific patterns whilst nuclear markers revealed contrasting patterns among markers, irrespective of region. In eastern Australia (EA) the species form a single mtDNA lineage, but the two species are reciprocally monophyletic in Western Australia (WA). This supports previous evidence for hybridisation between these species on the east coast. WA P. leopardus forms a sister relationship with the EA P. leopardus-maculatus clade while WA P. maculatus is more basal and sister to the P. leopardus lineages, indicating mtDNA does not suffer from incomplete lineage sorting for these species. In contrast, one of three nuclear markers (locus 7-90TG) differentiated the species into two reciprocally monophyletic clades, with no evidence of hybridisation or ancestral polymorphism. The remaining two nuclear markers (2-22 and ETS-2) did not separate these two species, while distinguishing other plectropomid species, suggesting incomplete lineage sorting at these nuclear loci. These results together with coalescence analyses suggest that P. leopardus females have hybridised historically with P. maculatus males and that P. maculatus mitochondria were displaced through introgressive hybridisation and fixation in the P. maculatus founder population on the Great Barrier Reef. The contrasting regional patterns of mtDNA structure may be attributed to Quaternary sea-level changes and shelf width differences driving different reef configurations on each coast. These reef configurations have provided opportunities for local scale interaction and reproduction among species on the narrower EA continental shelves, but not on the broader WA continental shelves.     

Phenotypic plasticity rather than locally adapted ecotypes allows the invasive alligator weed to colonize a wide range of habitats

Both phenotypic plasticity and locally adapted ecotypes may contribute to the success of invasive species in a wide range of habitats. Here, we conducted common garden experiments and molecular marker analysis to test the two alternative hypotheses in invasive alligator weed (Alternanthera philoxeroides), which colonizes both aquatic and terrestrial habitats. Ninety individuals from three pairs of aquatic versus terrestrial populations across southern China were analyzed, using inter simple sequence repeat (ISSR) marker, to examine population differentiation in neutral loci. Two common gardens simulating aquatic and terrestrial habitats were set up to examine population differentiation in quantitative traits. We found no evidence of population differentiation in both neutral loci and quantitative traits. Most individuals shared the same ISSR genotype. Meanwhile, plants from different habitats showed similar reaction norms across the two common gardens. In particular, plants allocated much more biomass to the belowground roots in the terrestrial environment, where alligator weed may lose part or all of the aboveground shoots because of periodical or accidental disturbances, than those in the aquatic environment. The combined evidence from molecular marker analysis and common garden experiments support the plasticity hypothesis rather than the ecotype hypothesis in explaining the adaptation of alligator weed in a wide range of habitats.