An integrated genomic and metabolomic framework for cell wall biology in rice

Down-regulation of OsMYB103L distinctively alters beta-1,4-glucan polymerization and cellulose microfibers assembly for enhanced biomass enzymatic saccharification in rice

BACKGROUND

As a major component of plant cell walls, cellulose provides the most abundant biomass resource convertible for biofuels. Since cellulose crystallinity and polymerization have been characterized as two major features accounting for lignocellulose recalcitrance against biomass enzymatic saccharification, genetic engineering of cellulose biosynthesis is increasingly considered as a promising solution in bioenergy crops. Although several transcription factors have been identified to regulate cellulose biosynthesis and plant cell wall formation, much remains unknown about its potential roles for genetic improvement of lignocellulose recalcitrance.

RESULTS

In this study, we identified a novel rice mutant (Osfc9/myb103) encoded a R2R3-MYB transcription factor, and meanwhile generated OsMYB103L-RNAi-silenced transgenic lines. We determined significantly reduced cellulose levels with other major wall polymers (hemicellulose, lignin) slightly altered in mature rice straws of the myb103 mutant and RNAi line, compared to their wild type (NPB). Notably, the rice mutant and RNAi line were of significantly reduced cellulose features (crystalline index/CrI, degree of polymerization/DP) and distinct cellulose nanofibers assembly. These alterations consequently improved lignocellulose recalcitrance for significantly enhanced biomass enzymatic saccharification by 10-28% at p < 0.01 levels (n = 3) after liquid hot water and chemical (1% H₂SO₄, 1% NaOH) pretreatments with mature rice straws. In addition, integrated RNA sequencing with DNA affinity purification sequencing (DAP-seq) analyses revealed that the OsMYB103L might specifically mediate cellulose biosynthesis and deposition by regulating OsCesAs and other genes associated with microfibril assembly.

CONCLUSIONS

This study has demonstrated that down-regulation of OsMYB103L could specifically improve cellulose features and cellulose nanofibers assembly to significantly enhance biomass enzymatic saccharification under green-like and mild chemical pretreatments in rice. It has not only indicated a powerful strategy for genetic modification of plant cell walls in bioenergy crops, but also provided insights into transcriptional regulation of cellulose biosynthesis in plants.

The Divergent Roles of the Rice bcl-2 Associated Athanogene (BAG) Genes in Plant Development and Environmental Responses

Bcl-2-associated athanogene (BAG), a group of proteins evolutionarily conserved and functioned as co-chaperones in plants and animals, is involved in various cell activities and diverse physiological processes. However, the biological functions of this gene family in rice are largely unknown. In this study, we identified a total of six BAG members in rice. These genes were classified into two groups, OsBAG1, -2, -3, and -4 are in group I with a conserved ubiquitin-like structure and OsBAG5 and -6 are in group Ⅱ with a calmodulin-binding domain, in addition to a common BAG domain. The BAG genes exhibited diverse expression patterns, with OsBAG4 showing the highest expression level, followed by OsBAG1 and OsBAG3, and OsBAG6 preferentially expressed in the panicle, endosperm, and calli. The co-expression analysis and the hierarchical cluster analysis indicated that the OsBAG1 and OsBAG3 were co-expressed with primary cell wall-biosynthesizing genes, OsBAG4 was co-expressed with phytohormone and transcriptional factors, and OsBAG6 was co-expressed with disease and shock-associated genes. β-glucuronidase (GUS) staining further indicated that OsBAG3 is mainly involved in primary young tissues under both primary and secondary growth. In addition, the expression of the BAG genes under brown planthopper (BPH) feeding, N, P, and K deficiency, heat, drought and plant hormones treatments was investigated. Our results clearly showed that OsBAGs are multifunctional molecules as inferred by their protein structures, subcellular localizations, and expression profiles. BAGs in group I are mainly involved in plant development, whereas BAGs in group II are reactive in gene regulations and stress responses. Our results provide a solid basis for the further elucidation of the biological functions of plant BAG genes.

