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

Disruption of a DUF247 Containing Protein Alters Cell Wall Polysaccharides and Reduces Growth in Arabidopsis

Plant cell wall biosynthesis is a complex process that requires proteins and enzymes from glycan synthesis to wall assembly. We show that disruption of At3g50120 (DUF247-1), a member of the DUF247 multigene family containing 28 genes in Arabidopsis, results in alterations to the structure and composition of cell wall polysaccharides and reduced growth and plant size. An ELISA using cell wall antibodies shows that the mutants also exhibit ~50% reductions in xyloglucan (XyG), glucuronoxylan (GX) and heteromannan (HM) epitopes in the NaOH fraction and ~50% increases in homogalacturonan (HG) epitopes in the CDTA fraction. Furthermore, the polymer sizes of XyGs and GXs are reduced with concomitant increases in short-chain polymers, while those of HGs and mHGs are slightly increased. Complementation using 35S:DUF247-1 partially recovers the XyG and HG content, but not those of GX and HM, suggesting that DUF247-1 is more closely associated with XyGs and HGs. DUF247-1 is expressed throughout Arabidopsis, particularly in vascular and developing tissues, and its disruption affects the expression of other gene members, indicating a regulatory control role within the gene family. Our results demonstrate that DUF247-1 is required for normal cell wall composition and structure and Arabidopsis growth.

A Survey of the Transcriptomic Resources in Durum Wheat: Stress Responses, Data Integration and Exploitation

Durum wheat (Triticum turgidum subsp. durum (Desf.) Husn.) is an allotetraploid cereal crop of worldwide importance, given its use for making pasta, couscous, and bulgur. Under climate change scenarios, abiotic (e.g., high and low temperatures, salinity, drought) and biotic (mainly exemplified by fungal pathogens) stresses represent a significant limit for durum cultivation because they can severely affect yield and grain quality. The advent of next-generation sequencing technologies has brought a huge development in transcriptomic resources with many relevant datasets now available for durum wheat, at various anatomical levels, also focusing on phenological phases and environmental conditions. In this review, we cover all the transcriptomic resources generated on durum wheat to date and focus on the corresponding scientific insights gained into abiotic and biotic stress responses. We describe relevant databases, tools and approaches, including connections with other "omics" that could assist data integration for candidate gene discovery for bio-agronomical traits. The biological knowledge summarized here will ultimately help in accelerating durum wheat breeding.

Toward kingdom-wide analyses of gene expression

Gene expression data for Archaeplastida are accumulating exponentially, with more than 300 000 RNA-sequencing (RNA-seq) experiments available for hundreds of species. The gene expression data stem from thousands of experiments that capture gene expression in various organs, tissues, cell types, (a)biotic perturbations, and genotypes. Advances in software tools make it possible to process all these data in a matter of weeks on modern office computers, giving us the possibility to study gene expression in a kingdom-wide manner for the first time. We discuss how the expression data can be accessed and processed and outline analyses that take advantage of cross-species analyses, allowing us to generate powerful and robust hypotheses about gene function and evolution.

Identification of growth regulators using cross-species network analysis in plants

With the need to increase plant productivity, one of the challenges plant scientists are facing is to identify genes that play a role in beneficial plant traits. Moreover, even when such genes are found, it is generally not trivial to transfer this knowledge about gene function across species to identify functional orthologs. Here, we focused on the leaf to study plant growth. First, we built leaf growth transcriptional networks in Arabidopsis (Arabidopsis thaliana), maize (Zea mays), and aspen (Populus tremula). Next, known growth regulators, here defined as genes that when mutated or ectopically expressed alter plant growth, together with cross-species conserved networks, were used as guides to predict novel Arabidopsis growth regulators. Using an in-depth literature screening, 34 out of 100 top predicted growth regulators were confirmed to affect leaf phenotype when mutated or overexpressed and thus represent novel potential growth regulators. Globally, these growth regulators were involved in cell cycle, plant defense responses, gibberellin, auxin, and brassinosteroid signaling. Phenotypic characterization of loss-of-function lines confirmed two predicted growth regulators to be involved in leaf growth (NPF6.4 and LATE MERISTEM IDENTITY2). In conclusion, the presented network approach offers an integrative cross-species strategy to identify genes involved in plant growth and development.

Multi-Omics Approaches and Resources for Systems-Level Gene Function Prediction in the Plant Kingdom

In higher plants, the complexity of a system and the components within and among species are rapidly dissected by omics technologies. Multi-omics datasets are integrated to infer and enable a comprehensive understanding of the life processes of organisms of interest. Further, growing open-source datasets coupled with the emergence of high-performance computing and development of computational tools for biological sciences have assisted in silico functional prediction of unknown genes, proteins and metabolites, otherwise known as uncharacterized. The systems biology approach includes data collection and filtration, system modelling, experimentation and the establishment of new hypotheses for experimental validation. Informatics technologies add meaningful sense to the output generated by complex bioinformatics algorithms, which are now freely available in a user-friendly graphical user interface. These resources accentuate gene function prediction at a relatively minimal cost and effort. Herein, we present a comprehensive view of relevant approaches available for system-level gene function prediction in the plant kingdom. Together, the most recent applications and sought-after principles for gene mining are discussed to benefit the plant research community. A realistic tabulation of plant genomic resources is included for a less laborious and accurate candidate gene discovery in basic plant research and improvement strategies.

CAGEs are Golgi-localized GT31 enzymes involved in cellulose biosynthesis in Arabidopsis

Cellulose is the main structural component in the plant cell walls. We show that two glycosyltransferase family 31 (GT31) enzymes of Arabidopsis thaliana, here named cellulose synthesis associated glycosyltransferases 1 and 2 (CAGE1 and 2), influence both primary and secondary cell wall cellulose biosynthesis. cage1cage2 mutants show primary cell wall defects manifesting as impaired growth and cell expansion in seedlings and etiolated hypocotyls, along with secondary cell wall defects, apparent as collapsed xylem vessels and reduced xylem wall thickness in the inflorescence stem. Single and double cage mutants also show increased sensitivity to the cellulose biosynthesis inhibitor isoxaben. The cage1cage2 phenotypes were associated with an approximately 30% reduction in cellulose content, an approximately 50% reduction in secondary cell wall CELLULOSE SYNTHASE (CESA) protein levels in stems and reduced cellulose biosynthesis rate in seedlings. CESA transcript levels were not significantly altered in cage1cage2 mutants, suggesting that the reduction in CESA levels was caused by a post-transcriptional mechanism. Both CAGE1 and 2 localize to the Golgi apparatus and are predicted to synthesize β-1,3-galactans on arabinogalactan proteins. In line with this, the cage1cage2 mutants exhibit reduced levels of β-Yariv binding to arabinogalactan protein linked β-1,3-galactan. This leads us to hypothesize that defects in arabinogalactan biosynthesis underlie the cellulose deficiency of the mutants.

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.

FASCICLIN-LIKE 18 Is a New Player Regulating Root Elongation in Arabidopsis thaliana

