Integrated Glycoproteomics Identifies a Role of N-Glycosylation and Galectin-1 on Myogenesis and Muscle Development.

Many cell surface and secreted proteins are modified by the covalent addition of glycans that play an important role in the development of multicellular organisms. These glycan modifications enable communication between cells and the extracellular matrix via interactions with specific glycan-binding lectins and the regulation of receptor-mediated signaling. Aberrant protein glycosylation has been associated with the development of several muscular diseases, suggesting essential glycan- and lectin-mediated functions in myogenesis and muscle development, but our molecular understanding of the precise glycans, catalytic enzymes, and lectins involved remains only partially understood. Here, we quantified dynamic remodeling of the membrane-associated proteome during a time-course of myogenesis in cell culture. We observed wide-spread changes in the abundance of several important lectins and enzymes facilitating glycan biosynthesis. Glycomics-based quantification of released N-linked glycans confirmed remodeling of the glycome consistent with the regulation of glycosyltransferases and glycosidases responsible for their formation including a previously unknown digalactose-to-sialic acid switch supporting a functional role of these glycoepitopes in myogenesis. Furthermore, dynamic quantitative glycoproteomic analysis with multiplexed stable isotope labeling and analysis of enriched glycopeptides with multiple fragmentation approaches identified glycoproteins modified by these regulated glycans including several integrins and growth factor receptors. Myogenesis was also associated with the regulation of several lectins, most notably the upregulation of galectin-1 (LGALS1). CRISPR/Cas9-mediated deletion of Lgals1 inhibited differentiation and myotube formation, suggesting an early functional role of galectin-1 in the myogenic program. Importantly, similar changes in N-glycosylation and the upregulation of galectin-1 during postnatal skeletal muscle development were observed in mice. Treatment of new-born mice with recombinant adeno-associated viruses to overexpress galectin-1 in the musculature resulted in enhanced muscle mass. Our data form a valuable resource to further understand the glycobiology of myogenesis and will aid the development of intervention strategies to promote healthy muscle development or regeneration. Crown

The bulk of skeletal muscle is composed of postmitotic multinucleated myofibers that form via the fusion of mononucleated progenitor myoblasts. Myofiber formation is achieved via myogenesis, a highly ordered process including differentiation, elongation, migration, cell adhesion, membrane alignment, and ultimately cell fusion of myoblasts (1). The initial differentiation of myoblasts is regulated by external growth factors, cytokines, steroid hormones, and signal transduction pathways that activate a series of muscle-specific and pleiotropic transcription factors (2). Elongation of myoblasts is achieved by extension of filopodia and lamellipodia to contact surrounding muscle cells. Myoblasts subsequently migrate to each other, which requires extracellular matrix (ECM) remodeling to facilitate cell motility before cell recognition and adherence. Here, interactions between multiple adherence molecules trigger integrin signaling and a rearrangement of the actin-cytoskeleton. This is coupled to the regulation of several GTPases and guanine nucleotide exchange factors that contribute to membrane remodeling and cell fusion via the ARP2/3, WASP, and WAVE protein complexes (3). In mammals, distinct phases of myogenesis contribute to the formation of mature skeletal muscle (4, 5). Muscle patterning is established by the fusion of embryonic myoblasts. The second phase involves fusion of fetal myoblasts followed by the formation of the basal lamina and the expansion of adult precursor satellite cells (muscle stem cells). Accompanying this second phase is innervation of the myofibers leading to the formation of the neuromuscular junctions. Finally, postnatal myogenesis is achieved via myoblasts derived from satellite cells that are responsible for growth and regeneration of mature skeletal muscle.

Myogenesis and muscle development involve the interaction of cell surfaces and ECM with hundreds of glycosylated proteins. It is therefore not surprising that defects in glycosylation have been associated with several developmental disorders. More than 50 congenital disorders of glycosylation (CDGs) have been identified in humans, and these typically present as abnormalities in development of the nervous system and/or skeletal muscle during infancy (6). The majority of CDGs are inherited defects in one or more of the multiple enzymes responsible for glycosylation of asparagine residues (N-linked glycosylation) that occur on membrane-associated, cell surface, and secreted proteins. For example, several loss-of-function mutations have been identified in the PMM2 gene involved in the synthesis of GDP-mannose, a nucleotide-sugar donor responsible for the transfer of mannose residues to N-glycans and thus critical for their maturation. Individuals with PMM2-CDG have hypo-N-glycosylation and display severe hypotonic muscles and underdeveloped cerebellum (7). Advances in next-generation DNA sequencing have pin-pointed additional CDGs in patients presenting with abnormal N-glycosylation and defects in muscle and/or nervous system development. This includes the identification of mutations in the STT3A and STT3B genes forming the catalytic subunits of the N-oligosaccharyl transferase complex (8), and MAN1B1, which are involved in the regulation of N-linked glycosylation (9). Despite the well-documented mutations in genes regulating glycosylation and their associations with poor muscle function, we know very little about the regulation of N-glycosylation during myogenesis and muscle development. Furthermore, the roles of glycan-binding proteins such as lectins on myogenesis and muscle development remain only partially understood. For example, loss of Lgals1 expression results in defects in muscle development in mice and fish (10, 11). The molecular mechanisms remain poorly defined, but in vitro experiments suggest galectin-1 binds to a variety of ECM proteins including laminins and integrins at the sarcolemma to regulate intracellular signaling (12). Muscular dystrophies are characterized by progressive muscle degeneration, sarcolemmal damage, and loss of muscle function. Analysis of immortalized healthy and dystrophic human muscle cells with lectin histochemistry has recently revealed several changes in lectin binding, suggesting discrete changes in glycosylation or glycoprotein abundance (13). Excitingly, administration of recombinant galectin-1 in a genetic mouse model of muscular dystrophy (newborn mdx mice) displays improved muscle function and sarcolemmal integrity, suggesting critical roles of galectin-1 in myogenesis (14). These experiments are encouraging as they highlight new potential treatment options for a variety of muscle diseases. However, further experiments are required to characterize the role of N-glycosylation and galectin-1 on myogenesis and muscle development.

Experimental Procedures

L6 Cell Culture

