Morphological and Molecular Characterization of Thamnocalamus Spathiflorus Subsp. Spathiflorus at Population Level

Plant Syst Evol (2009) 282:13–20 DOI 10.1007/s00606-008-0092-1 ORIGINAL ARTICLE Morphological and molecular characterization of Thamnocalamus spathiflorus subsp. spathiflorus at population level Samik Bhattacharya Æ Jayadri Sekhar Ghosh Æ Malay Das Æ Amita Pal Received: 29 January 2008 / Accepted: 12 August 2008 / Published online: 30 April 2009 Ó Springer-Verlag 2009 Abstract Comprehensive morphological characterization and an incidence of gregarious flowering in Thamnocalamus spathiflorus (Trin.) Munro subsp. spathiflorus are described. Twenty-eight key vegetative characters as well as reproductive morphology were studied. They are in gross agreement with prior taxonomic descriptions, yet more elaborate. Statistical analyses of the quantitative vegetative characters revealed significant high variability existing between populations, but not within populations. However, DNA fingerprinting analyses by applying highly polymorphic random primers could not detect any polymorphism either between populations or within populations. Insignificant within-population variability indicates possibility of clonal propagation from the donor(s) possessing similar genetic background and thus reducing genetic variability. Keywords Morphology Á Gregarious flowering Á DNA fingerprinting Á Genetic diversity Introduction Thamnocalamus spathiflorus (Trin.) Munro subsp. spathiflorus belongs to the subtribe Thamnocalaminae (=Arundinariinae) under subfamily Bambusoideae. The group includes several cold-tolerant members (http://www. bambooresearch.com/Specie/species-thamno.html) and is usually found in fragile ecosystems like canal banks, steep terrain and eroding hills. Hence, they play important role in stabilizing top soil by its widespread rhizomatous rooting system. On the basis of the flowering cycle, bamboos have been categorized in three major groups: annual flowering, sporadic or irregular flowering and gregarious flowering (Brandis 1899). In gregarious flowering, the onset of the reproductive phase occurs simultaneously in all the members of a cohort or species across a wide area and subsequently dies together. Records of bamboo flowering are limited in the literature (Siefriz 1950; Filgueiras and Pereira 1988; Ramanayake and Yakandawala 1998; Bhattacharya et al. 2006). In India, the past incidents of mast flowering in T. spathiflorus were recorded during 1818–1821 from North-Western Himalayas (Fig. 1, location A), during 1881–1882 from Jaunsar and Tehri Garhwal (Broun 1886; Fig. 1, location B), in 1942 from Chakrata (Fig. 1, location C) and in 2002 from Uttaranchal (Naithani et al. 2003; Fig. 1, location D). The estimated flowering period in T. spathiflorus is 60 years (Naithani et al. 2003). Here, we report the incidence of a gregarious flowering event in T. spathiflorus subsp. spathiflorus growing over an area of approximately 3.5 km2 in Sikkim, India. This incited us to re-evaluate the reproductive characteristics of the subspecies and compare with previous descriptions. Nevertheless, like most bamboo species, the identification of T. spathiflorus subsp. spathiflorus is usually dependant on vegetative Electronic supplementary material The online version of this article (doi:10.1007/s00606-008-0092-1) contains supplementary material, which is available to authorized users. S. Bhattacharya Á J. S. Ghosh Á A. Pal (&) Plant Molecular and Cellular Genetics Section, Bose Institute, P 1/12 CIT Scheme VIIM, Kolkata 700054, India e-mail: amita@bic.boseinst.ernet.in M. Das Department of Plant Pathology, Physiology and Weed Science, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA 123 14 S. Bhattacharya et al. Fig. 1 Map showing past flowering records in T. spathiflorus at various locations of India and Nepal (reported during A: 1818–1821; B: 1881–1882; C: 1942; D: 2002; E: 2006). An enlarged view of the flowering site (E) depicting nine locations surveyed in the present study; 1, Namchipong; 2, Khempur; 3, Kyongkhosla Falls; 4, 15 Mile; 5, 14 Mile; 6, 10 Mile; 7, 9 Mile; 8, 5 Mile; 9, Gangtok BSI. Gregarious flowering was observed only at locations 3, 4, 5, 6 and 7 characters, which often lead to taxonomic conflicts (Stapleton 1994). Similarly, in a number of instances, we have identified significant variability in vegetative features, which pose difficulty in distinguishing closely allied woody bamboo species (Bhattacharya et al. 2006). Hence, in this study, we attempted to assess the extent of population genetic diversity in T. spathiflorus by simultaneously utilizing a wide array of morphological characters and multilocus DNA fingerprinting profiles. In contrast to widespread interests on bamboo taxonomy and systematics, efforts on population level genetic diversity studies are in its infancy, till date (Das et al. 2008). Our previous study on B. tulda (Bhattacharya et al. 2006) indicated a low population genetic diversity existing in this species. Similar trend was reported in Phyllostachys pubescens from Taiwan (Lai and Hsiao 1997) and Guadua angustifolia from Colombia (Marulanda et al. 2002). It is likely that only few clones of individual species acted as the genetic donor within a particular geographic area resulting in low level among population genetic diversity. On the contrary, relatively higher clonal variation was found in Sasa senanensis from Japan (Suyama et al. 2000) and G. amplexifolia from Colombia (Marulanda et al. 2002). Most possibly, the differential reproductive systems have influenced the intensity of population genetic diversity, since it is expected that the allogamous species are usually more diverse than the autogamous ones. Table 1). Longitude, latitude and altitude of the collection sites were recorded with a MagellanÒ SportrackÒ GPS system, while temperature and rainfall data were collected from the local meteorological department. Vegetative characters Type of the rhizome was determined by digging the soil or from the already exposed soil at hill slopes. Among the selected 28 vegetative characters (qualitative plus quantitative; Das et al. 2007), 14 culm features were recorded at site, while 14 culm sheath characters were scored after detailed investigation. Each quantitative morphological character represents the mean of three randomly chosen culms from five individual clumps per each population. Flowering incidences, inflorescence and floral characters Flowering was recorded at five distant locations in East Sikkim, India, during 2006 (Fig. 1). Inflorescences were collected and key morphological characters of fresh spikelets were recorded. Five individual fertile florets were randomly selected to characterize each population. Genomic DNA isolation Genomic DNA was isolated from young leaves of five randomly selected culms as representatives of each population. Surface-sterilized leaf tissues weighing 0.1 g were sliced into small pieces and homogenized in liquid nitrogen using mortar and pestle. DNA was extracted using 25 ml warm CTAB extraction buffer following Doyle and Doyle (1987). After removal of RNA by RNAase treatment, DNA concentrations were determined by using an UV Materials and methods Data scoring and collection of materials Wild stands of Thamnocalamus spathiflorus (Trin.) Munro subsp. spathiflorus spanning across nine geographical locations in East Sikkim, India, were surveyed (Fig. 1, 123 Characterization of Thamnocalamus spathiflorus subsp. spathiflorus Table 1 Principal geographical features of the collection sites of T. spathiflorus subsp. spathiflorus Collection numbera AM/NMC/06/68 AM/KMP/06/74 JG/KKF/06/045 JG/15 M/06/056 AM/14 M/06/79 JG/10 M/06/062 JG/09 M/06/074 AM/05 M/06/85 JG/GBG/06/088 a b c 15 Place of collectionb Namchipong Khempur Kyongkhosla Fallsc 15 Milec 14 Mile c Latitude Longitude Altitude (m) Maximum temperature (°C) 11.5 9.0 9.0 8.5 8.5 8.0 9.5 10.5 11.5 Minimum temperature (°C) 2.5 2.0 1.5 1.5 1.5 1.0 1.5 2.5 3.5 Average annual rainfall (mm) 3,245.0 3,554.0 3,674.0 3,688.0 3,696.0 3,285.0 3,195.0 3,575.0 3,455.0 Soil type 27°15.1880 N 88°37.1790 E 1,321 27°22.8390 N 88°41.2580 E 2,945 27°22.3010 N 88°42.6990 E 3,164 27°22.0940 N 88°43.0760 E 3,200 27°22.1590 N 88°42.0130 E 3,289 27°22.210 N 88°41.1700 E 3,442 27°22.2290 N 88°41.0350 E 3,122 27°22.2270 N 88°40.6130 E 2,513 27°20.3290 N 88°36.6390 E 1,769 Silty clay loam Silty clay loam Loamy sand Loamy sand Loamy sand Loamy sand Loamy sand Silty clay loam Silty clay loam 10 Milec 9 Milec 5 Mile Gangtok BSI Collection number of representative individual of each population Locations under different districts of West Bengal, India Gregarious flowering observed spectrophotometer (Beckman-Coulter, DU-520) and purity was checked at A260 and A280. An OD of 1 at A260 is equivalent to 50 lg DNA, and pure DNA, free of protein, gave a value between 2.0 and 1.9 of A260:A280. Development of fingerprinting profiles by random primers Forty random decamer Operon primers (Operon Technologies, Germany) were initially screened to select 22 primers that generated distinct and highly reproducible amplified fragments (primer sequences are given in Table 4). The amplifications were carried out in a Perkin Elmer Cetus 2400 thermal cycler with 50 ll reaction mixture following the protocol of Das et al. (2005). The ability of these primers to generate reproducible and polymorphic banding patterns across a wide range of bamboo species had already been demonstrated in our laboratory (Das et al. 2007, 2008). Nevertheless, the ability of these primers to differentiate closely related congeners was further examined by carrying out amplification reactions using genomic DNA of Arundinaria maling Gamble [Synonym, Yushania maling (Gamble) Stapleton]. Among all the temperate, Sino-Himalayan, woody bamboo genera growing up to 3,000 m, A. maling is phylogenetically closest to Thamnocalamus (Guo et al. 2002). The thermal cycler programme was as follows: 4 min at 95°C followed by 35 cycles (45 s at 94°C, 45 s at 37°C and 1 min at 72°C) and finally 10 min at 72°C for elongation. The amplified PCR products were resolved by electrophoresis on 1.5% agarose gel in TAE buffer and recorded with a gel doc 1000 camera using the software programme Molecular Analyst version 1.5 (Bio Rad, USA). The exact molecular weights of the PCR-amplified fragments were calibrated from the standard 100 bp molecular weight marker (Genei, USA) by the same software. Statistical analysis One-way analysis of variance (ANOVA) was performed to estimate morphological variability between and within nine populations of T. spathiflorus subsp. spathiflorus. The quantitative characters considered for ANOVA analysis were height and diameter of culm, internode length, ratio of cavity versus culm diameter, ratio of culm sheath length versus breadth at base and ratio of culm sheath length to blade length. Coefficients of variation (CV%) between and within population for the said parameters were also computed. Results Vegetative characters The perennial, unicaespitose bamboo has short pachymorph rhizome. Culms were 5.8–6.9 m high, usually straight, terete, glaucous, distally pruinose, yellowish green in color (Fig. 2a); 22.0–28.0 mm in diameter; internodes 123 16 Fig. 2 a Vegetative culm and culm sheath of T. spathiflorus subsp. spathiflorus; b photograph of slightly asymmetric blade of the culm sheath with scanty hairs on abaxial surface; c wild stands of T. spathiflorus subsp. spathiflorus with gregarious flowering on stiff mountain slope; d close-up view of a flowering branch S. Bhattacharya et al. straight, 220–280 mm long; ratio of cavity diameter to culm diameter 0.51–0.59; nodes rarely swollen with sheath scars and whitish rings. Branching is typically found in the upper nodes. Culm sheaths were slightly asymmetrical in shape, deciduous, coriacious and almost glabrous, often with scanty hairs on abaxial surface (Fig. 2b), vary in size from base to top of the main culm; ratio of total length to breadth at the base 2.22–2.77; blade straight, ovate; ratio of total length to blade length 2.83–3.19; auricles absent; ligule fimbriate. Flowering incidence and floral morphology During August 2006, gregarious flowering of T. spathiflorus subsp. spathiflorus was noticed at five locations in East Sikkim, India (Fig. 2c, d) at altitudes of 3,122–3,442 m and stretched over an area of 3.5 km2 (Table 1). Inflorescence is a drooping panicle and branching pattern similar to that of a vegetative culm (Fig. 3a). Rachis was smooth, jointed, 100–120 mm with two to three spikelets (Fig. 3b) subtended by shining, chaffy, straw-colored sheathing 123 Characterization of Thamnocalamus spathiflorus subsp. spathiflorus Fig. 3 Reproductive morphology of T. spathiflorus subsp. spathiflorus: a inflorescence with fasciculate arrangement of spikelets; b a single spikelet; c–h individual floral parts: c empty glume; d lemma; e palea; f lodicule; g stamens; h pistil 17 bracts, 50–60 mm long; empty glumes 13–17 mm, six- to nine-nerved, ciliate at the tip and edges (Fig. 3c); spikelets pedicillate, 50–60 mm long, contains four to six fertile florets with distal floret imperfect (Fig. 3b). Lemma (Fig. 3d) was 20–25 mm long, acuminate with ciliate tip, 3–4 mm broad at base, glabrous, 11- to 13-nerved, strawcolored, overlapping with the palea (Fig. 3e); palea subtending a bisexual floret 10–15 mm long, 2–4 mm broad at the base, four-nerved, penicillate, two-keeled, bearing ciliae on the edges (Fig. 3e); lodicules 3, 2.5–3.0 mm long, ovate acute, hyaline, fimbriate, three-nerved (Fig. 3f); stamens 3; anthers 7.0–9.0 mm long, basifixed, purple, emarginated with linear dehiscence (Fig. 3g); ovary oblong, glabrous, 2–2.5 mm long; style short undivided, hairy; stigmas 3, 3–4 mm long, yellow, plumose (Fig. 3h). Seed was not found in any of the populations. Population variability based on vegetative features Analyses of variance was significant between population variations with respect to height and diameter of culm, internode length, ratios of culm cavity to culm diameter, 123 18 S. Bhattacharya et al. Table 2 Analysis of variance estimated within and between nine populations of T. spathiflorus subsp. spathiflorus for height of culm (H, in m), diameter of culm (D, in mm), internode length (I, in mm), ratio of culm cavity to culm diameter (CD), ratio of culm sheath length to breadth at base (LB) and ratio of culm sheath length to length of blade (LL) Source of variation df Sum of squares H Within population Between population Error Total Source of variation 4 8 32 44 df 0.061 2.872 0.779 3.712 D 0.105 0.774 0.211 1.090 Mean squares (MS) H Within population Between population Error Total Computed F Coefficient of variation (between population) a I 3.111 77.111 28.889 109.111 CD 0.0004 0.019 0.004 0.024 LB 0.002 0.893 0.046 0.942 LL 0.002 0.894 0.046 0.942 D 0.026 0.097 0.007 0.025 14.69198a 18.05% I 0.778 9.639 0.903 2.480 10.67692a 21.03% CD 0.0001 0.002 0.0001 0.0005 17.77686a 11.94% LB 0.0005 0.112 0.001 0.021 77.50478a 8.47% LL 0.0003 0.065 0.0004 0.012 144.5635a 3.99% 4 8 32 44 0.015 0.359 0.024 0.084 14.74501a 14.13% P B 0.001 level culm sheath length to breadth at base and culm sheath length to length of blade (Table 2). The most variable character was internode length having highest coefficient of variation between populations (21.03%) with significance tested at P \ 0.001, while within-population variation was reasonably low (CV 1.2–6.3%; Table 3). Few other quantitative characters also revealed significant (at P B 0.001 level) variations (Table 2) such as, height of culm (CV, 14.13%; F = 14.74501), culm diameter (CV, 18.05%; F = 14.69198), ratio of culm cavity to culm diameter (CV, 11.94%; F = 17.77686), culm sheath length to breadth at base (CV, 8.47%; F = 77.50478) and ratio of culm sheath length to length of blade (CV, 3.99%; F = 144.5635). However, within-population variability was significantly low (Table 3). Population genetic diversity based on DNA fingerprinting analyses Informative and reproducible DNA fingerprinting patterns were generated by 22 out of 40 random oligomers employed. Reproducibility of each fingerprinting pattern was confirmed by repetition (Fig. 4). Yet again, the ability of these primers to differentiate closely related congeners was confirmed by simultaneously carrying out amplification reactions with A. maling genomic DNA (Table 4, Figure S1). The approximate size range of the amplified products in T. spathiflorus subsp. spathiflorus was 134 bp to 1.38 kb, while in A. maling it was between 201 bp and 1.23 kb (Table 4). The highest number (7) of bands was obtained with OPOJ-01 and OPA-05, while the lowest (3) Table 3 Coefficient of variation (CV, in %) estimated within nine populations of T. spathiflorus subsp. spathiflorus for height of culm (H, in m), diameter of culm (D, in mm), internode length (I, in mm), ratio of culm cavity to culm diameter (CD), ratio of culm sheath length to breadth at base (LB) and ratio of culm sheath length to length of blade (LL) Populations Namchipong Khempur Kyongkhosla Falls 15 Mile 14 Mile 10 Mile 9 Mile 5 Mile Gangtok BSI a Ha 1.198486 2.817037 1.944039 1.792729 2.043637 1.798384 3.983609 3.245432 1.413277 Da 4.318846 2.828427 3.075956 3.515378 3.074377 3.193359 4.865076 3.521367 2.127738 Ia 2.828427125 6.271632338 3.63588289 3.19335888 2.226514461 1.98450202 1.644167631 4.488879626 1.150326657 CDa 1.543653185 1.818181818 0.76840824 1.543653185 2.413022997 1.437560183 3.064688742 3.45359356 0.649142659 LBa 0.879961083 3.039867379 0.958196814 1.941510436 1.348073007 1.275542344 0.887820616 0.75920487 0.318113295 LLa 0.518167777 1.066142487 0.404110847 0.54092787 0.903843182 0.563203717 0.914057667 0.463762764 0.198943755 P B 0.001 level 123 Characterization of Thamnocalamus spathiflorus subsp. spathiflorus 19 Fig. 4 A representative profile of DNA fingerprinting of T. spathiflorus subsp. spathiflorus. Lane 1 OPOJ-01; Lane 2 OPOJ-03; Lane 3 OPOJ-04; Lane 4 OPOJ-11; Lane 5 OPOJ-18; Lane 6 OPB-01; Lane 7 OPB-02; Lane 8 OPB-03; Lane 9 OPB-04; Lane 10 OPB-05; Lane 11 OPB-06; Lane 12 OPB-07; Lane 13 OPB-08; Lane 14 OPB-18; Lane 15 OPB-19; Lane 16 OPB-20; Lane 17 OPA-02; Lane 18 OPA-04; Lane 19 OPA-05; Lane 20 OPA-10; Lane 21 PW-01; Lane 22 PW-02; M1 100 bp DNA ladder marker, M2 1.5 kb ? 100 bp DNA ladder marker was scored with OPB-02, OPB-19 and PW-02 primers. No polymorphism could be detected among the nine populations studied (Fig. 5). Discussion Exclusively depending on the culm and culm sheath variations, three subspecies (spathiflorus, nepalensis and occidentalis) and two varieties (bhutanensis and crassinodous) had previously been described within T. spathiflorus (Stapleton 1994). The subspecies spathiflorus differs from Table 4 Total number and size range of amplified products generated from the genomic DNA of T. spathiflorus subsp. spathiflorus, using 22 random decamer primers others by having short pachymorph rhizome, short culm node diameter (\3 mm) and slightly asymmetric culm sheath covered with scanty abaxial hairs (Stapleton 1994; Fig. 2a, b). Although the nature of rhizome could not be critically assessed in our study, superficial observation of the clumping pattern confirms its short pachymorph nature. The vegetative and inflorescence/floral morphology of T. spathiflorus subsp. spathiflorus described here are in gross agreements with previous reports (Naithani et al. 2003; Clayton et al. 2006). However, the current study provides further details of individual floral parts of the subspecies spathiflorus, which had not been characterized before. We also illustrate the spatial and temporal patterns of gregarious flowering events in the subspecies reported till date from different climatic zones of India and Nepal (Fig. 1). Geographically isolated clones adopting variable flowering cycles have already been recorded in Dendrocalamus giganteus (Ramanayake and Yakandawala 1998). The flowering cycle reported was 25 years in South India, 40–45 years in North, East, Central India, while 45 years in Bangladesh and 65 years in West Indies. Another population level study enumerated that clones originating from the seeds of same flowering event may not necessarily have the same flowering interval (Isagi et al. 2004). Janzen (1976) was of opinion that, to resolve the issue, more Primers OPOJ-01 OPOJ-03 OPOJ-04 OPOJ-11 OPOJ-18 OPB-01 OPB-02 OPB-03 OPB-04 OPB-05 OPB-06 OPB-07 OPB-08 OPB-18 OPB-19 OPB-20 OPA-02 OPA-04 Nucleotide sequence (50 –30 ) CTGCTGGGAC GTAGACCCGT CCTTGACGCA TCCGCTCTGG AGGGAACGAG GTTTCGCTCC TGATCCCTGG CATCCCCCTG GGACTGGAGT TGCGCCCTTC TGCTCTGCCC GGTGACGCAG GTCCACACGG CCACAGCAGT ACCCCCGAAG GGACCCTTAC TGCCGAGCTG AATCGGGCTG AGGGGTCTTG GTGATCGCAG TGGTCACTGA TGGTCACTGA Total number of amplified fragmentsa 7 (5) 4 (2) 5 (6) 4 (3) 4 (3) 6 (1) 3 (5) 4 (5) 5 (4) 6 (4) 6 (5) 5 (4) 4 (3) 5 (4) 3 (5) 5 (5) 6 (5) 6 (5) 7 (5) 6 (6) 4 (5) 3 (5) Size ranges of amplified fragments (bp)a 451–1,079 (468–1,052) 416–888 (870–1,113) 211–691 (614–1,234) 134–846 (760–1,064) 260–790 (651–906) 387–1,135 (1,015) 148–790 (566–1,113) 380–818 (578–1,113) 310–910 (566–882) 529–1,185 (505–1,003) 225–790 (274–942) 366–980 (614–918) 564–825 (493–845) 167–1,022 (687–930) 451–888 (675–1,052) 345–1,164 (444–1,040) 176–1,382 (554–1,040) 416–1,086 (626–1,003) 366–966 (420–1,210) 402–839 (505–1,064) 451–980 (201–748) 543–959 (505–1,076) a Numbers and size range of amplified products generated from the genomic DNA of A. maling are given within parentheses OPA-05 OPA-10 PW-01 PW-02 123 20 S. Bhattacharya et al. and Dr. P. Singh of Botanical Survey of India, Sibpur, Howrah, for the identification of bamboo material. The help rendered by Dr. Aparajita Mitra during collection of material is sincerely acknowledged. References Bhattacharya S, Das M, Bar R, Pal A (2006) Morphological and molecular characterization of Bambusa tulda with a note on flowering. Ann Bot 98:529–535 Brandis D (1899) Biological notes on Indian bamboos. Ind For 25:1– 25 Broun AF (1886) Seeding of bamboos. Ind For 12:413–414 Clayton WD, Harman KT, Williamson H (2006) GrassBase—the online world grass flora. http://www.kew.org/data/grasses-db. html Das M, Bhattacharya S, Pal A (2005) Generation and characterization of SCARs by cloning and sequencing of RAPD products: a strategy for species specific marker development in bamboo. Ann Bot 95:835–841 Das M, Bhattacharya S, Basak J, Pal A (2007) Phylogenetic relationships among the bamboo species as revealed by morphological characters and polymorphism analyses. Biol Pl 51:667–672 Das M, Bhattacharya S, Singh P, Filgueiras TS, Pal A (2008) Bamboo taxonomy and diversity in the era of molecular markers. Adv Bot Res 47:225–268 Doyle JJ, Doyle JL (1987) A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem Bull 19:11–15 Filgueiras TS, Pereira BAS (1988) On the flowering of Actinocladum verticillatum (Gramineae: Bambusoideae). Biotropica 20:164– 166 Guo ZH, Chen YY, Li DZ (2002) Phylogenetic studies on the Thamnocalamus group and its allies (Gramineae: Bambusoideae) based on ITS sequence data. Mol Phylogenet Evol 22:20– 30 Isagi Y, Shimada K, Kushima H, Tanaka N, Nagao A, Ishikawa T, Ono Dera H, Watanabe S (2004) Clonal structure and flowering traits of a bamboo [Phyllostachys pubescens (Mazel) Ohwi] stand grown from a simultaneous flowering as revealed by AFLP analysis. Mol Ecol 13:2017–2021 Janzen DH (1976) Why bamboos wait so long to flower. Annu Rev Ecol Syst 7:347–391 Lai CC, Hsiao JY (1997) Genetic variation of Phyllostachys pubescens (Bambusoideae, Poaceae) in Taiwan based on DNA polymorphisms. Bot Bull Acad Sin 38:145–152 Marulanda ML, Marquez P, Londono X (2002) AFLP analysis of Guadua angustifolia (Poaceae: Bambusoideae) in Columbia with emphasis on the coffee region. Bamboo Sci Cult 16(1):32–42 Naithani HB, Pal M, Lepcha STS (2003) Gregarious flowering of Thamnocalamus spathiflorus and T. falconeri, bamboos from Uttaranchal, India. Ind For 129:517–526 Ramanayake SMSD, Yakandawala K (1998) Incidence of flowering, death and phenology of development in the giant bamboo Dendrocalamus giganteus Wall. ex Munro. Ann Bot 82:779–785 Siefriz W (1950) Gregarious flowering of Chusquea. Nature 22:635– 636 Stapleton CMA (1994) The bamboos of Nepal and Bhutan Part II. Edinb J Bot 51:275–295 Suyama Y, Obayashi K, Hayashi I (2000) Clonal structure in a dwarf bamboo (Sasa senanensis) population inferred from amplified fragment length polymorphism (AFLP) fingerprints. Mol Ecol 9:901–906 Fig. 5 Monomorphic fingerprinting profiles of one representative for each of the nine populations of T. spathiflorus subsp. spathiflorus, in the order as shown in Table 1, generated using random decamer primers a OPOJ- 04; b OPA- 04; M, 100 bp DNA ladder in a and 1500 bp DNA ladder marker in b attention should be paid to ‘habitat variation’ and ‘genetic versus phenotypic variation.’ The current incidence of gregarious flowering in T. spathiflorus subsp. spathiflorus at five sites out of nine sites surveyed in East Sikkim perhaps suggests the presence of two or more cohorts with variable intermast period. This could be verified only when the remaining four populations flower. The flowering cycle in a particular bamboo species/genetic background is predetermined. Hence, the prediction of mast flowering in a particular region is possible only when we know the clonal history of the donor plants and genetic structure of the populations. Statistical analysis of the quantitative morphological features revealed significant variation between populations (Table 3). However, these variations did not affect the identity of the subspecies, since no significant variation was observed with the subspecies determining culm sheath features. Insignificant within-population variability (Table 4) may suggest clonal propagation from the donor(s) possessing similar genetic background and thus reducing genetic variability. It is also apparent from the fingerprinting studies (representative Fig. 5a, b). Further studies involving a wide range of bamboo species from diverse eco-geographical locations are required to gain a better understanding of the population genetic diversity existing in bamboos. Acknowledgments We acknowledge Council of Scientific and Industrial Research, New Delhi, India, for the financial support [Grant No. 38 (1062)/03/EMR-II], CSIR-NET fellowships to the author (SB) 123