Abstract

Sorghum (Sorghum bicolor [L.] Moench) has excellent agronomic traits and biological properties, such as heat and drought-tolerance. It is a C4 grass and potential bioenergy-producing plant, which makes it an important crop worldwide. With the sorghum genome sequence released, it is essential to establish a sorghum functional genomics data mining platform. We collected genomic data and some functional annotations to construct a sorghum functional genomics database (SorghumFDB). SorghumFDB integrated knowledge of sorghum gene family classifications (transcription regulators/factors, carbohydrate-active enzymes, protein kinases, ubiquitins, cytochrome P450, monolignol biosynthesis related enzymes, R-genes and organelle-genes), detailed gene annotations, miRNA and target gene information, orthologous pairs in the model plants Arabidopsis, rice and maize, gene loci conversions and a genome browser. We further constructed a dynamic network of multidimensional biological relationships, comprised of the co-expression data, protein–protein interactions and miRNA-target pairs. We took effective measures to combine the network, gene set enrichment and motif analyses to determine the key regulators that participate in related metabolic pathways, such as the lignin pathway, which is a major biological process in bioenergy-producing plants.

Database URL:http://structuralbiology.cau.edu.cn/sorghum/index.html.

Introduction

Sorghum originated from Africa and is the 5th major cereal in terms of production and acreage throughout the world. It occupies 8 million hectares of farmland and provides food, feed, fiber and fuel. According to the Consultative Group for International Agricultural Research’s report in 2010, the USA was the world’s largest producer of sorghum (8.8 million metric tons annually), followed by India (7.0), Mexico (6.9), Nigeria (4.8) and Argentina (3.6). The excellent agronomic traits and biological properties, such as heat and drought-tolerant, of sorghum, a C4 grass and potential bioenergy-producing plant, make it an important crop around the world. Sorghum can be classified into several landraces with different levels of genetic diversity. Among them, three uses, grain, forage and sweet sorghum, based on work by Chromatin, Inc. (http://www.chromatininc.com/), attract a great number of sorghum experts. Grain sorghum has a high content of starch and can supply nutrition to humans and animals; forage sorghum is highly digestible and used as animal feed or as a cover crop; Sweet sorghum has a high content of sugar. These three kinds of sorghum can be used as raw materials for the industrial production of ethanol or other specialty chemicals.

In addition to maize and millet, sorghum is a common bioenergy-producing plant, as well as a C4 plant. Biomass is an organic carrier that transfers solar energy to biomass energy through photosynthesis (1, 2). Bioenergy-producing plants have characteristics like stress resistance, large size, quick growth and high productivity, and also higher contents of cellulose, hemicellulose, lignin, sugar, starch and grease. Sorghum is considered a source of biomass for the production of bioenergy and biofuels, as well as for use in the chemical industry (3). Nonetheless, the depth of research and yields in sorghum were more limited than in other grain crops, such as maize, rice and wheat. As an important bioenergy-producing plant, sorghum provides large amount of cellulose, which can be transformed to bioethanol through biological, physical and chemical methods (4, 5). However, because the mixture of lignin and cellulose reduces the conversion efficiency, lowering the lignin content or changing the lignin structure can enhance biomass availability (1, 6). At present, bioengineering methods have been applied to changing the activity of related enzymes, such as cinnamyl alcohol dehydrogenase (CAD), caffeic acid O-methyltransferase and ferulic acid/coniferaldehyde/coniferyl alcohol 5-hydroxylase, to alter the content or structure of lignin, making it more easily degradable (6, 7). However, the whole lignin pathway is still unclear. Discoveries and modifications of suitable enzymes cannot only accelerate biomass usage, but also enhance plant growth (2). Thus, the identification and functional analyses of related regulatory genes through the whole genome network may aid studies on sorghum.

The Sorghum bicolor BTx623 genome sequencing result was published in 2009 (8). Fundamental annotations of gene structures and other information on that assembly were later improved (version 1.0, 1.4 and 2.1), with more accurate parameters, comprehensive methods and a newly integrated dataset. Recent studies have made great progress in the bioenergy production, genetic variation, regulatory factors and metabolic pathways of sorghum under the sequenced genome’s background. For example, some synthesis databases have added sorghum information, including Phytozome (9), Gramene (10), National Center for Biotechnology Information (NCBI), PLAZA (11) and PlantsDB (12). MOROKOSHI (13), a sorghum transcriptome database, integrated functional annotations and used specific RNA-seq data to construct a co-expression network and further study expression profile variations.

Since the development of the microarray and next generation sequencing technology, more and more transcriptome data has become available. Prior to December 2015, the Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/) in NCBI had collected 17 series, 16 platforms and 177 samples for sorghum, including GSE50464 (14), GSE54705 (15) and GSE49879 (16), and their data types covered mRNA-seq, ncRNA-seq, microarrays and others. GSE49879 contains 78 samples of microarray data, which was from six genotypes (AR2400, Atlas, Fremont, PI152611 and PI455230) and four tissues (leaf, root, shoot and stem). Thus, it was possible to construct a whole genome co-expression network. In addition, the predicted sorghum protein–protein interactions (PPIs) from the experimentally validated plant data were also useful for predicting sorghum gene function annotations.

Due to the information provided by related databases and public papers, it was necessary to build a sorghum database for research workers using a more comprehensive search criteria and analysis, similar to the popular single species functional genomic databases in agricultural, such as TIGR, MaizeGDB and SIFGD. Driven by this need, we built a comprehensive platform for the genome functional annotation of sorghum, which was named as the sorghum genomics functional database (SorghumFDB). It contains eight gene family categories, super families for transcription factors/regulators (TFs/TRs), carbohydrate-active enzymes (CAZymes), protein kinases (PKs), ubiquitins (UBs), cytochrome P450 members (CYPs), monolignol biosynthesis (MBs) related protein coding genes, R-genes and organelle-genes. In addition, detailed gene annotations, miRNA and target mRNA information, orthologous relationships with Arabidopsis, maize and rice, gene loci conversions and a genome browser, as well as some analysis tools, such as gene set enrichment analysis (GSEA), motif significance analysis and pattern set tools. Users can visit SorghumFDB at the URL http://structuralbiology.cau.edu.cn/sorghum/. With the accessible information and online tools mentioned, we hope that SorghumFDB will aid sorghum functional genomics analyses and become useful for sorghum and other bioenergy plant-related studies.

Data source

