Plant glycoside hydrolase family 9 (GH9) comprises typical endo-β-1,4-glucanase (EGases, EC3.2.1.4). Although GH9A (KORRIGAN) family genes have been reported to be involved in cellulose biosynthesis in plants, much remains unknown about other GH9 subclasses. In this study, we observed a global gene co-expression profiling and conducted a correlation analysis between OsGH9 and OsCESA among 66 tissues covering most periods of life cycles in 2 rice varieties. Our results showed that OsGH9A3 and B5 possessed an extremely high co-expression with OsCESA1, 3, and 8 typical for cellulose biosynthesis in rice. Using two distinct rice non-GH9 mutants and wild type, we performed integrative analysis of gene expression level by qRT-PCR, cellulase activities in situ and in vitro, and lignocellulose crystallinity index (CrI) in four internodes of stem tissues. For the first time, OsGH9B1, 3, and 16 were characterized with the potential role in lignocellulose crystallinity alteration in rice, whereas OsGH9A3 and B5 were suggested for cellulose biosynthesis. In addition, phylogenetic analysis and gene co-expression comparison revealed GH9 function similarity in Arabidopsis and rice. Hence, the data can provide insights into GH9 function in plants and offer the potential strategy for genetic manipulation of plant cell wall using the five aforementioned novel OsGH9 genes.