BioMed Research International

Volume 2013 (2013), Article ID 289285, 11 pages

http://dx.doi.org/10.1155/2013/289285

Research Article

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

Sayanika Devi Waikhom,¹ Bengyella Louis,1,2,3 Chandradev K. Sharma,¹ Pushpa Kumari,⁴ Bharat G. Somkuwar,¹ Mohendro W. Singh,¹ and Narayan C. Talukdar¹

¹Institute of Bioresources and Sustainable Development (IBSD), Takyelpat, Imphal, Manipur 795001, India

²Department of Biochemistry, University of Yaoundé I, BP 812, Yaoundé, Cameroon

³Department of Biotechnology, Burdwan University, Golapbag More, West Bengal 713104, India

⁴AJC Bose Indian Botanic Garden, Botanical Survey of India, Botanic Garden, Howrah 711103, India

Received 18 June 2013; Accepted 3 October 2013

Academic Editor: Kuo-Chen Chou

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.

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, 7–9], 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].

Table 1: Principal geographical coordinates, mean soil pH, and identified edible bamboo species.

GenBank accession/bamboo species
1
2
3
4
5
Voucher 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, ²mean pH of the soil at the collection site in the months of July-August of 2009 to 2011, ³latitude (N), ⁴longitude (E), and ⁵altitude 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-it² 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 AICc₁ and AICc₂, 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 [4–6]. 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% Na₂CO₃) 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 (IC₅₀) is defined as the amount of extract required to reduce free scavenging activity by 50%. The IC₅₀ 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 (HNO₃ : HClO₄ : H₂SO₄). 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.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

(m)

(n)

(o)

Figure 1: Morphological features of 30-day-old bamboo shoots imaged with Nikon Coolpix S6200 and the nomenclature as submitted in the Genbank NCBI nucleotide database: (a) KC013282/C. callosa, (b) KC013285/B. cacharensis, (c) JX564900/B. manipureana, (d) JX564901/B. nutans, (e) JX507132/B. tulda, (f) JX507131/B. oliveriana, (g) JX564902/D. giganteus, (h) JX564903/D. hamiltonii, (i) JX564906/Bambusa sp., (j) JX564904/D. hookeri, (k) JX564905/D. manipureanus, (l) JX564907/Bambusa sp., (m) KC013288/B. tuldoides, (n) JX507133/M. baccifera, and (o) JX507134/S. dullooa.

Figure 2: A dendrogram based on morphological descriptors showing the relationship between 15 edible bamboo species generated in NTSYS-PC software computed based on simple matching coefficient [11].

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.

Figure 3: A dendrogram based on Jaccard’s similarity coefficient obtained from RAPD data showing the relationship between 15 edible bamboo species.

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.

Figure 4: A maximum likelihood tree based on Hasegawa et al. [20] nucleotide substitution model. The tree is drawn to scale, with branch lengths measured in the number of nucleotide substitutions per site. Evolutionary analyses were conducted in MEGA5 [21].

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).

Table 2: Total cyanide content (TCC) expressed in part per million (ppm) for different portions of bamboo shoots and their net antioxidant activity expressed as inhibition concentration (IC₅₀) of DPPH.

Samples
Total cyanide content (ppm)*
IC₅₀ (mg L⁻¹)

Tip
Middle
Base

KC013282/Chimonobambusa callosa
300a
210a
199a
0.09ab

KC013285/Bambusa cacharensis
1533d
1221f
735efgh
0.58h

JX564900/Bambusa manipureana
1007b
515b
761fgh
0.46fg

JX564901/Bambusa nutans
1001b
624c
267ab
1.83i

JX507132/Bambusa tulda
1579d
1406hi
779gh
0.30ef

JX507131/Bambusa oliveriana
1280c
1079e
543c
0.57gh

JX564902/Dendrocalamus giganteus
2604g
2243k
920i
0.60h

JX564903/Dendrocalamus hamiltonii
1897e
766d
654cdef
0.14bc

JX564904/Dendrocalamus hookeri
1595d
1322fgh
360b
0.45fg

JX564905/Dendrocalamus manipureanus
1838e
1270fg
600cd
0.32de

JX507133/Melocanna baccifera
1548d
484b
216a
0.41ef

JX507134/Schizostachyum dullooa
1521d
1279fg
727efgh
0.61h

JX564906/Bambusa sp.
1582d
1354gh
669defg
0.24cd

JX564907/Bambusa sp.
2063f
1557j
635cde
0.25cd

KC013288/Bambusa tuldoides
2528g
1511ij
825hi
0.45fg

LSD ()
85.78
106.53
109.31
0.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 (IC₅₀) 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 IC₅₀ value of 0.60. L-Ascorbic acid investigated under the same conditions had an IC₅₀ value of 0.003 mg/mL; implying that bamboo shoots have a moderate antioxidant activity (Table 2). Based on IC₅₀, 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).

Figure 5: Antioxidant activity of 15 bamboo shoot extracts against DPPH free radical with respect to L-ascorbic acid as reference antioxidant agent.

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).

Table 3: The data represented are obtained from atomic absorption spectrometry analysis showing the different mineral concentrations in milligram per 100 g dry weight (d.w) of bamboo shoots.

Samples
N
P
Ce
K
Na
Mg
Ca
Fe
Cu
Zn