The plasticity of root development represents a key trait that enables plants to adapt to diverse environmental cues. The pattern of cell wall deposition, alongside other parameters, affects the extent, and direction of root growth. In this study, we report that FASCICLIN-LIKE ARABINOGALACTAN PROTEIN 18 (FLA18) plays a role during root elongation in Arabidopsis thaliana. Using root-specific co-expression analysis, we identified FLA18 to be co-expressed with a sub-set of genes required for root elongation. FLA18 encodes for a putative extra-cellular arabinogalactan protein from the FLA-gene family. Two independent T-DNA insertion lines, named fla18-1 and fla18-2, display short and swollen lateral roots (LRs) when grown on sensitizing condition of high-sucrose containing medium. Unlike fla4/salt overly sensitive 5 (sos5), previously shown to display short and swollen primary root (PR) and LRs under these conditions, the PR of the fla18 mutants is slightly longer compared to the wild-type. Overexpression of the FLA18 CDS complemented the fla18 root phenotype. Genetic interaction between either of the fla18 alleles and sos5 reveals a more severe perturbation of anisotropic growth in both PR and LRs, as compared to the single mutants and the wild-type under restrictive conditions of high sucrose or high-salt containing medium. Additionally, under salt-stress conditions, fla18sos5 had a small, chlorotic shoot phenotype, that was not observed in any of the single mutants or the wild type. As previously shown for sos5, the fla18-1 and fla18-1sos5 root-elongation phenotype is suppressed by abscisic acid (ABA) and display hypersensitivity to the ABA synthesis inhibitor, Fluridon. Last, similar to other cell wall mutants, fla18 root elongation is hypersensitive to the cellulose synthase inhibitor, Isoxaben. Altogether, the presented data assign a new role for FLA18 in the regulation of root elongation. Future studies of the unique vs. redundant roles of FLA proteins during root elongation is anticipated to shed a new light on the regulation of root architecture during plant adaptation to different growth conditions.

Fungi.guru: Comparative genomic and transcriptomic resource for the fungi kingdom

The fungi kingdom is composed of eukaryotic heterotrophs, which are responsible for balancing the ecosystem and play a major role as decomposers. They also produce a vast diversity of secondary metabolites, which have antibiotic or pharmacological properties. However, our lack of knowledge of gene function in fungi precludes us from tailoring them to our needs and tapping into their metabolic diversity. To help remedy this, we gathered genomic and gene expression data of 19 most widely-researched fungi to build an online tool, fungi.guru, which contains tools for cross-species identification of conserved pathways, functional gene modules, and gene families. We exemplify how our tool can elucidate the molecular function, biological process and cellular component of genes involved in various biological processes, by identifying a secondary metabolite pathway producing gliotoxin in Aspergillus fumigatus, the catabolic pathway of cellulose in Coprinopsis cinerea and the conserved DNA replication pathway in Fusarium graminearum and Pyricularia oryzae. The tool is available at www.fungi.guru.

LSTrAP-Crowd: prediction of novel components of bacterial ribosomes with crowd-sourced analysis of RNA sequencing data

BACKGROUND

Bacterial resistance to antibiotics is a growing health problem that is projected to cause more deaths than cancer by 2050. Consequently, novel antibiotics are urgently needed. Since more than half of the available antibiotics target the structurally conserved bacterial ribosomes, factors involved in protein synthesis are thus prime targets for the development of novel antibiotics. However, experimental identification of these potential antibiotic target proteins can be labor-intensive and challenging, as these proteins are likely to be poorly characterized and specific to few bacteria. Here, we use a bioinformatics approach to identify novel components of protein synthesis.

RESULTS

In order to identify these novel proteins, we established a Large-Scale Transcriptomic Analysis Pipeline in Crowd (LSTrAP-Crowd), where 285 individuals processed 26 terabytes of RNA-sequencing data of the 17 most notorious bacterial pathogens. In total, the crowd processed 26,269 RNA-seq experiments and used the data to construct gene co-expression networks, which were used to identify more than a hundred uncharacterized genes that were transcriptionally associated with protein synthesis. We provide the identity of these genes together with the processed gene expression data.

CONCLUSIONS

We identified genes related to protein synthesis in common bacterial pathogens and thus provide a resource of potential antibiotic development targets for experimental validation. The data can be used to explore additional vulnerabilities of bacteria, while our approach demonstrates how the processing of gene expression data can be easily crowd-sourced.

Flux balance analysis of metabolism during growth by osmotic cell expansion and its application to tomato fruits

Cell expansion is a significant contributor to organ growth and is driven by the accumulation of osmolytes to increase cell turgor pressure. Metabolic modelling has the potential to provide insights into the processes that underpin osmolyte synthesis and transport, but the main computational approach for predicting metabolic network fluxes, flux balance analysis, often uses biomass composition as the main output constraint and ignores potential changes in cell volume. Here we present growth-by-osmotic-expansion flux balance analysis (GrOE-FBA), a framework that accounts for both the metabolic and ionic contributions to the osmotica that drive cell expansion, as well as the synthesis of protein, cell wall and cell membrane components required for cell enlargement. Using GrOE-FBA, the metabolic fluxes in dividing and expanding cells were analysed, and the energetic costs for metabolite biosynthesis and accumulation in the two scenarios were found to be surprisingly similar. The expansion phase of tomato fruit growth was also modelled using a multiphase single-optimization GrOE-FBA model and this approach gave accurate predictions of the major metabolite levels throughout fruit development, as well as revealing a role for transitory starch accumulation in ensuring optimal fruit development.

LSTrAP-Cloud: A User-Friendly Cloud Computing Pipeline to Infer Coexpression Networks

As genomes become more and more available, gene function prediction presents itself as one of the major hurdles in our quest to extract meaningful information on the biological processes genes participate in. In order to facilitate gene function prediction, we show how our user-friendly pipeline, the Large-Scale Transcriptomic Analysis Pipeline in Cloud (LSTrAP-Cloud), can be useful in helping biologists make a shortlist of genes involved in a biological process that they might be interested in, by using a single gene of interest as bait. The LSTrAP-Cloud is based on Google Colaboratory, and provides user-friendly tools that process quality-control RNA sequencing data streamed from the European Nucleotide Archive. The LSTRAP-Cloud outputs a gene coexpression network that can be used to identify functionally related genes for any organism with a sequenced genome and publicly available RNA sequencing data. Here, we used the biosynthesis pathway of Nicotiana tabacum as a case study to demonstrate how enzymes, transporters, and transcription factors involved in the synthesis, transport, and regulation of nicotine can be identified using our pipeline.

Expression Atlas of Selaginella moellendorffii Provides Insights into the Evolution of Vasculature, Secondary Metabolism, and Roots

Selaginella moellendorffii is a representative of the lycophyte lineage that is studied to understand the evolution of land plant traits such as the vasculature, leaves, stems, roots, and secondary metabolism. However, only a few studies have investigated the expression and transcriptional coordination of Selaginella genes, precluding us from understanding the evolution of the transcriptional programs behind these traits. We present a gene expression atlas comprising all major organs, tissue types, and the diurnal gene expression profiles for S. moellendorffii We show that the transcriptional gene module responsible for the biosynthesis of lignocellulose evolved in the ancestor of vascular plants and pinpoint the duplication and subfunctionalization events that generated multiple gene modules involved in the biosynthesis of various cell wall types. We demonstrate how secondary metabolism is transcriptionally coordinated and integrated with other cellular pathways. Finally, we identify root-specific genes and show that the evolution of roots did not coincide with an increased appearance of gene families, suggesting that the development of new organs does not coincide with increased fixation of new gene functions. Our updated database at conekt.plant.tools represents a valuable resource for studying the evolution of genes, gene families, transcriptomes, and functional gene modules in the Archaeplastida kingdom.