Association mapping identifies quantitative trait loci (QTL) for digestibility in rice straw

BACKGROUND

The conversion of lignocellulosic biomass from agricultural waste into biofuels and chemicals is considered a promising way to provide sustainable low carbon products without compromising food security. However, the use of lignocellulosic biomass for biofuel and chemical production is limited by the cost-effectiveness of the production process due to its recalcitrance to enzymatic hydrolysis and fermentable sugar release (i.e., saccharification). Rice straw is a particularly attractive feedstock because millions of tons are currently burned in the field each year for disposal. The aim of this study was to explore the underlying natural genetic variation that impacts the recalcitrance of rice (Oryza sativa) straw to enzymatic saccharification. Ultimately, we wanted to investigate whether we could identify genetic markers that could be used in rice breeding to improve commercial cultivars for this trait. Here, we describe the development and characterization of a Vietnamese rice genome-wide association panel, high-throughput analysis of rice straw saccharification and lignin content, and the results from preliminary genome-wide association studies (GWAS) of the combined data sets. We identify both QTL and plausible candidate genes that may have an impact on the saccharification of rice straw.

RESULTS

We assembled a diversity panel comprising 151 rice genotypes (Indica and Japonica types) from commercial, historical elite cultivars, and traditional landraces grown in Vietnam. The diversity panel was genotyped using genotype by sequencing (GBS) methods yielding a total of 328,915 single nucleotide polymorphisms (SNPs). We collected phenotypic data from stems of these 151 genotypes for biomass saccharification and lignin content. Using GWAS on the indica genotypes over 2 years we identified ten significant QTL for saccharification (digestibility) and seven significant QTL for lignin. One QTL on chromosome 11 occurred in both GWAS for digestibility and for lignin. Seven QTL for digestibility, on CH2, CH6, CH7, CH8, and CH11, were observed in both years of the study. The QTL regions for saccharification include three potential candidate genes that have been previously reported to influence digestibility: OsAT10; OsIRX9; and OsMYB58/63-L.

CONCLUSIONS

Despite the difficulties associated with multi-phasic analysis of complex traits in novel germplasm, a moderate resolution GWAS successfully identified genetic associations encompassing both known and/or novel genes involved in determining the saccharification potential and lignin content of rice straw. Plausible candidates within QTL regions, in particular those with roles in cell wall biosynthesis, were identified but will require validation to confirm their value for application in rice breeding.

AtCSLD3 and GhCSLD3 mediate root growth and cell elongation downstream of the ethylene response pathway in Arabidopsis

CSLD3, a gene of the cellulose synthase-like D family, affects root hair elongation, but its interactions with ethylene signaling and phosphate-starvation are poorly understood. Here, we aim to understand the role of CSLD3 in the context of the ethylene signaling and phosphate starvation pathways in Arabidopsis plant growth. Therefore, we performed a comparative analysis of the csld3-1 mutant, CSLD3-overexpressing lines, and ethylene-response mutants, such as the constitutive ethylene-response mutant i-ctr1. We found that CSLD3 overexpression enhanced root and hypocotyl growth by increasing cell elongation, and that the root growth was highly sensitive to ethylene treatment (1 µM ACC), in particular under phosphate starvation. However, the CSLD3-mediated hypocotyl elongation occurred independently of the ethylene signaling pathway. Notably, the typical induction of root hair and root elongation by ethylene and phosphate-starvation was completely abolished in the csld3-1 mutant. Furthermore, i-ctr1 csld3-1 double-mutants were hairless like the csld3-1 parent, confirming that CSLD3 acts downstream of the ethylene signaling pathway during root growth. Moreover, the CSLD3 levels positively correlated with cellulose levels, indicating a role of CSLD3 in cellulose synthesis, which may explain the observed growth effects. Our results establish how CSLD3 works in the context of the ethylene signaling and phosphate-starvation pathways during root hair growth, cell elongation, and cell wall biosynthesis.

