- Open Access
In the wild hybridization of annual Datura species as unveiled by morphological and molecular comparisons
Journal of Biological Research-Thessaloniki volume 21, Article number: 11 (2014)
The Erratum to this article has been published in Journal of Biological Research-Thessaloniki 2014 21:18
The present work aimed to verify whether intermediate variants were natural crosses between Datura species (D. stramonium forms and D. ferox). Their existence has been long ago insinuated but has not been studied using morphological features and molecular tools. The variants differed in stem coloring, upper bearing forks, and fruit characters.
Principal Components Analysis of 11 morphological characteristics showed that D. ferox and D. stramonium (forms stramonium and tatula) were quite different and the putative hybrids were intermittent. The D. ferox × D. stramonium f. tatula was closer to the latter of its parents. Sequencing analysis revealed identical amplified trnL intron in all variants and a 100% homology with D. stramonium accession number EU580984.1 suggested that this plastid cannot discern Datura variants. However, genomic analysis with URP markers indicated that the hybrids had >60% genetic makeup similarity with both parents suggesting that the intermediate variants were putative inter-specific hybrids. Moreover, the dendrogram stemmed from cluster analysis of the fingerprint profile of variants placed D. stramonium and D. ferox in different branches indicating their genetic differentiation from each other as well as from their hybrids.
The findings suggest that the natural hybridization of annual Datura species occurs. Extrapolating, this hybridization could be the first step for speciation. More possibly, it can alter population composition, its weediness and adaptability to local conditions.
Genus Datura, family Solanaceae, consists of nine (annual and tree) species, originating from the New and Old World . In Greece, Datura species are considered to be invasive , known since antiquity for their narcotic and medicinal actions [3, 4]. Nowadays in Greece, D. stramonium L. forms (f. stramonium and f. tatula), D. ferox L., and D. innoxia L. coexist in mixed populations in various combinations and relative ratios. Datura stramonium f. stramonium is the most common variant found as spring weed in fields, roadsides and dumps and usually coexists with the recently identified to occur in Greece D. stramonium f. tatula L. . Fierce thornapple (D. ferox L.) is the dominant Datura species in some sites in northern Greece where it shares habitat with D. stramonium. Finally, D. innoxia, commonly used as ornamental, is found as feral, usually, in dumps .
Datura stramonium is a predominantly self-fertilized species but cross-pollination is feasible to some extent by insects like hawkmoths and honeybees . The predominance of self-fertilization is ascribed to anther-stigma overlapping and results in inbreeding [7, 8], which found to reduce vigor and increase herbivory damages . However, within Datura populations, there are plants showing herkogamy (anther-stigma separation), which permits outcrossing at low rates ranging from 1.3% to 8.5% [6, 10]. The existence of anthocyanin in the purple-stemmed, violet-flowered D. stramonium f. tatula does not affect outcrossing rates .
In tree Datura species, inter-specific crossing is feasible (D. aurea × D. candida) descending hybrids with increased alkaloid content, which could be of economic interest . Regarding annual species, Husaini & Iwo  reported the existence of cytological compatibility between D. stramonium and D. ferox. Weaver & Warwick , reporting the findings of Rietsema , stated that these two Datura species are the only annual species, which gave identified hybrids in the wild, collected in South America. Since then, no relative report exists, making morphological and molecular confirmation of analogous findings necessary. Given that crosses between the annual Datura species yielded viable seeds , possible outcrossing in mixed communities would diversify the existing populations to new ones, with unpredicted features regarding competitiveness, resistance to herbivores, alkaloid content and herbicide tolerance.
During September 2011, in mixed Datura swards in northern Greece, specimens were found showing morphological features intermediate to those of the coexisting D. stramonium forms and D. ferox. Thus, the aim of this study was to verify naturally occurring Datura crosses and unmask the possible hybridization by morphological features and molecular tools.
The five variants differed significantly for all the characteristics determined (Table 1). Regarding stem color (SC), the putative D. ferox × f. stramonium hybrid had green stem resembling f. stramonium, while the putative D. ferox × f. tatula hybrid followed the stem coloring of f. tatula (Table 1).
The use of PCA with varimax rotation of the 11 morphological characteristics revealed two significant components which extracted the 70% of the total variance (Table 2). The component PC1 with eigenvalue 4.4 extracted the 40% of the total variance and the second component, PC2, with eigenvalue 3.3, the 30%. According to the results of parallel analysis, the critical value for significant eigenvalues was 2.1. Table 2 presents the loadings of the 11 morphological characteristics on the two components.
The PC1 was mainly correlated (positively or negatively) with long spine length, medium spine length, spines/capsule, short spine length, capsule length, and leaf length (Table 2). The PC2 was mainly loaded (positively or negatively) by the purple color, leaf width, length/width, and capsule width. The green color had almost the same negative loading on both components.
The projection of the 30 plant individuals on the 1 × 2 factorial plane derived from PCA indicated that the first component PC1 clearly separated the D. ferox from f. stramonium individuals while D. ferox × f. stramonium plants were projected around the origin almost in the middle of their two parents (Figure 1). The second component PC2 distinguished D. ferox and f. stramonium individuals from the f. tatula and D. ferox × f. tatula plants. The D. ferox × f. tatula variants were most similar with one of their parents, f. tatula. The individuals of D. ferox, f. stramonium, and D. ferox × f. stramonium showed greater variability along the PC1 axis than the f. tatula and D. ferox × f. tatula plants, which had very little variability on the same axis; in contrast they showed some variability along the PC2 axis.
The PCR fingerprinting, using all the primers and DNA, detected 63 clear and scorable bands in the samples including D. stramonium, D. ferox, and the two putative inter-specific hybrids (Figure 2). The band patterns of the D. stramonium forms were identical and only results from one form are shown. The URP primers produced multiple bands in all variants, varying in size from about 100 bp to more than 3000 bp. Polymorphic as well as monomorphic bands were revealed with 11 from the 12 URP primers. Only primer URP13R (number 7 in Figure 2) produced a single monomorphic band in all samples. Out of 63 scorable bands, 37 bands (58.73%) were found to be polymorphic (present in one to three variants) while 26 bands (41.27%) were monomorphic (present in all variants). Totally, 20 polymorphic bands (31.75%) were present in only one variant, while six were present in two (9.5%) and 11 in three variants (17.5%).
The 12 different primers generated various banding patterns, ranging from 1 to 9. The maximum number of scorable bands (9) was observed in primers 1F and 6R. Primer 1F also produced the highest number (8) of polymorphic bands (88.88%) and thus, it showed the higher level of polymorphism. On the other hand, the primer 13R produced a single monomorphic band. The numbers of total, polymorphic, and monomorphic bands indicated an average of 5.25 bands per primer of which an average of 3.08 were polymorphic and 2.16 monomorphic (Table 3). The URP band patterns were used to determine the genetic distances between Datura variants. The genetic similarity for pairs of variants was calculated using the Jaccard’s coefficient. The similarity matrix based on all possible pairs had a similarity range from 46% to 86% (Table 4). The lower similarity value of 46% indicating the higher distance was between D. stramonium and D. ferox, while the higher similarity value of 86% was between the two putative inter-specific hybrids indicating the closer relationship (smaller genetic distance).
The dendrogram obtained from the cluster analysis of the examined fingerprint profile of the variants showed grouping of the putative hybrids in one cluster (Figure 3). The position of D. stramonium and D. ferox in different branches of the tree indicates their genetic differentiation from each other as well as from the putative hybrids.
The sequencing analysis revealed identical amplified trnL intron in all variants. BLAST similarity search of the GenBank for homologous sequences revealed that a sequence of D. stramonium with accession number EU580984.1 had 100% homology over the 505 nucleotide amplified region (data not shown).
In the wild, inter-specific hybridization is a common phenomenon especially among plant species and may become the first step in the process of speciation [16, 17]. Sympatric coexistence, intermediate characters, inter-fertility, and biochemical additivity are criteria determining the feasibility of inter-specific hybridization , and several of them were met in the present study.
The first clue for this work was the intermediate characters of specimens found to coexist with their putative parents (D. ferox and D. stramonium). Despite the high percentage of self-fertilization, annual Datura species are cytological compatible , show herkogamy [6, 10] and attract insect pollinators . In particular, under Greek conditions, f. tatula is highly attractive for honeybees .
The possible existence of hybrids of annual Datura species in the wild goes back to 1950’s . Nevertheless, exempting its insinuation, no proof has been provided. Nowadays, a combination of morphological and molecular characters is used to provide evidence for the existence of a putative hybrid [19–23]. Selecting an adequate number, in statistical terms (six specimens per variant), of a coherent plant cohort, we found out that the Datura variants were mainly different in stem coloring, upper bearing forks, and fruit characters that were determined. PCA revealed that D. ferox and f. stramonium were different and their putative hybrid was intermediate. The other putative hybrid, D. ferox × f. tatula, was closer to the latter of its putative parents. However, morphological intermediates can also be derived from convergent evolution or environmental selection and this makes ambiguous the confirmation of putative hybrid individuals based solely on morphological evidence. Use of molecular methods provides a number of advantages over morphological analysis, among others being the large number of available markers, their apparent selective neutrality, and the low levels of non-heritable variation . Ideally, additive molecular markers displaying co-dominant inheritance (e.g. SSR), combined with morphological analysis, should be applied to decipher putative hybridization. The lack of such SSR markers for Datura species makes their development cumbersome and time-consuming. Nevertheless, genetic analyses by dominant markers, like RAPDs and AFLPs, solitary or flanked by morphological analysis, have been routinely employed with ultimate success in inter-specific hybrid verification [25–27].
In the present work, genetic distances of the examined variants were estimated using URP primers, which have been proved valuable in genomic fingerprinting of a variety of organisms including plants, animals, and microorganisms . The URP markers were selected for their advantages over RAPDs (high annealing temperature leading to high specificity and reproducibility) and AFLPs (less laborious and equally specific) in genetic analysis. The applied URP primers were employed for first time in Datura species genetic fingerprinting, and proved to be successful in amplification of polymorphic fragments that are sufficient to discriminate their phylogenetic origin.
The fact that the sequencing analysis revealed identical amplified trnL intron in all variants and a 100% homology with D. stramonium accession number EU580984.1  indicates that the typically used plant plastid phylogenetic marker trnL intron is not useful for discrimination, since it is identical between D. stramonium and D. ferox and, consequently, in their putative inter-specific hybrids.
The molecular characterization of the different Datura accessions provided the most compelling evidence on the genetic relationships and genetic makeup of the putative inter-specific hybrids. Using Jaccard’s coefficient for comparison, the similarities determined between the bands of the examined variants bolstered the possibility of putative hybrids (D. ferox × f. stramonium and D. ferox × f. tatula) to be indeed hybrids between annual Datura species. Datura ferox was equally related to both of the putative hybrids with 62% similarity, while D. stramonium was somewhat more similar to D. ferox × f. stramonium (67%) than to D. ferox × f. tatula (60%). The two species were similar at a level of 46% while D. ferox × f. stramonium and D. ferox × f. tatula share 86% similarity. In accordance, Dymshakova et al. , using genetic distances estimated by AFLPs, showed that F1 hybrids were intermediate between the parentals Saxifraga sibirica and S. cernua. It is clear from the results that the putative hybrids were highly similar to each other, meaning that share common ancestry, and a high percentage of their genetic makeup is equally similar to both D. stramonium and D. ferox at a level higher than 60%, and hence those were possibly their progenitors. This evidence supports the hypothesis that the two intermediate accessions may be inter-specific hybrids between D. ferox and D. stramonium.
The natural hybridization of annual Datura species confers putative implications in biological and agronomic terms. In extreme, it could be the first step for speciation but more possibly, it could change local population composition, which in turn could affect Datura weed competitiveness and its susceptibility to chemical or mechanical control. Finally, this hybridization could raise alkaloid content leading to commercial interest of its extraction.
Site description and morphological measurements
During September 2011, mixed swards of Datura species were identified at the locales of Eukarpia (40°32′N, 22°60′E, 14 m a.s.l.) and Agios Demetrios (40°53′N, 23°41′E, 66 m a.s.l.) in Serres region, northern Greece. The predominant D. ferox was coexisting with the typical green-stemmed, white-flowered D. stramonium f. stramonium (f. stramonium) and the purple-stemmed, violet-flowered D. stramonium f. tatula (f. tatula). The swards emerged postharvest, after erratic summer rains, in winter cereal fields. The preceding crop was cotton, the spring crop where D. ferox mainly resides. The soil of the fields in the wider area was derived from alluvial deposits and the climate is Mediterranean, with dry and hot summers.
Determinations were conducted in situ at Eukarpia site, where Datura population cohort was more homogenous. Six specimens, typical of each variant (Figure 4), were selected for measurements, which were conducted on the upper full expanded leaves (length, width, and length/width ratio), flowers (corolla and calyx length) and capsules (length, width, number of spines, short spine length, medium spine length, and long spine length). Lengths were measured using a HOREX electronic digital caliper (Helios-Preisser, Gammertingen, Germany). Stem color of the specimens was also recorded.
Plant material sampling and molecular determinations
The plant samples included Datura ferox (F), D. stramonium forms (S), and their putative inter-specific hybrids (D. ferox × f. tatula: A and D. ferox × f. stramonium: B).
Upper forks of each specimen were selected, sealed in a plastic bag, put in a portable refrigerator and transferred to the laboratory for DNA extraction, which was performed from 100 mg ground leaf tissue using the cetyl trimethylammonium bromide (CTAB) method according to the protocol outlined in the NucleoSpin®Plant II kit (Macherey Nagel GmbH & Co. KG, Düren, Germany).
Oligonucleotide primers and PCR conditions for fingerprinting and analysis of the cpDNA trnL intron
A KAPA Taq PCR kit was used to perform the PCR reactions (Kapa Biosystems Inc., Wilmington, USA). The PCR mixture in a total volume of 25 μl contained 1× PCR buffer A (containing MgCl2), separate MgCl2 solution, so that total Mg concentration was fixed at 3.0 mM, 0.2 mM of each dNTP (New England Biolabs, Ipswich, USA), 2.4 μM of each primer (forward and reverse) and 1.2 U KAPA Taq polymerase. The amount of genomic DNA added in the PCR mixture was about 250 ng. PCR amplification was carried out in a Veriti 96Well Thermal Cycler (Applied Biosystems, Foster City, USA) using the following profile: a first step of 3 min at 94°C; a second step of 10 cycles of 30 sec at 94°C, 45 sec at 56°C (touchdown: -1°C per cycle), 3 min at 72°C; a third step of 30 cycles of 30 sec at 94°C, 45 sec at 47°C, 3 min at 72°C; and a step of a final extension for 11 min at 72°C. For fingerprinting the different genotypes, 12 universal rice primers (URP) were used . URPs are 20-mer primers that were designed for fingerprinting rice genomes, based on the sequence of a rice-specific CACTA-like transposon element. The sequence of this repeated element (named pKRD) was not found in other plant, animal or fungal genomes . Despite that, sequences homologous to URP over a continuous range of 15 bases were widely observed on diverse genomes including Arabidopsis, human and bacteria . Thus, URP primers were used for fingerprinting other organisms , and were proved a useful tool for genetic characterization and grouping of most eukaryotic or prokaryotic genomes, especially at inter- and intraspecific levels [31, 32]. The URP are listed in Table 5. DNA fragments were detected by staining with ethidium bromide on a 3% MetaPhor™ agarose (Cambrex Bio Science, Copenhagen, Denmark) gel in TBE buffer.
For amplification of the cpDNA trnL intron, the PCR primes P1 and P2 were used (Table 5). PCR setup and instruments were as above. Amplification profile was: one cycle of 2 min at 94°C; 10 cycles of 30 sec at 94°C, 30 sec at 58°C (touchdown: -0.5°C per cycle), 50 sec at 72°C; 30 cycles of 30 sec at 94°C, 30 sec at 53°C, 50 sec at 72°C; one cycle of a final extension for 10 min at 72°C. DNA fragments were detected by staining with ethidium bromide on a 3% MetaPhor™ agarose (Cambrex Bio Science, Copenhagen, Denmark) gel in TBE buffer. The bands of the trnL intron were excised from the gel, DNA was recovered using a modified freeze-squeeze method  and sent for direct sequencing with the primers used in PCR.
Statistical analysis of morphological characteristics
With the exception of corolla and calyx lengths in D. ferox × f. tatula for which measurements were not taken due to the lack of appropriate samples, morphological trait means and standard errors were computed for summarizing the distributions of the corresponding variables. On the basis of morphological data (except stem colour), the five variants were compared with the Analysis of Variance (ANOVA) method. The Duncan’s multiple range test was used for means’ comparisons.
Principal Components Analysis (PCA) with varimax rotation was applied on the correlation matrix between the morphological variables (except corolla and calyx lengths) in order to study the groupings, similarities, and differences between the individuals of the five variants. Significant components were determined by the parallel analysis (PA) method . Since the stem colour was a nominal categorical variable, with three categories (purple, green, and grey-green), two dummy variables with binary coding (0, 1) were entered in the PCA; one variable for the purple colour and one for the green. A third dummy variable for the grey-green colour is redundant since it would be dependent on and negatively correlated with the other two dummy variables. Generally, the number of dummy-coded variables needed is one less than the number of modalities of the corresponding categorical variable. Since PCA is not a modelling but only a descriptive variance summarizing method, the use of binary along with scale variables is legitimate . The significance level for all hypotheses testing procedures was preset at p < 0.05. IBM SPSS package v. 20 (IBM Corp., New York, USA) was used for the analyses. Parallel analysis was conducted with the RanEigen v. 2.0 software .
Analysis of molecular data
Gel photographs were scored using the gel image analysis software GelAnalyzer 2010a (http://www.gelanalyzer.com/). The bands were binary coded with 1 or 0 for their presence or absence in each genotype, respectively, and the coded data were subjected to statistical analysis. Estimates of similarity among all genotypes were calculated from the Jaccard’s similarity coefficient using the IBM SPSS package v. 20 (IBM Corp., New York, USA). Hierarchical cluster analysis based on the Jaccard similarity matrix with the unweighted pair group method based on arithmetic averages (UPGMA)  was conducted using the software MEGA5 .
Amplified fragment length polymorphism
Analysis of variance
Cetyl trimethylammonium bromide
Long spike length
Medium spike length
Principal components analysis
Polymerase chain reaction
Random amplified polymorphic DNA
Short spike length
Simple sequence repeat
Unweighted pair group method based on arithmetic averages.
Geeta R, Gharaibeh W: Historical evidence for a pre-Columbian presence of Datura in the Old World and implications for a first millennium transfer from the New World. J Biosci 2007, 32: 1227–1244.
Arianoutsou M, Bazos I, Delipetrou P, Kokkoris Y: The alien flora of Greece: taxonomy, life traits and habitat preferences. Biol Invasions 2010, 12: 3525–3549.
Dafni A, Yaniv Z: Solanaceae as medicinal plants in Israel. J Ethnopharmacol 1994, 44: 11–18.
Higgins AJ: From ancient Greece to modern Athens: 3000 years of doping in competition horses. J Vet Pharmacol Therap 2006, 29: 4–8.
Tsialtas I, Eleftherohorinos I, Zarkos G, Christodoulou V, Tan K: New floristic records in the Balkans: 21. Report 130. Datura stramonium f. tatula (L.) Geerinck & Walravens. Phytol Balcan 2013, 19: 149–151.
van Kleunen M, Fischer M, Johnson SD: Reproductive assurance through self-fertilization does not vary with population size in the alien invasive plant Datura stramonium . Oikos 2007, 116: 1400–1412.
Motten AF, Stone JL: Heritability of stigma position and the effect of stigma-anther separation on outcrossing in a predominantly self-fertilizing weed, Datura stramonium (Solanaceae). Am J Bot 2000, 87: 339–347.
Stone JL, Motten AF: Anther-stigma separation is associated with inbreeding depression in Datura stramonium , a predominantly self-fertilizing annual. Evolution 2002, 56: 2187–2195.
Bello-Bedoy R, Núñez-Farfán J: Cost of inbreeding in resistance to herbivores in Datura stramonium . Ann Bot 2010, 105: 747–753.
Motten AF, Antonovics J: Determinants of outcrossing rate in predominantly self-fertilizing weed, Datura stramonium (Solanaceae). Am J Bot 1992, 79: 419–427.
Stone JL: Does anthocyanin affect outcrossing rates in Datura stramonium (Solanaceae)? Am J Bot 2000, 87: 348–354.
El-Dabbas SW, Evans WC: Alkaloids of genus Datura , section Brugmansia X, alkaloid content of Datura hybrids. Planta Med 1982, 44: 184–185.
Husaini SWH, Iwo GA: Cytomorphological studies in some species of the family Solanaceae from Jos Plateau, Nigeria. Feddes Repert 1990, 101: 41–47.
Weaver SE, Warwick SI: The biology of Canadian weeds. 64. Datura stramonium L. Can J Plant Sci 1984, 64: 979–991.
Rietsema J: Barriers to crossability: prefertilization. In The Genus Datura. Volume 20. Edited by: Avery AG, Satina S, Rietsema J. New York: Ronald Press Co; 1959:235–244.
Mallet J: Hybrid speciation. Nature 2005, 446: 279–283.
Mallet J: Hybridization as an invasion of the genome. Trends Ecol Evol 2007, 20: 229–237.
Gottlieb LD: Leaves of confidence in the analysis of hybridization in plants. Ann Mo Bot Gard 1972, 59: 435–446.
Baumel A, Ainouche ML, Misset MT, Gourret J-P, Bayer RJ: Genetic evidence for hybridization between the native Spartina maritima and the introduced Spartina alterniflora (Poaceae) in South-West France: Spartina × neyrautii re-examined. Plant Syst Evol 2003, 237: 87–97.
Lexer C, Welch ME, Raymond O, Rieseberg LH: The origin of ecological divergence in Helianthus paradoxus (Asteraceae): selection on transgressive characters in a novel hybrid habitat. Evolution 2003, 57: 1989–2000.
Zhang J-L, Zhang C-Q, Gao L-M, Yang J-B, Li H-T: Natural hybridization origin of Rhododendron agastum (Ericaceae) in Yunnan, China: inferred from morphological and molecular evidence. J Plant Res 2007, 120: 457–463.
Hayakawa H, Hamachi H, Matsuyama K, Muramatsu Y, Minamiya Y, Ito K, Yokoyama J, Fukuda T: Interspecific hybridization between Arisaema sikokianum and A. serratum (Araceae) confirmed through nuclear and chloroplast DNA comparisons. Am J Plant Sci 2011, 2: 521–526.
Du G-H, Zhang Z-Q, Li Q-J: Morphological and molecular evidence for natural hybridization in sympatric population of Roscoea humeana and R. cautleoides (Zingiberaceae). J Plant Res 2012, 125: 595–603.
Rieseberg LH, Linder CR, Seiler GJ: Chromosomal and genetic barriers to introgression in Helianthus. Genetics 1995, 141: 1163–1171.
Pooler MR, Riedel LGH, Bentz SE, Townsend AM: Molecular markers used to verify interspecific hybridization between hemlock ( Tsuga ) species. J Am Soc Hort Sci 2002, 127: 623–627.
Lee SY, Yae BW, Kim KS: Segregation patterns of several morphological characters and RAPD markers in interspecific hybrids between Dianthus giganteus and D. carthusianorum . Sci Hortic 2005, 105: 53–64.
Dymshakova OS, Semerikov VL, Lascoux M: AFLP analysis to estimate the genetic contribution of parents to progeny from hybridization between Saxifraga sibirica L. and S. cernua L. Russ J Ecol 2012, 43: 347–351.
Kang H-W, Park D-S, Go S-J, Eun M-Y: Fingerprinting of diverse genomes using PCR with universal rice primers generated from repetitive sequence of Korean weedy rice. Mol Cells 2002, 13: 281–287.
Olmstead RG, Bohs L, Migid HA, Santiago-Valentin E, Garcia VF, Collier SM: A molecular phylogeny of the Solanaceae. Taxon 2008, 57: 1159–1181.
Kang H-W, Kang K-K: Genomic characterization of Oryza species-specific CACTA-like transposon element and its application for genomic fingerprinting of rice varieties. Mol Breeding 2008, 21: 283–292.
An HS, Park JY, Myeong JI, An CM: Genetic relationships of Pacific abalone (Haliotidae) species determined using universal rice primer-polymerase chain reaction fingerprinting. Genet Mol Res 2013, 12: 6309–6318.
Jana TK, Singh NK, Koundal KR, Sharma TR: Genetic differentiation of charcoal rot pathogen, Macrophomina phaseolina , into specific groups using URP-PCR. Can J Microbiol 2005, 51: 159–164.
Kaloumenos NS, Chatzilazaridou SL, Mylona PV, Polidoros AN, Eleftherohorinos IG: Target-site mutation associated with cross-resistance to ALS-inhibiting herbicides in late watergrass ( Echinochloa oryzicola Vasing.). Pest Manag Sci 2013, 69: 865–873.
Lautenschlager GJ: A comparison of alternatives to conducting Monte Carlo analyses for determining parallel analysis criteria. Multivar Behav Res 1989, 24: 365–395.
Jolliffe IT: Principal Components Analysis. 2nd edition. New York: Springer; 2002.
Enzmann D: RanEigen: a program to determine the parallel analysis criterion for the number of principal components. Appl Psych Meas 1997, 21: 232.
Sneath PHA, Sokal RR: Numerical Taxonomy. San Francisco: Freeman; 1973.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 2011, 28: 2731–2739.
The authors declare that they have no competing interests.
Assistant Professor ITT: he found along with Professor IGE the Datura variants. He conducted morphological measurements, authored part of the text and synthesized the parts of the different contributors. Mrs EP, Dr PVM and Assistant Professor AP: decided and conducted the molecular analyses and authored the pertinent parts of the manuscript. Dr NK: he conducted morphological measurements. Assistant Professor GM: he was responsible for statistics. Professor IGE: he had the general supervision of the work, authored part of the text and corrected its final version. All authors read and approved the final manuscript.
An erratum to this article is available at http://dx.doi.org/10.1186/2241-5793-21-18.