The formation of β1,3-linkages on animal glycoconjugates is catalyzed by a subset of β1,3-glycosyltransferases grouped in the CAZy family GT31. This family represents an extremely diverse set of β1,3-N-acetylglucosaminyltransferases (B3GNTs and FNGs), β1,3-N-acetylgalactosaminyltransferases (B3GALNTs), β1,3-galactosyltransferases (B3GALTs and C1GALTs), β1,3-glucosaminyltransferase (B3GLCT), and β1,3-glucuronyl acid transferases (B3GLCATs or CHs). The mammalian enzymes were particularly well-studied and shown to use a large variety of sugar donors and acceptor substrates leading to the formation of β1,3-linkages in various glycosylation pathways. In contrast, there are only a few studies related to other metazoan and lower vertebrates GT31 enzymes and the evolutionary relationships of these divergent sequences remain obscure. In this study, we used bioinformatics approaches to identify more than 920 of putative GT31 sequences in Metazoa, Fungi, and Choanoflagellata revealing their deep ancestry. Sequence-based analysis shed light on conserved motifs and structural features that are signatures of all the GT31. We leverage pieces of evidence from gene structure, phylogenetic and sequence-based analyses to identified two major subgroups of GT31 named FR and BGR and demonstrate the existence of ten orthologue groups in the Urmetazoa, the hypothetical last common ancestor of all animals. Finally, synteny and paralogy analysis unveiled the existence of 30 subfamilies in vertebrates, among which 5 are new and were named C1GALT2, C1GALT3, B3GALT8, B3GNT10, and B3GNT11. Altogether, these various approaches enabled us to propose the first comprehensive analysis of the metazoan GT31 disentangling their evolutionary relationships.