Evolution, gene expression profiling and 3D modeling of CSLD proteins in cotton

BACKGROUND

Among CESA-like gene superfamily, the cellulose synthase-like D (CSLD) genes are most similar to cellulose synthase genes and have been reported to be involved in tip-growing cell and stem development. However, there has been no genome-wide characterization of this gene subfamily in cotton. We thus sought to analyze the evolution and functional characterization of CSLD proteins in cotton based on fully sequenced cotton genomes.

RESULTS

A total of 23 full-length CSLD proteins were identified in Gossypium raimondii, Gossypium arboreum and Gossypium hirsutum. The phylogenetic tree divided the CSLD proteins into five clades with strong support: CSLD1, CSLD2/3, CSLD4, CSLD5 and CSLD6. The total expression of GhCSLD genes was the highest in androecium & gynoecium (mostly contributed by CSLD1 and CSLD4) compared with other CSL genes. CSLD1 and CSLD4 were only highly expressed in androecium & gynoecium (A&G), and showed tissue-specific expression. The total expression of CSLD2/3, 5 and 6 was highest in the specific tissues. These results suggest that CSLD genes showed the different pattern of expression. Cotton CSLD proteins were subjected to different evolutionary pressures, and the CSLD1 and CSLD4 proteins exhibited episodic and long-term shift positive selection. The predicted three-dimensional structure of GrCSLD1 suggested that GrCSLD1 belongs to glycosyltransferase family 2. The amino acid residues under positive selection in the CSLD1 lineage are positioned in a region adjacent to the class-specific region (CSR), β1-strand and transmembrane helices (TMHs) in the GrCSLD1structure.

CONCLUSION

Our results characterized the CSLD proteins by an integrated approach containing phylogeny, transcriptional profiling and 3D modeling. The study added to the understanding about the importance of the CSLD family and provide a useful reference for selecting candidate genes and their associations with the biosynthesis of the cell wall in cotton.

Processes Underlying a Reproductive Barrier in indica-japonica Rice Hybrids Revealed by Transcriptome Analysis

In rice (Oryza sativa), hybrids between indica and japonica subspecies are usually highly sterile, which provides a model system for studying postzygotic reproductive isolation. A killer-protector system, S5, composed of three adjacent genes (ORF3, ORF4, and ORF5), regulates female gamete fertility of indica-japonica hybrids. To characterize the processes underlying this system, we performed transcriptomic analyses of pistils from rice variety Balilla (BL), Balilla with transformed ORF5+ (BL5+) producing sterile female gametes, and Balilla with transformed ORF3+ and ORF5+ (BL3+5+) producing fertile gametes. RNA sequencing of tissues collected before (MMC), during (MEI), and after (AME) meiosis of the megaspore mother cell detected 19,269 to 20,928 genes as expressed. Comparison between BL5+ and BL showed that ORF5+ induced differential expression of 8,339, 6,278, and 530 genes at MMC, MEI, and AME, respectively. At MMC, large-scale differential expression of cell wall-modifying genes and biotic and abiotic response genes indicated that cell wall integrity damage induced severe biotic and abiotic stresses. The processes continued to MEI and induced endoplasmic reticulum (ER) stress as indicated by differential expression of ER stress-responsive genes, leading to programmed cell death at MEI and AME, resulting in abortive female gametes. In the BL3+5+/BL comparison, 3,986, 749, and 370 genes were differentially expressed at MMC, MEI, and AME, respectively. Large numbers of cell wall modification and biotic and abiotic response genes were also induced at MMC but largely suppressed at MEI without inducing ER stress and programed cell death , producing fertile gametes. These results have general implications for the understanding of biological processes underlying reproductive barriers.

CmMYB19 Over-Expression Improves Aphid Tolerance in Chrysanthemum by Promoting Lignin Synthesis

