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

Mapping of flumioxazin tolerance in a snap bean diversity panel leads to the discovery of a master genomic region controlling multiple stress resistance genes

Introduction

Effective weed management tools are crucial for maintaining the profitable production of snap bean (Phaseolus vulgaris L.). Preemergence herbicides help the crop to gain a size advantage over the weeds, but the few preemergence herbicides registered in snap bean have poor waterhemp (Amaranthus tuberculatus) control, a major pest in snap bean production. Waterhemp and other difficult-to-control weeds can be managed by flumioxazin, an herbicide that inhibits protoporphyrinogen oxidase (PPO). However, there is limited knowledge about crop tolerance to this herbicide. We aimed to quantify the degree of snap bean tolerance to flumioxazin and explore the underlying mechanisms.

Methods

We investigated the genetic basis of herbicide tolerance using genome-wide association mapping approach utilizing field-collected data from a snap bean diversity panel, combined with gene expression data of cultivars with contrasting response. The response to a preemergence application of flumioxazin was measured by assessing plant population density and shoot biomass variables.

Results

Snap bean tolerance to flumioxazin is associated with a single genomic location in chromosome 02. Tolerance is influenced by several factors, including those that are indirectly affected by seed size/weight and those that directly impact the herbicide's metabolism and protect the cell from reactive oxygen species-induced damage. Transcriptional profiling and co-expression network analysis identified biological pathways likely involved in flumioxazin tolerance, including oxidoreductase processes and programmed cell death. Transcriptional regulation of genes involved in those processes is possibly orchestrated by a transcription factor located in the region identified in the GWAS analysis. Several entries belonging to the Romano class, including Bush Romano 350, Roma II, and Romano Purpiat presented high levels of tolerance in this study. The alleles identified in the diversity panel that condition snap bean tolerance to flumioxazin shed light on a novel mechanism of herbicide tolerance and can be used in crop improvement.

TF-Prioritizer: a Java pipeline to prioritize condition-specific transcription factors

BACKGROUND

Eukaryotic gene expression is controlled by cis-regulatory elements (CREs), including promoters and enhancers, which are bound by transcription factors (TFs). Differential expression of TFs and their binding affinity at putative CREs determine tissue- and developmental-specific transcriptional activity. Consolidating genomic datasets can offer further insights into the accessibility of CREs, TF activity, and, thus, gene regulation. However, the integration and analysis of multimodal datasets are hampered by considerable technical challenges. While methods for highlighting differential TF activity from combined chromatin state data (e.g., chromatin immunoprecipitation [ChIP], ATAC, or DNase sequencing) and RNA sequencing data exist, they do not offer convenient usability, have limited support for large-scale data processing, and provide only minimal functionality for visually interpreting results.

RESULTS

We developed TF-Prioritizer, an automated pipeline that prioritizes condition-specific TFs from multimodal data and generates an interactive web report. We demonstrated its potential by identifying known TFs along with their target genes, as well as previously unreported TFs active in lactating mouse mammary glands. Additionally, we studied a variety of ENCODE datasets for cell lines K562 and MCF-7, including 12 histone modification ChIP sequencing as well as ATAC and DNase sequencing datasets, where we observe and discuss assay-specific differences.

CONCLUSION

TF-Prioritizer accepts ATAC, DNase, or ChIP sequencing and RNA sequencing data as input and identifies TFs with differential activity, thus offering an understanding of genome-wide gene regulation, potential pathogenesis, and therapeutic targets in biomedical research.

Transcriptome-wide analysis reveals core transcriptional regulators associated with culm development and variation in Dendrocalamus sinicus, the strongest woody bamboo in the world

Transcription factors (TFs) play indispensable roles in plant development and stress responses. As the largest woody bamboo species in the world, Dendrocalamus sinicus is endemic to Yunnan Province, China, and possesses two natural variants characterized by culm shape, namely straight or bent culms. Understanding the transcriptional regulation network of D. sinicus provides a unique opportunity to clarify the growth and development characteristics of woody bamboos. In this study, 10,236 TF transcripts belonging to 57 families were identified from transcriptome data of two variants at different developmental stages, from which we constructed a transcriptional regulatory network and unigene-coding protein-TFs interactive network of culm development for this attractive species. Gene function enrichment analysis revealed that hormone signaling and MAPK signaling pathways were two most enriched pathways in TF-regulated network. Based on PPI analysis, 50 genes interacting with nine TFs were screened as the core regulation components related to culm development. Of them, 18 synergistic genes of seven TFs, including nuclear cap-binding protein subunit 1, transcription factor GTE9-like, and ATP-dependent DNA helicase DDX11 isoform X1, involved in culm-shape variation. Most of these genes would interact with MYB, C3H, and ARF transcription factors. Six members with two each from ARF, C3H, and MYB transcription factor families and six key interacting genes (IAA3, IAA19, leucine-tRNA ligase, nuclear cap-binding protein subunit 1, elongation factor 2, and coiled-coil domain-containing protein 94) cooperate with these transcription factors were differentially expressed at development stage of young culms, and were validated by quantitative PCR. Our results represent a crucial step towards understanding the regulatory mechanisms of TFs involved in culm development and variation of D. sinicus.

Transcriptome Analysis of Developing Grains from Wheat Cultivars TAM 111 and TAM 112 Reveal Cultivar-Specific Regulatory Networks

Wheat flour's end-use quality is tightly linked to the quantity and composition of storage proteins in the endosperm. TAM 111 and TAM 112 are two popular cultivars grown in the Southern US Great Plains with significantly different protein content. To investigate regulatory differences, transcriptome data were analyzed from developing grains at early- and mid-filling stages. At the mid-filling stage, TAM 111 preferentially upregulated starch metabolism-related pathways compared to TAM 112, whereas amino acid metabolism and transporter-related pathways were over-represented in TAM 112. Elemental analyses also indicated a higher N percentage in TAM 112 at the mid-filling stage. To explore the regulatory variation, weighted correlation gene network was constructed from publicly available RNAseq datasets to identify the modules differentially regulated in TAM 111 and TAM 112. Further, the potential transcription factors (TFs) regulating those modules were identified using graphical least absolute shrinkage and selection operator (GLASSO). Homologs of the OsNF-Y family members with known starch metabolism-related functions showed higher connectivities in TAM 111. Multiple TFs with high connectivity in TAM 112 had predicted functions associated with ABA response in grain. These results will provide novel targets for breeders to explore and further our understanding in mechanisms regulating grain development.

Antisense Transcription in Plants: A Systematic Review and an Update on cis-NATs of Sugarcane