Transcriptional Regulation of Sorghum Stem Composition: Key Players Identified Through Co-expression Gene Network and Comparative Genomics Analyses

Most sorghum biomass accumulates in stem secondary cell walls (SCW). As sorghum stems are used as raw materials for various purposes such as feed, energy and fiber reinforced polymers, identifying the genes responsible for SCW establishment is highly important. Taking advantage of studies performed in model species, most of the structural genes contributing at the molecular level to the SCW biosynthesis in sorghum have been proposed while their regulatory factors have mostly not been determined. Validation of the role of several MYB and NAC transcription factors in SCW regulation in Arabidopsis and a few other species has been provided. In this study, we contributed to the recent efforts made in grasses to uncover the mechanisms underlying SCW establishment. We reported updated phylogenies of NAC and MYB in 9 different species and exploited findings from other species to highlight candidate regulators of SCW in sorghum. We acquired expression data during sorghum internode development and used co-expression analyses to determine groups of co-expressed genes that are likely to be involved in SCW establishment. We were able to identify two groups of co-expressed genes presenting multiple evidences of involvement in SCW building. Gene enrichment analysis of MYB and NAC genes provided evidence that while NAC SECONDARY WALL THICKENING PROMOTING FACTOR NST genes and SECONDARY WALL-ASSOCIATED NAC DOMAIN PROTEIN gene functions appear to be conserved in sorghum, NAC master regulators of SCW in sorghum may not be as tissue compartmentalized as in Arabidopsis. We showed that for every homolog of the key SCW MYB in Arabidopsis, a similar role is expected for sorghum. In addition, we unveiled sorghum MYB and NAC that have not been identified to date as being involved in cell wall regulation. Although specific validation of the MYB and NAC genes uncovered in this study is needed, we provide a network of sorghum genes involved in SCW both at the structural and regulatory levels.

Gene expression analysis of Cyanophora paradoxa reveals conserved abiotic stress responses between basal algae and flowering plants

The glaucophyte Cyanophora paradoxa represents the most basal member of the kingdom Archaeplastida, but the function and expression of most of its genes are unknown. This information is needed to uncover how functional gene modules, that is groups of genes performing a given function, evolved in the plant kingdom. We have generated a gene expression atlas capturing responses of Cyanophora to various abiotic stresses. The data were included in the CoNekT-Plants database, enabling comparative transcriptomic analyses across two algae and six land plants. We demonstrate how the database can be used to study gene expression, co-expression networks and gene function in Cyanophora, and how conserved transcriptional programs can be identified. We identified gene modules involved in phycobilisome biosynthesis, response to high light and cell division. While we observed no correlation between the number of differentially expressed genes and the impact on growth of Cyanophora, we found that the response to stress involves a conserved, kingdom-wide transcriptional reprogramming, which is activated upon most stresses in algae and land plants. The Cyanophora stress gene expression atlas and the tools found in the https://conekt.plant.tools/ database thus provide a useful resource to reveal functionally related genes and stress responses in the plant kingdom.

Phenolic cross-links: building and de-constructing the plant cell wall

Covering: Up to 2019Phenolic cross-links and phenolic inter-unit linkages result from the oxidative coupling of two hydroxycinnamates or two molecules of tyrosine. Free dimers of hydroxycinnamates, lignans, play important roles in plant defence. Cross-linking of bound phenolics in the plant cell wall affects cell expansion, wall strength, digestibility, degradability, and pathogen resistance. Cross-links mediated by phenolic substituents are particularly important as they confer strength to the wall via the formation of new covalent bonds, and by excluding water from it. Four biopolymer classes are known to be involved in the formation of phenolic cross-links: lignins, extensins, glucuronoarabinoxylans, and side-chains of rhamnogalacturonan-I. Lignins and extensins are ubiquitous in streptophytes whereas aromatic substituents on xylan and pectic side-chains are commonly assumed to be particular features of Poales sensu lato and core Caryophyllales, respectively. Cross-linking of phenolic moieties proceeds via radical formation, is catalyzed by peroxidases and laccases, and involves monolignols, tyrosine in extensins, and ferulate esters on xylan and pectin. Ferulate substituents, on xylan in particular, are thought to be nucleation points for lignin polymerization and are, therefore, of paramount importance to wall architecture in grasses and for the development of technology for wall disassembly, e.g. for the use of grass biomass for production of 2^nd generation biofuels. This review summarizes current knowledge on the intra- and extracellular acylation of polysaccharides, and inter- and intra-molecular cross-linking of different constituents. Enzyme mediated lignan in vitro synthesis for pharmaceutical uses are covered as are industrial exploitation of mutant and transgenic approaches to control cell wall cross-linking.

Analysis of Changes in Plant Cell Wall Composition and Structure During Cold Acclimation

The cell wall has a crucial influence on the mechanical properties of plant cells. It therefore has a strong impact on the freezing behavior and very likely also the freezing tolerance of plants. However, not many studies have addressed the question how cell wall composition and structure impact plant freezing tolerance and cold acclimation. In this chapter, we describe a comprehensive workflow to extract total cell wall material from leaves of Arabidopsis thaliana and to separate this material into fractions enriched in crystalline cellulose, pectins, and hemicelluloses by sequential fractionation. We further describe methods for the analysis of chemical structure, monosaccharide composition, and cellulose and uronic acid contents in the total cell wall material and the fractions in response to cold acclimation. Structural properties of cell wall material are analyzed by attenuated total reflectance-Fourier-transform infrared spectrometry (ATR-FTIR) and monosaccharide composition by gas chromatography-mass spectrometry (GC-MS) after isolation of alditol acetate derivatives of the sugars.

Inferring biosynthetic and gene regulatory networks from Artemisia annua RNA sequencing data on a credit card-sized ARM computer