The gene encoding the MYB (v-myb avian myeloblastosis vira l oncogene homolog) transcription factor CmMYB19 was isolated from chrysanthemum. It encodes a 200 amino acid protein and belongs to the R2R3-MYB subfamily. CmMYB19 was not transcriptionally activated in yeast, while a transient expression experiment conducted in onion epidermal cells suggested that the CmMYB19 product localized to the nucleus. CmMYB19 transcription was induced by aphid (Macrosiphoniella sanborni) infestation, and the abundance of transcript was higher in the leaf and stem than in the root. The over-expression of CmMYB19 restricted the multiplication of the aphids. A comparison of transcript abundance of the major genes involved in lignin synthesis showed that CmPAL1 (phenylalanine ammonia lyase 1), CmC4H (cinnamate4 hydroxylase), Cm4CL1 (4-hydroxy cinnamoyl CoA ligase 1), CmHCT (hydroxycinnamoyl CoA-shikimate/quinate hydroxycinnamoyl transferase), CmC3H1 (coumarate3 hydroxylase1), CmCCoAOMT1 (caffeoyl CoA O-methyltransferase 1) and CmCCR1 (cinnamyl CoA reductase1) were all upregulated, in agreement with an increase in lignin content in CmMYB19 over-expressing plants. Collectively, the over-expression of CmMYB19 restricted the multiplication of the aphids on the host, mediated by an enhanced accumulation of lignin.

A precise and consistent assay for major wall polymer features that distinctively determine biomass saccharification in transgenic rice by near-infrared spectroscopy

BACKGROUND

The genetic modification of plant cell walls has been considered to reduce lignocellulose recalcitrance in bioenergy crops. As a result, it is important to develop a precise and rapid assay for the major wall polymer features that affect biomass saccharification in a large population of transgenic plants. In this study, we collected a total of 246 transgenic rice plants that, respectively, over-expressed and RNAi silenced 12 genes of the OsGH9 and OsGH10 family that are closely associated with cellulose and hemicellulose modification. We examined the wall polymer features and biomass saccharification among 246 transgenic plants and one wild-type plant. The samples presented a normal distribution applicable for statistical analysis and NIRS modeling.

RESULTS

Among the 246 transgenic rice plants, we determined largely varied wall polymer features and the biomass enzymatic saccharification after alkali pretreatment in rice straws, particularly for the fermentable hexoses, ranging from 52.8 to 95.9%. Correlation analysis indicated that crystalline cellulose and lignin levels negatively affected the hexose and total sugar yields released from pretreatment and enzymatic hydrolysis in the transgenic rice plants, whereas the arabinose levels and arabinose substitution degree (reverse xylose/arabinose ratio) exhibited positive impacts on the hexose and total sugars yields. Notably, near-infrared spectroscopy (NIRS) was applied to obtain ten equations for predicting biomass enzymatic saccharification and seven equations for distinguishing major wall polymer features. Most of the equations exhibited high R²/R²_cv/R²_ev and RPD values for a perfect prediction capacity.

CONCLUSIONS

Due to large generated populations of transgenic rice lines, this study has not only examined the key wall polymer features that distinctively affect biomass enzymatic saccharification in rice but has also established optimal NIRS models for a rapid and precise screening of major wall polymer features and lignocellulose saccharification in biomass samples. Importantly, this study has briefly explored the potential roles of a total of 12 OsGH9 and OsGH10 genes in cellulose and hemicellulose modification and cell wall remodeling in transgenic rice lines. Hence, it provides a strategy for genetic modification of plant cell walls by expressing the desired OsGH9 and OsGH10 genes that could greatly improve biomass enzymatic digestibility in rice.

Cell Wall Composition and Candidate Biosynthesis Gene Expression During Rice Development