Initially, natural antisense transcripts (NATs, natRNAs, or asRNAs) were considered repressors; however, their functions in gene regulation are diverse. Positive, negative, or neutral correlations to the cognate gene expression have been noted. Although the first studies were published about 50 years ago, there is still much to be investigated regarding antisense transcripts in plants. A systematic review of scientific publications available in the Web of Science databases was conducted to contextualize how the studying of antisense transcripts has been addressed. Studies were classified considering three categories: “Natural antisense” (208), artificial antisense used in “Genetic Engineering” (797), or “Natural antisense and Genetic Engineering”-related publications (96). A similar string was used for a systematic search in the NCBI Gene database. Of the 1132 antisense sequences found for plants, only 0.8% were cited in PubMed and had antisense information confirmed. This value was the lowest when compared to fungi (2.9%), bacteria (2.3%), and mice (54.1%). Finally, we present an update for the cis-NATs identified in Saccharum spp. Of the 1413 antisense transcripts found in different experiments, 25 showed concordant expressions, 22 were discordant, 1264 did not correlate with the cognate genes, and 102 presented variable results depending on the experiment.

NetREx: Network-based Rice Expression Analysis Server for abiotic stress conditions

Recent focus on transcriptomic studies in food crops like rice, wheat and maize provide new opportunities to address issues related to agriculture and climate change. Re-analysis of such data available in public domain supplemented with annotations across molecular hierarchy can be of immense help to the plant research community, particularly co-expression networks representing transcriptionally coordinated genes that are often part of the same biological process. With this objective, we have developed NetREx, a Network-based Rice Expression Analysis Server, that hosts ranked co-expression networks of Oryza sativa using publicly available messenger RNA sequencing data across uniform experimental conditions. It provides a range of interactable data viewers and modules for analysing user-queried genes across different stress conditions (drought, flood, cold and osmosis) and hormonal treatments (abscisic and jasmonic acid) and tissues (root and shoot). Subnetworks of user-defined genes can be queried in pre-constructed tissue-specific networks, allowing users to view the fold change, module memberships, gene annotations and analysis of their neighbourhood genes and associated pathways. The web server also allows querying of orthologous genes from Arabidopsis, wheat, maize, barley and sorghum. Here, we demonstrate that NetREx can be used to identify novel candidate genes and tissue-specific interactions under stress conditions and can aid in the analysis and understanding of complex phenotypes linked to stress response in rice.

Database URL: https://bioinf.iiit.ac.in/netrex/index.html

Transcriptome changes in the developing sugarcane culm associated with high yield and early-season high sugar content

Sugarcane, with its exceptional carbon dioxide assimilation, biomass and sugar yield, has a high potential for the production of bio-energy, bio-plastics and high-value products in the food and pharmaceutical industries. A crucial challenge for long-term economic viability and environmental sustainability is also to optimize the production of biomass composition and carbon sequestration. Sugarcane varieties such as KQ228 and Q253 are highly utilized in the industry. These varieties are characterized by a high early-season sugar content associated with high yield. In order to investigate these correlations, 1,440 internodes were collected and combined to generate a set of 120 samples in triplicate across 24 sugarcane cultivars at five different development stages. Weighted gene co-expression network analysis (WGCNA) was used and revealed for the first time two sets of co-expressed genes with a distinct and opposite correlation between fibre and sugar content. Gene identification and metabolism pathways analysis was used to define these two sets of genes. Correlation analysis identified a large number of interconnected metabolic pathways linked to sugar content and fibre content. Unsupervised hierarchical clustering of gene expression revealed a stronger level of segregation associated with the genotypes than the stage of development, suggesting a dominant genetic influence on biomass composition and facilitating breeding selection. Characterization of these two groups of co-expressed key genes can help to improve breeding program for high fibre, high sugar species or plant synthetic biology.

Isolation of Promoters and Transcription Factors Involved in the Regulation of Lignin Biosynthesis in Saccharum Species

Sugarcane bagasse has received attention as a raw material for the production of second-generation ethanol (E2G). However, its use is limited because of the cell wall recalcitrance, mostly conferred by lignin. Recently our knowledge of the genes coding for the enzymes of the lignin biosynthesis pathway has increased; however, still little is known about the transcription factors controlling the expression of these genes in sugarcane. Here we describe protocols to optimize the isolation of the promoters of the lignin biosynthetic genes ShCAD8, ShCOMT and ShF5H and the transcription factors (TFs) ShMYB85 and ShMYB58/63 in Saccharum species. To confirm whether these TFs are able to activate the target promoters, a transactivation assay in BY2 protoplasts of Nicotiana tabacum is also detailed.

Visualization of Suberization and Lignification in Sugarcane

Cell wall biopolymers are major factors responsible for the high recalcitrance of sugarcane biomass. The study of suberization and lignification mechanisms in sugarcane and of the networks that control biosynthesis of these polymers will contribute to the biotechnological improvement of this crop. Here, we describe experiments that allow the visualization of the suberization and lignification mechanism in response to mechanical injury in sugarcane.

Gene Co-expression Network Analysis

Gene co-expression analysis is a data analysis technique that helps identify groups of genes with similar expression patterns across several different conditions. By means of these techniques, different groups have been able to assign putative metabolic pathways and functions to understudied genes and to identify novel metabolic regulation networks for different metabolites. Some groups have even used network comparative studies to understand the evolution of these networks from green algae to land plants. In this chapter, we will go over the basic definitions required to understand network topology and gene module identification. Additionally, we offer the reader a walk-through a standard analysis pipeline as implemented in the package WGCNA that takes as input raw fastq files and obtains co-expressed gene clusters and representative gene expression patterns from each module for downstream applications.

Planting Season Impacts Sugarcane Stem Development, Secondary Metabolite Levels, and Natural Antisense Transcription

To reduce the potentially irreversible environmental impacts caused by fossil fuels, the use of renewable energy sources must be increased on a global scale. One promising source of biomass and bioenergy is sugarcane. The study of this crop's development in different planting seasons can aid in successfully cultivating it in global climate change scenarios. The sugarcane variety SP80-3280 was field grown under two planting seasons with different climatic conditions. A systems biology approach was taken to study the changes on physiological, morphological, agrotechnological, transcriptomics, and metabolomics levels in the leaf +1, and immature, intermediate and mature internodes. Most of the variation found within the transcriptomics and metabolomics profiles is attributed to the differences among the distinct tissues. However, the integration of both transcriptomics and metabolomics data highlighted three main metabolic categories as the principal sources of variation across tissues: amino acid metabolism, biosynthesis of secondary metabolites, and xenobiotics biodegradation and metabolism. Differences in ripening and metabolite levels mainly in leaves and mature internodes may reflect the impact of contrasting environmental conditions on sugarcane development. In general, the same metabolites are found in mature internodes from both "one-year" and "one-and-a-half-year sugarcane", however, some metabolites (i.e., phenylpropanoids with economic value) and natural antisense transcript expression are only detected in the leaves of "one-year" sugarcane.

