Genome-wide survey of the phosphofructokinase family in cassava and functional characterization in response to oxygen-deficient stress

Abstract Background Glycolytic pathway is common in all plant organs, especially in oxygen-deficient tissues. Phosphofructokinase (PFK) is a rate-limiting enzyme in the glycolytic pathway and catalyses the phosphorylation of fructose-6-phosphate to fructose-1,6-bisphosphate. Cassava (M. esculenta) root is a huge storage organ with low amount of oxygen. However, less is known about the functions of PFK from M. esculenta (MePFK). We conducted a systematic analysis of MePFK genes to explore the function of the MePFK gene family under hypoxic stress. Results We identified 13 MePFK genes and characterised their sequence structure. The phylogenetic tree divided the 13 genes into two groups: nine were MePFKs and four were pyrophosphate-fructose-6-phosphate phosphotransferase (MePFPs). We confirmed by green fluorescent protein fusion protein expression that MePFK03 and MePFPA1 were localised in the chloroplast and cytoplasm, respectively. The expression profiles of the 13 MePFKs detected by quantitative reverse transcription polymerase chain reaction revealed that MePFK02, MePFK03, MePFPA1, MePFPB1 displayed higher expression in leaves, root and flower. The expression of MePFK03, MePFPA1 and MePFPB1 in tuber root increased gradually with plant growth. We confirmed that hypoxia occurred in the cassava root, and the concentration of oxygen was sharply decreasing from the outside to the inside root. The expression of MePFK03, MePFPA1 and MePFPB1 decreased with the decrease in the oxygen concentration in cassava root. Waterlogging stress treatment showed that the transcript level of PPi-dependent MePFP and MeSuSy were up-regulated remarkably and PPi-dependent glycolysis bypass was promoted. Conclusion A systematic survey of phylogenetic relation, molecular characterisation, chromosomal and subcellular localisation and cis-element prediction of MePFKs were performed in cassava. The expression profiles of MePFKs in different development stages, organs and under waterlogging stress showed that MePFPA1 plays an important role during the growth and development of cassava. Combined with the transcriptional level of MeSuSy, we found that pyrophosphate (PPi)-dependent glycolysis bypass was promoted when cassava was under waterlogging stress. The results would provide insights for further studying the function of MePFKs under hypoxic stress