Prediction of gene function and gene regulatory networks is one of the most active topics in bioinformatics. The accumulation of publicly available gene expression data for hundreds of plant species, together with advances in bioinformatical methods and affordable computing, sets ingenuity as one of the major bottlenecks in understanding gene function and regulation. Here, we show how a credit card-sized computer retailing for <50 USD can be used to rapidly predict gene function and infer regulatory networks from RNA sequencing data. To achieve this, we constructed a bioinformatical pipeline that downloads and allows quality-control of RNA sequencing data; and generates a gene co-expression network that can reveal enzymes and transcription factors participating and controlling a given biosynthetic pathway. We exemplify this by first identifying genes and transcription factors involved in the biosynthesis of secondary cell wall in the plant Artemisia annua, the main natural source of the anti-malarial drug artemisinin. Networks were then used to dissect the artemisinin biosynthesis pathway, which suggest potential transcription factors regulating artemisinin biosynthesis. We provide the source code of our pipeline (https://github.com/mutwil/LSTrAP-Lite) and envision that the ubiquity of affordable computing, availability of biological data and increased bioinformatical training of biologists will transform the field of bioinformatics. This article is part of a Special Issue entitled: Transcriptional Profiles and Regulatory Gene Networks edited by Dr. Dr. Federico Manuel Giorgi and Dr. Shaun Mahony.

Co-expression networks for plant biology: why and how

Co-expression network analysis is one of the most powerful approaches for interpretation of large transcriptomic datasets. It enables characterization of modules of co-expressed genes that may share biological functional linkages. Such networks provide an initial way to explore functional associations from gene expression profiling and can be applied to various aspects of plant biology. This review presents the applications of co-expression network analysis in plant biology and addresses optimized strategies from the recent literature for performing co-expression analysis on plant biological systems. Additionally, we describe the combined interpretation of co-expression analysis with other genomic data to enhance the generation of biologically relevant information.

CoNekT: an open-source framework for comparative genomic and transcriptomic network analyses

The recent accumulation of gene expression data in the form of RNA sequencing creates unprecedented opportunities to study gene regulation and function. Furthermore, comparative analysis of the expression data from multiple species can elucidate which functional gene modules are conserved across species, allowing the study of the evolution of these modules. However, performing such comparative analyses on raw data is not feasible for many biologists. Here, we present CoNekT (Co-expression Network Toolkit), an open source web server, that contains user-friendly tools and interactive visualizations for comparative analyses of gene expression data and co-expression networks. These tools allow analysis and cross-species comparison of (i) gene expression profiles; (ii) co-expression networks; (iii) co-expressed clusters involved in specific biological processes; (iv) tissue-specific gene expression; and (v) expression profiles of gene families. To demonstrate these features, we constructed CoNekT-Plants for green alga, seed plants and flowering plants (Picea abies, Chlamydomonas reinhardtii, Vitis vinifera, Arabidopsis thaliana, Oryza sativa, Zea mays and Solanum lycopersicum) and thus provide a web-tool with the broadest available collection of plant phyla. CoNekT-Plants is freely available from http://conekt.plant.tools, while the CoNekT source code and documentation can be found at https://github.molgen.mpg.de/proost/CoNekT/.

PhytoNet: comparative co-expression network analyses across phytoplankton and land plants

Phytoplankton consists of autotrophic, photosynthesizing microorganisms that are a crucial component of freshwater and ocean ecosystems. However, despite being the major primary producers of organic compounds, accounting for half of the photosynthetic activity worldwide and serving as the entry point to the food chain, functions of most of the genes of the model phytoplankton organisms remain unknown. To remedy this, we have gathered publicly available expression data for one chlorophyte, one rhodophyte, one haptophyte, two heterokonts and four cyanobacteria and integrated it into our PlaNet (Plant Networks) database, which now allows mining gene expression profiles and identification of co-expressed genes of 19 species. We exemplify how the co-expressed gene networks can be used to reveal functionally related genes and how the comparative features of PhytoNet allow detection of conserved transcriptional programs between cyanobacteria, green algae, and land plants. Additionally, we illustrate how the database allows detection of duplicated transcriptional programs within an organism, as exemplified by two putative DNA repair programs within Chlamydomonas reinhardtii. PhytoNet is available from www.gene2function.de.

Metabolism of polysaccharides in dynamic middle lamellae during cotton fibre development

Evidence is presented that cotton fibre adhesion and middle lamella formation are preceded by cutin dilution and accompanied by rhamnogalacturonan-I metabolism. Cotton fibres are single cell structures that early in development adhere to one another via the cotton fibre middle lamella (CFML) to form a tissue-like structure. The CFML is disassembled around the time of initial secondary wall deposition, leading to fibre detachment. Observations of CFML in the light microscope have suggested that the development of the middle lamella is accompanied by substantial cell-wall metabolism, but it has remained an open question as to which processes mediate adherence and which lead to detachment. The mechanism of adherence and detachment were investigated here using glyco-microarrays probed with monoclonal antibodies, transcript profiling, and observations of fibre auto-digestion. The results suggest that adherence is brought about by cutin dilution, while the presence of relevant enzyme activities and the dynamics of rhamnogalacturonan-I side-chain accumulation and disappearance suggest that both attachment and detachment are accompanied by rhamnogalacturonan-I metabolism.

Gene regulatory networks for lignin biosynthesis in switchgrass (Panicum virgatum)

Cell wall recalcitrance is the major challenge to improving saccharification efficiency in converting lignocellulose into biofuels. However, information regarding the transcriptional regulation of secondary cell wall biogenesis remains poor in switchgrass (Panicum virgatum), which has been selected as a biofuel crop in the United States. In this study, we present a combination of computational and experimental approaches to develop gene regulatory networks for lignin formation in switchgrass. To screen transcription factors (TFs) involved in lignin biosynthesis, we developed a modified method to perform co-expression network analysis using 14 lignin biosynthesis genes as bait (target) genes. The switchgrass lignin co-expression network was further extended by adding 14 TFs identified in this study, and seven TFs identified in previous studies, as bait genes. Six TFs (PvMYB58/63, PvMYB42/85, PvMYB4, PvWRKY12, PvSND2 and PvSWN2) were targeted to generate overexpressing and/or down-regulated transgenic switchgrass lines. The alteration of lignin content, cell wall composition and/or plant growth in the transgenic plants supported the role of the TFs in controlling secondary wall formation. RNA-seq analysis of four of the transgenic switchgrass lines revealed downstream target genes of the secondary wall-related TFs and crosstalk with other biological pathways. In vitro transactivation assays further confirmed the regulation of specific lignin pathway genes by four of the TFs. Our meta-analysis provides a hierarchical network of TFs and their potential target genes for future manipulation of secondary cell wall formation for lignin modification in switchgrass.

Xylan in the Middle: Understanding Xylan Biosynthesis and Its Metabolic Dependencies Toward Improving Wood Fiber for Industrial Processing

Lignocellulosic biomass, encompassing cellulose, lignin and hemicellulose in plant secondary cell walls (SCWs), is the most abundant source of renewable materials on earth. Currently, fast-growing woody dicots such as Eucalyptus and Populus trees are major lignocellulosic (wood fiber) feedstocks for bioproducts such as pulp, paper, cellulose, textiles, bioplastics and other biomaterials. Processing wood for these products entails separating the biomass into its three main components as efficiently as possible without compromising yield. Glucuronoxylan (xylan), the main hemicellulose present in the SCWs of hardwood trees carries chemical modifications that are associated with SCW composition and ultrastructure, and affect the recalcitrance of woody biomass to industrial processing. In this review we highlight the importance of xylan properties for industrial wood fiber processing and how gaining a greater understanding of xylan biosynthesis, specifically xylan modification, could yield novel biotechnology approaches to reduce recalcitrance or introduce novel processing traits. Altering xylan modification patterns has recently become a focus of plant SCW studies due to early findings that altered modification patterns can yield beneficial biomass processing traits. Additionally, it has been noted that plants with altered xylan composition display metabolic differences linked to changes in precursor usage. We explore the possibility of using systems biology and systems genetics approaches to gain insight into the coordination of SCW formation with other interdependent biological processes. Acetyl-CoA, s-adenosylmethionine and nucleotide sugars are precursors needed for xylan modification, however, the pathways which produce metabolic pools during different stages of fiber cell wall formation still have to be identified and their co-regulation during SCW formation elucidated. The crucial dependence on precursor metabolism provides an opportunity to alter xylan modification patterns through metabolic engineering of one or more of these interdependent pathways. The complexity of xylan biosynthesis and modification is currently a stumbling point, but it may provide new avenues for woody biomass engineering that are not possible for other biopolymers.

Both cold and sub-zero acclimation induce cell wall modification and changes in the extracellular proteome in Arabidopsis thaliana

Cold acclimation (CA) leads to increased plant freezing tolerance during exposure to low, non-freezing temperatures as a result of many physiological, biochemical and molecular changes that have been extensively investigated. In addition, many plant species, such as Arabidopsis thaliana, respond to a subsequent exposure to mild, non-damaging freezing temperatures with an additional increase in freezing tolerance referred to as sub-zero acclimation (SZA). There is comparatively little information available about the molecular basis of SZA. However, previous transcriptomic studies indicated that cell wall modification may play an important role during SZA. Here we show that CA and SZA are accompanied by extensive changes in cell wall amount, composition and structure. While CA leads to a significant increase in cell wall amount, the relative proportions of pectin, hemicellulose and cellulose remained unaltered during both CA and SZA. However, both treatments resulted in more subtle changes in structure as determined by infrared spectroscopy and monosaccharide composition as determined by gas chromatography-mass spectrometry. These differences could be related through a proteomic approach to the accumulation of cell wall modifying enzymes such as pectin methylesterases, pectin methylesterase inhibitors and xyloglucan endotransglucosylases/hydrolases in the extracellular matrix.

Secondary Wall Regulating NACs Differentially Bind at the Promoter at a CELLULOSE SYNTHASE A4 Cis-eQTL

Arabidopsis thaliana CELLULOSE SYNTHASE A4/7/8 (CESA4/7/8) are three non-redundant subunits of the secondary cell wall cellulose synthase complex. Transcript abundance of these genes can vary among genotypes and expression quantitative trait loci (eQTL) were identified in a recombinant population of the accessions Bay-0 and Shahdara. Genetic mapping and analysis of the transcript levels of CESAs between two distinct near isogenic lines (NILs) confirmed a change in CESA4 expression that segregates within that interval. We sequenced the promoters and identified 16 polymorphisms differentiating CESA4^Sha and CESA4^Bay . In order to determine which of these SNPs could be responsible for this eQTL, we screened for transcription factor protein affinity with promoter fragments of CESA4^Bay, CESA4^Sha , and the reference genome CESA4^Col . The wall thickening activator proteins NAC SECONDARY WALL THICKENING PROMOTING FACTOR2 (NST2) and NST3 exhibited a decrease in binding with the CESA4^Sha promoter with a tracheary element-regulating cis-element (TERE) polymorphism. While NILs harboring the TERE polymorphisms exhibited significantly different CESA4 expression, cellulose crystallinity and cell wall thickness were indistinguishable. These results suggest that the TERE polymorphism resulted in differential transcription factor binding and CESA4 expression; yet A. thaliana is able to tolerate this transcriptional variability without compromising the structural elements of the plant, providing insight into the elasticity of gene regulation as it pertains to cell wall biosynthesis and regulation. We also explored available DNA affinity purification sequencing data to resolve a core binding site, C(G/T)TNNNNNNNA(A/C)G, for secondary wall NACs referred to as the VNS element.

A New Calmodulin-Binding Protein Expresses in the Context of Secondary Cell Wall Biosynthesis and Impacts Biomass Properties in Populus

A greater understanding of biosynthesis, signaling and regulatory pathways involved in determining stem growth and secondary cell wall chemistry is important for enabling pathway engineering and genetic optimization of biomass properties. The present study describes a new functional role of PdIQD10, a Populus gene belonging to the IQ67-Domain1 family of IQD genes, in impacting biomass formation and chemistry. Expression studies showed that PdIQD10 has enhanced expression in developing xylem and tension-stressed tissues in Populus deltoides. Molecular dynamics simulation and yeast two-hybrid interaction experiments suggest interactions with two calmodulin proteins, CaM247 and CaM014, supporting the sequence-predicted functional role of the PdIQD10 as a calmodulin-binding protein. PdIQD10 was found to interact with specific Populus isoforms of the Kinesin Light Chain protein family, shown previously to function as microtubule-guided, cargo binding and delivery proteins in Arabidopsis. Subcellular localization studies showed that PdIQD10 localizes in the nucleus and plasma membrane regions. Promoter-binding assays suggest that a known master transcriptional regulator of secondary cell wall biosynthesis (PdWND1B) may be upstream of an HD-ZIP III gene that is in turn upstream of PdIQD10 gene in the transcriptional network. RNAi-mediated downregulation of PdIQD10 expression resulted in plants with altered biomass properties including higher cellulose, wall glucose content and greater biomass quantity. These results present evidence in support of a new functional role for an IQD gene family member, PdIQD10, in secondary cell wall biosynthesis and biomass formation in Populus.

Insights into ABA-mediated regulation of guard cell primary metabolism revealed by systems biology approaches

Despite the fact that guard cell abscisic acid (ABA) signalling pathway is well documented, our understanding concerning how and to which extent ABA regulates guard cell metabolism remains fragmentary. Here we have adopted different systems approaches to investigate how ABA modulates guard cell central metabolism by providing genes that are possibly ABA-regulated. By using previous published Arabidopsis guard cell transcript profiling data, we carried out an extensive co-expression network analysis using ABA-related genes and those related to the metabolism and transport of sugars, starch and organic acids. Next, we investigated the presence of ABA responsive elements (ABRE) in the promoter of genes that are highly expressed in guard cells, responsive to ABA and co-expressed with ABA-related genes. Together, these analyses indicated that 44 genes are likely regulated by ABA and 8 of them are highly expressed in guard cells in both the presence and absence of ABA, including genes of the tricarboxylic acid cycle and those related to sucrose and hexose transport and metabolism. It seems likely that ABA may modulate both sucrose transport through guard cell plasma membrane and sucrose metabolism within guard cells. In this context, genes associated with sucrose synthase, sucrose phosphate synthase, trehalose-6-phosphate, invertase, UDP-glucose epimerase/pyrophosphorylase and different sugar transporters contain ABRE in their promoter and are thus possibly ABA regulated. Although validation experiments are required, our study highlights the importance of systems biology approaches to drive new hypothesis and to unravel genes and pathways that are regulated by ABA in guard cells.

Cellulose synthase 'class specific regions' are intrinsically disordered and functionally undifferentiated

Cellulose synthases (CESAs) are glycosyltransferases that catalyze formation of cellulose microfibrils in plant cell walls. Seed plant CESA isoforms cluster in six phylogenetic clades, whose non-interchangeable members play distinct roles within cellulose synthesis complexes (CSCs). A 'class specific region' (CSR), with higher sequence similarity within versus between functional CESA classes, has been suggested to contribute to specific activities or interactions of different isoforms. We investigated CESA isoform specificity in the moss, Physcomitrella patens (Hedw.) B. S. G. to gain evolutionary insights into CESA structure/function relationships. Like seed plants, P. patens has oligomeric rosette-type CSCs, but the PpCESAs diverged independently and form a separate CESA clade. We showed that P. patens has two functionally distinct CESAs classes, based on the ability to complement the gametophore-negative phenotype of a ppcesa5 knockout line. Thus, non-interchangeable CESA classes evolved separately in mosses and seed plants. However, testing of chimeric moss CESA genes for complementation demonstrated that functional class-specificity is not determined by the CSR. Sequence analysis and computational modeling showed that the CSR is intrinsically disordered and contains predicted molecular recognition features, consistent with a possible role in CESA oligomerization and explaining the evolution of class-specific sequences without selection for class-specific function.

Co-expression and Transcriptome Analysis of Marchantia polymorpha Transcription Factors Supports Class C ARFs as Independent Actors of an Ancient Auxin Regulatory Module

We performed differential gene expression (DGE) and co-expression analyses with genes encoding components of hormonal signaling pathways and the ∼400 annotated transcription factors (TFs) of M. polymorpha across multiple developmental stages of the life cycle. We identify a putative auxin-related co-expression module that has significant overlap with transcripts induced in auxin-treated tissues. Consistent with phylogenetic and functional studies, the class C ARF, MpARF3, is not part of the auxin-related co-expression module and instead is associated with transcripts enriched in gamete-producing gametangiophores. We analyze the Mparf3 and MpmiR160 mutant transcriptomes in the context of coexpression to suggest that MpARF3 may antagonize the reproductive transition via activating the MpMIR11671 and MpMIR529c precursors whose encoded microRNAs target SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) transcripts of MpSPL1 and MpSPL2. Both MpSPL genes are part of the MpARF3 co-expression group corroborating their functional significance. We provide evidence of the independence of MpARF3 from the auxin-signaling module and provide new testable hypotheses on the role of auxin-related genes in patterning meristems and differentiation events in liverworts.

