Glycan Profiling by Sequencing to Uncover Multicellular Communication: Launching Glycobiology in Single Cells and Microbiomes
Lalhaba Oinam
Hiroaki Tateno- Cellular and Molecular Biotechnology Research Institute, National Institute of Advanced Industrial Science and Technology, Ibaraki, Japan
Glycans are essential building blocks of life that are located at the outermost surface of all cells from mammals to bacteria and even viruses. Cell surface glycans mediate multicellular communication in diverse biological processes and are useful as “surface markers” to identify cells. Various single-cell sequencing technologies have already emerged that enable the high-throughput analysis of omics information, such as transcriptome and genome profiling on a cell-by-cell basis, which has advanced our understanding of complex multicellular interactions. However, there has been no robust technology to analyze the glycome in single cells, mainly because glycans with branched and heterogeneous structures cannot be readily amplified by polymerase chain reactions like nucleic acids. We hypothesized that the generation of lectins conjugated with DNA barcodes (DNA-barcoded lectins) would enable the conversion of glycan information to gene information, which may be amplified and measured using DNA sequencers. This technology will enable the simultaneous analysis of glycan and RNA in single cells. Based on this concept, we developed a technology to analyze glycans and RNA in single cells, which was referred to as scGR-seq. Using scGR-seq, we acquired glycan and gene expression profiles of individual cells constituting heterogeneous cell populations, such as tissues. We further extended Glycan-seq to the profiling of the surface glycans of bacteria and even gut microbiota. Glycan-seq and scGR-seq are new technologies that enable us to elucidate the function of glycans in cell–cell and cell–microorganism communication, which extends glycobiology to the level of single cells and microbiomes.
Introduction
The human body is composed of 37 trillion cells, which are covered in a dense layer of glycans with a diameter of ∼30 nm (Martinez-Palomo et al., 1969). Cell surface glycans have been determined to play important roles in diverse biological processes, including development (Haltiwanger and Lowe, 2004), pluripotency (Alisson-Silva et al., 2014; Nishihara, 2018), tumorigenesis (Dube and Bertozzi, 2005; Lau and Dennis, 2008), and immune escape (van Kooyk and Rabinovich, 2008). Cell surface glycans are known to vary depending upon the cell type and state, such as the degree of differentiation and tumorigenesis (Tateno et al., 2011; Hasehira et al., 2012; Taniguchi and Kizuka, 2015; Shimomura et al., 2018). Therefore, glycans are often referred to as a “cell signature” that reflects cellular characteristics (Naoki et al., 2013; Krištić et al., 2014). Most glycans are attached to proteins or lipids to form glycoconjugates, such as glycoproteins, glycolipids, proteoglycans, and glycosylphosphatidylinositol (GPI)-anchored proteins (Gahmberg and Tolvanen, 1996; Kinoshita et al., 1997; Maccioni et al., 2002). Glycans are composed of approximately ten different monosaccharides (Glc, GlcNAc, Gal, GalNAc, Man, Xyl, GlcA, Fuc, NeuAc, and IdoA), which are linked by glycosidic bonds (Cummings and Pierce, 2014; Seeberger, 2015). Anomers and linkage isomers exist in the glycan structure; thus a variety of isomers are possible (Laine, 1994; Werz et al., 2007). Unlike genes and proteins existing as a linear sequence, glycans typically exhibit a branched structure, which is regulated by the substrate specificity of the glycosyltransferases (Varki, 2017). Glycan structures are evolving rapidly, and they are different depending on the organismal species (Varki, 2006; Varki, 2017). Glycans are the secondary products of genes, which are synthesized by the activity of various glycogenes (as of a total >200), such as glycosyltransferases, glycosidases, sugar-nucleotide transporter synthases, and sugar-nucleotide transporters (Yarema and Bertozzi, 2001). As a result, glycan structures may be influenced by intrinsic and extrinsic environmental changes and cannot be easily predicted simply from gene expression profiles. Therefore, it is necessary to develop technologies to analyze cell surface glycans directly.
Lectin-Based Glycan Profiling
Various strategies have been developed to analyze the glycome, which represents the total set of glycans expressed in a cell or tissue. These techniques include mass spectrometry (MS), high-performance liquid chromatography (HPLC), nuclear magnetic resonance (NMR), and capillary electrophoresis (CE) (Haslam et al., 2006; Yamaguchi and Kato, 2010; Nakano et al., 2011). In 2005, a lectin-based glycan profiling technology, known as lectin microarray, emerged. Since then, it has been noted to play a pivotal role in surveying and mapping the structure of complex glycans in various biological samples (Kuno et al., 2005; Pilobello et al., 2005; Narimatsu et al., 2010; Narimatsu et al., 2018; Hirabayashi et al., 2013; Ribeiro and Mahal, 2013). In lectin microarrays, specific lectins with various glycan-binding specificities, which can discriminate structural isomers such as anomers and linkage isomers of glycans, are immobilized onto glass slides. Glycoproteins extracted from cells or tissues are incubated with the lectin microarray, wherein the binding profiles of lectins can be acquired. Because the specificity of each lectin is known, the glycan profiles and their differences in cells and tissues can be predicted. Lectin microarrays have also been applied to the analysis of not only glycoproteins but also small vesicles, such as exosomes (Gerlach et al., 2013; Shimoda et al., 2017; Shimoda et al., 2022; Saito et al., 2018); viruses (Stevens et al., 2006; Hiono et al., 2019); whole cells, such as live mammalian cells (Zheng et al., 2005; Lee et al., 2006; Tateno et al., 2007; Tao et al., 2008); and bacterial cells (Hsu et al., 2006; Gao et al., 2010; Yasuda et al., 2011). However, there are limitations to lectin microarrays as well as these other analytical methods. For example, 1) glycans cannot be analyzed at the single-cell level, 2) the glycan profile of each cell type in a mixed cell population cannot be obtained without prior separation, and 3) the relationship between the glycome and transcriptome in single cells cannot be determined.
Concept of Glycan Profiling by Sequencing
High-throughput single-cell sequencing has been transformative for the identification and study of complex cell populations (Stuart and Satija, 2019). Recently, simultaneous profiling of multiple types of molecules within a single cell has been developed in order to establish a more comprehensive molecular view of the cell (Peterson et al., 2017; Stoeckius et al., 2017; Stuart and Satija, 2019). However, there are no techniques yet to simultaneously analyze the glycome and transcriptome in single cells. One reason is that, unlike DNA and RNA, glycans cannot be amplified by methods such as the polymerase chain reaction (PCR). Therefore, we hypothesized that the generation of lectins conjugated with a DNA sequence (DNA-barcoded lectin) would enable the transformation of glycan information to gene information, which may be amplified and further measured using a DNA sequencer, such as a next-generation sequencer (Figure 1A). In addition, the simultaneous analysis of the glycome with other molecular profiles, such as the transcriptome, may be realized.
FIGURE 1. Single-cell glycan and RNA sequencing (scGR-seq). (A) Principle of converting glycan information into gene information by DNA-barcoded lectins. (B) Schematic experimental workflow of scGR-seq. (C) hiPSCs after 0- (red), 4- (green), and 11-day differentiation (blue) into neural progenitor cells were analyzed by scGlycan-seq (left panel), flow cytometry (middle panel), and principal component analysis (right panel). (D) Dimensional reduction and clustering. UMAP visualization based on only the scRNA-seq data (left panel), only the scGlycan-seq (middle panel), both scRNA-seq and scGlycan-seq (scGR-seq, right panel) data of hiPSCs (n = 53, red), and NPCs (n = 43, green). (E) PLS regression. A heatmap showing the association between each lectin and each component inferred by PLS regression. Rows represent lectins, and columns represent components. (F) Correlation between lectin signal and glycosyltransferase gene expression. A heatmap showing the association of glycogenes and lectins inferred by PLS regression. Rows represent genes, and columns represent lectins. Figures are reprinted from Minoshima et al. (2021) (Minoshima et al., 2021).
Glycan-seq
In 2016, we began developing a technology to analyze glycans in single cells based on the aforementioned concept as a Japan Science and Technology Agency (JST) PRESTO project. Lectins with known specificity were conjugated to DNA oligonucleotides containing a barcode sequence for the identification of the specific lectin, thus enabling the lectins to be specifically identified by sequence analysis (Figure 1A) (Minoshima et al., 2021). The lectins were conjugated via their amino groups with the photocleavable dibenzocyclooctyne-N-hydroxysuccinimidyl ester (DBCO-NHS), which allowed the efficient conjugation with 5′-azide-modified oligonucleotides. The oligonucleotides were released from the lectin following ultraviolet (UV) exposure (Odaka et al., 2022). We prepared a panel of 39 DNA-barcoded lectins that covered various glycans, such as sialylated, galactosylated, GlcNAcylated, mannosylated, and fucosylated glycans that are present in glycoconjugates (Minoshima et al., 2021), whereas DNA-barcoded mouse and goat IgG were used as negative controls. In total, 41 DNA-barcoded proteins were incubated with 1 × 105 cells, and the unbound lectins were removed by washing (Odaka et al., 2022) (Figure 1B). Then, the bulk or single cells were separated into a PCR tube and exposed to UV light. After centrifugation, the supernatants containing the released DNA barcodes were recovered, amplified by PCR, and analyzed by a next-generation sequencer to count the DNA barcodes (Odaka et al., 2022). We refer to this method as Glycan-seq.
We evaluated the ability of Glycan-seq in a comparative analysis of the following bulk samples: human-induced pluripotent stem cells (hiPSCs) vs. human dermal fibroblasts (hFibs), Chinese hamster ovary cells vs. glycosylation-defective Chinese hamster ovary mutants, and hiPSCs vs. hiPSC-derived neural progenitor cells (NPCs) (Minoshima et al., 2021). The results were compared by flow cytometry using fluorescence-labeled lectins as the gold standard. Essentially, the Glycan-seq data were consistent with the flow cytometry data (Minoshima et al., 2021). Therefore, bulk Glycan-seq can capture distinct and quantitative differences in glycan profiles in various cell populations as confirmed via flow cytometry.
Next, we tested the applicability of Glycan-seq in single cells, which we termed single-cell Glycan-seq (scGlycan-seq) (Minoshima et al., 2021) (Figure 1C). We applied scGlycan-seq for comparative analysis of hiPSCs and hFibs and hiPSCs before and after differentiation into NPCs. The relative quantitative differences in the rBC2LCN signal for hiPSCs before (day 0) and after differentiation to NPCs (days 4 and 11) observed by flow cytometry were also captured by scGlycan-seq (Minoshima et al., 2021) (Figure 1C, left and middle panels). The principal component analysis clearly separated single cells on days 0, 4, and 11, and the cells were clearly ordered with respect to the progression of differentiation (Minoshima et al., 2021) Figure 1C, right panel). Therefore, scGlycan-seq enabled glycan profiling in single cells and revealed cellular heterogeneity in the glycan profiles.
scGR-seq
scGlycan-seq was then combined with scRNA-seq for the simultaneous analysis of glycan and RNA profiles in single cells (scGR-seq) (Minoshima et al., 2021) (Figure 1B). For scRNA-seq, we used a plate-based method known as RamDA-seq, which is a full-length single-cell total RNA-sequencing method (Hayashi et al., 2018). We performed scGR-seq on human-induced pluripotent stem cells (hiPSCs) and hiPSC-derived NPCs (11-days differentiation). Using UMAP, a nonlinear dimensional clustering based on only the mRNA or glycan data, the two cell types (hiPSCs and NPCs) were partially separated (Minoshima et al., 2021) (Figure 1D). In contrast, when we performed UMAP based on both the mRNA and glycan data, the two cell types were clearly separated (Minoshima et al., 2021) (Figure 1D). Therefore, the combination of mRNA and glycan profiling techniques has further characterized the cell identities. Simultaneous transcriptome and glycome profiles can associate genes with glycans at the single-cell level. A PLS regression analysis was able to identify a group of mRNAs and lectins that were associated with one another differently per component (Minoshima et al., 2021) (Figure 1E). This analysis allowed us to infer each glycan’s potential function and role as a marker through the set of genes associated with the glycan. We also established the overall relationship between lectins and glycosylation-related genes (Minoshima et al., 2021) (Figure 1F). Therefore, scGR-seq is useful for finding potential relevance between the transcriptome and glycome profiles.
Glycan-seq of the Gut Microbiota
The gut microbiota is known to be populated with diverse microbial communities, of which the bacterial communities are present in high numbers (Costello et al., 2009; Neish, 2009). It is estimated that the human gastrointestinal tract is home to approximately 100 trillion (1014) microbes, including 1,000 species of bacteria, which is similar to the number of cells in the entire human body (Sender et al., 2016). The interactions between these microbial communities and the host provide important physiological functions that can affect human health (Sekirov et al., 2010; Kho and Lal, 2018; Valdes et al., 2018; Fan and Pedersen, 2021). Mammalian and bacterial cells are coated with glycans, which serve as an interface for crosstalk with the host (Mescher et al., 1974; Comstock and Kasper, 2006). Bacterial cell surface glycans are highly complex and quite different from those of eukaryotes: Gram-positive bacterial cells are enclosed by a single membrane covered by a thick peptidoglycan layer and lipoteichoic acids (Whitfield, 1988; Rajagopal and Walker, 2017), whereas Gram-negative bacterial cells are covered by two cell membranes (inner and outer membranes) separated by a periplasm containing a thin peptidoglycan layer and the outer membrane consisting of lipopolysaccharides (Figure 2A) (Beveridge, 1999; Casey et al., 2008). Understanding the cell surface glycans of gut bacteria may provide better insight into the interactions between the host and the gut microbiota. Previously, we and others used lectin microarrays to demonstrate that the cell surface glycans of bacteria are different from strain to strain and with culture conditions (Hsu et al., 2006; Gao et al., 2010; Yasuda et al., 2011). Utilizing lectin microarrays for glycan analysis of the gut microbiota is challenging for the following reasons: 1) it is difficult to fluorescently label all of the bacteria comprising the microbiota with the same intensity, 2) a large number of cells is required, and 3) bacterial cells are frequently washed out during the washing steps. Therefore, there has been no robust technology for glycan profiling of the gut microbiota.
FIGURE 2. Glycan profiling of the gut microbiota and 16S rRNA sequencing. (A) Illustration of the typical cell wall architecture of Gram-positive and negative bacteria. (B) Schematic experimental workflow of glycomic profiling, bacterial composition analysis, and the identification of lectin-reactive bacteria. (C) Hierarchical clustering heatmap of the gut microbiota of mouse pups (n = 3) and adult mice (n = 3) obtained from the Glycan-seq data. The column shows the pups and adult mouse sample separation, and the row shows the name of the lectins used in the Glycan-seq analysis. (D) The stacked bar graph represents the differential abundance of the bacterial family identified by 16S rRNA sequencing from each sample. Each colored bar represents the bacterial family identified. Figures are reprinted from Oinam et al. (2022) (Oinam et al., 2022).
We began developing technology for glycan profiling of the gut microbiota as an AMED-Prime project in 2018. For this purpose, the Glycan-seq technique described earlier was modified and implemented for the analysis of bacterial surface glycans. We first evaluated the applicability of Glycan-seq to bacterial glycan profiling in cultured bacteria: Gram-positive (Deinococcus radiodurans) vs. Gram-negative (Escherichia coli) bacteria (Oinam et al., 2022). After incubation of bacteria with DNA-barcoded lectins, the barcodes were released from the lectins bound to the bacteria by UV irradiation. The recovered DNA barcodes were amplified by PCR, and the number of barcodes derived from each lectin was counted by a next-generation sequencer, which corresponds to the resulting bacterial glycan profiles (Oinam et al., 2022). The Glycan-seq analysis revealed that the glycans were different between Gram-positive and Gram-negative bacteria. To confirm the Glycan-seq data, lectins differentially detected in the cultured bacteria were subjected to flow cytometry analysis, the gold standard for cellular lectin staining. The results indicated that the identified lectins are bound to the bacteria, which is in agreement with the Glycan-seq data and suggested that Glycan-seq data may be applied to bacterial cell surface glycan analysis (Oinam et al., 2022). The specific lectin binding was confirmed using a competition assay with saccharide inhibitors and glycosidase treatment, followed by flow cytometry.
After confirming the applicability of glycan analysis to cultured bacteria, Glycan-seq was used to analyze the gut microbiota along with genomic profiling of the bulk gut microbiota from pups (14–20 days old) and adult (12-month old) mice (Figure 2B) (Oinam et al., 2022). Glycan analysis revealed that the glycans of mouse pups and adult mice were different (Figure 2B) and lectins identifying sialylated glycans were higher in pups (Figure 2C). The genome content was isolated from the bacterial cells, and the V3-V4 region of the 16S rRNA was amplified by PCR (Oinam et al., 2022). The 16S rRNA was then sequenced to determine the bacterial composition of the gut microbiota (Oinam et al., 2022). The composition of the two microbiota was different as bacterial families belonging to Lactobacillaceae, Enterobacteriaceae, Pseudomonadaceae, and Gemellaceae were more abundant in pups (Figure 2D) (Oinam et al., 2022). Sialylated bacteria were then enriched by lectin pulldown using Sia-binding lectins and subjected to 16S rRNA sequencing. The sialylated bacteria identified were Lactobacillaceae, Lachnospiraceae, Enterobacteriaceae, and Muribaculaceae (Oinam et al., 2022). Using Glycan-seq, we were able to perform a glycan analysis of the gut microbiota along with the bacterial composition using the same sequencing instrument. Therefore, Glycan-seq analysis comprehensively revealed differences in the glycan profile of the gut microbiota of pups and adult mice and identified more sialylated bacteria in the mouse pups.
Concluding Remarks and Future Perspective
scGR-seq provides lectin-based glycan and gene expression profiles for individual cells, making it possible to obtain detailed glycan information and discriminate structural isomers of glycans on single cells constituting a tissue. These data will provide insight into complex multicellular communication networks, including tumor microenvironments and neural networks based on lectin-receptor interactions. scGR-seq can also be applied to the development of drug targets for rare cells, such as cancer stem cells and circulating tumor cells. However, there are limitations to the current Glycan-seq and scGR-seq techniques. Similar to flow cytometry and lectin microarray, absolute amounts of glycans cannot be determined. Another limitation of this current system is the throughput. Since scGR-seq is a plate-based platform, the processing is currently limited to hundreds of cells, whereas it can perform full-length total RNA sequencing. In contrast, droplet-based methods such as 10x Genomics (CITE-seq) can sequence thousands of cells at once but only target the 3′ ends of poly(A) transcripts (Baran-Gale et al., 2018). Because of this difference, scGR-seq will complement the study of single cells in complex biological systems. To resolve this limitation, we plan to improve scGR-seq and adapt it to a droplet-based high-throughput single-cell technology (B Rosenberg et al., 2018). Accordingly, SUrface-protein Glycan and RNA-seq (SUGAR-seq) based on the 10x Genomics platform was also recently reported for the detection of a lectin-binding signal together with the analysis of extracellular epitopes and the transcriptome at the single-cell level (Kearney et al., 2022), although SUGAR-seq detects only one lectin binding to a single cell.
We have also adopted Glycan-seq to approach the untapped glycomics of the gut microbiota, which mediates the direct crosstalk with the host. Bacteria containing a particular glycan surface marker may represent a novel diagnostic and therapeutic target of the disease. For glycomic profiling of the gut microbiota, only bulk analysis by Glycan-seq is currently available. The development of a technique for the simultaneous analysis of the glycome and genome in single cells is needed to fully realize bacterial glycomic profiling in single cells.
In conclusion, we have developed a lectin-based glycan profiling technique by sequencing and applied this technique to the joint analysis of glycan and RNA in single cells and the glycomic profiling of the gut microbiota. Glycan-seq and scGR-seq have the potential to advance our understanding of cellular heterogeneity and the biological role of glycans across diverse multicellular systems across species and lead to the launch of glycobiology in single cells and microbiomes.
Author Contributions
LO and HT wrote the paper.
Funding
This work was supported by JST PRESTO (JPMJPR16F6), JST A-STEP (JPMJTR20UD), AMED-Prime (21gm6010018h0004), JST-CREST (JPMJCR17H2), Mizutani Foundation for Glycoscience, Yakult Bio-Science Foundation, fellowship and Grant-in-Aid from the Tokyo Biochemical Research Foundation, and The Kurata Grants from the Hitachi Global Foundation.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank Fumi Minoshima, Haruki Odaka, Jinko Murakami, Sayoko Saito, Kayo Kiyoi, Keiko Hiemori, Kayo Suzuki, and Sunanda Keisham at the National Institute of Advanced Industrial Science and Technology and Haruka Ozaki at the University of Tsukuba for their contribution to the development of Glycan-seq and scGR-seq.
Abbreviations
CE, capillary electrophoresis; HPLC, high-performance liquid chromatography; MS, mass spectrometry; NMR, nuclear magnetic resonance; NPC, neural progenitor cells; and PCR, polymerase chain reaction.
References
Alisson-Silva, F., de Carvalho Rodrigues, D., Vairo, L., Asensi, K. D., Vasconcelos-dos-Santos, A., Mantuano, N. R., et al. (2014). Evidences for the Involvement of Cell Surface Glycans in Stem Cell Pluripotency and Differentiation. Glycobiology 24, 458–468. doi:10.1093/glycob/cwu012
Baran-Gale, J., Chandra, T., and Kirschner, K. (2018). Experimental Design for Single-Cell RNA Sequencing. Brief. Funct. Genomics 17, 233–239. doi:10.1093/bfgp/elx035
Beveridge, T. J. (1999). Structures of Gram-Negative Cell Walls and Their Derived Membrane Vesicles. J. Bacteriol. 181, 4725–4733. doi:10.1128/JB.181.16.4725-4733.1999
Huang, K. C., Mukhopadhyay, R., Wen, B., Gitai, Z., and Wingreen, N. S. (2008). Cell Shape and Cell-Wall Organization in Gram-Negative Bacteria. Proc. Natl. Acad. Sci. U.S.A. 105, 19282–19287. doi:10.1073/pnas.0805309105
Comstock, L. E., and Kasper, D. L. (2006). Bacterial Glycans: Key Mediators of Diverse Host Immune Responses. Cell 126, 847–850. doi:10.1016/j.cell.2006.08.021
Costello, E. K., Lauber, C. L., Hamady, M., Fierer, N., Gordon, J. I., and Knight, R. (2009). Bacterial Community Variation in Human Body Habitats Across Space and Time. Science 326, 1694–1697. doi:10.1126/science.1177486
Cummings, R. D., and Pierce, J. M. (2014). The Challenge and Promise of Glycomics. Chem. Biol. 21, 1–15. doi:10.1016/j.chembiol.2013.12.010
Dube, D. H., and Bertozzi, C. R. (2005). Glycans in Cancer and Inflammation - Potential for Therapeutics and Diagnostics. Nat. Rev. Drug Discov. 4, 477–488. doi:10.1038/nrd1751
Fan, Y., and Pedersen, O. (2021). Gut Microbiota in Human Metabolic Health and Disease. Nat. Rev. Microbiol. 19, 55–71. doi:10.1038/s41579-020-0433-9
Gahmberg, C. G., and Tolvanen, M. (1996). Why Mammalian Cell Surface Proteins Are Glycoproteins. Trends Biochem. Sci. 21, 308–311. doi:10.1016/S0968-0004(96)10034-7
Gao, J., Liu, D., and Wang, Z. (2010). Screening Lectin-Binding Specificity of Bacterium by Lectin Microarray with Gold Nanoparticle Probes. Anal. Chem. 82, 9240–9247. doi:10.1021/ac1022309
Gerlach, J. Q., Krüger, A., Gallogly, S., Hanley, S. A., Hogan, M. C., Ward, C. J., et al. (2013). Surface Glycosylation Profiles of Urine Extracellular Vesicles. PLoS One 8, e74801. doi:10.1371/journal.pone.0074801
Haltiwanger, R. S., and Lowe, J. B. (2004). Role of Glycosylation in Development. Annu. Rev. Biochem. 73, 491–537. doi:10.1146/annurev.biochem.73.011303.074043
Hasehira, K., Tateno, H., Onuma, Y., Ito, Y., Asashima, M., and Hirabayashi, J. (2012). Structural and Quantitative Evidence for Dynamic Glycome Shift on Production of Induced Pluripotent Stem Cells. Mol. Cell. Proteomics 11, 1913–1923. doi:10.1074/mcp.M112.020586
Haslam, S., North, S., and Dell, A. (2006). Mass Spectrometric Analysis of N- and O-Glycosylation of Tissues and Cells. Curr. Opin. Struct. Biol. 16, 584–591. doi:10.1016/J.SBI.2006.08.006
Hayashi, T., Ozaki, H., Sasagawa, Y., Umeda, M., Danno, H., and Nikaido, I. (2018). Single-cell Full-Length Total RNA Sequencing Uncovers Dynamics of Recursive Splicing and Enhancer RNAs. Nat. Commun. 9, 619. doi:10.1038/s41467-018-02866-0
Hiono, T., Matsuda, A., Wagatsuma, T., Okamatsu, M., Sakoda, Y., and Kuno, A. (2019). Lectin Microarray Analyses Reveal Host Cell-specific Glycan Profiles of the Hemagglutinins of Influenza A Viruses. Virology 527, 132–140. doi:10.1016/j.virol.2018.11.010
Hirabayashi, J., Yamada, M., Kuno, A., and Tateno, H. (2013). Lectin Microarrays: Concept, Principle and Applications. Chem. Soc. Rev. 42, 4443–4458. doi:10.1039/C3CS35419A
Hsu, K.-L., Pilobello, K. T., and Mahal, L. K. (2006). Analyzing the Dynamic Bacterial Glycome with a Lectin Microarray Approach. Nat. Chem. Biol. 2, 153–157. doi:10.1038/nchembio767
Kearney, C. J., Vervoort, S. J., Ramsbottom, K. M., Todorovski, I., Lelliott, E. J., Zethoven, M., et al. (2021). SUGAR-seq Enables Simultaneous Detection of Glycans, Epitopes, and the Transcriptome in Single Cells. Sci. Adv. 7, eabe3610. doi:10.1126/sciadv.abe3610
Kho, Z. Y., and Lal, S. K. (2018). The Human Gut Microbiome - A Potential Controller of Wellness and Disease. Front. Microbiol. 9, 1835. doi:10.3389/fmicb.2018.01835
Kinoshita, T., Ohishi, K., and Takeda, J. (1997). GPI-Anchor Synthesis in Mammalian Cells: Genes, Their Products, and a Deficiency. J. Biochem. 122, 251–257. doi:10.1093/oxfordjournals.jbchem.a021746
Krištić, J., Vučković, F., Menni, C., Klarić, L., Keser, T., Beceheli, I., et al. (2014). Glycans are a Novel Biomarker of Chronological and Biological Ages. Journals Gerontol. - Ser. A Biol. Sci. Med. Sci. 69, 779–789. doi:10.1093/gerona/glt190
Kuno, A., Uchiyama, N., Koseki-Kuno, S., Ebe, Y., Takashima, S., Yamada, M., et al. (2005). Evanescent-field Fluorescence-Assisted Lectin Microarray: a New Strategy for Glycan Profiling. Nat. Methods 2, 851–856. doi:10.1038/nmeth803
Laine, R. A. (1994). Invited Commentary: A Calculation of All Possible Oligosaccharide Isomers Both Branched and Linear Yields 1.05 × 1012 Structures for a Reducing Hexasaccharide: the Isomer Barrier to Development of Single-Method Saccharide Sequencing or Synthesis Systems. Glycobiology 4, 759–767. doi:10.1093/glycob/4.6.759
Lau, K. S., and Dennis, J. W. (2008). N-Glycans in Cancer Progression. Glycobiology 18, 750–760. doi:10.1093/glycob/cwn071
Lee, M.-r., Park, S., and Shin, I. (2006). Protein Microarrays to Study Carbohydrate-Recognition Events. Bioorg. Med. Chem. Lett. 16, 5132–5135. doi:10.1016/j.bmcl.2006.07.028
Maccioni, H. J. F., Giraudo, C. G., and Daniotti, J. L. (2002). Understanding the Stepwise Synthesis of Glycolipids. Neurochem. Res. 27, 629–636. doi:10.1023/A:1020271932760
Martinez-Palomo, A., Braislovsky, C., and Bernhard, W. (1969). Ultrastructural Modifications of the Cell Surface and Intercellular Contacts of Some Transformed Cell Strains. Cancer Res. 29, 925–937.
Mescher, M. F., Strominger, J. L., and Watson, S. W. (1974). Protein and Carbohydrate Composition of the Cell Envelope of Halobacterium Salinarium. J. Bacteriol. 120, 945–954. doi:10.1128/jb.120.2.945-954.1974
Minoshima, F., Ozaki, H., Odaka, H., and Tateno, H. (2021). Integrated Analysis of Glycan and RNA in Single Cells. iScience 24, 102882. doi:10.1016/j.isci.2021.102882
Nakano, M., Kakehi, K., Taniguchi, N., and Kondo, A. (2011). in Capillary Electrophoresis and Capillary Electrophoresis–Mass Spectrometry for Structural Analysis of N-Glycans Derived from Glycoproteins BT - Capillary Electrophoresis of Carbohydrates: From Monosaccharides to Complex Polysaccharides. Editor N. Volpi (Totowa, NJ: Humana Press), 205–235. doi:10.1007/978-1-60761-875-1_9
Fujitani, N., Furukawa, J.-i., Araki, K., Fujioka, T., Takegawa, Y., Piao, J., et al. (2013). Total Cellular Glycomics Allows Characterizing Cells and Streamlining the Discovery Process for Cellular Biomarkers. Proc. Natl. Acad. Sci. U.S.A. 110, 2105–2110. doi:10.1073/pnas.1214233110
Narimatsu, H., Sawaki, H., Kuno, A., Kaji, H., Ito, H., and Ikehara, Y. (2010). A Strategy for Discovery of Cancer Glyco-Biomarkers in Serum Using Newly Developed Technologies for Glycoproteomics. FEBS J. 277, 95–105. doi:10.1111/j.1742-4658.2009.07430.x
Narimatsu, H., Kaji, H., Vakhrushev, S. Y., Clausen, H., Zhang, H., Noro, E., et al. (2018). Current Technologies for Complex Glycoproteomics and Their Applications to Biology/Disease-Driven Glycoproteomics. J. Proteome Res. 17, 4097–4112. doi:10.1021/ACS.JPROTEOME.8B00515
Neish, A. S. (2009). Microbes in Gastrointestinal Health and Disease. Gastroenterology 136, 65–80. doi:10.1053/j.gastro.2008.10.080
Nishihara, S. (2018). Glycans in Stem Cell Regulation: from Drosophila Tissue Stem Cells to Mammalian Pluripotent Stem Cells. FEBS Lett. 592, 3773–3790. doi:10.1002/1873-3468.13167
Odaka, H., Ozaki, H., and Tateno, H. (2022). scGR-Seq: Integrated Analysis of Glycan and RNA in Single Cells. Star. Protoc. 3, 101179. doi:10.1016/j.xpro.2022.101179
Oinam, L., Minoshima, F., and Tateno, H. (2022). Glycan Profiling of the Gut Microbiota by Glycan-Seq. Isme Commun. 2, 1. doi:10.1038/s43705-021-00084-2
Peterson, V. M., Zhang, K. X., Kumar, N., Wong, J., Li, L., Wilson, D. C., et al. (2017). Multiplexed Quantification of Proteins and Transcripts in Single Cells. Nat. Biotechnol. 35, 936–939. doi:10.1038/nbt.3973
Pilobello, K. T., Krishnamoorthy, L., Slawek, D., and Mahal, L. K. (2005). Development of a Lectin Microarray for the Rapid Analysis of Protein Glycopatterns. ChemBioChem 6, 985–989. doi:10.1002/cbic.200400403
Rajagopal, M., and Walker, S. (2015). Envelope Structures of Gram-Positive Bacteria. Curr. Top. Microbiol. Immunol. 404, 1–44. doi:10.1007/82_2015_5021
Ribeiro, J. P., and Mahal, L. K. (2013). Dot by Dot: Analyzing the Glycome Using Lectin Microarrays. Curr. Opin. Chem. Biol. 17, 827–831. doi:10.1016/J.CBPA.2013.06.009
Rosenberg, A. B., Roco, C. M., Muscat, R. A., Kuchina, A., Sample, P., Yao, Z., et al. (2018). Single-cell Profiling of the Developing Mouse Brain and Spinal Cord with Split-Pool Barcoding. Science 360, 176–182. doi:10.1126/science.aam8999
Saito, S., Hiemori, K., Kiyoi, K., and Tateno, H. (2018). Glycome Analysis of Extracellular Vesicles Derived from Human Induced Pluripotent Stem Cells Using Lectin Microarray. Sci. Rep. 8, 3997. doi:10.1038/s41598-018-22450-2
Seeberger, P. H. (2015). Monosaccharide Diversity. 3rd ed. Cold Spring Harbor (NY): Cold Spring Harbor Laboratory Press. Available at: http://europepmc.org/abstract/MED/28876863. doi:10.1101/glycobiology.3e.002
Sekirov, I., Russell, S. L., Antunes, L. C. M., and Finlay, B. B. (2010). Gut Microbiota in Health and Disease. Physiol. Rev. 90, 859–904. doi:10.1152/physrev.00045.2009
Sender, R., Fuchs, S., and Milo, R. (2016). Revised Estimates for the Number of Human and Bacteria Cells in the Body. PLOS Biol. 14, e1002533. doi:10.1371/journal.pbio.1002533
Shimoda, A., Tahara, Y., Sawada, S.-i., Sasaki, Y., and Akiyoshi, K. (2017). Glycan Profiling Analysis Using Evanescent-Field Fluorescence-Assisted Lectin Array: Importance of Sugar Recognition for Cellular Uptake of Exosomes from Mesenchymal Stem Cells. Biochem. Biophysical Res. Commun. 491, 701–707. doi:10.1016/j.bbrc.2017.07.126
Shimoda, A., Miura, R., Tateno, H., Seo, N., Shiku, H., Sawada, S. i., et al. (2022). Assessment of Surface Glycan Diversity on Extracellular Vesicles by Lectin Microarray and Glycoengineering Strategies for Drug Delivery Applications. Small methods 6, 2100785. doi:10.1002/smtd.202100785
Shimomura, O., Oda, T., Tateno, H., Ozawa, Y., Kimura, S., Sakashita, S., et al. (2018). A Novel Therapeutic Strategy for Pancreatic Cancer: Targeting Cell Surface Glycan Using rBC2LC-N Lectin-Drug Conjugate (LDC). Mol. Cancer Ther. 17, 183–195. doi:10.1158/1535-7163.MCT-17-0232
Stevens, J., Blixt, O., Paulson, J. C., and Wilson, I. A. (2006). Glycan Microarray Technologies: Tools to Survey Host Specificity of Influenza Viruses. Nat. Rev. Microbiol. 4, 857–864. doi:10.1038/nrmicro1530
Stoeckius, M., Hafemeister, C., Stephenson, W., Houck-Loomis, B., Chattopadhyay, P. K., Swerdlow, H., et al. (2017). Simultaneous Epitope and Transcriptome Measurement in Single Cells. Nat. Methods 14, 865–868. doi:10.1038/nmeth.4380
Stuart, T., and Satija, R. (2019). Integrative Single-Cell Analysis. Nat. Rev. Genet. 20, 257–272. doi:10.1038/s41576-019-0093-7
Taniguchi, N., and Kizuka, Y. (2015). “Chapter Two - Glycans and Cancer: Role of N-Glycans in Cancer Biomarker, Progression and Metastasis, and Therapeutics” in Glycosylation and Cancer. Editors R. R. Drake, and L. E Ball (cambridge: Academic Press), 11–51. doi:10.1016/bs.acr.2014.11.001
Tao, S.-C., Li, Y., Zhou, J., Qian, J., Schnaar, R. L., Zhang, Y., et al. (2008). Lectin Microarrays Identify Cell-specific and Functionally Significant Cell Surface Glycan Markers. Glycobiology 18, 761–769. doi:10.1093/glycob/cwn063
Tateno, H., Uchiyama, N., Kuno, A., Togayachi, A., Sato, T., Narimatsu, H., et al. (2007). A Novel Strategy for Mammalian Cell Surface Glycome Profiling Using Lectin Microarray. Glycobiology 17, 1138–1146. doi:10.1093/glycob/cwm084
Tateno, H., Toyota, M., Saito, S., Onuma, Y., Ito, Y., Hiemori, K., et al. (2011). Glycome Diagnosis of Human Induced Pluripotent Stem Cells Using Lectin Microarray. J. Biol. Chem. 286, 20345–20353. doi:10.1074/jbc.M111.231274
Valdes, A. M., Walter, J., Segal, E., and Spector, T. D. (2018). Role of the Gut Microbiota in Nutrition and Health. BMJ 361, k2179. doi:10.1136/bmj.k2179
van Kooyk, Y., and Rabinovich, G. A. (2008). Protein-glycan Interactions in the Control of Innate and Adaptive Immune Responses. Nat. Immunol. 9, 593–601. doi:10.1038/ni.f.203
Varki, A. (2006). Nothing in Glycobiology Makes Sense, except in the Light of Evolution. Cell 126, 841–845. doi:10.1016/j.cell.2006.08.022
Werz, D. B., Ranzinger, R., Herget, S., Adibekian, A., von der Lieth, C.-W., and Seeberger, P. H. (2007). Exploring the Structural Diversity of Mammalian Carbohydrates (“Glycospace”) by Statistical Databank Analysis. ACS Chem. Biol. 2, 685–691. doi:10.1021/cb700178s
Whitfield, C. (1988). Bacterial Extracellular Polysaccharides. Can. J. Microbiol. 34, 415–420. doi:10.1139/m88-073
Yamaguchi, Y., and Kato, K. (2010). “Chapter Fifteen - Dynamics and Interactions of Glycoconjugates Probed by Stable-Isotope-Assisted NMR Spectroscopy,” in Glycomics. Editor M. Fukuda (cambridge: Academic Press), 305–322. doi:10.1016/S0076-6879(10)78015-0
Yarema, K. J., and Bertozzi, C. R. (2001). Characterizing Glycosylation Pathways. Genome Biol. 2, reviews0004. doi:10.1186/gb-2001-2-5-reviews0004
Yasuda, E., Tateno, H., Hirabarashi, J., Iino, T., and Sako, T. (2011). Lectin Microarray Reveals Binding Profiles of Lactobacillus Casei Strains in a Comprehensive Analysis of Bacterial Cell Wall Polysaccharides. Appl. Environ. Microbiol. 77, 4539–4546. doi:10.1128/AEM.00240-11
Keywords: glycan, single cell, glycobiology, microbiome, sequencing, multicellular communication, glycan profiling
Citation: Oinam L and Tateno H (2022) Glycan Profiling by Sequencing to Uncover Multicellular Communication: Launching Glycobiology in Single Cells and Microbiomes. Front. Cell Dev. Biol. 10:919168. doi: 10.3389/fcell.2022.919168
Received: 13 April 2022; Accepted: 02 May 2022;
Published: 27 May 2022.
Edited by:
Kayvon Pedram, Janelia Research Campus, United StatesReviewed by:
Toshio Takahashi, Suntory Foundation for Life Sciences, JapanCopyright © 2022 Oinam and Tateno. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hiroaki Tateno, h-tateno@aist.go.jp