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

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