Systems Identification and Characterization of Cell Wall Reassembly and Degradation Related Genes in Glycine max (L.) Merill, a Bioenergy Legume

Soybean is a promising biomass resource for generation of second-generation biofuels. Despite the utility of soybean cellulosic biomass and post-processing residues in biofuel generation, there is no comprehensive information available on cell wall loosening and degradation related gene families. In order to achieve enhanced lignocellulosic biomass with softened cell walls and reduced recalcitrance, it is important to identify genes involved in cell wall polymer loosening and degrading. Comprehensive genome-wide analysis of gene families involved in cell wall modifications is an efficient stratagem to find new candidate genes for soybean breeding for expanding biofuel industry. We report the identification of 505 genes distributed among 12 gene families related to cell wall loosening and degradation. 1262 tandem duplication events contributed towards expansion and diversification of studied gene families. We identified 687 Simple Sequence Repeat markers and 5 miRNA families distributed on 316 and 10 genes, respectively. Publically available microarray datasets were used to explore expression potential of identified genes in soybean plant developmental stages, 68 anatomical parts, abiotic and biotic stresses. Co-expression networks revealed transcriptional coordination of different gene families involved in cell wall loosening and degradation process.

Expression atlas and comparative coexpression network analyses reveal important genes involved in the formation of lignified cell wall in Brachypodium distachyon