Rat L6 myoblasts were maintained in an α-minimum essential medium containing 5.5-mM glucose (Gibco) and 10% fetal bovine serum (FBS; Hyclone Laboratories) in a 10% CO2 incubator. L6 myoblasts were differentiated into myotubes over 7 days with 2% FBS when myoblasts reached ∼90 to 95% confluency. For CRISPR/Cas9 experiments targeting complete disruption of Lgals1, myoblasts were transfected at 60% confluency with DNA constructs expressing CMV-Cas9(D10A) and paired U6-gRNAs (5’-GTTGTTGCTGTCTTTCCCCAGG and 5’- ACCCCCGCTTCAACGCCCATGG) (Sigma Aldrich) using TransIT-X2 reagent (Mirus Bio). After 48 h, cells were trypsinized, counted, and serial-diluted to ∼0.8 cells per 30 μl in fresh conditioned α-minimum essential medium media containing 10% FBS. Cells were seeded into 384-well plates with media replaced every 2 to 3 days and single clones selected. After 10 days, cells were expanded and screened for LGALS1 ablation using Western blot analysis.

L6 Cell Lysis for Mass Spectrometry Analysis

Cells were washed twice with ice-cold PBS and lysis in ice-cold 100-mM sodium carbonate containing protease inhibitor cocktail (Roche) by tip-probe sonication. Lysates were rotated at 4 °C for 1 h and then centrifuged at 150,000g for 60 min at 4 °C to pellet microsomal protein fraction. The pellet was resuspended in 6 M urea, 2 M thiourea, 1% SDS containing 25 mM triethylammonium bicarbonate (TEAB), pH 8.0, and protein precipitated with chloroform:methanol:water (1:3:4). Precipitated protein was washed with methanol and resuspended in 6 M urea and 2 M thiourea containing 25-mM TEAB, pH 8.0. The protein concentration was determined via Qubit (Invitrogen), normalized and stored at −80 °C.

(Glyco)proteomic Sample Preparation

Peptides were prepared essentially as described previously (15). Briefly, 35 μg of membrane-enriched protein from each of the eight time points (0–7 days) and each biological replicate (n = 3) were reduced with 10-mM DTT for 1 h at RT followed by alkylation with 25-mM iodoacetamide for 30 min at RT in the dark. The reaction was quenched to 25-mM DTT and digested with 0.7-μg of sequencing-grade LysC (Wako Chemicals) for 3 h at RT. The digest was diluted 5-fold with 25-mM TEAB and digested with 0.7 μg of sequencing-grade trypsin (Sigma) overnight at 30 °C. Samples were acidified to a final concentration of 2% formic acid, centrifuged at 13,000g for 10 min at RT and desalted using 30 mg of hydrophilic–lipophilic balance–solid phase extraction material in a 96-well plate (Waters) using a vacuum manifold. The plate was washed with 5% acetonitrile (MeCN) containing 0.1% TFA. The peptides were eluted with 50% MeCN containing 0.1% TFA and dried by vacuum centrifugation. Peptides were resuspended in 20 μl of 100-mM Hepes, pH 8.0, and isotopically labeled with 73 μg of 10-plex tandem mass tags (TMTs) in a final concentration of 33% MeCN for 2 h at RT. The sample channels were labeled as follows: 126 = day 0, 127N = day 1, 127C = day 2, 128N = day 3, 128C = day 4, 129N = day 5, 129C = day 6, 130N = day 7. The reactions were quenched with 0.8% hydroxylamine for 15 min at RT, and then acidified and diluted to a final concentration of 0.1% TFA and 5% MeCN. The labeled peptides were pooled and purified by hydrophilic–lipophilic balance–solid phase extraction as described above. Peptides were resuspended in 90% MeCN containing 0.1% TFA and 10% of the peptide separated directly into 12 fractions for total proteome analysis using amide-80 hydrophilic interaction liquid chromatography as described previously (16). Glycopeptides were enriched from 90% of the remaining peptide using zwitterionic-hydrophilic interaction liquid chromatography microcolumns as previously described (17, 18).

(Glyco)proteomic Mass Spectrometry

The analysis of each fraction for proteome quantification was performed on a Dionex 3500RS coupled to an Orbitrap Fusion (Thermo Scientific) operating in a positive polarity mode. Peptides were separated using an in-house packed 75-μm × 50-cm pulled column (1.9-μm particle size, C18AQ; Dr Maisch, Germany) with a gradient of 2 to 40% MeCN containing 0.1% formic acid over 150 min at 250 nl/min at 55 °C. MS1 scans were acquired from 350 to 1400 m/z (120,000 resolution, 5e5 AGC, 50 ms injection time) followed by tandem mass spectrometry (MS/MS) data-dependent acquisition of the 10 most intense ions with collision-induced dissociation (CID) and detection in the ion-trap (rapid scan rate, 1e4 AGC, 70-ms injection time, 30% NCE, 1.6 m/z isolation width). Synchronous precursor selection was enabled with multinotch isolation of the 10 most abundant fragment ions excluding precursor window of 40 m/z and loss of TMT reporter ions for MS3 analysis by higher collisional dissociation (HCD) and detection in the Orbitrap (60,000 resolution, 1e5 AGC, 300-ms injection time, 100–500 m/z, 55% NCE, 2 m/z) (19).

The analysis of glycopeptides was performed as single-shot analysis without fractionation on the identical system as described above. MS1 scans were acquired from 550 to 1750 m/z (120,000 resolution, 5e5 AGC, 100-ms injection time) followed by MS/MS data–dependent acquisition of the 7 most intense ions and highest charge state with HCD and detection in the Orbitrap (60,000 resolution, 2e5 AGC, 200-ms injection time, 40 NCE, 2 m/z quadrupole isolation width). The acquisition strategy included a product ion triggered reisolation of the precursor ion if HexNAc oxonium ions (138.0545 and 204.0867 m/z) were detected amongst the top 20 fragment ions of the HCD-MS/MS spectrum. The reisolated precursor ions were subjected to both electron transfer dissociation with higher collisional dissociation supplemental activation (EThcD)-MS/MS and CID-MS/MS analysis (20, 21, 22). EThcD-MS/MS analysis was detected in the Orbitrap (60,000 resolution, 3e5 AGC, 250 ms injection time, calibrated charge dependent ETD reaction times [2+ 121 ms; 3+ 54 ms; 4+ 30 ms; 5+ 20 ms; 6+ 13 ms; 7+; 10 ms], 2 m/z quadrupole isolation width), and CID-MS/MS analysis was detected in the ion trap (rapid scan rate, 1e4 AGC, 70-ms injection time, 35% NCE, 2 m/z quadrupole isolation width). Data are available via ProteomeXchange with identifier PXD019372 (23). Username: reviewer24022@ebi.ac.uk and Password: KGAgOAkh

