Perturbation ofBrachypodium distachyon CELLULOSE SYNTHASE A4or7results in abnormal cell walls
© Handakumbura et al.; licensee BioMed Central Ltd. 2013
Received: 20 March 2013
Accepted: 5 September 2013
Published: 11 September 2013
Cellulose is an integral component of the plant cell wall and accounts for approximately forty percent of total plant biomass but understanding its mechanism of synthesis remains elusive. CELLULOSE SYNTHASE A (CESA) proteins function as catalytic subunits of a rosette-shaped complex that synthesizes cellulose at the plasma membrane. Arabidopsis thaliana and rice (Oryza sativa) secondary wall CESA loss-of-function mutants have weak stems and irregular or thin cell walls.
Here, we identify candidates for secondary wall CESAs in Brachypodium distachyon as having similar amino acid sequence and expression to those characterized in A. thaliana, namely CESA4/7/8. To functionally characterize BdCESA4 and BdCESA7, we generated loss-of-function mutants using artificial microRNA constructs, specifically targeting each gene driven by a maize (Zea mays) ubiquitin promoter. Presence of the transgenes reduced BdCESA4 and BdCESA7 transcript abundance, as well as stem area, cell wall thickness of xylem and fibers, and the amount of crystalline cellulose in the cell wall.
These results suggest BdCESA4 and BdCESA7 play a key role in B. distachyon secondary cell wall biosynthesis.
KeywordsArtificial microRNA Brachypodium distachyon Crystallinity Secondary cell wall Sum-frequency-generation vibration spectroscopy Xylem
With continued consumption of fossil fuels, humankind faces a growing challenge of finding renewable sources of energy. Although plant derived biomass contains appreciable energy, to serve as a fuel for transportation it must be chemically or biologically liquefied, a conversion that is the subject of strenuous efforts to make economical. Success will presumably require not only advances in chemistry such as better catalysts, but also improved plants to serve as feedstock. As an ideal input, grasses are receiving considerable attention because certain species grow to great density, are perennial, and can require little if any fertilizer or irrigation . But their very size, longevity, and genome complexity makes these species difficult subjects to study and breed.
For studying grasses, whether as sources for biofuel or for any other grass-specific question, an emerging model is Brachypodium distachyon. This species is closely related to cereals and temperate grasses, and is about the same size as Arabidopsis thaliana and has about the same generation time. Studies on B. distachyon benefit from rapidly developing community resources, such as a completely sequenced genome, mutant collections (both chemically derived and sequence-indexed DNA insertions), and efficient protocols for transformation and crossing [2–5]. Of particular note for biofuels research, B. distachyon has cell walls that are similar compositionally to that of other grasses, like wheat (Tritium aestavum), barley (Hordeum vulgare), and Miscanthus[6–8]. Insofar as the cell wall constitutes almost the entirety of the input for converting biomass to biofuel, this similarity, along with the genetic attributes of this small grass, emphasize its suitability as a model for grass-related, biomass crop research.
Within the cell wall, as a target for optimization, cellulose is pre-eminent. Cellulose is the most abundant of any single wall component and is made exclusively of glucose, a tractable and energy rich molecule. Cellulose comprises long polymers of (1–4) β-linked glucose that are synthesized at the plasma membrane and associate laterally into a microfibril. Because of the configuration of the glucose residues, hydrogen bonds form at great density both within and between chains, a density that gives cellulose an elastic modulus rivaling that of steel but makes the structure impervious to degradation whether chemical or enzymatic.
Within the plasma membrane, the structure synthesizing cellulose is called a “terminal complex” [9–11]. In land plants and related green algae, the terminal complex as seen in the electron microscope comprises six subunits with hexagonal symmetry and is called a “rosette” . The primary constituents of the rosette are CELLULOSE SYNTHASE A (CESA) proteins. These proteins belong to processive glycosyltransferase family 2 and are thought to be the catalytic subunits for polymerizing the glucose chain. In angiosperms, CESAs usually comprise a small gene family with around ten members [13, 14].
Identification of CESA proteins and characterization of their function has greatly benefited from the facile genetics of A. thaliana. From this work, it emerged that certain CESA proteins synthesize the primary cell wall whereas others synthesize the secondary cell wall [15, 16]. Furthermore, it appeared that a given cell must express three distinct CESA proteins to produce cellulose at optimal levels. For the secondary cell wall, a screen based on collapsed xylem cells led to the identification of several irregular xylem (irx) lines, three of which, irx5, irx3, and irx1, harbor lesions in AtCESA4, 7, and 8, respectively [17–19]. Supporting the hypothesis of non-redundancy, these genes are expressed at similar levels in similar cell types, and the null mutants have indistinguishable phenotypes, including weak stems, collapsed xylem, and thin secondary cell walls that are deficient in cellulose.
Identification of a trio of CESA genes primarily responsible for synthesizing cellulose in secondary cell walls has been supported by work in other systems, including grasses. First, AtCESA4, 7, and 8 are represented in some angiosperm species by a single sequence each, and orthologs are more closely related than homologs (i.e., CESA4s of various species resemble each other more closely than do CESA4, 7, and 8 of a single species) [13, 20, 21]. Second, in grasses, a mutant of barley, fragile stem 2, with brittle stems and low cellulose content in the mature plant was attributed to a lesion in the barley ortholog of AtCESA4 . In addition, “brittle culm” mutants in rice (Oryza sativa) have been mapped to the three orthologs of AtCEA4, 7 and 8 and again these mutants have similar phenotypes, including modest dwarfism, thinner and weaker culms, and reduced cellulose content . However, it is not understood why three distinct proteins are needed nor is it known which of the loss-of-function phenotypes result directly from the lost protein and which result as a consequence of cumulative effects.
Here we describe the CESA gene family in B. distachyon with the aid of gene expression profiling and phylogeny. Furthermore a detailed analysis of the candidate secondary CESAs were performed by functionally characterizing mutants generated with specific artificial microRNA constructs.
Brachypodium distachyonCESA gene family
In view of our interest in the secondary cell wall, we examined the B. distachyon secondary cell wall sequences at greater depth (Additional file 2: Figure S2). In the cystine rich RING-type zinc finger domain, the sequences differ in the spacing between the second and third cystine. BdCESA4 has the canonical 15 amino acids whereas BdCESA7 has a single amino acid insertion and BdCESA8 has an eight amino acid deletion. Such spacing variations have apparently not been previously reported for A. thaliana secondary CESAs. However, the rice secondary CESA homologous to BdCESA8 exhibits the same eight amino acid deletion . The RING-type zinc finger domain is the distinguishing feature of CESA proteins seen in the first portion of the N-terminus and is thought to be involved in CESA protein dimerization [10, 11].
Brachypodium distachyonsecondary cell wall CESA gene expression
Localization of putative secondary cell wall CESA transcripts
Artificial microRNAs targeting BdCESA4 and BdCESA7
BdCESA4 and BdCESA7 knock-down lines and the structure of the stem
To examine the effects of the artificial microRNA constructs on cell wall structure, we measured cell wall thickness in the toluidine blue-stained sections. For both targets, the constructs reduced the thickness of cell walls modestly, but significantly (Figure 7D-G). The reduction was similar for xylem as well as for interfascicular fibers. This observation is consistent with the phenotypes of secondary CESA mutants characterized in other grass species [22, 23, 25].
Knock-down of BdCESA4 and BdCESA7 and crystalline cellulose content
Second, we evaluated the amount of crystalline cellulose by sum-frequency-generation (SFG) vibration spectroscopy . In this method, the sample is irradiated simultaneously by visible laser pulses at 532 nm and by infrared laser pulses with a tunable frequency. Among reflected and scattered lights, there are photons whose frequency is the sum of two input laser frequencies, which can be filtered and recorded separately. Due to the symmetry requirements, this frequency summation can be caused by crystalline cellulose but not by amorphous non-crystalline cell wall components . For this analysis, lines were grown at the same time. In the spectra, there are three prominent peaks indicative of crystalline cellulose Iβ: namely 2850 cm-1 ascribed to symmetric CH2 stretching, 2944 cm-1 ascribed to asymmetric CH2 stretching, and 3320 cm-1 ascribed to the intra-chain hydrogen-bonded hydroxyl stretch (Figure 8C) [26, 27]. Taking Avicel, a model cellulose Iβ with a high crystallinity, as a standard, the intensity at 2944 cm-1 can be used to estimate the crystalline cellulose amount by means of the previously determined calibration curve . The comparison of the intensities recorded at 2944 cm-1 indicated that the Avicel-equivalent crystalline cellulose content tended to be reduced in amiR-CESA7 and was significantly reduced in amiR-CESA4 (Figure 8D). These results are comparable to those from X-ray diffraction.
For B. distachyon, this is apparently the first study detailing the CESA gene family and functionally characterizing BdCESAs involved in secondary cell wall synthesis. Members of the CESA family are best characterized in A. thaliana, where they are assigned a role for cellulose synthesis in either primary or secondary cell wall. Other vascular plants have a similar gene family structure . The B. distachyon CESA gene family includes the uncharacterized grass-specific CESA clade (BdCESA10), which has been identified in all grass genomes sequenced to date, and all the other previously described clades [13, 24, 28].
Using amino acid sequences, the ten BdCESA genes were categorized into two groups: primary and secondary. Each of the three secondary CESA clades contains a single A. thaliana protein and a single rice protein, which have all been functionally characterized [17, 18, 23]. Fittingly, a single BdCESA protein was present in each of these clades, suggesting they have not expanded since the time of eudicot and monocot divergence 140–150 million years ago . This differs for the primary CESAs, which have differentially expanded between eudicots and monocots. There are different numbers of proteins in the CESA1 and CESA3 clades in A. thaliana, B. distachyon, and rice, and the CESA6 clade has diverged into separate eudicot and monocot clades, referred to as CESA6A and CESA6B, respectively . Among the three secondary B. distachyon CESAs, BdCESA8 was the most divergent, as was that clade member (OsCESA4) in rice . Moreover, BdCESA8 and OsCESA4 both lack eight amino acids in the first RING type zinc finger motif. However, A. thaliana CESA8 neither lacks those eight amino acids nor is notably divergent.
Confirming the deduction from phylogeny, expression of the secondary BdCESAs was generally enriched in stems, which are abundant in secondary cell walls, and transcripts were specifically abundant in stem vascular tissue and the surrounding mechanical tissue including sub-epidermal cell layers, tissues that all make secondary cell walls. This is consistent with the tissue-specific patterns of secondary CESA expression observed in rice, maize, A. thaliana, and barley [19, 24, 30, 31]. Interestingly, in developing stems of barley, HvCESA8 expression was two-fold greater than that of HvCESA4, which was also observed here for BdCESA8 and BdCESA4. Taken together, we conclude that BdCESA4, 7, and 8 encode cellulose synthase catalytic subunits for the secondary cell wall.
To our knowledge, this is the first report of artificial microRNA-mediated gene silencing in B. distachyon. We characterized the loss-of-function lines generated for BdCESA4 or BdCESA7 for changes in morphology and cellulose crystallinity. While these lines had a small but significant delay in inflorescence emergence and a reduction in stature, whole plant morphology was similar to wild-type plants. Similarly, delayed flowering and reduced stature was observed in the A. thaliana irx3 mutant and in rice brittle culm mutants [25, 32]. Sections of the stem revealed that the knock-down lines had a small but significant reduction in stem internode transverse cross-sectional area, a reduction that resembles that observed in rice brittle culm mutants and barley fragile stem 2[23, 33].
The loss-of-function lines for BdCESA4 and 7 had cell walls in xylem and sclerenchyma that were modestly but significantly thinner than wild-type. Thinner cell walls are expected from a defect in secondary cell wall cellulose synthesis, and have been observed consistently in secondary cell wall cesa mutants [22, 23, 25]. In A. thaliana, the defining phenotype of secondary CESA mutants is irregular or collapsed xylem [17, 18]. However, xylem contours appeared to be regular here, as they do in the secondary CESA mutants of rice . These differences probably reflect distinctions between eudicots and monocots in xylem tension or cell wall composition rather than differences in protein function [34, 35].
Cellulose is a linear glucan polymer made of β (1–4) linked glucose molecules. Thirty-six such glucan chains are predicted to form a cellulose microfibril . Hydrogen bonding between microfibrils crystallizes multiple cellulose chains together providing the physical properties necessary to maintain rigid plant cell walls. Crystallinity also allows cellulose to be readily quantified. We took advantage of spectroscopic techniques (X-ray diffraction and sum frequency generation) and polarized-light microscopy to assess the status of cellulose in the loss-of-function lines.
For the amiR-CESA7 mutants there was no significant decrease in crystallinity seen spectroscopically and the small decrease in birefringent retardance might reflect the thinner cell walls. Given that BdCESA7 expression was reduced by less than a factor of two, a generally wild-type cellulose status is perhaps not unexpected. In contrast, the amiR-CESA4 mutants had substantially decreased retardance and significantly decreased crystallinity from both diffraction and spectroscopic methods. The most parsimonious explanation for these data is that the cell walls of this line simply contained less cellulose than those of the wild-type. A nuanced explanation is that the cellulose synthesized in this line contained more defects, such as amorphous regions, or comprised microfibrils with fewer than the usual number (e.g., 36) of glucose chains. These alternatives are difficult to distinguish. Be that as it may, the significant decrease in crystalline cellulose amount or quality in this line matches the nearly 10-fold reduction in BdCESA4 expression. This result implies that the function of BdCESA4 is at least partially non-redundant and it is noteworthy that the expression of neither BdCESA7 nor 8 was elevated in an attempt at compensation.
Interestingly, even though the two lines differed in the strength of the knock-down and in the consequent loss of cellulose, the morphological defects in the plants were similar and modest. This implies that small decreases in CESA activity are sufficient to impact the plant, possibly through a feedback system analogous to the one A. thaliana invoking THESEUS , but that progressively larger changes to the plant need not follow from larger decreases in CESA activity. In this connection, it will be interesting to observe the phenotype of a complete loss of function mutant for BdCESA4.
Plant material and growth
Brachypodium distachyon (L.) line Bd 21–3 was used throughout. Seeds were imbibed in moist paper towels for seven days at 6°C, planted on potting mix (#2; Conrad Fafard Inc. Agawam, MA), and grown in a growth chamber at 20°C with 20 h light: 4 h dark, at a fluence rate of 220 μmol.m-2.s-1 and relative humidity of ~68%. For plate-grown plants, seeds were de-hulled and imbibed in water for 2 h with shaking. Then, seeds were treated with 70% ethanol for 20 s, rinsed with sterile water, and soaked in 1.3% NaClO for 4 min at room temperature while shaking. Seeds were subsequently rinsed three times with sterile water and stored in the dark at 4°C for a minimum of 2 days in a sterile Petri dish with filter paper. Seedlings were grown for seven days on 0.5 strength MS medium adjusted to a pH of 5.8 with KOH and containing 0.7% agar (Difco “Bacto agar”).
Identification of CELLULOSE SYNTHASE Agenes
Complete, translated amino acid sequences of the ten A. thaliana CESA proteins and members of the cellulose synthase-like families were used as queries in BLASTP of Phytozome v8.0 and NCBI databases to identify homologous B. distachyon CESA genes and identified proteins were named according to the A. thaliana counterparts, where possible. Multiple sequence alignments were performed using ClustalW program and a phylogenetic trees were generated by the neighbor-joining method using MEGA5 software with 1000 bootstrap permutations .
Measurements of transcript abundance
For microarray analysis, different growth regimes were used for leaves and stems, versus roots. For leaves and stems, approximately 30 days following germination and growth on soil, total leaf and stem were collected when the inflorescence began to emerge from the flag leaf. Leaves were separated from the stems with a curved-tip probe. Nodes and internodes from the second leaf junction to the internode below the inflorescence were frozen in liquid nitrogen. For roots, seven-day-old whole seedlings were flash frozen in liquid nitrogen and then the roots were snapped off into a sterile culture tube. For all organs, material was harvested six times during the day (2, 6, 10, 14, 18, and 22 h circadian time). Three plants were dissected for each time point and in triplicate for each tissue type. Samples were stored in liquid nitrogen or at −80°C until RNA extraction. Tissue was ground with mortar and pestle in liquid nitrogen. RNA was extracted using a kit (Plant RNaeasy, Qiagen, Valencia, CA) according to the manufacturer’s instructions. For hybridization, cDNA probes were synthesized using a kit (WT Ambion Santa Clara, CA).
The probes were applied to the B. distachyon BradiAR1b520742 whole genome tiling array (Affymetrix, Santa Clara, CA). The array contains ~6.5 million unique 25-mer oligonucleotide features, both the forward and reverse strand sequence. The complete genome sequence is tiled with an average of 30 bases between each array feature; 1.6 million features correspond to exons and introns and 4.9 million features between gene models (Todd Mockler, Donald Danforth Plant Science Center, personal communication). Approximately ~95% (~26,670) of the genes have at least five corresponding exon array features and from those a summary value was calculated for each gene model. Probeset values were calculated using gcRMA .
For obtaining RNA from the transgenics, stems were collected at the same developmental stage when the inflorescence was just visible from the flag leaf. First, second, and third nodes and internodes of the tallest stem were frozen in liquid nitrogen, homogenized, and RNA extracted as described above.
For RT-QPCR, on-column DNA digestion was performed using RNase-free DNase I (Qiagen). First strand cDNA was synthesized from 1 μg of total RNA using Superscript III reverse transcriptase with oligo dT primers (Invitrogen, Grand Island, NY). Samples were diluted three-fold with RNase free water (Qiagen) and 1 μL from each cDNA sample was used for RT-QPCR to check for genomic DNA contamination using GapDH primers. Triplicate quantitative PCR reactions were performed using 20 μL reaction volumes with 1 μL of diluted cDNA in each reaction with the QuantiFast SYBR Green PCR Kit (Qiagen). The reactions were conducted in an Eppendorff Realplex2 Mastercycler using the following conditions: 95°C for 2 min, followed by 40 cycles of 95°C for 15 s, 60°C for 15 s and 68°C for 20 s. As reference genes for normalization, BdUBC18 (ubiquitin-conjugating enzyme 18) and Bd5g25870 (Belonging to nuclear hormone receptor binding category) were used . QuantiPrime primer design tool was used for qPCR primer design .
RNA in situ hybridization
RNA in situ hybridization was performed using methods described previously [42, 43]. Briefly, the 3′ end of BdCESA4 and BdCESA7 were cloned into the pGEM-T Easy vector and used as the template to generate labeled sense and anti-sense ribo probes using a kit (DIG labeling kit, Roche, Indianapolis, IN). Three week old stem sections were frozen at −80°C and fixed in 4% paraformaldehyde/ethanol:acetic acid (3:1) overnight at 4°C . Fixed tissue was encased in 4% agarose and sectioned using a Vibratome. Sections were removed from agarose using forceps and washed in phosphate-buffered saline (PBS; 33 mM Na2HPO4, 1.8 mM NaH2PO4 and 140 mM NaCl, pH 7.2) and post-fixed in PBS containing 3.7% (v/v) formaldehyde and 2 mg/mL glycine for 20 min at room temperature. Sections were washed in PBS and dehydrated in graded ethanol series and pre-hybridized for 1 h at 65°C with hybridization solution comprising 20X SSC (3M NaCl, 0.3M Na citrate) containing 20% SDS, 3.7% formamide, and 10 mg/mL yeast tRNA made up in DEPC-treated water. Stems were then hybridized with DIG-labeled sense and anti-sense probes overnight at 65°C and washed with a series of SSC and SDS containing solutions. Anti-DIG alkaline phosphatase fab fragments (Roche) were used at a dilution of 1:1000 and incubated over night at 4°C. Alkaline phosphatase was detected using nitroblue tetrazolium and 5-bromo-4-chloro-3-indoyl-phosphate (Roche) in 0.1 M Tris (pH 9.5), 0.1 M NaCl and imaged with a PixeLINK camera attached to a Nikon eclipse 200 microscope.
Artificial microRNA constructs
Artificial microRNA sequences were designed on the Web MicroRNA Designer platform (http://wmd3.weigelworld.org) based on JGI B. distachyon genome annotation version 1.0 . The amiR-CESA4 construct targets AAGGGACCCATTCTTAAGCCA with hybridization energy of −37.74 kcal/mol and the amiR-CESA7 construct targets ACGCCCACCATTGTCATCATC with hybridization energy of −37.37 kcal/mol. Constructs were engineered from the pNW55 plasmid to replace the targeting regions of the native rice microRNA precursor osaMIR528. MicroRNA targets were PCR amplified (Additional file 5: Table S2) according to Warthmann et al.  and cloned into pENTR/D-TOPO (Invitrogen). Sequence confirmed clones were recombined with a modified version of the destination vector pOL001 , pOL001-ubigate-ori1 and transformed into Agrobaterium tumefaciens strain AGL1 via electroporation.
Transformation was carried out according to Vogel et al. . Briefly, seeds were collected from six to seven week old plants and deglumed. Seeds were surface sterilized with a 1.3% NaClO solution containing 0.01% Triton-X100 for 4 min. Embryos were dissected and placed on callus initiation medium under aseptic conditions. Calli were propagated for seven weeks with two subsequent subcultures at four and six weeks following dissection. Seven-week-old calli were immersed in an A. tumefaciens suspension for 5 min and dried on filter paper. Next, they were co-cultivated on dry filter paper for three days at 22°C in the dark. Following co-cultivation, calli were moved to selective plates containing 40 mg/L hygromycin and 200 mg/L timentin for four weeks in the dark at 28°C. Following selection, calli were moved to Linsmaier and Skoog media for regeneration at 28°C under constant light and next onto Murashige and Skoog media for root establishment under same conditions. Next they were transplanted to soil and grown as described above.
Genomic DNA extraction and genotyping
Genomic DNA was extracted from leaves according to Csaikl et al.  with slight modification. Briefly, leaves were frozen in liquid nitrogen and ground using 3.2 mm diameter stainless steel metal balls (Biospec Products, Bartlesville, OK) in a ball mill (Mixer Mill, MM400, Retsch, Newtown, PA). Ground samples were incubated in DNA extraction buffer (100 mM NaCl, 50 mM Tris, 25 mM EDTA, 1% SDS, 10 mM 2-mercaptoethanol) for 10 min at 65°C. Next, they were mixed with 5 M potassium acetate and incubated on ice for 20 min and centrifuged for 10 min. DNA was pelleted by mixing the supernatant with isopropanol and centrifuging at maximum speed for 10 min followed by a 70% ethanol wash. Pelleted DNA was resuspended in 1X TE and the integrity of the samples was measured using a spectrophotometer (NanoDrop 1000, Themo Scientific, Waltham, MA). Genotyping was carried out by PCR for the Hgyromycin locus using the extracted genomic DNA as template with following cycler conditions; Initial denaturation at 98°C for 2 min, flowed by 30 cycles of 98°C for 30 s, 59°C for 30 s, 72°C for 55 s and a final extension at 72°C for 7 min. PCR confirmed positive transformants were used in subsequent experiments.
For histochemical analysis, stems were hand sectioned using a razor blade and stained with 0.002% Toluidine blue in water for 30 s. Stained sections were mounted with water and observed under an Eclipse E200MV R microscope (Nikon) and imaged using a PixeLINK 3 MP camera. Images were captured at 4 X magnification and stem area was measured by freehand tracing of a perimeter in ImageJ (http://rsb.info.nih.gov/ij/). Images captured at 20 X magnification were used for cell wall thickness measurements.
For polarized-light microscopy, internode segments were fixed in 2% glutaraldehyde in 50 mM Na2PO4 buffer (pH 7.2) for a minimum of 2 h at room temperature and then handled and embedded in Spurr’s resin by standard techniques. Semi-thin sections (0.5 μm) were cut on an ultra-microtome, mounted in immersion oil, and observed through a microscope equipped with the LC-PolScope (CRI Cambridge MA) as described previously . Briefly, this instrument uses circularly polarized light and computer-controlled liquid crystal compensators to obtain four images with known compensator settings and from them calculates a fifth image in which the intensity of each pixel is proportional to birefringent retardation and a sixth image (not used here) in which the intensity of each pixel represents the orientation of the crystal’s optical axis within the specimen plane . Several sections of each genotype were observed and the images shown are representative.
X-ray diffraction profiles and sum-frequency-generation vibration spectroscopy
X-ray diffraction and calculation of a crystallinity index were done as described by Ruland et al.  with slight modifications. Fully senesced stems were ground as described above for genomic DNA extraction. Diffraction was analyzed on an X’Pert Pro powder X-ray diffractometer (PANanalytical BV, The Netherlands) operated at 45 kV and 45 mA using CuKα radiation at both Kα1 (λ = 1.5406 Å) and Kα2 (λ = 1.5444 Å). The diffraction profile was acquired from 5 to 50˚ in 0.0167 steps, with 66 s per step. The crystallinity index was calculated using the amorphous subtraction method, which determines crystallinity by subtracting the amorphous contribution from the diffraction profile obtained from xylan (Aldrich) measured in parallel. A scale factor was applied to the xylan profile to avoid negative values in the subtracted profiles. For each group, eight to twelve individuals were analyzed. Sum-frequency vibration spectroscopy was done as previously described . Intact first internodes were excised from completely senesced plants, and each group contained four to six individuals with ten measurements per stem. The amount of crystalline cellulose was estimated by comparing the 2944 cm-1 intensity of the sample with that of Avicel, as previously described .
For each measurement, three to twelve independent plants were sampled from three or four different T3 families for each transgene. Student’s t-tests were performed in R (v 2.15.0). Significance was set a P < 0.05. No significant differences were observed among the different independent transgenic lines for a given construct and were thus pooled.
Availability of supporting data
The microarray raw data (CEL files) have been deposited at PLEXdb [Accession no: BD3] .
This research was supported by the Division of Chemical Sciences, Geosciences, and Biosciences, Office of Basic Energy Sciences of the U.S. Department of Energy through Grant DE-FG-03ER15421 to TIB and by the Office of Science (BER) Department of Energy Grant DE-FG02-08ER64700DE to SPH. The spectroscopic analysis of cellulose crystallinity was supported by The Center for Lignocellulose Structure and Formation, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, and Office of Basic Energy Sciences under Award Number DE-SC0001090.
- Heaton EA, Flavell RB, Mascia PN, Thomas SR, Dohleman FG, Long SP: Herbaceous energy crop development: recent progress and future prospects. Curr Opin Biotechnol. 2008, 19: 202-10.1016/j.copbio.2008.05.001.PubMedView ArticleGoogle Scholar
- Brkljacic J, Grotewold E, Scholl R, Mockler T, Garvin DF, Vain P, Brutnell T, Sibout R, Bevan M, Budak H, et al: Brachypodium as a model for the grasses: today and the future. Plant Physiol. 2011, 157 (1): 3-13. 10.1104/pp.111.179531.PubMedPubMed CentralView ArticleGoogle Scholar
- Bragg JN, Wu J, Gordon SP, Guttman ME, Thilmony R, Lazo GR, Gu YQ, Vogel JP: Generation and characterization of the Western Regional Research Center Brachypodium T-DNA insertional mutant collection. PLoS ONE. 2012, 7 (9): e41916-10.1371/journal.pone.0041916.PubMedPubMed CentralView ArticleGoogle Scholar
- d'Yvoire MB, Bouchabke-Coussa O, Voorend W, Antelme S, Cézard L, Legée F, Lebris P, Legay S, Whitehead C, McQueen-Mason SJ, et al: Disrupting the cinnamyl alcohol dehydrogenase 1 gene (BdCAD1) leads to altered lignification and improved saccharification in Brachypodium distachyon. Plant J. 2012, 73: 496-508.View ArticleGoogle Scholar
- International Brachypodium Initiative: Genome sequencing and analysis of the model grass Brachypodium distachyon. Nature. 2010, 463 (7282): 763-768. 10.1038/nature08747.View ArticleGoogle Scholar
- Gomez L, Bristow J, Statham E, McQueen-Mason S: Analysis of saccharification in Brachypodium distachyon stems under mild conditions of hydrolysis. Biotechnol Biofuels. 2008, 1: 15-10.1186/1754-6834-1-15.PubMedPubMed CentralView ArticleGoogle Scholar
- Christensen U, Alonso-Simon A, Scheller HV, Willats WGT, Harholt J: Characterization of the primary cell walls of seedlings of Brachypodium distachyon - potential model plant for temperate grasses. Phytochemistry. 2010, 71 (1): 62-69. 10.1016/j.phytochem.2009.09.019.PubMedView ArticleGoogle Scholar
- Guillon F, Bouchet B, Jamme F, Robert P, Quemener B, Barron C, Larre C, Dumas P, Saulnier L: Brachypodium distachyon grain: characterization of endosperm cell walls. J Exp Bot. 2011, 62 (3): 1001-1015. 10.1093/jxb/erq332.PubMedView ArticleGoogle Scholar
- Somerville C: Cellulose synthesis in higher plants. Annu Rev Cell Dev Biol. 2006, 22: 53-10.1146/annurev.cellbio.22.022206.160206.PubMedView ArticleGoogle Scholar
- Haigler CH, Brown RM: Transport of rosettes from the Golgi apparatus to the plasma membrane in isolated mesophyll cells of Zinnia elegans during differentiation to tracheary elements in suspension culture. Protoplasma. 1986, 134: 111-10.1007/BF01275709.View ArticleGoogle Scholar
- Kimura S, Laosinchai W, Itoh T, Cui XJ, Linder CR, Brown RM: Immunogold labeling of rosette terminal cellulose-synthesizing complexes in the vascular plant Vigna angularis. Plant Cell. 2075, 1999: 11-Google Scholar
- Doblin MS, Kurek I, Jacob-Wilk D, Delmer DP: Cellulose biosynthesis in plants: from genes to rosettes. Plant Cell Physiol. 2002, 43: 1407-10.1093/pcp/pcf164.PubMedView ArticleGoogle Scholar
- Carroll A, Specht C: Understanding plant cellulose synthases through a comprehensive investigation of the cellulose synthase family sequences. Front Plant Sci. 2011, 2: 5-PubMedPubMed CentralView ArticleGoogle Scholar
- Yin YB, Huang JL, Xu Y: The cellulose synthase superfamily in fully sequenced plants and algae. BMC Plant Biol. 2009, 9: 99-113. 10.1186/1471-2229-9-99.PubMedPubMed CentralView ArticleGoogle Scholar
- Desprez T, Juraniec M, Crowell EF, Jouy H, Pochylova Z, Parcy F, Hofte H, Gonneau M, Vernhettes S: Organization of cellulose synthase complexes involved in primary cell wall synthesis in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2007, 104: 15572-15577. 10.1073/pnas.0706569104.PubMedPubMed CentralView ArticleGoogle Scholar
- Persson S, Paredez A, Carroll A, Palsdottir H, Doblin M, Poindexter P, Khitrov N, Auer M, Somerville CR: Genetic evidence for three unique components in primary cell-wall cellulose synthase complexes in Arabidopsis. Proc Natl Acad Sci U S A. 2007, 104: 15566-15571. 10.1073/pnas.0706592104.PubMedPubMed CentralView ArticleGoogle Scholar
- Taylor NG, Scheible WR, Cutler S, Somerville CR, Turner SR: The irregular xylem3 locus of Arabidopsis encodes a cellulose synthase required for secondary cell wall synthesis. Plant Cell. 1999, 11: 769-PubMedPubMed CentralView ArticleGoogle Scholar
- Taylor NG, Laurie S, Turner SR: Multiple cellulose synthase catalytic subunits are required for cellulose synthesis in Arabidopsis. Plant Cell. 2000, 12 (12): 2529-2540.PubMedPubMed CentralView ArticleGoogle Scholar
- Taylor NG, Howells RM, Huttly AK, Vickers K, Turner SR: Interactions among three distinct CesA proteins essential for cellulose synthesis. Proc Natl Acad Sci U S A. 2003, 100: 1450-1455. 10.1073/pnas.0337628100.PubMedPubMed CentralView ArticleGoogle Scholar
- Li A, Xia T, Xu W, Chen T, Li X, Fan J, Wang R, Feng S, Wang Y, Wang B, et al: An integrative analysis of four CESA isoforms specific for fiber cellulose production between Gossypium hirsutum and Gossypium barbadense. Planta. 2013, 237 (6): 1585-1597. 10.1007/s00425-013-1868-2.PubMedView ArticleGoogle Scholar
- Song D, Shen J, Li L: Characterization of cellulose synthase complexes in Populus xylem differentiation. New Phytol. 2010, 187 (3): 777-790. 10.1111/j.1469-8137.2010.03315.x.PubMedView ArticleGoogle Scholar
- Burton RA, Ma G, Baumann U, Harvey AJ, Shirley NJ, Taylor J, Pettolino F, Bacic A, Beatty M, Simmons CR, et al: A customized gene expression microarray reveals that the brittle stem phenotype fs2 of barley Is attributable to a retroelement in the HvCesA4 cellulose synthase gene. Plant Physiol. 2010, 153 (4): 1716-1728. 10.1104/pp.110.158329.PubMedPubMed CentralView ArticleGoogle Scholar
- Tanaka K, Murata K, Yamazaki M, Onosato K, Miyao A, Hirochika H: Three distinct rice cellulose synthase catalytic subunit genes required for cellulose synthesis in the secondary wall. Plant Physiol. 2003, 133 (1): 73-83. 10.1104/pp.103.022442.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang L, Guo K, Li Y, Tu Y, Hu H, Wang B, Cui X, Peng L: Expression profiling and integrative analysis of the CESA/CSL superfamily in rice. BMC Plant Biol. 2010, 10 (1): 282-10.1186/1471-2229-10-282.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang B, Deng L, Qian Q, Xiong G, Zeng D, Li R, Guo L, Li J, Zhou Y: A missense mutation in the transmembrane domain of CESA4 affects protein abundance in the plasma membrane and results in abnormal cell wall biosynthesis in rice. Plant Mol Biol. 2009, 71 (4): 509-524.PubMedView ArticleGoogle Scholar
- Barnette AL, Bradley LC, Veres BD, Schreiner EP, Park YB, Park J, Park S, Kim SH: Selective detection of crystalline cellulose in plant cell walls with sum-frequency-generation (SFG) vibration spectroscopy. Biomacromolecules. 2011, 12 (7): 2434-2439. 10.1021/bm200518n.PubMedView ArticleGoogle Scholar
- Barnette AL, Lee C, Bradley LC, Schreiner EP, Park YB, Shin H, Cosgrove DJ, Park S, Kim SH: Quantification of crystalline cellulose in lignocellulosic biomass using sum frequency generation (SFG) vibration spectroscopy and comparison with other analytical methods. Carbohydr Polym. 2012, 89 (3): 802-809. 10.1016/j.carbpol.2012.04.014.PubMedView ArticleGoogle Scholar
- Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, Haberer G, Hellsten U, Mitros T, Poliakov A, et al: The Sorghum bicolor genome and the diversification of grasses. Nature. 2009, 457 (7229): 551-556. 10.1038/nature07723.PubMedView ArticleGoogle Scholar
- Chaw S-M, Chang C-C, Chen H-L, Li W-H: Dating the monocot–dicot divergence and the origin of core eudicots using whole chloroplast genomes. J Mol Evol. 2004, 58 (4): 424-441. 10.1007/s00239-003-2564-9.PubMedView ArticleGoogle Scholar
- Appenzeller L, Doblin M, Barreiro R, Wang H, Niu X, Kollipara K, Carrigan L, Tomes D, Chapman M, Dhugga KS: Cellulose synthesis in maize: isolation and expression analysis of the cellulose synthase (CesA) gene family. Cellulose. 2004, 11 (3): 287-299.View ArticleGoogle Scholar
- Burton RA, Shirley NJ, King BJ, Harvey AJ, Fincher GB: The CesA gene family of barley (Hordeum vulgare): quantitative analysis of transcripts reveals two groups of co-expressed genes. Plant Physiol. 2004, 134: 224-236. 10.1104/pp.103.032904.PubMedPubMed CentralView ArticleGoogle Scholar
- Turner SR, Somerville CR: Collapsed xylem phenotype of Arabidopsis identifies mutants deficient in cellulose deposition in the secondary cell wall. Plant Cell. 1997, 9 (5): 689-701.PubMedPubMed CentralView ArticleGoogle Scholar
- Kokubo A, Sakurai N, Kuraishi S, Takeda K: Culm brittleness of barley (Hordeum vulgare L.) mutants is caused by smaller number of cellulose molecules in cell wall. Plant Physiol. 1991, 97 (2): 509-514. 10.1104/pp.97.2.509.PubMedPubMed CentralView ArticleGoogle Scholar
- Vogel J: Unique aspects of the grass cell wall. Curr Opin Plant Biol. 2008, 11 (3): 301-307. 10.1016/j.pbi.2008.03.002.PubMedView ArticleGoogle Scholar
- Handakumbura PP, Hazen SP: Transcriptional regulation of grass secondary cell wall biosynthesis: playing catch-up with Arabidopsis thaliana. Front Plant Sci. 2012, 3: 74-PubMedPubMed CentralView ArticleGoogle Scholar
- Taylor NG: Cellulose biosynthesis and deposition in higher plants. New Phytol. 2008, 178 (2): 239-252. 10.1111/j.1469-8137.2008.02385.x.PubMedView ArticleGoogle Scholar
- Hématy K, Sado P-E, Van Tuinen A, Rochange S, Desnos T, Balzergue S, Pelletier S, Renou J-P, Höfte H: A receptor-like kinase mediates the response of Arabidopsis cells to the inhibition of cellulose synthesis. Curr Biol. 2007, 17 (11): 922-931. 10.1016/j.cub.2007.05.018.PubMedView ArticleGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28 (10): 2731-2739. 10.1093/molbev/msr121.PubMedPubMed CentralView ArticleGoogle Scholar
- Wu Z, Irizarry RA, Gentleman R, Martinez-Murillo F, Spencer F: A model-based background adjustment for oligonucleotide expression arrays. J Am Statis Assoc. 2004, 99 (468): 909-917. 10.1198/016214504000000683.View ArticleGoogle Scholar
- Hong S-Y, Seo P, Yang M-S, Xiang F, Park C-M: Exploring valid reference genes for gene expression studies in Brachypodium distachyon by real-time PCR. BMC Plant Biol. 2008, 8 (1): 112-10.1186/1471-2229-8-112.PubMedPubMed CentralView ArticleGoogle Scholar
- Arvidsson S, Kwasniewski M, Riano-Pachon D, Mueller-Roeber B: QuantPrime - a flexible tool for reliable high-throughput primer design for quantitative PCR. BMC Bioinforma. 2008, 9 (1): 465-10.1186/1471-2105-9-465.View ArticleGoogle Scholar
- Harding SA, Leshkevich J, Chiang VL, Tsai C-J: Differential substrate inhibition couples kinetically distinct 4-Coumarate:coenzyme A Ligases with spatially distinct metabolic roles in quaking aspen. Plant Physiol. 2002, 128 (2): 428-438. 10.1104/pp.010603.PubMedPubMed CentralView ArticleGoogle Scholar
- Kao Y-Y, Harding SA, Tsai C-J: Differential expression of two distinct phenylalanine ammonia-lyase genes in condensed tannin-accumulating and lignifying cells of quaking aspen. Plant Physiol. 2002, 130 (2): 796-807. 10.1104/pp.006262.PubMedPubMed CentralView ArticleGoogle Scholar
- Harrington M, Hong E, Fasanmi O, Brewster R: Cadherin-mediated adhesion regulates posterior body formation. BMC Dev Biol. 2007, 7 (1): 130-10.1186/1471-213X-7-130.PubMedPubMed CentralView ArticleGoogle Scholar
- Warthmann N, Chen H, Ossowski S, Weigel D, Hervre P: Highly specific gene silencing by artificial miRNAs in rice. PLoS ONE. 2008, 3 (3): e1829-10.1371/journal.pone.0001829.PubMedPubMed CentralView ArticleGoogle Scholar
- Vogel J, Garvin D, Leong O, Hayden D: Agrobacterium-mediated transformation and inbred line development in the model grass Brachypodium distachyon. Plant Cell Tiss Org. 2006, 84 (2): 100179-100191.View ArticleGoogle Scholar
- Vogel J, Hill T: High-efficiency Agrobacterium-mediated transformation of Brachypodium distachyon inbred line Bd21-3. Plant Cell Rep. 2008, 27 (3): 471-478. 10.1007/s00299-007-0472-y.PubMedView ArticleGoogle Scholar
- Csaikl UM, Bastian H, Brettschneider R, Gauch S, Meir A, Schauerte M, Scholz F, Sperisen C, Vornam B, Ziegenhagen B: Comparative analysis of different DNA extraction protocols: a fast, universal maxi-preparation of high quality plant DNA for genetic evaluation and phylogenetic studies. Plant Mol Biol Rep. 1998, 16 (1): 69-86. 10.1023/A:1007428009556.View ArticleGoogle Scholar
- Baskin TI, Beemster GTS, Judy-March JE, Marga F: Disorganization of cortical microtubules stimulates tangential expansion and reduces the uniformity of cellulose microfibril alignment among cells in the root of Arabidopsis. Plant Physiol. 2004, 135: 2279-10.1104/pp.104.040493.PubMedPubMed CentralView ArticleGoogle Scholar
- Oldenbourg R, Mei G: New polarized light microscope with precision universal compensator. J Microsc. 1995, 180 (2): 140-147. 10.1111/j.1365-2818.1995.tb03669.x.PubMedView ArticleGoogle Scholar
- Ruland W: X-ray determination of crystallinity and diffuse disorder scattering. Acta Cryst. 1961, 14: 1180-1185. 10.1107/S0365110X61003429.View ArticleGoogle Scholar
- Dash S, Van Hemert J, Hong L, Wise RP, Dickerson JA: PLEXdb: gene expression resources for plants and plant pathogens. Nucleic Acids Res. 2011, 40: D1194-D1201.PubMedPubMed CentralView ArticleGoogle Scholar