Biomass Production of Three Biofuel Energy Plants ’ Use of a New Carbon Resource by Carbonic Anhydrase in Simulated Karst Soils : Mechanism and Capacity

To determine whether the bicarbonate in karst limestone soil could be used as a new carbon resource for biomass production by the catalysis of carbonic anhydrase (CA), a simulative karst drought stress experiment was designed and performed. Three plants used for biofuel energy, Orychophragmus violaceus L. (Ov), Brassica juncea L. (Bj), and Euphorbia lathyris L. (El), were grown under simulated karst drought stress. In response to drought stress, the photosynthesis of the three energy plants was inhibited, but their CA activity increased. The hypothesis was confirmed by plant physiological and stable isotope techniques. The obtained results showed that plant biomass was produced with atmospheric CO2 as well as bicarbonate under drought stress. Bicarbonate use was proportional to the CA activity of the plants. With high CA activity over a long period, El had the highest proportional bicarbonate use compared to Ov and Bj, reaching 26.95%. Additionally, a new method is proposed for the screening of plants grown for energy in karst habitats.


Introduction
Currently, a deficiency in feedstock remains a principal barrier in the biomass industry.Due to limited land availability for cultivation, few countries can grow plants for biofuel without other considerations.Although marginal lands are not suited for agricultural production due to poor soil quality, harsh climates, rough terrain, and high levels of pollution, these lands are recommended for the production of biomass feedstock, which has garnered widespread attention [1][2][3][4][5][6][7][8].As an important marginal land resource, karst landscapes occupy approximately 11.2% of the surface of the earth.Additionally, the karst landscapes in Southwestern China represent more than 30% of the total potentially available land in China [9,10].In karst habitats, plants must tolerate stresses from high concentrations of bicarbonate, high pH, frequent arid-humid alternation, high concentrations of Ca 2+ and Mg 2+ , and phosphorus deficiencies.Although these factors adversely affect the survival and growth of plants, the growth of some plants was improved by the regulation of carbonic anhydrase Energies 2017, 10, 1370 2 of 14 (CA, EC 4.2.1.1)gene expression, including pH regulation, ion exchange, respiration, CO 2 transfer, and photosynthetic CO 2 fixation.These plants are classified as karst-adapted and can grow across the karst regions in limestone soils [9,11].
In terms of the accumulation of biomass in plants, most researchers have only focused on the incorporation of atmospheric CO 2 in photosynthesis.However, bicarbonate is another inorganic carbon source that can be used in photosynthesis, which is absorbed through the roots.In fact, plant growth in karst habitats is almost always adversely affected by high concentrations of HCO 3 − in the range of 3 mmol•kg −1 to 10 mmol•kg −1 , and high pH values ranging from 7.4 to 8.3 [10].In a recent study, the addition of an appropriate concentration of bicarbonate increased the biofuel production of microalgae (Scenedesmus sp.CCNM 1077) [12].As a class of key regulatory and detoxifying enzymes, CA helps plants to mitigate the adverse effects pH and HCO 3 − with rapid catalysis of the reversible reaction HCO 3 − + H + CO 2 + H 2 O.The resulting CO 2 and H 2 O are used in photosynthesis.Additionally, in the guard cells of plant leaves, CA effectively regulates gas exchange between plants and the atmosphere [11,[13][14][15].Both lower plants and higher plants, CA has been attracting attention in the biomass field.Using the catalytic activity of CA, algal reactors (Chlorella sp. and Spirulina platensis) also increased the reversible conversion between bicarbonate and CO 2 to develop biomass in a small raceway pond [16].Although different from lower plants (algage), the CA of higher plants was also found have the similar function.According to Gu et al. [17], the activities of phosphoenolpyruvate carboxylase and CA in two transgenic rice genotypes increased three-to five-fold compared to those in an untransformed wild type cultivar, with the same level of drought stress.Notably, the biomass production of the transgenic plants was 7%-15% higher than that of the wild type cultivar.Jiang et al. also studied the effect of CA activity on biomass accumulation.They found that the expression of mitochondria-localized CA increased when the respiration rate of leaves decreased.Moreover, with an increase in the expression of CA activity, shoot biomass accumulated in Arabidopsis thaliana [18].In our previous report [19], we found that microalgae promote the dissolution of limestone.Using a stable carbon isotope technique, we deduced that the CA in microalgae catalyzed the dissolution of bicarbonate from limestone for biomass production, with a three-to seven-fold increase in biomass production compared to those with an inhibited CA activity treatment.In the study of Swarnalatha et al. [20], biomass productivity and carbon dioxide bio-fixation increased two-to four-fold in treated cultures of indigenous freshwater microalgae compared to controlled cultures after 16 days.They attributed these results to the high CA activity in treatment cultures, which increased with the increase in pH and bicarbonate concentration, due to the high concentration of CO 2 in a closed photo-bioreactor.
Plants in the harsh karst habitats of Southwestern China also face the strain of spring drought, which typically occurs in March and April and lasts several weeks.In response to drought stress, karst-adapted plants close stomata and these plants use the bicarbonate in limestone soils as an alternative inorganic carbon source to maintain continuous photosynthesis [15,21].To distinguish between atmospheric CO 2 and exogenous bicarbonate sources of carbon, a stable carbon isotope technique was successfully used to trace the metabolism of inorganic carbon sources during photosynthesis [22].Using a method based on stable carbon isotopes, we quantified the proportion of bicarbonate use by plants [11].Although the conversion to biomass from bicarbonate by the activity of CA has been proposed in studies of plant biomass production for energy [17,18,20], the mechanism and the capacity for bicarbonate use in plants in karst habitats are not fully understood.In this study, we hypothesized that CA activity of plants increased under drought stress, which produced H 2 O and CO 2 from the bicarbonate (HCO 3 − ) in limestone dissolution of karst soils for use in photosynthesis and biomass production.With the integration of plant physiological and stable carbon isotope techniques, the effect of simulated drought stress in karst soil on plant physiological characteristics, biomass production, and capacity for bicarbonate use were investigated for three plant species used for biomass in biofuel energy production-Euphorbia lathyris L. (El), Orychophragmus violaceus L. (Ov), and Brassica juncea L. (Bj)-among which El belonged to Euphorbiaceae family (C4, annual herbs) that could thrive in drought, frost, and arid soils; while Ov and Bj belonged to Cruciferae family (C3, annual herbs), two kinds of common plants in karst area.We then attempted to elucidate the relationship between CA activity and bicarbonate use capacity of the three plants.Additionally, a new method is proposed to screen plants grown for energy in karst soils.

Plant Materials and Substrate
The El, Ov, and Bj seeds were provided by the Biomass Energy Research Institute of Guizhou Province, China, the Institute of Geochemistry, Chinese Academy of Sciences, China, and the Rapeseed Institute of Guizhou Province, China, respectively.The substrate used in plastic pots was a 1:1 quartz sand to vermiculite mixture (m/m) which was fertilized with modified Hoagland's nutrient solution (pH 8.1 ± 0.1).To quantify the proportion of bicarbonate use of the three species, the plants in each treatment level of drought stress were divided into two groups that received additions of NaHCO 3 with two different δ 13 C values, −27.28 Pee Dee Belemnite (PDB) or −18.81 PDB.The δ 13 C values of NaHCO 3 were determined by gas isotope ratio mass spectrometry.
After air drying, the substrate contained 6 mmol•kg  [10,14], the concentration of NaHCO 3 in the substrate was 10 mmol•kg −1 , and the pH value was 8.1.

Experimental Design and Growth Conditions
A systematic experiment was conducted in a growth chamber under natural sunlight in the Institute of Geochemistry, Chinese Academy of Sciences, Guizhou Province, China (26.35 • N, 106.42 • E).The experiment was divided into growth and stress phases.During the growth phase from 21 September 2013 (autumn) to 18 February 2014 (spring), totaling 150 days, the room temperature and relative humidity were balanced with outside conditions using ventilators and were recorded with an intelligent temperature and humidity recorder at 60-minute intervals.The three species were treated similarly; for example, El seeds were surface sterilized with 5% (v/v) H 2 O 2 for 30 min, washed with deionized water and sown with eight seeds per plastic pot (height: 24 cm; volume: 5 L), and filled with substrate.At 15 days after sowing, the pots were thinned to four uniform seedlings with good growth.A total of 256 plants in 64 pots, which were divided into two groups, were grown for 180 days in well-watered cultures.Based on gravimetric measurements of the pots, the moisture content of the substrate was maintained at 20 ± 1% with additions of NaOH solution (pH 8.1) every five days.Following the volumetric procedure used by Zuo et al. [23], the concentration of bicarbonate in the substrate at 10 cm was determined and modified 1/10 Hoagland's nutrient solution with the NaHCO 3 at one of two different δ 13 C values was added at 10-day intervals to maintain a concentration of 10 mmol•kg −1 .During the growth phase, the room temperature ranged from 268.9 K to 303.9 K and the relative humidity ranged from 48.7% to 76.2%, with average values of 284.4 K and 63.1%.There were 23 days below 273 K.
Following the growth phase of 150 days, a 21-day stress phase occurred from 19 February to 10 March 2014.Room temperature ranged from 279.1 K to 305.9 K, relative humidity from 63.7% to 75.2%, and the average values were 290.9K and 71.1%.Based on literature values for karst soils in drought conditions [4], four drought treatments were applied: well-watered (control, 20% moisture), mild drought stress (D 1 , 17% moisture), moderate drought stress (D 2 , 14% moisture), and severe drought stress (D 3 , 11% moisture).Therefore, the 64 pots were classified into 4 treatments.For each level of moisture, the 16 pots were further divided into 2 groups that were watered with NaOH solution that contained NaHCO 3 at one of two different δ 13 C values, −27.28 PDB or −18.81 PDB, every two days.The bicarbonate concentration of the substrate was maintained at 10 mmol•kg −1 for the entire phase of the experiment.

Physiological Characteristics
On days 1, 7, 14, and 21 of the drought stress phase, malondialdehyde (MDA), net photosynthetic rate P n (µmol•m 2 •s −1 ), and stomatal conductance G s (mmol•m 2 •s −1 ) were determined.To obtain accurate data for these biological characteristics, the fourth fully expanded leaf from the top of the plant was used.For MDA determinations, 1 g of leaf tissue obtained from 10 randomly selected plants was homogenized and analyzed by UV spectrophotometry [24].The P n and G s were determined between 9:30 and 11:00 a.m. using a portable LI-6400XT photosynthesis measurement system (LI-COR Inc., Lincoln, NE, USA).The photosynthetically active radiation, temperature, and CO 2 concentration measurements taken were 600 µmol•m −2 •s −1 , 298 K, and 400 µmol•mol −1 , respectively, which were determined from 10 random plants in different pots.

Biomass Production
At the end of the stress phase, 256 plants were harvested according to two differently labeled groups.Then were separated into four parts: roots, stems, leaves, and seeds.To determine dry biomass production (g), the samples were dried under a vacuum (335 K, 0.08 MPa) for 4 h.The oil content (wt %) and acid value (mg KOH/g) of seeds were determined according to ISO 659:1998 and ISO 660:1996, respectively.After methylation according to ISO 5509: 2001, the fatty acid composition (FAC) of oils was obtained using the external standard method and a Shimadzu GC/MS-QP2010E (Shimazhu, Kyoto, Japan) equipped with an HP-Innowax column (30 m × 0.32 mm, 0.5 µm) (Agilent, Palo Alto, CA, USA) and a flame ionization detector (Agilent, Palo Alto, CA, USA) [25].Furthermore, the fuel properties were measured by cetane number (CN), higher heating value (HHV), and kinematic viscosity (KV) from each oil sample which were determined using the equations of Ramos et al. [26,27].

CA Activity
CA activity (WAU g −1 (FW)) of the leaf tissue was determined every seven days during the stress phase using the method reported previously [11].The fourth youngest fully expanded leaf from the top of the plant was used.Immediately, leaf tissue samples, weighing 0.1 g-0.2 g, from each treatment were ground in liquid nitrogen with 3 mL of extraction buffer (0.01 M sodium barbitone with 0.05 M mercaptoethanol, pH 8.3).Homogenates were centrifuged at 10,000× g for 5 min at 273 K followed by standing for 20 min at 277 K.Then, 0.4 mL of homogenate was added to 4.5 mL of 0.02 M barbitone buffer (5,5-diethylbarbituric acid; pH 8.3).Using an electrochemical instrument (CHI 615C; Chenhua, Shanghai, China), the reaction was initialized with the addition of 3 mL of CO 2 -saturated deionized water.The CA activity was calculated using Equation (1) CA activity (WAU where t 0 and t are the reaction times (s) for the pH to change from 8.2 to 7.2 with buffer alone (t 0 ) and with the sample (t).

Stable Carbon-Isotope Ratio Analysis
Plants were cultivated on two kinds of simulate karst substrate, which containing different NaHCO 3 with different δ13C values (−27.28PDB or −18.81 PDB).On days 1, 7, and 14 of the drought stress phase, the first fully expanded leaves from the tops of 10 random plants of each treatment were collected.On day 21 of the stress phase, seed samples were collected.The seeds were immediately deactivated at 383 K for 20 min and were dried at 353 K for 8 h.Dried samples were ground to pass through a 100-mesh sieve.Powdered samples, weighing 3.0 mg-3.2 mg, were enclosed in a quartz tube and converted to CO 2 by combustion.The stable carbon isotope values (δ 13 C) were Energies 2017, 10, 1370 5 of 14 determined using a gas isotope ratio mass spectrometer (Mat-252; Finnigan MAT, San Jose, CA, USA).Stable carbon isotope values (δ) were expressed relative to the C reference materials (PDB) as and were calculated using Equation ( 2)

Proportion of Bicarbonate Use
In this study, the exogenous dissolved inorganic carbon from added NaHCO 3 (NaHCO 3 labels) was the focus.To minimize the contribution of other sources of isotope fractionation, such as photosynthetic fractionation and stomatal conductance, a bidirectional stable carbon isotope tracer technology was employed under the same culture conditions.Based on the results of stable carbon isotope ratio analyses, the proportion of bicarbonate use was quantified for the three plant species under drought stress.A bivariate isotope-mixture model was used, shown in Equation (3), which is described in detail in our previous report [11] where f B is the proportion of bicarbonate use, δ T1 and δ T2 are the δ 13 C values of samples from the two differently labeled groups at identical moisture levels, and δ C1 and δ C2 are the δ 13 C values of the two stable isotope-labeled NaHCO 3 additions.

Statistical Analyses
All data are expressed as the mean values ± standard deviation, with the exception of biomass and seed production.Statistical analyses were performed using the SPSS 13.0 statistical software package (SPSS Inc., Chicago, IL, USA).Differences amongst means were compared with Duncan's Multiple Range Tests (DMRT) at a significance level of p < 0.05.

MDA Content of Plants Leaves Under Drought Stress
Malondialdehyde (MDA) is the main product of membrane lipid peroxidation, and its content is proportional to the degree of damage to cells.Therefore, it is frequently used as the indicator in studies on stress of plants.In this study, the MDA content in the leaves of the three energy plant species was determined every seven days.As shown in Figure 1, under drought stress (D1: 17% moisture, D2: 14% moisture, and D3: 11% moisture), the MDA content in the leaves of the three species increased significantly compared with that in the respective control.On day seven, compared to each controlled treatment, the MDA content in Ov leaves in the D 2 and D 3 treatments increased by 246.17% and 354.23%, respectively, which were larger increases than those in the leaves of Bj (140.47% and 172.41%, respectively) and El (44.43% and 160.61%, respectively).However, as the experiment continued, the MDA content of the Bj leaves in the three drought stress treatments increased dramatically by 252.29% to 272.13% on day 21 compared to the controls, whereas the increase for Ov and El leaves ranged from 56.50% to 208.51%.For the entire stress phase, the lowest accumulation of MDA occurred in El leaves, compared with those of Ov and Bj for the identical stress treatments.Based on the strong contrast in accumulation of MDA in the leaves of the three species, we speculated that Bj was more vulnerable to simulated drought stress in karst soils than Ov and El.
increased dramatically by 252.29% to 272.13% on day 21 compared to the controls, whereas the increase for Ov and El leaves ranged from 56.50% to 208.51%.For the entire stress phase, the lowest accumulation of MDA occurred in El leaves, compared with those of Ov and Bj for the identical stress treatments.Based on the strong contrast in accumulation of MDA in the leaves of the three species, we speculated that Bj was more vulnerable to simulated drought stress in karst soils than Ov and El.

Photosynthesis of Plant Leaves under Drought Stress
Photosynthesis plays an important role in plant growth and biomass production, and stomatal closure is a common mechanism to avoid loss of water through transpiration when plants suffer from drought stress.Thus, drought-induced stomatal closure could inhibit photosynthesis.Net photosynthetic rate (P n (µmol m 2 •s −1 )), and stomatal conductance (G s (mmol m 2 s −1 )) were employed as direct evaluation indexes [15].Two important parameters of photosynthesis, G s and P n , for the three plants are displayed in Figure 2. Compared to those in the control, the values of G s and P n of the three energy plants declined at different rates with continuous drought stress.After drought stress was initiated, the P n values of El were higher than Ov and Bj at the same level of stress.Bj had the lowest P n value of the three plant species, except on the seventh day.A similar phenomenon was also observed in the G s of the three plant species.At the beginning of the stress phase on day seven, P n and G s values of both Ov and Bj reduced significantly under the D 1 , D 2 , and D 3 treatments (p < 0.05); However, the P n and G s of El did not decrease significantly (Figure 2C,F) under mild stress (D 1 ) when compared to the controls during the entire stress phase.For Bj on day seven, under mild drought stress (D 1 ), G s and P n decreased 53.69% and 31.77%,respectively, compared to the controls, and the parameters for Ov declined by 42.98% and 22.43%, respectively.Under the mild stress treatment at the end of the experiment on day 21, these two indices for Bj decreased by 85.08% and 83.06%, and those of Ov decreased 47.96% and 81.70%, respectively; whereas those of El declined only 1.14% and 1.78%, respectively.On day 21 of severe stress (D 3 ), both G s and P n values of Bj and Ov decreased sharply by more than 96%; by contrast, those values for El decreased 89.14% and 89.15%, respectively.In contrast to Ov and Bj, El is in the Euphorbiaceae family, which can thrive drought and saline environments [25].When combined with the results for MDA accumulation (Figure 1), we concluded that the drought resistance of El was superior to that of Ov and Bj under simulated drought stress in karst soil.only 1.14% and 1.78%, respectively.On day 21 of severe stress (D3), both Gs and Pn values of Bj and Ov decreased sharply by more than 96%; by contrast, those values for El decreased 89.14% and 89.15%, respectively.In contrast to Ov and Bj, El is in the Euphorbiaceae family, which can thrive drought and saline environments [25].When combined with the results for MDA accumulation (Figure 1), we concluded that the drought resistance of El was superior to that of Ov and Bj under simulated drought stress in karst soil.

Biomass Production
Plants are an important biomass feedstock because of their biochemical composition and their biofuel-like compounds such as cellulose, hemicelluloses, lignin, and lipids, from root, stem, leaf, and seed, and so on.In this study, dry biomass production (g), oil content (wt %), and acid value of the oil samples (mg KOH/g) were determined, as shown in Table 1.In this experiment, the photosynthesis of three plants was inhibited when drought stress was imposed (Figure 2).Therefore, differences were observed in dry biomass production (g), oil content (wt %), acid value (mg KOH/g) of oil samples and seed production, especially in seed production (Table 1).Under the three drought treatments, Bj (111.38 g) and Ov (127.55 g) had higher average total biomass production than that of El (81.65 g).Compared to each control, El showed the lowest rate of decline (5.38%) compared with those of Bj (12.87%) and Ov (8.86%).When severe drought (D 3 ) was imposed on Bj, seed production decreased sharply to 20.68% of the control, whereas those of Ov and El were 49.50% and 61.69%, respectively.Based on seed production, we concluded that El had the highest biomass production among the three plant species under drought stress, which we attributed to the higher P n of El under all drought stress treatments (Figure 1).However, the total accumulation of biomass of Ov and Bj was higher during the growth phase than that of El.The accumulation of oil in the seeds of the three species was also inhibited at varying degrees of drought stress.However, the oil content of El seeds was not significantly affected in the D 1 (29.13 wt % ± 1.51 wt %) and D 2 (28.83 wt % ± 1.19 wt %) treatments compared to the controls (29.83 wt % ± 1.33 wt %).
Acid value is an important indicator in the evaluation of biodiesel feedstock and seed ripeness, which is proportional to the free fatty acid content of oil.As shown in the Table 1, the acid values for all the plant seed samples were affected by drought stress.Similar to MDA content and P n values (Figures 1 and 2), the acid values of Ov and Bj, under all levels of drought, were significantly different (P < 0.05).However, the acid values of El seed oil were significantly affected only under the D 2 and D 3 treatments (P < 0.05) compared to controlled and D 1 tretments.Thus, at the same levels of drought stress, the ripeness of Ov and Bj seeds was inhibited more than that of El seeds.Additionally, the acid values of all the oil samples were higher than 3.0 mg KOH/g, which prevents the application of these materials in an alkali-catalyzed biodiesel process [25].
Table 1.Effect of different levels (Con: 25% moisture, D 1 : 17% moisture, D 2 : 14% moisture, and D 3 : 11% moisture) of drought stress on the biomass production of three plant species.Different letters indicate a significant difference between treatments of the same plant (p < 0.05) and different lowercase letters in the same row denote significant differences among four drought treatments (p < 0.05).Drought stress at different levels also significantly affected the FAC of the three energy plants, compared with the respective controls (Tables A1 and A2).Overall, in the seed oil samples of the three species, with the increase in stress, monounsaturation (Cn:1) and saturation (Cn:0) increased and polyunsaturation (Cn:2,3) and the degree of unsaturation decreased.For example, oleic acid (C18:1) and linoleic acid (C18:2) were the primary fatty acids in Ov seed oil samples, and in drought-stressed treatments, oleic acid increased 5.47%, from 22.47% in the controls to 27.94% in the D 3 treatment, and linoleic acid decreased 8.64% from 42.58% in controls to 33.94% in the D 3 treatment.In some previous studies, recommendations for ideal biodiesel included a high proportion of monounsaturated fatty acids (C16:1 and C18:1), a reduced proportion of polyunsaturated acids, and a controlled content of saturated acids [26,27].CN is related to the ignition delay time and combustion quality.The higher the cetane number is, the better it is in its ignition properties [27].HHV also regarded as the gross calorific value or gross energy, it is usually used to define the energy content of fuels.KV is a property of the fluid which opposes the relative motion between the two surfaces of the fluid in a fluid, low KV is better [28].As shown in Table A2, the CN and HHV (MJ/kg) of the oil samples from the three species improved with the increase in stress level as a whole.However, the KVs (mm 2 /s; 40 • C) deteriorated.Based on the decrease in the degree of unsaturation, we predicted an improvement in the oxidation stability of all resulting biodiesel samples obtained from stressed cultures.

CA Activity and Proportion of Bicarbonate Use
Carbonic anhydrases are key regulatory and detoxifying enzymes for plants that are upstream regulators of CO 2 -controlled stomatal movements and are also an important tool for plants to assimilate HCO 3 − [12,14].In this experiment, the CA activity profile of the three species was assayed at seven-day intervals during the drought stress phase.Drought stress clearly modulated the CA activity of the three energy plants (Figure 3).On day 1, the CA activity of the three plant species was in the range of 4580 WAU g −1 (FW)-5520 WAU g −1 (FW), and Bj had the highest CA activity at 5520.07 WAU g −1 (FW) ± 178.29 WAU g −1 (FW).On day seven, CA activity increased in all cultures under drought stress, with the exception of that of Bj in the D 3 treatment.In sharp contrast to the control, CA activity was 8994 WAU g −1 (FW) ± 195 WAU g −1 (FW) for Bj in the D 2 treatment, which was an increase of 38.62%, whereas the increase was 30.88% and 16.73% for Ov and El, respectively.However, CA activity was inhibited when Ov and Bj suffered from severe stress (D 3 ); for example, the CA Energies 2017, 10, 1370 9 of 14 activity of Bj decreased 81.83% in the D 3 treatment compared to the control.As the experiment continued, the CA activity for each of the three plant species was different in response to continuous stress.For example, the CA activity of Ov on day 14 was 6012 WAU g −1 (FW) ± 139 WAU g −1 (FW) (control), 6570 WAU g −1 (FW) ± 235 WAU g −1 (FW) (D 1 ), 7950 WAU g −1 (FW) ± 297 WAU g −1 (FW) (D 2 ), and 2904 WAU g −1 (FW) ± 190 WAU g −1 (FW) (D 3 ).However, at the end of the stress phase on day 21, only the CA activity of the D 1 treatment was higher than that of the control at 6201 WAU g −1 (FW) ± 75 WAU g −1 (FW) vs. 5760 WAU g −1 (FW) ± 151 WAU g −1 (FW), with the other two treatments, D 2 and D 3 , showing a decrease in activity.When combined with the results for MDA content (Figure 1), we attributed the inhibition of activity to oxidative damage to plant cells from drought stress.Notably, the CA activity of El leaves in all treatments was higher than that in controls, with the exception of the leaves from the D 3 treatment on day 21.
Energies 2017, 10, 1370 9 of 14 To investigate the capacity of the three energy plants for bicarbonate use, a stable isotope technique was used.The δ 13 C values of the leaves from the three species were determined regularly during the stress phase, and are listed in the Table A3.We observed that the δ 13 C values of the leaves showed a difference under different levels of drought.Then, the proportion of bicarbonate use by the three energy plants from all treatments was calculated using Equation (3), with the results shown in Figure 3.After the 180-day growth phase, the proportion of bicarbonate use in leaves of the three species was different: 11.45% for Ov, 10.39% for Bj, and 9.44% for El.After a 21-day drought stress phase, the highest proportion of bicarbonate use was reached for each species: 26.72% for Ov in D3, 22.10% for Bj in D2, and 26.95% for El in D3.Based on these results, under continuous drought stress, the accumulation of inorganic carbon derived from bicarbonate increased compared to the controls.Notably, El showed the highest capacity for bicarbonate use compared to Ov and Bj.
In this study, several distinct phenomena were observed.(1) MDA content (μmol/kg FW) in the leaves of the three species increased with the severity of stress and with time during the stress phase (Figure 1), which indicated that cell damage increased in the leaves of these three species with the increase in drought stress.(2) drought-induced stomatal closure in the three plants led to the inhibition of photosynthesis (Figure 2).( 3) With the inhibition of photosynthesis, biomass production was effected, particularly that of seed production (Table 1), and compared to the respective controls, the decrease in seed production in descending order was Bj, Ov, and El.Additionally, the acid value and the degree of saturation of seed oils increased with the severity of drought stress.( 4) CA activity of the leaves of the three species was clearly affected by drought stress, but the activity increased under drought stress, compared to the controls (Figure 3).( 5) The To investigate the capacity of the three energy plants for bicarbonate use, a stable isotope technique was used.The δ 13 C values of the leaves from the three species were determined regularly during the stress phase, and are listed in the Table A3.We observed that the δ 13 C values of the leaves showed a difference under different levels of drought.Then, the proportion of bicarbonate use by the three energy plants from all treatments was calculated using Equation (3), with the results shown in Figure 3.After the 180-day growth phase, the proportion of bicarbonate use in leaves of the three species was different: 11.45% for Ov, 10.39% for Bj, and 9.44% for El.After a 21-day drought stress phase, the highest proportion of bicarbonate use was reached for each species: 26.72% for Ov in D 3 , 22.10% for Bj in D 2 , and 26.95% for El in D 3 .Based on these results, under continuous drought stress, the accumulation of inorganic carbon derived from bicarbonate increased compared to the controls.Notably, El showed the highest capacity for bicarbonate use compared to Ov and Bj.
In this study, several distinct phenomena were observed.(1) MDA content (µmol/kg FW) in the leaves of the three species increased with the severity of stress and with time during the stress phase (Figure 1), which indicated that cell damage increased in the leaves of these three species with the increase in drought stress.(2) drought-induced stomatal closure in the three plants led to the inhibition of photosynthesis (Figure 2).( 3) With the inhibition of photosynthesis, biomass production was effected, particularly that of seed production (Table 1), and compared to the respective controls, the decrease in seed production in descending order was Bj, Ov, and El.Additionally, the acid value and the degree of saturation of seed oils increased with the severity of drought stress.(4) CA activity of the leaves of the three species was clearly affected by drought stress, but the activity increased under drought stress, compared to the controls (Figure 3).(5) The bicarbonate use capacity of the three species was also different, and the highest proportion of bicarbonat use for Ov, Bj, and El occurred in those treatments with the highest CA activity for that species.
When these phenomena were examined collectively, most notably, the decreases were disproportionate for photosynthesis (P n and G s ) and biomass production.In sharp contrast to the respective controls on day 21, the P n values of Bj in the D 2 and D 3 treatments decreased by 98.97% and 99.41% (Figure 1 and Table 1), with decreases in seed production of 71.97% and 79.31%, respectively.For Ov with the same treatments, the decreases were 96.23% and 99.12% for P n , and 25.20% and 50.50% for seed production, respectively; and for El, the decreases were 89.49% and 94.67% for P n , and 13.48% and 38.80% for seed production, respectively.Both Bj and Ov are members of the Brassicaceae, and although the decrease in P n was approximately similar in the D 2 and D 3 treatments, the decreases in seed production were significantly different compared to the respective controls.We hypothesized that this result could be explained, at least in part, by the high CA activity of Ov and Bj under drought stress compared to the respective controls.For example, the CA activity of Bj in the D 2 treatment gradually decreased from day 7 to 21 and was lower than that in the control from day 14, indicating no further capacity to regulate photosynthesis.However, the CA activity of Ov in the D 2 treatment was higher than that in the control on days 7 and 14 and decreased gradually after day 14.Combined with a high proportion of bicarbonate use at 22.10% (Figure 3), the primary conclusion was that, with the regulation and catalytic activity of CA, the conversion of bicarbonate in the substrate into biomass increased.Thus, the decrease in Ov seed production was lower than that of Bj at the same stress level.Similarly, compared to the controls, the CA activity of Bj in the D 1 treatment and El in the three drought stress treatments were maintained at a higher level, with the exception of El in the D 3 treatment on day 21.The leaves of Bj in the D 1 and El in the D 3 treatment had higher proportion of bicarbonate use at the end of the drought stress compared to Ov, reaching 26.72% and 26.95%.Moreover, the lowest decrease in seed production under drought stress occurred with El (Table 1).In summary, through the regulation and catalytic activity of CA, the three species of biofuel energy plants could continue to produce biomass from atmospheric CO 2, but when under drought stress, alternatively, could do the same with bicarbonate.Additionally, the capacity for bicarbonate use of the three species was proportional to the CA activity under simulated drought stress in karst soil.Here, the present work just establishes the relationship between CA activity and bicarbonate-use capacity of the three energy plants.The mechanism and pathway for bicarbonate use in plants in karst habitats are not fully understood.We continue research by means of molecular biological technique in the future.

Conclusions
In karst habitats, high soil concentrations of bicarbonate are almost universally an adverse factor for the growth of plants.Drought is another stress factor for plants in the karst region of Southwestern China.In this study, three species of biofuel energy plants, Orychophragmus violaceus L. (Ov), Brassica juncea L. (Bj), and Euphorbia lathyris L. (El), were grown in soil simulating a karst habitat.As a result of integrating plant physiological and stable isotope techniques, the following results were obtained: (1) as a new carbon resource for biomass production, the use of bicarbonate in karst limestone soil was directly confirmed; (2) photosynthesis of the three plants was inhibited by drought stress, but the activity of CA increased; thus, the three energy plants could continue photosynthesis for biomass production with the alternate use of atmospheric CO 2 and bicarbonate through the regulation and catalysis of CA; and (3) the capacity for bicarbonate use of the three energy plants was proportional to the activity of CA.Among the three species, the highest bicarbonate use proportion was reached by El, and given the high CA activity, the bicarbonate use proportion reached 26.95% under severe drought stress.To our knowledge, few studies have focused on the contribution of this carbon resource in karst soils toward biomass production or on the mechanism for this use.Additionally, a new method is proposed for the screening of energy plants for use in karst soils.It is worth considering that the adapted plants with high carbonic anhydrase activity and high bicarbonate-use capacity were used for energy plants in karst soils.

Table A1 .
Fatty acid composition (FAC) of three energy plant species oil samples.
a Not detected because the seed sample of Bj under D 3 was too small.

Table A2 .
Effect of different levels (Con: 25% moisture, D 1 : 17% moisture, D 2 : 14% moisture, and D 3 : 11% moisture) of drought stress on the fatty acid composition and fuel properties of the oil samples from the three plant species.Calculated according to the report of Romas et al.; b Not detected because the seed sample of Bj under D 3 was too small. a