While Brachypodium distachyon (Brachypodium) is an emerging model for grasses, no expression atlas or gene coexpression network is available. Such tools are of high importance to provide insights into the function of Brachypodium genes. We present a detailed Brachypodium expression atlas, capturing gene expression in its major organs at different developmental stages. The data were integrated into a large-scale coexpression database ( www.gene2function.de), enabling identification of duplicated pathways and conserved processes across 10 plant species, thus allowing genome-wide inference of gene function. We highlight the importance of the atlas and the platform through the identification of duplicated cell wall modules, and show that a lignin biosynthesis module is conserved across angiosperms. We identified and functionally characterised a putative ferulate 5-hydroxylase gene through overexpression of it in Brachypodium, which resulted in an increase in lignin syringyl units and reduced lignin content of mature stems, and led to improved saccharification of the stem biomass. Our Brachypodium expression atlas thus provides a powerful resource to reveal functionally related genes, which may advance our understanding of important biological processes in grasses.

Genome and transcriptome-wide analyses of cellulose synthase gene superfamily in soybean

The plant cellulose synthase gene superfamily belongs to the category of type-2 glycosyltransferases, and is involved in cellulose and hemicellulose biosynthesis. These enzymes are vital for maintaining cell-wall structural integrity throughout plant life. Here, we identified 78 putative cellulose synthases (CS) in the soybean genome. Phylogenetic analysis against 40 reference Arabidopsis CS genes clustered soybean CSs into seven major groups (CESA, CSL A, B, C, D, E and G), located on 19 chromosomes (except chromosome 18). Soybean CS expansion occurred in 66 duplication events. Additionally, we identified 95 simple sequence repeat makers related to 44 CSs. We next performed digital expression analysis using publically available datasets to understand potential CS functions in soybean. We found that CSs were highly expressed during soybean seed development, a pattern confirmed with an Affymatrix soybean IVT array and validated with RNA-seq profiles. Within CS groups, CESAs had higher relative expression than CSLs. Soybean CS models were designed based on maximum average RPKM values. Gene co-expression networks were developed to explore which CSs could work together in soybean. Finally, RT-PCR analysis confirmed the expression of 15 selected CSs during all four seed developmental stages.

Cellulosic fibres of flax recruit both primary and secondary cell wall cellulose synthases during deposition of thick tertiary cell walls and in the course of graviresponse

Cellulose synthesising complex consists of cellulose synthase (CESA) subunits encoded by a multigene family; different sets of CESA genes are known to be expressed during primary and secondary cell wall formation. We examined the expression of LusCESAs in flax (Linum usitatissimum L.) cellulosic fibres at various stages of development and in the course of graviresponse by means of RNA-Seq and quantitative PCR. Transcripts for both primary and secondary cell wall-related CESAs were abundant in fibres depositing highly cellulosic tertiary cell walls. Gravistimulation of flax plants temporally increased the abundance of CESA transcripts, specifically in phloem fibres located at the pulling stem side. Construction of coexpression networks for LusCESAs revealed that both primary and secondary cell wall-related CESAs were involved in the joint coexpression group in fibres depositing tertiary cell walls, as distinct from other tissues, where these genes were within separate groups. The obtained data suggest that fibres depositing tertiary cell walls have a specific mechanism of cellulose biosynthesis and a specific way of its regulation.

Phylogenomic analysis of gene co-expression networks reveals the evolution of functional modules

Molecular evolutionary studies correlate genomic and phylogenetic information with the emergence of new traits of organisms. These traits are, however, the consequence of dynamic gene networks composed of functional modules, which might not be captured by genomic analyses. Here, we established a method that combines large-scale genomic and phylogenetic data with gene co-expression networks to extensively study the evolutionary make-up of modules in the moss Physcomitrella patens, and in the angiosperms Arabidopsis thaliana and Oryza sativa (rice). We first show that younger genes are less annotated than older genes. By mapping genomic data onto the co-expression networks, we found that genes from the same evolutionary period tend to be connected, whereas old and young genes tend to be disconnected. Consequently, the analysis revealed modules that emerged at a specific time in plant evolution. To uncover the evolutionary relationships of the modules that are conserved across the plant kingdom, we added phylogenetic information that revealed duplication and speciation events on the module level. This combined analysis revealed an independent duplication of cell wall modules in bryophytes and angiosperms, suggesting a parallel evolution of cell wall pathways in land plants. We provide an online tool allowing plant researchers to perform these analyses at http://www.gene2function.de.