Cell wall formation pathways are differentially regulated in sugarcane contrasting genotypes associated with endophytic diazotrophic bacteria

Differences in cell wall components between two BNF-contrasting sugarcane genotypes might result from genetic variations particular to the genotype and from the efficiency in diazotrophic bacteria association. Sugarcane is a plant of the grass family (Poaceae) that is highly cultivated in Brazil, as an important energy resource. Commercial sugarcane genotypes may be successfully associated with beneficial endophytic nitrogen-fixing bacteria, which can influence several plant metabolic pathways, such as cell division and growth, synthesis of hormones, and defense compounds. In this study, we investigated how diazotrophic bacteria associated with sugarcane plants could be involved in the regulation of cell wall formation pathways. A molecular and structural characterization of the cell wall was compared between two genotypes of sugarcane with contrasting rates of Biological Nitrogen Fixation (BNF): SP70-1143 (high BNF) and Chunee (low BNF). Differentially expressed transcripts were identified in transcriptomes generated from SP70-1143 and Chunee. Expression profiles of cellulose and lignin genes, which were more expressed in SP70-1134, and callose genes, which were more expressed in Chunee, were validated by RT-qPCR and microscopic analysis of cell wall components in tissue sections. A similar expression profile in both BNF-contrasting genotypes was observed in naturally colonized plants and in plants inoculated with G. diazotrophicus. Cell walls of the high BNF genotype have a greater cellulose content, which might contribute to increase biomass. In parallel, callose was concentrated in the vascular tissues of the low BNF genotype and could possibly represent a barrier for an efficient bacterial colonization and dissemination in sugarcane tissues. Our data show a correlation between the gene profiles identified in the BNF-contrasting genotypes and a successful association with endophytic diazotrophic bacteria.

Transcriptional reprogramming of xylem cell wall biosynthesis in tension wood

Tension wood (TW) is a specialized xylem tissue developed under mechanical/tension stress in angiosperm trees. TW development involves transregulation of secondary cell wall genes, which leads to altered wood properties for stress adaptation. We induced TW in the stems of black cottonwood (Populus trichocarpa, Nisqually-1) and identified two significantly repressed transcription factor (TF) genes: class B3 heat-shock TF (HSFB3-1) and MYB092. Transcriptomic analysis and chromatin immunoprecipitation (ChIP) were used to identify direct TF-DNA interactions in P. trichocarpa xylem protoplasts overexpressing the TFs. This analysis established a transcriptional regulatory network in which PtrHSFB3-1 and PtrMYB092 directly activate 8 and 11 monolignol genes, respectively. The TF-DNA interactions were verified for their specificity and transactivator roles in 35 independent CRISPR-based biallelic mutants and overexpression transgenic lines of PtrHSFB3-1 and PtrMYB092 in P. trichocarpa. The gene-edited trees (mimicking the repressed PtrHSFB3-1 and PtrMYB092 under tension stress) have stem wood composition resembling that of TW during normal growth and under tension stress (i.e., low lignin and high cellulose), whereas the overexpressors showed an opposite effect (high lignin and low cellulose). Individual overexpression of the TFs impeded lignin reduction under tension stress and restored high levels of lignin biosynthesis in the TW. This study offers biological insights to further uncover how metabolism, growth, and stress adaptation are coordinately regulated in trees.

Combined Analysis of Volatile Terpenoid Metabolism and Transcriptome Reveals Transcription Factors Related to Terpene Synthase in Two Cultivars of Dendrobium officinale Flowers

Dendrobium officinale is a kind of traditional Chinese herbal medicine. Its flowers could be used as health care tea for its aroma flavor and medicinal value. Most recent studies demonstrated that terpenoids are the main components of the aromatic compounds in the flowers, but the biosynthesis of terpenoids is poorly understood in D. officinale. In the experiment, the flowers from two cultivars of D. officinale with different smells were collected. The transcriptome analysis and combined volatile terpenoids determination were performed to identify the genes related to the biosynthesis of the terpenoids. The results showed that the different products of volatile terpenoids are α-thujene, linalool, α-terpineol, α-phellandrene, γ-muurolene, α-patchoulene, and δ-elemene in two cultivar flowers. The transcriptome analysis detected 25,484 genes in the flowers. And 18,650 differentially expressed genes were identified between the two cultivars. Of these genes, 253 genes were mapped to the terpenoid metabolism pathway. Among these genes, 13 terpene synthase (TPS) genes may have correlations with AP2/ERF, WRKY, MYB, bHLH, and bZIP transcription factors by weighted gene co-expression network analysis (WGCNA). The transcription factors have regulatory effects on TPS genes. These results may provide ideas for the terpenoid biosynthesis and regulatory network of D. officinale flowers.

Applying Molecular Phenotyping Tools to Explore Sugarcane Carbon Potential

Sugarcane (Saccharum spp.), a C₄ grass, has a peculiar feature: it accumulates, gradient-wise, large amounts of carbon (C) as sucrose in its culms through a complex pathway. Apart from being a sustainable crop concerning C efficiency and bioenergetic yield per hectare, sugarcane is used as feedstock for producing ethanol, sugar, high-value compounds, and products (e.g., polymers and succinate), and bioelectricity, earning the title of the world's leading biomass crop. Commercial cultivars, hybrids bearing high levels of polyploidy, and aneuploidy, are selected from a large number of crosses among suitable parental genotypes followed by the cloning of superior individuals among the progeny. Traditionally, these classical breeding strategies have been favoring the selection of cultivars with high sucrose content and resistance to environmental stresses. A current paradigm change in sugarcane breeding programs aims to alter the balance of C partitioning as a means to provide more plasticity in the sustainable use of this biomass for metabolic engineering and green chemistry. The recently available sugarcane genetic assemblies powered by data science provide exciting perspectives to increase biomass, as the current sugarcane yield is roughly 20% of its predicted potential. Nowadays, several molecular phenotyping tools can be applied to meet the predicted sugarcane C potential, mainly targeting two competing pathways: sucrose production/storage and biomass accumulation. Here we discuss how molecular phenotyping can be a powerful tool to assist breeding programs and which strategies could be adopted depending on the desired final products. We also tackle the advances in genetic markers and mapping as well as how functional genomics and genetic transformation might be able to improve yield and saccharification rates. Finally, we review how "omics" advances are promising to speed up plant breeding and reach the unexplored potential of sugarcane in terms of sucrose and biomass production.

