ORIGINAL RESEARCH article
Down-Regulation of KORRIGAN-Like Endo-β-1,4-Glucanase Genes Impacts Carbon Partitioning, Mycorrhizal Colonization and Biomass Production in Populus
- 1BioEnergy Science Center and Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, USA
- 2BioEnergy Science Center, School of Chemistry and Biochemistry, Institute of Paper Science and Technology, Georgia Institute of Technology, Atlanta, GA, USA
- 3The Biosciences Center, National Renewable Energy Laboratory, Golden, CO, USA
- 4Oak Ridge National Laboratory – Department of Chemical and Biomolecular Engineering and Department of Forestry, Wildlife and Fisheries, University of Tennessee, Knoxville, TN, USA
A greater understanding of the genetic regulation of plant cell wall remodeling and the impact of modified cell walls on plant performance is important for the development of sustainable biofuel crops. Here, we studied the impact of down-regulating KORRIGAN-like cell wall biosynthesis genes, belonging to the endo-β-1,4-glucanase gene family, on Populus growth, metabolism and the ability to interact with symbiotic microbes. The reductions in cellulose content and lignin syringyl-to-guaiacyl unit ratio, and increase in cellulose crystallinity of cell walls of PdKOR RNAi plants corroborated the functional role of PdKOR in cell wall biosynthesis. Altered metabolism and reduced growth characteristics of RNAi plants revealed new implications on carbon allocation and partitioning. The distinctive metabolome phenotype comprised of a higher phenolic and salicylic acid content, and reduced lignin, shikimic acid and maleic acid content relative to control. Plant sustainability implications of modified cell walls on beneficial plant-microbe interactions were explored via co-culture with an ectomycorrhizal fungus, Laccaria bicolor. A significant increase in the mycorrhization rate was observed in transgenic plants, leading to measurable beneficial growth effects. These findings present new evidence for functional interconnectedness of cellulose biosynthesis pathway, metabolism and mycorrhizal association in plants, and further emphasize the consideration of the sustainability implications of plant trait improvement efforts.
Plant cell walls play integral roles in determining plant form and function. Changes in the composition and structural properties of cell walls can have significant implications on the developmental and physiological processes relevant to water transport (Chen et al., 2005; Liang et al., 2010; Zhu et al., 2010), abiotic and biotic stress responses (Zhao and Dixon, 2014), and plant fitness (Yang et al., 2007; Goubet et al., 2009; Park et al., 2010; Thévenin et al., 2011). The relationship between plant cell wall integrity and biotic stress response, especially to pathogenic microbes, is well documented (Cantu et al., 2008; Bellincampi et al., 2014; Miedes et al., 2014). Genetic modification of the biosynthesis pathways of the major cell wall components; cellulose (Caño-Delgado et al., 2003; Hernández-Blanco et al., 2007), lignin (Tronchet et al., 2010), hemicellulose (Vega-Sánchez et al., 2012), and pectin (Bethke et al., 2016), have been shown to result in altered plant response to pathogen attack.
The impact of modified cell walls on biotic interactions with non-pathogenic symbiotic microbes or microbiomes is less understood. Emerging evidence suggests that genetically modifying the lignin biosynthesis pathway in plants can significantly alter the composition of the associated microbiome (Beckers et al., 2016). The functional impacts of modified cellulose properties on beneficial plant-microbe inter actions have received far less attention. Here, we investigated the functional impact of down-regulating AtKOR/KORRIGAN-like cellulose biosynthesis pathway genes, PdKOR1 and PdKOR2, on cell wall composition, growth, metabolism and beneficial microbial interactions of Populus deltoides.
KORRIGAN (KOR), an endo-1,4-β-glucanase (EGase) belonging to the glycosyl hydrolase family 9 (GH9) (Henrissat, 1991; Urbanowicz et al., 2007), has been shown to be associated with the Cellulose Synthase Complex, CSC (Vain et al., 2014; Mansoori et al., 2014; Worden et al., 2015; Zhang et al., 2016), and to be integral to the pathways controlling content, crystallinity and degree of polymerization of cellulose in Arabidopsis (kor-1 and kor-2 and irx2) (Szyjanowicz et al., 2004; Takahashi et al., 2009). Modification of endo-1,4-β-glucanase has been shown to alter the defense response to Pseudomonas syringae (López-Cruz et al., 2014) in Arabidopsis. Previous studies in Populus have studied the phenotypic impact of simultaneous RNAi targeting of several KOR paralogs. Such an approach led to a reduction in cellulose content and secondary wall thickness with no apparent impacts on plant morphology in P. euramericana (Yu et al., 2014), and a contrasting phenotype of higher cellulose crystallinity, curled leaves and reduced height study in P. alba x grandidentata RNAi plants (Maloney and Mansfield, 2010).
Driven by the hypothesis that genetic modification of host cell wall can have cascading and quantifiable phenotypic impacts on its secondary metabolome and interaction with beneficial microbes, we investigated the consequence of a cellulose biosynthesis pathway gene modification on Populus metabolism and its ability to interact with a well-characterized fungal symbiont, Laccaria bicolor (Martin et al., 2008). Our results reveal significant phenotypic impacts, beyond the altered cellulose properties, to carbon partitioning (primary vs. secondary metabolites) and allocation (leaves vs. stem) pathways, and a distinctive secondary metabolome phenotype in PdKOR RNAi plants. Considering the biological, ecological, and economic significance of plant-microbe interactions, our results further emphasize the need for holistic considerations of plant trait modifications in crop improvement efforts.
Materials and Methods
Phylogeny and Protein Sequence Alignment
Protein sequences of the 26 Populus GH9 were collected from Phytozome v9.1 P. trichocarpa (Tuskan et al., 2006). The P. trichocarpa (Ptr) sequences, with significant sequence similarity to AtKOR, were labeled PtrKOR1-5. Correspondingly, P. deltoides (Pd) sequences were designated PdKOR1-5 and were used as a phylogenetic/sequence reference in the design of qRT-PCR primers and the analyses of P. deltoides libraries and constructs. Protein sequences from other species were collected from NCBI. A phylogenetic tree was developed based on the maximum likelihood method in the MEGA5 software program (Tamura et al., 2011). Bootstrap values were calculated from 500 independent bootstrap runs. Protein sequence alignment and percent sequence similarity/identity were generated using the GeneDoc program (Nicholas et al., 1997).
RNA extraction, using the Spectrum Plant Total RNA Kit (Sigma), and cDNA synthesis and qRT-PCR assays, were performed according to methods described previously (Payyavula et al., 2014). Primers are listed in Supplementary Table S1. The efficiency of these primers for use in P. deltoides was verified by evaluating PCR products from serial dilutions of P. deltoides cDNA libraries. Specificity was confirmed by performing a melt curve analysis on qRT-PCR reaction products and verifying the expected band size (amplicon sizes ranged from 60 to 240 bp) on an agarose gel.
Plant Transformation, Growth and Sampling
A ∼200-bp fragment from the gene-specific 3′UTR of the PdKOR1 and PdKOR2 genes from P. deltoides was cloned in the binary vector pAGSM552 in both sense and antisense orientations to form a hairpin loop with a chalcone synthase intron and expressed under the control of a Ubiquitin promoter. Agrobacterium-mediated plant transformation of the wild-type P. deltoides ‘WV94’ line was performed at the ArborGen LLC facility in Ridgeville, SC. The gene-specific RNAi constructs were each used to generate a minimum of 20 independent transformation events or lines. Five clonal replicates (i.e., ramets) for each PdKOR transgenic line, along with equal numbers of clonal replicates for empty vector transformed control plants, were propagated at Oak Ridge National Laboratory greenhouses and maintained at a constant day/night temperature of 25°C with 16 h photoperiod. In a preliminary study, the empty vector transformed control plants performed similar to the wild-type, non-transformed plants (data not shown). As a result, we used the averaged data from three independent empty-vector lines as the control in all subsequent experiments. All plants were initially propagated in small tubes (up to 50 cm in height) and then moved to 6-L pots. After ∼6 months of growth, plant height and stem diameter were measured and stem samples were collected and air dried for cell wall chemistry analyses. Initial phenotypic assessments included 14 and 9 independent RNAi lines for PdKOR1 and PdKOR2, respectively, and three replicated phenotyping cycles. In a second study, two RT-PCR-confirmed independent lines were propagated from greenwood cuttings to generate four to five ramets each. Estimates of leaf area were then calculated for five individual leaves between leaf plastochron index (LPI)-1 through LPI-5 for each ramet by multiplying midrib length by maximum width. At harvest time, the leaf at LPI-6, its petioles, phloem (i.e., excised bark), xylem (scraping from the stem) and remaining stem were collected. The phloem samples were frozen with dry ice. The stem samples were air dried and the remaining tissues were frozen in liquid nitrogen and stored at -80°C. Root weight was measured after thorough washing and drying at 70°C.
For co-culture studies, the ectomycorrhizal fungal strain, L. bicolor S238N, was established and maintained on Pachlewski’s medium (Pachlewski and Pachlewska, 1974). Mycelia were produced in a peat-vermiculite nutrient mixture, then placed in glass jars and grown in darkness for 2 months at 25°C., They were subsequently stored at 4°C just prior to host inoculation (Duponnois and Garbaye, 1991). Co-culture was performed in the greenhouse with eight replicate plants for each line, according to prior published methods (Labbe et al., 2011). Soil inoculation was carried out with a pure fungal culture mixed in peat-vermiculite at a concentration of 10% (1:9, v:v). Control (no inoculation) plants were established in autoclaved peat-vermiculite soil. The percentage of mycorrhizal colonization was determined under a stereomicroscope as previously described (Labbe et al., 2011). In brief, washed roots were cut into 1 cm segments, and for each root system, 100 apices were randomly examined and assessed as either mycorrhizal or non-mycorrhizal. Fresh weights of leaves, stems and roots were recorded for each of the eight replicates per line.
Non-structural Carbohydrate Analysis
Soluble sugars were estimated as previously described (Payyavula et al., 2012), with minor modifications. Briefly, freeze-dried leaf and air-dried stem samples (20 mg) were each extracted three times with 80% ethanol at 80°C. The supernatants were then combined and re-extracted with 50 mg of activated charcoal. Next, the pigment-free extract was dried overnight at 50°C to eliminate ethanol, resuspended in either 800 or 250 μl of water (leaf or stem, respectively) and then employed in sucrose and glucose estimation assays using commercial kits (Sigma, St. Louis, MO, USA). Starch was estimated from the plant sample fraction leftover after ethanol extract (pellet remaining after soluble sugar extraction). The pellet was treated with 720 μl of 0.1 N NaOH at 50°C for 30 min and was neutralized in 900 μl of 0.1 N acetic acid. Pellet starch was digested using 1 U each of α-amylase (from Aspergillus oryzae, Sigma, St. Louis, MO, USA) and amyloglucosidase (from A. niger, Sigma), as described elsewhere (Chow and Landhäusser, 2004; Payyavula et al., 2011). Starch quantities were calculated from the glucose standard calibration curve.
Cellulose Content, and Lignin Content and Syringyl-to-Guaiacyl Ratio
Cellulose content was estimated by the anthrone method, as previously described (Updegraff, 1969), with modifications for a 96-well plate format (Payyavula et al., 2014). Lignin and the syringyl-to-guaiacyl ratio (S/G ratio) were analyzed at the National Renewable Energy Laboratory in Golden, CO by pyrolysis molecular beam mass spectrometry (MBMS), as described elsewhere (Tuskan et al., 2006; Sykes et al., 2009). Briefly, 4 mg stem samples were pyrolyzed at 500°C in a quartz reactor using a Frontier py2020 auto-sampler. The pyrolysis vapors were then analyzed using a custom Extrel single-quadrupole molecular beam mass spectrometer.
Structural Carbohydrate Analysis
After sugar and starch extraction, the pellet was dried and used for structural carbohydrate analysis. The extract-free sample (5 mg) was hydrolyzed with 50 μl of H2SO4 (72%, v:v) at 37°C for 1 h then diluted with 1.45 ml of water and autoclaved for 60 min. After cooling, an aliquot was neutralized with CaCO3 and the carbohydrate composition was then analyzed as described elsewhere (Yee et al., 2012) on a high performance liquid chromatograph (LaChrom Elite® system, Hitachi High Technologies America, Inc.) equipped with a refractive index detector (model L-2490) and a UV–Vis detector (model L-2420). Individual sugars were quantified based on the standard curve developed for each compound.
Cellulose Crystallinity and Molecular Weight Measurements
Populus stem samples were ground and extracted in a dichloromethane solvent with a Foss Soxtec unit (SoxtecTM2050). Lignin and hemicellulose were removed from the extractive-free samples and used in a 13C cross-polarization magic angle spinning (CP/MAS) nuclear magnetic resonance (NMR) cellulose analysis by acetic acid holocellulose pulping and pre-acid hydrolysis, following published methodologies (Foston et al., 2011). We isolated α-cellulose from the holocellulose sample with an alkaline hemicellulose extraction procedure, as according to published literature (Foston et al., 2011), and then subjected it to a gel permeation chromatography (GPC) cellulose analysis.
In preparation for a 13C-CPMAS NMR analysis, 4 mm cylindrical ceramic MAS rotors were filled with isolated α-cellulose. Solid-state NMR measurements of 13C were carried out in a Bruker Avance-400 spectrometer operating at 100.55 MHz for with a Bruker double-resonance MAS probe spinning at 10 kHz. CP/MAS experiments utilized a 5 μs (90°) proton pulse and a 1.5 ms contact pulse with a 4 s recycle delay and 4–8 K scans. All spectra were recorded on equilibrated moisture samples (∼35% water content). Crystallinity was determined via a 2-peak integration of the C4 crystalline carbon region (δ 85–92 ppm) over the integral of the entire C4 region (δ 80–92 ppm), as published earlier (Wickholm et al., 1998; Foston et al., 2011).
The number-average molecular weight (Mn) and weight-average molecular weight (Mw) were determined by GPC after tricarbanilation of cellulose. Isolated cellulose from each sample was derivatized with phenyl isocyanate, following a published procedure (Wood et al., 1986; Foston et al., 2011). Prior to GPC analysis, derivatized cellulose was dissolved in tetrahydrofuran (THF, 1 mg ml-1), then filtered through a 0.45 mm filter and placed in a 2 ml auto-sampler vial. The molecular weight distributions of the cellulose tricarbanilate samples were analyzed by an Agilent GPC Security 1200 system equipped with four Waters Styragel columns (HR1, HR2, HR4, HR5), an Agilent refractive index detector and an Agilent UV detector (270 nm) using THF as the mobile phase with injection volumes of 20 ml. A calibration curve was constructed based on eight narrow polystyrene standards ranging in molecular weight from 1.5 × 103 to 3.6 × 106 g mol-1. Data collection and processing were performed using the Polymer Standards Service WinGPC Unity software (Build 6807). Molecular weights (Mn and Mw) were calculated using software relative to the universal polystyrene calibration curve. Number-average degree of polymerization (DPn) and weight-average degree of polymerization (DPw) were both obtained by dividing Mn and Mw by 519 g mol-1, the molecular weight of the tricarbanilated cellulose repeat unit. The polydispersity index was calculated by dividing Mw by Mn.
In preparation for Gas Chromatography (GC)-Mass Spectrometry (MS) metabolite profiling, approximately 50–75 mg (fresh weight) of mature leaf, phloem and xylem tissues were each extracted in 2.5 ml of 80% ethanol. The extracts were then combined in a 1.0 ml aliquot and dried in a nitrogen stream. Sorbitol was added before extraction and used as an internal standard. Dried extracts were dissolved in acetonitrile, followed by trimethylsilylation (TMS) and a GC-MS analysis was performed after 2 days, as described elsewhere (Jung et al., 2009; Li et al., 2012; Payyavula et al., 2014).
Metabolites of interest to be quantified were identified by database matching with the Wiley Registry 10th Edition combined with NIST 2014 mass spectral database and a large user-created database (>2300 spectra) of mass spectral EI fragmentation patterns of TMS-derivatized compounds. Peaks were integrated and analyzed using a key selected ion, a characteristic relatively unique m/z fragment for each peak, rather than the total ion chromatogram, to minimize integrating co-eluting metabolites. The extracted peaks of known metabolites were scaled back up to the total ion current using predetermined scaling factors. Unidentified metabolites used the scaling factor for the internal standard (sorbitol) and were denoted by their retention time (RT; min) as well as key mass-to-charge (m/z) ratios. Peaks were quantified by area integration and the concentrations were normalized to the quantity of the internal standard recovered, volume of sample extracted, derivatized, and injected, as described elsewhere (Tschaplinski et al., 2012). Fold-change values of metabolite levels in leaf, phloem and xylem tissues of RNAi plants were generated relative to control tissues. Student’s t-test was applied to identify significantly increased or decreased metabolites (p < 0.05; n ≥ 3).
For phenolic colorimetric analysis, 5 mg of freeze-dried leaf sample was extracted with 1.5 ml of 80% EtOH using a Tissue Lyzer (Qiagen, Valencia, CA, USA) for 5 min at 30 Hz. After centrifuging, the supernatant was analyzed for phenolics, chlorophyll and carotenoids, the pellet for tannins. Total phenolics were estimated in a 96-well plate by the Folin–Ciocalteu (FC) method, as previously described (Singleton and Rossi, 1965; Payyavula et al., 2014). Tannins from the pellet were estimated as described elsewhere (Porter et al., 1985). Specifically, 200 μl of 80% EtOH and 600 μl of butanol:HCl (95:5, v:v) containing 25 μl of 2% ferric-ammonium-sulfate (Sigma) in 2N HCl were combined with the pellet and heated at 95°C for 20 min. The sample’s absorbance was then recorded at A550-700 nm.
Phylogenetic and Expression Analyses of AtKOR-Like, GH9A Family Endo-β-1,4-Glucanase Genes in Populus
To understand the phylogenetic relationship and extent of sequence similarity among GH9 family members in Populus, we carried out sequence analyses as follows. Analysis of the amino acid sequences of all the annotated GH9 glycosyl hydrolase, endo -β-1,4-glucanase family members in the genomes of P. trichocarpa and Arabidopsis generated a tree with three distinct clades (Figure 1), designated, A, B, and C (Urbanowicz et al., 2007). Group A, which includes AtKOR, contains five Populus PdGH9A or PdKOR-like isoforms. At the amino acid level, PdKOR1 and PdKOR2 have a high degree of similarity with each other relative to the three other Populus paralogs in group A: PdKOR3, PdKOR4, and PdKOR5 (Supplementary Table S2). An analysis of the chromosomal positions showed that PdKOR1 and PdKOR2 lie within homeologous blocks (Supplementary Figure S1), which were created during the Salicoid whole-genome duplication event estimated to have occurred approximately 65 Mya (Tuskan et al., 2006; Kalluri et al., 2007).
FIGURE 1. Phylogenetic analyses of endo-β-1,4-glucanases from Populus, Arabidopsis and selected KOR-like genes from other species using MEGA5 program. The evolutionary history was inferred by using the Maximum Likelihood method based on the JTT matrix-based model. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (500 replicates) are shown next to the branches. The three distinct branches are represented by (A–C) subgroups. KORs from this study are in bold. The analysis involved 55 predicted protein sequences and the corresponding accession numbers are: At4g39000; At4g38990; At4g39010; At1g70710; At1g23210; At4g02290; At1g02800; At1g22880; At1g71380; At4g23560; At4g09740; At3g43860; At1g64390; At4g11050; At2g32990; At1g48930; At2g44570; At2g44560; At2g44540; At2g44550; At1g75680; At1g19940; At1g65610/AtKOR2; At5g49720/AtKOR/AtGH9A1; At4g24260; PaxgKOR (Populus alba x grandidentata): ADB82903.1; PtGH (P. tremuloides): AAS45400.1); PtrKOR1/PdKOR1 (P. trichocarpa, Potri.003G151700); PtrKOR2/PdKOR2 (Potri.001G078900); PtrKOR3/PdKOR3 (Potri.008G079500); PtrKOR4/PdKOR4 (Potri.010G177300); PtrKOR5/PdKOR5 (Potri.005G188500); PtrCel9A6 (Potri.005G237700); PtrGH9A7 (Potri.002G023900);PtrGH9A8 (Potri.005G115400); PtrGH9A9 (Potri.002G225200); PtrCel9B1/PopCel1 (Potri.001G083200); PtrGH9B3 (Potri.019G069300); PtrGH9B4 (Poptri.014G126900); PtrGH9B5 (Potri.002G202400); PtrGH9B6 (Potri.010G109200); PtrGH9B7 (Potri.008G132700); PtrGH9B8 (Potri.009G123900); PtrGH9B9 (Potri.004G162200); PtrGH9B10 (Potri.015G12800); PtrGH9B11 (Potri.014G157600); PtrGH9B12 (Potri.001G356000); PtrGH9C1 (Potri.007G071200); PtrGH9C2 (Potri.003G139600); PtrGH9C3 (Potri.001G092200); PtGH9_4 (Potri.001G098800); PtrGH9_16 (Potri.006G219700); PtrGH9_18 (Potri.015G127900); PttCel9A (P. tremula x tremuloides):AAT75041.1); SlGH (Solanum lycopersicum): AAC49704.1
Quantitative RT-PCR analyses showed that the relative expression level of PdKOR1 gene was highest in xylem among various Populus tissue types measured (Figure 2). The expression level of PdKOR2 gene were observed to be weaker across all tissues relative to PdKOR1 gene. Expression levels of PdKOR3 and PdKOR4 genes were at significantly lower, at detection limit, in several tissues, while that of PdKOR5 gene was below detection limit. An analysis of the cDNA sequences of the PdKOR1, PdKOR2, PdKOR3, PdKOR4, and PdKOR5 paralogs in Populus showed the non-coding 3′ UTR region to be least conserved and hence this region was targeted in the design of gene-specific RNAi constructs (Supplementary Figures S2 and S3).
FIGURE 2. Relative expression of KOR-like genes in various Populus tissue types. Relative expression (arbitrary units) was based on difference in critical threshold (CT) values between target gene and average of housekeeping genes, Ubiquitin-conjugating enzyme E2 and 18S ribosomal RNA. Expression of PdKOR5 (KOR5) was below quantifiable limits. YL, young leaf; ML, mature leaf; YS, young stem; MS, mature stem; PET, petiole; PHL, phloem; XYL, xylem; RT, root. Data represent means ± SE (n = 3).
Morphology and Growth of PdKOR Down-Regulated RNAi Plants
PdKOR1 and PdKOR2 transgenic RNAi plants displayed a highly similar phenotype of slower growth and differential biomass accumulation relative to control plants. Independent lines from both PdKOR1 and PdKOR2 RNAi plants had reduced plant height, stem diameter and number of internodes, and larger leaves relative to control plants (Figure 3). A significant impact on growth and biomass properties of PdKOR RNAi plants was also observed. RNAi plants displayed a significant reduction in stem height and diameter, and number of internodes (Figures 3A–C), and an increase in leaf area and petiole length (Figure 3D).
FIGURE 3. Phenotypic characterization of PdKOR1 RNAi, PdKOR2 RNAi, and control plants. (A) Image of a two-month old control and transgenic RNAi plants propagated from stem cuttings. Measurement of (B) Plant height and stem diameter, (C) Number of internodes, and (D) Leaf area and petiole length. Please refer to left-y-axis for plant height and leaf area, and to right-axis for stem girth and petiole length. Data represent mean ± SE (n = 4–5). ∗indicates statistically significant, p ≤ 0.05 based on Student’s t-tests.
We conducted RT-PCR analyses to study the specificity of the RNAi construct. RT-PCR assays showed a significant reduction in PdKOR1 as well as PdKOR2 transcript levels in either PdKOR1 or PdKOR2 RNAi background (Supplementary Figure S4A). The impact of RNAi down-regulation on transcript levels of PdKOR3 and PdKOR4 was found to be insignificant (Supplementary Figures S4B,C). PdKOR5 transcripts were found to be below detection limit as in the case of native tissue expression profiling.
Given the down-regulation of both closely related Populus KOR paralogs, PdKOR1 and PdKOR2, and similar phenotypes (Figure 3; Supplementary Figures S4A and S5), we henceforth refer to the transgenic RNAi lines as PdKOR lines and present results from representative lines of PdKOR1.
Down-Regulation of PdKOR Alters Cellulose, Cell Wall and Biomass Properties of Populus
In order to corroborate whether the sequence homologs of AtKOR have similar functional roles in cell wall biosynthesis in P. deltoides, we analyzed the cell walls of RNAi plants. KOR transgenic plants were characterized by reduced cellulose content (Figure 4A) and increased cellulose crystallinity (Figure 4B). Degree of polymerization (DP) measured by GPC technique yielded inconclusive results with a statistically significant reduction observed in only one of the PdKOR RNAi lines (Figure 4C). The effects of PdKOR down-regulation were measureable beyond cellulose to that of other cell wall properties. Statistically significant reductions in lignin content and S/G ratio were also evident (Figures 4D,E).
FIGURE 4. Cell wall composition of control and PdKOR1 downregulated lines. (A) percentage (%) cellulose content, (B) NMR-based cellulose crystallinity, (C) degree of polymerization of cellulose, (D) lignin content, and (E) S/G (syringyl-to-guaiacyl) ratio are represented as means ± SE (n = 3–5). (A–E) represent values from debarked stem samples. ∗indicates statistically significant, p ≤ 0.05 based on Student’s t-tests.
Down-Regulation of PdKOR Impacts Carbon Metabolism, Carbon Partitioning and Allocation
The differential biomass observed for RNAi leaf, stem, and root relative to the control shows that carbon allocation is impacted in transgenic PdKOR plants. In order to understand the impact of modification of KOR genes on carbon partitioning, we complemented cell wall polymer analysis with analyses of soluble non-structural carbohydrates and phenolics.
Our colorimetric analyses showed a relative decrease in the non-structural carbohydrate or sugar content in transgenic PdKOR RNAi plants. As such, the content of glucose, sucrose and starch in mature RNAi plant leaves was reduced by 30–80% relative to the control (Figure 5). In addition to confirming the colorimetric observation of reduced glucose content in RNAi plant tissues, MS analysis revealed a reduction in fructose, galactose and raffinose levels relative to the control.
FIGURE 5. Colorimetric estimation of non-structural carbohydrates in mature leaves. Data are represented as means ± SE (n = 4). ∗indicates statistically significant, p ≤ 0.05 based on Student’s t-tests.
We found that phenolic compounds were elevated in the RNAi plant leaves as well as phloem. Some of the predominantly phenolics were caffeic acid derivatives, including caffeoyl-glycosides, caffeoyl-shikimate conjugates, and other caffeoyl-conjugates (Supplementary Table S3). Two novel caffeoyl-glycoside conjugates, eluted at 14.14 and 19.31 min, were uniquely found in the PdKOR transgenic leaves, but not in the control. Syringin levels were reduced in nearly all of the lines. These changes in monolignol glucosides occur in parallel with an observed reduction of the S/G ratio in cell walls of RNAi plants. The observed reduction in carbon in the form of soluble sugars and increase in carbon represented as soluble phenolics suggests increased carbon partitioning to secondary metabolite pathways in RNAi plants (Table 1 and Supplementary Table S4).
The concentration of salicin, salicortin, and 6-hydroxy-2-cyclohexenone-1-carboxylic acid, a putative precursor of the latter, were elevated in all lines, coupled with a decline in 6-hydroxy-2-cyclohexenone (enol), 6-hydroxy-2-cyclohexenone alcohol and catechol, which are all putative breakdown products of salicortin. In the xylem of transgenic PdKOR lines, there was a significant decrease (>50%) in the concentrations of shikimic, maleic and threonic acid, and an increase (50–150%) in cis-aconitic acid (Supplementary Table S3).
Down-Regulation of PdKOR Impacts Mycorrhizal Association with Laccaria bicolor
In order to test whether the specific genetic modification of the host cell wall and its cascading effects on secondary metabolome has a quantifiable phenotypic impact on the ability to interact with microbes, we grew both transgenic RNAi and control plants with and without the addition of the known beneficial fungal symbiont, L. bicolor, and observed the mycorrhizal association.
Laccaria bicolor was found to have a higher colonization rate (i.e., mycorrhization rates) on the transgenic PdKOR plants relative to the control (Figure 6A). The greater mycorrhizal association also correlated with the significant gain in stem biomass in PdKOR RNAi lines relative to its non-inoculated plants following a two-month co-culture period in the greenhouse (Figure 6B). While an increase in root biomass was observed in inoculated vs. non-inoculated plants (Figure 6C), the gain in stem biomass in inoculated vs. non-inoculated plants was quantitative greater for PdKOR lines (74–90%) relative the control (61%) lines.
FIGURE 6. (A) Root colonization rate of the Populus transgenic PdKOR1 downregulated and control lines by Laccaria bicolor. (B) Stem biomass of the plants with or without L. bicolor inoculation. (C) Root biomass weight of the plants with or without L. bicolor inoculation. ∗indicates statistically significant, p ≤ 0.05 based on Student’s t-tests.
Delineating the genetic underpinnings and physiological implications of plant cell wall properties in bioenergy relevant crops such as Populus is critical to the goals of lignocellulosic bioenergy research. A cascading phenotypic effect of cellulose modification on carbon partitioning was distinctively supported by the metabolic profile of transgenic plants. The large diversity and fold-change ratios of higher-order caffeoyl-conjugates suggest that the reduced carbon flux to the lignin pathway resulted in a diversion of caffeic acid, via caffeoyl-quinates and caffeoyl-shikimates, to these higher-order caffeoyl-conjugates and may be symptomatic of dysfunctional cell wall assembly (Payyavula et al., 2014). The results suggest the phloem contains the large diversity of higher-order hydroxycinnamate-conjugates and is the organ that serves in both the storage and transfer of the carbon that was destined for lignin synthesis, but remained unassimilated. Several higher-order salicylates, including salicin and salicortin, are overrepresented in RNAi plants, owing to the diversion of carbon from structural lignin to these soluble phenolic microbial defense metabolites.
Previous studies using distinct species of Populus have showed that down-regulation of PaxgKOR (Maloney and Mansfield, 2010) and KOR-like gene from P. euramericana (Yu et al., 2014) resulted in modification in cellulose properties and a reduction of cell wall glucose content. Our PdKOR knockdown lines showed an increase in cellulose crystallinity, which is consistent with the results reported from down-regulation of PaxgKOR, (Maloney and Mansfield, 2010). However, this is contrary to the report of reduced crystallinity for Arabidopsis loss-of-function irx2 mutants and Arabidopsis plants overexpressing a Populus PttCel9A1 gene (Takahashi et al., 2009), as well as the report of unaffected crystallinity for Picea KOR (Maloney et al., 2012). Furthermore, the observed reduction in sucrose content of PdKOR RNAi leaves relative to the control is consistent with the results reported for PaxgKOR down-regulated plants (Maloney and Mansfield, 2010; Yu et al., 2014), however, the observed reduction in glucose levels is in contrast. It is important to note that the RNAi constructs in the two other studies are designed to target in conserved coding region, so as it is the case in this study both PdKOR1 and PdKOR2-like homologs were impacted. Taken together, these findings suggest that KOR’s role as a cellulose and cell wall modifying factor can differ qualitatively among species. The functional variations reported among studies may be attributable to differences in transgene technology applied, species-specific variations in protein isoforms and/or growth habit.
The enhanced mycorrhizal association, resulting in pronounced increases in co-culture induced root biomass, may be related to the differential metabolic and cell wall features of the PdKOR RNAi plants. The indirect effects on metabolome of RNAi plants potentially resulted in increased responsiveness, or priming of the root system, to its ectomycorrhizal partner (Tschaplinski et al., 2014). The altered lignin (Beckers et al., 2016), and shikimic (Tschaplinski et al., 2014), and phenolic and salicylic acid content (Lebeis et al., 2015), observed in PdKOR RNAi plants, have been previously reported as significant metabolites underlying plant-microbe signaling and interactions.
An alternate explanation for the enhanced mycorrhization may be drawn from the microbial biocontrol activity reported for cell wall degrading enzymes such as glucanases (Chet, 1987; Lam and Gaffney, 1993; Lorito et al., 1994; Lo, 1998; Harman, 2000). In this scenario, it may be speculated that PdKOR endo-β-1,4-glucanases have a direct functional role in defense and that down-regulation of the gene leads to a muted defense response phenotype, which promotes establishment of the plant-microbe relationship.
While there is evidence in literature for potential beneficial effects of mycorrhizae on plant growth and a suggestion of the importance of cell wall chemistry on mycorrhization (Goulao et al., 2011; Danielsen et al., 2013; Balestrini and Bonfante, 2014), a direct link to the beneficial effect of Laccaria on growth-compromised transgenic Populus plants and an underlying increase in mycorrhization rate has not been previously reported. It is plausible that increased nutrient acquisition, potentially nitrogen and phosphorus as a resulting from enhanced mycorrhization (Szuba, 2015) contributed to improved carbon assimilation and measurable growth benefits in transgenic Populus RNAi plants. These results provide an additional strategy to improve growth in transgenic plants that would otherwise produce below average biomass. The precise structural, molecular, and biochemical bases of improved mycorrhization rates, as well as the significance of the mutualistic relationship as a means of improving transgenic plant performance, merit further studies.
The present study examined and validated the hypothesis that genetic modification of the host cell wall can have cascading and quantifiable phenotypic impacts on its secondary metabolome and its interaction with symbiotic microbes. It remains to be examined whether modifying endo-β-1,4-glucanase or a distinct cell wall pathway gene has an impact on functional microbiome associated with plants under field settings. It also remains to be clarified whether the pleiotropic phenotype is the result of a direct or an indirect, feed-back or feed-forward, cascade of changes triggered by the gene modification. Evaluation of additional plant mutants, genotypes or species, and plant metabolites for their interactions with microbes or microbiome are needed to delineate the specific structural and secondary metabolome characteristics influencing plant-microbe interactions.
UK, JL, and RP conceived the work, conducted analyses and wrote the manuscript. GT and TT conducted analyses and wrote the manuscript. RS, GB, TT, NE, SJ, AR, and MD, conducted analyses. UK and RP contributed equally to this work.
This research was supported by funding from U.S. Department of Energy, Office of Science, Biological and Environmental Research to the BioEnergy Science Center (BESC) and Plant Microbe Interfaces (PMI) projects. BESC is a U.S. DOE Energy Bioenergy Research Center supported by the Office of Biological and Environmental Research in the DOE Office of Science. Oak Ridge National Laboratory is managed by UT-Battelle, LLC for the U.S. Department of Energy under Contract Number DE–AC05–00OR22725. Oak Ridge National Laboratory is managed by UT-Battelle, LLC for the DOE under Contract Number DE–AC05–00OR22725.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We are grateful to Arborgen LLC for undertaking Populus transformation and tissue culture, Miguel Rodriguez Jr. for HPLC validation of processed cell wall samples. ORNL greenhouse staff for plant care and maintenance, and Daniel Smallwood for assistance with formatting the manuscript.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.01455
Beckers, B., Op De Beeck, M., Weyens, N., Van Acker, R., Van Montagu, M., Boerjan, W., et al. (2016). Lignin engineering in field-grown poplar trees affects the endosphere bacterial microbiome. Proc. Natl. Acad. Sci. U.S.A. 113, 2312–2317. doi: 10.1073/pnas.1523264113
Bethke, G., Thao, A., Xiong, G., Li, B., Soltis, N. E., Hatsugai, N., et al. (2016). Pectin biosynthesis is critical for cell wall integrity and immunity in Arabidopsis thaliana. Plant Cell 2016, 537–556. doi: 10.1105/tpc.15.00404
Caño-Delgado, A., Penfield, S., Smith, C., Catley, M., and Bevan, M. (2003). Reduced cellulose synthesis invokes lignification and defense responses in Arabidopsis thaliana. Plant J. 34, 351–362. doi: 10.1046/j.1365-313X.2003.01729.x
Cantu, D., Vicente, A. R., Labavitch, J. M., Bennett, A. B., and Powell, A. L. (2008). Strangers in the matrix: plant cell walls and pathogen susceptibility. Trends Plant Sci. 13, 610–617. doi: 10.1016/j.tplants.2008.09.002
Chen, Z., Hong, X., Zhang, H., Wang, Y., Li, X., Zhu, J. K., et al. (2005). Disruption of the cellulose synthase gene, AtCesA8/IRX1, enhances drought and osmotic stress tolerance in Arabidopsis. Plant J. 43, 273–283.
Chet, I. (1987). “Trichoderma–application, mode of action, and potential as a biocontrol agent of soilborne plant pathogenic fungi,” in Innovative Approaches to Plant Disease Control, ed. I. Chet (New York, NY: John Wiley & Sons), 137–160.
Danielsen, L., Lohaus, G., Sirrenberg, A., Karlovsky, P., Bastien, C., Pilate, G., et al. (2013). Ectomycorrhizal colonization and diversity in relation to tree biomass and nutrition in a plantation of transgenic poplars with modified lignin biosynthesis. PLoS ONE 8:e59207. doi: 10.1371/journal.pone.0059207
Duponnois, R., and Garbaye, J. (1991). Effect of dual inoculation of Douglas fir with the ectomycorrhizal fungus Laccaria laccata and mycorrhization helper bacteria (MHB) in two bare-root forest nurseries. Plant Soil 138, 169–176. doi: 10.1007/BF00012243
Goubet, F., Barton, C. J., Mortimer, J. C., Yu, X., Zhang, Z., Miles, G. P., et al. (2009). Cell wall glucomannan in Arabidopsis is synthesised by CSLA glycosyltransferases, and influences the progression of embryogenesis. Plant J. 60, 527–538. doi: 10.1111/j.1365-313X.2009.03977.x
Goulao, L. F., Vieira-Silva, S., and Jackson, P. A. (2011). Association of hemicellulose- and pectin-modifying gene expression with Eucalyptus globulus secondary growth. Plant Physiol. Biochem. 49, 873–881. doi: 10.1016/j.plaphy.2011.02.020
Hernández-Blanco, C., Feng, D. X., Hu, J., Sánchez-Vallet, A., Deslandes, L., Llorente, F., et al. (2007). Impairment of cellulose synthases required for Arabidopsis secondary cell wall formation enhances disease resistance. Plant Cell 19, 890–903. doi: 10.1105/tpc.106.048058
Labbe, J., Jorge, V., Kohler, A., Vion, P., Marcais, B. T., Bastien, C., et al. (2011). Identification of quantitative trait loci affecting ectomycorrhizal symbiosis in an interspecific F1 poplar cross and differential expression of genes in ectomycorrhizas of the two parents: Populus deltoides and Populus trichocarpa. Tree Genetics & Genomes 7, 617–627. doi: 10.1007/s11295-010-0361-3
Lebeis, S. L., Paredes, S. H., Lundberg, D. S., Breakfield, N., Gehring, J., McDonald, M., et al. (2015). PLANT MICROBIOME. Salicylic acid modulates colonization of the root microbiome by specific bacterial taxa. Science 349, 860–864. doi: 10.1126/science.aaa8764
Li, Y., Tschaplinski, T., Engle, N., Hamilton, C., Rodriguez, M. Jr., Liao, J. C., et al. (2012). Combined inactivation of the Clostridium cellulolyticum lactate and malate dehydrogenase genes substantially increases ethanol yield from cellulose and switchgrass fermentations. Biotechnol. Biofuels 5:2. doi: 10.1186/1754-6834-5-2
Liang, Y. K., Xie, X., Lindsay, S. E., Wang, Y. B., Masle, J., Williamson, L., et al. (2010). Cell wall composition contributes to the control of transpiration efficiency in Arabidopsis thaliana. Plant J. 64, 679–686. doi: 10.1111/j.1365-313X.2010.04362.x
López-Cruz, J., Finiti, I., Fernández-Crespo, E., Crespo-Salvador, O., García-Agustín, P., and González-Bosch, C. (2014). Absence of endo-1,4-(-glucanase KOR1 alters the jasmonate-dependent defence response to Pseudomonas syringae in Arabidopsis. J. Plant Physiol. 171, 1524–1532. doi: 10.1016/j.jplph.2014.07.006
Lorito, M., Hayes, C. K., Di Pietro, A., Woo, S. L., and Harman, G. E. (1994). Purification, characterization, and synergistic activity of a glucan 1,3-glucosidase and an N-acetyl-beta-glucosaminidase from Trichoderma harziarnium. Phytopathology 84, 398–405. doi: 10.1094/Phyto-84-398
Maloney, V. J., and Mansfield, S. D. (2010). Characterization and varied expression of a membrane-bound endo-β-1,4-glucanase in hybrid poplar. Plant Biotechnol. J. 8, 294–307. doi: 10.1111/j.1467-7652.2009.00483.x
Maloney, V. J., Samuels, A. L., and Mansfield, S. D. (2012). The endo-1,4-(-glucanase Korrigan exhibits functional conservation between gymnosperms and angiosperms and is required for proper cell wall formation in gymnosperms. New Phytol. 193, 1076–1087. doi: 10.1111/j.1469-8137.2011.03998.x
Mansoori, N., Timmers, J., Desprez, T., Alvim-Kamei, C. L., Dees, D. C. T., Vincken, J.-P., et al. (2014). KORRIGAN1 interacts specifically with integral components of the cellulose synthase machinery. PLoS ONE 9:e112387. doi: 10.1371/journal.pone.0112387
Martin, F., Aerts, A., Ahrén, D., Brun, A., Danchin, E. G., Duchaussoy, F., et al. (2008). The genome of Laccaria bicolor provides insights into mycorrhizal symbiosis. Nature 452, 88–92. doi: 10.1038/nature06556
Pachlewski, R., and Pachlewska, J. (1974). Studies on Symbiotic Properties of Mycorrhizal Fungi of Pine (Pinus silvestris L.) with the Aid of the Method of Mycorrhizal Synthesis in Pure Cultures on Agar. Warsaw: Forest Research Institue.
Park, J. I., Ishimizu, T., Suwabe, K., Sudo, K., Masuko, H., Hakozaki, H., et al. (2010). UDP-glucose pyrophosphorylase is rate limiting in vegetative and reproductive phases in Arabidopsis thaliana. Plant Cell Physiol. 51, 981–996. doi: 10.1093/pcp/pcq057
Payyavula, R. S., Navarre, D., Kuhl, J., Pantoja, A., and Pillai, S. (2012). Differential effects of environment on potato phenylpropanoid and carotenoid expression. BMC Plant Biol. 12:39. doi: 10.1186/1471-2229-12-39
Payyavula, R. S., Tay, K. H. C., Tsai, C.-J., and Harding, S. A. (2011). The sucrose transporter family in Populus: the importance of a tonoplast PtaSUT4 to biomass and carbon partitioning. Plant J. 65, 757–770. doi: 10.1111/j.1365-313X.2010.04463.x
Payyavula, R. S., Tschaplinski, T. J., Jawdy, S. S., Sykes, R. W., Tuskan, G. A., and Kalluri, U. C. (2014). Metabolic profiling reveals altered sugar and secondary metabolism in response to UGPase overexpression in Populus. BMC Plant Biol. 14, doi: 10.1186/s12870-014-0265-8
Sykes, R., Yung, M., Novaes, E., Kirst, M., Peter, G., and Davis, M. (2009). High-throughput screening of plant cell-wall composition using pyrolysis molecular beam mass spectroscopy. Methods Mol. Biol. 581, 169–183.
Szyjanowicz, P. M. J., McKinnon, I., Taylor, N. G., Gardiner, J., Jarvis, M. C., and Turner, S. R. (2004). The irregular xylem 2 mutant is an allele of korrigan that affects the secondary cell wall of Arabidopsis thaliana. Plant J. 37, 730–740. doi: 10.1111/j.1365-313X.2003.02000.x
Takahashi, J., Rudsander, U. J., Hedenström, M., Banasiak, A., Harholt, J., Amelot, N., et al. (2009). KORRIGAN1 and its aspen homolog PttCel9A1 decrease cellulose crystallinity in Arabidopsis stems. Plant Cell Physiol. 50, 1099–1115. doi: 10.1093/pcp/pcp062
Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M., and Kumar, S. (2011). MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 28, 2731–2739. doi: 10.1093/molbev/msr121
Thévenin, J., Pollet, B., Letarnec, B., Saulnier, L., Gissot, L., Maia-Grondard, A., et al. (2011). The simultaneous repression of CCR and CAD, two enzymes of the lignin biosynthetic pathway, results in sterility and dwarfism in Arabidopsis thaliana. Mol. Plant 4, 70–82. doi: 10.1093/mp/ssq045
Tronchet, M., Balagué, C., Kroj, T., Jouanin, L., and Roby, D. (2010). Cinnamyl alcohol dehydrogenases-C and D, key enzymes in lignin biosynthesis, play an essential role in disease resistance in Arabidopsis. Mol. Plant Pathol. 11, 83–92. doi: 10.1111/j.1364-3703.2009.00578.x
Tschaplinski, T. J., Plett, J. M., Engle, N. L., Deveau, A., Cushman, K. C., Martin, M. Z., et al. (2014). Populus trichocarpa and Populus deltoides exhibit different metabolomic responses to colonization by the symbiotic fungus Laccaria bicolor. Mol. Plant Microbe Interact. 27, 546–556. doi: 10.1094/MPMI-09-13-0286-R
Tschaplinski, T. J., Standaert, R. F., Engle, N. L., Martin, M. Z., Sangha, A. K., Parks, J. M., et al. (2012). Down-regulation of the caffeic acid O-methyltransferase gene in switchgrass reveals a novel monolignol analog. Biotechnol. Biofuels 5:71. doi: 10.1186/1754-6834-5-71
Urbanowicz, B. R., Bennett, A. B., del Campillo, E., Catalá, C., Hayashi, T., Henrissat, B., et al. (2007). Structural organization and a standardized nomenclature for plant endo-1,4-(-glucanases (cellulases) of glycosyl hydrolase family 9. Plant Physiol. 144, 1693–1696. doi: 10.1104/pp.107.102574
Vain, T., Crowell, E. F., Timpano, H., Biot, E., Desprez, T., Mansoori, N., et al. (2014). The cellulase KORRIGAN is part of the cellulose synthase complex. Plant Physiol. 165, 1521–1532. doi: 10.1104/pp.114.241216
Vega-Sánchez, M. E., Verhertbruggen, Y., Christensen, U., Chen, X., Sharma, V., Varanasi, P., et al. (2012). Loss of Cellulose synthase-like F6 function affects mixed-linkage glucan deposition, cell wall mechanical properties, and defense responses in vegetative tissues of rice. Plant Physiol. 159, 56–69. doi: 10.1104/pp.112.195495
Wickholm, K., Larsson, P. T., and Iversen, T. (1998). Assignment of non-crystalline forms in cellulose I by CP/MAS 13C NMR spectroscopy. Carbohydr. Res. 312, 123–129. doi: 10.1016/S0008-6215(98)00236-5
Wood, B. F., Conner, A. H., and Hill, C. G. (1986). The effect of precipitation on the molecular weight distribution of cellulose tricarbanilate. J. Appl. Polym. Sci. 32, 3703–3712. doi: 10.1002/app.1986.070320225
Worden, N., Wilkop, T. E., Esteve, V. E., Jeannotte, R., Lathe, R., Vernhettes, S., et al. (2015). CESA TRAFFICKING INHIBITOR inhibits cellulose deposition, and interferes with the trafficking of cellulose synthase complexes, and their associated proteins KORRIGAN1, and POM2/CELLULOSE SYNTHASE INTERACTIVE PROTEIN1. Plant Physiol. 167, 381–393. doi: 10.1104/pp.114.249003
Yang, C., Xu, Z., Song, J., Conner, K., Vizcay Barrena, G., and Wilson, Z. A. (2007). Arabidopsis MYB26/MALE STERILE35 regulates secondary thickening in the endothecium and is essential for anther dehiscence. Plant Cell 19, 534–548. doi: 10.1105/tpc.106.046391
Yee, K. L., Rodriguez, M. Jr., Tschaplinski, T. J., Engle, N. L., Martin, M. Z., Fu, C., et al. (2012). Evaluation of the bioconversion of genetically modified switchgrass using simultaneous saccharification and fermentation and a consolidated bioprocessing approach. Biotechnol. Biofuels 5:81. doi: 10.1186/1754-6834-5-81
Zhang, Y., Nikolovski, N., Sorieul, M., Vellosillo, T., McFarlane, H. E., Dupree, R., et al. (2016). Golgi-localized STELLO proteins regulate the assembly and trafficking of cellulose synthase complexes in Arabidopsis. Nat Commun. 9:11656. doi: 10.1038/ncomms11656
Keywords: plant cell wall, endo-β-1, 4-glucanase, carbon partitioning and allocation, cellulose, biomass, mycorrhiza, microbe, Populus
Citation: Kalluri UC, Payyavula RS, Labbé JL, Engle N, Bali G, Jawdy SS, Sykes RW, Davis M, Ragauskas A, Tuskan GA and Tschaplinski TJ (2016) Down-Regulation of KORRIGAN-Like Endo-β-1,4-Glucanase Genes Impacts Carbon Partitioning, Mycorrhizal Colonization and Biomass Production in Populus. Front. Plant Sci. 7:1455. doi: 10.3389/fpls.2016.01455
Received: 06 June 2016; Accepted: 12 September 2016;
Published: 04 October 2016.
Edited by:Chang-Jun Liu, Brookhaven National Laboratory, USA
Reviewed by:Henrik Scheller, Lawrence Berkeley National Laboratory, USA
Charles T. Anderson, Pennsylvania State University, USA
Copyright © 2016 Kalluri, Payyavula, Labbé, Engle, Bali, Jawdy, Sykes, Davis, Ragauskas, Tuskan and Tschaplinski. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Udaya C. Kalluri, firstname.lastname@example.org
†These authors have contributed equally to this work.