(Glyco)proteomic Data Analysis

The identification and quantification of peptides for proteomic analysis was performed with Proteome Discoverer (v2.1.0.801) using Sequest (24) against the rat UniProtKB database (November 2015; 31,095 entries). The precursor mass tolerance was set to 20 ppm with a maximum of two full trypsin miss cleavages while the CID-MS/MS fragment mass tolerance was set to 0.6 Da. The peptides were searched with oxidation of methionine set as variable modifications, whereas carbamidomethylation of cysteine and TMT of peptide N-terminus and lysine was set as a fixed modification. All data were searched as a single batch with peptide spectral match (PSM) and peptide false discovery rate (FDR) filtered to 1% using Percolator (25) and protein level FDR set to 1% using the Protein FDR Validator node. Quantification was performed using the Reporter Ion Quantifier node with integration set to 20 ppm and coisolation threshold set to 50%, and reporter ions were required in all channels. Peptides were grouped in each replicate based on unique sequence and unique modifications, and the median reporter ion areas calculated. Data were expressed as Log2 fold change to day 0 for each replicate and normalized to a median of 0. Total proteome from biological replicates across time points were batch effect corrected using an empirical Bayes model implemented in the sva R package (26). As recommended, the parametric shrinkage adjustment was applied to the data. The quality of the data after batch effect correction was assessed using the principal component analysis, and the corrected data were used for subsequent analysis. Significantly regulated glycopeptides were determined using ANOVA with permutation-based FDR set at 5% with Tukey’s post hoc test.

The identification and quantification of glycopeptides was performed with Proteome Discoverer (v2.1.0.801) using the Byonic node (v2.3.5) (27) against the rat UniProtKB database (November 2015; 31,095 entries). The precursor, HCD-MS/MS, and EThcD-MS/MS mass tolerance were set to 20 ppm with a maximum of two full trypsin miss cleavages. The peptides were searched with oxidation of methionine and N-glycan modification of asparagine (309 possible glycan compositions without sodium adducts available within Byonic) set as variable modifications, whereas carbamidomethylation of cysteine and TMT of peptide N-terminus and lysine were considered as fixed modifications. A precursor isotope off set was enabled (narrow) to account for incorrect precursor monoisotopic identification (±1.0 Da). All data were searched as a single batch with PSM FDR set to 1% using the PSM validator node, and a minimum Byonic score of 100 was applied, which has previously been shown to provide a good balance between accuracy and coverage (28). Only HCD-MS/MS spectra containing HexNAc oxonium ions (138.05–138.06 and 204.08–204.09 m/z) were annotated in the final list of glycopeptides as previously described (29). All identified glycopeptides were quantified based on HCD-MS/MS data using the Reporter Ions Quantifier node with integration set to 20 ppm and coisolation threshold set to 75%. Reporter ions were required in all channels. Peptides were grouped in each replicate based on unique sequence and unique modifications, and the median reporter ion areas were calculated. Data were expressed as Log2 fold change to day 0 for each replicate and normalized to a median of 0. Significantly regulated glycopeptides were determined using ANOVA with permutation-based FDR set at 5% with Tukey’s post hoc test.

Glycomic Sample Preparation

N-glycome profiling was performed essentially as described previously (30). Briefly, 10-μg protein extracts were dot-blotted onto polyvinylidene fluoride membranes and allowed to dry overnight. The membranes were stained briefly with direct blue 71 in 40% ethanol containing 10% acetic acid and washed with water. Immobilized proteins were excised and the membrane blocked with 1% polyvinyl pyrrolidine 4000 for 5 min followed by washing with water. N-glycans were released with 2.5 U PNGase F (Roche, Australia) for 16 h at 37 °C. Released glycans were collected, incubated with 100-mM ammonium acetate, pH 5.0, for 1 h at 23 °C and dried by vacuum centrifugation. N-glycans were reduced with 1 M NaBH4 in 50 mM KOH for 3 h at 50 °C and then desalted and enriched offline using AG 50W-X8 (Bio-Rad, Australia) strong cation exchange followed by porous graphitized carbon (PGC) solid-phase extraction microcolumns (Grace). For the determination of galactose linkages, aliquots of released glycans were incubated with combinations of 20 U broad-specificity sialidase (P0722S, α2-3,6,8,9), 8 U broad-specificity α-galactosidase (P0747S, α1-3,4,6), 10 U β-galactosidase (P0726S, β1-3), and 8 U β-galactosidase (P0746S, β1-3,4) (all from New England BioLabs). All reactions were performed in a final volume of 10 μl in 50-mM sodium acetate buffer, 5 mM CaCl2, pH 5.5, for 16 h at 37 °C.

Glycomic Mass Spectrometry

PGC-LC-ESI-MS/MS experiments were performed on a 3D ion trap using an Agilent 1100 capillary LC system (Agilent Technologies) interfaced with an Agilent 6330 LC-MSD 3D Trap XCT ultra. A PGC LC column (3 μm, 100 mm × 0.18 mm, Hypercarb, Thermo Scientific) was maintained at RT and at 50 °C. 10-mM ammonium bicarbonate aqueous solution (solvent A) and 10-mM ammonium bicarbonate aqueous solution with 45% MeCN (solvent B) were used as mobile phases at 2 μl/min with the following gradient: 0 min, 2% B; linear increase up to 35% B for 53 min; linear increase up to 100% B for 20 min; held constant for 5 min; and then equilibrated at 2% B for 5 min before the next injection—giving a total LC run time of 83 min. Glycans were analyzed using the following ESI-MS conditions: source voltage −3.2 kV, MS1 scan range 350 to 2200 m/z, 5 microscans, 0.13 m/z resolution (FWHM), 8e4 ion current control (ICC), and 200-ms accumulation time. Ion trap CID-MS/MS conditions were as follows: 0.13 m/z resolution (FWHM), 8e4 ICC, 200 ms accumulation time, 4 m/z isolation width and data-dependent acquisition of the three most abundant glycan precursors in each scan. The CID-MS/MS used ultrapure helium as the collision cell gas. Fragmentation amplitude was set to 1 V with Smart-Frag–enabled ramping from 30 to 200% of the fragmentation amplitude for CID-MS/MS with an activation time of 40 ms.

