Servicios Personalizados
Revista
Articulo
Indicadores
- Citado por SciELO
- Accesos
Links relacionados
- Similares en SciELO
Compartir
Agrociencia
versión On-line ISSN 2521-9766versión impresa ISSN 1405-3195
Agrociencia vol.50 no.7 Texcoco oct./nov. 2016
Biotechnology
Genetic diversity of the chloroplast (TrnL-F) region among populations of Pholisma culiacanum Y.
1 Universidad Autónoma de Sinaloa. Facultad de Agronomía. Culiacán de Rosales, km.17.5, Carretera Culiacán-Eldorado, 80000. Sinaloa. México. (orlandomer@hotmail.com), (spenap@hotmail.com).
2 Universidad Nacional Autónoma de México. Instituto de ecología. ERNO. Avenida Luis Donaldo Colosio M. esquina Madrid S/N, Campus UNISON, Apartado Postal 1039, 83000. Hermosillo, Sonora. México. (freaner@unam.mx).
Pholisma culiacanum is a root holoparasitic plant, endemic to the Mexican states of Sinaloa and Sonora that has cultural, nutritious, biological, and evolutionary importance. This plant is an adequate agricultural model as an alternative source of nutrients and as a tool for understanding photosynthesis’ molecular changes, from autotrophism to parasitism. The aim of this study was to analyze the genetic diversity of seven populations of P. culiacanum in Sinaloa and Sonora. Sequences from a chloroplast DNA (TrnL-F) fragment were analyzed in 70 samples, in order to calculate the genetic variation within (86.51 %) and between (13.49 %) populations. Correlation and simple linear regression analysis were performed for the genetic and geographic distances data. Haplotype and nucleotide diversity was detected at 0.85 and 0.27; 11 haplotypes were identified as well. The relationship between genetic and geographic distances was not significant (p>0.05); therefore, no evidence of isolation due to distance was found. Based on the analyzed information, areas of in situ exploitation and genetic resources conservation are proposed.
Key words: Holoparasitic; Pholisma culiacanum; genetic diversity; chloroplast; TrnL-F
Pholisma culiacanum es una holoparásita de raíz, endémica de los estados de Sinaloa y Sonora, México, y tiene importancia cultural, alimenticia, biológica y evolutiva. Esta planta es un modelo agrícola adecuado como fuente alternativa de nutrientes y para entender los cambios moleculares de la fotosíntesis del autotrofismo al parasitismo. El objetivo de este estudio fue analizar la diversidad genética de siete poblaciones de P. culiacanum de Sinaloa y Sonora. En 70 muestras se analizaron secuencias de un fragmento de ADN de cloroplasto (Trnt-TrnF) para calcular la variación genética dentro (86.51 %) y entre (13.49 %) poblaciones. Para los datos de distancias genéticas y geográficas se realizaron análisis de correlación y regresión lineal simple. La diversidad haplotípica y nucleotídica detectada fue 0.85 y 0.27; además, se identificaron 11 haplotipos. La red y filogenia de haplotipos mostraron similitud con tres grupos. La relación entre las distancias genéticas y geográficas no fue significativa, (p>0.05); por lo tanto, no hay evidencia de aislamiento por distancia. Con la información analizada se proponen zonas de explotación in situ y de conservación de recursos genéticos.
Palabras clave: Holoparásita; Pholisma culiacanum; diversidad genética; cloroplasto; TrnL-F
Introduction
Pholisma culiacanum is a root holoparasitic plant of great agricultural, nutritious, economic, cultural, biological, and evolutionary importance. Besides it represents an appropriate model to study changes in the photosynthesis process, from autotrophism to parasitism (Bungard, 2004).
Parasitic plants represent 1 % of angiosperms. There are four thousand different species worldwide, gathered in 265 genera, 36 of which are found in Mexico. There are six genera of parasitic plants in Sinaloa -including Pholisma, which can be found from southern Sinaloa, all the way to California and Arizona in the U.S. Out of the three species that belong to the genus, P. culiacanum is the one endemic to the states of Sonora and Sinaloa, and the area surrounding Culiacan -Sinaloa’s capital city- is likely to be the place of origin of the species and genus (Yatskievych et al., 1986).
The aim of this research was to determine the genetic variability in populations of P. culiacanum with genetic sequences of the TrnL-F region of the chloroplast DNA and to identify areas with good potential for both in situ exploitation and genetic resources conservation.
Materials and methods
Sampling
In order to locate P. culiacanum populations, different data bases were consulted: GIFB (Global Information of Biodiversity, 2013), Missouri Botanical Garden, University of Arizona, the Biology School and the Faculty of Agriculture of the Autonomous University of Sinaloa, and the Sonora Regional Station of the National Autonomous University of Mexico (ERNO-UNAM). Areas that comply with the distribution characteristics were also explored. In each population, 10 to 30 individuals were collected (with at least 10 m between them), they were georeferenced, and stored at -80 °C, until DNA purification was performed.
Amplification
For the amplification, the extraction methods proposed by Doyle (1991), Porebski et al. (1997), Sá et al. (2011), Cota-Sánchez et al. (2006), Krizman (2006), Sánchez-Hernández (2006), and Amani et al. (2011) were evaluated. The appropriate methodology was based on that described by Azmat et al. (2012).
Due to the lack of previous information on P. culiacanum’s genetic variability, several chloroplast DNA regions reported in other parasitic plants were evaluated, where sufficient variability to perform population genetics studies were found (Taberlet et al., 1991, de Pamphilis et al., 1997; Zuber et al., 2000; Heinze, 2007; Amico et al., 2009; Nickrent et al., 2009). Based on those results, c-Tab F and f-Tab R from Taberlet et al. (1991) were selected as strand primers, which amplify the intergenic region of Trn L-F. A touchdown-like PCR amplification protocol was used (Amico et al., 2007), with a modification in alignment temperatures: from 48 to 50 °C, and from 52 to 54 °C. Amplification reactions were performed in a Techne Unit Genius® thermocycler, visualized in a 1 % agarose gel, and stained with ethidium bromide (0.003 %) in an UV transilluminator.
The amplified product was sequenced in the HighThroughput Sequencing Solutions of the Genomic Science Department of the University of Washington, Seattle, U.S. To edit and align the sequences, two software were used: CodonCode Aligner (CodonCode Corporation, 2009), Finch TV 1.4.0 (Geospiza Inc., 2004-2006), and ClustalW multiple alignment from the BioEdit software suite (Hall, 1999).
Genetic variability analysis
DnaSP (Librado and Rosas, 2010) and ARLEQUIN 3.1 (Excoffier, 2005) software were used to calculate genetic diversity indexes: 1) nucleotide diversity (π) (Nei, 1987, Hedrick, 2005, Hartl and Clark, 2007); 2) haplotype diversity (Hd), with and without indels as informative sites (Nei, 1987); 3) variation components due to differences between and within populations (AMOVA). In order to evaluate if the populations fit in the isolation by distance model, correlation and linear regression analysis of the genetic and geographic distances were performed, using the Maximum Composite Likelihood model, without taking indels into consideration, in MEGA 4.0, Tamura et al. (2007), and Microsoft 2007.
Haplotype network
Haplotype sequences were obtained with the DnaSP software, excluding the aligned sites with gaps, missing data, and monomorphic sites. The Network 4.6.1.1 (Fluxus Technology, 2004-2012) software was used to build haplotype networks (Bandelt et al., 1999; Polzin et al., 2003) with Lennoa sequences as base group.
Haplotype phylogeny
The phylogenetic tree was generated with MEGA 4.0 (Tamura et al., 2007), using the Neighbor-Joining method, with a bootstrap of one thousand iterations. Genetic distances were calculated with the Maximum Composite Likelihood method, without taking into consideration indels. Lennoa sp. sequences were used as an outgroup.
Results and discussion
Sampling
Our P. culiacanum distribution area scouting found seven populations: five in Sinaloa and two in Sonora. Out of all the collected individuals, 10 of each population were selected for DNA purification. In some reported distribution sites (GIF, 2013), no populations were found, due to changes in land use - such as primary forests transformed into agricultural or forest exploitation areas (Ranchería Los Álamos, Cerro Siete Gotas, in Culiacán, San Blas)-, to the construction of highways -highway in south Culiacán, Los Mochis-Obregón highway, AltataDautillos highway-, or to tourism and real estate development projects (Nuevo Altata costal area), and the consequent ecosystem fragmentation.
Amplification
Out of all the sequenced samples, 70 belonged to P.culiacanum and two belonged to Lennoa madreporoides. The size of the amplified fragment was 735 nucleotides, 261 had gaps or were missing data, 151 were monomorphic sites and 323 were polymorphic sites. Two-hundred fifty-one of those polymorphic sites presented two types of substitutions, 69 presented three types, and 3 presented four types. The length of the amplification product, the unamplified regions (rps432-trnL [UAA], trnT [GGU]-trnE [UUC], atpB-rbcL), and the amplified regions with low variation levels (trnT [UGU] - trnL [UAA], rpl32F - trnL [UAG] y trnT [UGU] - trnL [UAA]) all suggest that in P. culiacanum, as in other parasitic plants, the size of the chloroplast was reduced, in comparison with the original size in tobacco plants (1560 nucleotides) (Bungard, 2004).
Variability analysis
The seven analyzed populations are polymorphic, with 11 different haplotypes (2-4 per population), without taking indels into consideration as informative sites (Figure 1).
Nucleotide total diversity (π) -without taking indels into consideration- was 0.2737 ± 0.0242; Masiaca and Tosalibampo had the highest and lowest values (0.2990, and 0.0007, respectively). Overall haplotype diversity -without taking indels (Hd) into consideration- was 0.855 ± 0.024; Altata and Cerro de Álamos shared the highest value (0.7330), while Tosalibampo had the lowest (0.4670) (Table 1).
Localidad | h | Hd | hi | Hdi | π |
Altata | 3 | 0.73300 | 10 | 1.00000 | 0.25957 |
Culiacán | 4 | 0.72200 | 8 | 0.97222 | 0.12660 |
Tosalibampo | 2 | 0.46700 | 9 | 0.97778 | 0.00069 |
San Blas | 3 | 0.60000 | 9 | 0.97778 | 0.28244 |
Heraclio SANB | 3 | 0.66700 | 6 | 0.84444 | 0.26463 |
Masiaca | 2 | 0.53600 | 8 | 1.00000 | 0.29903 |
Cerro de Alamos | 4 | 0.73300 | 10 | 1.00000 | 0.2960 |
Total | 11 | 0.855 ± 0.024 | 57 | 0.99367 | 0.27366 ± 0.02425 |
Number of haplotypes (h), haplotype diversity (Hd), number of haplotypes, taking indels into consideration (hi), haplotype diversity, taking indels into consideration (Hdi), nucleotide diversity (π), standard deviation (±).
First value (0.1349, p≤0.05) indicates a moderate level of genetic differentiation among studied populations. AMOVA showed that most of the observed genetic variation is mainly due to differences within populations (86.51 %), rather than differences between populations (13.49 %). This result is consistent with those of other parasitic plants studies (Jerome et al., 2002; Amico et al., 2009) (Table 2). This was to be expected, because the populations’ Hd values are over 0.5 (except in Tosalibampo); high levels of molecular variance within populations result in high haplotipic diversity. High variation within populations can be interpreted as a consequence of gene flow between individuals of the same population (García-Franco et al., 1998). This contrasts with other studies on parasitic and dependent plants, in which most genetic variation was found at population level. These species are characterized by: fragmented or patchy distribution; dependence on their hosts (de Vega et al., 2008); the life history of both the parasite and the host; suffering geographic and stochastic events (bottle necks and founder events). However, there are species that are wind-pollinated or that use the wind to spread their seeds (such as P. culiacanum). These species usually have a lower genetic diversity index between populations (Allendorf et al., 2014).
AMOVA Método de distancia: Tajima & Nei | ||||
Variación | G.L | Suma de cuadrados | Varianza | Porciento de varianza |
Entre poblaciones | 6 | 2247.257 | 23.44503 Va | 13.49 |
Dentro de poblaciones | 60 | 9021.418 | 150.3569 Vb | 86.51 |
Total | 66 | 11268.675 | 173.80201 |
Genetic differentiation index Fst=0.13490. † Significance test 1023 permutations. ᾽Va and FST: P value= 0.0137±0.0034. §Weir, B.S. and Cockerham, C.C. 1984. Excoffier, L., Smouse, P., and Quattro, J. 1992. Weir, B. S., 1996.
Although population structure had low levels of genetic differentiation among populations, the Pearson correlation coefficient and the linear regression model for genetic distances showed no evidence of isolation by distance (p>0.05 y r2=0.0026), unlike in other parasitic plants (Holzapfel et al., 2002; Jerome et al., 2002a; Mutikainen et al., 2002; deVega et al., 2008). This allowed us to infer that there is a gene flow amongst the studied populations and that overall they behave as a single evolutionary population unit (García-Franco et al., 1998; Allendorf et al., 2014), which is large enough to avoid the effects of genetic drift. However, in rare species -such as P. culiacanum-, where variation is equal or larger than in common species, evaluating other non-ecologic factors -such as life cycle strategies, reproductive attributes, demography, interaction with other organisms (the host in this case), and others factors- is advised (Eguiarte et al. 2007).
Haplotype distribution
Haplotype 5 had the highest frequency (19), and it was found in Tosalibampo, Masiaca, San Blas and Cerro de Álamos; haplotype 6 had the highest distribution rate, and it was found in Culiacán, San Blas, Heraclio Sb, Masiaca and Cerro de Álamos (Figure 2).
Populations with unique haplotypes were Altata (haplotype 1), Culiacán (haplotypes 3, 9, and 11), Tosalibampo (haplotype 10), and Cerro de Álamos (haplotypes 2 and 8). The other three populations (Masiaca, San Blas, and Heraclio Sb) were made up entirely by haplotypes present in other populations (haplotypes 4, 5, 6, and 7). Unshared haplotypes can be the result of mutations related to the colonization of new ecological niches (hosts, ecosystems, or both). This P. culiacanum ability to parasitize different plant genera in different ecosystems was documented by Yatskievych (1985). It can also be a result of the expansion of hosts to new territories and the subsequent colonization of the parasite as proposed by deVega et al. (2008).
Haplotype networks
In Figure 2, each haplotype is represented by a black circle the size of which is proportional to the frequency of the sequences that it is made up by. Mean vectors are represented in white and were generated by the software to link distant haplotypes. Mean vectors can be mutational steps, extinct, or not sampled subspecies. The analysis detected three major haplotype groups: 1) haplotypes 1 and 4, which are 51 mutational steps away from Lennoa (H12); 2) a group formed by haplotypes 2, 3, 5, 9, 10 and 11 is 160 steps away ; and 3) another group formed by haplotypes 6, 7 and 8 is 207 steps away. The methods used in haplotype networks are designed to model the evolutionary past via a stochastic process (coalescence), which is based on the idea that the alleles in a population can be traced back in time to the point where they coalesce with an ancestral allele. Therefore, the haplotypes furthest from the root sequence (6, 7, 8, 3, 11 y 9) are the most recent and the haplotypes closest to the root sequence (1 and 4) are the oldest (Eguiarte et al., 2007).
Haplotype phylogeny
The phylogenetic tree of the 11 haplotypes showed four groups: 1) the outgroup group (Lennoa); 2) haplotypes 1 and 4, as the basal group of P. culiacanum and the closest to Lennoa; 3) a branch which bifurcates into two major groups (80 % support), on the one hand, haplotypes close to the basal group 6, 8 and 7 (100 % support); and (4) a branch (100 % support) with haplotypes 2, 3, 5, 9, 10, 11. The general distribution matched the haplotype network diagram, including the proximity or distance from the root sequences (Figure 3).
Populations’ morphology varies geographically. There were two major differentiated groups. The first one includes the populations of Altata and Heraclio Sb, and the second is made up by the populations of Culiacán, Tosalibampo, San Blas, Masiaca, and Cerro de Álamos. The former has characteristic thin or branched stems originating from a single haustorium, similar to those of Lennoa. The second group has one single thick stem, over 1 cm in diameter. These observed characteristics matched the molecular evidence that resulted from the phylogenetic analysis and the haplotype network, which located most of group 1 haplotypes far apart from group 2 haplotypes.
The haplotype networks showed several mutational steps between those haplotypes that are mostly formed by individuals from Altata and Heraclio Sb (haplotypes 1, 4, 6, 7 and 8) and those from all other populations. Haplotype phylogeny also located these morphologically different populations’ haplotypes close to each other and to Lennoa. Other relevant characteristics were the different kind of soils in which the populations were found, ranging from coastal dunes (Altata) and sandy soil (Heraclio Sb), to clay soil in deciduous forests and the ecosystems (thorn forest, xeric shrublands, deciduous and semideciduous forests). Pholisma arenarium and Pholisma sonorae share these types of soil (Nabhan, 1979; Felger, 1980; Yatskievych, 1985; Yatskievych and Mason, 1986; AGFD Plant Abstract, 2004), which might suggests that these P. culiacanum populations belong to a different species of the Pholisma genus or to a P. culiacanum subspecies. For this reason, a comparative molecular analysis of the three species of the genus is necessary to determine the reason for these differences.
Given the high levels of haplotype genetic diversity, it is advisable to establish genetic resources conservation areas in the Culiacán, Cerro de Álamos and Altata ranges. These would include almost all the genetic diversity of the studied populations and most of the unique haplotypes would be found in the preserved areas. This is important because the aforementioned ranges are the most affected by anthropogenic activities, such as urban and touristic developments, agriculture, forestry, and land-use changes.
The Tosalibampo range is recommended for in situ exploitation, because the genetic diversity and haplotype indexes were the lowest; only two haplotypes were found and it had the most abundant population.
Literatura citada
Amani, J., K. Roohallah, R. A. Ali, and H. S. Ali. 2011. A simple and rapid leaf genomic DNA extraction method for polymerase chain reaction analysis. Iran. J. Biotechnol. 9: 69-71. [ Links ]
Amico G. C., R. Vidal-Russell, and D. L Nickrent. 2007. Phylogenetic relationships and ecological speciation in the mistletoe Tristerix (Loranthaceae): The influence of pollinators, dispersers, and hosts. Am. J. Bot. 94: 558-567. [ Links ]
Amico, G. C., and D. L Nickrent. 2009. Population structure and phylogeography of the mistletoes Tristerix corymbosus and T. aphyllus (Loranthaceae) using chloroplast DNA sequence variation. Am. J. Bot. 96: 1571-1580. [ Links ]
Allendorf, F. A., G. Luikart, and S. N. Aitken. 2014. Conservation and the Genetics of Populations. Wiley-Blackwell. 2nd ed. 624 p. [ Links ]
Arizona Game and Fish Department. 2004. Lennoaceae. PDLNN02020. Unpublished abstract compiled and edited by the Heritage Data Management System, Arizona Game and Fish Department, Phoenix, AZ. 3 p. [ Links ]
Azmat, M. A., A. K. Iqrar, M. N. Hafiza C., A. R. Ishtiaq, S. K. Ahmad, and A. H. Asif. 2012. Extraction of DNA suitable for PCR applications from mature leaves of Mangifera Indica L. J. Zhejiang University-Science B. Biomed. Biotechnol. 13: 239-243. [ Links ]
Bandelt, H., P. Forster, and A. Röhl. 1999. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 16: 37-48. [ Links ]
Bungard, R. A. 2004. Photosynthetic evolution in parasitic plants: insight from the chloroplast genome. BioEssays 26: 235-247. [ Links ]
CodonCode Corporation, 2009. CodonCode Aligner versión 3.5. http://www.codoncode.com/ . (Consulta: Febrero 2010). [ Links ]
Cota-Sánchez, H. J., K. Remarchuk, and U. A. Kumary. 2006. Ready-to-use DNA extracted with a CTAB method adapted for herbarium specimens and mucilaginous plant tissue. Plant Mol. Biol. Rep. 24: 161-167. [ Links ]
dePamphilis C. W., D. Y. Nelson, and A. D. Wolfe. 1997. Evolution of plastid gene rps2 in a lineage of hemiparasitic and holoparasitic plants: Many losses of photosynthesis and complex patterns of rate variation. PNAS 94: 7367-7372. [ Links ]
deVega, C., R. Berjano, M. Arista, P. L. Ortiz, T. Salvador, and T. F. Stuessy. 2008. Genetic races associated with the genera and sections of host species in the holoparasitic plant Cytinus (Cytinaceae) in the western mediterranean basin. New Phytol. 178: 875-887. [ Links ]
Doyle, J. J. 1991. DNA Protocols for plants. In: Hewitt, G., A. W. B. Johnson, and J. P. W. Young (eds), Molecular Techniques In Taxonomy. NATO ASI Series H, Cell Biol. 57: 283-293. [ Links ]
Dressler, R. L. 1968. A second species of Ammobroma (Lennoaceae) in Sinaloa, Mexico. Madroño 19: 179-182. [ Links ]
Excoffier, L., and S. Schneider. 2005. Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evol. Bioinformatics Online 1: 47-50. [ Links ]
Funk, H. T., S. Berg, K. Krupinska, U. G. Maier, and K. Krause. 2007. Complete DNA sequences of the plastid genomes of two parasitic flowering plant species, Cuscuta reflexa and Cuscuta gronovii. Plant Biol. 45: 1-12. [ Links ]
Garcia-Franco, J. G., V. Souza, L. E. Eguiarte, and V. Rico-Gray. 1998. Genetic variation, genetic structure and effective population size in tropical holoparasitic endophyte Bdallophyton bambusarum (Rafflesiaceae). Plant Syst. Evol. 210: 271-288. [ Links ]
Hall, T. A. 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/ NT. Nucleic Acids Symp. Ser 41: 95-98. [ Links ]
Hartl, D. L., and A. G. Clark. 2007. Principles of Population Genetics. 4ta ed. Ed.Sinauer Associates, Inc. Publishers. Sunderland, Massachusetts. USA. 652 p. [ Links ]
Hedrick, P. W. 2005. Genetic of Populations. Jones and Bartlett. Third Edition. USA. 725 p. [ Links ]
Heinze, B. 2007. A database of PCR primers for the chloroplast genomes of higher plants. Plant Method 3: 1-7. [ Links ]
Holzapfel, S., M. Z. Faville, and C. E. Gemmill. 2002. Genetic variation of the endangered holoparasite Dactylanthus taylorii (Balanophoraceae). J. Biogeogr. 29: 663-676. [ Links ]
Jerome, C. A., and B. A. Ford. 2002-a. The discovery of three genetic races of the dwarf mistletoe Arceuthobium americanum (Viscaceae) provides insight into the evolution of parasitic angiosperms. Mol. Ecol. 11: 387-405. [ Links ]
Krizman, M., J. Jernej, B. Dea, J. Branka, and P. Mirko. 2006. Robust CTAB-activated harcoal protocol for plant DNA extraction. Acta Agriculturae Slovenica 87: 427-433. [ Links ]
Librado, P., and J. Rosas. 2009. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25: 1451-1452. [ Links ]
Mutikainen, P., and T. Koskela. 2002. Population structure of a parasitic plant and its perennial host. Heredity 89: 318-324. [ Links ]
Nabhan G. P. 1980. The ecology of floodwater farming in southwestern North America. Agroecosy. 5: 245-55. [ Links ]
Nickrent, D. L. and M. A. García. 2009. On the brink of holoparasitism: plastome evolution in dwarf mistletoes (Arceuthobium, Viscaceae). J. Mol. Evol. 68: 603-615. [ Links ]
Zuber, D., and A. Widmer. 2000. Genetic evidence for host specificity in the hemi-parasitic Viscum album L. (Viscaceae). Mol. Ecol. 9: 1069-1073. [ Links ]
Polzin, T., and S. V. Daneschmand. 2003. On Steiner trees and minimum spanning trees in hypergraphs. Oper. Res. Lett. 31:12-20. [ Links ]
Porebski, S., Grant B., and R. B. Bernard. 1997. Modification of a CTAB DNA extraction protocol for plants containing high polysaccharide and polyphenol components. Plant Mol. Biol. Rep. 15: 8-15. [ Links ]
Sá, O., J. A. Pereira and P. Baptista. 2011. Optimization of DNA extraction for RAPD and ISSR analysis of Arbutus unedo L. leaves. Int. J. Mol. Sci. 12: 4156-4164. [ Links ]
Sánchez-Hernández, C., and J. C. Gaytán-Oyarzún. 2006. Two mini-preparation protocols to DNA extraction from plants with high polysaccharide and secondary metabolites. African J. Biotech. 5: 1864-1867. [ Links ]
Taberlet, P., L. Gielly, G. Pautou, and J. Bouvet. 1991. Universal primers for amplification of three non-coding regions o chloroplast DNA. Plant Mol. Biol. 17:1105-1109. [ Links ]
Tamura, K., J. Dudley, M. Nei, and S. Kumar. 2007. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol. Biol. Evol. 24: 1596-1599. [ Links ]
Yatskievych, G. 1985. Notes on the biology of the Lennoaceae. Cactus Succulent 57: 73-79. [ Links ]
Yatskievych, G. and C. T Mason. 1986. A revision of the Lennoaceae. Syst. Bot. 11: 531-584. [ Links ]
Received: June 2015; Accepted: March 2016