Large-scale co-expression approach to dissect secondary cell wall formation across plant species

Development of genetic markers in Eucalyptus species by target enrichment and exome sequencing

The advent of next-generation sequencing has facilitated large-scale discovery, validation and assessment of genetic markers for high density genotyping. The present study was undertaken to identify markers in genes supposedly related to wood property traits in three Eucalyptus species. Ninety four genes involved in xylogenesis were selected for hybridization probe based nuclear genomic DNA target enrichment and exome sequencing. Genomic DNA was isolated from the leaf tissues and used for on-array probe hybridization followed by Illumina sequencing. The raw sequence reads were trimmed and high-quality reads were mapped to the E. grandis reference sequence and the presence of single nucleotide variants (SNVs) and insertions/ deletions (InDels) were identified across the three species. The average read coverage was 216X and a total of 2294 SNVs and 479 InDels were discovered in E. camaldulensis, 2383 SNVs and 518 InDels in E. tereticornis, and 1228 SNVs and 409 InDels in E. grandis. Additionally, SNV calling and InDel detection were conducted in pair-wise comparisons of E. tereticornis vs. E. grandis, E. camaldulensis vs. E. tereticornis and E. camaldulensis vs. E. grandis. This study presents an efficient and high throughput method on development of genetic markers for family– based QTL and association analysis in Eucalyptus.

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

Background

Plant cell walls are complex structures that full-fill many diverse functions during plant growth and development. It is therefore not surprising that thousands of gene products are involved in cell wall synthesis and maintenance. However, functional association for the majority of these gene products remains obscure. One useful approach to infer biological associations is via transcriptional coordination, or co-expression of genes. This approach has proved useful for several biological processes. Nevertheless, combining co-expression with other large-scale measurements may improve the biological inferences.

Results

In this study, we used a combined approach of co-expression and cell wall metabolomics to obtain new insight into cell wall synthesis in rice. We initially created a weighted gene co-expression network from publicly available datasets, and then established a comprehensive cell wall dataset by determining cell wall compositions from 29 tissues that almost cover the whole life cycle of rice. We subsequently combined the datasets through the conversion of co-expressed gene modules into eigen-vectors, representing expression profiles for the genes in the modules, and performed comparative analyses against the cell wall contents. Here, we made three major discoveries. First, we confirmed our approach by finding primary and secondary wall cellulose biosynthesis modules, respectively. Second, we found co-expressed modules that strongly correlated with re-organization of the secondary cell walls and with modifications and degradation of hemicellulosic structures. Third, we inferred that at least one module is likely to play a regulatory role in the production of G-rich lignification.

Conclusions

Here, we integrated transcriptomic associations and cell wall metabolism and found that certain co-expressed gene modules are positively correlated with distinct cell wall characteristics. We propose that combining multiple data-types, such as coordinated transcription and cell wall analyses, may be a useful approach to glean new insight into biological processes. The combination of multiple datasets, as illustrated here, can further improve the functional inferences that typically are generated via a single type of datasets. In addition, our data extend the typical co-expression approach to allow deeper insight into cell wall biology in rice.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-596) contains supplementary material, which is available to authorized users.

Cell wall-related proteins of unknown function: missing links in plant cell wall development

Lignocellulosic biomass is an important feedstock for the pulp and paper industry as well as emerging biofuel and biomaterial industries. However, the recalcitrance of the secondary cell wall to chemical or enzymatic degradation remains a major hurdle for efficient extraction of economically important biopolymers such as cellulose. It has been estimated that approximately 10-15% of about 27,000 protein-coding genes in the Arabidopsis genome are dedicated to cell wall development; however, only about 130 Arabidopsis genes thus far have experimental evidence validating cell wall function. While many genes have been implicated through co-expression analysis with known genes, a large number are broadly classified as proteins of unknown function (PUFs). Recently the functionality of some of these unknown proteins in cell wall development has been revealed using reverse genetic approaches. Given the large number of cell wall-related PUFs, how do we approach and subsequently prioritize the investigation of such unknown genes that may be essential to or influence plant cell wall development and structure? Here, we address the aforementioned question in two parts; we first identify the different kinds of PUFs based on known and predicted features such as protein domains. Knowledge of inherent features of PUFs may allow for functional inference and a concomitant link to biological context. Secondly, we discuss omics-based technologies and approaches that are helping identify and prioritize cell wall-related PUFs by functional association. In this way, hypothesis-driven experiments can be designed for functional elucidation of many proteins that remain missing links in our understanding of plant cell wall biosynthesis.