Glycomic Data Analysis

Data were analyzed as described previously (31). Briefly, a candidate list of glycans from relevant spectra of each sample was extracted using the ESI compass, v1.3, Bruker Daltonics Software (Bruker DALTONIK GmbH). The extracted monoisotopic precursor masses were searched against GlycoMod (http://www.expasy.ch/tools/glycomod) to identify putative monosaccharide compositions. Interpretation and validation of glycan identities, assisted by GlycoWorkBench v2.1 (available from https://code.google.com/archive/p/glycoworkbench), were based on the existence of A-, B-, C-, X-, Y-, and Z-product ions consistently found across the majority of the CID-MS/MS scans over the elution times of the respective precursor ions. Glycan annotation nomenclature was presented as described previously (32).

Relative glycan quantitation based on glycan precursor intensity was performed using Skyline. The relative abundance of each glycan was determined by the relative AUC peak areas of the monoisotopic m/z for relevant precursors divided by the AUC peak area sum of all glycans identified. The precursor mass analyzer was set to a quadrupole ion trap with resolution of 0.5 m/z, and precursor ions within the range 50 to 2000 m/z were included. Data are freely available on GlycoPOST (33) with the identifier GPST000079 (raw files fhttps://glycopost.glycosmos.org/preview/7030709905ed7e73ea0928, PIN: 5219) and Panorama Public (34) (skyline assays, https://panoramaweb.org/RatMuscleGlyco.url).

Cell and Tissue Lysis for Western Blot Analysis

Cells were washed twice with ice-cold PBS and lysis in 2% SDS by tip-probe sonication, whereas frozen muscle was sonicated directly in 2% SDS. Lysates were then centrifuged at 16,000g for 10 min at RT. Protein concentrations were determined via the BCA (Thermo Fisher), normalized, and diluted with Laemmli buffer. Samples were incubated at 65 °C, and 10 μg of protein was separated by SDS-PAGE using 10% gels. Proteins were transferred to polyvinylidene fluoride membrane and subjected to immunoblot analysis with anti-LGALS1 (#12936; Cell Signaling Technologies) and pan 14-3-3 loading control (#8312; Cell Signaling Technologies). Detection was achieved with HRP-conjugated donkey anti-rabbit secondary antibody (Jackson ImmunoResearch) and imaged by the Bio-Rad ChemiDoc Gel Imager system.

Immunofluorescence Microscopy

Cells were grown on Matrigel coverslips and analyzed 3 days after differentiation. Coverslips were washed twice with ice-cold PBS and fixed with 4% paraformaldehyde in PBS for 10 min at RT. The cells were washed three times with PBS and permeabilized in 0.2% triton X-100 in PBS for 10 min at RT. Three additional washes were performed and then blocked with 3% BSA in PBS for 1 h at RT. Cells were stained with MY-32 anti-myosin primary antibody (M4276; Sigma) at 1:400 dilution in PBS containing 1% BSA overnight at 4 °C. The coverslips were washed tree times with PBS and stained with goat anti-mouse secondary antibody Alexa Fluor 488 (A28175; Invitrogen) at 1:500 dilution in PBS containing 1% BSA for 1 h at RT. The cells were washed three times and stained with Hoechst 33342 (62249; Thermo Fisher Scientific) at 1 μg/ml in PBS for 5 min at RT. The coverslips were washed three times and mounted on glass slides and imaged on a Nikon C2 confocal microscope. Data was processed in ImageJ using Analyze Particle function to count nuclei (35).

Generation of Recombinant Adeno-Associated Viral Serotype-6 Vectors

DNA constructs expressing mouse Lgals1 were subcloned into AAV:multiple cloning site (MCS)-SV40pA plasmid (VectorBuilder). Recombinant adeno-associated viral vectors were generated by cotransfection of 10 μg of plasmids containing cDNA constructs with 20 μg of pDGM6 packaging plasmid into HEK293 cells (seeded 16 h prior at a density equivalent to 3.2–3.8 × 106 cells per P-100 tissues culture plate) using the calcium phosphate precipitate method to generate type-6 pseudotyped viral vectors as previously described (36). After 72 h, cells and culture media were collected, subjected to cell lysis by freeze-thawing, and clarified using a 0.22-μm filter (EMD Millipore). Virus was purified by heparin affinity chromatography on an Akta fast protein liquid chromatography (HiTrap, GE Healthcare), ultracentrifuged overnight, and resuspended in sterile Ringer’s solution. Vector concentrations were determined using qPCR (Applied Biosystems).

Intramuscular Injections of AAV6 Vectors

Intramuscular injections of AAV6 vectors were administered to either 2-day or 42-day-old C57BL/6J male mice maintained in a 12:12 h light:dark cycle at 22 °C with ad libitum access to normal chow diet and water. Mice were placed under general anesthesia (2% isoflurane in O2) and maintained on 37 °C heat pads. For injections into 2-day-old mice, the left hindlimb was injected with recombinant adeno-associated viral serotype-6 (rAAV6):MCS (1 × 1010 vector genomes per 10 μl) as a control, and the right contralateral limb was administered rAAV6:LGALS1 (1 × 1010 vector genomes per 10 μl). For injections in tibialis anterior (TA) muscle of 42-day-old mice, the left muscle was injected with rAAV6:MCS (1 × 109–1 × 1010 vector genomes per 30 μl) as a control, and the right contralateral limb was administered rAAV6:LGALS1 (1 × 109–1 × 1010 vector genomes per 30 μl). At experimental endpoint, mice were humanely killed by cervical dislocation and tissues rapidly excised. Ethical approval for mouse procedures was obtained from The University of Queensland (AEC approval #SBMS/101/13/NHMRC or #SBMS/AIBN/138/16/NHMRC/NHF), QIMR Berghofer Medical Research Institute (AEC approval #A18603M), and the University of Melbourne (AEC approval #1914940).

Experimental Design and Statistical Rationale

All presented data show a minimum of three independent biological replicates. For proteomic, glycoproteomic, and glycomic analysis of L6 cells, three biological replicates were performed. For glycomic analysis of skeletal muscle, three individual mice per timepoint were analyzed. ANOVA with permutation-based FDR set at 5% with Tukey’s post hoc test was used to determine significantly regulated proteins, glycans, or glycopeptides. For generation of clonal CRISPR cells, four independent WT and KO clones were analyzed, and three coverslips analyzed and compared with two-way ANOVA. For in vivo AAV delivery, five individual mice were analyzed and compared with paired t-test.

Results

Dynamic Remodeling of the Membrane-Associated Proteome During In Vitro Myogenesis

To identify membrane-associated proteins potentially important for myogenesis, we performed a quantitative proteomic analysis of L6 myoblasts during differentiation into multinucleated myotubes in 2D cell culture. Proteins were extracted daily over a 7-day time course of differentiation, and membrane-associated proteins were enriched with ultracentrifugation. Proteins were digested with trypsin followed by multiplexed stable isotope labeling with TMTs and analysis by 2D liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) (n = 3 biological replicates). A total of 5898 proteins groups were quantified with at least two peptides represented by 5801 unique genes in all biological replicates and time points (supplemental Table S1). More than half of the proteins contained membrane Gene Ontology cellular compartment annotation, and broad coverage of proteins associated within the Golgi apparatus and endoplasmic reticulum was also quantified (Fig. 1A). Principal component analysis revealed clustering of the biological replicates from each of the time points with a rapid and dramatic redistribution of the membrane-associated proteome during differentiation (Fig. 1B). Remarkably, more than 2700 proteins were differentially regulated in 7-day myotubes compared with undifferentiated myoblasts (±1.5-fold and q < 0.05 ANOVA with Tukey’s post hoc test) (Fig. 1C). Induction of well-characterized myogenic markers showed distinct profiles. For example, we observed rapid upregulation of various troponins (TNNI1/2, TNNC1, and TNNT3) within 1 to 2 days, whereas other proteins associated with the transcriptional regulation such as cysteine and glycine-rich protein-3 (CSRP3) and proteins associated with the contractile apparatus such as dystrophin (DMD) and myosin-binding proteins (MYOM1 and MYBPC1) progressively increased throughout the entire differentiation process (Fig. 1D). The kyoto encyclopedia of genes and genomes pathway enrichment revealed the regulation of several pathways associated with amino acid, fatty acid, and glucose metabolism in addition to other important signaling pathways and proteins associated with N-glycosylation (Fig. 1E). We next performed a clustering analysis of various anabolic and catabolic enzymes associated with glycosylation (Fig. 1F). We observed upregulation of various N-acetyl-β-hexosaminidases (HEXA/B/D) and bifunctional UDP-N-acetylglucosamine transferases (ALG13/14) along with the rapid increase of β-galactoside α-2,6-sialyltransferase-1 (ST6GAL1). Conversely, we observed decreased expression of mannosyl glycoprotein acetylglucosaminyl-transferases (MGAT1/2), five isoforms of β-galactoside α-2,3-sialyltransferases (ST3GAL1-5), and various galactosyltransferases (B3GALT6 and B4GALT1/4/5). Taken together, these data will be a rich resource for the study of membrane-associated proteins involved in myogenesis and suggest that the glycosylation machinery, and by extension protein glycosylation, is dynamically regulated.

Fig. 1

Dynamic membrane-associated proteomic analysis of myogenesis during L6 myotube formation.A, gene Ontology cellular compartments analysis. B, principal component analysis. C, the number of regulated protein (ANOVA with permutation-based FDR set at 5% with Tukey’s post hoc test). D, quantification of myogenic markers. E, KEGG pathway enrichment analysis. F, temporal clustering of glycosylation-associated proteins. FDR, false discovery rate.

Integrated and Dynamic N-Glycomic and N-Glycoproteomic Analysis Identifies Changes in Digalactosylation to α-2,3-Sialylation on Cell Adhesion and ECM Proteins During Myogenesis

We next investigated the regulation of protein glycosylation by first performing a glycomic analysis of released N-glycans with PGC LC-MS/MS over the 7-day time course of L6 myoblast differentiation (n = 3 biological replicates). We identified 45 unique monosaccharide compositions that were confirmed by MS/MS and quantified in all biological replicates and time points corresponding to a total of 68 different glycan structures (supplemental Table S2). The N-glycome was made up of paucimannosidic- (Man1–3GlcNAc2Fuc0–1), oligomannosidic- (Man5–9GlcNAc2), hybrid-type glycans capped by a single α-2,3- or α-2,6-NeuAc or di-Gal moiety with or without core fucosylation, and finally complex-type glycans with monoantennae, diantennae, or triantennae with one to three terminal α-2,3-NeuAc, α-2,6-NeuAc, or di-Gal with or without core fucosylation and minor amounts of bisecting β-1,4-GlcNAc (supplemental Fig. S1). Analysis of glycans treated with broad-specific sialidase only, a combination of broad-specific sialidase, and β-1,4-galactosidase or sialidase and α-galactosidase revealed that the di-Gal comprised a galactose residue α-linked to β-1,4-galactose (supplemental Fig. S2). Importantly, none of the di-Gal moieties were found to be capped by sialic acids, suggesting that the α-linked Gal is not a substrate for sialyltransferases in these rat-derived cells. A total of 34 glycans were differentially regulated during differentiation (q < 0.05 ANOVA with Tukey’s post hoc test) (Fig. 2A). There was progressive decrease in N-glycans containing α-2,3-NeuAc (Fig. 2B) along with a rapid increase in N-glycans containing α-2,6-NeuAc (Fig. 2C) (supplemental Fig. S3, A–B). The temporal glycan expression matched the expression of sialyltransferases catalyzing these glycoforms confirming that myoblast-to-myotube differentiation is associated with sialic acid linkage switching (Fig. 2D). We also observed a progressive decrease in N-glycans containing terminal di-Gal epitopes (Fig. 2E) with a concomitant decrease in the expression of galactosyltransferases during differentiation (Fig. 2F). This included several B4GALT enzymes catalyzing the formation of LacNAc moieties that are required substrates for di-Gal extensions and sialic acid capping. Consistent with the downregulation of these enzymes, glycans containing LacNAc moieties were also decreased accompanied with an increase in nongalactosylated β-linked GlcNAc termini (supplemental Fig. S3, C–D). Furthermore, glycans containing α-1,3-Gal were dramatically downregulated; however, our proteomic analysis did not quantify N-acetyllactosaminide alpha-1,3-galactosyltransferase (GGTA1), which catalyzes the transfer of α-1,3-Gal to β-1,4-linked Gal of LacNAc moieties (supplemental Fig. S3E). Therefore, we investigated previously published transcriptomics data during myogenesis in mouse C2C12 cells (37) revealing a rapid downregulation Ggta1 in the immediate hours of differentiation (supplemental Fig. S3F). As expected, the decreased available LacNAc substrates resulted in reduced total amount of capping features in particular the α-2,3-sialylation and α-1,3-galactosylation capping efficiency (supplemental Fig. S3G). Interestingly, paucimannosylation rapidly increased during differentiation along with the expression of the catalyzing enzyme N-acetyl-β-hexosaminidase-A (HEXA) while HEXB only subtly increased (Fig. 2G). Grouping the glycan types revealed both oligomannosylation and paucimannosylation are elevated, whereas complex- and hybrid-type glycans are weakly reduced in advanced myotube formation (supplemental Fig. S3, H–I)

Fig. 2

Dynamic A, the heat map showing the average fold change of N-glycans relative to predifferentiation at day 0. B, quantification of α-2,3-NeuAc–containing N-glycans. C, quantification of α-2,6-NeuAc–containing N-glycans. D, quantification of sialyltransferases. E, quantification of N-glycans containing terminal di-Gal. F, quantification of galactosyltranferases. G, quantification of paucimannosidic N-glycans and N-acetyl-β-hexosaminidase A/B (HEXA/B). ∗q < 0.05 ANOVA with permutation-based FDR. FDR, false discovery rate.

To identify the glycoproteins modified during differentiation, we next performed quantitative N-glycoproteomics. Membrane-associated proteins were isolated over a 7-day time course of differentiation (n = 3) and digested with trypsin followed by multiplexed stable isotope labeling with TMT. Glycopeptides were enriched with zwitterionic-hydrophilic interaction liquid chromatography and analyzed by LC-MS/MS using data-dependent fragmentation decisions of HCD and product-dependent oxonium ion–triggered electron transfer dissociation with HCD supplemental activation (pdEThcD). A total of 2751 unique glycopeptides were identified with 2238 identified by HCD and 1518 identified by pdEThcD (1005 with both fragmentation types) (supplemental Table S3). A change in the abundance of a specific glycopeptide can arise from either alteration in the (i) site occupancy, (ii) type of glycosylation at a given site or, (iii) a change in the expression of the glycoprotein itself. Figure 3, A–B displays changes in the abundance of glycopeptides plotted against changes in their total protein levels at day 1 and day 7 of differentiation relative to day 0, respectively. The green box highlights glycopeptides from integrin α-3 (ITGA3), which reveals increased protein expression on the y-axis and a distribution of glycoforms on the x-axis. The annotated HCD-MS/MS spectrum indicated with the arrow is shown in Figure 3C, identifying a glycopeptide containing Asn-605 in the N-linked glycosylation motif (NxS/T/C) conjugated with HexNAc(4)Hex(6)Fuc(1) corresponding to a glycan composition with biantennary terminal di-Gal moieties. These data are consistent with the global glycomic data showing downregulation of di-Gal–containing N-glycans and the observed downregulation of galactosyltransferases in the proteomics data. Figure 3D plots the relative changes of this glycopeptide and the total protein levels of ITGA3 across the differentiation time course. ITGA3 was more than 6-fold upregulated after 7 days of differentiation at the protein level, whereas the indicated glycopeptide was almost 2-fold downregulated. Normalizing these data revealed this specific glycosylation was reduced by more than 10-fold. A total of 931 unique glycopeptides were quantified in at least two biological replicates in all time points, and 91 were regulated across the differentiation time course (q < 0.05 ANOVA with Tukey’s post hoc test) (supplemental Table S4). Several additional cell surface and extracellular proteins contained decreased di-Gal–containing N-glycans such as other integrins (ITGA7 and ITGB1), receptors (CD63, ATRN, NOMO1, IL6ST), and ECM-associated proteins (FN1 and EMILIN1). Unlike the glycomic analysis of released N-glycans, our analysis of intact glycopeptides is not able to accurately differentiate glycopeptides containing different sialic acid linkages (α-2,3-NeuAc versus α-2,6-NeuAc). Despite this, we observed the regulation of several glycopeptides modified with sialic acid–containing N-glycans, suggesting changes in the amount of sialylation. These include several proteins also containing digalactosylation and additional proteins such as insulin-like growth factor-2 receptor. Many of the proteins with changes in glycosylation including integrins (38) and insulin-like growth factor-2 receptor (39) have been previously implicated in the regulation of myogenesis, and here, we provide support for and expand on these findings by showing precise regulation of di-Gal and sialylation, which may play important roles in either protein function or interaction with lectins during the various stages of myogenesis including differentiation, elongation, migration, cell adhesion, membrane alignment, or cell fusion.

Fig. 3

Dynamic Quantification of glycopeptides plotted against changes in their total protein levels at day 1 (A) and day 7 (B) of differentiation relative to day 0. The green box highlights glycopeptides from integrin α-3 (ITGA3). The arrow indicates annotated HCD-MS/MS spectrum in panel C identifying a glycopeptide containing Asn605 in the N-glycosylation motif (NxS/T/C) with the glycan HexNAc(4)Hex(6)Fuc(1) corresponding to a glycan composition with biantennary terminal di-Gal. D, Relative changes of the indicated glycopeptide and the total protein levels of ITGA3 across the differentiation time course. HCD, higher collisional dissociation; MS/MS, tandem mass spectrometry.

Galectin-1 Is Required for Myoblast Differentiation In Vitro

Many glycoproteins form glycan–protein interactions with glycan-binding proteins (40). The regulation of galactosylation during myogenesis is interesting as these glycans can be recognized by various galectins, forming interactions that play important roles in several biological processes (41). Mice lacking various galectins display alterations in smooth muscle, cardiac muscle, and skeletal muscle development or regeneration/remodeling after injury (10, 42, 43). However, galectins are ubiquitously expressed in many cell types and their roles in the various stages of myogenesis and muscle development are incompletely understood. Our proteomic analysis identified time-dependent regulation of several galectins during the differentiation of L6 myoblasts in vitro including the upregulation of galectin-1 (LGALS1) and downregulation of galectin-3 (LGALS3) (Fig. 4A). We also observed an initial subtle increase in galectin-8 (LGALS8), which declined after 2 days whereas galectin-3 binding protein (LGALS3BP) initially decreased and then increased after 3 days of differentiation. Here, we focus on galectin-1 and confirm time-dependent increases in protein expression during differentiation of L6 myoblasts by Western blot analysis (Fig. 4B). We next genetically ablated Lgals1 from L6 myoblasts using paired gRNAs and CRISPR/Cas9(D10A) (Fig. 4C). Clonal cells were isolated, and four independent KO lines were validated by Western blot analysis (Fig. 4D). We also isolated four WT clonal cells from the original population transfected with the gRNAs and CRISPR/Cas9(D10A) that did not undergo editing to serve as controls. After 3 days of differentiation, WT myoblasts formed elongated multinucleated and myosin-positive myotubes, whereas KO cells mostly remained mononucleated with a small number of multinucleated cells which stained for myosin (Fig. 4E). Quantification of the differentiation index showed a marked decrease in the percentage of nuclei in myosin-positive cells, indicating loss of Lgals1 significantly reduced myoblast differentiation (Fig. 4F). Our results support previous in vivo data showing mice lacking Lgals1 have reduced myofiber diameter and delayed muscle regeneration (10), and we propose impaired myotube formation is likely mediated by an early defect in the myogenic program and failure to initiate differentiation.

Fig. 4

Galectin-1 is upregulated and required for L6 myotube formation.A, quantification of galectins during myogenesis. B, western blot analysis of galectin-1 during myogenesis. C, CRISPR/Cas9(D10A) paired gRNA strategy to KO Lgals1. D, western blot analysis of galectin-1 and 14-3-3 loading control in WT and KO L6 myoblasts. E, immunofluorescence microscopy of WT and KO cells at day 3 of differentiation (blue = Hoechst; green = myosin; white scale bar = 50 μm). F, quantification of the differentiation index expressed as the percentage of nuclei in myosin-positive cells. ∗q < 0.05 ANOVA with permutation-based FDR. #p < 0.05 two-way ANOVA. FDR, false discovery rate.

Galectin-1 and N-glycosylation Are Regulated During Mouse Muscle Development In Vivo

We next investigated the regulation of protein N-glycosylation and galectin-1 expression during postnatal mouse skeletal muscle development. Male mice were sacrificed at 1, 7, 14, and 21 days after birth and skeletal muscle collected from lower hindlimbs (n = 3 biological replicates). Western blot analysis revealed an increase in galectin-1 throughout development consistent with the in vitro myogenesis data (Fig. 5A). Next, released N-glycans were quantified in total protein lysates throughout postnatal development with PGC-LC-MS/MS. We identified 54 unique monosaccharide compositions that were confirmed by MS/MS and quantified in all biological replicates and developmental time points corresponding to a total of 104 different glycan structures (supplemental Table S5). The N-glycome was made up of paucimannosidic-, oligomannosidic-, and hybrid-type glycans capped with a single α-2,3- or α-2,6-NeuAc or di-Gal moiety with or without core fucosylation, and also complex-type glycans displaying monoantennae, diantennae, or triantennae with one or two terminal α-2,3-NeuAc, α-2,6-NeuAc, α-2,3-NeuGc, α-2,6-NeuGc, or di-Gal with or without core fucosylation (supplemental Fig. S4). There was a progressive decrease in N-glycans containing α-2,3-NeuAc that was significant 21 days after birth (Fig. 5C). Furthermore, a rapid reduction in α-2,6-NeuAc was also observed 7 days after birth, but this gradually increased throughout the next 2 weeks of development. The seemingly inverse relationship between these two sialic acid linkages is consistent with the in vitro myogenesis model, and the biological consequences of these changes on the myogenic program warrants further investigation. Also in line with the in vitro data, the N-glycans containing terminal di-Gal were reduced at day 7, 14, and 21 days after birth, paucimannosylation was significantly increased 14 days after birth, and by 21 days, the levels were >2.5-fold greater than newborn mice. Taken together, our data suggest that N-glycan remodeling is a strong feature associated with in vitro and in vivo rodent models of myogenesis.

Fig. 5

Quantification of LGALS1 and released A, postnatal time points. B, western blot analysis of galecin-1 and 14-3-3 loading control in skeletal muscle development. C, quantification of selected N-glycan features during skeletal muscle development. ∗q < 0.05 ANOVA with permutation-based FDR.

In Vivo Galectin-1 Gain-Of-Function Increases Muscle Mass

Administration of recombinant galectin-1 has recently been shown to improve muscle function in a mouse model of muscular dystrophy (14). Therefore, we next investigated the ability of galectin-1 to promote postnatal muscle development using an in vivo gene therapy gain-of-function mouse model. Hindlimbs of newborn mice were injected with rAAV6 using a paired experimental design where the left legs were treated with viral vectors encoding an empty MCS to serve as a control and the right legs received constructs expressing Lgals1 to transduce the entire lower limb musculature (n = 5 biological replicates) (Fig. 6A). An initial pilot dose response was performed in the TA muscle of 42-day-old mice revealing injections of ∼1 × 1010 vector genomes (vg) produced ∼3-fold average overexpression by Western blot analysis (data not shown). Two days after birth, mice were injected with rAAV6, and after 42 days, the animals were sacrificed and muscles dissected. Western blot analysis of gastrocnemius muscle revealed a variable but significant ∼2.2-fold average overexpression of galectin-1 (Fig. 6, B–C). The TA muscles injected with rAAV6:LGALS1 all increased in mass with an average of 5.3%, but there was considerable variation likely because of variation in expression (Fig. 6D). Our data are the first to reveal that enhanced muscle is mass associated with galecin-1 expression in early skeletal muscle development.

Fig. 6

A, strategy for in vivo gain-of-function mouse model using recombinant adeno-associated virus serotype-6 (rAAV6) to express an empty multiple cloning site (MCS) control or LGALS1. B, western blot analysis of galectin-1 and 14-3-3 loading control of skeletal muscle from left leg (L) injected with rAAV6:MCS control and right leg (R) injected with rAAV6:LGALS1. C, quantification of western blot. D, mass of tibialis anterior (TA) muscle normalized to body weight injected with rAAV6:MCS control or rAAV6:LGALS1. ∗p < 0.05 paired t-test.

Discussion

The integration of proteomics, glycoproteomics, and glycomics is proving to be a valuable approach for holistic characterization of glycoproteome regulation (44). No single technology can comprehensively unravel the complexities of glycosylation and each platform offers unique insights, proteomics enables quantification of glycoproteins, lectins, and glycosyltransferase/glycosidase abundance, glycoproteomics enables quantification of site-specific changes in glycan composition, and glycomics enables detailed structural characterization of the glycome. Here, we integrate these technologies to study dynamic regulation of protein N-glycosylation, enabling us to define changes in the abundance of galectins, galactosyltransferases, sialyltransferases, and N-acetyl-β-hexosaminidases that correlate with specific glycan structures on functionally important receptors and adhesion molecules. These integrated data substantially improve our confidence in defining the mechanisms of glycosylation changes during myogenesis. Specifically, we show a switch in terminal di-galactosylation and sialylation linkages mediated by changes in glycosyltransferase/sialyltransferase expression. The is interesting because these enzymes are colocalized in the Golgi apparatus and will compete for the generation of galacto- or sialo-epitopes presented for reactivity. It is important to note that digalactosylation is not a common feature in the human glycome, meaning the use of these glycoepitopes in myogenesis may be different to rodent models analyzed in the present study. Our data also revealed an increase in paucimannosylation during terminal differentiation, which is surprising as this modification has previously been associated with proliferation in cancer (45) and neuronal progenitor cells (46) and is downregulated during macrophage differentiation (47).

Recent advances in glycoproteomics are dramatically increasing our ability to identify thousands of glycopeptides (48, 49, 50, 51, 52, 53, 54, 55), and future studies are likely to see an expansion of methods to analyze large sample cohorts and further enhance identification confidence or structural characterization (56, 57, 58). One missing piece of information in our study is the assessment of glycan heterogeneity and site-specific occupancy on intact protein species. New methods focused on top-down glycoproteomics to characterize intact glycoproteins will further enhance our ability to achieve this level of detail on a proteome-wide scale (59). Our study focused on N-glycosylation, but it is important to note that other types of protein glycosylation such as O-, GPI-, or C-linked may also be important during myogenesis. Furthermore, additional glycoconjugates such as glycosphingolipids, glycosaminoglycans, and proteoglycans are likely to play important roles in the development of skeletal muscle as mutations in the catalytic enzymes responsible for their biosynthesis cause developmental neurological defects (60), and pharmacological inhibition modulates adult skeletal muscle metabolic function (61).

Our proteomic analysis quantified an increase in galecin-1 abundance and a reciprocal downregulation of galectin-3. This is interesting because both galectins have similar functions such as the regulation of Ras-mediated signaling in cancer (62). Reciprocal regulation of galectin-1 and galectin-3 has also been observed in sera from patients with rheumatoid arthritis undergoing various treatments (63), suggesting possible compensatory mechanisms. Galectins form a family of 15 members that bind to galactose via conserved carbohydrate-binding domains (64). Galectin-1 is a prototypical galectin containing a single carbohydrate-binding domain that form dimers, whereas galectin-3 contains an extended N-terminal domain with the ability to form pentamers (65). Administration of recombinant galectin-1 to the mdx mouse model of muscle dystrophy increased utrophin and integrin expression and improved skeletal muscle function (14). Furthermore, enhanced muscle function has also been observed with exogenous galectin-3 expression in a similar mdx mouse model (66). Our data further support these observations, revealing overexpression of galectin-1 increases muscle mass during postnatal muscle development in mice. Galectins have been shown to bind galecto-epitopes on various cell surface receptors including integrins to promote signaling (64). Surprisingly, we show that various integrins display decreased di-galactosylation during myogenesis, suggesting reduced galectin binding ability. This is somewhat contradictory; however, the presence of α-2,3-NeuAc versus α-2,6-NeuAc can further influence the binding of various galectins to galactose epitopes. For example, α-2,6-NeuAc inhibits galectin-1 binding to galactose but has no effect on galectin-3 binding while both galectin-1 and -3 can bind galactose containing α-2,3-NeuAc-linked glycans (67). The differential binding of various galectins to glycans containing different sialic acid linkage isomers further supports a complex interplay between these glycoforms in the various stages of myogenesis. Further experiments are required to precisely quantify the glycan–protein interactions during myogenesis and muscle growth. Galectins are secreted via nonclassical secretory mechanisms and may also play a role in endocytosis (68, 69). Therefore, galectins have been primarily studied in the context of extracellular lattice formation; however, there is growing evidence to suggest galectins may regulate intracellular signaling via lysosomal trafficking (70). For example, LGALS8 has recently been shown to inhibit mTOR signaling at the lysosome (71). This is interesting because the activation of mTOR has well established roles in skeletal muscle growth (72). Our data revealed a decrease in LGALS8 expression during myogenesis, suggesting a relief in mTOR inhibition required for cell growth and further experiments are required to investigate this in skeletal muscle. We anticipate our integrated analysis of proteomics, glycomics, and glycoproteomics will be a valuable resource to further our understanding of the role of glycosylation on muscle development.

Data Availability

Supplementary data are available online. Proteomic and glycoproteomic data are available via ProteomeXchange with identifier PXD021006 (23). Glycomic data are freely available on GlycoPOST (32) with the identifier GPST000079 and Panorama Public (33) (skyline assays, https://panoramaweb.org/RatMuscleGlyco.url).

Conflict of interest

The authors declare no competing interests.