Abstract
The hydrolysis of sugar-containing compounds by glycoside hydrolases (GHs) plays essential roles in many major biological processes, but to date our systematic understanding of the functional diversity and evolution of GH families remains largely limited to a few well-studied terrestrial animals. Molluscs represent the largest marine phylum in the animal kingdom, and many of them are herbivorous that utilize algae as a main nutritional source, making them good subjects for studying the functional diversity and adaptive evolution of GH families. In the present study, we conducted genome-wide identification and functional and evolutionary analysis of all GH families across major molluscan lineages. We revealed that the remarkable expansion of the GH9, GH10, GH18 and GH20 families and the wide adoption of carbohydrate-binding modules in molluscan expanded GH families likely contributed to the efficient hydrolysis of marine algal polysaccharides and were involved in the consolidation of molluscan algae-feeding habits. Gene expression and network analysis revealed the hepatopancreas as the main organ for the prominent expression of approximately half of the GH families (well corresponding to the digestive roles of the hepatopancreas) and key or hub GHs in the coexpression gene network with potentially diverse functionalities. We also revealed the evolutionary signs of differential expansion and functional divergence of the GH family, which possibly contributed to lineage-specific adaptation. Systematic analysis of GH families at both genomic and transcriptomic levels provides important clues for understanding the functional divergence and evolution of GH gene families in molluscs in relation to their algae-feeding biology.