Analysis of Uruguayan weedy rice genetic diversity using AFLP molecular markers María
Teresa Federici* Duncan
Vaughan Norihiko
Tomooka Akita
Kaga Xin
Wang Wang Koji
Doi Marta
Francis Gonzalo Zorrilla Nestor
Saldain * Corresponding author Financial support: Japan Cooperation Agency (JICA)- Individual Training. Keywords: clusters,
dendrograms, polymorphic bands, red rice, population dynamics.
The term "weedy rice" refers to populations of usually annual Oryza species that diminish farmer income both quantitatively through yield reduction, and qualitatively, through lowered commodity value at harvest. Weedy rice occurs in all the major rice growing areas in the tropics, being a particular problem in direct seeded rice agriculture of Latin and North America, the Caribbean, Africa, South and South East Asia (Oka, 1988). Weedy rice occurs in temperate regions where no wild relatives of rice occur in these areas, such as the USA and Europe. Weedy rice occurs where conditions of direct seeding and intensive production prevail (Bres-Patry et al. 2001). Weedy rice may have indica or japonica characteristics and show weedy traits such as phenotypic plasticity, a high seed dispersal ability and seed dormancy. Weedy rice seeds usually have a red pericarp (thus commonly called "red rice"), earlier tillering plus flowering than cultivated rice, and anthocyanin pigmentation of different plant parts, such as, collar, ligule, grain apiculus, stigma and awns (Cho et al. 1995; Suh et al. 1997). Several studies on weedy rice genetic diversity have been reported. Weedy rice strains also appear to be differentiated into indica and japonica types based on morphological and physiological traits, isozymes, restriction fragment length polymorphism (RFLP) and random amplified polymorphic DNA (RAPD) markers (Cho et al. 1995; Suh et al. 1997). DNA fingerprinting using simple sequence repeats (SSR or microsatellites) has also been reported indicating differentiation among weedy rice populations (awnless types; straw hull and black hull awned types; cultivated rice and hybrids between cultivars and each weedy type respectively (Gealy, 2000, personal communication). Recently, amplified fragment length polymorphism (AFLP) methodology has been found to be a useful tool to study relationships within members of various taxa (i.e. Caicedo et al. 1999; Yoon et al. 2000). The AFLP technique is a novel polymerase chain reaction (PCR)- based molecular marker assay (Vos et al. 1995) that has the capacity to detect a higher number of loci and thus a higher rate of polymorphism in a single assay than RFLPs or RAPDs (Powell et al. 1996). AFLP methodology has been employed, because of their effectiveness and reliability for studies on diversity, phylogeny, genomic linkage mapping, and identification of varieties (Becker et al. 1995; Maughan et al. 1996; Xu et al. 2000). In this study genetic diversity among weedy and cultivated rice was assessed using the AFLP technique. The objectives of the study were to use the AFLP technique to fingerprint 26 biotypes of weedy rice and 6 cultivars from Uruguay and analyze genetic diversity in these materials; determine genetic relationships between weedy rice and Uruguayan cultivars and find out if this provides insights into the origin of weedy rice; and relate diversity analysis to strategies to control weedy rice in Uruguayan fields.
Plant materials. Panicles of different red rice biotypes were collected from rice fields across Uruguay from 1996 to 1998 (Table 1, Figure 1). Individual plants of each type were transplanted in the field in 1999. Records of morphological traits used as variety descriptors were taken on a single plant basis. In March 2000, seeds were harvested from these plants. The seeds were grown in the tropical greenhouse of NIAR, Tsukuba, Japan for DNA extraction. DNA analysis. Total genomic DNA was extracted from young leaves of 20 days old plants using the CTAB method (Saghai-Maroof et al. 1984). Briefly, 3 g of fresh leaf tissue were finely ground after immersion in liquid nitrogen (at this step material can be stored -80ºC). 10% CTAB solution containing 0.2% beta-mercaptoethanol, was added to the finely ground leaves and incubated at 65ºC in a water bath with shaking for 20 min. The solution was mixed with 20 ml. of chloroform/ isoamylalcohol (24:1), mixed by inversion at room temperature during 20 min, and then centrifuged at 5000 rpm at RT during another 20 min. The supernatant was then mixed with 2 ml of pre- heated 10% CTAB solution and incubated for 3-5 min at 65ºC. After that, an equal volume of chloroform /isoamylalcohol was added, mixed at room temperature for 20 min and centrifuged at 5000 rpm during 20 min (at room temperature). Then, 1.5 volumes of CTAB precipitation buffer (1% CTAB, 5 mM Tris HCl, 10 mM EDTA pH 8.0) was added and fully mixed. After storing at room temperature 10 min, the mixture was centrifuged at 3000 rpm for 20 min (also at room temperature). After that, samples were incubated with 10 ml of buffer (1M NaCl, 10 mM Tris HCl, 1 mM EDTA) at 65ºC to dissolve pellet completely and treated with 8 to 10 ml of RNAse (810 mg/ml) at 37ºC during 30 min. Then, the enzyme was inactivated by incubating samples at 65ºC during 10 min. After ethanol precipitation, the pellet was washed twice with 70% ethanol during 7 min followed by another 5 min wash with 99.5% ethanol, dissolved in 0.1X TE and stored at -20ºC for use. DNA concentrations were estimated in 1% agarose gels and adjusted to 50 ng/µl. PCRs were performed as described by Vos et al. (1995) with minor modifications. The procedures followed were: 500 ng of DNA from each sample were digested with 5 units of EcoRI (recognition sequence 5`GAATTC-3`) and 5 units of Mse I (5' TTAA 3') in a final volume of 40 µl. This mixture was incubated at 37ºC during 3:30 hours. The reaction buffer was 5X RL-buffer (50 mM Tris. HAc pH 7.5, 50 mM MgAc, 250 mM KAc, 25 mM DTT, 250 ng/µl BSA). To the double digested DNA sample, a total volume of 10 µL of ligation reaction mixture including 5 pmoles of EcoRI adapter and 50 pmoles of Mse I adaptor, 10 mM ATP, 1 unit of T4-DNA ligase, and the same reaction buffer as the above mentioned were added and then incubated overnight at 37ºC. Adapter sequences used are listed (Table 2). Digested/ ligated DNA samples were diluted 10 times with TE. A 20 µl pre-selective PCR mixture contained 5 µl of diluted template mixture, 75 ng of both MseI 00 primer (MseI adaptor sequence) and EcoRI 00 primer (EcoRI adaptor (Table 2), 200 µM dNTPs, 10 mM Tris-HCl (pH 8.3), 1.5 mM MgCl2, 50 mM KCl, and 1 unit Taq polymerase. The cycle profile used was: 94ºC 30 sec, 65ºC 30 sec, 72ºC 60 sec (cycle 1), then annealing temperature is lowered 0.7ºC each cycle during next 12 cycles and 94ºC 30 sec, 56ºC 30 sec and 72ºC 60 sec from cycles 14 to 36. The pre-amplification product was diluted 20 times with TE (10 mM Tris HCl, pH 8.0, 0.1 mM and EDTA- Etylene-diamine-tetra-acetic acid), and then used as a template for selective amplification to generate AFLPs. The selective amplification reaction was conducted in a final volume of 20 µl containing 5 µL of diluted preamplification reaction, 30 ng of each MseI+3 and EcoRI+3 primers (Table 2), 200 µM dNTPs, 10 mM Tris-HCl (pH 8.3), 1,5 mM MgCl2, 50 mM KCl, and 1 unit Taq polymerase. Amplifications were conducted on a 9600 Perkin Elmer DNA thermocycler (Norwalk, CT). The PCR was programmed with the same cycle profile above mentioned.A total of 5 µl of the AFLP selective amplification product were mixed with 2.5 µl of loading buffer (98% formamide, 10 mM EDTA, 0,25% each of bromophenol blue and xylene cyanol FF), then denatured at 95ºC for 2 min and place immediately on ice. Electrophoresis was performed on an EC 160 standard sequencing unit (E-C Apparatus Corporation, USA). About 3 µl of mixture were loaded onto 6% (w/v) polyacrylamide denaturing gel with 1X TBE buffer after a pre-run electrophoresis at 50 W for 1 hour, and then were run at 50 W for about 2 hours. After electrophoresis, the gel was silver- stained using the silver sequence DNA staining kit (Promega cat # Q4132, Madison, WI) following the instructions provided. AFLP bands were visually scored with the aid of digital pictures of the gels and adobe Photoshop computer software. AFLP bands were scored as present (1) and absent (0) and only bands showing unambiguous polymorphism were entered into a data matrix that was used for the following analyses. Genetic similarity was calculated with the simple matching coefficients (Sokal and Michener, 1958) by the NTSYS-pc version 1.70 computer software (Rohlf, 2000). Then, the similarity matrices were subjected to cluster analysis by the unweighted pair-group method with the arithmetic averages (UPGMA) cluster method and Plus analysis (Innan et al. 1999).
Each of the 10 primer combinations was able to generate an informative AFLP fingerprint across all 32 DNA samples including 26 weedy accessions and 6 cultivars (Table 1, Figure 2). Each primer pair produced between 12 and 44 distinguishable AFLPs bands with an average of (25.6) per primer pair (Table 3). About 130 unambiguous polymorphic AFLP bands were scored across the 32 samples, of which 127 bands were polymorphic among weedy material and 71 were polymorphic among cultivars. Associations among the accessions were investigated by UPGMA cluster analysis using NTSYS software (Figure 3). Samples formed three clusters based on genetic distances, two clusters (A and C) consisted of weedy samples only and one cluster (B) included all six cultivars analyzed and five weedy samples. Within cluster B there were two subgroups. One included two cultivars (INIA Tacuari, INIA Cuaro) and three weedy accessions (Tacuarembo-6, Unknown-2 and Rio Branco-5). The other included 4 cultivars (INIA Caraguata, Bluebelle, INIA Zapata, EEA 404) and two weedy accessions (Zapata-1, Unknown-3). Principal component analysis (PCA) based on genetic similarity (Figure 4) shows the relationships among accessions in two dimensional space. This shows that several weedy samples appear intermediate between groups, such as Tacuarembó-4, Vichadero-3 and unknown-3 and 4 (Figure 4). Associations among accessions were also analyzed by the Neighboring Joining Method (NJ) using AFLPDIST software to construct a matrix based on inferred nucleotide diversity, using a modified version of software developed by Innan et al. 1999 (Figure 5, Figure 6). Dendrograms with bootstrap values showing genetic relationships among accessions were constructed using Phylip software (Felsenstein, 1995) (Figure 5). The different analyses gave similar results. The samples in each of the AFLP groups were compared to geographic origin and morphological characteristics (Figure 1, Table 1). The weedy rice samples came from two major geographical regions of Uruguay, the east around Trent y Tres and about 80 km to the north around Tacuarembó department. While most samples in Group A are from the east (9 out of 12 samples) and most samples from the north were group C (4 out of 7) there was no clear trend. Based on morphological characteristics from AFLP group A, 75% (9 samples) have black hull and long or short awns. In group B about 90% of samples have straw hull and short or no awn. Group C consists only of samples with straw hull and no awn.
The AFLP technique using silver staining to visualize bands has been successfully applied to the molecular characterization of Uruguayan weedy and cultivated rice. Two types of dendrograms were constructed to analyze genetic diversity between samples: one based on phenotypic data (NTSYS; Felsenstein, 1995) (Figure 3, Figure 4) and the other using phenotypic data to generate a genetic matrix that was subsequently analysed (AFLPDIST; Innan et al. 1999) (Figure 5, Figure 6). Both dendrograms showed similar results. Three clusters, A, B and C, were differentiated by both sets of analyses. Clusters A and B consist of only weedy rice samples whereas cluster B consists of all the cultivars and some weedy rice samples. No clear trend was detected between geographical origin of samples and genetic diversity. This may reflect the rather small sample number for which location was known and the limited geographic area from which samples were collected. It seems from the results that in northern and eastern rice growing areas of Uruguay all three genetic groups detected by AFLP analysis can be found. There was a clear relationship between AFLP groups and morphological characters. Wild type traits, such as, black hull and long awn was found in most of the samples in group A, whereas traits found in cultivated rice were found in the weedy rice of group B and C. It appears that weedy rice in Uruguay has two strategies for survival. The one is to adapt to the natural environment and the other is to mimic the cultigen. A hypothesis to explain the evolution of weedy rice in Uruguay is presented in Figure 7a. This hypothesis suggests that at different times there has been selection for shattering and subsequently disruptive selection for either wild traits (Group A) or domestication traits (Group C). Recently evolved weedy types that still are genetically similar to present day cultivars are represented by weedy samples in group B. At the same time hybridization-differentiation cycles between weedy rice and cultivars results in intermediate types that were found in this study such as Tacuarembó-4, Vichadero-3, unknown 3 and 4 (Figure 7b) (Harlan, 1975). The results of this study suggest that the control of weedy rice that have evolved from rice should consider the following:
This research was conducted while the author was a JICA fellow of the Japanese government in Crop Evolutionary Dynamics Laboratory, Genetic Resources 2, National Institute of Agrobiological Resources, Tsukuba, Japan. I express my thanks to Dr. Duncan A. Vaughan, head of the Crop Evolutionary Dynmics Laboratory for receiving me and kindly supervising my work; to staff of the Crop Evolutionary Dynamics Laboratory particularly Dr. Xin Wang Wang for technical advice and Dr. Koji Doi and Dr. Akito Kaga for their assistance in computer processing. I express my thanks to Dra. Marta Francis (Unit of Biotechnology, INIA Las Brujas, Uruguay) for supporting my fellowship; and to Dr. Gonzalo Zorrilla and Ing. Nestor Saldain (INIA Treinta y Tres, Uruguay) for providing access and information on morphological characteristics of the samples. I am also grateful to Dr. Jordi García (IRTA, Barcelona) for corrections and suggestions about my work.
Bres- Patry, C.; Lorieux, M.; Clement, G.; Bangratz, M. and Gesquiere, A. (2001). Heredity and genetic mapping of domestication- related traits in a temperate japonica weedy rice. Theoretical Applied Genetics 102:118-126. Caicedo, A. L.; Caitan, E.; Duque, M.C., Toro Chica, O.; Debouck, D.G. and Tohme, J. (1999). AFLP fingerprinting of Phaseolus lunatus L. and related wild species form South America. Crop Science 39:1497-1509. Cho, Y.G., Blair, M.W.; Panaud, O. and McCouch, S.R. (1995). Cloning and mapping of variety specific rice genomic DNA sequences amplified length fragment polymorphisms (AFLP) from silver-stained polyacrylamide gels. Genome 39:373-378. Felsenstein, J. (1995). Phylogenies from restriction sites, a maximum likelihood approach. Evolution 46:557-574. Harlan, J.R. (1975). Crops and man. American Society of Agronomy, Madison, Wisconsin. Innan, H.; Terauchi, R.; Kahl, G. and Tajima, F. (1999). A method for estimating nucleotide diversity from AFLP data. Genetics 151:1157- 1164. Maughan, P.J.; Saghai Maroof, M.A.; Buss, G.R. and Huestis, G.M. (1996). Amplified fragment length polymorphism (AFLP) in soybean: Species diversity, inheritance and near isogenetic line analysis. Theoretical Applied Genetics 94:255-263. Oka, H.I. (1988). Origin of cultivated rice. Japan Scientific Societies Press, Tokyo. 254 pp. Powell, W.; Morgante, M.; Andre, C.; Hanafey, M.; Vogel, J.; Tingey, S. and Rafalski, A. (1996). The comparison of RFLP, RAPD, AFLP and SSR (microsatellite) markers for germplasm analysis. Molecular Breeding 2:225-238. Rohlf, F.J. (2000). NTSYS-pc: Numerical taxonomy and multivariate analysis system. Version 2.1 Exeter Publications, New York, USA. Saghai-Maroof, K.; Soliman, M.; Jorgensen, R.A. and Allard, R.W. (1984). Ribosomal DNA Spacer-Length Polymorphisms in Barley: Mendelian Inheritance, Chromosomal Location, and Population Dynamics. Proceedings of the National Academy of Sciences of the United States of America 81:8014-8018. Sokal, R.R. and Michener, C.D. (1958). A statistical method for evaluating systematic relationships. Univ. Kan. Sci. Bull. 38:1409-1438. Suh, H.S.; Sato, Y.I. and Morishima, H. (1997). Genetic characterization of weedy rice (Oryza sativa L.) based on morpho-physiology, isozymes and RAPD markers. Theoretical Applied Genetics 94:316-321. Vaughan, D.A.; Watanabe, H.; HilleRisLambers, D.; Abdullah Zain, M. and Tomooka, N. (1999). Weedy rice complexes in direct seeding rice cultures. Proc. Int. Symp. "World Food Security", Kyoto, Japan. pp. 277-280. Vos, P.; Hogers, R.; Bleeker, M.; Rijans, M.; Van de Lee, T.; Hornes, M.; Frijters, A.; Pot, J.; Peleman, J.; Kuiper M. and Zabeau, M. (1995). AFLP: A new technique for DNA fingerprinting. Nucleic Acids Research 23:4407-4414. Xu, R.Q; Tomooka, N. and Vaughan, D.A. (2000). AFLP Markers for characterizing the azuki bean complex. Crop Science 40:808-815. Yoon, M.S.; Doi, K.; Kaga, A.; Tomooka, N. and Vaughan, D.A. (2000). Analysis of the genetic diversity in the Vigna minima complex and related species in East Asia. Journal of Plant Research 113:375-386. |
|
Note: Electronic Journal of Biotechnology is not responsible if on-line references cited on manuscripts are not available any more after the date of publication. |
||