Beyond Genomics: Studying Evolution with Gene Coexpression Networks

Understanding how genomes change as organisms become more complex is a central question in evolution. Molecular evolutionary studies typically correlate the appearance of genes and gene families with the emergence of biological pathways and morphological features. While such approaches are of great importance to understand how organisms evolve, they are also limited, as functionally related genes work together in contexts of dynamic gene networks. Since functionally related genes are often transcriptionally coregulated, gene coexpression networks present a resource to study the evolution of biological pathways. In this opinion article, we discuss recent developments in this field and how coexpression analyses can be merged with existing genomic approaches to transfer functional knowledge between species to study the appearance or extension of pathways.

PlaNet: Comparative Co-Expression Network Analyses for Plants

Functional relations between genes can be represented as networks. These networks have been successfully used to infer gene function and to mediate transfer of functional knowledge between species. Transcriptionally coordinated or co-expressed genes tend to be functionally related, which combined with availability of transcriptomic data for multiple plant species make the co-expression networks a useful resource for the plant community. In this chapter, we describe PlaNet ( www.gene2function.de ), a database that includes comparative analyses for co-expression networks of 11 plant species. We exemplify how the tools included in PlaNet can be used to predict gene function, transfer knowledge, and discover conserved and multiplied gene modules.

Using gene networks in EvoDevo analyses

An intricate web of regulatory relationships between DNA, RNA, proteins and metabolites regulates how organisms achieve form and function. Genome sequencing combined with computational methods has allowed us to look at diverse readouts and generate a comprehensive framework for how molecules generate morphological phenotypes. RNAseq has evolved and proved useful for identifying links between transcription factor activity and transcript abundance, and for the generation of transcriptomes in non-model species through de novo assembly. Gene coexpression networks, combined with differential correlation analysis, offer the most versatile gene interaction exploratory tools, using gene expression patterns to determine potential associations and modularity of gene interactions and the differential associations between networks. Networks incorporating different types of biological data can help reduce the complexity of the entirety of biological information into discrete and interpretable pieces of data that can be the starting point for hypothesis construction and testing.

Co-expression network analysis reveals transcription factors associated to cell wall biosynthesis in sugarcane

Sugarcane is a hybrid of Saccharum officinarum and Saccharum spontaneum, with minor contributions from other species in Saccharum and other genera. Understanding the molecular basis of cell wall metabolism in sugarcane may allow for rational changes in fiber quality and content when designing new energy crops. This work describes a comparative expression profiling of sugarcane ancestral genotypes: S. officinarum, S. spontaneum and S. robustum and a commercial hybrid: RB867515, linking gene expression to phenotypes to identify genes for sugarcane improvement. Oligoarray experiments of leaves, immature and intermediate internodes, detected 12,621 sense and 995 antisense transcripts. Amino acid metabolism was particularly evident among pathways showing natural antisense transcripts expression. For all tissues sampled, expression analysis revealed 831, 674 and 648 differentially expressed genes in S. officinarum, S. robustum and S. spontaneum, respectively, using RB867515 as reference. Expression of sugar transporters might explain sucrose differences among genotypes, but an unexpected differential expression of histones were also identified between high and low Brix° genotypes. Lignin biosynthetic genes and bioenergetics-related genes were up-regulated in the high lignin genotype, suggesting that these genes are important for S. spontaneum to allocate carbon to lignin, while S. officinarum allocates it to sucrose storage. Co-expression network analysis identified 18 transcription factors possibly related to cell wall biosynthesis while in silico analysis detected cis-elements involved in cell wall biosynthesis in their promoters. Our results provide information to elucidate regulatory networks underlying traits of interest that will allow the improvement of sugarcane for biofuel and chemicals production.

Cellulose synthase gene expression profiling of Physcomitrella patens

The cellulose synthase (CESA) gene family of seed plants comprises six clades that encode isoforms with conserved expression patterns and distinct functions in cellulose synthesis complex (CSC) formation and primary and secondary cell wall synthesis. In mosses, which have rosette CSCs like those of seed plants but lack lignified secondary cell walls, the CESA gene family diversified independently and includes no members of the six functionally distinct seed plant clades. There are seven CESA isoforms encoded in the genome of the moss Physcomitrella patens. However, only PpCESA5 has been characterised functionally, and little information is available on the expression of other PpCESA family members. We have profiled PpCESA expression through quantitative RT-PCR, analysis of promoter-reporter lines, and cluster analysis of public microarray data in an effort to identify expression and co-expression patterns that could help reveal the functions of PpCESA isoforms in protein complex formation and development of specific tissues. In contrast to the tissue-specific expression observed for seed plant CESAs, each of the PpCESAs was broadly expressed throughout most developing tissues. Although a few statistically significant differences in expression of PpCESAs were noted when some tissues and hormone treatments were compared, no strong co-expression patterns were observed. Along with CESA phylogenies and lack of single PpCESA mutant phenotypes reported elsewhere, broad overlapping expression of the PpCESAs indicates a high degree of inter-changeability and is consistent with a different pattern of functional specialisation in the evolution of the seed plant and moss CESA families.

Tools of the trade: studying molecular networks in plants

Driven by recent technological improvements, genes can be now studied in a larger biological context. Genes and their protein products rarely operate as a single entity and large-scale mapping by protein-protein interactions can unveil the molecular complexes that form in the cell to carry out various functions. Expression analysis under multiple conditions, supplemented with protein-DNA binding data can highlight when genes are active and how they are regulated. Representing these data in networks and finding strongly connected sub-graphs has proven to be a powerful tool to predict the function of unknown genes. As such networks are gradually becoming available for various plant species, it becomes possible to study how networks evolve. This review summarizes currently available network data and related tools for plants. Furthermore we aim to provide an outlook of future analyses that can be done in plants based on work done in other fields.

Expression and functional analyses of a putative phenylcoumaran benzylic ether reductase in Arabidopsis thaliana

A candidate gene for phenylcoumaran benzylic ether reductase in Arabidopsis thaliana encodes a peptide with predicted functional activity and plays a crucial role in secondary metabolism. Phenylcoumaran benzylic ether reductase (PCBER) is thought to be an enzyme crucial in the biosynthesis of 8-5'-linked neolignans. Genes of the enzyme have been isolated and characterized in several plant species. In this study, we cloned cDNA and the 5'-untranslated region of one PCBER candidate gene (At4g39230, designated AtPCBER1) from Arabidopsis thaliana. At the amino acid level, AtPCBER1 shows high sequence identity (64-71 %) with PCBERs identified from other plant species. Expression analyses of AtPCBER1 by reverse transcriptase-polymerase chain reaction and histochemical analysis of transgenic plants harboring the 5'-untranslated region of AtPCBER1 linked with gus coding sequence indicate that expression is induced by wounding and is expressed in most tissues, including flower, stem, leaf, and root. Catalytic analysis of recombinant AtPCBER1 with neolignan and lignans in the presence of NADPH suggests that the protein can reduce not only the 8-5'-linked neolignan, dehydrodiconiferyl alcohol, but also 8-8' linked lignans, pinoresinol, and lariciresinol, with lower activities. To investigate further, we performed metabolomic analyses of transgenic plants in which the target gene was up- or down-regulated. Our results indicate no significant effects of AtPCBER1 gene regulation on plant growth and development; however, levels of some secondary metabolites, including lignans, flavonoids, and glucosinolates, differ between wild-type and transgenic plants. Taken together, our findings indicate that AtPCBER1 encodes a polypeptide with PCBER activity and has a critical role in the biosynthesis of secondary metabolites in A. thaliana.

