Weidong Li
Kristof De Schutter
Els J. M. Van Damme
Guy Smagghe- 1Laboratory of Agrozoology, Department of Plants and Crops, Faculty of Bioscience Engineering, Ghent University, Ghent, Belgium
- 2Laboratory of Biochemistry and Glycobiology, Department of Biotechnology, Faculty of Bioscience Engineering, Ghent University, Ghent, Belgium
Protein glycosylation is one of the most common and most important post-translational modifications. Despite the growing knowledge on N-glycosylation, the research on O-glycosylation is lagging behind. This study investigates the importance of O-glycosylation in the post-embryonic development of insects using the red flour beetle, Tribolium castaneum, as a model. We identified 28 O-glycosylation-related genes (OGRGs) in the genome of the red flour beetle. 14 OGRGs were selected for functional analysis based on their involvement in the initial attachment of the carbohydrate in the different O-glycosylation pathways or the further elongation of the most abundant O-glycans and, in addition, showing severe RNAi-induced phenotypes in Drosophila melanogaster. The expression profile of these OGRGs was mapped throughout the developmental stages of the insect and in the different tissues of the pupa and adult. Subsequently, these genes were silenced using RNA interference (RNAi) to analyze their role in development. A broad spectrum of phenotypes was observed: from subtle effects and disrupted wing formation when silencing the genes involved in O-mannosylation, to blockage of pupation and high mortality after silencing of the genes involved in O-GalNAc and core 1 O-glycan (O-GalNAc-Gal) synthesis. RNAi experiments were also performed to assess the effects of blocking multiple pathways of O-glycosylation. However, the observed phenotypes induced by multiple RNAi were similar to those of the single gene RNAi experiments. The silencing of OGRGs often resulted in high mortality and wing phenotypes, indicating the importance of O-glycosylation for the survival of the insect and the formation of wings during the post-embryonic development of T. castaneum.
Introduction
The decoration of proteins with carbohydrates is one of the most important post-translational modifications. One of the major types is O-glycosylation, in which a sugar is transferred to the oxygen atom of a hydroxyl group of mainly serine and threonine, but in some cases hydroxyproline (plant), hydroxylysine (collagens) or tyrosine (prokaryotes) residues can also be O-glycosylated (Spiro et al., 1971; Kieliszewski, 2001; Zarschler et al., 2010). O-glycosylation is found in all kingdoms of life, including mammals, insects, fungi, archaea, plants and bacteria (Ten Hagen, 2015; Eichler and Koomey, 2017). O-glycosylation plays a crucial role in many biological processes, including, stress response and immunization through cell-to-cell interaction and protein-to-protein interaction, and influences protein stability, folding, localization and secretion (Fukuda, 2002; Schjoldager et al., 2011; Ashwood et al., 2017; Chaffey et al., 2017a, b; Martinez et al., 2017). Previous studies have shown that unusual O-glycosylation is associated with a variety of diseases and cancers (Yang et al., 2017). In contrast to N-glycosylation, where an oligosaccharide (glycan precursor) is transferred to the nascent protein in the endoplasmic reticulum (ER), O-glycans are synthesized by the stepwise addition of monosaccharides in the ER, Golgi or cytosol (Spiro, 2002; De Pourcq et al., 2010; Dell et al., 2010; Figure 1). In insects, O-glycans consist mainly of eight monosaccharides –N-acetylgalactosamine (GalNAc), galactose (Gal), N-acetylglucosamine (GlcNAc), fucose (Fuc), mannose (Man), glucose (Glc), glucuronic acid (GlcA) and xylose (Xyl)– which can be assembled into oligosaccharides or polysaccharides. The attached residue in O-glycosylation can be GalNAc, GlcNAc, Fuc, Man, Glc, Xyl (Haltom and Jafar-Nejad, 2015; Walski et al., 2017), while the attached monosaccharide residue in N-glycosylation is always a GlcNAc. This heterogeneity in the attachment residue makes O-glycosylation more diverse and complex than N-glycosylation. This complexity makes the study of O-glycosylation in insects challenging. So far there are only a couple of studies focused on O-glycosylation in insects such as Drosophila melanogaster (Aoki et al., 2008), mosquito (Kurz et al., 2015) and silkworm (Soya et al., 2016). These studies investigated the O-glycan profile during specific developmental stages (embryos, larvae) or in specific tissues (nervous system). In Drosophila embryos, the O-GalNAc-Gal (core 1 structure or T antigen) is the most prevalent O-glycan structure representing 55% of the total O-glycan pool; O-GalNAc (Tn antigen) and O-GlcNAc account for 19% of the total profile; O-GalNAc-Gal-GlcA and O-Fuc(GlcNAc)-GlcA occupy about 10% each; while O-GalNAc(Gal)-GlcA, O-Man/O-Glc, and O-GalNAc(HexNAc)-GlcA each account for around 1% of the total O-glycan profile (HexNAc, N-acetylhexosamine, including GalNAc and GlcNAc). Other O-glycans are much less abundant (Aoki et al., 2008). Although the O-glycan profiles of a specific tissue or at a specific stage yield valuable information, they do not provide a general understanding of O-glycosylation in the insect.
Figure 1. O-glycan biosynthesis pathway. O-glycans are attached to serine and threonine. The known enzymes responsible for each attachment are indicated beside the arrow. Ofut1 is responsible for transferring monosaccharide specifically to EGF repeats; Ofut2 specifically transfers Fuc to TSRs. EOGT and Rumi specifically transfer sugars to EGF repeats. POMT1 and POMT2 cooperate in the transfer of Man. pgants, N-acetylgalactosaminyltransferases; POMTs, protein O-mannosyltransferases; Ofuts, O-fucosyltransferases; OGTs, O-GlcNAc transferases; C1GalT, core 1 β1,3-galactosyltransferase; fng, Fuc-specific O-GlcNAc transferase; Rumi, protein O-glucosyltransferase; GalNAc, N-acetylgalactosamine; Gal, galactose; GlcNAc, N-acetylglucosamine; Fuc, fucose; Man, mannose; Glc, glucose; GlcA, glucuronic acid; Xyl, xylose.
O-glycosylation involves several different pathways, leading to the attachment of O-GalNAc, O-GlcNAc, O-Fuc, O-Man, O-Glc type O-glycans and the OGRGs encode the enzymes responsible for the synthesis of these diverse O-glycan structures (Figure 1). In Drosophila, 12 polypeptide N-acetylgalactosaminyltransferases (pgants) are responsible for transferring GalNAc to protein substrates, and 8 core 1 β1,3-galactosyltransferases (C1GalTs) are responsible for transferring Gal to the O-linked GalNAc. These pgant genes are characterized by a distinct spatial and temporal expression pattern during development (Tian and Hagen, 2006), for example, pgant3 is located in the wing disk and is responsible for wing adhesion, loss of pgant3 causes wing blistering in Drosophila (Zhang et al., 2008, 2010). Pgant5 and pgant35A have an influence on gut acidification and the respiratory system, respectively (Tran et al., 2012) and mutations or silencing of these genes causes complete lethality in the fruit fly (Yamamoto-Hino et al., 2015). In D. melanogaster, many proteins were shown to carry mucin-type O-glycans (Schwientek et al., 2007; Zhang et al., 2008; Walski et al., 2017) and knockdown of these proteins resulted in lethality (Mummery-Widmer et al., 2009; Fraichard et al., 2010; Kong et al., 2010; Maynard et al., 2010; Neely et al., 2010; Nagarkar-Jaiswal et al., 2015). The O-GlcNAc transferases OGT is responsible for transferring GlcNAc to intracellular target proteins (Ogawa et al., 2015), while EOGT specifically glycosylates Notch and its ligands (Li et al., 2020), and interacts with Notch signaling, required for many developmental processes (Müller et al., 2013). O-mannosylation of Dystroglycan (a major non-integrin extracellular matrix receptor) occurs through the cooperation of the protein O-mannosyltransferases (POMT) POMT1 and POMT2 (Ichimiya et al., 2004; Nakamura et al., 2010). Since this protein is required for proper muscle formation, its mutation caused major mortality and locomotor behavior defective phenotypes in Drosophila (Shcherbata et al., 2007; Yatsenko et al., 2014). In Drosophila, Ofut1 is an O-fucosyltransferase (Ofut) specific for EGF repeats while Ofut2 transfers Fuc to thrombospondin type-1 repeats (TSRs) (Luo et al., 2006). The Fucose-specific O-GlcNAc transferase (Fng) encodes a glycosyltransferase responsible for transferring GlcNAc specifically to the O-fucose of EGF repeats (Matsumoto et al., 2015), and mutation of the gene in fruit fly is lethal. O-Glc structures represent less than 1% of the identified O-glycans in Drosophila embryos. However, in silkworms, monosaccharide profiling revealed that Glc is the most abundant monosaccharide in the nervous system during development (Aoki et al., 2008; Soya et al., 2016). Insects encode one glucosyltransferase, Rumi, which catalyzes the transfer of Glc residues to EGF repeats of secreted proteins and plays a redundant role with Ofut1 in the trafficking of Notch out of the ER and in the activation of Notch signaling (Ishio et al., 2015; Matsumoto et al., 2016).
In D. melanogaster RNAi of these OGRGs indicated that O-glycosylation was essential for viability and movement, and plays vital roles in the respiratory, immune and nervous system (Ichimiya et al., 2004; Okajima et al., 2005; Tian and Hagen, 2006; Zhang et al., 2008; Sinclair et al., 2009; Sanders et al., 2010; Zhang et al., 2010; Park et al., 2011; Haltom and Jafar-Nejad, 2015; Ten Hagen, 2015). However, the specific function of O-glycosylation in the development of insects remains unclear.
This project addressed the importance of O-glycosylation in the post-embryonic development of the holometabolic model insect, Tribolium castaneum or red flour beetle. Therefore, we first identified the genes encoding enzymes involved in O-glycosylation in the genome of T. castaneum and analyzed the changes in transcription levels of the OGRGs involved in the synthesis of the main O-glycan structures and their precursor monosaccharides throughout the different life stages and body parts of the red flour beetle. These genes include the UDP-Galactose 4-epimerase (GALE, an enzyme catalyzing the interconversion of UDP-Gal and UDP-Glc in the final step of the Leloir pathway of galactose metabolism, and also the interconversion of UDP-GalNAc and UDP-GlcNAc), and the OGRGs involved in the first two steps of the different O-glycosylation pathways (pgants, POMTs, Ofuts, OGTs, C1GalT, fng and Rumi). In addition, the role of O-glycosylation in development was studied using RNAi of the OGRGs by injecting dsRNA in the larval stage. We discovered that most O-glycan transferases reached a peak of gene expression in the pupal stage, and a decrease in OGRG transcription levels resulted in mortality as well as abnormal wings and antennae. Therefore, we hypothesized that O-glycosylation plays an important role in metamorphosis and wing development of T. castaneum.
Materials and Methods
Insect Culture
The wild type strain T. castaneum GA-1 was used in all experiments. Insects were reared on whole wheat flour containing 5% yeast powder as daily diet under standard conditions (30°C and humidity 60% in darkness) as previously described (Walski et al., 2016).
Identification of OGRGs and Phylogenetic Analyses
Protein sequences encoding the OGRGs from D. melanogaster were used to perform BLAST searches against the T. castaneum genome through Tribolium BLAST1 (E-value < 0.0001). Protein sequences were downloaded from2 and screened for the presence of conserved protein domains using interproscan 5 (Jones et al., 2014). The program was downloaded3 and locally installed. Only sequences with the required domains were retained. Phylogenetic trees were constructed including known sequences from D. melanogaster, Leptinotarsa decemlineata, Aedes aegypti, Bombus terrestris, Bombyx mori, Nilaparvata lugens, Acyrthosyphon pisum and Homo sampiens in MEGA X using Maximum Likelihood method. The ML substitution matrix was analyzed and the best fitting model was selected based on the BIC value. The phylogeny was tested with 1000 bootstrap replicates. Protein sequences from all species were identified in BLAST searches against the respective genomes available from NCBI4 (E-value < 0.0001). All identified OGRGs of T. castaneum are listed in Table 1. All protein sequences used in the phylogenetic analysis were collected in Supplementary File 1.
Table 1. Overview of putative genes involved in O-glycosylation of T. castaneum, included in this study.
RNA Extraction and cDNA Synthesis
Samples were collected from egg to adult stage of T. castaneum. First, eggs that were laid within a period of 24 h were collected using an 800 Mic sieve. After hatching, insects were collected every three days for second to sixth instar and prepupa. Then samples of pupa and late pupa were collected at 20 and 22 days post hatching. Insects collected at 24 and 27 days post hatching represent the early adult and mature adult samples, respectively. Ten milligrams of fresh insects were taken as a sample for RNA extraction. Additionally, head, thorax, abdomen, wings (forewings and hindwings) and legs were dissected under a microscope, and then separately collected from 10 pupae or 20 adults. Total RNA was extracted using the RNeasy mini kit (Qiagen) according to the manufacturer’s instructions. The concentration of total RNA was measured by Nanodrop ND-1000 (Thermo Scientific) and the quantity was checked on a 1.5% agarose gel after removing genomic DNA using the TURBO DNA-freeTM Kit (Ambion/Life Technologies). cDNA synthesis was conducted following the manufacturer’s protocol using SuperScriptTM II Reverse Transcriptase (Invitrogen/Life Technologies). One microgram of total RNA was used in cDNA synthesis with an Oligo(dT)12 primer. Finally, the cDNA solution was stored at −20°C for subsequent experiments.
Primers
All the primers were designed by Primer3Plus5. The sequences of the primers for the reference genes were obtained from publications (Lord et al., 2010; Toutges et al., 2010; Vuerinckx et al., 2011; Feliciello et al., 2015). Primer specificity was checked using PrimerBlast in NCBI and on gel. All primer sequences are listed in the Supplementary Tables 1, 2.
Reverse Transcription Quantitative PCR
Reverse transcription quantitative PCR (RT-qPCR) was performed using a CFX96TM Real-Time PCR Detection system (BioRad) with GoTaq® qPCR Master Mix (Promega). A standard curve with a serial dilution of cDNA was made for each primer pair to determine the amplification efficiency, followed by an analysis of the melting curve to ensure specificity. Each qPCR reaction was performed in duplicate, which contained 10 μl of GoTaq® qPCR Master Mix, 0.4 μl of 10 μM forward primer, 0.4 μl of 10 μM reverse primer (Invitrogen), 8.2 μl of water and 1 μl of cDNA. All reactions were performed with 95°C for 5 min, followed by 40 cycles of 95°C for 15 s, 60°C for 60 s. No-template controls and no-reverse transcriptase controls were performed to exclude contamination. Ubiquitin and RPS18 were chosen as reference genes based on previous studies and our optimization with qBase + software (Biogazelle). All the gene expressions were analyzed by qBase + software (Biogazelle).
Larva RNAi Experiments
DsRNAs were synthesized following the instructions of the MEGAscript® RNAi Kit (Life Technologies) with primers containing the T7 promoter sequence. The concentration of dsRNA was measured by Nanodrop ND-1000 (Thermo Scientific). The quality was checked on a 1.5% agarose gel. Fourth-instar larvae (with an individual fresh weight of 1.2 ± 0.2 mg) were injected with approximately 200 nl of dsRNA at 1 μg/μl using a FemtoJet® injector (Eppendorf) after 4 min diethyl ether-anesthetization. Insects in the control group were injected with an equal amount of dsGFP (200 ng) in the same volume. About 60 larvae were injected for each gene in two biological replicates. In the double and triple RNAi experiments, fourth-instar larvae were injected with the same amount of dsRNA for each gene as in the single RNAi experiments.
After 3 h of recovery, the insects that died because of injection or physical injury were discarded. Larvae were kept in an incubator under standard conditions and checked every second day for mutant phenotype, mortality, time of pupation as well as time of adult eclosion. Images were taken for each insect using a Leica SD6 stereomicroscope. Elytra length and elytra gap of pupa was measured in Fiji software6. RNAi efficiency was inspected by RT-qPCR at larva [6 days post-injection (dpi)], pupa (11 dpi) and adult (16 dpi) stage. For measurement of adult motility at 20 dpi, adults were put in the middle of a square cardboard box and filmed for at least 30 s using a digital camera after a short adaption period. The recorded videos were transformed into images at one-second intervals using Free Video to JPG Converter. Then adult motility was measured in Fiji software. In brief, images were combined to the extended depth of field pictures as previously optimized (Walski et al., 2016) and measured using a segmented line. Data analysis was performed using independent t-tests in SPSS 23.
Results
Expression of Putative Genes Related to O-Glycosylation During Development
O-glycans are generated through an extended set of enzymes encoded by OGRGs. The identification of these genes and their expression patterns throughout post-embryonic development could shed light on the role of O-glycosylation during insect development. OGRGs in the genome of T. castaneum were identified through BLAST searches. Using the protein sequences for the OGRGs from D. melanogaster (Adams et al., 2000) as a query, a total of 28 putative OGRG sequences were identified in the red flour beetle genome (Supplementary Table 4). Phylogenetic trees were constructed to confirm the true orthology of the identified sequences (Supplementary Figure 1). Overall, the identified Tribolium sequences show close homology with the orthologs from L. decemlineata. Similarly, the dipterian (D. melanogaster and A. aegypti) and hemipteran (N. lugens and A. pisum) sequences cluster together in most of the trees. Remarkable was the divergence of the coleopteran Rumi cluster from the other insects, forming a separate branch in the phylogenetic tree (Supplementary Figure 1F). The clustering of pgants 2, 5, 7, and 9 suggest an ancestral evolution, with the former three even including human orthologs. In contrast, the orthologs of pgants 4, 6, 8, and 10 group in one cluster and suggest a more species specific development (Supplementary Figure 1G).
Based on their position in the O-glycosylation pathway, their involvement in the synthesis of the most abundant O-glycans or with orthologs showing clear mutant phenotypes in D. melanogaster, 19 OGRGs (Table 1) were selected for an in-depth analysis of their expression patterns (Aoki et al., 2008; Ten Hagen et al., 2009; Haltom and Jafar-Nejad, 2015). Transcription levels were quantified for the different developmental stages from egg to adult insects (Figure 2). Statistically significant differences were observed in the transcription levels for all genes (ANOVA, P < 0.05). Most of the OGRGs showed an increase in expression during the pupal stages (four in prepupa, ten in pupa and one in late pupa) (Figure 2). In contrast, mRNA levels of TcGALE1 were high during the larval stages, with the peak of expression at the sixth-instar larval stage, and a decrease during the pupal stage and early adult stage. Tcpgant6 transcription levels are low during the larval and early pupal stages, and increase during the late pupal stage, reaching a peak at the early adult stage. Several genes showed high transcription levels at the egg stage, with 5 genes showing their peak expression at this stage. While the expression profiles for genes coding for paralogous enzymes mostly showed similar patterns (pgants, POMTs, Ofuts, and OGTs), the expression profiles for TcGALE1 and TcGALE2 showed opposite patterns (Figure 2).
Figure 2. Relative expression of the OGRGs during development of Tribolium castaneum. Eggs samples (E) were collected within a period of 24 h. second to sixth-instar larvae (L2-L6) and prepupae (PP) were collected every three days after egg-hatching. At 20 and 22 days post egg-hatching, samples of pupa (P) and late pupa (LP) were collected. Early adult (EA) and mature adult (MA) insects were collected at 24 and 27 days post egg-hatching. Relative gene expression was measured by RT-qPCR and normalized using TcUbiquitin and TcRPS18 as reference genes. Lines connect the mean values of two biological replicates.
To further characterize the expression profile of the OGRGs in the red flour beetle, the transcription levels for these genes were analyzed in different tissues (head, thorax, abdomen, wings (forewings and hindwings) and legs) at the pupal and adult stage. The results showed that during the pupal stage TcGALE2, Tcpgant1, Tcpgant3, Tcpgant5, Tcpgant7, Tcpgant9, Tcpgant35A, TcOGT, TcOfut1, TcPOMT1, TcPOMT2, Tcfng, and TcRumi are highly expressed in the wings (Figure 3A). In addition, Tcpgant6, TcGALE1, TcC1GalTA, TcOfut1,TcOfut2, and TcRumi exhibited their highest expression in the thorax. Interestingly, the unique gene showing topmost expression in the legs is Tcfng. Except for Tcpgant7, TcC1GalTA and TcOfut2, most genes showed a lower expression in the abdomen during the pupal stage (Figure 3A). As observed in the temporal expression profile, most paralogous enzymes (pgants, POMTs, Ofuts, and OGTs) showed a similar pattern over the different tissues of the pupa; except for the GALEs that show opposite patterns. Moreover, POMTs, Ofuts and OGTs are nearly constitutively expressed in all parts of the pupal body (Figure 3A). At the adult stage, the tissue specificity of some genes changed (Figure 3B). The peak in Tcpgant7 expression shifted from the wings in the pupal stage to the abdomen in the adult stage. Also for TcRumi, a shift in expression to the abdomen was observed. Similarly, Tcpgant1 expression became dominant in the thorax and decreased in the wings. Expression of Tcpgant2 peaked in the head at the adult stage. GALE2, Tcpgant3, Tcpgant5, Tcpgant6, Tcpgant9, Tcpgant35A, TcC1GalTA, TcOGT, TcOfut1, TcPOMT1 and TcPOMT2 were highly expressed in the wings. In the abdomen, Tcpgant7, TcOfut2, TcPOMT2, and TcRumi showed the highest expression, while Tcpgant2, Tcpgant5, Tcpgant9, TcC1GalTA, TcOGT, TcOfut1 and TcPOMT1 showed their lowest expression levels (Figure 3B). In addition, Tcpgant2, TcC1GalTA, TcOGT, TcOfut1, TcPOMT1, Tcfng showed their peak expression in the head. TcGALE2 and TcEOGT were expressed at the highest level in the legs. Tcpgant1 and TcGALE1 showed the highest expression in the thorax. Although their spatial expression pattern was similar in the pupa, OGTs, Ofuts and POMTs showed an inverse expression pattern at the adult stage (Figure 3B). The results of the expression analysis indicated that several pgant genes have distinct temporal and spatial expression patterns (Figures 2, 3).
Figure 3. Tissue specificity of the OGRGs at pupa and adult stage of Tribolium castaneum. OGRGs expression in head, thorax, abdomen, wings (forewings and hindwings) and legs was examined at the pupa and adult stage of the red flour beetle. The different tissues were dissected under a microscope, and separately collected from 10 pupae or 20 adults. Relative gene expression was measured by RT-qPCR and normalized to reference genes (TcUbiquitin and TcRPS18). The bars represent the mean values of two biological replicates. (A) The tissue specificity of OGRGs at the pupal stage of the red flour beetle. (B) The tissue specificity of OGRGs at the adult stage.
Silencing of Genes Related to Core 1 O-Glycan Synthesis Affects Pupal Development
The core 1 O-glycan structure is the most prevalent O-glycan structure in insects and consists of a single GalNAc extended with a single Gal residue. These sugar moieties are transferred to the protein through the sequential action of pgants and core 1 β1,3-galactosyltransferases (C1GalTs) (Figure 1). DsRNA injection targeting the selected genes resulted in an efficient silencing of gene transcription (Supplementary Figure 2). As shown in Figure 4, RNAi of Tcpgant3 decreased adult mobility by 10% (P = 0.011). Silencing of Tcpgant5 caused a 33% larger pupa elytra gap compared to control (P = 0.007), and malformed adult elytra (Figure 4). RNAi of Tcpgant35A prevented pupation and eclosion, and yielded a reduction in pupa elytra length (13%, P < 0.001) and an increased elytra gap (67%, P = 0.045). About 20% of the injected insects died during the larval stages, 40% of the insects died during the transition from larva to pupa, 20% died during the transition from pupa to adult, and in total the mortality was about 90%. Pupae were malformed, and the dead insects showed appearance with adult thorax and pupa abdomen (Figure 4). RNAi of TcC1GalTA impaired adult motility (32%, P < 0.001) (Figure 5).
Figure 4. RNAi phenotypes of single RNAi experiments. RNAi of Tcpgant5 affected elytra development; RNAi of Tcpgant35A caused severe phenotypes and prevented pupation as well as eclosion. RNAi of TcOGT and TcEOGT caused strong phenotypes. RNAi of TcPOMT1 and TcPOMT2 resulted in adult wing malformation. RNAi of TcOfut1 destroyed elytra development. Scale bar is 1 mm.
Figure 5. RNAi effects on adult motility. RNAi of TcC1GalTA decreased adult mobility (n = 28, 38, 51, 57, from top to down). Error bars represent SEM of two biological replicates, *P < 0.05, **P < 0.01, ***P < 0.001.
O-GlcNAc Glycan Is Required for Pupa and Adult Development
The Tribolium genome contains two genes encoding O-GlcNAc transferases potentially involved in O-GlcNAcylation (Figure 1): OGT, transferring GlcNAc to common O-glycosylation sites of the peptide, and EOGT, which specifically transfers GlcNAc moieties to O-glycosylation sites in the EGF domain of secreted and membrane proteins (Varshney and Stanley, 2017). After injection of 4th-instar larvae with specific dsRNAs targeting either TcOGT or TcEOGT, the transcription levels of these genes were reduced by 70% (Supplementary Figure 2). Silencing of either of these genes caused conspicuous phenotypes (Figure 4) and adult mortality (Table 2). Silencing of TcEOGT resulted in antenna misorientation, a larger elytra gap compared to the control group (57%, P = 0.014), and wing blistering at pupal stage in 5% of the insects (Figure 4). The wings of adults were not closed, and the molting was not completely finished. Upon injection of dsRNAs targeting TcOGT, a misorientation of the antennas was observed in the pupa and the adult insects (Figure 4). The elytra of the pupa showed a reduction in size by 10% (P < 0.001) and a 5-fold increase in the gap distance between the pupal elytra (P < 0.001) (Figure 4) as shown in Supplementary Table 3. As was observed for TcEOGT silencing, silencing of TcOGT resulted in the inability to close the wings in the adult insects.
Table 2. Summary of phenotypes caused by larval RNAi.
Decreased Expression of O-Mannosylation Genes Affects Adult Wing Formation
Insects encode two protein O-mannosyltransferases which cooperate to the transfer of Man residues to proteins: POMT1 and POMT2 (Figure 1). Silencing of either of these two genes in T. castaneum larvae caused abnormal elytra (Figure 4) and adult mortality (Table 2). The length of the pupal elytra was found to be significantly reduced (P < 0.001) after injection of dsRNAs targeting TcPOMT1 and TcPOMT2 (10% and 9%, respectively) compared to the control treatment. In addition, adult wings were malformed and not closed. Silencing of these genes did not have an impact on the pupal elytra gap.
Knockdown of Two Putative O-Fucosyltransferases Showed Entirely Different Effects
The red flour beetle genome has two putative O-fucosyltransferases: Ofut1, which transfers Fuc specifically to O-glycosylation sites in the EGF domain of secreted proteins (including Notch, Delta, and Serrate), and Ofut2, which is responsible for O-fucosylation of thrombospondin type I repeats (TSRs) in proteins (Ten Hagen et al., 2009; Figure 1). While silencing of TcOfut1 caused severe effects on pupa and adult development, RNAi of TcOfut2 did not cause any phenotypical effect (Figure 4). Injection of dsRNAs targeting TcOfut1 resulted in a significant increase in the gap between the pupal elytra (88% compared to control, P = 0.017) and a decrease of the pupal elytra length (5%, P < 0.001).
RNAi of Rumi and Fng Caused Minor Effects
Silencing of the transcription levels for Tcfng and TcRumi, which are orthologs of the Drosophila Fuc-specific O-GlcNAc transferase (fng) and protein O-glucosyltransferase (Rumi), caused minor effects. Fng is responsible for transferring GlcNAc specifically to O-Fuc residues on the Notch protein (Moloney et al., 2000). The function of Rumi is adding Glc to serine residues in certain EGF repeats of Notch. RNAi of TcRumi resulted in a significant reduction of the pupal elytra gap by 32% (P = 0.032) as shown in Supplementary Table 3. Knockdown of Tcfng decreased the pupal elytra gap by 40% (P < 0.001). No appreciable phenotypes were observed after injecting specific dsRNA for these genes to larvae.
RNAi of TcGALE1 Reduced Adult Mobility
TcGALE1 and TcGALE2 are two putative genes encoding UDP-galactose 4-epimerases. These two enzymes catalyze two distinct but similar reactions: the epimerization of UDP-Glc to UDP-Gal, and the epimerization of UDP-GlcNAc to UDP-GalNAc. These nucleotide-activated sugars can function as sugar donors in O-glycosylation. Therefore, we assumed that knockdown of GALE can significantly reduce the related glycans simultaneously. To verify this assumption, we injected specific dsRNAs targeting either of the two GALEs into 4th-instar larvae. The induced RNAi effects persisted throughout the development and into the adult stage, and resulted in an approximately 75% reduction of transcript abundance (Supplementary Figure 2). Knockdown of either TcGALE1 or TcGALE2 caused no effects on morphology (Figure 6). However, silencing of TcGALE1 expression resulted in a decrease in the mobility of the adult insect (Figure 5), since the movement speed showed a 13% decrease compared to the control insect (P = 0.028).
Figure 6. Silencing of GALE genes blocks pupation. Individual RNAi of TcGALE1 or TcGALE2 caused minor effects, Silencing of TcGALE1 and TcGALE2 simultaneously resulted in mortality at the prepupal stage. Red arrow indicates the split of dorsal cuticle. Scale bar is 1 mm.
Single, Double or Triple RNAi Yield Similar Phenotypes
In addition to the experiments where only one of the O-glycosylation pathways was silenced, double and triple RNAi experiments were conducted to assess the redundancy of genes involved in the same pathway and the effect of blocking multiple O-glycosylation pathways on the development and survival of the insects. Several combinations were selected: Ofut1, EOGT and Rumi are enzymes known to transfer O-glycans to the EGF repeats on the extracellular domain of secreted proteins in insects, and silencing of TcOGT, TcPOMT1 and TcPOMT2 simultaneously could block all the intracellular O-glycans, except for mucin (O-GalNAc)-type O-glycans. While co-silencing usually revealed a phenotype similar to that of the silencing of a single gene, the co-silencing of GALE1 and GALE2 showed additive effects. Co-silencing of TcGALE1 and TcGALE2 simultaneously resulted in a 77%-increase in lethality before the pupal stage (Table 3). The insects that died showed a kraurotic abdomen and a split of the dorsal cuticle (Figure 6). This phenotype was similar to the one caused by RNAi of Tcpgant35A. Injection of dsRNAs targeting TcOGT-TcEOGT, Tcpgant3-TcOGT, Tcpgant35A-TcOGT, and TcOGT-TcPOMT1-TcPOMT2 caused similar phenotypes (Supplementary Figure 3). In all cases, a decrease in pupal elytra size was observed together with a misorientation of the antennae in the pupa and adult insects. At the adult stage, the insects additionally showed a bending of the body and these were unable to close their wings (Supplementary Figure 3). Also for RNAi of Tcpgant3-TcEOGT and TcOfut1-TcEOGT-TcRumi similar phenotypes were observed, co-silencing of these genes caused an increase in the elytra gap, a bending of the adult body, and a misorientation of the antennae (Supplementary Figure 3). Co-silencing of Tcpgant5 and TcOGT impaired adult wing formation and decreased larval growth as well as pupal elytra length (8%, P < 0.001). RNAi of Tcpgant3-Tcpgant5-Tcpgant35A caused pupa and adult elytra to be apart (Supplementary Figure 3). Double RNAi of TcOfut1 and TcOfut2 resulted in the malformation of the adult elytra (Supplementary Figure 3). RNAi of TcPOMT1 and TcPOMT2 caused malformation of the adult wing and a reduction in length of pupal elytra (5%, P = 0.011). In most cases, double or triple RNAi resulted in partial larval mortality and partial adult mortality (Table 3).
Table 3. Summary of phenotypes caused by double and triple RNAi.
Discussion
Identification and Expression Analysis of OGRGs in T. castaneum
A bioinformatics study resulted in the identification of 28 putative OGRGs in the genome of T. castaneum. To investigate the role of O-glycosylation in the post-embryonic development of insects, we studied the transcription profiles of 19 of the OGRGs in the red flour beetle during its development from egg to adult and in different tissues (including head, thorax, abdomen, wings and legs) at the pupa and adult stage of the red flour beetle. As observed in D. melanogaster (Hagen et al., 2003), the overall temporal expression profile showed the highest expression at the pupal stage and a higher expression at the adult stage compared to the larval stage. This profile is similar to the one observed for N-glycosylation related genes (NGRGs) in T. castaneum (Walski et al., 2016). Their increase in expression levels at the onset of metamorphosis suggested a role for N-glycosylation in the transition from larva to adult. Similarly, we hypothesized a role for the OGRGs in the post-embryonic development of insects.
The spatial profile revealed that most of the OGRGs were highly expressed in the wings both at the pupal and adult stage, suggesting that O-glycosylation is required for wing development. This is consistent with previous studies in Drosophila showing several OGRGs with a strong expression in the wing disc of the third-instar larvae (Tian and Hagen, 2006) (Gelbart and Emmert, 2013).
Functional Analysis of OGRGs During Post-embryonic Development of T. castaneum
To investigate the role of O-glycosylation in the metamorphosis and wing development of T. castaneum, we silenced the expression of the OGRG through RNAi. The results showed a broad range of phenotypes from blockage of pupation and high mortality to more mild effects on the development of wings and adult mobility (Figure 7). While many phenotypes were consistent with previous silencing and mutation experiments in D. melanogaster (reviewed in Li et al., 2020) and T. castaneum (Dönitz et al., 2015), several differences were observed with the previously described phenotypes. Although the non-lethal phenotype after silencing of Tcpgant3 is consistent with the observations in D. melanogaster and mutations in T. castaneum (Zhang et al., 2008; Dönitz et al., 2015; Yamamoto-Hino et al., 2015), differences were observed with the results in the iBeetle-database of Tribolium RNAi phenotypes (Dönitz et al., 2015). While in our hands silencing of Tcpgant3 only causes a slight reduction in adult motility, iBeetle reports partial mortality (up to 50%) and a delay or block of eclosion (Dönitz et al., 2015). These differences with the RNAi experiments in the iBeetle database could be related to the small sample size of the experiments in iBeetle where only 10 insects were injected for each target gene. The phenotypes of the other OGRGs in the synthesis of mucin-type O-glycans are consistent with the results in D. melanogaster, with RNAi of Tcpgant5 and Tcpgant35A causing severe lethality (Tran et al., 2012; Yamamoto-Hino et al., 2015) and silencing of C1GalTA causing a potential defect in the nervous system (Lin et al., 2008; Yoshida et al., 2008). Interesting, knock-down of pgant35A transcription levels resulted in major lethality during the transition from larva to pupa and from pupa to adult, with the insects that died between larva and pupa showing a split in the dorsal cuticle and the beetles that died just before the adult stage showing a cephalothorax in adult appearance while their abdomen was in pupal appearance. This indicated that in the RNAi insects pupation and eclosion were initiated but could not be completed.
Figure 7. Overview of rnai effects of O-glycosylation-related genes in the post-embryonic development of the red flour beetle Tribolium castaneum. (A) Break of the interconversion of UDP-GalNAc and UDP-Gal, UDP-GlcNAc and UDP-Glc blocked pupation. (B) RNAi of N-acetylgalactosaminyltransferases impaired pupation and eclosion. (C) Knockdown of O-GlcNAc transferases caused malformation of wing and antenna. (D) RNAi of protein O-mannosyltransferases resulted in malformation of adult wing. (E) Silencing of protein O-glucosyltransferase caused minor effects. (F) Knockdown of O-fucosyltransferases conduced to disruption of adult wing development.
While the major lethality (about 80%) and decreased larval growth, observed upon silencing of OGT are consistent with previous studies in D. melanogaster (Sinclair et al., 2009; Park et al., 2011), silencing in T. castaneum also caused a strong phenotype on wings, antennae and body, which is in agreement with the data from iBeetle-Base (Dönitz et al., 2015), being described as bent antennae and irregular wing orientation. However, some differences in the timing of the phenotype were observed for RNAi of TcEOGT compared to D. melanogaster. The pupal wing blistering in T. castaneum was also observed in D. melanogaster but at adult stage and at a high rate (Müller et al., 2013). In addition, loss of EOGT caused larval lethality in fruit flies but resulted in adult lethality in the red flour beetle (Cox and Spradling, 2006; Buszczak et al., 2007; Neely et al., 2010; Leonardi et al., 2011; Müller et al., 2013). Most of these beetles could not finish their molting completely, suggesting a novel role for EOGT during eclosion.
In contrast with the twisted abdomen observed in fruit flies with a mutant POMT (Ichimiya et al., 2004), RNAi-induced silencing of POMT in T. castaneum caused a malformation of the adult wings, suggesting novel roles for O-mannosylation in Coleoptera. In addition, RNAi of the POMTs in Tribolium caused major lethality (about 65%), whereas complete lethality and behavior defective phenotypes (flight and locomotor behavior defective) were reported in the fruit fly (Ichimiya et al., 2004; Ueyama et al., 2010). In agreement with a previous study in Drosophila (Ichimiya et al., 2004), RNAi of each or both of the genes resulted in a similar phenotype, suggesting that POMT1 and POMT2 are dependent on each other.
The effects of RNAi-induced silencing of TcOfut1 and TcOfut2 showed contrasting results. While RNAi of TcOfut1 resulted in a severe phenotype with disrupted wing development and high mortality, no phenotypes were observed after injecting dsRNAs specific to TcOfut2. The potential reason for this could be the different targets for the two enzymes (Luo et al., 2006).
Previous studies in Drosophila (Ishio et al., 2015; Matsumoto et al., 2016) suggested a redundant function for Ofut1 and Rumi, which catalyze the addition of Glc to serine residues in certain EGF repeats of Notch. However, our results of single RNAi showed that RNAi of TcOfut1 caused severe effects, while RNAi of TcRumi had a minor impact. Specifically, while the injection of dsOfut1 caused a malformation of the elytra and a high adult mortality, we only observed a smaller pupal elytra gap upon silencing of Rumi, which indicated a unidirectional complementation of Ofut1 to Rumi in T. castaneum. Combining the results of double RNAi, we expect that Ofut2 plays a role in controlling the complementation between Ofut1 and Rumi. Co-silencing Ofut2 and Ofut1 activated the complementation of Rumi to Ofut1. That could be an explanation for the low mortality in double RNAi of TcOfut1 and TcOfut2.
Interestingly, we did not observe any clear RNAi phenotypes and less than 10% mortality after the larval injection of dsTcfng, whereas in the iBeetle-Base, pupal injection resulted in 90% mortality and blocked eclosion (Dönitz et al., 2015). A reason for this discrepancy could be found in the differences in experimental setup. In the iBeetle-Base experiments, female pupa were injected with the dsTcfng (Dönitz et al., 2015), while in our experiments 4th instar larvae were treated. This emphasizes the need for comparative research to understand these differences in phenotypes after silencing the same gene but in different developmental stages.
Disruption of GALE, the enzyme catalyzing the interconversion of UDP-Gal and UDP-Glc, and UDP-GalNAc and UDP-GlcNAc, affects multiple O-glycosylation pathways as it causes a dependency not only on exogenous Gal but also on exogenous GalNAc being a necessary precursor for the synthesis of O-glycans. The high larval mortality when silencing the expression of TcGALE, is consistent with previous studies showing embryonic and larval lethality in D. melanogaster (Sanders et al., 2010; Daenzer et al., 2012). In humans, mutations in GALE result in epimerase-deficiency galactosemia, a disease characterized by liver damage, early-onset cataracts, deafness and cognitive disability (Timson, 2006). Interestingly, the two putative genes of GALE (TcGALE1, TcGALE2) identified in the genome of T. castaneum, showed an opposite expression pattern. Whereas RNAi of either TcGALE1 or TcGALE2 caused almost no effects on viability, co-silencing of both TcGALE1 and TcGALE2 significantly reduced the viability at the prepupal stage. This suggests that the proteins coded by these two genes have, at least partially overlapping functions. In addition, the split on the dorsal cuticle indicates that pupation was initiated but not completed, suggesting that GALE or O-glycans consisting out of these four sugars play an important role in pupation.
The O-glycosylation is complicated and diverse. For example, the insect genome encodes more than 12 pgant genes, that are functionally redundant in transferring GalNAc to proteins. In this case, double or triple RNAi could be a more powerful strategy to study gene function. However, from an overall perspective, most double and triple RNAi experiments resulted in almost identical effects compared to single RNAi. This was against our expectations. When performing the double or triple RNAi experiments, we injected the same amount of the individual dsRNAs as we did in the single RNAi experiments and similar RNAi efficiencies were obtained for each gene in the triple RNAi experiments as in the single RNAi experiments. The results observed in the triple RNAi experiments are consistent over independent biological repeats. A reasonable explanation for the phenotypes observed in the triple RNAi of the pgants could be that these pgants or their targets probably function in a particular order in biological processes. For instance, if Pgant5 plays a role upstream in the biological processes, then silencing of Tcpgant5 already disrupted the pathway, and consequently, this triple RNAi showed a similar phenotype as RNAi of Tcpgant5. Another possible explanation is that due to the increased pressure by silencing three genes a compensation effect takes place. A similar observation was made when silencing the catalytic subunit of the N-glycosylation oligosaccharyl transferase complex, after silencing of one isoform, the expression of the other isoform increased (De Schutter et al., 2019). Nevertheless, it is interesting that all double or triple RNAi experiments involving TcOGT resulted in a significantly higher larval mortality than the single RNAi, while the larval mortality was only minor in the single RNAi of TcOGT. In contrast, RNAi of Tcpgant35A and TcOGT caused 22% lower larval mortality than the single RNAi of Tcpgant35A. These results indicate a complicated interplay between different OGRGs.
In conclusion, we demonstrated that O-glycosylation is essential for survival, metamorphosis and wing development of the red flour beetle, T. castaneum. Moreover, we report that the core 1 O-glycan synthesized by pgant35A is required for pupation and eclosion. Our results provide an insight into the functions of O-glycosylation during the post-embryonic development of insects.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
WL performed the gene expression analysis, RNAi experiments, and data analyses. KDS and WL performed the phylogenetic analysis. EVD, GS, and KDS were involved in data analysis and discussion and manuscript corrections. All authors contributed to manuscript and approved the final manuscript.
Funding
This work is supported by the Special Research Fund from Ghent University (Belgium). WL is a recipient of a Chinese Scholarship Council (CSC) doctoral grant.
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2021.629682/full#supplementary-material
Supplementary Figure 1 | Phylogeny of putative proteins involved in O-glycosylation in Tribolium castaneum. Phylogenetic tree was constructed in MEGA X using Maximum Likelihood method. The best fitting model LG + G (A–E,G) or WAG + F + I (F,H,I) was determined based on the BIC score. The phylogeny was tested with 1000 bootstrap replicates. Bootstrap values are indicated at the branches. Hs, Homo sapiens; Dm, Drosophila melanogaster; Tc, Tribolium castaneum; Nl, brown planthopper (Nilaparvata lugens); Ld, colorado potato beetle (Leptinotarsa decemlineata); Bt, large earth bumblebee (Bombus terrestris); Ap, pea aphid (Acyrthosiphon pisum); Bm, silk moth (Bombyx mori). (A) GALE orthologs. (B) C1GalT orthologs. (C) OGT and EOGT orthologs. (D) Ofut orthologs. (E) POMT orthologs. (F) orthologs of fng and Rumi. (G) pgant orthologs. (H) orthologs of proteins involved in synthesis of heparan sulfate and chondroitin sulfate. (I) all OGRG orthologs from (A–H).
Supplementary Figure 2 | Efficiency of gene silencing. Fourth-instar larvae were injected with approx. 200 ng dsRNA. DsGFP was injected in control insects. Total RNA was extracted at 6 d (larva), 11 d (pupa), 16 d (adult) after injection. Relative gene expression was measured by qPCR using Ubiquitin and RPS18 as reference genes. Error bar is SD (n = 2). ∗, the gene expression was not examined because of high mortality.
Supplementary Figure 3 | Phenotypes caused by double and triple RNAi in Tribolium castaneum. In the double and triple RNAi experiments, fourth-instar larvae were injected with a same amount of dsRNA for each gene as in the single RNAi experiments. DsGFP was injected in the control. Scale bar is 1 mm.
Supplementary Table 1 | Primers used for gene expression analysis.
Supplementary Table 2 | Primers used for dsRNA synthesis. taatacgactcactatagggaga – T7 promoter sequence.
Supplementary Table 3 | Pupal elytra length and gap after dsRNA injection at larva stage.
Supplementary Table 4 | List of all identified genes.
Supplementary Table 5 | Gene silencing efficiency of double and triple RNAi.
Footnotes
- ^ http://ibeetle-base.uni-goettingen.de/blast/
- ^ https://www.ncbi.nlm.nih.gov/guide/proteins/
- ^ https://www.ebi.ac.uk/interpro/download/
- ^ https://blast.ncbi.nlm.nih.gov
- ^ http://www.bioinformatics.nl/cgi-bin/primer3plus/primer3plus.cgi
- ^ https://fiji.sc/
References
Adams, M. D., Celniker, S. E., Holt, R. A., Evans, C. A., Gocayne, J. D., Amanatides, P. G., et al. (2000). The genome sequence of Drosophila melanogaster. Science 287, 2185–2195. doi: 10.1126/science.287.5461.2185
Aoki, K., Porterfield, M., Lee, S. S., Dong, B., Nguyen, K., Mcglamry, K. H., et al. (2008). The diversity of O-linked glycans expressed during Drosophila melanogaster development reflects stage- and tissue-specific requirements for cell signaling. J. Biol. Chem. 283, 30385–30400. doi: 10.1074/jbc.M804925200
Ashwood, C., Abrahams, J. L., Nevalainen, H., and Packer, N. H. (2017). Enhancing structural characterisation of glucuronidated O−linked glycans using negative mode ion trap−HCD mass spectrometry. Rapid Commun. Mass Spectrom. 31, 851–858. doi: 10.1002/rcm.7851
Buszczak, M., Paterno, S., Lighthouse, D., Bachman, J., Planck, J., Owen, S., et al. (2007). The carnegie protein trap library: a versatile tool for Drosophila developmental studies. Genetics 175, 1505–1531. doi: 10.1534/genetics.106.065961
Chaffey, P. K., Guan, X., Chen, C., Ruan, Y., Wang, X., Tran, A. H., et al. (2017a). Structural insight into the stabilizing effect of O-glycosylation. Biochemistry 56, 2897–2906. doi: 10.1021/acs.biochem.7b00195
Chaffey, P. K., Guan, X., Wang, X., Ruan, Y., and Biochemistry, L.-Y. (2017b). Quantitative effects of O-linked glycans on protein folding. Biochemistry 56, 4539–4548. doi: 10.1021/acs.biochem.7b00483
Cox, R. T., and Spradling, A. C. (2006). Milton controls the early acquisition of mitochondria by Drosophila oocytes. Development 133, 3371–3377. doi: 10.1242/dev.02514
Daenzer, J. M., Sanders, R. D., Hang, D., and Fridovich-Keil, J. L. (2012). UDP-galactose 4’-epimerase activities toward UDP-Gal and UDP-GalNAc play different roles in the development of Drosophila melanogaster. PLoS Genet. 8:e1002721. doi: 10.1371/journal.pgen.1002721
De Pourcq, K., De Schutter, K., and Callewaert, N. (2010). Engineering of glycosylation in yeast and other fungi: current state and perspectives. Appl. Microbiol. Biotechnol. 87, 1617–1631. doi: 10.1007/s00253-010-2721-1
De Schutter, K., Chen, P., Shen, Y., Van Damme, E. J. M., and Smagghe, G. (2019). The OST-complex as target for RNAi-based pest control in Nilaparvata lugens. Arch. Insect. Biochem. Physiol. 101:e21555. doi: 10.1002/arch.21555
Dell, A., Galadari, A., Sastre, F., and Hitchen, P. (2010). Similarities and differences in the glycosylation mechanisms in prokaryotes and eukaryotes. Int. J. Microbiol. 2010:148178. doi: 10.1155/2010/148178
Dönitz, J., Schmitt-Engel, C., Grossmann, D., Gerischer, L., Tech, M., Schoppmeier, M., et al. (2015). iBeetle-Base: a database for RNAi phenotypes in the red flour beetle Tribolium castaneum. Nucleic Acids Res. 43, D720–D725. doi: 10.1093/nar/gku1054
Eichler, J., and Koomey, M. (2017). Sweet new roles for protein glycosylation in prokaryotes. Trends Microbiol. 25, 662–672. doi: 10.1016/j.tim.2017.03.001
Feliciello, I., Akrap, I., and Ugarkovic, D. (2015). Satellite DNA modulates gene expression in the beetle Tribolium castaneum after heat stress. PLoS Genet. 11:e1005547. doi: 10.1371/journal.pgen.1005466
Fraichard, S., Bouge, A. L., Kendall, T., Chauvel, I., Bouhin, H., and Bunch, T. A. (2010). Tenectin is a novel alphaPS2betaPS integrin ligand required for wing morphogenesis and male genital looping in Drosophila. Dev. Biol. 340, 504–517. doi: 10.1016/j.ydbio.2010.02.008
Fukuda, M. (2002). Roles of mucin-type O-glycans in cell adhesion. Biochim. Biophys. Acta 1573, 394–405. doi: 10.1016/S0304-4165(02)00409-9
Gelbart, W., and Emmert, D. (2013). Flybase High Throughput Expression Pattern Data. FlyBase Analysis. Availabe at: https://flybase.org/reports/FBrf0212041.html (accessed October 29, 2013).
Hagen, T. K. G., Tran, D. T., Gerken, T. A., and Stein, D. S. (2003). Functional characterization and expression analysis of members of the UDP-GalNAc: polypeptide N-acetylgalactosaminyltransferase family from Drosophila melanogaster. J. Biol. Chem. 278, 35039–35048. doi: 10.1074/jbc.M303836200
Haltom, A. R., and Jafar-Nejad, H. (2015). “O-Linked glycans in Drosophila development: overview,” in Glycoscience: Biology and Medicine, eds N. Taniguchi, T. Endo, G. Hart, P. Seeberger, and C. H. Wong (Tokyo: Springer), 809–815.
Ichimiya, T., Manya, H., Ohmae, Y., Yoshida, H., Takahashi, K., Ueda, R., et al. (2004). The twisted abdomen phenotype of Drosophila POMT1 and POMT2 mutants coincides with their heterophilic protein O-mannosyltransferase activity. J. Biol. Chem. 279, 42638–42647. doi: 10.1074/jbc.M404900200
Ishio, A., Sasamura, T., Ayukawa, T., Kuroda, J., Ishikawa, H. O., Aoyama, N., et al. (2015). O-fucose monosaccharide of Drosophila Notch has a temperature-sensitive function and cooperates with O-glucose glycan in Notch transport and Notch signaling activation. J. Biol. Chem. 290, 505–519. doi: 10.1074/jbc.M114.616847
Jones, P., Binns, D., Chang, H. Y., Fraser, M., Li, W., Mcanulla, C., et al. (2014). InterProScan 5: genome-scale protein function classification. Bioinformatics 30, 1236–1240. doi: 10.1093/bioinformatics/btu031
Kieliszewski, M. J. (2001). The latest hype on Hyp-O-glycosylation codes. Phytochemistry 57, 319–323. doi: 10.1016/S0031-9422(01)00029-2
Kong, E. C., Allouche, L., Chapot, P. A., Vranizan, K., Moore, M. S., Heberlein, U., et al. (2010). Ethanol-regulated genes that contribute to ethanol sensitivity and rapid tolerance in Drosophila. Alcohol. Clin. Exp. Res. 34, 302–316. doi: 10.1111/j.1530-0277.2009.01093.x
Kurz, S., Aoki, K., Jin, C. S., Karlsson, N. G., Tiemeyer, M., Wilson, I. B. H., et al. (2015). Targeted release and fractionation reveal glucuronylated and sulphated N- and O-glycans in larvae of dipteran insects. J. Proteomics 126, 172–188. doi: 10.1016/j.jprot.2015.05.030
Leonardi, J., Fernandez-Valdivia, R., Li, Y.-D., Simcox, A. A., and Jafar-Nejad, H. (2011). Multiple O-glucosylation sites on Notch function as a buffer against temperature-dependent loss of signaling. Development 138, 3569–3578. doi: 10.1242/dev.068361
Li, W., De Schutter, K., Van Damme, E. J. M., and Smagghe, G. (2020). Synthesis and biological roles of O-glycans in insects. Glycoconj. J. 37, 47–56. doi: 10.1007/s10719-019-09867-1
Lin, Y. R., Reddy, B., and Irvine, K. D. (2008). Requirement for a core 1 galactosyltransferase in the Drosophila nervous system. Dev. Dyn. 237, 3703–3714. doi: 10.1002/dvdy.21775
Lord, J. C., Hartzer, K., Toutges, M., and Oppert, B. (2010). Evaluation of quantitative PCR reference genes for gene expression studies in Tribolium castaneum after fungal challenge. J. Microbiol. Methods 80, 219–221. doi: 10.1016/j.mimet.2009.12.007
Luo, Y., Koles, K., Vorndam, W., Haltiwanger, R. S., and Panin, V. M. (2006). Protein O-fucosyltransferase 2 adds O-fucose to thrombospondin type 1 repeats. J. Biol. Chem. 281, 9393–9399. doi: 10.1074/jbc.M511975200
Martinez, M. R., Dias, T. B., Natov, P. S., and Zachara, N. E. (2017). Stress-induced O-GlcNAcylation: an adaptive process of injured cells. Biochem. Soc. Trans. 45, 237–249. doi: 10.1042/BST20160153
Matsumoto, K., Ayukawa, T., Ishio, A., Sasamura, T., Yamakawa, T., and Matsuno, K. (2016). Dual roles of O-glucose glycans redundant with monosaccharide O-fucose on Notch in Notch trafficking. J. Biol. Chem. 291, 13743–13752. doi: 10.1074/jbc.M115.710483
Matsumoto, K., Ishio, A., and Matsuno, K. (2015). “O-fucose glycan in Drosophila Notch signaling,” in Glycoscience: Biology and Medicine, eds N. Taniguchi, T. Endo, G. Hart, P. Seeberger, and C. H. Wong (Tokyo: Springer), 841–847.
Maynard, J. C., Pham, T., Zheng, T. L., Jockheck-Clark, A., Rankin, H. B., Newgard, C. B., et al. (2010). Gp93, the Drosophila GRP94 ortholog, is required for gut epithelial homeostasis and nutrient assimilation-coupled growth control. Dev. Biol. 339, 295–306. doi: 10.1016/j.ydbio.2009.12.023
Moloney, D. J., Panin, V. M., Johnston, S. H., Chen, J., Shao, L., Wilson, R., et al. (2000). Fringe is a glycosyltransferase that modifies Notch. Nature 406, 369–375. doi: 10.1038/35019000
Müller, R., Jenny, A., and Stanley, P. (2013). The EGF repeat-specific O-GlcNAc-transferase Eogt interacts with notch signaling and pyrimidine metabolism pathways in Drosophila. PLoS One 8:e62835. doi: 10.1371/journal.pone.0062835
Mummery-Widmer, J. L., Yamazaki, M., Stoeger, T., Novatchkova, M., Bhalerao, S., Chen, D., et al. (2009). Genome-wide analysis of Notch signalling in Drosophila by transgenic RNAi. Nature 458, 987–992. doi: 10.1038/nature07936
Nagarkar-Jaiswal, S., Lee, P. T., Campbell, M. E., Chen, K., Anguiano-Zarate, S., Gutierrez, M. C., et al. (2015). A library of MiMICs allows tagging of genes and reversible, spatial and temporal knockdown of proteins in Drosophila. eLife 4:e05338. doi: 10.7554/eLife.05338
Nakamura, N., Stalnaker, S. H., Lyalin, D., Lavrova, O., Wells, L., and Panin, V. M. (2010). Drosophila Dystroglycan is a target of O-mannosyltransferase activity of two protein O-mannosyltransferases, Rotated Abdomen and Twisted. Glycobiology 20, 381–394. doi: 10.1093/glycob/cwp189
Neely, G. G., Hess, A., Costigan, M., Keene, A. C., Goulas, S., Langeslag, M., et al. (2010). A genome-wide Drosophila screen for heat nociception identifies alpha 2 delta 3 as an evolutionarily conserved pain gene. Cell 143, 628–638.
Ogawa, M., Sawaguchi, S., Furukawa, K., and Okajima, T. (2015). N-acetylglucosamine modification in the lumen of the endoplasmic reticulum. Biochim. Biophys. Acta 1850, 1319–1324. doi: 10.1016/j.bbagen.2015.03.003
Okajima, T., Xu, A., Lei, L., and Irvine, K. D. (2005). Chaperone activity of protein O-fucosyltransferase 1 promotes notch receptor folding. Science 307, 1599–1603. doi: 10.1126/science.1108995
Park, S., Park, S.-H., Baek, J. Y., Jy, Y. J., Kim, K. S., Roth, J., et al. (2011). Protein O-GlcNAcylation regulates Drosophila growth through the insulin signaling pathway. Cell. Mol. Life Sci. 68, 3377–3384. doi: 10.1007/s00018-011-0640-7
Sanders, R. D., Sefton, J. M. I., Moberg, K. H., and Fridovich-Keil, J. L. (2010). UDP-galactose 4 epimerase (GALE) is essential for development of Drosophila melanogaster. Dis. Models Mech. 3, 628–638. doi: 10.1242/dmm.005058
Schjoldager, K. T., Vester-Christensen, M. B., Goth, C. K., Petersen, T. N., Brunak, S., Bennett, E. P., et al. (2011). A systematic study of site-specific GalNAc-type O-glycosylation modulating proprotein convertase processing. J. Biol. Chem. 286, 40122–40132. doi: 10.1074/jbc.M111.287912
Schwientek, T., Mandel, U., Roth, U., Muller, S., and Hanisch, F. G. (2007). A serial lectin approach to the mucin-type O-glycoproteome of Drosophila melanogaster S2 cells. Proteomics 7, 3264–3277. doi: 10.1002/pmic.200600793
Shcherbata, H. R., Yatsenko, A. S., Patterson, L., Sood, V. D., Nudel, U., Yaffe, D., et al. (2007). Dissecting muscle and neuronal disorders in a Drosophila model of muscular dystrophy. EMBO J. 26, 481–493. doi: 10.1038/sj.emboj.7601503
Sinclair, D. A., Syrzycka, M., Macauley, M. S., Rastgardani, T., Komljenovic, I., Vocadlo, D. J., et al. (2009). Drosophila O-GlcNAc transferase (OGT) is encoded by the Polycomb group (PcG) gene, super sex combs (sxc). Proc. Natl. Acad. Sci. U.S.A. 106, 13427–13432. doi: 10.1073/pnas.0904638106
Soya, S., Şahar, U., and Karaçalı, S. (2016). Monosaccharide profiling of silkworm (Bombyx mori L.) nervous system during development and aging. Invert. Neurosci. 16:8. doi: 10.1007/s10158-016-0191-6
Spiro, R. G. (2002). Protein glycosylation: nature, distribution, enzymatic formation, and disease implications of glycopeptide bonds. Glycobiology 12, 43R–56R. doi: 10.1093/glycob/12.4.43R
Spiro, R. G., Lucas, F., and Rudall, K. M. (1971). Glycosylation of hydroxylysine in collagens. Nat. New Biol. 231, 54–55. doi: 10.1038/newbio231054a0
Ten Hagen, K. G. (2015). “O-glycosylation and development,” in Glycoscience: Biology and Medicine, eds N. Taniguchi, T. Endo, G. Hart, P. Seeberger, and C. H. Wong (Tokyo: Springer), 1065–1072.
Ten Hagen, K. G., Zhang, L., Tian, E., and Zhang, Y. (2009). Glycobiology on the fly: developmental and mechanistic insights from Drosophila. Glycobiology 19, 102–111. doi: 10.1093/glycob/cwn096
Tian, E., and Hagen, K. G. (2006). Expression of the UDP-GalNAc: polypeptide N-acetylgalactosaminyltransferase family is spatially and temporally regulated during Drosophila development. Glycobiology 16, 83–95. doi: 10.1093/glycob/cwj051
Timson, D. J. (2006). The structural and molecular biology of type III galactosemia. IUBMB Life 58, 83–89. doi: 10.1080/15216540600644846
Toutges, M. J., Hartzer, K., and Lord, J. (2010). Evaluation of reference genes for quantitative polymerase chain reaction across life cycle stages and tissue types of Tribolium castaneum. J. Agric. Food Chem. 58, 8948–8951. doi: 10.1021/jf101603j
Tran, D. T., Zhang, L., Zhang, Y., Tian, E., Earl, L. A., and Ten Hagen, K. G. (2012). Multiple members of the UDP-GalNAc: polypeptide N-acetylgalactosaminyltransferase family are essential for viability in Drosophila. J. Biol. Chem. 287, 5243–5252. doi: 10.1074/jbc.M111.306159
Ueyama, M., Akimoto, Y., Ichimiya, T., Ueda, R., Kawakami, H., Aigaki, T., et al. (2010). Increased apoptosis of myoblasts in Drosophila model for the Walker-Warburg syndrome. PLoS One 5:e11557. doi: 10.1371/journal.pone.0011557
Varshney, S., and Stanley, P. (2017). EOGT and O-GlcNAc on secreted and membrane proteins. Biochem. Soc. Trans. 45, 401–408. doi: 10.1042/BST20160165
Vuerinckx, K., Verlinden, H., Lindemans, M., Broeck, J. V., and Huybrechts, R. (2011). Characterization of an allatotropin-like peptide receptor in the red flour beetle, Tribolium castaneum. Insect Biochem. Mol. Biol. 41, 815–822. doi: 10.1016/j.ibmb.2011.06.003
Walski, T., De Schutter, K., Van Damme, E. J. M., and Smagghe, G. (2017). Diversity and functions of protein glycosylation in insects. Insect Biochem. Mol. Biol. 83, 21–34. doi: 10.1016/j.ibmb.2017.02.005
Walski, T., Van Damme, E. J. M., Smargiasso, N., Christiaens, O., De Pauw, E., and Smagghe, G. (2016). Protein N-glycosylation and N-glycan trimming are required for postembryonic development of the pest beetle Tribolium castaneum. Sci. Rep. 6:35151. doi: 10.1038/srep35151
Yamamoto-Hino, M., Yoshida, H., and Ichimiya, T. (2015). Phenotype−based clustering of glycosylation−related genes by RNAi−mediated gene silencing. Genes Cells 20, 521–542. doi: 10.1111/gtc.12246
Yang, S., Höti, N., Yang, W., Liu, Y., Chen, L., Li, S., et al. (2017). Simultaneous analyses of N-linked and O-linked glycans of ovarian cancer cells using solid-phase chemoenzymatic method. Clin. Proteomics 14:3. doi: 10.1186/s12014-017-9137-1
Yatsenko, A. S., Marrone, A. K., and Shcherbata, H. R. (2014). miRNA-based buffering of the cobblestone-lissencephaly-associated extracellular matrix receptor dystroglycan via its alternative 3′-UTR. Nat. Commun. 5:4906. doi: 10.1038/ncomms5906
Yoshida, H., Fuwa, T. J., Arima, M., Hamamoto, H., Sasaki, N., Ichimiya, T., et al. (2008). Identification of the Drosophila core 1 1, 3-galactosyltransferase gene that synthesizes T antigen in the embryonic central nervous system and hemocytes. Glycobiology 18, 1094–1104. doi: 10.1093/glycob/cwn094
Zarschler, K., Janesch, B., Pabst, M., Altmann, F., Messner, P., and Schaffer, C. (2010). Protein tyrosine O-glycosylation-a rather unexplored prokaryotic glycosylation system. Glycobiology 20, 787–798. doi: 10.1093/glycob/cwq035
Zhang, L., Tran, D. T., and Ten Hagen, K. G. (2010). An O-glycosyltransferase promotes cell adhesion during development by influencing secretion of an extracellular matrix integrin ligand. J. Biol. Chem. 285, 19491–19501. doi: 10.1074/jbc.M109.098145
Keywords: O-glycosylation, oligosaccharide, gene expression, RNAi, Tribolium castaneum, metamorphosis
Citation: Li W, De Schutter K, Van Damme EJM and Smagghe G (2021) RNAi-Mediated Silencing of Pgants Shows Core 1 O-Glycans Are Required for Pupation in Tribolium castaneum. Front. Physiol. 12:629682. doi: 10.3389/fphys.2021.629682
Received: 15 November 2020; Accepted: 01 March 2021;
Published: 24 March 2021.
Edited by:
Peng He, Guizhou University, ChinaReviewed by:
Jia Fan, State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection (CAAS), ChinaBin Li, Nanjing Normal University, China
Yoonseong Park, Kansas State University, United States
Copyright © 2021 Li, De Schutter, Van Damme and Smagghe. 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: Guy Smagghe, guy.smagghe@ugent.be