Differentiation of Tracheary Elements in Sugarcane Suspension Cells Involves Changes in Secondary Wall Deposition and Extensive Transcriptional Reprogramming

Plant lignocellulosic biomass, mostly composed of polysaccharide-rich secondary cell walls (SCWs), provides fermentable sugars that may be used to produce biofuels and biomaterials. However, the complex chemical composition and physical structure of SCWs hinder efficient processing of plant biomass. Understanding the molecular mechanisms underlying SCW deposition is, thus, essential to optimize bioenergy feedstocks. Here, we establish a xylogenic culture as a model system to study SCW deposition in sugarcane; the first of its kind in a C4 grass species. We used auxin and brassinolide to differentiate sugarcane suspension cells into tracheary elements, which showed metaxylem-like reticulate or pitted SCW patterning. The differentiation led to increased lignin levels, mainly caused by S-lignin units, and a rise in p-coumarate, leading to increased p-coumarate:ferulate ratios. RNAseq analysis revealed massive transcriptional reprogramming during differentiation, with upregulation of genes associated with cell wall biogenesis and phenylpropanoid metabolism and downregulation of genes related to cell division and primary metabolism. To better understand the differentiation process, we constructed regulatory networks of transcription factors and SCW-related genes based on co-expression analyses. Accordingly, we found multiple regulatory modules that may underpin SCW deposition in sugarcane. Our results provide important insights and resources to identify biotechnological strategies for sugarcane biomass optimization.

Genetic linkage map construction and quantitative trait loci mapping of agronomic traits in Gloeostereum incarnatum

Gloeostereum incarnatum is an edible medicinal mushroom widely grown in China. Using the whole genome of G. incarnatum, simple sequence repeat (SSR) markers were developed and synthetic primers were designed to construct its first genetic linkage map. The 1,048.6 cm map is composed of 10 linkage groups and contains 183 SSR markers. In total, 112 genome assembly sequences were anchored, representing 16.43 Mb and covering 46.41% of the genome. Selfing populations were used for quantitative trait loci (QTL) targeting, and the composite interval mapping method was used to co-localize the mycelium growth rate (potato dextrose agar and sawdust), growth period, yield and fruiting body length, and width and thickness. The 14 QTLs of agronomic traits had LOD values of 3.20-6.51 and contribution rates of 2.22-13.18%. No linkage relationship was found between the mycelium growth rate and the growth period, but a linkage relationship was observed among the length, width and thickness of the fruiting bodies. Using NCBI's BLAST alignment, the genomic sequences corresponding to the QTL regions were compared, and a TPR-like protein candidate gene was selected. Using whole-genome data, 138 candidate genes were found in four sequence fragments of two SSR markers located in the same scaffold. The genetic map and QTLs established in this study will aid in developing selective markers for agronomic traits and identifying corresponding genes, thereby providing a scientific basis for the further gene mapping of quantitative traits and the marker-assisted selection of functional genes in G. incarnatum breeding programs.

The sugarcane ShMYB78 transcription factor activates suberin biosynthesis in Nicotiana benthamiana

KEY MESSAGE

A sugarcane MYB present in the culm induces suberin biosynthesis and is involved both with fatty acid and phenolics metabolism. Few transcription factors have been described as regulators of cell wall polymers deposition in C4 grasses. Particularly, regulation of suberin biosynthesis in this group of plants remains poorly understood. Here, we showed that the sugarcane MYB transcription factor ShMYB78 is an activator of suberin biosynthesis and deposition. ShMYB78 was identified upon screening genes whose expression was upregulated in sugarcane internodes undergoing suberization during culm development or triggered by wounding. Agrobacterium-mediated transient expression of ShMYB78 in Nicotiana benthamiana leaves induced the ectopic deposition of suberin and its aliphatic and aromatic monomers. Further, the expression of suberin-related genes was induced by ShMYB78 heterologous expression in Nicotiana benthamiana leaves. ShMYB78 was shown to be a nuclear protein based on its presence in sugarcane internode nuclear protein extracts, and protoplast transactivation assays demonstrated that ShMYB78 activates the promoters of the sugarcane suberin biosynthetic genes β-ketoacyl-CoA synthase (ShKCS20) and caffeic acid-O-methyltransferase (ShCOMT). Our results suggest that ShMYB78 may be involved in the transcriptional regulation of suberin deposition, from fatty acid metabolism to phenylpropanoid biosynthesis, in sugarcane internodes.

Differential expression in leaves of Saccharum genotypes contrasting in biomass production provides evidence of genes involved in carbon partitioning

BACKGROUND

The development of biomass crops aims to meet industrial yield demands, in order to optimize profitability and sustainability. Achieving these goals in an energy crop like sugarcane relies on breeding for sucrose accumulation, fiber content and stalk number. To expand the understanding of the biological pathways related to these traits, we evaluated gene expression of two groups of genotypes contrasting in biomass composition.

RESULTS

First visible dewlap leaves were collected from 12 genotypes, six per group, to perform RNA-Seq. We found a high number of differentially expressed genes, showing how hybridization in a complex polyploid system caused extensive modifications in genome functioning. We found evidence that differences in transposition and defense related genes may arise due to the complex nature of the polyploid Saccharum genomes. Genotypes within both biomass groups showed substantial variability in genes involved in photosynthesis. However, most genes coding for photosystem components or those coding for phosphoenolpyruvate carboxylases (PEPCs) were upregulated in the high biomass group. Sucrose synthase (SuSy) coding genes were upregulated in the low biomass group, showing that this enzyme class can be involved with sucrose synthesis in leaves, similarly to sucrose phosphate synthase (SPS) and sucrose phosphate phosphatase (SPP). Genes in pathways related to biosynthesis of cell wall components and expansins coding genes showed low average expression levels and were mostly upregulated in the high biomass group.

CONCLUSIONS

Together, these results show differences in carbohydrate synthesis and carbon partitioning in the source tissue of distinct phenotypic groups. Our data from sugarcane leaves revealed how hybridization in a complex polyploid system resulted in noticeably different transcriptomic profiles between contrasting genotypes.

Identification of genes from the general phenylpropanoid and monolignol-specific metabolism in two sugarcane lignin-contrasting genotypes