Cell walls of grasses, including cereal crops and biofuel grasses, comprise the majority of plant biomass and intimately influence plant growth, development and physiology. However, the functions of many cell wall synthesis genes, and the relationships among and the functions of cell wall components remain obscure. To better understand the patterns of cell wall accumulation and identify genes that act in grass cell wall biosynthesis, we characterized 30 samples from aerial organs of rice (Oryza sativa cv. Kitaake) at 10 developmental time points, 3-100 d post-germination. Within these samples, we measured 15 cell wall chemical components, enzymatic digestibility and 18 cell wall polysaccharide epitopes/ligands. We also used quantitative reverse transcription-PCR to measure expression of 50 glycosyltransferases, 15 acyltransferases and eight phenylpropanoid genes, many of which had previously been identified as being highly expressed in rice. Most cell wall components vary significantly during development, and correlations among them support current understanding of cell walls. We identified 92 significant correlations between cell wall components and gene expression and establish nine strong hypotheses for genes that synthesize xylans, mixed linkage glucan and pectin components. This work provides an extensive analysis of cell wall composition throughout rice development, identifies genes likely to synthesize grass cell walls, and provides a framework for development of genetically improved grasses for use in lignocellulosic biofuel production and agriculture.

A Novel FC116/BC10 Mutation Distinctively Causes Alteration in the Expression of the Genes for Cell Wall Polymer Synthesis in Rice

We report isolation and characterization of a fragile culm mutant fc116 that displays reduced mechanical strength caused by decreased cellulose content and altered cell wall structure in rice. Map-based cloning revealed that fc116 was a base substitution mutant (G to A) in a putative beta-1,6-N-acetylglucosaminyltransferase (C2GnT) gene (LOC_Os05g07790, allelic to BC10). This mutation resulted in one amino acid missing within a newly-identified protein motif "R, RXG, RA." The FC116/BC10 gene was lowly but ubiquitously expressed in the all tissues examined across the whole life cycle of rice, and slightly down-regulated during secondary growth. This mutant also exhibited a significant increase in the content of hemicelluloses and lignins, as well as the content of pentoses (xylose and arabinose). But the content of hexoses (glucose, mannose, and galactose) was decreased in both cellulosic and non-cellulosic (pectins and hemicelluloses) fractions of the mutant. Transcriptomic analysis indicated that the typical genes in the fc116 mutant were up-regulated corresponding to xylan biosynthesis, as well as lignin biosynthesis including p-hydroxyphenyl (H), syringyl (S), and guaiacyl (G). Our results indicate that FC116 has universal function in regulation of the cell wall polymers in rice.

Genetic modification of plant cell walls to enhance biomass yield and biofuel production in bioenergy crops

Plant cell walls represent an enormous biomass resource for the generation of biofuels and chemicals. As lignocellulose property principally determines biomass recalcitrance, the genetic modification of plant cell walls has been posed as a powerful solution. Here, we review recent progress in understanding the effects of distinct cell wall polymers (cellulose, hemicelluloses, lignin, pectin, wall proteins) on the enzymatic digestibility of biomass under various physical and chemical pretreatments in herbaceous grasses, major agronomic crops and fast-growing trees. We also compare the main factors of wall polymer features, including cellulose crystallinity (CrI), hemicellulosic Xyl/Ara ratio, monolignol proportion and uronic acid level. Furthermore, the review presents the main gene candidates, such as CesA, GH9, GH10, GT61, GT43 etc., for potential genetic cell wall modification towards enhancing both biomass yield and enzymatic saccharification in genetic mutants and transgenic plants. Regarding cell wall modification, it proposes a novel groove-like cell wall model that highlights to increase amorphous regions (density and depth) of the native cellulose microfibrils, providing a general strategy for bioenergy crop breeding and biofuel processing technology.

An Integrated Bioinformatics Analysis Reveals Divergent Evolutionary Pattern of Oil Biosynthesis in High- and Low-Oil Plants

