Print version ISSN 0100-204X
Pesq. agropec. bras. vol.47 no.8 Brasília Aug. 2012
Genetic diversity of cultivated accessions and wild species of rubber tree using EST‑SSR markers
Diversidade genética de acessos cultivados e espécies silvestres de seringueira por meio de marcadores EST‑SSR
Juliana Morini Küpper Cardoso PerseguiniI; Lineu Roberto de Castro RomãoI; Boris BriñezI; Erivaldo José Scaloppi JuniorII; Paulo de Souza GonçalvesIII; Luciana Lasry BenchimolIII
IUniversidade Estadual de Campinas, Campus Universitário Zeferino Vaz, s/no, Cidade Universitária, CEP 13083‑970 Campinas, SP, Brazil, E‑mail: firstname.lastname@example.org, email@example.com, firstname.lastname@example.org
IIAgência Paulista de Tecnologia dos Agronegócios, Rodovia Péricles Belini, Km 121, Caixa Postal 61, CEP 15500‑970 Votuporanga, SP, Brazil, E‑mail: email@example.com
IIIInstituto Agronômico de Campinas, Avenida Barão de Itapura, no 1.481, Jardim Guanabara, Caixa Postal 28, CEP 13012‑970 Campinas, SP, Brazil, E‑mail: firstname.lastname@example.org, email@example.com
The objective of this work was to evaluate the efficiency of EST‑SSR markers in the assessment of the genetic diversity of rubber tree genotypes (Hevea brasiliensis) and to verify the transferability of these markers for wild species of Hevea. Forty‑five rubber tree accessions from the Instituto Agronômico (Campinas, SP, Brazil) and six wild species were used. Information provided by modified Roger's genetic distance were used to analyze EST‑SSR data. UPGMA clustering divided the samples into two major groups with high genetic differentiation, while the software Structure distributed the 51 clones into eight groups. A parallel could be established between both clustering analyses. The 30 polymorphic EST‑SSRs showed from two to ten alleles and were efficient in amplifying the six wild species. Functional EST‑SSR microsatellites are efficient in evaluating the genetic diversity among rubber tree clones and can be used to translate the genetic differences among cultivars and to fingerprint closely related materials. The accessions from the Instituto Agronômico show high genetic diversity. The EST‑SSR markers, developed from Hevea brasiliensis, show transferability and are able to amplify other species of Hevea.
Index terms: fingerprinting analysis, genetic structure, functional molecular markers, polymorphism information content, transferability.
O objetivo deste trabalho foi avaliar a eficiência de marcadores EST‑SSR na determinação da diversidade genética de genótipos de seringueira e verificar a transferibilidade destes marcadores para espécies silvestres de Hevea. Foram utilizados 45 acessos de seringueira (H. brasiliensis) do Instituto Agronômico e seis espécies silvestres. As informações fornecidas pela distância genética de Roger modificada foram usadas para analisar os dados de EST‑SSR. O agrupamento UPGMA dividiu as amostras em dois grandes grupos com alta diferenciação genética, enquanto o programa Structure distribuiu os 51 clones em oito grupos. Foi possível traçar um paralelo entre ambos os métodos de agrupamento. Os 30 EST‑SSRs polimórficos mostraram de dois a dez alelos e foram eficientes em amplificar as seis espécies silvestres. Microssatélites funcionais EST‑SSR são eficientes na avaliação da diversidade genética entre clones de seringueira e podem ser usados para traduzir diferenças genéticas entre cultivares e para gerar perfis genéticos de materiais próximos. Os acessos do Instituto Agronômico apresentam elevada diversidade genética. Os marcadores EST‑SSR, desenvolvidos para Hevea brasilensis, apresentam transferabilidade e são capazes de amplificar outras espécies de Hevea.
Termos para indexação: análise de caracterização genética, estrutura genética, marcadores moleculares funcionais, conteúdo de informação polimórfica, transferabilidade.
The rubber tree [Hevea brasiliensis (Willd. ex. Adr. de Juss.) Muell. Arg.] is native to the Amazonian rainforest and the most widely cultivated tree species for producing natural rubber tree latex, which has been an undeniably beneficial commodity for the past 100 years (Priyadarshan & Gonçalves, 2003). Despite being the center of origin of the rubber tree and the main producer and exporter at the end of the 19th century, Brazil started to import this raw material at the beginning of the last century. Data from the International Rubber Study Group (2011) showed that, in 2010, the world production of natural rubber reached 10,291 thousand tons, of which Brazil only contributed with 132 thousand tons, i.e., about 1% of the total.
Cultivated rubber tree clones are characterized by low genetic variability. According to Yu et al. (2011), the genetic basis of Hevea breeding programs has to be expanded in order to improve the breeding efficiency and to accelerate the breeding of new cultivars. However, the genetic base is decreasing and little improvement in the productivity of new cultivars is being achieved because of the long period of time required for interspecific breeding among related species.
The rubber tree is a perennial cross‑pollinating and monoecious species of the Euphorbiaceae family, which has a long breeding process, with many selection cycles and difficulties in raising F2 progenies. These difficulties are attributed to several factors related to the reproductive biology of rubber, such as lack of synchrony in flowering among clones, seasonality in flowering, pollen longevity, low pollination success, and seed recalcitrance (Gonçalves & Marques, 2008). Therefore, conventional genetic analysis in this crop is difficult and time consuming (Saha et al., 2007).
Molecular markers have been used in several studies to assess genetic diversity among rubber tree genotypes (Roy et al., 2004; Saha et al., 2005). Microsatellites, or simple sequence repeats (SSRs), combine co‑dominance, high polymorphism, good genome coverage, and Mendelian inheritance. Up to now, they have also been capable of discriminating closely related individuals (Brondani et al., 1998). Reports confirm the successful use and application of SSRs in rubber tree breeding (Feng et al., 2009; Le Guen et al., 2009, 2011; Gouvêa et al., 2010).
The presence of SSRs in the transcripts of genes indicates that they might have a role in gene expression or function; however, it remains to be seen whether any unusual phenotypic variation might be associated with the length of SSRs in coding regions (Varshney et al., 2005). This fact, together with the increasing availability of sequences in databases, mainly expressed sequence tags (ESTs), has enabled the development of functional microsatellite markers or EST‑SSRs. The use of these markers has been reported in several species, including Hevea sp. (Feng et al., 2009). The fact that these markers are derived from genes and that, consequently, their sequences are more conserved, make it extremely robust to identify polymorphisms in closely related species (Varshney et al., 2005).
The objective of this work was to evaluate the efficiency of EST‑SSR markers in assessing the genetic diversity of 51 rubber trees and to verify the transferability of these markers for six wild species of Hevea.
Materials and Methods
A total of 45 rubber tree genotypes (H. brasiliensis) from the Instituto Agronômico (IAC, Campinas, SP, Brazil) and six wild species of Hevea (H. guianensis, H. rigidifolia, H. nitida, H. pauciflora, H. benthamiana, and H. camargoana), provided by Embrapa Amazônia Ocidental (Manaus, AM, Brazil), were used (Table 1). Total genomic DNA was extracted from powdered lyophilized young leaves using the CTAB method (Hoisington et al., 1994). The IAC genotypes were from Asiatic, African, Amazonian, and IAC series. Several of the Asiatic clones were derived from the Wickham collection originally introduced into Asia in 1876, which were known as the Wickham clones (Wycherley, 1968). The Amazonian clones were derived from the selection and crossings carried out in Brazil by the Ford Motor Company and the Instituto Agronômico do Norte (IAN), which is currently the Instituto de Pesquisa Agropecuária do Norte. The IAC genotypes resulted from controlled crossings and open pollinations performed by IAC.
A total of 30 EST‑SSRs (Table 2) were selected from Feng et al. (2009). These markers were derived from a cDNA library developed to analyze gene expression in laticifers by Chow et al. (2007). Amplifications were performed in a final volume of 25 µL containing 50 ng DNA, 1X buffer, 0.2 µmol L‑1 of each forward and reverse primer, 100 µmol L‑1 of each dNTP, 2.0 mmol L‑1 MgCl2, 10 mmol L‑1 Tris‑HCl (pH 8.0), 50 mmol L‑1 KCl, and 0.5 U Taq DNA polymerase. The following conditions were used for amplification: 1 min at 94ºC, then 30 cycles of 1 min at 94ºC, 1 min at the specific annealing temperature for each SSR, and 1 min at 72ºC, with a final extension of 5 min at 72ºC. The PCR products were separated with 6% silver‑stained polyacrylamide gel electrophoresis. Molecular size of the DNA fragments was estimated by comparison with standard 10‑bp DNA ladder (Invitrogen, São Paulo, SP, Brazil). A binary matrix was constructed from the gel readings, in which the individuals were genotypically characterized for presence (1) and absence (0) of bands. The percentage of polymorphism obtained with each primer was calculated from this matrix.
Genetic distances (GDs) were calculated from the EST‑SSR data for all possible inbred pairs using modified Roger's genetic distance (MRD) (Goodman & Stuber, 1983), version 1.3 (Miller, 1997). Cluster analyses were performed using the unweighted pair‑group method with arithmetic averages (UPGMA).
Polymorphism information content (PIC) values were calculated using the formula,
in which: n is the number of alleles and fi and fj are the frequencies of the ith and jth alleles, respectively (Lynch & Walsh, 1998). Discrimination power (DP) values for the kth primer were calculated using the formula,
in which: N is the number of individuals and pj is the frequency of the jth pattern (Tessier et al., 1999). PIC was used to measure the information of a given marker locus for the pool of genotypes, whereas DP was used to measure the efficiency of SSRs in order to identify varieties by taking into account the probability of two randomly chosen individuals having different patterns. In addition, Wright F statistics and molecular analysis of variance (AMOVA) were determined with Arlequin v. 3.11 (Excoffier et al., 2005).
Bayesian clustering was performed using the software Structure 2.2 (Pritchard et al., 2000). The number of clusters was defined from K = 2 to K = 20, and ten runs of each K were carried out using: the admixture model and correlated allele frequencies, a 200,000 burn‑in period, and 500,000 MCMC. Ad hoc statistics were related to rate changes in the log probability of data according to the number of Ks proposed by Evanno et al. (2005), with ΔK being used as a predictor of the ideal number of clusters.
Results and Discussion
All of the thirty EST‑SSRs were polymorphic (Table 2) and produced a total of 149 polymorphic alleles. The number of alleles ranged from two to ten, with a mean of 4.96 alleles per marker. The highest number of alleles was observed for the microsatellites HBE001, HBE117, HBE126, and HBE167.
The highest value of PIC found for EST‑SSRs was 0.86 (HBE117) and the lowest was 0.23 (HBE200), with a mean value of 0.59. Souza et al. (2009) assessed the genetic diversity of H. brasiliensis and six wild species of rubber tree using genomic SSR markers and reported a similar range for PIC values (0.13 to 0.88).
Cross‑species amplification with the 30 EST‑SSRs was successful for the six wild species evaluated in the present study. Saha et al. (2005) also observed that SSRs specifically developed for H. brasiliensis efficiently amplified H. benthamiana and H. spruceana. Other studies indicated that SSRs developed for H. brasiliensis were successfully transferable to Hevea wild species (Souza et al., 2009; Gouvêa et al., 2010; Le Guen et al., 2011).
The DP analyses (Table 2) showed values ranging from 0.07 (HBE017) to 0.90 (HBE021). According to the number of alleles, discrimination power, and polymorphic information content, the polymorphism estimated for the rubber tree with EST‑SSRs indicated that these loci were capable of generating sufficient information to ensure their use to establish genetic relationships among closely related materials. These markers could also be useful for a wide range of genetic investigations, such as linkage map construction and association mapping studies.
Genotypes (Figure 1) showed high genetic variability, as genetic distances varied from 0.35 to 0.77, and a high genetic structure was observed from the clustering pattern. Similar genetic distances (0.26 to 0.67) were found by Gouvêa et al. (2010), who reported high genetic variability among rubber tree genotypes. The 51 clones were separated into two major groups. The first group was composed mostly of IAC clones, developed within the IAC breeding program. Fewer Amazonian, Asian, and African clones were also included in this group. The second group consisted of Malaysian clones from the Rubber Research Institute of Malaysia (RRIM). This group showed a clear separation between wild Hevea species, in comparison to other H. brasiliensis genotypes. These samples have undergone several breeding cycles, as confirmed by the highest interspecific genetic distances observed. The smallest genetic distances were intraspecific, belonging to genotypes derived from the same breeding institution (Group 1: PB, Prang Besar; Group 2: PB and RRIM) (Figure 1). Moreover, IAC clones (Group 1) showed greater genetic distances than Malaysian clones (Group 2), indicating a greater genetic variability among Brazilian genotypes, which could be explained by the selection pressure present in each breeding program.
The overall observed heterozygosity (Ho) assessed for group 1 was 0.25883, and the expected heterozigosity (He) was 0.60264; for group 2, Ho was 0.27090 and He was 0.51366. The pair‑wise FST (Weir & Cockerham, 1984) was high (27.64%), indicating great variability among the dendrogram groups (FST). Furthermore, 72.36 % of variation was contained within the populations (FIS).
The cophenetic correlation coefficient was high (r=0.70, p<0.002), indicating that UPGMA clustering accurately represented the genetic distances among rubber tree clones. Gouvêa et al. (2010) observed a similar cophenetic value for the genomic SSR UPGMA clustering. Odong et al. (2011) found a relationship between the cophenetic correlation coefficient (CPCC) and subgroup differentiation (FST), and that a high CPCC with UPGMA clustering is an indication of a reliable population structure in the data.
The EST‑SSRs used in the present study were efficient in promoting the distinction of rubber tree accessions that have common male and female parents, as observed by the following clones: RRIM 710 and RRIM 713, whose parents were RRIM 605 and RRIM 71; RRIM 908 and RRIM 911, derived from the PB 551 and RRIM 623 intersection; and RRIM 915 and RRIM 919, whose parents were RRIM 605 and PB 5/51 (Figure 1).
The clones PB 235 (false) and PB 235 (true) were differentiated (Figure 1) by the EST‑SSRs, proving that these markers can be extremely useful for conducting effective fingerprinting analysis. In fact, fingerprinting was applied to verify the identity of parental clones during micropropagation and to certify the purity of F1 hybrids. The use of molecular markers to perform fingerprinting analysis was reported in other studies (Caruso et al., 2008; Van Treuren et al., 2010), with the same efficiency, to separate phenotypically similar genotypes, such as the two PB 235 clones used in the present study.
The EST‑SSRs, despite being developed from H. brasiliensis, were able to amplify other species of Hevea. This transferability is very important in the sense that it enables the use of microsatellite analysis in intra‑ and interspecific genetic diversity and gene flow evaluation, as well as phylogenetic relationships between cultivated clones of H. brasiliensis and of other species of wild Hevea. Souza et al. (2009) used genomic and microsatellite markers and also observed transferability of loci with wild species. Mantello et al. (2012) found a percentage of transferability from 82 to 87%.
According to Yu et al. (2011), there are homogeneous microsatellite sequences within species, between species, and across genes. All evidence indicates that, with large genetic distances, the ability to amplify microsatellite loci with the same primer pairs decreases. These authors found a good discriminatory power for genomic microsatellites used to investigate genetic diversity and that these loci were similar to the EST‑SSRs reported by Feng et al. (2009). Indeed, EST‑SSRs were efficiently used to access the genetic diversity of H. brasiliensis, being capable of dividing the genotypes into two major and clearly structured groups.
A total of eight groups were identified by ΔK as being the ideal number of groups, according to Evanno et al. (2005). Feng's SSRs were used to establish a correspondence between the dendrogram and Structure clustering (Figures 1 and 2). It was observed that group 1 of the dendrogram corresponded entirely to Structure groups I, II, IV, and VI. Group 2 of the dendrogram corresponded to Structure groups III and VII, whereas group 3 corresponded entirely to group VI of the Structure clustering. Therefore, Structure analysis was also efficient to access the genetic organization of Hevea species, although it tended to generate a deeper subgroup differentiation. The consistency in genotype clusters indicated a nonrandom distribution of alleles and their frequencies. A parallel could be established between both analyses, but UPGMA clustering was clearer. According to Odong et al. (2011), traditional cluster analysis, such as UPGMA, can provide an easier and more effective way of determining structure in germplasm collections, using molecular marker data, in comparison to Bayesian methods.
Microsatellites derived from genic sequences (EST‑SSR) have additional advantages because they can be developed from regions of the genome associated with a trait of interest and could be part of the gene controlling the character (Thiel et al., 2003). Another advantage is the fact that the gene regions are more conserved, when compared to anonymous ones. Therefore, these functional markers are likely to be more easily transferred to related species (Varshney et al., 2005).
1. The EST‑SSRs are efficient to evaluate genetic diversity among rubber tree accessions.
2. EST‑SSRs can be used to translate the genetic differences among cultivars and fingerprint closely related materials.
3. The EST‑SSRs developed from Hevea brasiliensis show transferability and are able to amplify other species of Hevea.
4. The accessions from Instituto Agronômico (IAC) show high genetic diversity.
To Fundação de Amparo à Pesquisa do Estado de São Paulo and to Conselho Nacional de Desenvolvimento Científico e Tecnológico, for scholarships granted.
BRONDANI, R.P.V.; BRONDANI, C.; TARCHINI, R.; GRATTAPAGLIA, D. Development, characterization and mapping of microsatellite markers in Eucalyptus grandis and E. uruphylla. Theoretical and Applied Genetics, v.97, p.816‑827, 1998. [ Links ]
CARUSO, M.; FEDERICI, C.T.; ROOSE, M.L. EST‑SSR markers for asparagus genetic diversity evaluation and cultivar identification. Molecular Breeding, v.21, p.195‑204, 2008. [ Links ]
CHOW, K.‑S.; WAN, K.‑L.; ISA, M.N.M.; BAHARI, A.; TAN, S.H.; HARIKRISHNA, K.; YEANG, H.Y. Insights into rubber biosynthesis from transcriptome analysis of Hevea brasiliensis latex. Journal of Experimental Botany, v.58, p.2429‑2440, 2007. [ Links ]
EVANNO, G.; REGNAUT, S.; GOUDET, J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Molecular Ecology, v.14, p.2611‑2620, 2005. [ Links ]
EXCOFFIER, L.; LAVAL, G.; SCHNEIDER, S. Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evolutionary Bioinformatics, v.1, p.47‑50, 2005. [ Links ]
FENG, S.P.; LI, W.G.; HUANG, H.S.; WANG, J.Y.; WU, Y.T. Development, characterization and cross‑species/genera transferability of EST‑SSR markers for rubber tree (Hevea brasiliensis). Molecular Breeding, v.23, p.85‑97, 2009. [ Links ]
GONÇALVES, O.S.; MARQUES, J.R.B. Melhoramento genético da seringueira: passado, presente e futuro. In: ALVARENGA, A. de P.; CARMO, C.A.F. de S. do (Coord.). Seringueira. Viçosa: EPAMIG, 2008. p.399‑498. [ Links ]
GOODMAN, M.M.; STUBER, C.W. Races of maize. VI. Isozyme variation among races of maize in Bolivia. Maydica, v.28, p.169‑187, 1983. [ Links ]
GOUVÊA, L.R.L.; RUBIANO, L.B.; CHIORATTO, A.F.; ZUCCHI, M.I.; GONÇALVES, P. de S. Genetic divergence of rubber tree estimated by multivariate techniques and microsatellite markers. Genetics and Molecular Biology, v.33, p.308‑318, 2010. [ Links ]
HOISINGTON, D.; KHAIRALLAH, M.; GONZALEZ DE LEÓN, D. Laboratory protocols: CIMMYT applied molecular genetics laboratory. 2nd ed. Mexico: CIMMYT, 1994. 102p. [ Links ]
LE GUEN, V.; DOARÉ, F.; WEBER, C.; SEGUIN, M. Genetic structure of Amazonian populations of Hevea brasiliensis is shaped by hydrographical network and isolation by distance. Tree Genetics and Genomes, v.5, p.673‑683, 2009. [ Links ]
LE GUEN, V.; GAY, C.; XIONG, T.C.; SOUZA, L.M.; RODIER‑GOUD, M.; SEGUIN, M. Development and characterization of 296 new polymorphic microsatellite markers for rubber tree (Hevea brasiliensis). Plant Breeding, v.130, p.294‑296, 2011. [ Links ]
LYNCH, M.; WALSH, J.B. Genetics and analysis of quantitative traits. Sunderland: Sinauer, 1998. 980p. [ Links ]
Mantello, C.C.; SUZUKI, F.I.; Souza, L.M.; Gonçalves, P. de S.; SOUZA, A.P. Microsatellite marker development for the rubber tree (Hevea brasiliensis): characterization and cross‑amplification in wild Hevea species. BMC Research Notes, v.5, 2012. Doi:10.1186/1756‑0500‑5‑329. [ Links ]
MILLER, M. Tools for population genetic analyses (TFPGA). Version 1.3. Flagstaff: Northern Arizona University, 1997. 33p. [ Links ]
ODONG, T.L.; HEERWAARDEN, J. van; JANSEN, J.; HINTUM, T.J.L. van; EEUWIJK, F.A. van. Determination of genetic structure of germplasm collections: are traditional hierarchical clustering methods appropriate for molecular marker data? Theoretical and Applied Genetics, v.123, p.195‑205, 2011. [ Links ]
PRITCHARD, J.K.; STEPHENS, M.; DONNELLY, P. Inference of population structure using multilocus genotype data. Genetics, v.155, p.945‑959, 2000. [ Links ]
PRIYADARSHAN, P.M.; GONÇALVES, P. de S. Hevea gene pool for breeding. Genetic Resources and Crop Evolution, v.50, p.101‑114, 2003. [ Links ]
ROY, C.B.; NAZEER, M.A.; SAHA, T. Identification of simple sequence repeats in rubber (Hevea brasiliensis). Current Science, v.87, p.807‑811, 2004. [ Links ]
SAHA, T.; ROY, C.B.; NAZEER, M.A. Microsatellite variability and its use in the characterization of cultivated clones of Hevea brasiliensis. Plant Breeding, v.124, p.86‑92, 2005. [ Links ]
SAHA, T.; ROY, C.B.; RAVINDRAN, M.; BINI, K.; NAZEER, M.A. Allelic diversity revealed through SSR polymorphisms at the locus encoding HMG‑CoA reductase in rubber (Hevea brasiliensis). Silvae Genetica, v.56, p.58‑65, 2007. [ Links ]
SOUZA, L.M.; MANTELLO, C.C.; SANTOS, M.O.; GONÇALVES, P. de S.; SOUZA, A.P. Microsatellites from rubber tree (Hevea brasiliensis) for genetic diversity analysis and cross‑amplification in six Hevea wild species. Conservation Genetics Resources, v.1, p.75‑79, 2009. [ Links ]
TESSIER, C.; DAVID, J.; THIS, P.; BOURSIQUOT, J.M.; CHARRIER, A. Optimization of the choice of molecular markers for varietal identification in Vitis vinifera L. Theoretical and Applied Genetics, v.98, p.171‑177, 1999. [ Links ]
THIEL, T.; MICHALEK, W.; VARSHNEY, R.K.; GRANER, A. Exploiting EST databases for the development and characterization of gene‑derived SSR‑markers in barley (Hordeum vulgare L.). Theoretical and Applied Genetics, v.106, p.411‑422, 2003. [ Links ]
VAN TREUREN, R.; KEMP, H.; ERNSTING, G.; JONGEJANS, B.; HOUTMAN, H.; VISSER, L. Microsatellite genotyping of apple (Malus x domestica Borkh.) genetic resources in the Netherlands: application in collection management and variety identification. Genetic Resources and Crop Evolution, v.57, p.853‑865, 2010. [ Links ]
VARSHNEY, R.K.; GRANER, A.; SORRELLS, M.E. Genic microsatellite markers in plants: features and applications. Trends in Biotechnology, v.23, p.48‑55, 2005. [ Links ]
VARSHNEY, R.K.; SIGMUND, R.; BORNER, A.; KORZUN, V.; STEIN, N.; SORRELLS, M.E.; LANGRIDGE, P.; GRANER, A. Interspecific transferability and comparative mapping of barley EST‑SSR markers in wheat, rye and rice. Plant Science, v.168, p.195‑202, 2005. [ Links ]
WEIR, B.S.; COCKERHAM, C.C. Estimating F‑statistics for the analysis of population structure. Evolution, v.38, p.1358‑1370, 1984. [ Links ]
WYCHERLEY, P.R. Introduction of Hevea to the Orient. The Planter, v.44, p.127‑137, 1968. [ Links ]
YU, F.; WANG, B.‑H.; FENG, S.‑P.; WANG, J.‑Y.; LI, W.‑G.; WU, Y.‑T. Development, characterization, and cross‑species/genera transferability of SSR markers for rubber tree (Hevea brasiliensis). Plant Cell Reports, v.30, p.335‑344, 2011. [ Links ]
Received on April 4, 2012 and accepted on July 16, 2012