The phenylpropanoid pathway is an important route of secondary metabolism involved in the synthesis of different phenolic compounds such as phenylpropenes, anthocyanins, stilbenoids, flavonoids, and monolignols. The flux toward monolignol biosynthesis through the phenylpropanoid pathway is controlled by specific genes from at least ten families. Lignin polymer is one of the major components of the plant cell wall and is mainly responsible for recalcitrance to saccharification in ethanol production from lignocellulosic biomass. Here, we identified and characterized sugarcane candidate genes from the general phenylpropanoid and monolignol-specific metabolism through a search of the sugarcane EST databases, phylogenetic analysis, a search for conserved amino acid residues important for enzymatic function, and analysis of expression patterns during culm development in two lignin-contrasting genotypes. Of these genes, 15 were cloned and, when available, their loci were identified using the recently released sugarcane genomes from Saccharum hybrid R570 and Saccharum spontaneum cultivars. Our analysis points out that ShPAL1, ShPAL2, ShC4H4, Sh4CL1, ShHCT1, ShC3H1, ShC3H2, ShCCoAOMT1, ShCOMT1, ShF5H1, ShCCR1, ShCAD2, and ShCAD7 are strong candidates to be bona fide lignin biosynthesis genes. Together, the results provide information about the candidate genes involved in monolignol biosynthesis in sugarcane and may provide useful information for further molecular genetic studies in sugarcane.

Assembly of the 373k gene space of the polyploid sugarcane genome reveals reservoirs of functional diversity in the world's leading biomass crop

BACKGROUND

Sugarcane cultivars are polyploid interspecific hybrids of giant genomes, typically with 10-13 sets of chromosomes from 2 Saccharum species. The ploidy, hybridity, and size of the genome, estimated to have >10 Gb, pose a challenge for sequencing.

RESULTS

Here we present a gene space assembly of SP80-3280, including 373,869 putative genes and their potential regulatory regions. The alignment of single-copy genes in diploid grasses to the putative genes indicates that we could resolve 2-6 (up to 15) putative homo(eo)logs that are 99.1% identical within their coding sequences. Dissimilarities increase in their regulatory regions, and gene promoter analysis shows differences in regulatory elements within gene families that are expressed in a species-specific manner. We exemplify these differences for sucrose synthase (SuSy) and phenylalanine ammonia-lyase (PAL), 2 gene families central to carbon partitioning. SP80-3280 has particular regulatory elements involved in sucrose synthesis not found in the ancestor Saccharum spontaneum. PAL regulatory elements are found in co-expressed genes related to fiber synthesis within gene networks defined during plant growth and maturation. Comparison with sorghum reveals predominantly bi-allelic variations in sugarcane, consistent with the formation of 2 "subgenomes" after their divergence ∼3.8-4.6 million years ago and reveals single-nucleotide variants that may underlie their differences.

CONCLUSIONS

This assembly represents a large step towards a whole-genome assembly of a commercial sugarcane cultivar. It includes a rich diversity of genes and homo(eo)logous resolution for a representative fraction of the gene space, relevant to improve biomass and food production.

Genomic resources for energy cane breeding in the post genomics era

Sugarcane is one of the most sustainable energy crops among cultivated crops presenting the highest tonnage of cultivated plants. Its high productivity of sugar, bioethanol and bioelectricity make it a promising green alternative to petroleum. Furthermore, the myriad of products that can be derived from sugarcane biomass has been driving breeding programs towards varieties with a higher yield of fiber and a more vigorous and sustainable performance: the energy cane. Here we provide an overview of the energy cane including plant description, breeding efforts, types, and end-uses. In addition, we describe recently published genomic resources for the development of this crop, discuss current knowledge of cell wall metabolism, bioinformatic tools and databases available for the community.

The lignin toolbox of the model grass Setaria viridis

The core set of biosynthetic genes potentially involved in developmental lignification was identified in the model C4 grass Setaria viridis. Lignin has been recognized as a major recalcitrant factor negatively affecting the processing of plant biomass into bioproducts. However, the efficient manipulation of lignin deposition in order to generate optimized crops for the biorefinery requires a fundamental knowledge of several aspects of lignin metabolism, including regulation, biosynthesis and polymerization. The current availability of an annotated genome for the model grass Setaria viridis allows the genome-wide characterization of genes involved in the metabolic pathway leading to the production of monolignols, the main building blocks of lignin. Here we performed a comprehensive study of monolignol biosynthetic genes as an initial step into the characterization of lignin metabolism in S. viridis. A total of 56 genes encoding bona fide enzymes catalyzing the consecutive ten steps of the monolignol biosynthetic pathway were identified in the S. viridis genome. A combination of comparative phylogenetic studies, high-throughput expression analysis and quantitative RT-PCR analysis was further employed to identify the family members potentially involved in developmental lignification. Accordingly, 14 genes clustered with genes from closely related species with a known function in lignification and showed an expression pattern that correlates with lignin deposition. These genes were considered the "core lignin toolbox" responsible for the constitutive, developmental lignification in S. viridis. These results provide the basis for further understanding lignin deposition in C4 grasses and will ultimately allow the validation of biotechnological strategies to produce crops with enhanced processing properties.

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.

Integration of Metabolite Profiling and Transcriptome Analysis Reveals Genes Related to Volatile Terpenoid Metabolism in Finger Citron (C. medica var. sarcodactylis)

Finger citron (Citrus medica var. sarcodactylis) is a popular ornamental tree and an important source of essential oils rich in terpenoids, but the mechanisms behind volatile formation are poorly understood. We investigated gene expression changes combined with volatile profiling of ten samples from three developing organs: flower, leaf, and fruit. A total of 62 volatiles were identified with limonene and γ-terpinene being the most abundant ones. Six volatiles were identified using partial least squares discriminant analysis (PLS-DA) that could be used as markers for distinguishing finger citron from other citrus species. RNA-Seq revealed 1,611,966,118 high quality clean reads that were assembled into 32,579 unigenes. From these a total of 58 terpene synthase (TPS) gene family members were identified and the spatial and temporal distribution of their transcripts was measured in developing organs. Transcript levels of transcription factor genes AP2/ERF (251), bHLH (169), bZIP (76), MYB (155), NAC (184), and WRKY (66) during finger citron development were also analyzed. From extracted subnetworks of three modules constructed by weighted gene co-expression network analysis (WGCNA), thirteen TPS genes and fifteen transcription factors were suggested to be related to volatile terpenoid formation. These results provide a framework for future investigations into the identification and regulatory network of terpenoids in finger citron.

QTL mapping for maize starch content and candidate gene prediction combined with co-expression network analysis