The SorghumFDB was mainly built on version 2.1, which contained 33 032 genes, 39 441 transcripts and other annotation information, including the eukaryote-specific version of the Clusters of Orthologous Groups (KOG) (17) and gene ontology (GO) (18). The NCBI version of the sorghum organelle genome annotation was composed of NC_008602.1 for the chloroplast and NC_008360.1 for the mitochondrial genome. We collected miRNAs from miRBase (19) and PMRD (20). Single nucleotide polymorphisms (SNPs), expressed sequence tags (ESTs) and motifs located in the nucleic genome sequence may regulate gene expression, directly or indirectly. Sorghum SNP data were collected from SorGSD (http://sorgsd.big.ac.cn/snp/index.jsp). EST and repeats were downloaded from Phytozome v10. Combining text-mining technology with data from plantCARE (21), PLACE (22) and AthMap (23), we collected 930 motifs with publication annotations. The orthologous pairs between sorghum and its related species (rice and maize) were derived from Gramene (10). We focused mainly on TFs/TRs and seven other super families, CAZymes, PKs, UBs, CYPs, MBs, R-genes and organelle-genes. The TFs/TRs were collected from PlantTFDB (24) and PlnTFDB (25). The cytochrome P450 family members were collected from Nelson et al. (26). The ubiquitin families were downloaded from UUCD (27). The monolignol biosynthesis families and R-genes were collected by text-mining from Shakoor et al. (16) and Mace et al. (28), respectively. To strengthen the functional annotations, we provided GO (18), Uniprot, Panther (29) and KOG (17) annotations for most of the genes. To further analyse the relationships between genes or protein pairs, the experimentally validated plant PPI data were collected from databases, IntAct (30), Biogrid (31), iMEX (32), TAIR (33), BAR (34), CCSB (35) and published papers (36). Because of the lack of sorghum PPI data, we first collected data from maize and rice, which have close relative relationships with sorghum, but the resulting gene pairs were too few to construct a network. Then, we chose the model plant Arabidopsis and adapted more stringent filter criteria. Together with the orthologous pairs this data was used to build the predicted PPI network. The transcriptome data GSE50464 (14) and GSE54705 (15) were used to create the expression profile tendency chart, and GSE49879 (16) was selected to be used to construct the co-expression network.

Construction

Functional annotation

In the gene detail page of the SorghumFDB, there are many functional module annotations, such as KOG (17), Panther (29), domain, SNP, GO (18), Uniprot and pathway. The UniProt database has abundant information and extensive protein resources. We determined the Uniprot IDs and annotations for the genes. Additionally, the Panther (29) classification system was designed to classify proteins based on family, pathway, molecular function and biological process. We downloaded the S. bicolor information from Panther (29) and incorporated it into our gene detail. The KOG (17) annotation was collected from Phytozome and linked to the NCBI conserved domain database, which consists of a collection of well-annotated multiple sequence alignment models for ancient domains and full-length proteins. The protein domains were predicted using the PfamScan software from the Pfam database (37). SNP sites located either in the 2 kb upstream of the gene transcript start site or gene body region could be linked to SorGDB (http://sorgsd.big.ac.cn/snp/index.jsp), a collection of S. bicolor SNPs, and visualized by Gbrowse. Different databases also use specialized gene identifiers or distinct versions to define the same gene sequence. Therefore, we transformed other gene names to the uniform version 2.1 for a broader annotation and user convenience. The KEGG pathway annotation of sorghum proteins contained 131 kinds of metabolic processes, but 433 enzymes could not be mapped to version 2.1. Then, we use the BLASTP algorithm to compare the protein sequences to version 2.1, while considering whether the two sequences have the same domains predicted by Pfam (37). Finally, we identified 4399 enzymes and 131 pathways from KEGG. PlantCyc (38) provides a broad network of plant metabolic pathway databases that contain curated information from the literature and a computational analysis of the genes, enzymes, compounds, reactions and pathways involved in plant primary and secondary metabolism. We downloaded the SorghumBicolorCyc 3.0 (38) dataset and gained 3377 annotated genes and 535 pathways.

Gene family classification

Because of the limited functional annotations, it is essential to predict potential gene functions for sorghum through comparative genomics and related algorithms. At present, we mainly perform the super family classifications for eight gene categories (39, 40), TFs/TRs, CAZymes, PKs, UBs, CYPs, MBs enzymes, R-genes and organelle-genes, which perform remarkable functions in biological processes. In addition to the eight important super families, we will add new gene families constantly as the family information is published, as revealed by new software and even for specific user requirements. Sorghum is a C4 model plant that can use solar energy with high-efficiency during photosynthesis. The chloroplast is part of main reaction site for C4 photosynthesis. Thus, integrating related carbohydrate-active enzymes and organelle-genes will aid research on the engineering of C4 photosynthesis and in C3 crops (41). In particular, the enzymes from the monolignol biosynthesis gene family have important functions in the lignin pathway. CYPs act on the production of structural components, light harvesting and hormone biosynthesis (26). They can be coupled directly to the photosynthetic energy output to obtain an environmentally friendly production of complex chemical compounds. In addition, they also participate in the biosynthesis of physiologically important compounds, such as fatty acids and steroid hormones (42). TFs/TRs, PKs and UBs participate in many signal transduction processes corresponding to biotic and abiotic stress, such as drought stress (43–45). Disease resistance is the main target of sorghum genetic improvement and contributes to high yield and quality. With the development of sorghum genome research, a large number of sorghum resistance genes have been found and located on the genetic linkage map (46).

We have three main analysis strategies: (i) comparative genomics based on the BLASTP algorithm and Inparanoid (47) software, (ii) annotation data collected directly from corresponding databases or published papers and (iii) predictions using certain software, such as iTAK (48) with the Hidden Markov Model (HMM) constructed by specific protein domain rules.

miRNA

We have collected 223 mature miRNAs with 241 precursor sequences (-5p and -3p) in miRBase (19) Release 21. Still 173 miRNAs in PMRD (20) and only three mature miRNA (sbi-miR165, sbi-miR169l and sbi-miR157) sequences cannot be found in miRBase. We used the GMAP version 2015-11-20 (49) software map to try and identify the three precursor sequences in the sorghum genome, but only sbi-miR169l matched. The locational information on the chromosome was well annotated. We also predicted the target mRNAs of miRNA using psRNAtarget (50), a plant small RNA target prediction toolkit using a proven scoring scheme that calculates unpaired energy (UPE) online. We obtained the expression profiles of 82 miRNAs using the miRNA-seq data GSE32458 from PNRD and 47 miRNAs from the microarray data series GSE49879 to determine the expression trends.

Network

The multidimensional network, including PPIs, co-expression relationships between genes (including miRNAs) and miRNA–target pairs (Figure 2E). The sorghum genome-wide network could help us understand the relationships between molecules.

(1) Predicted PPI network

The PPI network was based on all of the available experimentally verified PPI pairs in Arabidopsis, maize and rice, as well as sorghum, from IntAct (30), Biogrid (31), iMEX (32), TAIR (33), BAR (34), CCSB (35) and eTRAIN (36) owing to the deficiency in the predicted sorghum PPI network (Table 1). Although the relationship between Arabidopsis and sorghum is not close, the PPI pairs are too limited to construct the network only using maize and rice. Therefore, we chose relatively strict parameters, with up to 60% bootstrap support, produced by Inparanoid (47) software to determine orthologs. We used the comparative genomics method to map the PPIs of other plants to sorghum PPIs.

Table 1.

Summary of the data sources

CategoryDescriptionSource/methodDetails
Genomeversion 2.1Phytozome v1033032
version 1.4Phytozome v927608
version 1plantsDB34567
OrganelleNCBI194
miRNAmiRBase, PMRD242
Sequence elementSNPSorGSB∼62 million
Repeat maskerPhytozome v10639245
ESTPhytozome v10∼404650
MotifplantCARE, PLACE, AthMap930
TranscriptomeRNA-seqNCBI/GEOGSE50464, GSE54709, SRP008469
MicroarrayNCBI/GEOGSE49879
Interactionprotein–protein interactionBiogrid, iMEX, Intact TAIR, CCSB, ETRAIN3871 genes, 11621 pairs
co-expression interactionPCC,MR144,901 positive pairs
136,596 negative pairs
miRNA-targetspsRNAtarget4,376 pairs
Functional annotationKEGG pathwayKEGG4399 genes, 131 pathways
plantCyc pathwayplantCyc3377 genes, 535 pathways
KOGPhytozome v103150
PantherPanther2836
UniprotUniProt32796
GOPhytozome v10, Gramene1177
DomainPfam3832
OrthologArabidopsis111,461 pairs
rice21,176 pairs
maize19,691 pairs
Gene familyTranscription regulators/factorsplantTFDB, plnTFDB2541genes,85 families
Carbohydrate-active enzymesCAZy1166 genes, 6 classes, 99 families
Protein kinasesPCC database, iTAK1253genes,84 families
UbiquitinsUUCD1846 genes,35 families
Cytochrome P450Nelson.et al.176 genes, 81 families
Monolignol biosynthesis enzymesShakoor et al.36 genes,10 families
R-genesMace et al.231 genes,12 classes
Organelle-genesNCBI194 elements, 2 families
CategoryDescriptionSource/methodDetails
Genomeversion 2.1Phytozome v1033032
version 1.4Phytozome v927608
version 1plantsDB34567
OrganelleNCBI194
miRNAmiRBase, PMRD242
Sequence elementSNPSorGSB∼62 million
Repeat maskerPhytozome v10639245
ESTPhytozome v10∼404650
MotifplantCARE, PLACE, AthMap930
TranscriptomeRNA-seqNCBI/GEOGSE50464, GSE54709, SRP008469
MicroarrayNCBI/GEOGSE49879
Interactionprotein–protein interactionBiogrid, iMEX, Intact TAIR, CCSB, ETRAIN3871 genes, 11621 pairs
co-expression interactionPCC,MR144,901 positive pairs
136,596 negative pairs
miRNA-targetspsRNAtarget4,376 pairs
Functional annotationKEGG pathwayKEGG4399 genes, 131 pathways
plantCyc pathwayplantCyc3377 genes, 535 pathways
KOGPhytozome v103150
PantherPanther2836
UniprotUniProt32796
GOPhytozome v10, Gramene1177
DomainPfam3832
OrthologArabidopsis111,461 pairs
rice21,176 pairs
maize19,691 pairs
Gene familyTranscription regulators/factorsplantTFDB, plnTFDB2541genes,85 families
Carbohydrate-active enzymesCAZy1166 genes, 6 classes, 99 families
Protein kinasesPCC database, iTAK1253genes,84 families
UbiquitinsUUCD1846 genes,35 families
Cytochrome P450Nelson.et al.176 genes, 81 families
Monolignol biosynthesis enzymesShakoor et al.36 genes,10 families
R-genesMace et al.231 genes,12 classes
Organelle-genesNCBI194 elements, 2 families
Table 1.

Summary of the data sources

CategoryDescriptionSource/methodDetails
Genomeversion 2.1Phytozome v1033032
version 1.4Phytozome v927608
version 1plantsDB34567
OrganelleNCBI194
miRNAmiRBase, PMRD242
Sequence elementSNPSorGSB∼62 million
Repeat maskerPhytozome v10639245
ESTPhytozome v10∼404650
MotifplantCARE, PLACE, AthMap930
TranscriptomeRNA-seqNCBI/GEOGSE50464, GSE54709, SRP008469
MicroarrayNCBI/GEOGSE49879
Interactionprotein–protein interactionBiogrid, iMEX, Intact TAIR, CCSB, ETRAIN3871 genes, 11621 pairs
co-expression interactionPCC,MR144,901 positive pairs
136,596 negative pairs
miRNA-targetspsRNAtarget4,376 pairs
Functional annotationKEGG pathwayKEGG4399 genes, 131 pathways
plantCyc pathwayplantCyc3377 genes, 535 pathways
KOGPhytozome v103150
PantherPanther2836
UniprotUniProt32796
GOPhytozome v10, Gramene1177
DomainPfam3832
OrthologArabidopsis111,461 pairs
rice21,176 pairs
maize19,691 pairs
Gene familyTranscription regulators/factorsplantTFDB, plnTFDB2541genes,85 families
Carbohydrate-active enzymesCAZy1166 genes, 6 classes, 99 families
Protein kinasesPCC database, iTAK1253genes,84 families
UbiquitinsUUCD1846 genes,35 families
Cytochrome P450Nelson.et al.176 genes, 81 families
Monolignol biosynthesis enzymesShakoor et al.36 genes,10 families
R-genesMace et al.231 genes,12 classes
Organelle-genesNCBI194 elements, 2 families
CategoryDescriptionSource/methodDetails
Genomeversion 2.1Phytozome v1033032
version 1.4Phytozome v927608
version 1plantsDB34567
OrganelleNCBI194
miRNAmiRBase, PMRD242
Sequence elementSNPSorGSB∼62 million
Repeat maskerPhytozome v10639245
ESTPhytozome v10∼404650
MotifplantCARE, PLACE, AthMap930
TranscriptomeRNA-seqNCBI/GEOGSE50464, GSE54709, SRP008469
MicroarrayNCBI/GEOGSE49879
Interactionprotein–protein interactionBiogrid, iMEX, Intact TAIR, CCSB, ETRAIN3871 genes, 11621 pairs
co-expression interactionPCC,MR144,901 positive pairs
136,596 negative pairs
miRNA-targetspsRNAtarget4,376 pairs
Functional annotationKEGG pathwayKEGG4399 genes, 131 pathways
plantCyc pathwayplantCyc3377 genes, 535 pathways
KOGPhytozome v103150
PantherPanther2836
UniprotUniProt32796
GOPhytozome v10, Gramene1177
DomainPfam3832
OrthologArabidopsis111,461 pairs
rice21,176 pairs
maize19,691 pairs
Gene familyTranscription regulators/factorsplantTFDB, plnTFDB2541genes,85 families
Carbohydrate-active enzymesCAZy1166 genes, 6 classes, 99 families
Protein kinasesPCC database, iTAK1253genes,84 families
UbiquitinsUUCD1846 genes,35 families
Cytochrome P450Nelson.et al.176 genes, 81 families
Monolignol biosynthesis enzymesShakoor et al.36 genes,10 families
R-genesMace et al.231 genes,12 classes
Organelle-genesNCBI194 elements, 2 families

(2) Co-expression network

Since the development of next generation sequencing technologies, large amounts of ‘omics’ data, like RNA-seq and microarray, have become available on public platforms (e.g. GEO and SRA). A large number of probes arranged on a microarray can represent the whole or part of a genome. By RNA extraction, followed by reverse transcription to cDNA, and then fluorescence labeling, we obtained the expression quantity through the fluorescent signals. While the RNA-seq technology is much easier, it can reveal unknown transcripts. The RNA-seq dataset displays the gene expression trend in different tissues and after different treatments in our gene detail webpage. We only chose the microarray data GSE49879 (Supplementary Table S1) that contains 78 samples from a combination of four different tissues (shoot, root, leaf and stem), two dissected stem tissues (pith and rind) and six diverse lines (R159, Atlas, Fremont, PI152611, AR2400 and PI455230) among three landraces (sweet, grain and forage) to construct the co-expression network, as well as genotype-specific and tissue-preferential networks. The computational process (51) was described as follows: First, if a gene matched with more than one probe set, then the highest normalized intensity among those probe sets was selected as its expression value. The screened data was maintained for the next step. Second, we filtered the data with a strict cutoff [max–min = 0.58, which is log2(1.5)] to prevent housekeeping genes from making background noise. Third, Pearson’s correlation coefficient (PCC) scores between genes or miRNAs were computed to evaluate their expression relationships. To confirm the correction, we used mutual rank (MR), which showed a higher prediction efficiency than PCC values by taking a geometric average of the PCC rank (52, 53). Finally, we filtered the interaction pairs and divided them into three levels, the first level contained genes having the top three PCCs; the second level contained genes with the MR ≤ 5, and the third level kept the ones with the 5 < MR ≤ 30 (Supplementary Figure S3). For the network, we set a cut-off value of 4.56. The gene is specifically expressed if the expression value is higher than 4.56 in any of the repeats or is not expressed if the expression value is lower than 4.56 in all of the repeats under one condition.

PCC calculates the linear correlation between two variables, X and Y, using the following formula:

Here, X and Y represent the expression quantities of two genes. The i indicates different samples under various conditions or times (Supplementary Figure S2).

The MR values were calculated using the following computational formula:

For two genes, X and Y, selecting all of the pairs related to gene X and ranking the PCC from high to low, the position (starting from 1, with a step size of 1) of the Y gene is represented by a, while b represents the position of X in Y’s PCC list.

(3) miRNA-target network

The miRNA–target pairs were computed by psRNAtarget (50) through base pairing of miRNAs with their complementary mRNA targets. miRNAs play dominant roles in post-transcriptional gene regulation and have been experimentally proven to greatly affect crop plant productivity and quality (54).

(4) Functional annotation for each gene

We retrieved the single gene annotation from Phytozome, and it is displayed by clicking on the nodes in the network. A detailed gene list appears below, containing PCC, MR and PPI sources with links.

(5) Functional module prediction and annotation

After filtering false positive gene pairs based on the MR method, we had 281 497 remaining gene pairs. We used Cfinder (55), a free software based on the Clique Percolation Method, to find overlapping dense groups of nodes in the co-expression network. Using a balance between gene coverage and overlapping rates among modules, we selected the result when k = 5 (Supplementary Figure S4) as the threshold, which meant that a module must have at least five nodes based on the gene coverage degree, community distribution and overlap rate, and then we filtered the communities that contained more than 500 nodes. Finally, we selected 987 modules contained 3954 genes. Each network module was subjected to a GSEA analysis (P-value ≤ 0.05 and false discovery rate ≤ 0.05) to annotate the models.

Tools

Motifs (cis-regulatory elements) are a series of short conserved sequences in the promoter region. They can be recognized and bound by TFs/TR and then participate together in the regulation of downstream genes. With identified motifs and published annotations, we predicted TF/TR binding regions, which could benefit studies on the relationships between TFs and their target genes. In addition to scanning for motif sequences, we offered a Z-score method to calculate the significance of enriched motifs (39). All of the sorghum genes’ promoters have been calculated as background, and users can submit gene sequences or gene lists, in which they can search for related motifs and carry out motif enrichment analyses.

where Nmotifi indicates the number of occurrences of motifi in your promoters, meanmotifi indicates the average occurrences of motifi in these 1000 sets, stdevmotifi indicates their standard deviation and pnorm indicates the distribution function for the normal distribution in the R language (56).

The pattern set tool also uses similar algorithms as in the motif analysis. We included the transcriptome datasets GSE50464, GSE54705, GSE49879 and SRP008469, which contain samples from different tissues, treatments, and with different genotypes as background. Users can set specific expression levels, ranked as high, even and low in each sample, indicating specific expression patterns in a data series. Then, users can obtain the screening gene list based on their settings (51).
(a)
(b)

where x- represents the average expression value of the repetitions for one gene, x= represents the average expression value of all of the samples for one dataset, stdev represents the standard deviation of all of the samples for one dataset; and x-T represents the average expression value of the treated samples, x-C represents the average expression value of the control samples. The (a) indicates a tissue-specific data series, and (b) indicates a treatment-specific data series. Otherwise, the z-score is 0.

Additionally, we took the gene families, GO terms (18), miRNA targets and pathway information from our functional annotations, described before as knowledge background sets. Users can submit their gene lists, derived from high-throughput experiments or other resources, and obtain the biological functions or molecular activities in which the submitted genes were significantly involved based on the p-values of statistical tests and the adjusted p-value after the false discovery rate correction (57).

where N is the total number of sorghum genes (currently 33 032), and n is the number of genes in the query list. K represents the total number of genes in one gene set and k represents the number of overlapped genes. The default P-value is 0.05.

Visualization

Cytoscape (58) can be used to display the gene pairs by nodes and edges. We used different color and shape combinations to represent the nodes and edges. The yellow nodes indicate query genes or miRNAs, and the green nodes indicate the interaction elements. The orange lines link the positive co-expression pairs, while the blue color lines indicate negative co-expression pairs. To further enhance the links from gene nodes to each genes annotation, we added a pop-up window, which displays a short annotation when a gene from the network is clicked on (Figure 2E).

The datasets (GSE50464 and GSE54705) from GEO were changed from the sra to fastq format using the sratool (sratoolkit.2.4.0-1-centos_linux64). We used FastQC (fastqc_v0.10.1) to determine the data quality. Sequence reads were aligned to the sorghum genome through tophat (tophat-2.0.12.Linux_x86_64) and the annotation file was referred to using cufflinks (cufflinks-2.2.1.Linux_x86_64), which computed the fragments per kilobase of transcript per million mapped reads with the default parameters (59, 60) (Supplementary Table S5). The microarray dataset (GSE49879) probes corresponded to sorghum gene version 2.1 ID. After computing the RNA-seq data series and microarray data series, we calculated three expression matrices for different genotypes and tissues. To view each gene expression pattern, we used the open-flesh-chart (http://www.highcharts.com/) to draw the bar graph. The expression value was also included behind the bar graph (Figure 2C).

Gbrowse (v2.0) is a combination of databases and interactive web pages for manipulating and displaying annotations (like gene structure and location) on genomes. We have shown seven integrated genomic features, nuclear genome coding genes, miRNAs, ESTs, organelle genome information, repeat elements, mRNA-seq data and miRNA-seq data, in Gbrowse. Users can select the features that they are interested in, which will be displayed on the website (Supplementary Figure S1).

Results

SorghumFDB was constructed under the LAMP (Linux + Apache + Mysql + PHP) environment. The dataset was uploaded to MySQL, and the web interface was built by PHP and Perl (http://www.perl.org) scripts. The SorghumFDB (Figure 1) has eight main parts. In the (a) ‘Search’ part, we have four special searching boxes, gene ID search, KEGG pathway, gene ID switch and ortholog search, as well as a keywords search function located on the navigation bar. The orthologous relationships define 111 461 pairs between 24 651 sorghum genes and 21 070 Arabidopsis genes, 19 691 pairs between 16 233 sorghum genes and 15 113 rice genes, and 21 176 pairs between 17 485 sorghum genes and 16 904 maize genes. The metabolic pathway includes 4399 genes in 131 KEGG pathways and 3377 genes in 535 PlantCyc (38) pathways. For the functional annotation, there are 3150 KOGs (17), 2836 Panthers (29), 32 796 Uniprot, 1177 GO (18) terms and 3832 domains. The (b) ‘Gene family’ part, contains eight commonly used families, 2541 TRs/TFs, 1253 PKs, 1166 CAZymes, 176 CYPs, 1846 UBs, 36 MBs, 231 R-genes and 194 organelle elements. The (c) ‘miRNA’ part has 242 miRNAs and 1742 predicted target genes. In addition, owing to the limited amount of research, some of the genes could not be annotated. Thus, we offer the (d) ‘Network’ part to identify genes that have similar functions as the query genes or related up-stream or down-stream genes. In conclusion, the multi-dimensional network contains 11 621 predicted PPIs with 3871 proteins, a co-expression network with 144 901 positive and 136 596 negative correlation pairs between 17 665 genes or miRNAs, and miRNA-target network with 4376 pairs. In other words, the network included 240 miRNAs and 21 494 genes, which covered 65.1% of the whole genome (33 032 genes). Then, the co-expression network was split into 987 modules containing 3954 genes by Cfinder (55). For the (e) ‘Tool’ part, four useful analysis tools, GSEA, motif analysis, BLAST and pattern set, that enabled a statistical analysis of differentially expressed genes (RNA-seq or microarray), modification-enriched (TF binding, histone modification and DNA methylation) genes and network genes, which could illustrate the related biological phenomena. In the (f) ‘Generic genome browser’ (Gbrowse) part, the locational information on the seven tracks, including sequence elements and high throughput data, are accessible (Table 1). All of the information was integrated into the gene or miRNA detail webpage. The bar graphs of expression profiles were also added to the detail pages. Additionally, the (g) ‘Homepage’ and (h) ‘Help’ webpages will help users to understand the key points, construction and accessibility of SorghumFDB.

Overview of the sorghum database. The sorghum database has eight parts in total. The eight parts are shown in the picture as yellow rectangles, and the blue rectangles are the detailed classifications or annotations for the previous level. Finally, each part of the detailed information is linked to the gene annotation.
Figure 1.

Overview of the sorghum database. The sorghum database has eight parts in total. The eight parts are shown in the picture as yellow rectangles, and the blue rectangles are the detailed classifications or annotations for the previous level. Finally, each part of the detailed information is linked to the gene annotation.

We established some ways to introduce the uses of the SorghumFDB through the following: Sb02g024190 and Sb04g005950 (61) are two cinnamyl alcohol dehydrogenases (CADs), which are key enzymes in the lignin pathway that affect the biosynthesis of lignin monolignols (62). The search tool transformed the old version ID, Sb02g024190, to Sobic.002G195400 (the sorghum version 2.1). After that, a gene details search displayed integrated information on the gene. For example, (i) the structure and location of the gene (Gbrowse); (ii) the gene, transcript, coding DNA sequence and protein sequence; (iii) the functional domain predicted by Pfam (37), such as the ‘ADH_N’ or ‘ADH_zinc_N’ domains (Figure 2A); (iv) the GO (18) term annotations, such as ‘zinc ion binding’ or ‘plant-type hypersensitive response’; (v) the Uniprot ID ‘C5XC49’ and Panther (29) ID ‘PTHR11695:SF563’; and the (vi) orthologous pairs in Arabidopsis, rice and maize (Figure 2B). In the bar graphs it was evident that the gene expression level was significantly high in roots in the data series of GSE50464 and was still expressed after the ABA treatment. However, in the data series of GSE54705, the genes were more highly expressed in the genotype of Ch17 than in other genotypes, such as CK60 and TX623 (Figure 2C).

Detail information on the sorghum webpage. By searching with the gene Sobic.002G195400, we retrieved (A) the basic information on gene structure, location and sequence; (B) the GO, Uniprot, Panther and gene family annotations; (C) the expression patterns based on the data series GSE50464, GSE54705 and GSE49879, for different tissues, genotypes and treatments. (D) The secondary structure and complementary base pairing of miRNA sbi-miR396a. (E) The search results of Sobic.004G317000’s network and module 5.
Figure 2.

Detail information on the sorghum webpage. By searching with the gene Sobic.002G195400, we retrieved (A) the basic information on gene structure, location and sequence; (B) the GO, Uniprot, Panther and gene family annotations; (C) the expression patterns based on the data series GSE50464, GSE54705 and GSE49879, for different tissues, genotypes and treatments. (D) The secondary structure and complementary base pairing of miRNA sbi-miR396a. (E) The search results of Sobic.004G317000’s network and module 5.

Apart from the single gene annotation, searches for correlated relationships with other genes was provided. Here, we took Sobic.004G317000, a GRF transcript factor as an example. When we tried to analyse the network of the query gene, we could choose from several kinds of network, such as positive co-expression, negative co-expression, miRNA-targets and PPI networks (Figure 2E). To further analyse the expression specificity in different genotypes and tissues, we selected the sweet sorghum Atlas and stem to construct dynamic networks. In addition to the relationships with genes, relationships between functional modules could be displayed. For instance, Sobic.004G129600 was not well annotated using other bioinformatics methods, but after the module analysis several overlapping modules contained it (Figure 2E). Thus, we could hypothesize the function of the query genes based on the modules.

In the gene Sobic.004G317000’s network, the pink line represented the gene targeted by the miRNA sbi-miR396 family. In addition, the miRNA detail webpage was included in the database. For example, the miRNA sbi-miR396a detail page displayed the miRNA sequence of mature, precursor and miRNA secondary structures, as well as the expression bar chart (Figure 2D).

Finally, some analysis tools have been added to our database, such as BLAST, gene set enrichment (57) and motif significance analyses, as well as a pattern set tool.

Discussion

SorghumFDB contains single gene functional annotations and searches, as well as multiple gene analyses. Compared with SorghumFDB, multiple species synthesis databases, such as Phytozome and Gramene, may elaborate less data than that of single species functional genomics databases that possess analyses tools and gene family classifications involving species specificities and flexible updates. For the popular single species functional genomic databases in agricultural, such as TIGR, MaizeGDB and SIFGD, that integrates information from previous works. Except for basic functions that were referred to in other single species’ functional databases, we added new characteristics aimed at producing good agronomic traits and analyses tools.

SorghumFDB should be applied as a practical analysis tool. Biomass is derived from solar energy and a kind of renewable energy source. Lignocellulose biomass is the most important component. S. bicolor is a popular bioenergy plant that can relieve traditional fossil fuels shortages through the biosynthesis of cellulose. Sugarcane and corn, which are known to have a high biomass of sugar and starch, have become the first generation of biofuels, while sorghum and cotton became the new generation of bioenergy plants, owing to their abundant lignocellulosic biomass (63). Lignin is crucial for the structural integrity of the cell wall, and the stiffness and strength of the stem. However, a low lignin content is the prerequisite for enzymatic degradation of cellulose and hemicellulose, which can be further converted to the biofuel cellulosic ethanol. Nowadays, scientists dream of modifying lignin’s composition and content using bioengineering methods that improve the conversion efficiency from cellulose to ethyl alcohol (2).

Lignin is a kind of terrestrial biopolymer which rank only second to cellulose. Lignin is the product of phenylpropanoid pathway and contains primarily three types of subunits which are p-hydroxyphenyl (H), guaiacyl (G) and syringyl (S) lignin (6). Lignin biosynthesis pathway is mainly lignin monomer synthesis. P-coumaryl, coniferyl and sinapyl alcohols which experienced different degree of methoxylation (64, 65) are three dominant lignin monomers. After the synthesis of monolignol, they are transferred to deposition sites. In the end, through dehydrogenation and polymerization, they form lignin in the cell wall. The monolignol biosynthetic pathway was detailed introduced below in (Figure 3C) and the key enzymes were list in Supplementary Table S2 (6).

Network of the CAD gene Sobic.004G071000. (A) The CAD network. The query gene Sobic.004G071000 is highlighted by yellow, the blue line represents the negative co-expression of the conjoined gene, the orange line represents the positive co-expression of the conjoined gene, and the purple line indicates that the two nodes have a PPI relationship. The pink rhombus represents TFs. Here, one is MYB and the other is bHLH. The purple rectangles represent the carbohydrate-active enzyme family. The dark yellow hexagonal gene participated in phenylpropanoid biosynthesis pathway. (B) The transcriptome analysis of all of the CAD co-expression network genes. We selected the genotype type of PI455230 and different tissue, such as seed and leaf, and used cluster 3.0 software to illustrate the expression profile. We also highlighted the important genes of the network with the same colors as in Figure 3A. (C) The lignin pathway illustrates the biosynthesis process of the three kinds of lignin monomers, S-, G- and H-lignin. In the pathway, the related regulation factors, as found in published papers, are labeled with arrows. The up arrows represent up-regulation, while the down arrows represent down-regulation.
Figure 3.

Network of the CAD gene Sobic.004G071000. (A) The CAD network. The query gene Sobic.004G071000 is highlighted by yellow, the blue line represents the negative co-expression of the conjoined gene, the orange line represents the positive co-expression of the conjoined gene, and the purple line indicates that the two nodes have a PPI relationship. The pink rhombus represents TFs. Here, one is MYB and the other is bHLH. The purple rectangles represent the carbohydrate-active enzyme family. The dark yellow hexagonal gene participated in phenylpropanoid biosynthesis pathway. (B) The transcriptome analysis of all of the CAD co-expression network genes. We selected the genotype type of PI455230 and different tissue, such as seed and leaf, and used cluster 3.0 software to illustrate the expression profile. We also highlighted the important genes of the network with the same colors as in Figure 3A. (C) The lignin pathway illustrates the biosynthesis process of the three kinds of lignin monomers, S-, G- and H-lignin. In the pathway, the related regulation factors, as found in published papers, are labeled with arrows. The up arrows represent up-regulation, while the down arrows represent down-regulation.

Another CAD gene (Sobic.004G071000) has two transcripts and the coding protein embed ‘ADH_N’ and ‘ADH_Zinc_N’ domains. At the same time, it was annotated by ‘alcohol dehydrogenase activity’, ‘zinc-dependent’ GO annotations and participated in phenylpropanoid biosynthesis pathway.

The network analysis was our starting point. The construction of a multidimensional network, as introduced above, was based on the microarray data. In the public NCBI platform, the transcriptome data, including microarray and RNA-seq, both had certain accumulation levels. While the differences in experimental designs, analyses processes and the normalization methods make it difficult to combine RNA-seq and microarray datasets. Transcriptome data covering broad biological samples could contribute to the accuracy of the constructed co-expression network. In particular, in large sample size research, microarray data was processed faster and easier, while the larger coverage of various conditions is necessary to detect subtle functional connections (53). The comparison of RNA-seq datasets, which focus on different treatments and tissues, with microarray datasets main results in research on diverse tissues and genotypes. The published database for sorghum transcriptomes, named MOROKOSHI, displays a network based on the RNA-seq dataset. Then, we used the microarray data to build a network. We analysed the CAD gene in both the internal network, which was computed to the microarray data (Figure 3A), and the outside network was based on the RNA-seq data. After a network analysis, the outside network found 19 genes with 23 positive co-expression relationships after transferring version 1 to version 2.1 and filtering for incomplete pairs (Supplementary Figure S5B). Additionally, SorghumFDB produced 45 genes that had relationships with the query genes, including 25 positive co-expressers, 14 negative co-expressers and 6 PPIs.

We applied all 46 genes (Table 2) in the selected network to a GSEA analysis using default parameters (Fisher’s exact test, P-value ≤ 0.05). These genes were significantly involved in lignin biosynthesis pathways, such as flavonoid and phenylpropanoid biosynthesis. Meanwhile, some of the key elements were found, including MYB, bHLH transcription factors and the carbohydrate-active enzyme family (Supplementary Table S3). These elements were previously reported to play important roles in the lignin pathway (62,65). To explain and highlight the key functions of the CAD network, we replaced some gene IDs with the corresponding annotation information from the result of the GSEA analysis using different colors and shapes as prescribed by the Cytoscape software (Figure 3A). For example, the query gene Sobic.004G071000 was replaced with the CAD and highlighted in bright yellow. Meanwhile, the network genes from MOROKOSHI had two overlapping genes, except for the query gene, with our CAD network gene list (Supplementary Figure S5A). The specific genes (17) in the MOROKOSHI network had additional functions, such as CYP75A of the CYPs family, which has an unregulated function as seen in Figure 3C, which shows the GESA annotation analysis (66) (Supplementary Table S3). The co-expression network between the microarray and RNA-seq data may be complementary.

Table 2.

Details of the CAD gene Sobic.004G071000’s network

Protein AProtein BPCCMRRelationship
Sobic.004G273800Sobic.004G0710000.7420.49Positive
Sobic.001G316400Sobic.004G0710000.749.38Positive
Sobic.002G260000Sobic.004G0710000.776.928Positive
Sobic.007G115100Sobic.004G0710000.824.898Positive
Sobic.003G225300Sobic.004G0710000.814.898Positive
Sobic.005G180700Sobic.004G0710000.7720.39Positive
Sobic.003G244100Sobic.004G0710000.7312.72Positive
Sobic.003G343900Sobic.004G0710000.7419.74Positive
Sobic.004G145200Sobic.004G0710000.7519.49Positive
Sobic.004G071000Sobic.004G2203000.834.898Positive
Sobic.004G071000Sobic.001G3173000.6928.61Positive
Sobic.004G071000Sobic.001G4836000.794.898Positive
Sobic.004G071000Sobic.010G0522000.7922.44Positive
Sobic.004G071000Sobic.009G0336000.7727.14Positive
Sobic.004G071000Sobic.010G1924000.728.46Positive
Sobic.004G071000Sobic.007G0473000.843.162Positive
Sobic.004G071000Sobic.010G2096000.7713.49Positive
Sobic.004G071000Sobic.009G0272000.6715.87Positive
Sobic.004G071000Sobic.007G1799000.87.071Positive
Sobic.004G071000Sobic.008G0089000.7328.98Positive
Sobic.004G071000Sobic.006G0820000.7312.04Positive
Sobic.004G071000Sobic.004G2123000.7817.6Positive
Sobic.004G071000Sobic.001G0529000.82.449Positive
Sobic.004G071000Sobic.004G2329000.793Positive
Sobic.004G071000Sobic.010G0307000.775.656Positive
Sobic.004G273800Sobic.001G3164000.7217.54Positive
Sobic.004G273800Sobic.003G2253000.7618.33Positive
Sobic.004G273800Sobic.003G3439000.7321.63Positive
Sobic.004G273800Sobic.004G2203000.7524.49Positive
Sobic.004G273800Sobic.001G4836000.82.828Positive
Sobic.004G273800Sobic.010G2096000.7514.96Positive
Sobic.004G273800Sobic.007G1799000.841Positive
Sobic.004G273800Sobic.004G2329000.749.165Positive
Sobic.001G316400Sobic.003G2253000.748Positive
Sobic.001G316400Sobic.004G2203000.775.196Positive
Sobic.001G316400Sobic.001G4836000.7122.97Positive
Sobic.001G316400Sobic.007G0473000.7416.88Positive
Sobic.001G316400Sobic.004G2123000.7322.8Positive
Sobic.001G316400Sobic.004G2329000.6927.92Positive
Sobic.002G260000Sobic.003G2441000.712.84Positive
Sobic.002G260000Sobic.004G1452000.784.242Positive
Sobic.002G260000Sobic.001G4836000.748.124Positive
Sobic.002G260000Sobic.010G1924000.811Positive
Sobic.002G260000Sobic.009G0272000.6425.9Positive
Sobic.002G260000Sobic.004G2329000.728.366Positive
Sobic.002G260000Sobic.010G0307000.7211.48Positive
Sobic.007G115100Sobic.003G2253000.7729.34Positive
Sobic.007G115100Sobic.003G2441000.7320.97Positive
Sobic.007G115100Sobic.010G0522000.8122.58Positive
Sobic.007G115100Sobic.009G0336000.7928.58Positive
Sobic.007G115100Sobic.007G0473000.8115.71Positive
Sobic.007G115100Sobic.007G1799000.819.899Positive
Sobic.007G115100Sobic.001G0529000.7618.43Positive
Sobic.003G225300Sobic.005G1807000.841.414Positive
Sobic.003G225300Sobic.003G3439000.7810.9Positive
Sobic.003G225300Sobic.004G2203000.7915.49Positive
Sobic.003G225300Sobic.001G4836000.822.236Positive
Sobic.003G225300Sobic.007G0473000.836Positive
Sobic.003G225300Sobic.010G2096000.86.928Positive
Sobic.003G225300Sobic.004G2123000.7822.64Positive
Sobic.003G225300Sobic.010G0307000.775Positive
Sobic.005G180700Sobic.003G2441000.7223.83Positive
Sobic.005G180700Sobic.004G2203000.828.366Positive
Sobic.005G180700Sobic.007G0473000.7923.66Positive
Sobic.005G180700Sobic.010G2096000.842.449Positive
Sobic.003G244100Sobic.007G0473000.7513.26Positive
Sobic.003G244100Sobic.004G2123000.7227Positive
Sobic.003G244100Sobic.010G0307000.6917.97Positive
Sobic.003G343900Sobic.004G1452000.7326.26Positive
Sobic.003G343900Sobic.001G4836000.778.366Positive
Sobic.003G343900Sobic.007G0473000.833.741Positive
Sobic.003G343900Sobic.004G2123000.7622.24Positive
Sobic.004G145200Sobic.010G1924000.768.062Positive
Sobic.004G145200Sobic.007G0473000.816.782Positive
Sobic.004G145200Sobic.008G0089000.7329.98Positive
Sobic.004G220300Sobic.010G0522000.874Positive
Sobic.004G220300Sobic.007G0473000.847.416Positive
Sobic.004G220300Sobic.007G1799000.826.708Positive
Sobic.004G220300Sobic.004G2123000.854.898Positive
Sobic.004G220300Sobic.004G2329000.7127.92Positive
Sobic.001G317300Sobic.006G0820000.6729.98Positive
Sobic.001G317300Sobic.010G0307000.6915.87Positive
Sobic.001G483600Sobic.010G2096000.779.165Positive
Sobic.001G483600Sobic.004G2329000.754.898Positive
Sobic.001G483600Sobic.010G0307000.757.937Positive
Sobic.010G052200Sobic.009G0336000.894Positive
Sobic.010G052200Sobic.007G0473000.912.828Positive
Sobic.010G052200Sobic.004G2123000.8510Positive
Sobic.009G033600Sobic.007G0473000.911.732Positive
Sobic.009G033600Sobic.004G2123000.821.02Positive
Sobic.010G192400Sobic.009G0272000.6614.56Positive
Sobic.010G192400Sobic.010G0307000.756.708Positive
Sobic.007G047300Sobic.004G2123000.921Positive
Sobic.010G209600Sobic.010G0307000.7510.95Positive
Sobic.009G027200Sobic.001G0529000.733.741Positive
Sobic.007G179900Sobic.004G2329000.7318.33Positive
Sobic.008G008900Sobic.006G0820000.748.66Positive
Sobic.006G082000Sobic.001G0529000.746.633Positive
Sobic.007G207400Sobic.004G071000−0.6523.36Negative
Sobic.004G317300Sobic.004G071000−0.717Negative
Sobic.002G031000Sobic.004G071000−0.749.591Negative
Sobic.003G426700Sobic.004G071000−0.7421.09Negative
Sobic.004G071000Sobic.005G005500−0.7324.49Negative
Sobic.004G071000Sobic.009G254500−0.7310.58Negative
Sobic.004G071000Sobic.001G035800−0.7714.49Negative
Sobic.004G071000Sobic.009G191900−0.7229.25Negative
Sobic.004G071000Sobic.004G041500−0.649.433Negative
Sobic.004G071000Sobic.007G142200−0.7114.69Negative
Sobic.004G071000Sobic.010G128200−0.7518.41Negative
Sobic.004G071000Sobic.010G202700−0.79.165Negative
Sobic.004G071000Sobic.010G190700−0.7910.29Negative
Sobic.004G071000Sobic.009G127500−0.714.14Negative
Sobic.004G071000Sobic.010G232100Predicted PPI
Sobic.004G071000Sobic.003G403700Predicted PPI
Sobic.004G071000Sobic.004G049900Predicted PPI
Sobic.004G071000Sobic.004G050000Predicted PPI
Sobic.004G071000Sobic.004G050100Predicted PPI
Sobic.004G071000Sobic.010G239500Predicted PPI
Protein AProtein BPCCMRRelationship
Sobic.004G273800Sobic.004G0710000.7420.49Positive
Sobic.001G316400Sobic.004G0710000.749.38Positive
Sobic.002G260000Sobic.004G0710000.776.928Positive
Sobic.007G115100Sobic.004G0710000.824.898Positive
Sobic.003G225300Sobic.004G0710000.814.898Positive
Sobic.005G180700Sobic.004G0710000.7720.39Positive
Sobic.003G244100Sobic.004G0710000.7312.72Positive
Sobic.003G343900Sobic.004G0710000.7419.74Positive
Sobic.004G145200Sobic.004G0710000.7519.49Positive
Sobic.004G071000Sobic.004G2203000.834.898Positive
Sobic.004G071000Sobic.001G3173000.6928.61Positive
Sobic.004G071000Sobic.001G4836000.794.898Positive
Sobic.004G071000Sobic.010G0522000.7922.44Positive
Sobic.004G071000Sobic.009G0336000.7727.14Positive
Sobic.004G071000Sobic.010G1924000.728.46Positive
Sobic.004G071000Sobic.007G0473000.843.162Positive
Sobic.004G071000Sobic.010G2096000.7713.49Positive
Sobic.004G071000Sobic.009G0272000.6715.87Positive
Sobic.004G071000Sobic.007G1799000.87.071Positive
Sobic.004G071000Sobic.008G0089000.7328.98Positive
Sobic.004G071000Sobic.006G0820000.7312.04Positive
Sobic.004G071000Sobic.004G2123000.7817.6Positive
Sobic.004G071000Sobic.001G0529000.82.449Positive
Sobic.004G071000Sobic.004G2329000.793Positive
Sobic.004G071000Sobic.010G0307000.775.656Positive
Sobic.004G273800Sobic.001G3164000.7217.54Positive
Sobic.004G273800Sobic.003G2253000.7618.33Positive
Sobic.004G273800Sobic.003G3439000.7321.63Positive
Sobic.004G273800Sobic.004G2203000.7524.49Positive
Sobic.004G273800Sobic.001G4836000.82.828Positive
Sobic.004G273800Sobic.010G2096000.7514.96Positive
Sobic.004G273800Sobic.007G1799000.841Positive
Sobic.004G273800Sobic.004G2329000.749.165Positive
Sobic.001G316400Sobic.003G2253000.748Positive
Sobic.001G316400Sobic.004G2203000.775.196Positive
Sobic.001G316400Sobic.001G4836000.7122.97Positive
Sobic.001G316400Sobic.007G0473000.7416.88Positive
Sobic.001G316400Sobic.004G2123000.7322.8Positive
Sobic.001G316400Sobic.004G2329000.6927.92Positive
Sobic.002G260000Sobic.003G2441000.712.84Positive
Sobic.002G260000Sobic.004G1452000.784.242Positive
Sobic.002G260000Sobic.001G4836000.748.124Positive
Sobic.002G260000Sobic.010G1924000.811Positive
Sobic.002G260000Sobic.009G0272000.6425.9Positive
Sobic.002G260000Sobic.004G2329000.728.366Positive
Sobic.002G260000Sobic.010G0307000.7211.48Positive
Sobic.007G115100Sobic.003G2253000.7729.34Positive
Sobic.007G115100Sobic.003G2441000.7320.97Positive
Sobic.007G115100Sobic.010G0522000.8122.58Positive
Sobic.007G115100Sobic.009G0336000.7928.58Positive
Sobic.007G115100Sobic.007G0473000.8115.71Positive
Sobic.007G115100Sobic.007G1799000.819.899Positive
Sobic.007G115100Sobic.001G0529000.7618.43Positive
Sobic.003G225300Sobic.005G1807000.841.414Positive
Sobic.003G225300Sobic.003G3439000.7810.9Positive
Sobic.003G225300Sobic.004G2203000.7915.49Positive
Sobic.003G225300Sobic.001G4836000.822.236Positive
Sobic.003G225300Sobic.007G0473000.836Positive
Sobic.003G225300Sobic.010G2096000.86.928Positive
Sobic.003G225300Sobic.004G2123000.7822.64Positive
Sobic.003G225300Sobic.010G0307000.775Positive
Sobic.005G180700Sobic.003G2441000.7223.83Positive
Sobic.005G180700Sobic.004G2203000.828.366Positive
Sobic.005G180700Sobic.007G0473000.7923.66Positive
Sobic.005G180700Sobic.010G2096000.842.449Positive
Sobic.003G244100Sobic.007G0473000.7513.26Positive
Sobic.003G244100Sobic.004G2123000.7227Positive
Sobic.003G244100Sobic.010G0307000.6917.97Positive
Sobic.003G343900Sobic.004G1452000.7326.26Positive
Sobic.003G343900Sobic.001G4836000.778.366Positive
Sobic.003G343900Sobic.007G0473000.833.741Positive
Sobic.003G343900Sobic.004G2123000.7622.24Positive
Sobic.004G145200Sobic.010G1924000.768.062Positive
Sobic.004G145200Sobic.007G0473000.816.782Positive
Sobic.004G145200Sobic.008G0089000.7329.98Positive
Sobic.004G220300Sobic.010G0522000.874Positive
Sobic.004G220300Sobic.007G0473000.847.416Positive
Sobic.004G220300Sobic.007G1799000.826.708Positive
Sobic.004G220300Sobic.004G2123000.854.898Positive
Sobic.004G220300Sobic.004G2329000.7127.92Positive
Sobic.001G317300Sobic.006G0820000.6729.98Positive
Sobic.001G317300Sobic.010G0307000.6915.87Positive
Sobic.001G483600Sobic.010G2096000.779.165Positive
Sobic.001G483600Sobic.004G2329000.754.898Positive
Sobic.001G483600Sobic.010G0307000.757.937Positive
Sobic.010G052200Sobic.009G0336000.894Positive
Sobic.010G052200Sobic.007G0473000.912.828Positive
Sobic.010G052200Sobic.004G2123000.8510Positive
Sobic.009G033600Sobic.007G0473000.911.732Positive
Sobic.009G033600Sobic.004G2123000.821.02Positive
Sobic.010G192400Sobic.009G0272000.6614.56Positive
Sobic.010G192400Sobic.010G0307000.756.708Positive
Sobic.007G047300Sobic.004G2123000.921Positive
Sobic.010G209600Sobic.010G0307000.7510.95Positive
Sobic.009G027200Sobic.001G0529000.733.741Positive
Sobic.007G179900Sobic.004G2329000.7318.33Positive
Sobic.008G008900Sobic.006G0820000.748.66Positive
Sobic.006G082000Sobic.001G0529000.746.633Positive
Sobic.007G207400Sobic.004G071000−0.6523.36Negative
Sobic.004G317300Sobic.004G071000−0.717Negative
Sobic.002G031000Sobic.004G071000−0.749.591Negative
Sobic.003G426700Sobic.004G071000−0.7421.09Negative
Sobic.004G071000Sobic.005G005500−0.7324.49Negative
Sobic.004G071000Sobic.009G254500−0.7310.58Negative
Sobic.004G071000Sobic.001G035800−0.7714.49Negative
Sobic.004G071000Sobic.009G191900−0.7229.25Negative
Sobic.004G071000Sobic.004G041500−0.649.433Negative
Sobic.004G071000Sobic.007G142200−0.7114.69Negative
Sobic.004G071000Sobic.010G128200−0.7518.41Negative
Sobic.004G071000Sobic.010G202700−0.79.165Negative
Sobic.004G071000Sobic.010G190700−0.7910.29Negative
Sobic.004G071000Sobic.009G127500−0.714.14Negative
Sobic.004G071000Sobic.010G232100Predicted PPI
Sobic.004G071000Sobic.003G403700Predicted PPI
Sobic.004G071000Sobic.004G049900Predicted PPI
Sobic.004G071000Sobic.004G050000Predicted PPI
Sobic.004G071000Sobic.004G050100Predicted PPI
Sobic.004G071000Sobic.010G239500Predicted PPI
Table 2.

Details of the CAD gene Sobic.004G071000’s network

Protein AProtein BPCCMRRelationship
Sobic.004G273800Sobic.004G0710000.7420.49Positive
Sobic.001G316400Sobic.004G0710000.749.38Positive
Sobic.002G260000Sobic.004G0710000.776.928Positive
Sobic.007G115100Sobic.004G0710000.824.898Positive
Sobic.003G225300Sobic.004G0710000.814.898Positive
Sobic.005G180700Sobic.004G0710000.7720.39Positive
Sobic.003G244100Sobic.004G0710000.7312.72Positive
Sobic.003G343900Sobic.004G0710000.7419.74Positive
Sobic.004G145200Sobic.004G0710000.7519.49Positive
Sobic.004G071000Sobic.004G2203000.834.898Positive
Sobic.004G071000Sobic.001G3173000.6928.61Positive
Sobic.004G071000Sobic.001G4836000.794.898Positive
Sobic.004G071000Sobic.010G0522000.7922.44Positive
Sobic.004G071000Sobic.009G0336000.7727.14Positive
Sobic.004G071000Sobic.010G1924000.728.46Positive
Sobic.004G071000Sobic.007G0473000.843.162Positive
Sobic.004G071000Sobic.010G2096000.7713.49Positive
Sobic.004G071000Sobic.009G0272000.6715.87Positive
Sobic.004G071000Sobic.007G1799000.87.071Positive
Sobic.004G071000Sobic.008G0089000.7328.98Positive
Sobic.004G071000Sobic.006G0820000.7312.04Positive
Sobic.004G071000Sobic.004G2123000.7817.6Positive
Sobic.004G071000Sobic.001G0529000.82.449Positive
Sobic.004G071000Sobic.004G2329000.793Positive
Sobic.004G071000Sobic.010G0307000.775.656Positive
Sobic.004G273800Sobic.001G3164000.7217.54Positive
Sobic.004G273800Sobic.003G2253000.7618.33Positive
Sobic.004G273800Sobic.003G3439000.7321.63Positive
Sobic.004G273800Sobic.004G2203000.7524.49Positive
Sobic.004G273800Sobic.001G4836000.82.828Positive
Sobic.004G273800Sobic.010G2096000.7514.96Positive
Sobic.004G273800Sobic.007G1799000.841Positive
Sobic.004G273800Sobic.004G2329000.749.165Positive
Sobic.001G316400Sobic.003G2253000.748Positive
Sobic.001G316400Sobic.004G2203000.775.196Positive
Sobic.001G316400Sobic.001G4836000.7122.97Positive
Sobic.001G316400Sobic.007G0473000.7416.88Positive
Sobic.001G316400Sobic.004G2123000.7322.8Positive
Sobic.001G316400Sobic.004G2329000.6927.92Positive
Sobic.002G260000Sobic.003G2441000.712.84Positive
Sobic.002G260000Sobic.004G1452000.784.242Positive
Sobic.002G260000Sobic.001G4836000.748.124Positive
Sobic.002G260000Sobic.010G1924000.811Positive
Sobic.002G260000Sobic.009G0272000.6425.9Positive
Sobic.002G260000Sobic.004G2329000.728.366Positive
Sobic.002G260000Sobic.010G0307000.7211.48Positive
Sobic.007G115100Sobic.003G2253000.7729.34Positive
Sobic.007G115100Sobic.003G2441000.7320.97Positive
Sobic.007G115100Sobic.010G0522000.8122.58Positive
Sobic.007G115100Sobic.009G0336000.7928.58Positive
Sobic.007G115100Sobic.007G0473000.8115.71Positive
Sobic.007G115100Sobic.007G1799000.819.899Positive
Sobic.007G115100Sobic.001G0529000.7618.43Positive
Sobic.003G225300Sobic.005G1807000.841.414Positive
Sobic.003G225300Sobic.003G3439000.7810.9Positive
Sobic.003G225300Sobic.004G2203000.7915.49Positive
Sobic.003G225300Sobic.001G4836000.822.236Positive
Sobic.003G225300Sobic.007G0473000.836Positive
Sobic.003G225300Sobic.010G2096000.86.928Positive
Sobic.003G225300Sobic.004G2123000.7822.64Positive
Sobic.003G225300Sobic.010G0307000.775Positive
Sobic.005G180700Sobic.003G2441000.7223.83Positive
Sobic.005G180700Sobic.004G2203000.828.366Positive
Sobic.005G180700Sobic.007G0473000.7923.66Positive
Sobic.005G180700Sobic.010G2096000.842.449Positive
Sobic.003G244100Sobic.007G0473000.7513.26Positive
Sobic.003G244100Sobic.004G2123000.7227Positive
Sobic.003G244100Sobic.010G0307000.6917.97Positive
Sobic.003G343900Sobic.004G1452000.7326.26Positive
Sobic.003G343900Sobic.001G4836000.778.366Positive
Sobic.003G343900Sobic.007G0473000.833.741Positive
Sobic.003G343900Sobic.004G2123000.7622.24Positive
Sobic.004G145200Sobic.010G1924000.768.062Positive
Sobic.004G145200Sobic.007G0473000.816.782Positive
Sobic.004G145200Sobic.008G0089000.7329.98Positive
Sobic.004G220300Sobic.010G0522000.874Positive
Sobic.004G220300Sobic.007G0473000.847.416Positive
Sobic.004G220300Sobic.007G1799000.826.708Positive
Sobic.004G220300Sobic.004G2123000.854.898Positive
Sobic.004G220300Sobic.004G2329000.7127.92Positive
Sobic.001G317300Sobic.006G0820000.6729.98Positive
Sobic.001G317300Sobic.010G0307000.6915.87Positive
Sobic.001G483600Sobic.010G2096000.779.165Positive
Sobic.001G483600Sobic.004G2329000.754.898Positive
Sobic.001G483600Sobic.010G0307000.757.937Positive
Sobic.010G052200Sobic.009G0336000.894Positive
Sobic.010G052200Sobic.007G0473000.912.828Positive
Sobic.010G052200Sobic.004G2123000.8510Positive
Sobic.009G033600Sobic.007G0473000.911.732Positive
Sobic.009G033600Sobic.004G2123000.821.02Positive
Sobic.010G192400Sobic.009G0272000.6614.56Positive
Sobic.010G192400Sobic.010G0307000.756.708Positive
Sobic.007G047300Sobic.004G2123000.921Positive
Sobic.010G209600Sobic.010G0307000.7510.95Positive
Sobic.009G027200Sobic.001G0529000.733.741Positive
Sobic.007G179900Sobic.004G2329000.7318.33Positive
Sobic.008G008900Sobic.006G0820000.748.66Positive
Sobic.006G082000Sobic.001G0529000.746.633Positive
Sobic.007G207400Sobic.004G071000−0.6523.36Negative
Sobic.004G317300Sobic.004G071000−0.717Negative
Sobic.002G031000Sobic.004G071000−0.749.591Negative
Sobic.003G426700Sobic.004G071000−0.7421.09Negative
Sobic.004G071000Sobic.005G005500−0.7324.49Negative
Sobic.004G071000Sobic.009G254500−0.7310.58Negative
Sobic.004G071000Sobic.001G035800−0.7714.49Negative
Sobic.004G071000Sobic.009G191900−0.7229.25Negative
Sobic.004G071000Sobic.004G041500−0.649.433Negative
Sobic.004G071000Sobic.007G142200−0.7114.69Negative
Sobic.004G071000Sobic.010G128200−0.7518.41Negative
Sobic.004G071000Sobic.010G202700−0.79.165Negative
Sobic.004G071000Sobic.010G190700−0.7910.29Negative
Sobic.004G071000Sobic.009G127500−0.714.14Negative
Sobic.004G071000Sobic.010G232100Predicted PPI
Sobic.004G071000Sobic.003G403700Predicted PPI
Sobic.004G071000Sobic.004G049900Predicted PPI
Sobic.004G071000Sobic.004G050000Predicted PPI
Sobic.004G071000Sobic.004G050100Predicted PPI
Sobic.004G071000Sobic.010G239500Predicted PPI
Protein AProtein BPCCMRRelationship
Sobic.004G273800Sobic.004G0710000.7420.49Positive
Sobic.001G316400Sobic.004G0710000.749.38Positive
Sobic.002G260000Sobic.004G0710000.776.928Positive
Sobic.007G115100Sobic.004G0710000.824.898Positive
Sobic.003G225300Sobic.004G0710000.814.898Positive
Sobic.005G180700Sobic.004G0710000.7720.39Positive
Sobic.003G244100Sobic.004G0710000.7312.72Positive
Sobic.003G343900Sobic.004G0710000.7419.74Positive
Sobic.004G145200Sobic.004G0710000.7519.49Positive
Sobic.004G071000Sobic.004G2203000.834.898Positive
Sobic.004G071000Sobic.001G3173000.6928.61Positive
Sobic.004G071000Sobic.001G4836000.794.898Positive
Sobic.004G071000Sobic.010G0522000.7922.44Positive
Sobic.004G071000Sobic.009G0336000.7727.14Positive
Sobic.004G071000Sobic.010G1924000.728.46Positive
Sobic.004G071000Sobic.007G0473000.843.162Positive
Sobic.004G071000Sobic.010G2096000.7713.49Positive
Sobic.004G071000Sobic.009G0272000.6715.87Positive
Sobic.004G071000Sobic.007G1799000.87.071Positive
Sobic.004G071000Sobic.008G0089000.7328.98Positive
Sobic.004G071000Sobic.006G0820000.7312.04Positive
Sobic.004G071000Sobic.004G2123000.7817.6Positive
Sobic.004G071000Sobic.001G0529000.82.449Positive
Sobic.004G071000Sobic.004G2329000.793Positive
Sobic.004G071000Sobic.010G0307000.775.656Positive
Sobic.004G273800Sobic.001G3164000.7217.54Positive
Sobic.004G273800Sobic.003G2253000.7618.33Positive
Sobic.004G273800Sobic.003G3439000.7321.63Positive
Sobic.004G273800Sobic.004G2203000.7524.49Positive
Sobic.004G273800Sobic.001G4836000.82.828Positive
Sobic.004G273800Sobic.010G2096000.7514.96Positive
Sobic.004G273800Sobic.007G1799000.841Positive
Sobic.004G273800Sobic.004G2329000.749.165Positive
Sobic.001G316400Sobic.003G2253000.748Positive
Sobic.001G316400Sobic.004G2203000.775.196Positive
Sobic.001G316400Sobic.001G4836000.7122.97Positive
Sobic.001G316400Sobic.007G0473000.7416.88Positive
Sobic.001G316400Sobic.004G2123000.7322.8Positive
Sobic.001G316400Sobic.004G2329000.6927.92Positive
Sobic.002G260000Sobic.003G2441000.712.84Positive
Sobic.002G260000Sobic.004G1452000.784.242Positive
Sobic.002G260000Sobic.001G4836000.748.124Positive
Sobic.002G260000Sobic.010G1924000.811Positive
Sobic.002G260000Sobic.009G0272000.6425.9Positive
Sobic.002G260000Sobic.004G2329000.728.366Positive
Sobic.002G260000Sobic.010G0307000.7211.48Positive
Sobic.007G115100Sobic.003G2253000.7729.34Positive
Sobic.007G115100Sobic.003G2441000.7320.97Positive
Sobic.007G115100Sobic.010G0522000.8122.58Positive
Sobic.007G115100Sobic.009G0336000.7928.58Positive
Sobic.007G115100Sobic.007G0473000.8115.71Positive
Sobic.007G115100Sobic.007G1799000.819.899Positive
Sobic.007G115100Sobic.001G0529000.7618.43Positive
Sobic.003G225300Sobic.005G1807000.841.414Positive
Sobic.003G225300Sobic.003G3439000.7810.9Positive
Sobic.003G225300Sobic.004G2203000.7915.49Positive
Sobic.003G225300Sobic.001G4836000.822.236Positive
Sobic.003G225300Sobic.007G0473000.836Positive
Sobic.003G225300Sobic.010G2096000.86.928Positive
Sobic.003G225300Sobic.004G2123000.7822.64Positive
Sobic.003G225300Sobic.010G0307000.775Positive
Sobic.005G180700Sobic.003G2441000.7223.83Positive
Sobic.005G180700Sobic.004G2203000.828.366Positive
Sobic.005G180700Sobic.007G0473000.7923.66Positive
Sobic.005G180700Sobic.010G2096000.842.449Positive
Sobic.003G244100Sobic.007G0473000.7513.26Positive
Sobic.003G244100Sobic.004G2123000.7227Positive
Sobic.003G244100Sobic.010G0307000.6917.97Positive
Sobic.003G343900Sobic.004G1452000.7326.26Positive
Sobic.003G343900Sobic.001G4836000.778.366Positive
Sobic.003G343900Sobic.007G0473000.833.741Positive
Sobic.003G343900Sobic.004G2123000.7622.24Positive
Sobic.004G145200Sobic.010G1924000.768.062Positive
Sobic.004G145200Sobic.007G0473000.816.782Positive
Sobic.004G145200Sobic.008G0089000.7329.98Positive
Sobic.004G220300Sobic.010G0522000.874Positive
Sobic.004G220300Sobic.007G0473000.847.416Positive
Sobic.004G220300Sobic.007G1799000.826.708Positive
Sobic.004G220300Sobic.004G2123000.854.898Positive
Sobic.004G220300Sobic.004G2329000.7127.92Positive
Sobic.001G317300Sobic.006G0820000.6729.98Positive
Sobic.001G317300Sobic.010G0307000.6915.87Positive
Sobic.001G483600Sobic.010G2096000.779.165Positive
Sobic.001G483600Sobic.004G2329000.754.898Positive
Sobic.001G483600Sobic.010G0307000.757.937Positive
Sobic.010G052200Sobic.009G0336000.894Positive
Sobic.010G052200Sobic.007G0473000.912.828Positive
Sobic.010G052200Sobic.004G2123000.8510Positive
Sobic.009G033600Sobic.007G0473000.911.732Positive
Sobic.009G033600Sobic.004G2123000.821.02Positive
Sobic.010G192400Sobic.009G0272000.6614.56Positive
Sobic.010G192400Sobic.010G0307000.756.708Positive
Sobic.007G047300Sobic.004G2123000.921Positive
Sobic.010G209600Sobic.010G0307000.7510.95Positive
Sobic.009G027200Sobic.001G0529000.733.741Positive
Sobic.007G179900Sobic.004G2329000.7318.33Positive
Sobic.008G008900Sobic.006G0820000.748.66Positive
Sobic.006G082000Sobic.001G0529000.746.633Positive
Sobic.007G207400Sobic.004G071000−0.6523.36Negative
Sobic.004G317300Sobic.004G071000−0.717Negative
Sobic.002G031000Sobic.004G071000−0.749.591Negative
Sobic.003G426700Sobic.004G071000−0.7421.09Negative
Sobic.004G071000Sobic.005G005500−0.7324.49Negative
Sobic.004G071000Sobic.009G254500−0.7310.58Negative
Sobic.004G071000Sobic.001G035800−0.7714.49Negative
Sobic.004G071000Sobic.009G191900−0.7229.25Negative
Sobic.004G071000Sobic.004G041500−0.649.433Negative
Sobic.004G071000Sobic.007G142200−0.7114.69Negative
Sobic.004G071000Sobic.010G128200−0.7518.41Negative
Sobic.004G071000Sobic.010G202700−0.79.165Negative
Sobic.004G071000Sobic.010G190700−0.7910.29Negative
Sobic.004G071000Sobic.009G127500−0.714.14Negative
Sobic.004G071000Sobic.010G232100Predicted PPI
Sobic.004G071000Sobic.003G403700Predicted PPI
Sobic.004G071000Sobic.004G049900Predicted PPI
Sobic.004G071000Sobic.004G050000Predicted PPI
Sobic.004G071000Sobic.004G050100Predicted PPI
Sobic.004G071000Sobic.010G239500Predicted PPI

Furthermore, we executed a motif enrichment analysis using the 46 genes. Among the significant enrichment (P-value < 0.05) motifs, some MYB- and bHLH-related binding sites were found. Notably, most of the network genes, including the query gene Sobic.004G071000, had MYB TF binding sites.

Then, a CAD co-expression network gene (GSE49879) expression matrix was used to investigate expression patterns. The data of PI455230, which contains most of the samples and tissues of the six genotypes, was used. A heatmap was produced by Cluster 3.0 software (67), with colors from blue to orange corresponding to expression profiles from low to high (Figure 3B). In the heatmap, we can see that Sobic.004G071000 was highly expressed in rind and pind, instead of shoot and leaf. In addition, a similarity expression trend was seen in Sobic.004G232900, another lignin pathway-related gene. However, Sobic.002G031000, encoding a bHLH transcript factor, negatively regulated lignin biosynthesis, which had the opposite expression profile as the genes mentioned above.

Finally, lignin is mainly deposited on the stem, especially in more mature plants (68). Thus, we set an expression pattern threshold (high in the 70 DAP pith and rind, even in the 60 DAP pith and rind, and low in the last samples) to find more lignin pathway-related genes with the pattern set tool. As a result, we found 347 genes with a PCC threshold over 0.7, and the genes Sobic.009G186600 (PCC = 0.76), Sobic.010G245500 (PCC = 0.73) and Sobic.009G186600 (PCC = 0.73), participate in the phenylalanine metabolism pathway (Supplementary Table S4).

Overall, our analyses results corroborated those of other published paper. For example, Sobic.004G273800 is a MYB TF, and is up-regulated with the CAD gene (62), while Sobic.002g031000 is a bHLH TF that is down-regulated with the CAD gene (65). In the network, genes also have been reported, such as Sobic.004G220300 (phenylalanine ammonia lyase), Sobic.010G052200 (Caffeoyl-CoA 3-O-methyltransferase, Sobic.007G047300 (Caffeic acid O-methyltransferase) and Sobic.004G212300 (Hydroxyl-cinnamoyl CoA:shikimate/quinate-Hydroxyl-cinnamoyltransferase) (16).

In a summary, we integrated analysis and visualization methods into SorghumFDB to predict the lignin pathway-related regulators and found genes that had been identified by published papers, as well as some new genes, that may play roles in metabolic processes. However, the analysis flow to integrate the network analysis, GSEA tool, pattern set and motif analysis, may have been a good way to explore the gene functions and regulatory relationship. In the development of sorghum research, newer annotations and tools may be added to our database to comprehensively analyse other biological problems, such as C4 photosynthesis and drought stress.

Conclusion

At present, SorghumFDB is a platform for genome functional annotations and multi-dimensional network analyses. It encompassed most information, such as various annotations of whole genome assemblies, miRNA sequences and target genes, common gene families, network constructions using transcriptome data, PPI data and miRNA-target pairs, as well as multiple gene function annotation elements. Visualization tools (Gbrowse, Cytoscape and open-flash-chart) and four analysis-based tools, BLAST, GSEA, motif significance analysis and pattern set, were provided to determine the functional prediction. The SorghumFDB is online at http://structuralbiology.cau.edu.cn/sorghum/index.html. We hope that this will improve the accuracy and robustness of sorghum functional genomics analyses, and further help in understanding the gene regulatory networks involved in effective crop improvement.

Funding

This work was supported by grants from the Ministry of Science and Technology of China (31371291 and 2012CB215301).

Conflict of interest. None declared.

References

1

Shen
H.
Mazarei
M.
Hisano
H
. et al. . (
2013
)
A genomics approach to deciphering lignin biosynthesis in switchgrass
.
Plant Cell
,
25
,
4342
4361
.

2

Lockhart
J.
(
2015
)
Altering lignin composition to improve biofuel production
.
Plant Cell
,
27
,
2082.

3

Vandenbrink
F.A.F.a.J.P.
(
2012
)
Bioenergy grass feedstock: current options and prospects for trait improvement using emerging genetic, genomic, and systems biology toolkits
.
Biotechnol. Biofuels
.,
5
,
80
.

4

Kim
J.S.
Lee
Y.Y.
Kim
T.H.
(
2016
)
A review on alkaline pretreatment technology for bioconversion of lignocellulosic biomass
.
Bioresource Technol
.,
199
,
42
48
.

5

Peplow
M.
(
2014
)
Cellulosic ethanol fights for life
.
Nature
,
507
,
152
153
.

6

Li
X.
Weng
J.K.
Chapple
C.
(
2008
)
Improvement of biomass through lignin modification
.
Plant J. Cell Mol. Biol
.,
54
,
569
581
.

7

Park
S.H.
Ong
R.G.
Sticklen
M.
(
2015
)
Strategies for the production of cell wall-deconstructing enzymes in lignocellulosic biomass and their utilization for biofuel production
.
Plant Biotechnol. J
.,
14
,
1329
44
.

8

Paterson
A.H.
Bowers
J.E.
Bruggmann
R
. et al. . (
2009
)
The Sorghum bicolor genome and the diversification of grasses
.
Nature
,
457
,
551
556
.

9

Goodstein
D.M.
Shu
S.
Howson
R
. et al. . (
2012
)
Phytozome: a comparative platform for green plant genomics
.
Nucleic Acids Res
.,
40
,
D1178
D1186
.

10

Youens-Clark
K.
Buckler
E.
Casstevens
T
. et al. . (
2011
)
Gramene database in 2010: updates and extensions
.
Nucleic Acids Res
.,
39
,
D1085
D1094
.

11

Proost
S.
Van Bel
M.
Vaneechoutte
D
. et al. . (
2015
)
PLAZA 3.0: an access point for plant comparative genomics
.
Nucleic Acids Res
.,
43
,
D974
D981
.

12

Nussbaumer
T.
Martis
M.M.
Roessner
S.K
. et al. . (
2013
)
MIPS PlantsDB: a database framework for comparative plant genome research
.
Nucleic Acids Res
.,
41
,
D1144
D1151
.

13

Makita
Y.
Shimada
S.
Kawashima
M
. et al. . (
2015
)
MOROKOSHI: transcriptome database in Sorghum bicolor
.
Plant Cell Physiol
.,
56
,
e6.

14

Olson
A.
Klein
R.R.
Dugas
D.V
. et al. . (
2014
)
Expanding and vetting gene annotations through transcriptome and methylome sequencing
.
Plant Genome
,
7
,
0.

15

Gelli
M.
Duo
Y.
Konda
A.R
. et al. . (
2014
)
Identification of differentially expressed genes between sorghum genotypes with contrasting nitrogen stress tolerance by genome-wide transcriptional profiling
.
BMC Genomics
,
15
,
179.

16

Shakoor
N.
Nair
R.
Crasta
O.
et al. . (
2014
)
A Sorghum bicolor expression atlas reveals dynamic genotype-specific expression profiles for vegetative tissues of grain, sweet and bioenergy sorghums
.
BMC Bioinformatics
.,
14
,
35
.

17

Tatusov
R.L.
Fedorova
N.D.
Jackson
J.D
. et al. . (
2003
)
The COG database: an updated version includes eukaryotes
.
BMC Bioinformatics
,
4
,
41.

18

Ashburner
M.
Ball
C.A.
Blake
J.A
. et al. . (
2000
)
Gene ontology: tool for the unification of biology. The Gene Ontology Consortium
.
Nat. Genet
.,
25
,
25
29
.

19

Van Peer
G.
Lefever
S.
Anckaert
J
. et al. . (
2014
)
miRBase Tracker: keeping track of microRNA annotation changes
.
Database J. Biol. Datab. Curat
., 2014, bau080.

20

Zhang
Z.
Yu
J.
Li
D
. et al. . (
2010
)
PMRD: plant microRNA database
.
Nucleic Acids Res
.,
38
,
D806
D813
.

21

Lescot
M.
Dehais
P.
Thijs
G
. et al. . (
2002
)
PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences
.
Nucleic Acids Res
.,
30
,
325
327
.

22

Higo
K.
Ugawa
Y.
Iwamoto
M
. et al. . (
1999
)
Plant cis-acting regulatory DNA elements (PLACE) database: 1999
.
Nucleic Acids Res
.,
27
,
297
300
.

23

Hehl
R.
Bulow
L.
(
2014
)
AthaMap web tools for the analysis of transcriptional and posttranscriptional regulation of gene expression in Arabidopsis thaliana
.
Methods Mol. Biol
.,
1158
,
139
156
.

24

Jin
J.
Zhang
H.
Kong
L
. et al. . (
2014
)
PlantTFDB 3.0: a portal for the functional and evolutionary study of plant transcription factors
.
Nucleic Acids Res
.,
42
,
D1182
D1187
.

25

Riano-Pachon
D.M.
Ruzicic
S.
Dreyer
I
. et al. . (
2007
)
PlnTFDB: an integrative plant transcription factor database
.
BMC Bioinformatics
,
8
,
42.

26

Nelson
D.
Werck-Reichhart
D.
(
2011
)
A P450-centric view of plant evolution
.
Plant J. Cell Mol. Biol
.,
66
,
194
211
.

27

Gao
T.
Liu
Z.
Wang
Y
. et al. . (
2013
)
UUCD: a family-based database of ubiquitin and ubiquitin-like conjugation
.
Nucleic Acids Res
.,
41
,
D445
D451
.

28

Mace
E.
Tai
S.
Innes
D
. et al. . (
2014
)
The plasticity of NBS resistance genes in sorghum is driven by multiple evolutionary processes
.
BMC Plant Biol
.,
14
,
253.

29

Mi
H.
Muruganujan
A.
Thomas
P.D.
(
2013
)
PANTHER in 2013: modeling the evolution of gene function, and other gene attributes, in the context of phylogenetic trees
.
Nucleic Acids Res
.,
41
,
D377
D386
.

30

Hermjakob
H.
Montecchi-Palazzi
L.
Lewington
C
. et al. . (
2004
)
IntAct: an open source molecular interaction database
.
Nucleic Acids Res
.,
32
,
D452
D455
.

31

Stark
C.
Breitkreutz
B.J.
Reguly
T
. et al. . (
2006
)
BioGRID: a general repository for interaction datasets
.
Nucleic Acids Res
.,
34
,
D535
D539
.

32

Orchard
S.
Kerrien
S.
Abbani
S
. et al. . (
2012
)
Protein interaction data curation: the International Molecular Exchange (IMEx) consortium
.
Nat. Methods
,
9
,
345
350
.

33

Poole
R.L.
(
2007
)
The TAIR database
.
Methods Mol. Biol
.,
406
,
179
212
.

34

Provart
N.J.
Gil
P.
Chen
W
. et al. . (
2003
)
Gene expression phenotypes of Arabidopsis associated with sensitivity to low temperatures
.
Plant Physiol
.,
132
,
893
906
.

35

Arabidopsis Interactome Mapping
C.
(
2011
)
Evidence for network evolution in an Arabidopsis interactome map
.
Science
,
333
,
601
607
.

36

Lumba
S.
Toh
S.
Handfield
L.F
. et al. . (
2014
)
A mesoscale abscisic acid hormone interactome reveals a dynamic signaling landscape in Arabidopsis
.
Dev. Cell
,
29
,
360
372
.

37

Finn
R.D.
Bateman
A.
Clements
J
. et al. . (
2014
)
Pfam: the protein families database
.
Nucleic Acids Res
.,
42
,
D222
D230
.

38

Chae
L.
Kim
T.
Nilo-Poyanco
R
. et al. . (
2014
)
Genomic signatures of specialized metabolism in plants
.
Science
,
344
,
510
513
.

39

You
Q.
Zhang
L.
Yi
X
. et al. . (
2015
)
SIFGD: Setaria italica functional genomics database
.
Mol. Plant
,
8
,
967
970
.

40

Zhang
L.
Guo
J.
You
Q
. et al. . (
2015
)
GraP: platform for functional genomics analysis of Gossypium raimondii
.
Datab. J. Biol. Datab. Curat
.,
2015
,
bav047.

41

Sage
R.F.
Zhu
X.G.
(
2011
)
Exploiting the engine of C(4) photosynthesis
.
J. Exp. Bot
.,
62
,
2989
3000
.

42

Lassen
L.M.
Nielsen
A.Z.
Ziersen
B
. et al. . (
2014
)
Redirecting photosynthetic electron flow into light-driven synthesis of alternative products including high-value bioactive natural compounds
.
ACS Synth. Biol
.,
3
,
1
12
.

43

Yu
X.
Liu
Y.
Wang
S
. et al. . (
2015
)
CarNAC4, a NAC-type chickpea transcription factor conferring enhanced drought and salt stress tolerances in Arabidopsis
.
Plant Cell Rep
.,
35
,
613
27
.

44

Li
Z.
Wang
W.
Li
G
. et al. . (
2015
)
MAPK-mediated regulation of growth and essential oil composition in a salt-tolerant peppermint (Mentha piperita L.) under NaCl stress
.
Protoplasma
., DOI: 10.1007/s00709-015-0915-1.

45

Xu
J.
Xing
S.
Cui
H
. et al. . (
2015
)
Genome-wide identification and characterization of the apple (Malus domestica) HECT ubiquitin-protein ligase family and expression analysis of their responsiveness to abiotic stresses
.
Mol. Genet. Genomics MGG
.,
291
,
635
46
.

46

Gururani
M.A.
Venkatesh
J.
Upadhyaya
C.P
. et al. . (
2012
)
Plant disease resistance genes: Current status and future directions
.
Physiol. Mol. Plant Pathol
.,
78
,
51
65
.

47

Sonnhammer
E.L.
Ostlund
G.
(
2015
)
InParanoid 8: orthology analysis between 273 proteomes, mostly eukaryotic
.
Nucleic Acids Res
.,
43
,
D234
D239
.

48

Lehti-Shiu
M.D.
Shiu
S.H.
(
2012
)
Diversity, classification and function of the plant protein kinase superfamily
.
Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci
.,
367
,
2619
2639
.

49

Wu
T.D.
Watanabe
C.K.
(
2005
)
GMAP: a genomic mapping and alignment program for mRNA and EST sequences
.
Bioinformatics
,
21
,
1859
1875
.

50

Dai
X.
Zhao
P.X.
(
2011
)
psRNATarget: a plant small RNA target analysis server
.
Nucleic Acids Res
.,
39
,
W155
W159
.

51

Yu
J.
Zhang
Z.
Wei
J
. et al. . (
2014
)
SFGD: a comprehensive platform for mining functional information from soybean transcriptome data and its use in identifying acyl-lipid metabolism pathways
.
BMC Genomics
,
15
,
271.

52

Obayashi
T.
Kinoshita
K.
(
2009
)
Rank of correlation coefficient as a comparable measure for biological significance of gene coexpression
.
DNA Res. Int. J. Rapid Publ. Rep. Genes Genomes
,
16
,
249
260
.

53

Obayashi
T.
Hayashi
S.
Saeki
M
. et al. . (
2009
)
ATTED-II provides coexpressed gene networks for Arabidopsis
.
Nucleic Acids Res
.,
37
,
D987
D991
.

54

Li
C.
Zhang
B.
(
2016
)
MicroRNAs in control of plant development
.
J. Cell. Physiol
.,
231
,
303
313
.

55

Niklas
N.
Hafenscher
J.
Barna
A
. et al. . (
2015
)
cFinder: definition and quantification of multiple haplotypes in a mixed sample
.
BMC Res. Notes
,
8
,
422.

56

Nemhauser
J.L.
Mockler
T.C.
Chory
J.
(
2004
)
Interdependency of brassinosteroid and auxin signaling in Arabidopsis
.
PLoS Biol
.,
2
,
E258.

57

Yi
X.
Du
Z.
Su
Z.
(
2013
)
PlantGSEA: a gene set enrichment analysis toolkit for plant community
.
Nucleic Acids Res
.,
41
,
W98
103
.

58

Su
G.
Morris
J.H.
Demchak
B
. et al. . (
2014
)
Biological network exploration with cytoscape 3
.
Curr. Protoc. Bioinf
.,
47
,
8.13.1
24
.

59

Trapnell
C.
Roberts
A.
Goff
L
. et al. . (
2012
)
Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks
.
Nat. Protoc
.,
7
,
562
578
.

60

Ghosh
S.
Chan
C.K.
(
2016
)
Analysis of RNA-Seq data using TopHat and Cufflinks
.
Methods Mol. Biol
.,
1374
,
339
361
.

61

Sattler
S.E.
Saathoff
A.J.
Haas
E.J
. et al. . (
2009
)
A nonsense mutation in a cinnamyl alcohol dehydrogenase gene is responsible for the Sorghum brown midrib6 phenotype
.
Plant Physiol
.,
150
,
584
595
.

62

Li
J.
Wang
L.
Zhan
Q
. et al. . (
2013
)
Sorghum bmr6 mutant analysis demonstrates that a shared MYB1 transcription factor binding site in the promoter links the expression of genes in related pathways
.
Funct. Integr. Genomics
,
13
,
445
453
.

63

Mullet
J.
Morishige
D.
McCormick
R
. et al. . (
2014
)
Energy sorghum – a genetic model for the design of C4 grass bioenergy crops
.
J. Exp. Bot
.,
65
,
3479
3489
.

64

Dixon
R.A.
Chen
F.
Guo
D.
, et al. . (
2001
)
The biosynthesis of monolignols: a “metabolic grid”, or independent pathways to guaiacyl and syringyl units?
.
Phytochemistry
,
57
,
1069
1084
.

65

Yan
L.
Xu
C.
Kang
Y
. et al. . (
2013
)
The heterologous expression in Arabidopsis thaliana of sorghum transcription factor SbbHLH1 downregulates lignin synthesis
.
J. Exp. Bot
.,
64
,
3021
3032
.

66

Han
J.
Clement
J.M.
Li
J
. et al. . (
2010
)
The cytochrome P450 CYP86A22 is a fatty acyl-CoA omega-hydroxylase essential for Estolide synthesis in the stigma of Petunia hybrida
.
J. Biol. Chem
.,
285
,
3986
3996
.

67

de Hoon
M.J.
Imoto
S.
Nolan
J
. et al. . (
2004
)
Open source clustering software
.
Bioinformatics
,
20
,
1453
1454
.

68

Li
C.
Wang
X.
Ran
L
. et al. . (
2015
)
PtoMYB92 is a transcriptional activator of the lignin biosynthetic pathway during secondary cell wall formation in Populus tomentosa
.
Plant Cell Physiol
.,
56
,
2436
2446
.

Author notes

Citation details: Tian,T., You,Q., Zhang,L. et al. SorghumFDB: sorghum functional genomics database with multidimensional network analysis. Database, Vol. 2016: article ID baw099; doi: 10.1093/database/baw099

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Supplementary data