FamNet: A Framework to Identify Multiplied Modules Driving Pathway Expansion in Plants

Gene duplications generate new genes that can acquire similar but often diversified functions. Recent studies of gene coexpression networks have indicated that, not only genes, but also pathways can be multiplied and diversified to perform related functions in different parts of an organism. Identification of such diversified pathways, or modules, is needed to expand our knowledge of biological processes in plants and to understand how biological functions evolve. However, systematic explorations of modules remain scarce, and no user-friendly platform to identify them exists. We have established a statistical framework to identify modules and show that approximately one-third of the genes of a plant's genome participate in hundreds of multiplied modules. Using this framework as a basis, we implemented a platform that can explore and visualize multiplied modules in coexpression networks of eight plant species. To validate the usefulness of the platform, we identified and functionally characterized pollen- and root-specific cell wall modules that multiplied to confer tip growth in pollen tubes and root hairs, respectively. Furthermore, we identified multiplied modules involved in secondary metabolite synthesis and corroborated them by metabolite profiling of tobacco (Nicotiana tabacum) tissues. The interactive platform, referred to as FamNet, is available at http://www.gene2function.de/famnet.html.

Transcriptome responses in alfalfa associated with tolerance to intensive animal grazing

Tolerance of alfalfa (Medicago sativa L.) to animal grazing varies widely within the species. However, the molecular mechanisms influencing the grazing tolerant phenotype remain uncharacterized. The objective of this study was to identify genes and pathways that control grazing response in alfalfa. We analyzed whole-plant de novo transcriptomes from grazing tolerant and intolerant populations of M. sativa ssp. falcata subjected to grazing by sheep. Among the Gene Ontology terms which were identified as grazing responsive in the tolerant plants and differentially enriched between the tolerant and intolerant populations (both grazed), most were associated with the ribosome and translation-related activities, cell wall processes, and response to oxygen levels. Twenty-one grazing responsive pathways were identified that also exhibited differential expression between the tolerant and intolerant populations. These pathways were associated with secondary metabolite production, primary carbohydrate metabolic pathways, shikimate derivative dependent pathways, ribosomal subunit composition, hormone signaling, wound response, cell wall formation, and anti-oxidant defense. Sequence polymorphisms were detected among several differentially expressed homologous transcripts between the tolerant and intolerant populations. These differentially responsive genes and pathways constitute potential response mechanisms for grazing tolerance in alfalfa. They also provide potential targets for molecular breeding efforts to develop grazing-tolerant cultivars of alfalfa.

A Meta-Analysis Based Method for Prioritizing Candidate Genes Involved in a Pre-specific Function

The identification of genes associated with a given biological function in plants remains a challenge, although network-based gene prioritization algorithms have been developed for Arabidopsis thaliana and many non-model plant species. Nevertheless, these network-based gene prioritization algorithms have encountered several problems; one in particular is that of unsatisfactory prediction accuracy due to limited network coverage, varying link quality, and/or uncertain network connectivity. Thus, a model that integrates complementary biological data may be expected to increase the prediction accuracy of gene prioritization. Toward this goal, we developed a novel gene prioritization method named RafSee, to rank candidate genes using a random forest algorithm that integrates sequence, evolutionary, and epigenetic features of plants. Subsequently, we proposed an integrative approach named RAP (Rank Aggregation-based data fusion for gene Prioritization), in which an order statistics-based meta-analysis was used to aggregate the rank of the network-based gene prioritization method and RafSee, for accurately prioritizing candidate genes involved in a pre-specific biological function. Finally, we showcased the utility of RAP by prioritizing 380 flowering-time genes in Arabidopsis. The "leave-one-out" cross-validation experiment showed that RafSee could work as a complement to a current state-of-art network-based gene prioritization system (AraNet v2). Moreover, RAP ranked 53.68% (204/380) flowering-time genes higher than AraNet v2, resulting in an 39.46% improvement in term of the first quartile rank. Further evaluations also showed that RAP was effective in prioritizing genes-related to different abiotic stresses. To enhance the usability of RAP for Arabidopsis and non-model plant species, an R package implementing the method is freely available at http://bioinfo.nwafu.edu.cn/software.

A R2R3-MYB transcription factor that is specifically expressed in cotton (Gossypium hirsutum) fibers affects secondary cell wall biosynthesis and deposition in transgenic Arabidopsis

Secondary cell wall (SCW) is an important industrial raw material for pulping, papermaking, construction, lumbering, textiles and potentially for biofuel production. The process of SCW thickening of cotton fibers lays down the cellulose that will constitute the bulk (up to 96%) of the fiber at maturity. In this study, a gene encoding a MYB-domain protein was identified in cotton (Gossypium hirsutum) and designated as GhMYBL1. Quantitative real-time polymerase chain reaction (RT-PCR) analysis revealed that GhMYBL1 was specifically expressed in cotton fibers at the stage of secondary wall deposition. Further analysis indicated that this protein is a R2R3-MYB transcription factor, and is targeted to the cell nucleus. Overexpression of GhMYBL1 in Arabidopsis affected the formation of SCW in the stem xylem of the transgenic plants. The enhanced SCW thickening also occurred in the interfascicular fibers, xylary fibers and vessels of the GhMYBL1-overexpression transgenic plants. The expression of secondary wall-associated genes, such as CesA4, CesA7, CesA8, PAL1, F5H and 4CL1, were upregulated, and consequently, cellulose and lignin biosynthesis were enhanced in the GhMYBL1 transgenic plants. These data suggested that GhMYBL1 may participate in modulating the process of secondary wall biosynthesis and deposition of cotton fibers.

Pinoresinol reductase 1 impacts lignin distribution during secondary cell wall biosynthesis in Arabidopsis

Pinoresinol reductase (PrR) catalyzes the conversion of the lignan (-)-pinoresinol to (-)-lariciresinol in Arabidopsis thaliana, where it is encoded by two genes, PrR1 and PrR2, that appear to act redundantly. PrR1 is highly expressed in lignified inflorescence stem tissue, whereas PrR2 expression is barely detectable in stems. Co-expression analysis has indicated that PrR1 is co-expressed with many characterized genes involved in secondary cell wall biosynthesis, whereas PrR2 expression clusters with a different set of genes. The promoter of the PrR1 gene is regulated by the secondary cell wall related transcription factors SND1 and MYB46. The loss-of-function mutant of PrR1 shows, in addition to elevated levels of pinoresinol, significantly decreased lignin content and a slightly altered lignin structure with lower abundance of cinnamyl alcohol end groups. Stimulated Raman scattering (SRS) microscopy analysis indicated that the lignin content of the prr1-1 loss-of-function mutant is similar to that of wild-type plants in xylem cells, which exhibit a normal phenotype, but is reduced in the fiber cells. Together, these data suggest an association of the lignan biosynthetic enzyme encoded by PrR1 with secondary cell wall biosynthesis in fiber cells.