A major QTL Qsta9.1 was identified on chromosome 9, combined with GWAS, and co-expression network analysis showed that GRMZM2G110929 and GRMZM5G852704 are the potential candidates for association with maize kernel starch content. Increasing maize kernel starch content may not only lead to higher maize kernel yields and qualities, but also help meet industry demands. By using the intermated B73 × Mo17 population, QTLs were mapped for starch content in this study. A major QTL Qsta9.1 was detected in a 1.7 Mb interval on chromosome 9 and validated by allele frequency analysis in extreme tails of a newly constructed segregating population. According to genome-wide association study (GWAS) based on genotyping of a natural population, we identified a significant SNP for starch content within the ORF region of GRMZM5G852704_T01 colocalized with QTL Qsta9.1. Co-expression network analysis was also conducted, and 28 modules were constructed during six seed developmental stages. Functional enrichment was performed for each module, and one module showed the most possibility for the association with carbohydrate-related processes. In this module, one transcripts GRMZM2G110929_T01 located in the Qsta9.1 assigned 1.7 Mb interval encoding GLABRA2 expression modulator. Its expression level in B73 was lower than that in Mo17 across all seed developmental stages, implying the possibility for the candidate gene of Qsta9.1. Our studies combined GWAS, mRNA profiling, and traditional QTL analyses to identify a major locus for controlling seed starch content in maize.

Diel oscillations in cell wall components and soluble sugars as a response to short-day in sugarcane (Saccharum sp.)

BACKGROUND

Sugarcane is a tropical crop that can accumulate high concentration of sucrose in the stem as a storage carbohydrate. For that reason, sugarcane accounts for approximately 75% of all the sugar produced in the world and has become the main sugar source to produce first-generation bioethanol in Brazil. Daily rhythms cause plants to adapt and coordinate their metabolism to achieve maximum photosynthesis and carbohydrate production throughout the day. Circadian rhythms arise from the interaction of an internal oscillator and external stimuli, whereas diel rhythms occur in response to a light-dark cycle. Diel signalling contributes to synchronizing circadian rhythms to photoperiods, and levels of carbohydrates oscillate in a diel fashion. Under regular photoperiods, they are synthesized during the daytime and consumed throughout the night as an energy reserve. However, short days can induce higher rates of synthesis during daytime and lower rates of consumption in the dark. Cell wall carbohydrates are also diurnally regulated, and it has been shown that celluloses, hemicelluloses and pectin are deposited/degraded at different times of the day. To assess the diel carbohydrate profile in young sugarcane plants, we measured soluble sugars and cell wall components along a time course in plants subjected either to a regular day or short day.

RESULTS

Short-day influenced sucrose synthesis and cell wall components. In short-day a 44% increase in sucrose concentration was detected in the dark, but was stable during the day. Cellulose, hemicellulose and pectin also fluctuate within a 24 h interval when subjected to a short day. A 38% increase in leaf sheath cellulose was observed from the middle of the day to the first hour of the night. Leaf sheath pectin and hemicellulose also increased from the day to the night, while it decreased in leaves.

CONCLUSIONS

The presented data show diurnal patterns of soluble sugar metabolism together with temporal regulation of cell wall metabolism for a short day, suggesting that diel signalling has a role in how sugarcane manages sugar accumulation and partitioning. Understanding cell wall synthesis/degradation dynamics may help to improve the yield of sugarcane.

Deposition of lignin in four species of Saccharum

We used primers designed on conserved gene regions of several species to isolate the most expressed genes of the lignin pathway in four Saccharum species. S. officinarum and S. barberi have more sucrose in the culms than S. spontaneum and S. robustum, but less polysaccharides and lignin in the cell wall. S. spontaneum, and S. robustum had the lowest S/G ratio and a lower rate of saccharification in mature internodes. Surprisingly, except for CAD, 4CL, and CCoAOMT for which we found three, two, and two genes, respectively, only one gene was found for the other enzymes and their sequences were highly similar among the species. S. spontaneum had the highest expression for most genes. CCR and CCoAOMT B presented the highest expression; 4CL and F5H showed increased expression in mature tissues; C3H and CCR had higher expression in S. spontaneum, and one of the CADs isolated (CAD B) had higher expression in S. officinarum. The similarity among the most expressed genes isolated from these species was unexpected and indicated that lignin biosynthesis is conserved in Saccharum including commercial varieties Thus the lignin biosynthesis control in sugarcane may be only fully understood with the knowledge of the promotor region of each gene.

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.

The control of endopolygalacturonase expression by the sugarcane RAV transcription factor during aerenchyma formation

The development of lysigenous aerenchyma starts with cell expansion and degradation of pectin from the middle lamella, leading to cell wall modification, and culminating with cell separation. Here we report that nutritional starvation of sugarcane induced gene expression along sections of the first 5 cm of the root and between treatments. We selected two candidate genes: a RAV transcription factor, from the ethylene response factors superfamily, and an endopolygalacturonase (EPG), a glycosyl hydrolase related to homogalacturonan hydrolysis from the middle lamella. epg1 and rav1 transcriptional patterns suggest they are essential genes at the initial steps of pectin degradation during aerenchyma development in sugarcane. Due to the high complexity of the sugarcane genome, rav1 and epg1 were sequenced from 17 bacterial artificial chromosome clones containing hom(e)ologous genomic regions, and the sequences were compared with those of Sorghum bicolor. We used one hom(e)olog sequence from each gene for transactivation assays in tobacco. rav1 was shown to bind to the epg1 promoter, repressing β-glucuronidase activity. RAV repression upon epg1 transcription is the first reported link between ethylene regulation and pectin hydrolysis during aerenchyma formation. Our findings may help to elucidate cell wall degradation in sugarcane and therefore contribute to second-generation bioethanol production.

A Trihelix Family Transcription Factor Is Associated with Key Genes in Mixed-Linkage Glucan Accumulation

Mixed-linkage glucan (MLG) is a polysaccharide that is highly abundant in grass endosperm cell walls and present at lower amounts in other tissues. Cellulose synthase-like F (CSLF) and cellulose synthase-like H genes synthesize MLG, but it is unknown if other genes participate in the production and restructuring of MLG. Using Brachypodium distachyon transcriptional profiling data, we identified a B distachyon trihelix family transcription factor (BdTHX1) that is highly coexpressed with the B distachyon CSLF6 gene (BdCSLF6), which suggests that BdTHX1 is involved in the regulation of MLG biosynthesis. To determine the genes regulated by this transcription factor, we conducted chromatin immunoprecipitation sequencing (ChIP-seq) experiments using immature B distachyon seeds and an anti-BdTHX1 polyclonal antibody. The ChIP-seq experiment identified the second intron of BdCSLF6 as one of the most enriched sequences. The binding of BdTHX1 to the BdCSLF6 intron sequence was confirmed using electrophoretic mobility shift assays (EMSA). ChIP-seq also showed that a gene encoding a grass-specific glycoside hydrolase family 16 endotransglucosylase/hydrolase (BdXTH8) is bound by BdTHX1, and the binding was confirmed by EMSA. Radiochemical transglucanase assays showed that BdXTH8 exhibits predominantly MLG:xyloglucan endotransglucosylase activity, a hetero-transglycosylation reaction, and can thus produce MLG-xyloglucan covalent bonds; it also has a lower xyloglucan:xyloglucan endotransglucosylase activity. B distachyon shoots regenerated from transformed calli overexpressing BdTHX1 showed an abnormal arrangement of vascular tissue and seedling-lethal phenotypes. These results indicate that the transcription factor BdTHX1 likely plays an important role in MLG biosynthesis and restructuring by regulating the expression of BdCSLF6 and BdXTH8.

