The dynamic impact approach (DIA) represents an alternative to overrepresentation analysis (ORA) for functional analysis of time-course experiments or those involving multiple treatments. The DIA can be used to estimate the biological impact of the differentially expressed genes (DEGs) associated with particular biological functions, for example, as represented by the Kyoto encyclopedia of genes and genomes (KEGG) annotations. However, the DIA does not take into account the correlated dependence structure of the KEGG pathway hierarchy. We have developed herein a path analysis model (KEGG-PATH) to subdivide the total effect of each KEGG pathway into the direct effect and indirect effect by taking into account not only each KEGG pathway itself, but also the correlation with its related pathways. In addition, this work also attempts to preliminarily estimate the impact direction of each KEGG pathway by a gradient analysis method from principal component analysis (PCA). As a result, the advantage of the KEGG-PATH model is demonstrated through the functional analysis of the bovine mammary transcriptome during lactation.

BackgroundThe knowledge base-driven pathway analysis is becoming the first choice for many investigators, in that it not only can reduce the complexity of functional analysis by grouping thousands of genes into just several hundred pathways, but also can increase the explanatory power for the experiment by identifying active pathways in different conditions. However, current approaches are designed to analyze a biological system assuming that each pathway is independent of the other pathways.ResultsA decision analysis model is developed in this article that accounts for dependence among pathways in time-course experiments and multiple treatments experiments. This model introduces a decision coefficient—a designed index, to identify the most relevant pathways in a given experiment by taking into account not only the direct determination factor of each Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway itself, but also the indirect determination factors from its related pathways. Meanwhile, the direct and indirect determination factors of each pathway are employed to demonstrate the regulation mechanisms among KEGG pathways, and the sign of decision coefficient can be used to preliminarily estimate the impact direction of each KEGG pathway. The simulation study of decision analysis demonstrated the application of decision analysis model for KEGG pathway analysis.ConclusionsA microarray dataset from bovine mammary tissue over entire lactation cycle was used to further illustrate our strategy. The results showed that the decision analysis model can provide the promising and more biologically meaningful results. Therefore, the decision analysis model is an initial attempt of optimizing pathway analysis methodology.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-016-1285-1) contains supplementary material, which is available to authorized users.

Cinnamoyl-CoA reductase (CCR) is the first enzyme in the monolignol-specific branch of the lignin biosynthetic pathway. In this research, three sorghum CCR genes including SbCCR1, SbCCR2-1 and SbCCR2-2 were cloned and characterized. Analyses of the structure and phylogeny of the three CCR genes showed evolutionary conservation of the functional domains and divergence of function. Transient expression assays in Nicotiana benthamiana leaves demonstrated that the three CCR proteins were localized in the cytoplasm. The expression analysis showed that the three CCR genes were induced by drought. But in 48 h, the expression levels of SbCCR1 and SbCCR2-2 did not differ between CK and the drought treatment; while the expression level of SbCCR2-1 in the drought treatment was higher than in CK. The expression of the SbCCR1 and SbCCR2-1 genes was not induced by sorghum aphid [Melanaphis sacchari (Zehntner)] attack, but SbCCR2-2 was significantly induced by sorghum aphid attack. It is suggested that SbCCR2-2 is involved in the process of pest defense. Absolute quantitative real-time PCR revealed that the three CCR genes were mainly expressed in lignin deposition organs. The gene copy number of SbCCR1 was significantly higher than those of SbCCR2-1 and SbCCR2-2 in the tested tissues, especially in stem. The results provide new insight into the functions of the three CCR genes in sorghum.