KC013282/C. callosa
1153m
154g
4070.75ab
3377e
39bc
110f
63.43gh
9.25abcd
6.12g
7.2de

KC013285/B. cacharensis
815g
87ab
3157.82a
2093abc
45cde
74ab
29.60bc
13.15def
3.58de
5.63bcd

JX564900/B. manipureana
957k
94bc
5936.51abc
2343bcd
53e
175i
32.4bcd
16.5f
4.94f
12.80i

JX564901/B. nutans
825gh
104cd
5372.79abc
2090abc
32b
175i
53.5fg
11.6cde
4.16e
6.50cde

JX507132/B. tulda
721c
116de
3701.59ab
1453abc
51de
92d
42de
25.8g
3.19cd
15.0j

JX507131/B. oliveriana
763e
85ab
9972.34cd
2293abc
43cd
123g
60.95fg
10.1bcd
3.58de
4.88bc

JX564902/D. giganteus
673a
70a
7048.98abcd
3093de
22a
69a
37.9cd
7.83abc
2.47bc
10.7gh

JX564903/D. hamiltonii
747d
118de
4305.22ab
1910abc
50de
121g
24.4ab
16.82f
3.0cd
7.87ef

JX564904/D. hookeri
692b
75a
5697.05abc
2400cd
45cde
84cd
73.3hi
6.08a
2.64bc
3.03a

JX564905/D. manipureanus
894i
83ab
6784.58abc
3533e
42cd
85cd
14.98a
6.68ab
2.55bc
4.80b

JX507133/M. baccifera
1102l
135f
10126.98cd
1757abc
39bc
149h
51.88ef
13.37ef
1.96ab
11.93hi

JX507134/S. dullooa
837h
121ef
7098.87abcd
1310a
42cd
180i
56.95fg
14.42ef
1.55a
9.18fg

JX564906/Bambusa sp.
785f
73a
8306.12bcd
1373ab
32b
92d
40.1cd
11.5cde
2.54bc
11.33g

JX564907/Bambusa sp.
940j
129ef
11044.9d
1523abc
53e
112f
50.8ef
13def
7.85h
11.03j

KC013288/B. tuldoides
900i
122ef
4988.66ab
2060abc
31b
101e
78.15j
11.2cde
2.7c
11.0h

LSD ()
13.32
14.92
2858.21
873.06
8.51
8.28
10.17
3.48
0.64
1.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 [29–31], 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].

Figure 6: PCA polygonal biplot analysis showing the interrelatedness of nutritional attributes of 12 edible bamboos species. The numbers represent bamboo species and the vectors are biochemical traits: KC013285/B. cacharensis (1), JX564900/B. manipureana (2), JX564901/B. nutans (3), JX507132/B. tulda (4), JX507131/B. oliveriana (5), JX564902/D. giganteus (6), JX564903/D. hamiltonii (7), JX564906/Bambusa sp. (8), JX564904/D. hookeri (9), JX564905/D. manipureanus (10), JX564907/Bambusa sp. (11), and KC013288/B. tuldoides (12). AA: antioxidant activity, Ca: calcium, Ce: cellulose, Cu: copper, Fe: iron, K: potassium, Mg: magnesium, N: nitrogen, Na: sodium, P: phosphorous, TCC: total cyanide content, and Zn: zinc.

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.

References

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.

U. Schwarzmaier, “Cyanogenesis of Dendrocalamus: taxiphyllin,” Phytochemistry, vol. 16, no. 10, pp. 1599–1600, 1977.

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.

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.

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.

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.

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.

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.

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.

Y. P. Kalra, “Determination of pH soils by different methods: collaborative study,” Journal of AOAC International, vol. 78, no. 2, pp. 310–321, 1995.

F. J. Rohlf, NTSYS-pc Numerical Taxonomy and Multivariate Analysis System, Version 2.20e, Exeter Publication, New York, NY, USA, 2000.

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.

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.

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.

P. Jaccard, “Nouvelles recherché sur la distribution florale,” Bulletin de la Société Vaudoise de Sciences Naturelles, vol. 44, pp. 223–270, 1908.

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.

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.

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.

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.

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.

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.

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.

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.

AOAC, Official Method of Analysis, Association of Official Agriculture Chemicals, Washington, DC, USA, 10th edition, 1985.

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.

D. M. Updegraff, “Semimicro determination of cellulose inbiological materials,” Analytical Biochemistry, vol. 32, no. 3, pp. 420–424, 1969.

F. Tajima and M. Nei, “Estimation of evolutionary distance between nucleotide sequences,” Molecular Biology and Evolution, vol. 1, no. 3, pp. 269–285, 1984.

S. X. Lin and J. Lapointe, “Theoretical and experimental biology in one,” Journal of Biomedical Science and Engineering, vol. 6, pp. 435–442, 2013.

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.

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.

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.

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.

K. C. Chou, “Graphic rule for drug metabolism systems,” Current Drug Metabolism, vol. 11, no. 4, pp. 369–378, 2010.

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.

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.

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.

P. M. Kroonenberg, “Introduction to biplots for G × E tables,” Research Report 51, Department of Mathematics, University of Queensland, 1995.

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.

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.

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.

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.

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.

Institute of Medicine (IOM), “US National Academy of Sciences, dietary reference intakes for individuals, food and nutritional board,” 2004

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.