Optimization of RT-PCR reactions in studies with genes of lignin biosynthetic route in Saccharum spontaneum

Saccharum spontaneum has been used for the development of energy cane a crop aimed to be used for the production of second-generation ethanol, or lignocellulosic ethanol. Lignin is a main challenge in the conversion of cell wall sugars into ethanol. In our studies to isolate the genes the lignin biosynthesis in S. spontaneum we have had great difficulty in RT-PCR reactions. Thus, we evaluated the effectiveness of different additives in the amplification of these genes. While COMT and CCoAOMT genes did not need any additives for other genes there was no amplification (HCT, F5H, 4CL and CCR) or the yield was very low (CAD and C4H). The application of supplementary cDNA was enough to overcome the non-specificity and low yield for C4H and C3H, while the addition of 0.04% BSA + 2% formamide was effective to amplify 4CL, CCR, F5H and CCR. HCT was amplified only by addition of 0.04% BSA + 2% formamide + 0.1 M trehalose and amplification of PAL was possible with addition of 2% of DMSO. Besides optimization of expression assays, the results show that additives can act independently or synergistically.

Differential gene expression profiling through transcriptome approach of Saccharum spontaneum L. under low temperature stress reveals genes potentially involved in cold acclimation

Sugarcane (Saccharum sp.) is predominantly grown in both tropics and subtropics in India, and the subtropics alone contribute more than half of sugarcane production. Sugarcane active growth period in subtropics is restricted to 8-9 months mainly due to winter's low temperature stress prevailing during November to February every year. Being a commercial crop, tolerance to low temperature is important in sugarcane improvement programs. Development of cold tolerant sugarcane varieties require a deep knowledge on molecular mechanism naturally adapted by cold tolerant genotypes during low temperature stress. To understand gene regulation under low temperature stress, control and stressed (10 °C, 24 h) leaf samples of cold tolerant S. spontaneum IND 00-1037 collected from high altitude region in Arunachal Pradesh were used for transcriptome analysis using the Illumina NextSeq 500 platform with paired-end sequencing method. Raw reads of 5.1 GB (control) and 5.3 GB (stressed) obtained were assembled using trinity and annotated with UNIPROT, KEGG, GO, COG and SUCEST databases, and transcriptome was validated using qRT-PCR. The differential gene expression (DGE) analysis showed that 2583 genes were upregulated and 3302 genes were down-regulated upon low temperature stress. A total of 170 cold responsive transcriptional factors belonging to 30 families were differentially regulated. CBF6 (C-binding factor), a DNA binding transcriptional activation protein associated with cold acclimation and freezing tolerance was differentially upregulated. Many low temperature responsive genes involved in various metabolic pathways, viz. cold sensing through membrane fluidity, calcium and lipid signaling genes, MAP kinases, phytohormone signaling and biosynthetic genes, antioxidative enzymes, membrane and cellular stabilizing genes, genes involved in biosynthesis of polyunsaturated fatty acids, chaperones, LEA proteins, soluble sugars, osmoprotectants, lignin and pectin biosynthetic genes were also differentially upregulated. Potential cold responsive genes and transcriptional factors involved in cold tolerance mechanism in cold tolerant S. spontaneum IND 00-1037 were identified. Together, this study provides insights into the cold tolerance to low temperature stress in S. spontaneum, thus opening applications in the genetic improvement of cold stress tolerance in sugarcane.

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.

Host-Plant Resistance to Spittlebugs: A Life-Table Study With Mahanarva fimbriolata (Hemipera: Cercopidae) in Sugarcane Genotypes

Understanding how host-crop genotypes affect the life history of insect pests is important for developing and using varietal resistance as a pest control measure. Here we determined how wild and cultivated sugarcane genotypes affect key life-history traits and the demographic performance of the root spittlebug, Mahanarva fimbriolata (Stål; Hemiptera: Cercopidae), a major pest of sugarcane, which produces most of the world's sugar and bioethanol. In the greenhouse, plants of four sugarcane genotypes (two wild and two cultivated) were infested with newly-ecloded spittlebug nymphs. A longitudinal life-table experiment was conducted recording developmental time, survival, reproductive output, and longevity of the released spittlebug nymphs. One of the promising sugarcane genotypes was the accession IM76-229 (Saccharum robustum Brandes & Jesw. ex Grassl), which allowed only 20% of nymph survival to adulthood and reduced the spittlebug longevity by 10 d. Such effects on these life-history traits led to a negative intrinsic rate of population growth (rm = -0.002) of the insects, indicating that the test population would eventually go extinct on this sugarcane genotype. An opposite trend (i.e., population growth) was obtained for the insects raised on Saccharum spp. (SP81-3250 and SP80-1816 cultivars) and Erianthus arundinaceus (Retz) Jeswiet (Kawandang genotype). These results show that S. robustum IM76-229 genotype has resistance traits that drastically reduce spittlebug population growth and indicate that it can be used in breeding programs aiming to develop sugarcane cultivars with resistance to spittlebugs.

Association of variation in the sugarcane transcriptome with sugar content

Background

Sugarcane is a major crop of the tropics cultivated mainly for its high sucrose content. The crop is genetically less explored due to its complex polyploid genome. Sucrose synthesis and accumulation are complex processes influenced by physiological, biochemical and genetic factors, and the growth environment. The recent focus on the crop for fibre and biofuel has led to a renewed interest on understanding the molecular basis of sucrose and biomass traits. This transcriptome study aimed to identify genes that are associated with and differentially regulated during sucrose synthesis and accumulation in the mature stage of sugarcane. Patterns of gene expression in high and low sugar genotypes as well as mature and immature culm tissues were studied using RNA-Seq of culm transcriptomes.

Results

In this study, 28 RNA-Seq libraries from 14 genotypes of sugarcane differing in their sucrose content were used for studying the transcriptional basis of sucrose accumulation. Differential gene expression studies were performed using SoGI (Saccharum officinarum Gene Index, 3.0), SAS (sugarcane assembled sequences) of sugarcane EST database (SUCEST) and SUGIT, a sugarcane Iso-Seq transcriptome database. In total, about 34,476 genes were found to be differentially expressed between high and low sugar genotypes with the SoGI database, 20,487 genes with the SAS database and 18,543 genes with the SUGIT database at FDR < 0.01, using the Baggerley’s test. Further, differential gene expression analyses were conducted between immature (top) and mature (bottom) tissues of the culm. The DEGs were functionally annotated using GO classification and the genes consistently associated with sucrose accumulation were identified.

