BioMed Research International

BioMed Research International / 2013 / Article

Research Article | Open Access

Volume 2013 |Article ID 289285 | 11 pages | https://doi.org/10.1155/2013/289285

Grappling the High Altitude for Safe Edible Bamboo Shoots with Rich Nutritional Attributes and Escaping Cyanogenic Toxicity

Academic Editor: Kuo-Chen Chou
Received18 Jun 2013
Accepted03 Oct 2013
Published19 Nov 2013

Abstract

Consumption of bamboo species with high level of total cyanogenic content (TCC) in Asia by many ethnic groups is significantly associated with food poisoning and occasionally Konzo (a neurological disorder). Adequate characterization of edible bamboo species with low level of TCC and high nutritious attributes is required for consumer’s safety as well as for the conservation of the gene pool. Here, we employed morphological descriptors, atomic absorption spectrophotometer, RAPD, and trnL-F intergenic spacer to characterize 15 indigenous edible bamboo species of north-east India. The study indicates that morphologically and genetically evolved edible bamboo species having large and robust bamboo-shoot texture and growing at low altitude contain high level of TCC, low antioxidant properties, and low levels of beneficial macronutrients and micronutrients. Importantly, Dendrocalamus species are shown to be rich in TCC irrespective of the growing altitude while Bambusa species are found to have moderate level of TCC. The findings clearly demonstrated that Chimonobambusa callosa growing at high altitude represents safe edible bamboo species with nutritious attributes.

1. Introduction

Bamboo shoots are popular traditional food delicacies which are consumed as fresh, fermented, or canned in many South-East Asian countries. Cyanogenic glycosides are inherently produced in cyanogenic plants as defence arsenals and abundantly produced in bamboo shoots in the form of taxiphyllin [1, 2]. High intake of cyanogenic glycosides is life threatening and significantly associated with neurological disorder called Konzo [3]. Under optimal conditions, lactic acid from fermentation of bamboo shoot reduces total cyanide content (TCC) [4]. Nonetheless, better results can be achieved if the initial starting material is poor in TCC. Therefore, ingestion of fresh or inappropriately fermented bamboo shoots can lead to cyanide poisoning.

Irrespective of toxicity, some young edible bamboo shoots (≤30 days) possess enormous nutritious potentials such as high fibre content with antioxidant and antitumor properties [5, 6]. Consequently, the demand for bamboo shoots is high and farmers are often challenged to match supply in both quality and quantity. The magnitude of this increasing demand is exacerbated by nondomestication of edible bamboo species. As a result, nonedible bamboo (or poisonous) shoots are made available in the market which are harmful to consumers. In order to effectively outwit this anthropogenic pressure, a holistic approach is required to identify the species with rich nutritional attributes for safe human consumption and domestication.

Biochemical studies on antioxidant activity, cyanogenic glycoside content, and nutrient content have been reported for a few species [5, 79], without correlation to the altitudes of sample collection, equally lacking morphological and molecular characterization. Based on previous work, it is difficult to sort out edible bamboo species, requiring low-processing capital inputs to eliminate the poisonous components. Here, we exploit the geographic positioning of edible bamboo species in a dynamic biota and their morphological descriptors, RAPD and trnL-F intergenic spacer, to study the interrelatedness of coevolved species of genera Bambusa, Dendrocalamus, Chimonobambusa, Schizostachyum, and Melocanna. We also exploit the geographic positioning of same bamboo species to shed light on how altitude influences the valuable nutritional attributes of bamboo shoots.

2. Material and Method

2.1. Plant Material

Fifteen species of edible bamboo-shoots belonging to the genera Bambusa, Dendrocalamus, Chimonobambusa, Schizostachyum, and Melocanna were collected from different altitudes of Manipur, India (23°47′–25°41′ NL; 92°58′–94°47′ EL), during July-August of 2009-2010. This region often receives an average rainfall of 1320 ± 3 mm and temperature of 23 ± 3°C during the months of July-August. This group of samples was used for morphological characterization and authenticated by the Botanical Survey of India (BSI), Kolkata. The voucher specimens were deposited at the Central National Herbarium in BSI for future references (Table 1). For biochemical analysis, 30-day-old bamboo-shoots (from the day of emergence) from the field were collected in 2009, 2010, and 2011 following routine sampling. Subsequent to harvest, the bamboo-shoots were directly frozen in liquid nitrogen in a thermocol box and transported to the laboratory where they were stored at −80°C for downstream analysis. This group of samples was collected at sunset and the soil pH for the sites of collection was determined as previously described [10].


GenBank accession/bamboo species12345Voucher no.

KC013282/Chimonobambusa callosa Leimaram 24°44′93°43′1843 m (hill)IBSD/WS/019
KC013285/Bambusa cacharensis Arapti 24°44′93°51′807 m (hill)IBSD/WS/020
JX564900/Bambusa manipureana Khong Khang 24°21′94°11′240 m (valley)IBSD/WS/008
JX564901/Bambusa nutans Arapti 24°44′93°56′226 m (valley)IBSD/WS/023
JX507132/Bambusa tulda Arapti 24°44′93°56′226 m (valley)IBSD/WS/022
JX507131/Bambusa oliveriana Arapti 24°44′93°56′728 m (hill)IBSD/WS/010
JX564902/Dendrocalamus giganteus Arapti 24°44′93°56′803 m (hill)IBSD/WS/001
JX564903/Dendrocalamus hamiltonii Kwatha 24°19′94°16′358 m (valley)IBSD/WS/004
JX564904/Dendrocalamus hookeri Arapti 24°44′93°56′770 m (hill)IBSD/WS/005
JX564905/Dendrocalamus manipureanus Arapti 24°21′93°57′769 m (hill)IBSD/WS/002
JX507133/Melocanna baccifera Lokchao 24°20′94°14′497 m (hill)IBSD/WS/018
JX507134/Schizostachyum dullooa Arapti 24°42′93°57′728 m (hill)IBSD/WS/003
JX564906/Bambusa sp.Arapti 24°44′93°56′770 m (hill)IBSD/WS/024
JX564907/Bambusa sp.Andro 24°47′93°03′803 m (hill)IBSD/WS/007
KC013288/Bambusa tuldoides Andro 24°47′93°03′777 m (hill)IBSD/WS/006