ComPlEx: conservation and divergence of co-expression networks in A. thaliana, Populus and O. sativa

Background

Divergence in gene regulation has emerged as a key mechanism underlying species differentiation. Comparative analysis of co-expression networks across species can reveal conservation and divergence in the regulation of genes.

Results

We inferred co-expression networks of A. thaliana, Populus spp. and O. sativa using state-of-the-art methods based on mutual information and context likelihood of relatedness, and conducted a comprehensive comparison of these networks across a range of co-expression thresholds. In addition to quantifying gene-gene link and network neighbourhood conservation, we also applied recent advancements in network analysis to do cross-species comparisons of network properties such as scale free characteristics and gene centrality as well as network motifs. We found that in all species the networks emerged as scale free only above a certain co-expression threshold, and that the high-centrality genes upholding this organization tended to be conserved. Network motifs, in particular the feed-forward loop, were found to be significantly enriched in specific functional subnetworks but where much less conserved across species than gene centrality. Although individual gene-gene co-expression had massively diverged, up to ~80% of the genes still had a significantly conserved network neighbourhood. For genes with multiple predicted orthologs, about half had one ortholog with conserved regulation and another ortholog with diverged or non-conserved regulation. Furthermore, the most sequence similar ortholog was not the one with the most conserved gene regulation in over half of the cases.

Conclusions