Conclusions

The large number of DEGs may be due to the large number of genes that influence sucrose content or are regulated by sucrose content. These results indicate that apart from being a primary metabolite and storage and transport sugar, sucrose may serve as a signalling molecule that regulates many aspects of growth and development in sugarcane. Further studies are needed to confirm if sucrose regulates the expression of the identified DEGs or vice versa. The DEGs identified in this study may lead to identification of genes/pathways regulating sucrose accumulation and/or regulated by sucrose levels in sugarcane. We propose identifying the master regulators of sucrose if any in the future.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4302-5) contains supplementary material, which is available to authorized users.

Luxurious Nitrogen Fertilization of Two Sugar Cane Genotypes Contrasting for Lignin Composition Causes Changes in the Stem Proteome Related to Carbon, Nitrogen, and Oxidant Metabolism but Does Not Alter Lignin Content

Sugar cane is an important crop for sugar and biofuel production. Its lignocellulosic biomass represents a promising option as feedstock for second-generation ethanol production. Nitrogen fertilization can affect differently tissues and its biopolymers, including the cell-wall polysaccharides and lignin. Lignin content and composition are the most important factors associated with biomass recalcitrance to convert cell-wall polysaccharides into fermentable sugars. Thus it is important to understand the metabolic relationship between nitrogen fertilization and lignin in this feedstock. In this study, a large-scale proteomics approach based on GeLC-MS/MS was employed to identify and relatively quantify proteins differently accumulated in two contrasting genotypes for lignin composition after excessive nitrogen fertilization. From the ∼1000 nonredundant proteins identified, 28 and 177 were differentially accumulated in response to nitrogen from IACSP04-065 and IACSP04-627 lines, respectively. These proteins were associated with several functional categories, including carbon metabolism, amino acid metabolism, protein turnover, and oxidative stress. Although nitrogen fertilization has not changed lignin content, phenolic acids and lignin composition were changed in both species but not in the same way. Sucrose and reducing sugars increased in plants of the genotype IACSP04-065 receiving nitrogen.

Transcriptome analysis of genes involved in secondary cell wall biosynthesis in developing internodes of Miscanthus lutarioriparius

Miscanthus is a promising lignocellulosic bioenergy crop for bioethanol production. To identify candidate genes and regulation networks involved in secondary cell wall (SCW) development in Miscanthus, we performed de novo transcriptome analysis of a developing internode. According to the histological and in-situ histochemical analysis, an elongating internode of M. lutarioriparius can be divided into three distinct segments, the upper internode (UI), middle internode (MI) and basal internode (BI), each representing a different stage of SCW development. The transcriptome analysis generated approximately 300 million clean reads, which were de novo assembled into 79,705 unigenes. Nearly 65% of unigenes was annotated in seven public databases. Comparative profiling among the UI, MI and BI revealed four distinct clusters. Moreover, detailed expression profiling was analyzed for gene families and transcription factors (TFs) involved in SCW biosynthesis, assembly and modification. Based on the co-expression patterns, putative regulatory networks between TFs and SCW-associated genes were constructed. The work provided the first transcriptome analysis of SCW development in M. lutarioriparius. The results obtained provide novel insights into the biosynthesis and regulation of SCW in Miscanthus. In addition, the genes identified represent good candidates for further functional studies to unravel their roles in SCW biosynthesis and modification.

Combining transcription factor binding affinities with open-chromatin data for accurate gene expression prediction

The binding and contribution of transcription factors (TF) to cell specific gene expression is often deduced from open-chromatin measurements to avoid costly TF ChIP-seq assays. Thus, it is important to develop computational methods for accurate TF binding prediction in open-chromatin regions (OCRs). Here, we report a novel segmentation-based method, TEPIC, to predict TF binding by combining sets of OCRs with position weight matrices. TEPIC can be applied to various open-chromatin data, e.g. DNaseI-seq and NOMe-seq. Additionally, Histone-Marks (HMs) can be used to identify candidate TF binding sites. TEPIC computes TF affinities and uses open-chromatin/HM signal intensity as quantitative measures of TF binding strength. Using machine learning, we find low affinity binding sites to improve our ability to explain gene expression variability compared to the standard presence/absence classification of binding sites. Further, we show that both footprints and peaks capture essential TF binding events and lead to a good prediction performance. In our application, gene-based scores computed by TEPIC with one open-chromatin assay nearly reach the quality of several TF ChIP-seq data sets. Finally, these scores correctly predict known transcriptional regulators as illustrated by the application to novel DNaseI-seq and NOMe-seq data for primary human hepatocytes and CD4+ T-cells, respectively.

Characterization of sugarcane (Saccharum spp.) leaf senescence: implications for biofuel production

BACKGROUND

Second-generation ethanol (2G-bioethanol) uses lignocellulosic feedstocks for ethanol production. Sugarcane is one among the most suitable crops for biofuel production. Its juice is extracted for sugar production, while sugarcane bagasse, straw, and senescing leaves are considered industrial waste. Senescence is the age-dependent deterioration of plant cells, ultimately leading to cell death and completion of the plant life cycle. Because senescing leaves may also be used for biofuel production, understanding the process of natural senescence, including remobilization of nutrients and its effect on cell walls can provide useful information for 2G-bioethanol production from sugarcane leaves.

RESULTS

The natural senescence process in leaves of the commercial sugarcane cultivar RB867515 was investigated. Senescence was characterized by strong reduction in photosynthetic pigments content, remobilization of the nutrients N, P, K, B, Cu, Fe, and Zn, and accumulation of Ca, S, Mg, B, Mn, and Al. No significant changes in the cell-wall composition occurred, and only small changes in the expression of cell wall-related genes were observed, suggesting that cell walls are preserved during senescence. Senescence-marker genes, such as SAG12-like and XET-like genes, were also identified in sugarcane and found to be highly expressed.

CONCLUSIONS

Our study on nutrient remobilization under senescence in a vigorous sugarcane cultivar can contribute to the understanding on how nutrient balance in a high-yielding crop is achieved. In general, neutral monosaccharide profile did not change significantly with leaf senescence, suggesting that senescing leaves of sugarcane can be as a feedstock for biofuel production using pretreatments established for non-senescing leaves without additional efforts. Based on our findings, the potential biotechnological applications for the improvement of sugarcane cultivars are discussed.