Collection site, 2mean pH of the soil at the collection site in the months of July-August of 2009 to 2011, 3latitude (N), 4longitude (E), and 5altitude in meters showing either valley or hill where samples were collected.
2.2. Morphological Analysis

To evaluate the morphological interrelatedness among the edible bamboo species, we used 35 morphological descriptors based on 11 culm types, 13 culm-sheath types, and 11 leaf types (Table S1 available online at http://dx.doi.org/10.1155/2013/289285). The morphological data were analysed using NTSYS-PC version 2.2 [11]. Simple coefficient matching was performed using SIMQUAL option for generating dataset similarity matrix [12]. The best dendrogram was computed with Unweighted Pair-Grouped Method Arithmetic Averages (UPGMA) [13]. All the samples used in this study were scored in triplicate.

2.3. DNA Extraction

To characterize at molecular level, genomic DNA was isolated as described by Aras et al. [14]. The quantity and quality of DNA were checked on BioSpec NanoDrop spectrophotometer (Thermo Scientific, USA) and on 0.8% w/v agarose gel electrophoresis, respectively.

2.4. Random Amplified Polymorphic DNA (RAPD) Analysis

PCR was performed using standard RAPD PCR kit (Invitrogen, USA) in a volume of 50 μL in C1000 Touch Thermal Cycler (BIO-RAD, USA). The run was programmed as follows: initial denaturation at 95°C for 5 min, followed by 35 cycles of amplification (95°C for 1 min, 37°C for 1 min, and 72°C for 2 min) and a final extension at 72°C for 7 min. Amplicons were profiled on a 1.8% agarose gel electrophoresis and revealed with ethidium bromide in Gel Doc-it2 Imager (UVP Co., Ltd.). The primers used for this analysis are represented (Table S2). Polymorphism was scored as 1/0 (presence or absence) to generate a primary binary matrix. The primary binary matrix was used for producing similarity data using Jaccard’s similarity coefficient [15] and computed using UPGMA. This analysis was performed in NTSYS-PC software [11].

2.5. Analysis of trnL-F Intergenic Spacer

In order to authenticate the edible bamboo species, the trnL-F region was amplified using the primers set (forward: 5′-ggttcaagtccctctatccc-3′; reverse: 5′-atttgaactggtgacacgag-3′) as described in Taberlet et al. [16]. PCR products of about 350–400 bp were purified and sequenced in ABI370X1 Cycler Sequencer (ABI, USA) using the same set of primers. Sequences were automatically trimmed and assembled in DNAbaser 3.5.3 software (http://www.dnabaser.com/). Following annotation, sequences were assigned to molecular species based on 98–100% sequence similarity threshold in the GenBank (http://www.ncbi.nlm.nih.gov/) and in accordance with morphological descriptors. The sequences are available in GenBank as accessions JX564900 (Bambusa manipureana), JX564901 (Bambusa nutans), JX564902 (Dendrocalamus giganteus), JX564903 (Dendrocalamus hamiltonii), JX564904 (Dendrocalamus hookeri), JX564905 (Dendrocalamus manipureanus), JX507131 (Bambusa oliveriana), JX507132 (Bambusa tulda), JX507133 (Melocanna baccifera), JX507134 (Schizostachyum dullooa), KC013282 (Chimonobambusa callosa), KC013285 (Bambusa cacharensis), JX564906 (Bambusa sp.), JX564907 (Bambusa sp.), and KC013288 (Bambusa tuldoides), respectively.

Sequences were aligned using Clustal omega program [17]. The program BioEdit [18] was used to assess the nature of variability and the entropy of the alignment among the species. The program TOPALi v.2.5 [19] was used to select the best substitution model and the best method for phylogenetic tree reconstruction based on Akaike Information Criterion, corrected AICc1 and AICc2, and Bayesian Information Criterion (BIC). Phylogenetic analysis was performed using the maximum likelihood (ML) method. Hasegawa-Kishino-Yano nucleotide substitution model was used [20]. The initial trees for the heuristic search were obtained automatically as follows. When the number of common sites was <100 or less than one-fourth of the total number of sites, the maximum parsimony method was used; otherwise BIONJ method with Monte Carlo localization (MCL) distance matrix was used. A discrete Gamma (+G) distribution was used to model evolutionary rate differences among sites (5 categories (+G, parameter = 200.00)). The analysis was performed in MEGA5 [21].

2.6. Determination of Total Cyanide Content

To determine the poisonous potential of each species, the total cyanide content (TCC) was evaluated using the picrate method as described in Bradbury et al. [22] with some modifications. Firstly, bamboo shoot sheaths were removed and the innermost edible portion was measured using a ruler and a slide calliper. The full length was divided into three equal parts, that is, the tip, the middle, and the base. Woody bamboos generally grow rapidly and their shoots are often eaten young [46]. Moreover, based on preliminary findings using 10-, 20-, and 30-day-old bamboo shoots (data not shown), revealed that TCC and other nutritional parameters studied insignificantly varied with only 30-day-old sample. Cogently, only 30-day-old bamboo shoots were used for the biochemical analysis. The standard curve for determination of HCN was established using NaCN solution as follows: 5 mL of alkaline picrate solution (1.4 g of picric acid in 2.5% Na2CO3) and 5 mL of NaCN solution (181 mg of NaCN in 1 L sterile milli-Q water) were pooled together to obtain 100 μg HCN/mL and heated for 5 min in boiling water. Volumes of 0.1, 0.2, 0.4, 0.6, and 0.8 and 1 mL of the resultant NaCN alkaline picrate solution were adjusted to 5 mL with sterile milli-Q water to obtain 5, 10, 20, 30, 40, and 50 μg HCN, respectively.

2.7. Antioxidant Activity Estimation

100 g of sliced bamboo shoot was boiled in 300 mL of double distilled water for 2 h at 100°C. The crude extract was filtered through Whatman no. 42 filter paper and concentrated in a rotary evaporator at 100°C. The solid residues were stored at 4°C till used. The scavenging effect of 2,2-diphenyl-1-picrylhydrazyl (DPPH) free radical was assayed as previously described in Mensor et al. [23]. L-Ascorbic acid (Sigma, USA) was used as reference antioxidant control. One mL of 0.3 mM DPPH ethanolic solution was added to each sample at different concentrations of 20 μg, 50 μg, 100 μg, 200 μg, and 400 μg. The mixture was vortexed for 1 min and then left to stand at room temperature in the dark. After 30 min, absorbance was read at 517 nm in UV1700 spectrophotometer (Shimadzu, USA). The scavenging activity of DPPH free radical was calculated using the following equation: is the absorbance of the control reaction (containing all reagents except for the test compound) and is the absorbance of the test compound. The inhibition concentration (IC50) is defined as the amount of extract required to reduce free scavenging activity by 50%. The IC50 values were obtained from the inhibition curve by extrapolation.

2.8. Estimation of Macro- and Micronutrients

Total nitrogen contents were determined through digestion and distillation of dry bamboo shoots in Kel-Plus digestion system (Pelican, India) according to AOAC [24] protocol. Crude protein was calculated as Kjeldahl N × 6.25 based on the assumption that nitrogen (N) constitutes 16.00% of a protein. The element contents such as potassium (K), sodium (Na), calcium (Ca), magnesium (Mg), copper (Cu), iron (Fe), and zinc (Zn) were estimated using atomic absorption spectrophotometer (AAS) (Perkin Elmer, USA). 500 mg of dried weight (d.w) bamboo shoot powder was digested in 10 : 4 : 1 (HNO3 : HClO4 : H2SO4). The digested sample was appropriately diluted with sterile deionized water and filtered with Whatman no. 42 filter paper. All the elements were analyzed with appropriate multielement hollow cathode lamps (Lumina Lamp, Perkin Elmer) against a standard reference solution for AAS (Accutrace Reference Standard, USA). Phosphorus was estimated by absorbance measurement at 420 nm of the vanadomolybdophosphoric heteropoly complex formed in the digestate [25]. Cellulose (Ce) content was estimated at 620 nm using cold anthrone reagent method [26].

2.9. Statistical Analysis

One-way analysis of variance (ANOVA) was implemented to compare the means of different treatments. The differences between individual means were tested using the least significant difference (LSD). Computation was performed in SPSS software (version 22.0, SPSS Inc., Chicago, USA). The relationship between the different biochemical attributes and 13 edible bamboo species was analysed using principal component analysis (PCA). PCA on standardized data was performed in NTSYS-PC version 2.2 [11]. Principal components with eigenvalues were selected and correlation values were considered as relevant for the PCA. Three genera underrepresented in the study set of 15 were excluded from principal component analysis.

3. Results

3.1. Morphological Analysis

Of the 15 identified edible bamboo species studied, only 13 were morphologically identified at the species level and deposited in BSI, Kolkata (Table 1). The morphological characteristics of the young (<30 days) bamboo shoots are depicted in Figure 1. In the study set, bamboo shoots of JX564902/D. giganteus, JX564903/D. hamiltonii, JX564904/D. hookeri, and JX564905/D. manipureanus were generally deep green, broad based with robust texture (Figures 1(g), 1(h), 1(j), and 1(k)). An overall strong cophenetic correlation coefficient of 0.70 based on morphological characteristics (Table S1) was obtained, indicating a faithfully constructed dendrogram (Figure 2). The bamboo species clustered into two main clades (I and II) with JX507134/S. dullooa evolving in a polyphyletic pattern. The clustering pattern was significantly affected by the morphological characteristics such as colour, shape, and presence of hairs in culm sheaths covering the shoots.

3.2. RAPD Analysis

High level of polymorphism was observed among the 15 species based on the 9 primer sets (Table S2). A strong co-phenetic coefficient of 0.77 based on UPGMA analysis was obtained. A dendrogram (Figure 3) based on band differences revealed three clades (I, II, and III) with C. callosa forming an out group. In this analysis, RAPD evidence of JX564903/D. hamiltonii followed a polyphyletic evolutionary pattern (clade III). The gel profile (Figure S1) provides evidence of JX507131/B. oliveriana (lane 7) and JX564907/Bambusa sp. (lane 15) with higher numbers of dominant characters than those of other species as reflected in the banding pattern.

3.3. DNA Sequence Analysis

Based on DNA sequences, the estimated model parameters were base frequencies (A = 25%, T/U = 25%, C = 25%, and G = 25%) and substitution model , , , , , and . The estimated transition-transversion bias () ratio was at 0.76. The overall mean Tajima-Nei [27] evolutionary distance among the species was 0.51. In the sequence set, the entropy of the alignment (Figure S2) showed 164 patterns (out of a total of 491 sites) and 276 sites were without polymorphism (56.21%). A maximum likelihood tree with the highest log likelihood ) supported by 1000 bootstrap test of replicates showing two main clades (I and II) was generated (Figure 4). It was observed that Dendrocalamus spp. formed a close complex relationship with Bambusa spp. The tree without branch swapping evidence C. callosa has evolved differently from the rest of the edible bamboo species, thus, forming an out group.

3.4. Total Cyanide Content

The level of TCC in all the species varied from 300 to 2604 ppm (for the tip portion), 210 to 2243 ppm (for the middle portion), and 199 to 920 ppm (for the basal portion). These significant differences (at ) in toxicity level suggest that bamboo shoot tips are generally toxic. On the contrary, a low level TCC was observed in C. callosa collected from high altitude in all the studied segments, that is, the tip, the middle, and the base portions of the bamboo shoots (Table 2). Often, all bamboo shoots collected from low altitude were rich in TCC, except for M. baccifera and B. manipureana, with respect to C. callosa collected from high altitude. Overall, all Dendrocalamus species, growing at either high altitude (>700 m) or low altitude (<400 m), were rich in TCC in comparison with other genera (Tables 1 and 2).


SamplesTotal cyanide content (ppm)*IC50 (mg L−1)
TipMiddleBase

KC013282/Chimonobambusa callosa 300a210a199a0.09ab
KC013285/Bambusa cacharensis 1533d1221f735efgh0.58h
JX564900/Bambusa manipureana 1007b515b761fgh0.46fg
JX564901/Bambusa nutans 1001b624c267ab1.83i
JX507132/Bambusa tulda 1579d1406hi779gh0.30ef
JX507131/Bambusa oliveriana 1280c1079e543c0.57gh
JX564902/Dendrocalamus giganteus 2604g2243k920i0.60h
JX564903/Dendrocalamus hamiltonii 1897e766d654cdef0.14bc
JX564904/Dendrocalamus hookeri 1595d1322fgh360b0.45fg
JX564905/Dendrocalamus manipureanus 1838e1270fg600cd0.32de
JX507133/Melocanna baccifera 1548d484b216a0.41ef
JX507134/Schizostachyum dullooa 1521d1279fg727efgh0.61h
JX564906/Bambusa sp.1582d1354gh669defg0.24cd
JX564907/Bambusa sp.2063f1557j635cde0.25cd
KC013288/Bambusa tuldoides 2528g1511ij825hi0.45fg
LSD ( )85.78106.53109.310.11

ppm = mg HCN equivalents/kg bamboo shoots and each value is the mean of three replicates for 2009, 2010, and 2011. The same letter(s) associated with mean values within a column is (are) not significantly different at P ≤ 0.05 based on LSD.
3.5. Antioxidant Activity

Bamboo shoot extract of C. callosa showed the highest significant antioxidant activity (of 53.46%, ) at 400 μg/mL of DPPH (Figure 5), whereas the lowest antioxidant activity (of 2.90%, ) was obtained with extract of B. nutans. Akin to this pattern, the half-inhibition concentration (IC50) obtained by linear regression analysis showed a significant variation from 0.09 mg/L for C. callosa to 1.83 mg/mL for B. nutans (Table 2). Among the three Dendrocalamus species obtained at altitude (>700 m), D. giganteus possessed the least scavenging activity with an IC50 value of 0.60. L-Ascorbic acid investigated under the same conditions had an IC50 value of 0.003 mg/mL; implying that bamboo shoots have a moderate antioxidant activity (Table 2). Based on IC50, D. giganteus is 30-fold poorer in antioxidant activity than an equivalent weight of L-ascorbic acid. When put together, our results suggest that Dendrocalamus species collected at high or low altitude possess low antioxidant attributes (Tables 1 and 2).

3.6. Estimation of Macro- and Micronutrients

Atomic absorption spectrometry (AAS) data revealed that bamboo shoots are generally rich in nitrogen, phosphorous and contain moderate amount of calcium. Remarkably, the highest nitrogen (1153 mg/100 g d.w) and phosphorous (1154 mg/100 g d.w) content was found in C. callosa. By contrast, the lowest nitrogen (673 mg/100 g d.w) and phosphorous content (70 mg/100 g d.w) was observed in D. giganteus. Among the important macronutrients studied (N, P, and K), all the bamboo shoots were rich in potassium and poor in sodium (Table 3). Of all the microelements assessed, iron was significantly present in B. tulda at the rate of 25.80 mg/100 g dry weight (d.w) and was lowest in D. hookeri at the rate of 6.08 mg/100 g d.w. We observed a high level of zinc in B. tulda at the rate of 15 mg/100 g d.w and its lowest level in D. hookeri at the rate of 3.03 mg/100 g d.w. We found a significant amount of copper only in C. callosa at the rate of 6.12 mg/100 g d.w (Table 3).


SamplesNPCeKNaMgCaFeCuZn

KC013282/C. callosa 1153m154g4070.75ab3377e39bc110f63.43gh9.25abcd6.12g7.2de
KC013285/B. cacharensis 815g87ab3157.82a2093abc45cde74ab29.60bc13.15def3.58de5.63bcd
JX564900/B. manipureana 957k94bc5936.51abc2343bcd53e175i32.4bcd16.5f4.94f12.80i
JX564901/B. nutans 825gh104cd5372.79abc2090abc32b175i53.5fg11.6cde4.16e6.50cde
JX507132/B. tulda 721c116de3701.59ab1453abc51de92d42de25.8g3.19cd15.0j
JX507131/B. oliveriana 763e85ab9972.34cd2293abc43cd123g60.95fg10.1bcd3.58de4.88bc
JX564902/D. giganteus 673a70a7048.98abcd3093de22a69a37.9cd7.83abc2.47bc10.7gh
JX564903/D. hamiltonii 747d118de4305.22ab1910abc50de121g24.4ab16.82f3.0cd7.87ef
JX564904/D. hookeri 692b75a5697.05abc2400cd45cde84cd73.3hi6.08a2.64bc3.03a
JX564905/D. manipureanus 894i83ab6784.58abc3533e42cd85cd14.98a6.68ab2.55bc4.80b
JX507133/M. baccifera 1102l135f10126.98cd1757abc39bc149h51.88ef13.37ef1.96ab11.93hi
JX507134/S. dullooa 837h121ef7098.87abcd1310a42cd180i56.95fg14.42ef1.55a9.18fg
JX564906/Bambusa sp.785f73a8306.12bcd1373ab32b92d40.1cd11.5cde2.54bc11.33g
JX564907/Bambusa sp.940j129ef11044.9d1523abc53e112f50.8ef13def7.85h11.03j
KC013288/B. tuldoides 900i122ef4988.66ab2060abc31b101e78.15j11.2cde2.7c11.0h
LSD ( )13.3214.922858.21873.068.518.2810.173.480.641.55

The mean values are for replicates obtained in 2009, 2010, and 2011. The same letter(s) associated with mean values within column is (are) not significantly different at P ≤ 0.05 using LSD.
3.7. Principal Component Analysis (PCA)

Using graphical approaches to study biological problems can provide an intuitive picture or useful insights for analysing complicated relationship in large data set [28], as demonstrated by many previous studies on a series of important biological topics, such as enzyme-catalysed reactions [2931], inhibition of HIV-1 reverse transcriptase [32], drug metabolism systems [33], and using Wenxiang diagram [34] to study protein-protein interactions [35, 36]. The interrelatedness of all the biochemical traits studied and their relationship with the bamboo species based on PCA generated four principal components: PC-1, PC-2, PC-3, and PC-4. The four principal components with eigenvalues accounted for 72.25% of the total variation in nutritional attributes (Table S3). A positive correlation coefficient was observed in PC-1 among phosphorous, iron, sodium, copper, magnesium, zinc, nitrogen, and antioxidant activity, accountable for 32.19% () variation in nutritional content. Interestingly, we observed that PC-2 and PC-3 were highly associated with nitrogen, cellulose, magnesium, copper, and potassium (16.91% variation, ) and total cyanide content, calcium, cellulose, zinc, and copper (14.50% variation, ), respectively. In PC-2 and PC-3, no significant correlation was observed among nutritional attributes. PC-4 was associated with cellulose, antioxidant activity and accounted for 11.65% () net variation in nutritional attributes. The interrelatedness between the nutritional attributes is represented in a polygonal biplot (Figure 6) as previously described [37].

The polygonal biplot (Figure 6) is divided into four sectors by values 2, 3, 4, 5, 6, 7, 8, and 10 with the vertex representing the species. In this representation, the species are the best or the poorest contributing for some or all of the biochemical traits [38], depending on the length of the vector lines from the origin. The biplot indicates that B. manipureana (from 240 m altitude) and B. nutans (from 226 m altitude) had the highest magnesium (Mg), copper (Cu), and nitrogen (N) denoted by 2 and 3, respectively. Similarly, B. tulda (from 226 m altitude) and D. hamiltonii (from 358 m altitude) had the highest value for antioxidant activity (AA), iron (Fe), magnesium (Mg), and zinc (Zn) denoted by 7 and 4, respectively. On the other hand, JX564906/Bambusa sp. (from 770 m altitude) and D. giganteus (from 803 m altitude) had the highest total cyanide content (TCC) denoted by 6 and 10, respectively. The numbers 5 and 8 representing B. oliveriana (from 728 m altitude) and D. hookeri (from 770 m altitude) were found to be tightly linked with high level of calcium (Ca), cellulose (Ce), and potassium (K), respectively. In relative terms, PCA polygonal biplot intuitively evidences that D. giganteus (from 803 m altitude) is rich in TCC. C. callosa, M. baccifera, and S. dullooa were excluded from PCA biplot analysis for reason of underrepresentation of the genera.

4. Discussion

For any meaningful differentiation of edible bamboo species from poisonous ones, the delineating parameters must show a low level of conflicting signals. Using morphological descriptors as discerning tools to categorize the edible bamboo species, the clustering pattern was often significantly affected by dominant morphological characteristics such as colour, shape, and presence of hairs in the culm sheath covering shoots (Figure 1). Considering that most bamboo species and their shoots are green in colour, this produces conflicting signal and renders thought-provoking judgements. This is challenging because the homogeneity and conflicting heterogeneity among species are time consuming and require a high level of expertise. However, the benefit of morphological descriptors is that it permits the initial discernment of poisonous bamboo shoot (which is usually large and has robust-deep green colour) from less toxic bamboo shoots.

Although a less efficient technique, the benefit of RAPD analysis to delineate bamboo species is its rapidity and cost effectiveness. Nevertheless, without any public database for RAPD data, the technique was not effective to barcode our edible bamboo species. The dendrogram generated by RAPD and morphological descriptors had similarity: a cophenetic coefficient is greater than 0.5 but differs by the interspecies clustering pattern. Nonetheless, RAPD discriminated JX564903/D. hamiltonii as a polyphyletic evolving species in the data set. This might be a false positive grouping since RAPD markers are dominant and cannot discriminate if a DNA fragment is amplified from a homozygous locus or a heterozygous locus. By contrast, the species were discriminated by point mutations observed at the trnL-F intergenic spacer region (Figure S3). This phenomenon of low variation can possibly be explained by the unique asexual reproductive life cycle and also the high rate of transversion transition at the trnL-F locus. Thus, using Kumar and Gadagkar [39] disparity test of substitution pattern homogeneity rate, we showed that variations are generally low at the trnL-F intergenic spacer region (Figure S4). Nonetheless, based on two molecular tools used in this study, we found that C. callosa collected at high altitude (1843 m) had undergone a unique evolving pattern (Figures 3 and 4). Evidently, DNA sequence (Figure S3) shows that C. callosa has suffered less from point mutation in comparison with Bambusa species, Schizostachyum species, Dendrocalamus species, and Melocanna species.

Undoubtedly, since few bamboo species have been adapted for consumption through ancestral customs and not domesticated for extensive farming, no proper study has been conducted to identify low level TCC species vis-à-vis their geographical positions in a dynamic biota. In this era of high demand for bamboo shoots, irrespective of the health risk (such as Konzo), identifying species with nutritious attributes, nonetheless, with low level TCC is indispensable. A holistic approach encompassing morphological and molecular tools which are required, considering the three barcoding tools used in this study, showed different interspecies clustering pattern. From the present data, grappling the high altitude for safe edible bamboo species allowed for low input capital for eliminating TCC during bamboo shoot processing. Nonetheless, prolonged and acute intake of diet containing cyanogen derivatives could lead to death; exacerbate goitre, cretinism in iodine deficient regions, mental confusion, irreversible paralysis of legs, and ultimately the neurological disorder of Konzo [3, 39, 40]. Intriguingly, an optimal cooking procedure (at 98–102°C for 148–180 min) has shown to reduce cyanogens by 97% in D. giganteus [41, 42]. Since, only D. giganteus was studied [42], this cooking protocol might not be applicable for other bamboo species because the binding of taxiphyllin may vary from species to species and even at the genus level.

Accordingly, the safest procedure for consuming bamboo shoots might only be by using species with low TCC as starting material. Previously, Haque and Bradbury [9] showed that the total cyanide content decreased from the tip to the base of bamboo shoots which agreed with the present study, except for the Dendrocalamus species which further showed significant variability, irrespective of the altitude of collection. This suggests that the Dendrocalamus species should most probably not be consumed or considered at all as edible bamboo. Although cyanogenic, some succulent bamboo species are shown to possess a free radical scavenging activity which is considered as good food. Explicitly, some species showed appreciable amount of antioxidant activity except for the Dendrocalamus species, further confirming their low nutritious value as a whole. Recently, Park and Jhon [7] evaluated the antioxidant activities of Phyllostachys pubescens and Phyllostachys nigra only and concluded that the parameter varies between the species. Here, using a larger data set we showed that there is an interspecific variation in antioxidant properties of bamboo shoots collected at different altitudes.

The recommended dietary allowance (RDA) intake for magnesium (Mg) is 320 mg for female adults and 400–420 mg for male adults [43]. Consequently, 100 g of the bamboo shoots (except the Dendrocalamus species) can supply 28–75% of potassium, 1–4% of sodium, 10–22% of phosphorous and 22–56% of magnesium for females and 16–43% for adults, respectively. The RDA of iron for pregnant women is 27 mg, and thus 100 g of the bamboo shoots can provide 26–96% RDA for pregnant women. Equally, 100 g of the bamboo shoots also contained 28–136% of zinc and 2–7% of copper required as per RDA [43]. Based on our studies, C. callosa provides more of these beneficial RDA intakes of micro- and macronutrients as well as substantial antioxidant activity. Plants with antioxidant properties have been shown to be important natural antitumor and antimutagenic reserves [6, 44]. When put together, results of mineral content obtained in this study agreed with the previous quantification reports [8].

5. Conclusions

Considering the commercial implications, health risks associated with poisonous bamboo shoots, and nonavailability of optimized method for removal of TCC, we highlight the following for validating safe edible bamboo shoots prior to commercialization: (1) barcoding bamboo shoots using morphological descriptors and trnL-F intergenic spacer; (2) tagging the packed products with at least the genus name and beneficial nutritional components; (3) cogent tagging of genus Dendrocalamus with warning such as “bamboo shoots are injurious to health”; (4) refining the available edible bamboo shoots processing technique to minimize the residual quantity of TCC in the final products and normalizing the processing technique for stakeholders involved in the trade; (5) conserving germplasm, domesticating, and extensively farming C. callosa at high altitude of at least 1843 m from sea level to avoid extinction.

Conflict of Interests

The authors declare that there is no conflict of interests.

Acknowledgments

This research was supported by the Department of Biotechnology (DBT), Government of India, in the form of a Postdoctoral Fellowship to Sayanika Devi Waikhom. The author, Bengyella louis, is grateful to The Academy of Sciences for Developing World (TWAS) and DBT, Government of India (program no. 3240223450), for his fellowship. We extend our gratitude to D. K. Hore, D. Biswas and N. Mazumder for proofreading the manuscript.

Supplementary Materials

Table S1 shows qualitative and quantitative morphological descriptors. Table S2 shows list of primers use for performing random amplified polymorphic DNA analysis. Figure S1 shows RAPD gel profile using OPA-04 primer. Figure S2 shows sequence alignment entropy generated in BioEdit [19]. Table S3 shows principal component analysis for interrelatedness between nutritional attributes. Figure S3 shows comparative profile of point mutation at trnL-F intergenic spacer. Figure S4 shows test of homogeneity of substitution patterns between sequences.

  1. Supplementary Material

References

  1. F. Nartey, “Toxicological aspects of cyanogen’s is in tropical food stuff,” in Toxicology in the Tropics, R. L. Smith and E. A. Bababumni, Eds., pp. 53–73, Taylor and Francis, London, UK, 1980. View at: Google Scholar
  2. U. Schwarzmaier, “Cyanogenesis of Dendrocalamus: taxiphyllin,” Phytochemistry, vol. 16, no. 10, pp. 1599–1600, 1977. View at: Google Scholar
  3. H. Nzwalo and J. Cliff, “Konzo: from poverty, cassava, and cyanogen intake to toxico-nutritional neurological disease,” PLoS Neglected Tropical Diseases, vol. 5, no. 6, article e1051, 2011. View at: Publisher Site | Google Scholar
  4. M. D. Thiyam, S. S. Giri, D. Y. Pramoda, and D. G. A. Shantibala, “An investigation into elimination of cyanide from under process soibum,” The Bioscan, vol. 5, no. 4, pp. 639–642, 2010. View at: Google Scholar
  5. S. D. Waikhom, S. Ghosh, N. C. Talukdar, and S. S. Mandi, “Assessment of genetic diversity of landraces of Dendrocalamus hamiltonii using AFLP markers and association with biochemical traits,” Genetics and Molecular Research, vol. 11, no. 3, pp. 2107–2121, 2012. View at: Publisher Site | Google Scholar
  6. N. Chongtham, M. S. Bisht, and S. Haorongbam, “Nutritional properties of bamboo shoots: potential and prospects for utilization as a health food,” Comprehensive Reviews in Food Science and Food Safety, vol. 10, no. 3, pp. 153–168, 2011. View at: Publisher Site | Google Scholar
  7. E. J. Park and D. Y. Jhon, “The antioxidant, angiotensin converting enzyme inhibition activity, and phenolic compounds of bamboo shoot extracts,” LWT—Food Science and Technology, vol. 43, no. 4, pp. 655–659, 2010. View at: Publisher Site | Google Scholar
  8. C. Nirmala, E. David, and M. L. Sharma, “Changes in nutrient components during ageing of emerging juvenile bamboo shoots,” International Journal of Food Sciences and Nutrition, vol. 58, no. 8, pp. 612–618, 2007. View at: Publisher Site | Google Scholar
  9. M. R. Haque and J. H. Bradbury, “Total cyanide determination of plants and foods using the picrate and acid hydrolysis methods,” Food Chemistry, vol. 77, no. 1, pp. 107–114, 2002. View at: Publisher Site | Google Scholar
  10. Y. P. Kalra, “Determination of pH soils by different methods: collaborative study,” Journal of AOAC International, vol. 78, no. 2, pp. 310–321, 1995. View at: Google Scholar
  11. F. J. Rohlf, NTSYS-pc Numerical Taxonomy and Multivariate Analysis System, Version 2.20e, Exeter Publication, New York, NY, USA, 2000.
  12. P. H. A. Sneath and R. R. Sokal, Numerical Taxonomy: The Principles and Practice of Numerical Classification, W. H. Freeman, San Francisco, Calif, USA, 1st edition, 1973.
  13. R. R. Sokal and C. D. Michener, “A statistical method for evaluating systematic relationships,” University of Kansas Science Bulletin, vol. 28, pp. 1409–1438, 1958. View at: Google Scholar
  14. S. Aras, A. Duran, and G. Yenilmez, “Isolation of DNA for RAPD analysis from dry leaf material of some Hesperis L. Specimens,” Plant Molecular Biology Reporter, vol. 21, no. 4, pp. 461–462, 2003. View at: Google Scholar
  15. P. Jaccard, “Nouvelles recherché sur la distribution florale,” Bulletin de la Société Vaudoise de Sciences Naturelles, vol. 44, pp. 223–270, 1908. View at: Google Scholar
  16. P. Taberlet, L. Gielly, G. Pautou, and J. Bouvet, “Universal primers for amplification of three non-coding regions of chloroplast DNA,” Plant Molecular Biology, vol. 17, no. 5, pp. 1105–1109, 1991. View at: Publisher Site | Google Scholar
  17. F. Sievers, A. Wilm, D. Dineen et al., “Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega,” Molecular Systems Biology, vol. 7, article 539, 2011. View at: Publisher Site | Google Scholar
  18. T. A. Hall, “BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT,” Nucleic Acids Symposium Series, vol. 41, pp. 95–98, 1999. View at: Google Scholar
  19. I. Milne, F. Wright, G. Rowe, D. F. Marshall, D. Husmeier, and G. McGuire, “TOPALi: software for automatic identification of recombinant sequences within DNA multiple alignments,” Bioinformatics, vol. 20, no. 11, pp. 1806–1807, 2004. View at: Publisher Site | Google Scholar
  20. M. Hasegawa, H. Kishino, and T. Yano, “Dating of the human-ape splitting by a molecular clock of mitochondrial DNA,” Journal of Molecular Evolution, vol. 22, no. 2, pp. 160–174, 1985. View at: Google Scholar
  21. K. Tamura, D. Peterson, N. Peterson, G. Stecher, M. Nei, and S. Kumar, “MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods,” Molecular Biology and Evolution, vol. 28, no. 10, pp. 2731–2739, 2011. View at: Publisher Site | Google Scholar
  22. M. J. Bradbury, S. V. Egan, and J. H. Bradbury, “Determination of all forms of cyanogens in cassava roots and cassava products using picrate kits,” Journal of the Science of Food and Agriculture, vol. 79, no. 4, pp. 593–601, 1999. View at: Google Scholar
  23. L. L. Mensor, F. S. Menezes, G. G. Leitão et al., “Screening of Brazilian plant extracts for antioxidant activity by the use of DPPH free radical method,” Phytotherapy Research, vol. 15, no. 2, pp. 127–130, 2001. View at: Publisher Site | Google Scholar
  24. AOAC, Official Method of Analysis, Association of Official Agriculture Chemicals, Washington, DC, USA, 10th edition, 1985.
  25. R. A. Koenig and C. R. Johnson, “Colorimetric determination of phosphorus in biological materials,” Industrial and Engineering Chemistry, vol. 14, no. 2, pp. 155–156, 1942. View at: Google Scholar
  26. D. M. Updegraff, “Semimicro determination of cellulose inbiological materials,” Analytical Biochemistry, vol. 32, no. 3, pp. 420–424, 1969. View at: Google Scholar
  27. F. Tajima and M. Nei, “Estimation of evolutionary distance between nucleotide sequences,” Molecular Biology and Evolution, vol. 1, no. 3, pp. 269–285, 1984. View at: Google Scholar
  28. S. X. Lin and J. Lapointe, “Theoretical and experimental biology in one,” Journal of Biomedical Science and Engineering, vol. 6, pp. 435–442, 2013. View at: Publisher Site | Google Scholar
  29. G. P. Zhou and M. H. Deng, “An extension of Chou's graphic rules for deriving enzyme kinetic equations to systems involving parallel reaction pathways,” Biochemical Journal, vol. 222, no. 1, pp. 169–176, 1984. View at: Google Scholar
  30. K. C. Chou, “Graphic rules in steady and non-steady state enzyme kinetics,” Journal of Biological Chemistry, vol. 264, no. 20, pp. 12074–12079, 1989. View at: Google Scholar
  31. J. Andraos, “Kinetic plasticity and the determination of product ratios for kinetic schemes leading to multiple products without rate laws—New methods based on directed graphs,” Canadian Journal of Chemistry, vol. 86, no. 4, pp. 342–357, 2008. View at: Publisher Site | Google Scholar
  32. I. W. Althaus, A. J. Gonzales, J. J. Chou et al., “The quinoline U-78036 is a potent inhibitor of HIV-1 reverse transcriptase,” Journal of Biological Chemistry, vol. 268, no. 20, pp. 14875–14880, 1993. View at: Google Scholar
  33. K. C. Chou, “Graphic rule for drug metabolism systems,” Current Drug Metabolism, vol. 11, no. 4, pp. 369–378, 2010. View at: Publisher Site | Google Scholar
  34. K. C. Chou, W. Z. Lin, and X. Xiao, “Wenxiang: a web-server for drawing wenxiang diagrams,” Natural Science, vol. 3, no. 10, pp. 862–865, 2011. View at: Publisher Site | Google Scholar
  35. G. P. Zhou, “The disposition of the LZCC protein residues in wenxiang diagram provides new insights into the protein-protein interaction mechanism,” Journal of Theoretical Biology, vol. 284, no. 1, pp. 142–148, 2011. View at: Publisher Site | Google Scholar
  36. G. P. Zhou and R. B. Huang, “The pH-triggered conversion of the PrP(c) to PrP(sc.),” Current Topics in Medicinal Chemistry, vol. 13, no. 10, pp. 1152–1163, 2013. View at: Publisher Site | Google Scholar
  37. P. M. Kroonenberg, “Introduction to biplots for G × E tables,” Research Report 51, Department of Mathematics, University of Queensland, 1995. View at: Google Scholar
  38. W. Yan, M. S. Kang, B. Ma, S. Woods, and P. L. Cornelius, “GGE biplot versus AMMI analysis of genotype-by-environment data,” Crop Science, vol. 47, no. 2, pp. 643–655, 2007. View at: Publisher Site | Google Scholar
  39. S. Kumar and S. R. Gadagkar, “Disparity index: a simple statistic to measure and test the homogeneity of substitution patterns between molecular sequences,” Genetics, vol. 158, no. 3, pp. 1321–1327, 2001. View at: Google Scholar
  40. N. Mlingi, N. H. Poulter, and H. Rosling, “An outbreak of acute intoxications from consumption of insufficiently processed cassava in Tanzania,” Nutrition Research, vol. 12, no. 6, pp. 677–687, 1992. View at: Google Scholar
  41. O. S. A. Oluwole, A. O. Onabolu, I. A. Cotgreave, H. Rosling, A. Persson, and H. Link, “Incidence of endemic ataxic polyneuropathy and its relation to exposure to cyanide in a Nigerian community,” Journal of Neurology, Neurosurgery and Psychiatry, vol. 74, no. 10, pp. 1417–1422, 2003. View at: Publisher Site | Google Scholar
  42. V. L. Ferreira, K. Yotsuyanagi, and C. R. Carvalho, “Elimination of cyanogenic compounds from bamboo shoots Dendrocalamus giganteus Munro,” Tropical Science, vol. 35, no. 4, pp. 342–346, 1995. View at: Google Scholar
  43. Institute of Medicine (IOM), “US National Academy of Sciences, dietary reference intakes for individuals, food and nutritional board,” 2004. View at: Google Scholar
  44. M. Zahin, F. Aqil, F. M. Husain, and I. Ahmad, “Antioxidant capacity and antimutagenic potential of Murrayakoengii,” BioMed Research International, vol. 2013, Article ID 263509, 10 pages, 2013. View at: Publisher Site | Google Scholar

Copyright © 2013 Sayanika Devi Waikhom et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1500 Views | 808 Downloads | 4 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder