|AFLP genome scan to detect genetic structure and candidate loci under selection for local adaptation of the invasive weed Mikania micrantha.|
|Jump to Full Text|
|PMID: 22829939 Owner: NLM Status: MEDLINE|
|Why some species become successful invaders is an important issue in invasive biology. However, limited genomic resources make it very difficult for identifying candidate genes involved in invasiveness. Mikania micrantha H.B.K. (Asteraceae), one of the world's most invasive weeds, has adapted rapidly in response to novel environments since its introduction to southern China. In its genome, we expect to find outlier loci under selection for local adaptation, critical to dissecting the molecular mechanisms of invasiveness. An explorative amplified fragment length polymorphism (AFLP) genome scan was used to detect candidate loci under selection in 28 M. micrantha populations across its entire introduced range in southern China. We also estimated population genetic parameters, bottleneck signatures, and linkage disequilibrium. In binary characters, such as presence or absence of AFLP bands, if all four character combinations are present, it is referred to as a character incompatibility. Since character incompatibility is deemed to be rare in populations with extensive asexual reproduction, a character incompatibility analysis was also performed in order to infer the predominant mating system in the introduced M. micrantha populations. Out of 483 AFLP loci examined using stringent significance criteria, 14 highly credible outlier loci were identified by Dfdist and Bayescan. Moreover, remarkable genetic variation, multiple introductions, substantial bottlenecks and character compatibility were found to occur in M. micrantha. Thus local adaptation at the genome level indeed exists in M. micrantha, and may represent a major evolutionary mechanism of successful invasion. Interactions between genetic diversity, multiple introductions, and reproductive modes contribute to increase the capacity of adaptive evolution.|
|Ting Wang; Guopei Chen; Qijie Zan; Chunbo Wang; Ying-juan Su|
Related Documents :
|22677909 - The peptide nucleic acid-locked nucleic acid polymerase chain reaction clamp-based test...
22410809 - Association of ppar-γ gene polymorphisms with obesity and obesity-associated phenotype...
19653989 - Head morphometry and isoenzymatic profile of two triatoma infestans klůg, 1834 (hemipt...
12406239 - The genetic structure of the rare lagoonal sea anemone, nematostella vectensis stephens...
22703879 - Secondary variants in individuals undergoing exome sequencing: screening of 572 individ...
22334579 - Genome-wide selection on codon usage at the population level in the fungal model organi...
22677909 - The peptide nucleic acid-locked nucleic acid polymerase chain reaction clamp-based test...
18469049 - Genetic variation of plasma insulin-like growth factor-1 in young crossbred ewes and it...
18760859 - Impact of egfr mutation analysis in non-small cell lung cancer.
|Type: Journal Article; Research Support, Non-U.S. Gov't Date: 2012-07-19|
|Title: PloS one Volume: 7 ISSN: 1932-6203 ISO Abbreviation: PLoS ONE Publication Date: 2012|
|Created Date: 2012-07-25 Completed Date: 2012-12-07 Revised Date: 2013-07-12|
Medline Journal Info:
|Nlm Unique ID: 101285081 Medline TA: PLoS One Country: United States|
|Languages: eng Pagination: e41310 Citation Subset: IM|
|CAS Key Laboratory of Plant Germplasm Enhancement and Specialty Agriculture, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan, China.|
|APA/MLA Format Download EndNote Download BibTex|
Amplified Fragment Length Polymorphism Analysis / methods*
Mikania / genetics*
Plant Weeds / genetics*
Journal ID (nlm-ta): PLoS One
Journal ID (iso-abbrev): PLoS ONE
Journal ID (publisher-id): plos
Journal ID (pmc): plosone
Publisher: Public Library of Science, San Francisco, USA
Wang et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Received Day: 30 Month: 11 Year: 2011
Accepted Day: 25 Month: 6 Year: 2012
collection publication date: Year: 2012
Electronic publication date: Day: 19 Month: 7 Year: 2012
Volume: 7 Issue: 7
E-location ID: e41310
PubMed Id: 22829939
Publisher Id: PONE-D-11-24197
|AFLP Genome Scan to Detect Genetic Structure and Candidate Loci under Selection for Local Adaptation of the Invasive Weed Mikania micrantha Alternate Title:AFLP Genome Scan to Invasive Mikania micrantha|
|Dmitry A. Filatovedit1||
1CAS Key Laboratory of Plant Germplasm Enhancement and Specialty Agriculture, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan, China
2State Key Laboratory of Biocontrol, School of Life Sciences, Sun Yat-sen University, Guangzhou, China
3Shenzhen Wildlife Rescue and Rehabilitation Center, Shenzhen, China
|University of Oxford, United Kingdom
|Correspondence: * E-mail: email@example.com
Contributed by footnote: Conceived and designed the experiments: TW YS. Performed the experiments: TW GC QZ YS. Analyzed the data: TW YS. Wrote the paper: TW YS CW.
Why some species become successful invaders is an important issue in invasive biology. When species are introduced into a new region, they face two fates. Some species quickly go extinct, whereas others persist and finally become highly competitive invaders, posing a serious threat to native diversity and ecosystems , . Successful invasions involve three major phases: introduction, naturalization, and invasion , . In the initial introduction phase, invasive species often contain low levels of genetic diversity due to bottleneck and founder effects –. Then, invaders produce pre-adapted genotypes in response to the abrupt environmental changes during naturalization –. Finally, exotic species become broadly invasive in the extended period , . The rapid population expansion of invaders is expected to promote adaptive evolution, since it has been shown that the rapidly increasing population size is conducive to withstanding (and responding to) strong directional selection , . A substantial time lag is involved during the transition from introduction via naturalization to invasion . The occurrence of a lag phase allows populations to adapt to new environmental factors such as ecological niche, temperature, precipitation, soils, frost, and wind speed or growing season length , , –. On the other hand, neutral or deleterious alleles, which become favored in new ecological contexts, will contribute to adaptive changes of invasive populations . These changes may increase the survival rate of invasive species , making them become gradually dominant in the introduced range . Therefore, pre-adaptation to novel environments is often counted as a premise for successful invasion , , .
Genomic scans are useful to identify potential adaptive loci under selection at the genomic level , . All loci across the genome are anticipated to possess similar demography and neutral evolution history of populations, including genetic drift and gene dispersal . If variation of a locus is beyond the genomic pattern with an unusual frame of higher genetic differentiation, it is deemed an “outlier locus” under natural selection , . The outlier locus can be identified explicitly in the genes under selection and also in neutral flanking regions due to hitchhiking effects , . In model organisms for which whole genomic information is available, it is easy to track the “outlier locus” under selection . However, for non-model organism such as invasive species, it becomes difficult to identify candidate genes and pinpoint the evolutionary and genetic factors involved in invasiveness because of restricted genomic resources , . New methods, especially those based on the polymerase chain reaction (PCR) to obtain amplified polymorphisms, have frequently been used to scan genomes in non-model organisms , . Among these, amplified fragment length polymorphisms (AFLPs) are reported to be the most efficient approach to identify candidate genomic regions under selection (candidate loci or outlier loci) , , –. They can provide several hundred random loci scattered throughout the genome at less cost . So far, the most commonly used analysis approaches for AFLP genome scans are Dfdist and the hierarchical Bayesian method (Bayescan). Dfdist, which was originally developed by Beaumont & Nichols , employs a classical Wright's island model to generate the expected neutral distribution of FST estimates –. In contrast, Bayescan relies on a logistic regression model , representing an extended method devised by Beaumont & Balding . To decrease false positives, Dfdist mainly depends on the trimmed mean of the empirical FST distribution for the simulation , whereas Bayescan uses a likelihood ratio test to evaluate the most likely of two alternative models, one that includes the effect of selection and another that excludes it , . Bayescan has been suggested to be more efficient at detecting high selective loci with low false positives .
Mikania micrantha H. B. K. (Asteraceae) is one of the top ten worst weeds in the world . It is a many-branched climbing perennial vine native to tropical Central and South America where it is a weed of minor importance –. M. micrantha was deliberately introduced into Asia as early as 1900s . In humid sub-tropical China, M. micrantha is subjected to new selective pressures, including such abiotic and biotic stresses as new climate, soils, pathogens, herbivores, pollinators, and competitors . By adaptive evolutionary changes, M. micrantha survived and reproduced in new environments using strikingly different strategies from those employed in its original niche . It demonstrated the ability to outcompete native plants in utilizing such limited resources as soil nutrients and sunlight and releasing phytotoxic compounds to inhibit growth of neighboring plants –. Once established, M. micrantha became a dominant plant across most regions, influencing ecosystems, biological diversity, and natural communities . By the late 1980s and early 1990s, it had spread extensively in southern China, colonizing agricultural land, orchards, nurseries, lawns, mangroves, secondary forests, scrubland, waste ground, ponds, and seashore , . Unlike in its native range, in southern China M. micrantha grows on dry soils as well as shady sites . Its favorable growth conditions have changed to an average annual temperature higher than 21°C and soil moisture content over 15% . Since spreading at accelerating rates, it is reasonable to postulate that adaptive evolution provides a key mechanism allowing the success of M. micrantha in new environments , . Such adaptive processes may leave specific signatures in the genome of M. micrantha. Thus we expected to detect the “outlier locus” signature of such local adaptations under selection at the genomic level. It is particular important to find the “outlier locus” associated with the local adaptation because, firstly, the information is critical to understanding the introduction history and genetic consequences of introduced populations of M. micrantha, . Second, it may provide fresh insights into the evolutionary potential of M. micrantha populations, which is helpful to predict their adaptability in response to management practices , . Finally, the locus will offer a prime candidate for functional surveys targeting the linked gene and dissecting the molecular mechanism of invasiveness .
To test our prediction that outlier loci can be found in southern China populations of M. micrantha (Figure 1; Table 1) and to detect genetic structure, in the present study we performed a genome scan analysis based on a large number of AFLP polymorphisms. The goals of the study were (1) to identify candidate loci under selection for local adaptation in M. micrantha and (2) to understand contributions of interactions between genetic diversity, reproductive modes, bottlenecks, and multiple introductions to the adaptive evolution of M. micrantha.
AFLP analysis of M. micrantha was optimized by determining which selective primer sets produced the clearest DNA fragments based on the eight primer pair combinations. The eight primer pairs were polymorphic in each population. A total of 483 clear polymorphic fragments were detected. The proportion of polymorphic loci within a population ranged from 10.57% to 73.47%, with an average value of 41.84% (Table 2). The average number of fragments per individual was 338.39, and 67.75 per primer combination.
At the regional level, either Hong Kong or Zhuhai was consistently found to maintain the highest level of variation according to all the genetic diversity statistics (Table 2). However, populations such as MA106, SZ72, NLD10, NLD26, and NLD30 growing in regions of Macao, Shenzhen, and Neilingding were observed to possess lower amounts of genetic variation than others in the same region (Table 2).
We estimated levels of linkage disequilibrium (LD) and character compatibility among the 483 AFLP loci. Each region was observed to possess significant LD and matrix compatibility obtained by testing the index of association (IA and ) and the incompatibility excess ratio (IER) (Table 2). Likewise, significant LD and matrix compatibility were also detected at the population scale (Table 2).
Alleles were defined as rare when they occurred at a frequency of less than 0.05 in the examined populations . In total, 22 rare alleles were identified at the species level. Across regions, 14, 12, and 9 rare alleles were detected in Hong Kong, Shenzhen, and Neilingding, respectively, whereas no rare alleles were found in Macao, Zhuhai, or Dongguan. In addition, no rare allele was detected within any populations at the population level.
Using the stepwise mutation model (SMM) and the infinite allele model (IAM), we identified bottleneck signatures in each population with a heterozygosity excess/deficiency ratio that significantly deviated from the expected ratio (1∶1) at mutation-drift equilibrium (Table 3, P<0.05). Moreover, the entire range of M. micrantha appeared subject to a significant bottleneck (Table 3).
The genetic differentiation ΦST measured by AMOVA analysis was 0.3335; 66.65% of the variation was partitioned within populations, 28.61% was attributed to differences among populations within regions, and only 4.74% of the variation was due to regional differences (all three hierarchical levels were significant with P<0.001) (Table 4). When individual pairs of populations were compared, of the 378 pairwise ΦST values, 373 were significant (P<0.001), whereas only five values derived from SZ64 vs. SZ75, SZ64 vs. SZ77, SZ64 vs. DG91, SZ75 vs. SZ77, and SZ75 vs. DG91 were not (P>0.05), highlighting remarkable differences among populations.
The f-free model chosen based on the smallest mean DIC was more suitable than other models for the AFLP data set of the 28 populations. Therefore, θB = 0.2927 (95% credible interval: 0.2772–0.3058) was determined to be an unbiased Bayesian estimate of all population genetic differentiation.
No significant associations between geographical distances and pairwise estimates of θB were observed in the full AFLP data set (r = 0.0032, P = 0.452). A similar pattern was derived at the regional level. Hence no evidence of “isolation by distance”  was revealed by Mantel tests at either the regional or entire range level.
The ΔK criterion of Evanno et al. was applied to estimate the number of population clusters. The maximal value of ΔK was K = 24, indicating the need to divide the samples into 24 clusters. Bar plots showed varying extents of admixture among populations (Figure 2). Ten clusters including populations NLD10, NLD26, NLD30, NLD31, SZ35, SZ72, HK87, HK90, MA105, and MA106 were detected, although each population also contained a very few individuals from different regions. The substructure of the rest populations was weak and could only be resolved when information about the sampling locations was included. A mix of individuals from other clusters was often observed in each cluster. The results suggested that M. micrantha populations in southern China maintained a relatively weak genetic structure.
In addition, the UPGMA tree rooted with M. cordata revealed that populations of M. micrantha from different locations usually were intermingled (Figure 3), suggesting a lack of geographic pattern. HK87, which was documented as the earliest introduced population in China , formed an independent branch, one that first diverged from the other M. micrantha populations.
The AFLP data set of M. micrantha was used in a global analysis for outlier detection, both with Dfdist and Bayescan (Table 5). By using Dfdist on the 28 populations of M. micrantha, 23 out of 483 loci (4.8%) were identified as outlier loci under directional selection at the 99.5% probability level (Figure 4A). Bayescan analysis produced high differentiation loci at a threshold of log10 PO>2.0 (posterior probabilities higher than 0.99), corresponding to 7.87% of the 483 investigated loci. This approach identified 38 outlier loci potentially under selection or linked to a locus under selection (Table 5; Figure 5A). Notably, Bayescan found 24 outliers (loci 4, 5, 11, 14, 43, 54, 59, 153, 166, 167, 173, 209, 224, 228, 275, 276, 388, 420, 425, 434, 451, 455, 473, and 483) that were not detected by Dfdist. By contrast, Dfdist detected nine outliers (loci 35, 50, 56, 136, 231, 240, 324, 342, and 465) not identified by Bayescan. By pooling the results of the totally different two detection approaches, 14 outlier loci (3, 17, 25, 57, 116, 132, 177, 183, 184, 219, 325, 347, 381 and 426) were identified by Dfdist and Bayescan. The 14 outliers represented truly adaptive loci (not false positives) because the two approaches used different algorithms and assumptions, and very stringent significance criteria were considered (99.5% probability level for Dfdist, posterior probability >0.90 for Bayescan) in the study. Additionally, in the two populations of M. cordata, no outlier locus under selection was detected by either Dfdist or Bayescan (Figures 4B and 5B). Nevertheless, this lack of detection of any outliers may also be due to lower sample size.
To the best of our knowledge, this study is the first report on detection of candidate loci under selection by genome scan in the invasive weed M. micrantha since its introduction to southern China. To avoid spurious outlier loci, several authors advocated employment of two or more outlier detection methods and a conservative significance level , , –. In this study, the AFLP explorative genome scan has revealed 14 loci as under selection among a total of 483 loci in M. micrantha. These 14 loci are considered to possess high credibility because they were picked up by two complementary and exhaustive methods, Dfdist and Bayescan, applying very stringent significance criteria and a 30% trimmed mean FST. Meyer et al. noted that the power of the analysis is directly associated with the genome coverage . Since a high density of AFLP loci were examined, the selective loci detected here should prove to have good reliability. Taking these together, the 14 outlier loci we identified can be considered useful for understanding the successful invasion of M. micrantha in southern China.
Currently, there is increasing interest in identifying genes or outlier loci that underlie local adaptations in invasive species –. Understanding the process of invasive adaptive evolution is critical to the introduction history and genetic consequences on introduced populations, the planning of management strategies of invasive species, and for identifying the possible risks of introduced non-native species in the future , . Because selection for local adaptation is unlikely to simultaneously occur at many loci , outlier loci are limited to a minor part of the genome, which has been suggested in other species , , –. Our results indicate that only a small part (2.9%) of the genome of M. micrantha has been under directional selection during its invasive processes. The percentage of loci detected is slightly lower than the 5–10% reported in the generality of AFLP genome scans , . As AFLP loci are likely located in non-coding DNA, some of the outlier loci may only exhibit the signature of selection because they are linked to the actual target , . Although it is difficult to know the location and function of the loci involved in the adaptation to invasiveness, a genome scan of M. micrantha still offers unique opportunity to unravel the genetic basis of invasive adaptation without known phenotypes and whole genome sequences. In particular, the AFLP primers that were developed to amplify the outlier loci identified here can be directly used to construct a reduced representation library of the M. micrantha genome, which will allow efficient sequencing of the linked genomic regions by next-generation sequencing technology .
Genetic diversity is anticipated to increase the adaptive potential of invasive populations in the new environment , . In invasive species, the alleles from standing genetic variation have been supposed to control the adaptive evolution , , because favorable alleles compared to neutral or deleterious alleles are immediately available and often occur at a greater frequency in populations , . Compared with other non-invasive and invasive weedy species (Table 6), we have detected relatively high levels of genetic diversity in the introduced populations of M. micrantha in southern China. This substantial genetic variation enables to provide a large pool of raw material for adaptive evolution , .
High genetic diversity can be created by multiple introductions, which bring together large amounts of genetic variation and novel genetic combinations . Such admixtures elevate the expansive ability and pace of proliferation of invasive species and are able to enhance adaptive evolutionary responses to novel environmental selection , , , . Based on the UPGMA and STRUCTURE analysis results, closely related populations often are from different geographical locations (Figures 2 and 3), while individuals from the same location do not group together. Thus multiple introductions are inferred among populations of M. micrantha. In addition, the significant genetic differentiation detected among M. micrantha populations is also consistent with the expectation that each population has experienced severe bottleneck (Table 3). It is generally suggested that bottlenecks tend to increase genetic differentiation among populations by changing allelic frequencies ,  and they are perceived to reduce the potential for adaptive evolution as well , . Nevertheless, multiple introductions can ameliorate this effect and accelerate the speed of adaptive evolution . Such phenomena have been recorded in the populations of other invasive species including Heracleum mantegazzianum, Hieracium lepidulum, and Alliaria petiolata. Furthermore, no geographical signature is revealed in the AFLP variation pattern among the introduced populations of M. micrantha at either regional or entire range scale. Repeated or secondary introductions, which are induced by long-distance dispersals of copious amounts of fine and fluffy seeds, pollination or propagules mediated by humans, insects, and violent winds , may further complicate the geographic consequences of the population genetic variation of M. micrantha. Therefore, multiple introductions can add to the successful invasiveness of M. micrantha by increasing genetic variation, preventing genetic bottlenecks, and strengthening the capacity of adaptive evolution.
Adaptive evolution of invasive species is also determined by the reproductive system such as the amount of sexual and asexual reproduction and the patterns of mating . Mikania micrantha mainly propagates by seeds, but its local spread mostly results from vegetative propagation , , . When reproducing vegetatively, M. micrantha generates shoots from small stem fragments and rosettes , . In this study, significant LD and locus compatibility have been identified in a majority of M. micrantha populations (Table 2), indicating that asexual reproduction is the predominant propagation system. In addition, self-fertilization (selfing) has also been reported in the southern China populations of M. micrantha. The results are in accord with the findings that most invasive plants reproduce asexually or by selfing –, and reproductive systems may be subject to evolutionary modification during invasion such as switching to higher levels of inbreeding or vegetative reproduction , , –. For species that are able to reproduce both sexually and clonally, sexual reproduction would promote allele recombinations that result in high genetic variation, whereas clonal reproduction preserves successful combinations and maintains genetic variation by occasionally creating new alleles through somatic mutations –. The specific reproductive modes (asexual and selfing) of M. micrantha are likely to enhance its adaptive potential in southern China.
In summary, this study is the first report to identify outlier loci under selection in the genome of the invasive weed M. micrantha without reference to coding vs. noncoding sequences, or specific models of selection by explorative genomic scans . We find that 2.9% candidate loci are under selection, indicating that local adaptation has indeed occurred in M. micrantha since its introduction to southern China. This study also demonstrates high genetic variability, strong genetic differentiation, severe bottlenecks, and multiple introductions, as well as primarily asexual reproduction in the introduced populations. The interactions between genetic diversity, multiple introductions, and reproductive modes are revealed to be involved in the response of M. micrantha to the new local environment, allowing a more precise understanding of its successful invasion. In future studies, we will isolate and sequence the outlier AFLP bands to identify their genomic locations and neighboring genes to further detail the adaptive molecular mechanisms of M. micrantha.
This study was conducted in accordance with all People's Republic of China laws. No specific permits were required for the described field studies. No specific permissions were required for the locations/activities described in this study. The location is not privately owned or protected in any way. The field studies did not involve endangered or protected species.
For M. micrantha, individuals were sampled from 28 populations in six regions (Hong Kong, Macao, Shenzhen, Neilingding, Zhuhai, and Dongguan) representing different habitats (Figure 1A; Table 1). In each population, we collected fresh leaf material from 12 to 19 randomly selected individuals (Table 1). These samples represented a large number of individuals of M. micrantha throughout its introduced range in southern China. To further confirm that the outlier loci under selection in M. micrantha were caused by novel environments, we tested whether AFLP loci can demonstrate evidence of selection in the congeneric species M. cordata, which is native to sub-tropical China. For M. cordata, five to 12 individuals were collected from its two populations (HN1 and HN2) located on Hainan Island (Figure 1B; Table 1). Leaves were preserved in silica gel. Voucher specimens (M. micrantha: CGP1-409; M. cordata: CGP410-426) were deposited at the herbarium of Sun Yat-sen University (SYS), Guangzhou, China.
Total genomic DNA was extracted from ground tissue, following the modified CTAB protocol . The quality and quantity of the DNA were determined with a spectrophotometer (Pharmacia 2000 UV/Visible, Amersham Pharmacia Biotech, Piscataway, NJ) and on 0.8% agarose gels.
AFLP analyses were performed according to Vos et al. with the following modifications: genomic DNA (50 ng) was digested with 2 U EcoRI and 4 U MseI (New England Biolabs, Ipswitch, MA) for 3 h at 37°C in 1×NE buffer and incubated at 70°C for 20 min. To ligate the resulting fragments to the corresponding adapters, 10 µl of restriction products were added into a 20 µl reaction mixture containing 0.1 µM EcoRI adapter, 1 µM MseI adapter, and 60 U T4 DNA ligase (New England Biolabs). After being incubated at 20°C for 3 h, the samples were diluted 10 times with ddH2O, then 2 µl of each sample was used as a template to conduct the pre-amplification in a final volume of 20 µl that contained 1×PCR buffer, 100 nM each of EcoRI+A and MseI+C primers, 0.24 mM dNTPs, and 1.5 U Taq DNA polymerase. The preamplification reaction was carried out for 20 cycles of 30 s at 94°C, 1 min at 56°C and 1 min at 72°C. After that, preamplified product was diluted 1∶50 with ddH2O, and then used as template for the selective PCR amplifications to generate AFLPs.
Selective amplification was performed in a final volume of 10 µl containing 1× PCR buffer, 125 nM EcoRI primer, 125 nM 6-FAM-EcoRI primer, 250 nM MseI primer, 0.2 mM dNTPs, 0.75 U Taq DNA polymerase and 2 µl diluted pre-amplified DNA sample. The reaction was conducted for 13 cycles of 30 s at 94°C, 30 s at 65°C and 1 min at 72°C. The annealing temperature was reduced by 0.7°C per cycle. Then 23 cycles consisting of 30 s at 94°C, 30 s at 56°C, and 1 min at 72°C were performed.
An initial screening of 48 combinations of selective primers, in which six were EcoRI (EcoRI+AAG, +ACT, +AGC, +ACA, +AAC, +ACC) and eight were MseI (MseI+CAG, +CTC, +CTA, +CAA, +CAC, +CTG, +CAT, +CTT), were performed. Eight combinations were selected that generated clear and evenly distributed bands: EcoRI+ACT/MseI+CAG, EcoRI+AAC/MseI+CAG, EcoRI+ACT/MseI+CTA, EcoRI+AAC/MseI+CTC, EcoRI+AAC/MseI+CTA, EcoRI+ACT/MseI+CTC, EcoRI+AGC/MseI+CAC, and EcoRI+ACC/MseI+CAC. Therefore, the eight primer combinations were chosen for selective amplifications of all samples. PCR reactions were conducted on a PTC-100 Peltier Thermal Cycler (MJ Research, St. Bruno, Quebec). The selective PCR products were separated by electrophoresis on 6.5% polyacrylamide gels on an ABI 377 automated sequencer (Applied Biosystems, Carlsbad, CA). An ROX-500-labeled internal size standard (Applied Biosystems) was added to each sample to size fragments.
Software GeneScan 3.7 (Applied Biosystems) and genographer (version 1.6.0; http://hordeum.oscs.montana.edu/genographer) were utilized to collect and score raw fluorescent AFLP data. The presence (1) or absence (0) of data from unambiguous AFLP bands was used to establish the matrix of genetic identity of the sampled individuals. Genetic diversity statistics, including percentage of polymorphic loci and Nei's gene diversity , were calculated using POPGEN32 software . Shannon's index of phenotypic diversity was quantified as S = −∑ pilnpi where pi is the frequency of a given AFLP band in the population .
Analysis of Molecular Variance (AMOVA) based on a Euclidean squared distance matrix was hierarchically calculated to estimate the allocation of genetic variation among and within populations by ARLEQUIN 3.0 (available at http://cmpg.unibe.ch/software/arlequin3) . In this analysis, the AFLP data set was partitioned at three levels: regional, among- population, and within- population. One thousand random permutations were used to infer the significance of the variance components . Also using the same software, a Mantel test was performed to investigate the correlation between genetic differentiation and geographic distances (km) among populations. The matrix of genetic differentiation was composed by pairwise θB estimates.
Holsinger et al. proposed a Bayesian approach to estimating genetic structure for dominant and co-dominant markers. The nearly unbiased parameter estimations of heterozygosity, genetic distance, and population differentiation can be obtained using the Bayesian method . Its f and θB are equivalent to the inbreeding coefficient (FIS) and the fixation index (FST) of F-statistics, respectively. The posterior distributions of f and θB were estimated through Markov Chain Monte Carlo (MCMC) methods by HICKORY v1.0, with a burn-in of 50 000 iterations and a sampling run of 250 000 iterations from which every fiftieth sample was retained for posterior calculations . The analysis model was chosen based on the deviance information criterion (DIC). The “f-free” model, where f was not estimated but was chosen at random from the prior distribution, was decided due to its smaller DIC than for other models. Alleles were deemed rare if their frequencies were ≤0.05 in the sampled populations .
To gain further perspectives on genetic structure, we also employed the Bayesian clustering method to infer the pattern of population structure by implementing the STRUCTURE 2.2.3 program , . To determine the best number of clusters, 30 independent runs of K (K = 1 to 30) were performed with an admixture model at 100 000 MCMC iterations and a 10 000 burn-in period. We used ΔK, the second-order rate of change in ln P (X|K) for successive values of K to determine the number of clusters . The distribution map of STRUCTURE was plotted according to K value at the highest log likelihood. Moreover, based on Nei's genetic distance , an unweighted pair group method with arithmetic mean analysis (UPGMA) was used to generate a dendrogram of the relationships among the populations by TFPGA 1.3 . One thousand bootstrap replicates were permuted to assess the reliability of the UPGMA dendrogram.
We further tested whether populations have suffered a bottleneck using the program BOTTLENECK . The heterozygosity (Heq) expected at mutation-drift equilibrium was calculated based on both the stepwise mutation model (SMM) and the infinite allele model (IAM) , . The significance of heterozygosity excess was determined by the sign test .
Linkage disequilibrium, which is an association of alleles at different loci on chromosomes in a population, is mainly affect by population size . The dependency of population size can be removed by multilocus linkage disequilibrium with both indices of association IA– and its modified measure . These indices and their significance by randomization were computed using the software Multilocus v1.2 (http://www.bio.ic.ac.uk/evolve/software/multilocus/).
A character incompatibility analysis was carried out to probe the predominant mating system in the populations at the molecular level –. In a pair of binary character data, such as the presence or absence of AFLP bands at two loci, the presence of all four possible combinations of characters (0/0, 1/0, 0/1, 1/1) is more parsimoniously explained by sexual recombination than by three mutation events. Once all four character combinations arise, an incompatibility is said to occur, and can be used as a measure of recombination . The incompatibility excess ratio (IER) was calculated using PICA 4.0 by comparing the observed incompatibility count for the original data and mean incompatibility count for randomly permuted data . If none or a small fraction of incompatible loci occur, asexual reproduction is implied to have happened , .
To detect outlier loci under selection for local adaptation of M. micrantha and M. cordata, two complementary methods were applied to our AFLP data set. We first used Program Dfdist, which is the most popular software for detecting candidate loci . Dfdist software (http://www.rubic.rdg.ac.uk/~mab/stuff) was recently modified from Beaumont and Balding  to analyze dominant data, which is a hierarchical Bayesian approach based on summary statistics in a symmetrical island model , . Outlier loci were detected by comparing empirical FST values for each locus against a null distribution of FST values expected from a neutral drift model . A potential defect of the Dfdist approach for detecting outlier loci is the possibility of false positives –. The statistical power was enhanced to avoid false positives by setting a conservative constraint . A null distribution of FST close to the empirical distribution was acquired by 50 000 coalescent simulations. Simulations were computed with a mean FST similar to the trimmed mean FST, which was calculated by excluding 30% of the most extreme FST values observed in the empirical dataset , . The 4Nμ parameter value was set to 0.04 in all simulations. A global analysis was done for M. micrantha and M. micrantha using 28 populations and 2 populations, respectively. The threshold for outliers was set to the more conservative 0.005, estimated from simulated FST values to control for false positives , , . In addition, Dfdist assumes that populations are at migration-drift equilibrium, which does not often occur in natural populations . Hence, to cross-check the reliability of the outlier loci detected by Dfdist, we also ran Bayescan software (http://www-leca.ujf-grenoble.fr/logiciels.htm), which better handles dominant marker data by directly estimating the posterior probability of a given locus to be under selection . Assuming that allele frequencies within populations follow a Dirichlet distribution –, the Bayesian method not only permits for different demographic scenarios and different amounts of genetic drift between populations when estimating population-specific FST coefficients, but also considers all loci in the analyses . The Bayesian approach also disposes of the problem of multiple testing of a large number of genomic loci through prior distribution . In our genome scan, the log10 PO>2.0 was considered a threshold value for determining loci under selection according to Jeffreys' interpretation , which is a logarithmic scale for model choice as follows: log10 PO>0.5 (substantial); log10 PO>1.0 (strong); log10 PO>1.5 (very strong); and log10 PO>2.0 (decisive support for accepting a model) . We used 10 pilot runs of 5000 iterations to estimate model parameters. A burn-in of 50 000 iterations was employed to cover the MCMC. The sample size was set to 5000 and the thinning interval to 20, resulting in a total chain length of 150 000 iterations . The loci were ranked according to their estimated posterior probability and all loci with a value over 0.993 were retained as outliers. This corresponds to log10 PO>2.0, which provides decisive support for acceptation of the model. Outliers identified by both Dfdist and Bayescan are likely to be truly adaptive regions of the genome, because the two approaches differ in algorithms and assumptions .
Competing Interests: The authors have declared that no competing interests exist.
Funding: This work was supported by the National Natural Science Foundation of China (Grant Nos. 30970290, 31070594), the Opening Fund of Laboratory Sun Yat-sen University (2010), and Guangdong Key Laboratory of Plant Resources (plant01k13). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
We thank Zhanming Ying, Zhenyu Li, Junsheng Yu, Binghui Wei, and Mingzhao Yang of the School of Life Sciences, Sun Yat-sen University, for assistance with the collection of plant materials.
|1.||Holt RD. Year: 2009Up against the edge: invasive species as testbeds for basic questions about evolution in heterogeneous environments.Mol Ecol184347434819845861|
|2.||Xu CY,Julien MH,Fatemi M,Girod C,Van Klinken RD,et al. Year: 2010Phenotypic divergence during the invasion of Phyla canescens in Australia and France: evidence for selection-driven evolution.Ecol Lett13324419849709|
|3.||Richardson DM,Pyšek P,Rejmánek M,Barbour MG,Panetta FD,et al. Year: 2000Naturalization and invasion of alien plants: concepts and definitions.Divers Distrib693107|
|4.||Lachmuth S,Durka W,Schurr FM. Year: 2010The making of a rapid plant invader: genetic diversity and differentiation in the native and invaded range of Senecio inaequidens.Mol Ecol193952396720854275|
|5.||Barrett SCH,Husband B. Year: 1990The genetics of plant migration and colonization.Brown AHD,Clegg MT,Kahler AL,Weir BSPlant population genetics, breeding, and genetic resourcesSunderlandSinauer & Associates254277|
|6.||Novak SJ,Mack RN. Year: 1993Genetic variation in Bromus tectorum (Poaceae): comparison between native and introduced populations.Heredity71167176|
|7.||Milne RI,Abbott RJ. Year: 2000Origin and evolution of invasive naturalized material of Rhododendron ponticum L. in the British Isles.Mol Ecol954155610792698|
|8.||Dlugosch MK,Parker MI. Year: 2008Founding events in species invasions: genetic variation, adaptive evolution, and the role of multiple introductions.Mol Ecol1743144917908213|
|9.||Sakai AK,Allendorf FW,Holt JS,Lodge DM,Molofsky J,et al. Year: 2001The population biology of invasive species.Annu Rev Ecol Sys32305332|
|10.||Lee CE. Year: 2002Evolutionary genetics of invasive species.Trends Ecol Evol17386391|
|11.||Holt RD,Barfield M,Gomulkiewicz R. Year: 2005Theories of niche conservatism and evolution: could exotic species be potential tests?Sax DF,Stachowicz JJ,Gaines SDSpecies invasions: insights into ecology, evolution, and biogeographySunderlandSinauer & Associates259290|
|12.||Keller SR,Taylor DR. Year: 2008History, chance and adaptation during biological invasion: separating stochastic phenotypic evolution from response to selection.Ecol Lett1185286618422638|
|13.||Allan E,Pannell JR. Year: 2009Rapid divergence in physiological and life-history traits between northern and southern populations of the British introduced neo-species, Senecio squalidus.Oikos11810531061|
|14.||Reznick DN,Ghalambor CK. Year: 2001The population ecology of contemporary adaptations: what empirical studies reveal about the conditions that promote adaptive evolution.Genetica11218319811838765|
|15.||Whitney KD,Gabler CA. Year: 2008Rapid evolution in introduced species, ‘invasive traits’ and recipient communities: challenges for predicting invasive potential.Divers Distrib14569580|
|16.||Kaufman SR,Smouse PE. Year: 2001Comparing indigenous and introduced populations of Melaleuca quinquenervia (Cav.) Blake: response of seedlings to water and pH levels.Oecologia127487494|
|17.||Broennimann O,Treier UA,Müller-Scharer H,Thuiller W,Peterson AT,et al. Year: 2007Evidence of climatic niche shift during biological invasion.Ecol Lett1070170917594425|
|18.||Kanarek AR,Webb CT. Year: 2010Allee effects, adaptive evolution, and invasion success.Evol Appl3122135|
|19.||Barrett SCH,Colautti RI,Eckert CG. Year: 2008Plant reproductive systems and evolution during biological invasion.Mol Ecol1737338317868309|
|20.||Wang JY,Zhang HW,Huang RF. Year: 2007Expression analysis of low temperature responsive genes in Eupatorium adenophorum Spreng using cDNA-AFLP.Plant Mol Biol Rep253744|
|21.||Pyŝek P. Year: 1998Is there a taxonomic pattern to plant invasions?Oikos82282294|
|22.||Herrera CM,Bazaga P. Year: 2009Quantifying the genetic component of phenotypic variation in unpedigreed wild plants: tailoring genomic scan for within-population use.Mol Ecol182602261419457184|
|23.||Meyer CL,Vitalis R,Saumitou-Laprade P,Castric V. Year: 2009Genomic pattern of adaptive divergence in Arabidopsis halleri, a model species for tolerance to heavy metal.Mol Ecol182050206219434814|
|24.||Bonin A,Taberlet P,Miaud C,Pompanon F. Year: 2006Explorative genome scan to detect candidate loci for adaptation along a gradient of altitude in the common frog (Rana temporaria).Mol Biol Evol2377378316396915|
|25.||Murray MC,Hare MP. Year: 2006A genomic scan for divergent selection in a secondary contact zone between Atlantic and Gulf of Mexico oysters, Crassostrea virginica.Mol Ecol154229424217054515|
|26.||Schlötterer C. Year: 2003Hitchhiking mapping-functional genomics from the population genetics perspective.Trends Genet19323812493246|
|27.||Paris M,Boyer S,Bonin A,Collado A,David JP,et al. Year: 2010Genome scan in the mosquito Aedes rusticus: population structure and detection of positive selection after insecticide treatment.Mole Ecol19325337|
|28.||Galindo J,Morán P,Rolán-Alvaren E. Year: 2009Comparing geographical genetic differentiation between candidate and noncandidate loci for adaptation strengthens support for parallel ecological divergence in the marine snail Littorina saxatilis.Mol Ecol1891993019207246|
|29.||Prentis PJ,Wilson JRU,Dormontt EE,Richardson DM,Lowe AJ. Year: 2008Adaptive evolution in invasive species.Trends Plant Sci1328829418467157|
|30.||Prentis PJ,Woolfit M,Thomas-Hall SR,Ortiz-Barrientos, Pavasovic A,et al. Year: 2010Massively parallel sequencing and analysis of expressed sequence tags in a successful invasive plant.Ann Bot1061009101720929896|
|31.||Mealor BA,Hild AL. Year: 2006Potential selection in native grass populations by exotic invasion.Mol Ecol152291230016780441|
|32.||Campbell D,Bernatchez L. Year: 2004Generic scan using AFLP markers as a means to assess the role of directional selection in the divergence of sympatric whitefish ecotypes.Mol Biol Evol2194595615014172|
|33.||Gagnaire PA,Albert V,Jónsson B,Bernatchez L. Year: 2009Natural selection influences AFLP intraspecific genetic variability and introgression patterns in Atlantic eels.Mol Ecol181678169119302349|
|34.||Michalski SG,Durka W,Jentsch A,Kreyling J,Pompe S,et al. Year: 2010Evidence for genetic differentiation and divergent selection in an autotetraploid forage grass (Arrhenatherum elatius).Theor Appl Genet1201151116220049414|
|35.||Parisod C,Joost S. Year: 2010Divergent selection in trailing-versus leading-edge populations of Biscutella laevigata.Ann Bot10565566020181567|
|36.||Fischer MC,Foll M,Excoffier L,Heckel G. Year: 2011Enhanced AFLP genome scans detect local adaptation in high-altitude populations of a small rodent (Microtus arvalis).Mol Ecol201450146221352386|
|37.||Tollenaere C,Duplantier JM,Rahalison L,Ranjalahy M,Brouat C. Year: 2011AFLP genome scan in the black rat (Rattus rattus) from Madagascar: detecting genetic markers undergoing plague-mediated selection.Mol Ecol201026103820444082|
|38.||Beaumont MA,Nichols RA. Year: 1996Evaluating loci for use in the genetic analysis of population structure.Proc Roy Soc B-Biol Sci26316191626|
|39.||Wright S. Year: 1951The genetical structure of populations.Ann Eugen15323354|
|40.||Excoffier L,Hofer T,Foll M. Year: 2009Detecting loci under selection in a hierarchically structured population.Heredity10328529819623208|
|41.||Helyar SJ,Hemmer-Hansen J,Bekkevold D,Taylor MI,Ogden R,et al. Year: 2011Application of SNPs for population genetics of nonmodel organisms: new opportunities and challenges.Mol Ecol Resour11Suppl 112313621429169|
|42.||Foll M,Gaggiotti O. Year: 2008A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: a Bayesian perspective.Genetics18097799318780740|
|43.||Beaumont MA,Balding DJ. Year: 2004Identifying adaptive genetic divergence among populations from genome scans.Mol Ecol1396998015012769|
|44.||Butlin RK. Year: 2010Population genomics and speciation.Genetica13840941818777099|
|45.||Pérez-Figueroa A,García-Pereira MJ,Saura M,Rolán-Alvarez E,Caballero A. Year: 2010Comparing three different methods to detect selective loci using dominant markers.J Evol Biol232267227620796133|
|46.||Holm LG,Plucknett DL,Pancho JV,Herberger JP. Year: 1977The world's worst weeds: distribution and biologyHonoluluEast-West Center and University Press of Hawaii|
|47.||Wirjahar S. Year: 1976Autecological study of Mikania spp.Proceedings of fifth Asian-Pacific weed science society conference, 5–11, October, 1975TokyoAsian Weed Science Society7073|
|48.||Kong GH,Wu QG,Hu QM. Year: 2000Exotic weed Mikania micrantha H. B. K. appeared in south China.J Trop Subtrop Bot827|
|49.||Wang BS,Liao WB,Zan QJ,Li MG,Zhou XY,et al. Year: 2003The spreads of Mikania micrantha in China.Acta Scient Nat Univ Sunyatseni424750|
|50.||Ellison CA,Evans HC,Djeddour DH,Thomas SE. Year: 2008Biology and host range of the rust fungus Puccinia spegazzinii: A new classical biological control agent for the invasive, alien weed Mikania micrantha in Asia.Biol Control45133145|
|51.||Cock MJW,Ellison CA,Evans HC,Ooi PAC. Year: 2000Can failure be turned into success for biological control of mile-a-minute weed (Mikania micrantha)?Spencer NRProceedings of the X international symposium on biological control of weeds155167 Bozeman, Montana, USA.|
|52.||Huang ZL,Cao HL,Liang XD,Ye WH,Feng HL,et al. Year: 2000The growth and damaging effect of Mikania micrantha in different habitats.J Trop Subtrop Bot8131138|
|53.||Ismail BS,Chong TV. Year: 2002Effects of aqueous extracts and decomposition of Mikania micrantha H. B. K. debris on selected agronomic crops.Weed Biol Manag23138|
|54.||Ni GY,Song LY,Zhang JL,Peng SL. Year: 2006Effects of root extracts of Mikania micrantha H. B. K. on soil microbial community.Allelopathy J17247254|
|55.||Zan QJ,Wang YJ,Wang BS,Liao WB,Li MG. Year: 2000The distribution and harm of the exotic weed Mikania micrantha.Chinese J Ecol195861|
|56.||Zhang LY,Ye WH,Cao HL,Feng HL. Year: 2004Mikania micrantha H. B. K. in China – an overview.Weed Res444249|
|57.||Maron JL,Vilà M,Bommarco R,Elmendorf S,Beardsley P. Year: 2004Rapid evolution of an invasive plant.Ecol Monogr74261280|
|58.||Baker SA,Dyer RJ. Year: 2011Invasion genetics of Microstegium vimineum (Poaceae) within the James River Basin of Virginia, USA.Conserv Genet12793803|
|59.||Marshall DR,Brown ADH. Year: 1975Optimum sampling strategies for gene conservation.Frankel OH,Hawkes JGCrop genetic resources for today and tomorrowCambridgeCambridge University Press5380|
|60.||Wright S. Year: 1943Isolation by distance.Genetics2811413817247074|
|61.||Evanno G,Regnaut S,Goudet J. Year: 2005Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study.Mol Ecol142611262015969739|
|62.||Vasemägi A,Primmer CR. Year: 2005Challenges for identifying functionally important genetic variation: the promise of combining complementary research strategies.Mol Ecol143623364216202085|
|63.||Egan SP,Nosil P,Funk DJ. Year: 2008Selection and genomic differentiation during ecological speciation: isolating the contributions of host association via a comparative genome scan of Neochlamisus bebbianae leaf beetles.Evolution621162118118298645|
|64.||Williams MW,Oleksiak MF. Year: 2008Signatures of selection in natural populations adapted to chronic pollution.BMC Evol Biol828218847479|
|65.||Tice KA,Carlon DB. Year: 2011Can AFLP genome scans detect small islands of differentiation? The case of shell sculpture variation in the periwinkle Echinolittorina hawaiiensis.J Evol Biol241814182521605221|
|66.||Pimental D,Lach L,Zuniga R,Morrison D. Year: 2000Environmental and economic costs of nonindigenous species in the United States.BioScience505365|
|67.||Rieseberg LH,Kim SC,Randell RA,Whitney KD,Gross BL,et al. Year: 2007Hybridization and the colonization of novel habitats by annual sunflowers.Genetica12914916516955330|
|68.||Shimada Y,Shikano T,Merila J. Year: 2011A high incidence of selection on physiologically important genes in the three-spined stickleback, Gasterosteus aculeatus.Mol Biol Evol2818119320660084|
|69.||Kuchma O,Finkeldey R. Year: 2011Evidence for selection in response to radiation exposure: Pinus sylvestris in the Chernobyl exclusion zone.Environ Pollut1591606161221429640|
|70.||Nosil P,Egan SP,Funk DJ. Year: 2008Heterogeneous genomic differentiation between walking stick ecotypes: “Isolation by adaptation” and multiple roles for divergent selection.Evolution6231633617999721|
|71.||Smith TB,Milá B,Grether GF,Slabbekoorn H,Sepil I,et al. Year: 2008Evolutionary consequences of human disturbance in a rainforest bird species from Central Africa.Mol Ecol17587117868295|
|72.||Nunes VR,Beaumont MA,Butlin RK,Paulo OS. Year: 2011Multiple approaches to detect outliers in a genome scan for selection in ocellated lizards (Lacerta lepida) along an environmental gradient.Mol Ecol2019320521091562|
|73.||Nosil P,Funk DJ,Ortiz-Barrientos D. Year: 2009Divergent selection and heterogeneous genomic divergence.Mol Ecol1837540219143936|
|74.||Leroux S,Feve K,Vignoles F,Bouchez O,Klopp C,et al. Year: 2010Non PCR-amplified transcripts and AFLP® fragments as reduced representations of the quail genome for 454 Titanium sequencing.BMC Res Note3214|
|75.||Novak SJ,Mack RN. Year: 2005Genetic bottlenecks in alien plant species: influences of mating systems and introduction dynamics.Sax DF,Stachowicz JJ,Gaines SDSpecies invasions: insights into ecology, evolution, and biogeographySunderlandSinauer & Associates201228|
|76.||Hedrick PW. Year: 1999Perspective: highly variable loci and their interpretation in evolution and conservation.Evolution53313318|
|77.||Hedrick PW. Year: 2000Genetics of PopulationsSudburyJones and Bartlett Publishers|
|78.||Van Buskirk J,Willi Y. Year: 2006The change in quantitative genetic variation with inbreeding.Evolution602428243417263106|
|79.||Walker NF,Hulme PE,Hoelzel AR. Year: 2003Population genetics of an invasive species, Heracleum mantegazzianum: implications for the role of life history, demographics and independent introductions.Mol Ecol121747175612803628|
|80.||Chapman H,Robson B,Pearson ML. Year: 2004Population genetic structure of a colonising, triploid weed, Hieracium lepidulum.Heredity9218218814679390|
|81.||Durka W,Bossdorf O,Prati D,Auge H. Year: 2005Molecular evidence for multiple introductions of garlic mustard (Alliaria petiolata, Brassicaceae) to North America.Mol Ecol141697170615836643|
|82.||Zhou XM,Huang BQ. Year: 2001The invasion and control of Mikania micrantha.World Agr104243|
|83.||Swarmy PS,Ramakrishnan PS. Year: 1987Weed potential of Mikania micrantha H.B.K. and its control in fallows after shifting agriculture (Jhum) in north-east India.Agr Ecosyst Environ18195204|
|84.||Wen DZ,Ye WH,Feng HL,Cai CX. Year: 2000Comparison of basic photosynthetic characteristics between exotic invader weed Mikania micrantha and its companion species.J Trop Subtrop Bot8139146|
|85.||Hong L,Shen H,Ye WH,Cao HL,Wang ZM. Year: 2007Self-incompatibility in Mikania micrantha in South China.Weed Res47280283|
|86.||Price SC,Jain SK. Year: 1981Are inbreeders better colonisers?Oecologia49283286|
|87.||Husband BC,Barrett SCH. Year: 1991Colonization history and population genetic structure of Eichornia paniculata in Jamaica.Heredity66287296|
|88.||Amsellem L,Noyer JL,Hossaert-Mckey M. Year: 2001Evidence for a switch in the reproductive biology of Rubus alceifolius (Rosaceae) towards apomixis, between its native range and its area of introduction.Am J Bot882243225121669657|
|89.||Burdon JJ,Marshall DR. Year: 1981Biological control and the reproductive mode of weeds.J Appl Ecol18649658|
|90.||Barrett SCH,Richardson BJ. Year: 1986Genetic attributes of invading species.Groves RH,Burdon JJEcology of biological invasions: an Australian perspectiveCanberraAustralian Academy of Science2133|
|91.||Pellegrin D,Hauber DP. Year: 1999Isozyme variation among populations of the clonal species, Phragmites australis (Cav.) Trin. ex Steudel.Aquat Bot63241259|
|92.||Klekowski EJ Jr. Year: 1997Somatic mutation theory of clonality.de Kroon H,van Groenendael JThe ecology and evolution of clonal plantsLeidenBackhuys Publishers227241|
|93.||Houliston GJ,Chapman HM. Year: 2004Reproductive strategy and population variability in the facultative apomict Hieracium pilosella (Asteraceae).Am J Bot91374421653360|
|94.||Rottenberg A,Parker JS. Year: 2004Asexual populations of the invasive weed Oxalis pes-caprae are genetically variable.Proc Roy Soc B-Biol SciSuppl 271S206S208|
|95.||Lavergne S,Molofsky J. Year: 2007Increased genetic variation and evolutionary potential drive the success of an invasive grass.Proc Natl Acad Sci USA1043883388817360447|
|96.||Su YJ,Wang T,Zheng B,Jiang Y,Chen GP,et al. Year: 2005Genetic differentiation of relictual populations of Alsophila spinulosa in southern China inferred from cpDNA trnL–F noncoding sequences.Mol Phylogenet Evol3432333315619445|
|97.||Vos P,Hogers R,Bleeker M,Reijans M,van de Lee T,et al. Year: 1995AFLP: a new technique for DNA fingerprinting.Nuc Acids Res2344074414|
|98.||Nei M. Year: 1973Analysis of gene diversity in subdivided populations.Proc Natl Acad Sci USA70332133234519626|
|99.||Yeh FC,Yang R. Year: 1999POPGENE version 1.31. Dept of Renewable Resources, University of Alberta. Available: http://www.ualberta.ca/_/fyeh.|
|100.||Lewontin RC. Year: 1973The apportionment of human diversity.Evol Biol6381398|
|101.||Excoffier L,Laval G,Schneider S. Year: 2005Arlequin ver. 3.0: An integrated software package for population genetics data analysis.Evolutionary Bioinformatics Online1475019325852|
|102.||Holsinger KE,Lewis PO,Dey DK. Year: 2002A Bayesian approach to inferring population structure from dominant markers.Mol Ecol111157116412074723|
|103.||Pritchard JK,Stephens M,Donnelly P. Year: 2000Inference of population structure using multilocus genotype data.Genetics15594595910835412|
|104.||Falush D,Stephens M,Pritchard JK. Year: 2007Inference of population structure using multilocus genotype data: dominant markers and null alleles.Mol Ecol Notes757457818784791|
|105.||Nei M. Year: 1978Estimation of average heterozygosity and genetic distance from a small number of individuals.Genetics8958359017248844|
|106.||Miller MP. Year: 1997Tools for Populations Genetic Analyses (TFPGA) 1.3: A Windows program for the analysis of allozyme and molecular population genetic data. Department of Biological Sciences, Northern Arizona University.|
|107.||Luikart G,Cornuet JM. Year: 1999BOTTLENECK: a program for detecting recent effective population size reductions from allele data frequencies. Available: http://www.montpellier.inra.fr/URLB/bottleneck/bottleneck.html.|
|108.||Di Rienzo A,Peterson AC,Garza JC,Valdes AM,Slatkin M,et al. Year: 1994Mutational processes of simple sequence repeat loci in human populations.Proc Natl Acad Sci USA91316631708159720|
|109.||Godwin ID,Aitken EAB,Smith LW. Year: 1997Application of inter simple sequence repeat (ISSR) markers to plant genetics.Electrophoresis18152415289378115|
|110.||Cornuet JM,Luikart G. Year: 1996Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data.Genetics144200120148978083|
|111.||Kim YJ,Feng S,Zeng ZB. Year: 2008Measuring and partitioning the high-order linkage disequilibrium by multiple order Markov chains.Genet Epidemiol3230131218330903|
|112.||Brown AHD,Feldman MW,Nevo E. Year: 1980Multilocus structure of natural populations of Hordeum spontaneum.Genetics9652353617249067|
|113.||Maynard SJ,Smith NH,O'Rourke M,Spratt BG. Year: 1993How clonal are bacteria?Proc Natl Acad Sci USA90438443888506277|
|114.||Haubold B,Travisano M,Rainey PB,Hudson RR. Year: 1998Detecting linkage disequilibrium in bacterial populations.Genetics150134113489832514|
|115.||Agapow PM,Burt A. Year: 2001Indices of multilocus linkage disequilibrium.Mol Ecol Notes1101102|
|116.||Mes THM. Year: 1998Character compatibility of molecular markers to distinguish asexual and sexual reproduction.Mol Ecol717191727|
|117.||Van Der Hulst RG,Mes TH,Den Nijs JC,Bachmann K. Year: 2000Amplified fragment length polymorphism (AFLP) markers reveal that population structure of triploid dandelions (Taraxacum officinale) exhibits both clonality and recombination.Mol Ecol91810652071|
|118.||Wilkinson M. Year: 2001PICA 4.0: software and documentation distributed by the Department of Zoology, Natural History Museum, London, UK. Available: http://www.nhm.ac.uk/research-curation/projects/software/mwphylogeny.html.|
|119.||Hassel K,Ståstad SM,Gunnarsson U,Söderström L. Year: 2005Genetic variation and structure in the expanding moss Pogonatum dentatum (Polytrichaceae) in its area of origin and in a recently colonized area.Am J Bot921684169021646085|
|120.||Caballero A,Quesada H,Rolán-Alvarez E. Year: 2008Impact of amplified fragment length polymorphism size homoplasy on the estimation of population gene diversity and the detection of selective loci.Genetics17953955418493070|
|121.||Luikart G,England PR,Tallmon D,Jordan S,Taberlet P. Year: 2003The power and promise of population genomics: from genotyping to genome typing.Nat Rev Genet498199414631358|
|122.||Bonin A,Ehrich D,Manel S. Year: 2007Statistical analysis of amplified fragment length polymorphism data: a toolbox for molecular ecologists and evolutionists.Mol Ecol163737375817850542|
|123.||Herrera CM,Bazaga P. Year: 2008Population-genomic approach reveals adaptive floral divergence in discrete populations of a hawk moth pollinated violet.Mol Ecol175378539019121004|
|124.||Manel S,Conord C,Després L. Year: 2009Genome scan to assess the respective role of host-plant and environmental constraints on the adaptation of a widespread insect.BMC Evol Biol928820003282|
|125.||Miller N,Ciosi M,Sappington TW,Ratcliffe ST,Spencer JL,et al. Year: 2007Genome scan of Diabrotica virgifera virgifera for genetic variation associated with crop rotation tolerance.J Appl Entomol131378385|
|126.||Balding DJ,Nichols RA. Year: 1995A method for quantifying differentiation between populations at multi-allelic loci and its implications for investigating identity and paternity.Genetica963127607457|
|127.||Rannala B,Hartigan JA. Year: 1996Estimating gene flow in island populations.Genet Res671471588801187|
|128.||Balding DJ. Year: 2003Likelihood-based inference for genetic correlation coefficients.Theor Popul Biol6322123012689793|
|129.||Jeffreys H. Year: 1961Theory of probability (third edition)OxfordOxford University Press|
|130.||Geleta M,Bryngelsson T,Bekele E,Dagne K. Year: 2007AFLP and RAPD analyses of genetic diversity of wild and/or weedy Guizotia (Asteraceae) from Ethiopia.Hereditas144536217567442|
|131.||Bahulikar RA,Stanculescu D,Preston CA,Baldwin IT. Year: 2004ISSR and AFLP analysis of the temporal and spatial population structure of the post-fire annual, Nicotiana attenuata, in SW Utah.BMC Ecology41215350209|
|132.||Nielsen LR. Year: 2004Molecular differentiation within and among island populations of the endemic plant Scalesia affinis (Asteraceae) from the Galápagos Islands.Heredity9343444215280895|
|133.||Fjellheim S,Rognli OA. Year: 2005Molecular diversity of local Norwegian meadow fescue (Festuca pratensis Huds.) populations and Nordic cultivars-consequences for management and utilisation.Theor Appl Genet11164065016034583|
|134.||Haldimann P,Steinger T,Müller-Schärer H. Year: 2003Low genetic differentiation among seasonal cohorts in Senecio vulgaris as revealed by amplified fragment length polymorphism analysis.Mol Ecol122541255112969460|
|135.||Schöswetter P,Paun O,Tribsch A,Niklfeld H. Year: 2003Out of the Alps: colonization of Northern Europe by East Alpine populations of the Glacier Buttercup Ranunculus glacialis L. (Ranunculaceae).Mol Ecol123373338114629353|
|136.||Paun O,Greilhuber J,Temsch EM,Hörandl E. Year: 2006Patterns, sources and ecological implications of clonal diversity in apomictic Ranunculus carpaticola (Ranunculus auricomus complex, Ranunculaceae).Mol Ecol1589791016599955|
|137.||Jørgensen S,Mauricio R. Year: 2004Neutral genetic variation among wild North American populations of the weedy plant Arabidopsis thaliana is not geographically structured.Mol Ecol133403341315487999|
|138.||Treuren RV,Bas N,Goossens PJ,Jansen J,Van Soest LJM. Year: 2005Genetic diversity in perennial ryegrass and white clover among old Dutch grasslands as compared to cultivars and nature reserves.Mol Ecol14395215643949|
|139.||Menchari Y,Délye C,Le Corre V. Year: 2007Genetic variation and population structure in black-grass (Alopecurus myosuroides Huds.), a successful, herbicide-resistant, annual grass weed of winter cereal fields.Mol Ecol163161317217651194|
[Figure ID: pone-0041310-g001]
Figure 1 Map of sampled populations.
(A) Mikania micrantha. (B) Mikania cordata.
[Figure ID: pone-0041310-g002]
Figure 2 Bayesian assignment proportions for K = 28 clusters determined in STRUCTURE 2.3.3.
Each vertical bar represents one individual.
[Figure ID: pone-0041310-g003]
Figure 3 UPGMA dendrogram derived from AFLP data by Nei's  unbiased genetic distances.
It shows the relationships among 28 examined populations of Mikania micrantha. Populations of Mikania cordata are used as the root. Numbers above branches indicate bootstrap values (% of 1000 replicates). Only values larger than 40% are displayed. Branch lengths are proportional to genetic distances (see scale at the bottom of figure).
[Figure ID: pone-0041310-g004]
Figure 4 Results of the simulations with Dfdist for outlier detection.
Plots representing FST values are against heterozygosity. Each dot indicates an AFLP locus. The lower, intermediate, and higher lines represent the 0.5%, 50%, and 99.5% confidence intervals, respectively. Loci above the 99.5% line are regarded as outlier loci. (A) The result of Mikania micrantha. The 23 outlier loci under selection are represented by red dots accompanied by the locus number. (B) The result of Mikania cordata. No outlier locus under selection is detected.
[Figure ID: pone-0041310-g005]
Figure 5 Genomic scan to identify outlier loci subject to selection by Bayescan approach.
Each point corresponds to an AFLP locus. FST is plotted against the log10 of the posterior odds (PO), which provides evidence whether the locus is subject to selection or not. The vertical dashed line shows the decisive threshold value (log10 PO = 2.0) used for identifying outlier loci. (A) The result of Mikania micrantha. Shown as red dots, the 38 outlier loci with the corresponding locus number are candidates for being under positive selection. (B) The result of Mikania cordata. No outlier locus under selection is detected.
Table 1 Locations of Mikania micrantha and Mikania cordata populations surveyed in this study.
|Species||Region||Population||Location||Latitude||Longitude||Sample size||Altitude (m)|
|M. micrantha||Hong Kong||HK81||Hong Kong Island, Stubbs Road, in shrubs||22° 16′ 08″ N||114° 10′ 49″ E||15||74|
|HK82||Hong Kong Island, Barker Road, in tussock||22° 16′ 09″ N||114° 09′ 50″ E||13||300|
|HK84||Hong Kong Island, Mount Gough, in tussock||22° 16′ 05″ N||114° 09′ 43″ E||16||315|
|HK85||Hong Kong Island, Victoria Peak, in tussock||22° 16′ 30″ N||114° 08′ 57″ E||15||503|
|HK87||Hong Kong Zoological and Botanical Gardens, under bamboo forest||22° 16′55″N||114° 09′23″ E||12||372|
|HK88||New Territories, Hok Tau, in tussock||22° 29′ 53″ N||114° 10′ 49″ E||15||46|
|HK89||New Territories, Luk Keng, Pat Sin Leng Country Park, rivulet-side||22° 31′ 26″ N||114° 12′ 57″ E||12||8|
|HK90||Kowloon, Hong Kong Baptist University, slope||22° 20′ 07″ N||114° 10′ 57″ E||19||25|
|Macao||MA101||Hác Sá Beach, wasteland||22° 07′ 08″ N||113° 34′ 04″ E||15||3|
|MA105||Hác Sá Village, roadside||22° 07′ 02″ N||113° 34′ 03″ E||15||1|
|MA106||Hác Sá Beach, building site||22° 07′16″N||113°34′10″E||14||8|
|Shenzhen||SZ34||Fairy Lake Botanical Garden, Liang Yi Ting, in tussock||22° 34′ 49″ N||114° 10′ 08″ E||14||88|
|SZ35||Fairy Lake Botanical Garden, Desert Plant Section, roadside||22°35′10″N||114°10′26″E||16||42|
|SZ64||The office of Mangrove Natural Reserve, roadside||22° 31′ 58″ N||114° 00′ 01″ E||11||19|
|SZ72||Meilin Park, Huanshan Road, slope||22° 33′ 59″ N||114° 01′ 27″ E||16||83|
|SZ75||Lotus Hill Park, the Kite Square, in tussock||22° 33′ 14″ N||114° 03′ 22″ E||14||38|
|SZ77||Lotus Hill Park, under the Eucalypt forest||22° 33′ 25″ N||114° 02′ 57″ E||16||46|
|Neilingding||NLD1||Management station, roadside||22° 23′ 48″ N||113° 49′ 09″ E||15||3|
|NLD10||Dong Jiao Zui Wan, slope, in shrubs||22° 24′ 06″ N||113° 48′ 44″ E||14||51|
|NLD15||Dong Jiao Zui Wan, slope, under Spiny date palms||22°24′21″N||113°48′42″E||17||145|
|NLD26||Nan Wan, in shrubs||22° 23′ 41″ N||113° 48′ 52″ E||14||3|
|NLD29||Dong Wan, ravine||22° 23′ 48″ N||113° 49′ 38″ E||15||15|
|NLD30||Bei Wan Ma Guan, in shrubs||22° 25′ 12″ N||113° 47′ 17″ E||15||3|
|NLD31||Management station East, ravine||22° 24′ 01″N||113° 49′ 59″ E||14||8|
|Zhuhai||ZH43||Qi'ao Island, roadside||22° 24′ 37″ N||113° 38′ 38″ E||15||2|
|ZH50||Qi'ao Island, No Jia Le, roadside||22° 24′ 03″ N||113° 37′ 40″ E||12||2|
|Dongguan||DG91||Da Ling Shan forestry centre, Shan Zhu Wo||22° 51′ 51″ N||113° 46′ 21″ E||16||174|
|DG92||Da Ling Shan forestry centre, Chang Keng Kou||22° 51′ 38″ N||113° 46′ 27″ E||15||114|
|M. cordata||Hainan Island||HN1||Xinglong, in arbors||18° 42′04″ N||110° 13′ 23″ E||12||49|
|HN2||Xinglong, in tussock||18° 42′03″ N||110° 13′ 22″ E||5||44|
Table 2 Estimates of genetic diversity, and test for linkage disequilibrium and character compatibility in populations of Mikania micrantha.
|Population||Number of loci||Number of polymorphic loci||Percentage of polymorphic loci||Nei's gene diversity||IA||IER|
IA, index of association; , modified index of association; IER, incompatibility excess ratio.
Table 3 Genetic bottleneck of Mikania micrantha populations from six introduced regions in southern China.
P values are determined by a sign test under the stepwise mutation model (SMM) and the infinite allele model (IAM). He/Hd, the heterozygosity excess/deficiency ratio.
Table 4 Analysis of molecular variance (AMOVA) of 483 AFLP loci for 28 Mikania micrantha populations from six introduced regions in southern China.
|Variance components||Percentage of total variation||P||Φ statistics|
|Among regions||2.884||4.74||<0.001||ΦCT = 0.0474|
|Among populations within regions||17.404||28.61||<0.001||ΦSC = 0.3004|
|Within populations||40.538||66.65||<0.001||ΦST = 0.3335|
The P-value was calculated by a permutation procedure based on 1023 replicates.
Table 5 Comparison of outlier loci of Mikania micrantha under selection using Bayescan, Dfdist, and both with Dfdist and Bayescan, respectively.
|Approach||Outiler loci identified|
|Bayescan||3, 4, 5, 11, 14, 17, 25, 43, 54, 57, 59, 116, 132, 153, 166, 167, 173, 177, 183, 184, 209, 219, 224, 228, 275, 276, 325, 347, 381, 388, 420, 425, 426, 434, 451, 455, 473, 483|
|Dfdist||3, 17, 25, 35, 50, 56, 57, 116, 132, 136, 177, 183, 184, 219, 231, 240, 324, 325, 342, 347, 381, 426, 465|
|Bayescan and Dfdist||3, 17, 25, 57, 116, 132, 177, 183, 184, 219, 325, 347, 381, 426|
Table 6 Percentage of polymorphic loci (PL), Nei's total gene diversity (HT), GST, ΦST, and θB obtained from the populations of different weedy species based on AFLP data.
|Guizotia scabra ssp. scabrai||0.8455||0.32||0.18||0.22||Geleta et al. |
|Guizotia scabra ssp. schimperi||0.9002||0.32||0.19||0.17||Geleta et al. |
|Guizotia villosa||0.8393||0.33||0.19||0.26||Geleta et al. |
|Nicotiana attenuata||0.961||0.114||0.0549||Bahulikar et al. |
|Scalesia affinis||0.549||0.435||Nielsen |
|Oxalis pes-caprae*||0.8843||0.23||Rottenberg & Parker |
|Festuca pratensis||0.9320||0.31||Fjellheim & Rognli |
|Senecio vulgaris||0.34||0.1808||Haldimann et al. |
|Ranunculus glacialis||0.99||Schönswetter et al. |
|Ranunculus carpaticola||0.8594||0.7442||Paun et al. |
|Arabidopis thaliana||1.00||0.179||0.279||0.23||Jørgensen & Mauricio |
|Lolium perenne||0.886||Treuren et al. |
|Alopecurus myosuroides||0.905||Menchari et al. |
|Mikania micrantha||1.00||0.3376||0.4736||0.3335||0.2927||This Study|
Previous Document: A novel source of methylglyoxal and glyoxal in retina: implications for age-related macular degenera...
Next Document: Effect of intravitreal anti-vascular endothelial growth factor therapy on the risk of arterial throm...