We have provided a comprehensive analysis of gene regulation evolution in plants and built a web tool for Comparative analysis of Plant co-Expression networks (ComPlEx, http://complex.plantgenie.org/). The tool can be particularly useful for identifying the ortholog with the most conserved regulation among several sequence-similar alternatives and can thus be of practical importance in e.g. finding candidate genes for perturbation experiments.

Elucidating gene function and function evolution through comparison of co-expression networks of plants

The analysis of gene expression data has shown that transcriptionally coordinated (co-expressed) genes are often functionally related, enabling scientists to use expression data in gene function prediction. This Focused Review discusses our original paper (Large-scale co-expression approach to dissect secondary cell wall formation across plant species, Frontiers in Plant Science 2:23). In this paper we applied cross-species analysis to co-expression networks of genes involved in cellulose biosynthesis. We showed that the co-expression networks from different species are highly similar, indicating that whole biological pathways are conserved across species. This finding has two important implications. First, the analysis can transfer gene function annotation from well-studied plants, such as Arabidopsis, to other, uncharacterized plant species. As the analysis finds genes that have similar sequence and similar expression pattern across different organisms, functionally equivalent genes can be identified. Second, since co-expression analyses are often noisy, a comparative analysis should have higher performance, as parts of co-expression networks that are conserved are more likely to be functionally relevant. In this Focused Review, we outline the comparative analysis done in the original paper and comment on the recent advances and approaches that allow comparative analyses of co-function networks. We hypothesize that in comparison to simple co-expression analysis, comparative analysis would yield more accurate gene function predictions. Finally, by combining comparative analysis with genomic information of green plants, we propose a possible composition of cellulose biosynthesis machinery during earlier stages of plant evolution.

Genome-wide association mapping for wood characteristics in Populus identifies an array of candidate single nucleotide polymorphisms

Establishing links between phenotypes and molecular variants is of central importance to accelerate genetic improvement of economically important plant species. Our work represents the first genome-wide association study to the inherently complex and currently poorly understood genetic architecture of industrially relevant wood traits. Here, we employed an Illumina Infinium 34K single nucleotide polymorphism (SNP) genotyping array that generated 29,233 high-quality SNPs in c. 3500 broad-based candidate genes within a population of 334 unrelated Populus trichocarpa individuals to establish genome-wide associations. The analysis revealed 141 significant SNPs (α ≤ 0.05) associated with 16 wood chemistry/ultrastructure traits, individually explaining 3-7% of the phenotypic variance. A large set of associations (41% of all hits) occurred in candidate genes preselected for their suggested a priori involvement with secondary growth. For example, an allelic variant in the FRA8 ortholog explained 21% of the total genetic variance in fiber length, when the trait's heritability estimate was considered. The remaining associations identified SNPs in genes not previously implicated in wood or secondary wall formation. Our findings provide unique insights into wood trait architecture and support efforts for population improvement based on desirable allelic variants.

Network analysis for gene discovery in plant-specialized metabolism

Recent omics technologies provide information on multiple components of biological networks. Web-based data mining tools are continuously being developed. Because genes involved in specialized (secondary) metabolism are often co-ordinately regulated at the transcriptional level, a number of gene discovery studies have been successfully conducted using network analysis, especially by integrating gene co-expression network analysis and metabolomic investigation. In addition, next-generation sequencing technologies are currently utilized in functional genomics investigations of Arabidopsis and non-model plant species including medicinal plants. Systems-based approaches are expected to gain importance in medicinal plant research. This review discussed network analysis in Arabidopsis and gene discovery in plant-specialized metabolism in non-model plants.

Identification of molecular processes needed for vascular formation through transcriptome analysis of different vascular systems

Background

Vascular system formation has been studied through molecular and genetic approaches in Arabidopsis, a herbaceous dicot that is used as a model system. Different vascular systems have developed in other plants such as crops and trees. Uncovering shared mechanisms underlying vascular development by transcriptome analysis of different vascular systems may help to transfer knowledge acquired from Arabidopsis to other economically important species.

Results

Conserved vascular genes and biological processes fundamental to vascular development were explored across various plants. Through comparative transcriptome analysis, 226 genes from Arabidopsis, 217 genes from poplar and 281 genes from rice were identified as constituting 107 conserved vascular gene groups. These gene groups are expressed mainly in vascular tissues and form a complex coexpression network with multiple functional connections. To date, only half of the groups have been experimentally investigated. The conserved vascular gene groups were classified into 9 essential processes for vascular development. 18 groups (17%) lack of annotations were classified as having unknown functions.

Conclusion

The study provides a map of fundamental biological processes conserved across different vascular systems. It identifies gaps in the experimental investigation of pathways active in vascular formation, which if explored, could lead to a more complete understanding of vascular development.

Engineering secondary cell wall deposition in plants

Lignocellulosic biomass was used for thousands of years as animal feed and is now considered a great sugar source for biofuels production. It is composed mostly of secondary cell walls built with polysaccharide polymers that are embedded in lignin to reinforce the cell wall structure and maintain its integrity. Lignin is the primary material responsible for biomass recalcitrance to enzymatic hydrolysis. During plant development, deep reductions of lignin cause growth defects and often correlate with the loss of vessel integrity that adversely affects water and nutrient transport in plants. The work presented here describes a new approach to decrease lignin content while preventing vessel collapse and introduces a new strategy to boost transcription factor expression in native tissues. We used synthetic biology tools in Arabidopsis to rewire the secondary cell network by changing promoter-coding sequence associations. The result was a reduction in lignin and an increase in polysaccharide depositions in fibre cells. The promoter of a key lignin gene, C4H, was replaced by the vessel-specific promoter of transcription factor VND6. This rewired lignin biosynthesis specifically for vessel formation while disconnecting C4H expression from the fibre regulatory network. Secondly, the promoter of the IRX8 gene, secondary cell wall glycosyltransferase, was used to express a new copy of the fibre transcription factor NST1, and as the IRX8 promoter is induced by NST1, this also created an artificial positive feedback loop (APFL). The combination of strategies-lignin rewiring with APFL insertion-enhances polysaccharide deposition in stems without over-lignifying them, resulting in higher sugar yields after enzymatic hydrolysis.

Regulation of secondary wall synthesis and cell death by NAC transcription factors in the monocot Brachypodium distachyon

In several dicotyledonous species, NAC transcription factors act as master switches capable of turning on programmes of secondary cell-wall synthesis and cell death. This work used an oestradiol-inducible system to overexpress the NAC transcription factor BdSWN5 in the monocot model Brachypodium distachyon. This resulted in ectopic secondary cell-wall formation in both roots and shoots. Some of the genes upregulated in the process were a secondary cell-wall cellulose synthase (BdCESA4), a xylem-specific protease (BdXCP1) and an orthologue of AtMYB46 (BdMYB1). While activation of BdMYB1 may not be direct, this study showed that BdSWN5 is capable of transactivating the BdXCP1 promoter through two conserved binding sites. In the course of Brachypodium development, the BdXCP1 promoter was observed to be active in all types of differentiating tracheary elements. Together, these results suggest that Brachypodium SWNs can act as switches that turn on secondary cell-wall synthesis and programmed cell death.

Global identification of multiple OsGH9 family members and their involvement in cellulose crystallinity modification in rice

Plant glycoside hydrolase family 9 (GH9) comprises typical endo-β-1,4-glucanase (EGases, EC3.2.1.4). Although GH9A (KORRIGAN) family genes have been reported to be involved in cellulose biosynthesis in plants, much remains unknown about other GH9 subclasses. In this study, we observed a global gene co-expression profiling and conducted a correlation analysis between OsGH9 and OsCESA among 66 tissues covering most periods of life cycles in 2 rice varieties. Our results showed that OsGH9A3 and B5 possessed an extremely high co-expression with OsCESA1, 3, and 8 typical for cellulose biosynthesis in rice. Using two distinct rice non-GH9 mutants and wild type, we performed integrative analysis of gene expression level by qRT-PCR, cellulase activities in situ and in vitro, and lignocellulose crystallinity index (CrI) in four internodes of stem tissues. For the first time, OsGH9B1, 3, and 16 were characterized with the potential role in lignocellulose crystallinity alteration in rice, whereas OsGH9A3 and B5 were suggested for cellulose biosynthesis. In addition, phylogenetic analysis and gene co-expression comparison revealed GH9 function similarity in Arabidopsis and rice. Hence, the data can provide insights into GH9 function in plants and offer the potential strategy for genetic manipulation of plant cell wall using the five aforementioned novel OsGH9 genes.

Identification of novel transcription factors regulating secondary cell wall formation in Arabidopsis

The presence of lignin in secondary cell walls (SCW) is a major factor preventing hydrolytic enzymes from gaining access to cellulose, thereby limiting the saccharification potential of plant biomass. To understand how lignification is regulated is a prerequisite for selecting plant biomass better adapted to bioethanol production. Because transcriptional regulation is a major mechanism controlling the expression of genes involved in lignin biosynthesis, our aim was to identify novel transcription factors (TFs) dictating lignin profiles in the model plant Arabidopsis. To this end, we have developed a post-genomic approach by combining four independent in-house SCW-related transcriptome datasets obtained from (1) the fiber cell wall-deficient wat1 Arabidopsis mutant, (2) Arabidopsis lines over-expressing either the master regulatory activator EgMYB2 or (3) the repressor EgMYB1 and finally (4) Arabidopsis orthologs of Eucalyptus xylem-expressed genes. This allowed us to identify 502 up- or down-regulated TFs. We preferentially selected those present in more than one dataset and further analyzed their in silico expression patterns as an additional selection criteria. This selection process led to 80 candidates. Notably, 16 of them were already proven to regulate SCW formation, thereby validating the overall strategy. Then, we phenotyped 43 corresponding mutant lines focusing on histological observations of xylem and interfascicular fibers. This phenotypic screen revealed six mutant lines exhibiting altered lignification patterns. Two of them [Bel-like HomeoBox6 (blh6) and a zinc finger TF] presented hypolignified SCW. Three others (myb52, myb-like TF, hb5) showed hyperlignified SCW whereas the last one (hb15) showed ectopic lignification. In addition, our meta-analyses highlighted a reservoir of new potential regulators adding to the gene network regulating SCW but also opening new avenues to ultimately improve SCW composition for biofuel production.

Network analysis: tackling complex data to study plant metabolism

Incomplete knowledge of biochemical pathways makes the holistic description of plant metabolism a non-trivial undertaking. Sensitive analytical platforms, which are capable of accurately quantifying the levels of the various molecular entities of the cell, can assist in tackling this task. However, the ever-increasing amount of high-throughput data, often from multiple technologies, requires significant computational efforts for integrative analysis. Here we introduce the application of network analysis to study plant metabolism and describe the construction and analysis of correlation-based networks from (time-resolved) metabolomics data. By investigating the interactions between metabolites, network analysis can help to interpret complex datasets through the identification of key network components. The relationship between structural and biological roles of network components can be evaluated and employed to aid metabolic engineering.

Metabolic module mining based on Independent Component Analysis in Arabidopsis thaliana

Independent Component Analysis (ICA) has been introduced as one of the useful tools for gene-functional discovery in animals. However, this approach has been poorly utilized in the plant sciences. In the present study, we have exploited ICA combined with pathway enrichment analysis to address the statistical challenges associated with genome-wide analysis in plant system. To generate an Arabidopsis metabolic platform, we collected 4,373 Affy-metrix ATH1 microarray datasets. Out of the 3,232 metabolic genes and transcription factors, 99.47% of these genes were identified in at least one component, indicating the coverage of most of the metabolic pathways by the components. During the metabolic pathway enrichment analysis, we found components that indicate an independent regulation between the isoprenoid biosynthesis pathways. We also utilized this analysis tool to investigate some transcription factors involved in secondary cell wall biogenesis. This approach has identified remarkably more transcription factors compared to previously reported analysis tools. A website providing user-friendly searching and downloading of the entire dataset analyzed by ICA is available at http://kimjy.gnu.ac.kr/ICA.files/slide0002.htm . ICA combined with pathway enrichment analysis might provide a powerful approach for the extraction of the components responsible for a biological process of interest in plant systems.

4-O-methylation of glucuronic acid in Arabidopsis glucuronoxylan is catalyzed by a domain of unknown function family 579 protein

The hemicellulose 4-O-methyl glucuronoxylan is one of the principle components present in the secondary cell walls of eudicotyledonous plants. However, the biochemical mechanisms leading to the formation of this polysaccharide and the effects of modulating its structure on the physical properties of the cell wall are poorly understood. We have identified and functionally characterized an Arabidopsis glucuronoxylan methyltransferase (GXMT) that catalyzes 4-O-methylation of the glucuronic acid substituents of this polysaccharide. AtGXMT1, which was previously classified as a domain of unknown function (DUF) 579 protein, specifically transfers the methyl group from S-adenosyl-L-methionine to O-4 of α-D-glucopyranosyluronic acid residues that are linked to O-2 of the xylan backbone. Biochemical characterization of the recombinant enzyme indicates that GXMT1 is localized in the Golgi apparatus and requires Co(2+) for optimal activity in vitro. Plants lacking GXMT1 synthesize glucuronoxylan in which the degree of 4-O-methylation is reduced by 75%. This result is correlated to a change in lignin monomer composition and an increase in glucuronoxylan release during hydrothermal treatment of secondary cell walls. We propose that the DUF579 proteins constitute a previously undescribed family of cation-dependent, polysaccharide-specific O-methyl-transferases. This knowledge provides new opportunities to selectively manipulate polysaccharide O-methylation and extends the portfolio of structural targets that can be modified either alone or in combination to modulate biopolymer interactions in the plant cell wall.

Genome-scale identification of cell-wall related genes in Arabidopsis based on co-expression network analysis

BACKGROUND

Identification of the novel genes relevant to plant cell-wall (PCW) synthesis represents a highly important and challenging problem. Although substantial efforts have been invested into studying this problem, the vast majority of the PCW related genes remain unknown.

RESULTS

Here we present a computational study focused on identification of the novel PCW genes in Arabidopsis based on the co-expression analyses of transcriptomic data collected under 351 conditions, using a bi-clustering technique. Our analysis identified 217 highly co-expressed gene clusters (modules) under some experimental conditions, each containing at least one gene annotated as PCW related according to the Purdue Cell Wall Gene Families database. These co-expression modules cover 349 known/annotated PCW genes and 2,438 new candidates. For each candidate gene, we annotated the specific PCW synthesis stages in which it is involved and predicted the detailed function. In addition, for the co-expressed genes in each module, we predicted and analyzed their cis regulatory motifs in the promoters using our motif discovery pipeline, providing strong evidence that the genes in each co-expression module are transcriptionally co-regulated. From the all co-expression modules, we infer that 108 modules are related to four major PCW synthesis components, using three complementary methods.

CONCLUSIONS

We believe our approach and data presented here will be useful for further identification and characterization of PCW genes. All the predicted PCW genes, co-expression modules, motifs and their annotations are available at a web-based database: http://csbl.bmb.uga.edu/publications/materials/shanwang/CWRPdb/index.html.

Transcriptional Regulation of Grass Secondary Cell Wall Biosynthesis: Playing Catch-Up with Arabidopsis thaliana

Secondary cell wall synthesis occurs in specialized cell types following completion of cell enlargement. By virtue of mechanical strength provided by a wall thickened with cellulose, hemicelluloses, and lignin, these cells can function as water-conducting vessels and provide structural support. Several transcription factor families regulate genes encoding wall synthesis enzymes. Certain NAC and MYB proteins directly bind to the SNBE and AC elements upstream of structural genes and other transcription factors. The most detailed model of this regulatory network is established predominantly for a eudicot, Arabidopsis thaliana. In grasses, both the patterning and the composition of secondary cell walls are distinct from that of eudicots. These differences suggest transcriptional regulation is similarly distinct. Putative rice and maize orthologs of several eudicot cell wall regulators genetically complement mutants of A. thaliana or result in wall defects when constitutively overexpressed; nevertheless, aside from a maize, ZmMYB31, and a switchgrass protein, PvMYB4, function has not been tested in a grass. Similar to the seminal work conducted in A. thaliana, gene expression profiling in maize, rice, and other grasses implicates additional genes as regulators. Characterization of these genes will continue to elucidate the relationship between the transcription regulatory networks of eudicots and grasses.

Co-expression and co-responses: within and beyond transcription

Whole genome sequencing, the relative ease of transcript profiling by the use of microarrays and latterly RNA sequencing approaches have facilitated the capture of vast amounts of transcript data. However, despite the enormous progress made in gene annotation a substantial proportion of genes remain to be annotated at the functional level. Considerable progress has, however, been made by searching for transcriptional coordination between genes of known function and non-annotated genes on the premise that such co-expressed genes tend to be functionally related. Here we review progress made following this approach as well as its expansion to include phenotypic information from other levels of cellular organization such as proteomic and metabolomic data as well as physiological and developmental phenotypes.

Co-expression of cell-wall related genes: new tools and insights

Global transcript analyses based on publicly available microarray dataset have revealed that genes with similar function tend to be transcriptionally coordinated. Indeed, many genes involved in the formation of cellulose, hemicelluloses, and lignin have been identified using co-expression approaches in Arabidopsis. To facilitate these transcript analyses, several web-based tools have been developed that allow researchers to investigate co-expression relationships of their gene(s) of interest. In addition, several tools now also provide the possibility of comparative transcriptional analyses across species, which potentially increases the predictive power. In this short review, we describe recent developments and updates of plant-related co-expression tools, and summarize studies that have successfully used expression profiling in cell wall research. Finally, we illustrate the value of comparative co-expression relationships across species using genes involved in lignin biosynthesis.