Seed oils provide a renewable source of food, biofuel and industrial raw materials that is important for humans. Although many genes and pathways for acyl-lipid metabolism have been identified, little is known about whether there is a specific mechanism for high-oil content in high-oil plants. Based on the distinct differences in seed oil content between four high-oil dicots (20~50%) and three low-oil grasses (<3%), comparative genome, transcriptome and differential expression analyses were used to investigate this mechanism. Among 4,051 dicot-specific soybean genes identified from 252,443 genes in the seven species, 54 genes were shown to directly participate in acyl-lipid metabolism, and 93 genes were found to be associated with acyl-lipid metabolism. Among the 93 dicot-specific genes, 42 and 27 genes, including CBM20-like SBDs and GPT2, participate in carbohydrate degradation and transport, respectively. 40 genes highly up-regulated during seed oil rapid accumulation period are mainly involved in initial fatty acid synthesis, triacylglyceride assembly and oil-body formation, for example, ACCase, PP, DGAT1, PDAT1, OLEs and STEROs, which were also found to be differentially expressed between high- and low-oil soybean accessions. Phylogenetic analysis revealed distinct differences of oleosin in patterns of gene duplication and loss between high-oil dicots and low-oil grasses. In addition, seed-specific GmGRF5, ABI5 and GmTZF4 were predicted to be candidate regulators in seed oil accumulation. This study facilitates future research on lipid biosynthesis and potential genetic improvement of seed oil content.

Silica distinctively affects cell wall features and lignocellulosic saccharification with large enhancement on biomass production in rice

Rice is a typical silicon-accumulating crop with enormous biomass residues for biofuels. Silica is a cell wall component, but its effect on the plant cell wall and biomass production remains largely unknown. In this study, a systems biology approach was performed using 42 distinct rice cell wall mutants. We found that silica levels are significantly positively correlated with three major wall polymers, indicating that silica is associated with the cell wall network. Silicon-supplied hydroculture analysis demonstrated that silica distinctively affects cell wall composition and major wall polymer features, including cellulose crystallinity (CrI), arabinose substitution degree (reverse Xyl/Ara) of xylans, and sinapyl alcohol (S) proportion in three typical rice mutants. Notably, the silicon supplement exhibited dual effects on biomass enzymatic digestibility in the mutant and wild type (NPB) after pre-treatments with 1% NaOH and 1% H2SO4. In addition, silicon supply largely enhanced plant height, mechanical strength and straw biomass production, suggesting that silica rescues mutant growth defects. Hence, this study provides potential approaches for silicon applications in biomass process and bioenergy rice breeding.

High-level hemicellulosic arabinose predominately affects lignocellulose crystallinity for genetically enhancing both plant lodging resistance and biomass enzymatic digestibility in rice mutants

Rice is a major food crop with enormous biomass residue for biofuels. As plant cell wall recalcitrance basically decides a costly biomass process, genetic modification of plant cell walls has been regarded as a promising solution. However, due to structural complexity and functional diversity of plant cell walls, it becomes essential to identify the key factors of cell wall modifications that could not much alter plant growth, but cause an enhancement in biomass enzymatic digestibility. To address this issue, we performed systems biology analyses of a total of 36 distinct cell wall mutants of rice. As a result, cellulose crystallinity (CrI) was examined to be the key factor that negatively determines either the biomass enzymatic saccharification upon various chemical pretreatments or the plant lodging resistance, an integrated agronomic trait in plant growth and grain production. Notably, hemicellulosic arabinose (Ara) was detected to be the major factor that negatively affects cellulose CrI probably through its interlinking with β-1,4-glucans. In addition, lignin and G monomer also exhibited the positive impact on biomass digestion and lodging resistance. Further characterization of two elite mutants, Osfc17 and Osfc30, showing normal plant growth and high biomass enzymatic digestion in situ and in vitro, revealed the multiple GH9B candidate genes for reducing cellulose CrI and XAT genes for increasing hemicellulosic Ara level. Hence, the results have suggested the potential cell wall modifications for enhancing both biomass enzymatic digestibility and plant lodging resistance by synchronically overexpressing GH9B and XAT genes in rice.