PLoS Genetics
Public Library of Science
image
A new neuropeptide insect parathyroid hormone iPTH in the red flour beetle Tribolium castaneum
DOI 10.1371/journal.pgen.1008772, Volume: 16, Issue: 5,
Abstract

Vertebrate parathyroid hormone (PTH) and its receptors have been extensively studied with respect to their function in bone remodeling and calcium metabolism. Insect parathyroid hormone receptors (iPTHRs) have been previously described as counterparts of vertebrate PTHRs, however, they are still orphan receptors for which the authentic ligands and biological functions remain unknown. We describe an insect form of parathyroid hormone (iPTH) by analyzing its interactions with iPTHRs. Identification of this new insect peptidergic system proved that the PTH system is an ancestral signaling system dating back to the evolutionary time before the divergence of protostomes and deuterostomes. We also investigated the functions of the iPTH system in a model beetle Tribolium castaneum by using RNA interference. RNA interference of iPTHR resulted in defects in wing exoskeleton maturation and fecundity. Based on the differential gene expression patterns and the phenotype induced by RNAi, we propose that the iPTH system is likely involved in the regulation of exoskeletal cuticle formation and fecundity in insects.

Xie, Sang, Song, Zhang, Kim, Veenstra, Park, Li, and Schoofs: A new neuropeptide insect parathyroid hormone iPTH in the red flour beetle Tribolium castaneum

Introduction

Discoveries of neuropeptides and their receptors in various taxa in the postgenomic era have provided nearly comprehensive lists that provide crucial information for understanding their evolutionary processes and physiological functions. In addition, the development of new molecular techniques has rapidly expanded the knowledge of their functions in numerous cases of ancestral taxa. Homology-based searches in the existing genomic information are a powerful method; however, they may fail in cases where relatively rapid evolution occurs, thus leaving gaps in the knowledge due to punctuated equilibria in evolution. Functional studies of ancestral bilaterian neuropeptides have been successful by starting from the sequence similarities, e.g., vasopressin and thyrotropin-releasing hormone [14]

Vertebrate parathyroid hormone (PTH), the most important regulator of calcium ion homeostasis, and its receptors have been extensively studied in bone remodeling and calcium metabolism [5, 6]. Multiple PTH receptors in different phyla of vertebrates are known to be the consequence of multiple gene duplications and losses [7, 8]. Likewise, the ligand PTH also has undergone gene duplication and losses in vertebrates [9]. The PTH receptor lineage is traceable to before the time of the deuterostome-protostome split in the basal Bilateria [1012], but there are no obvious homologous PTH ligands in the basal protostome lineages. Therefore, the PTH receptors in the basal lineages of Eukaryotes are still orphan receptors for which the authentic ligands and biological functions remain unknown. In insects, likewise, our previous study has described two receptors in the red flour beetle, Tribolium castaneum , as an orthologous group similar to the PTH-receptors of vertebrates in the G protein-coupled receptor (GPCR) B group [13, 14]. The ligand for these receptors has not yet been identified [15].

In this study, we have uncovered the authentic ligand for the insect PTH receptor (iPTHR). Our strategy for deorphanization of this ancestral GPCR utilized a genomic survey of numerous species in insects and further basal lineages. We looked for the similar patterns of taxonomic distributions of iPTHR and of neuropeptide-like sequences. We found that the taxonomic distribution of the neuropeptide previously described as PXXXamide (where X is any amino acid) in the cuttlefish Sepia officinalis (Mollusca) [16], is found in a mirror image of the distribution pattern of iPTHR. Specifically, the loss of iPTHR in Diptera and Lepidoptera coincided with the absence of this specific neuropeptide in the same taxa. We demonstrated that this previously uncharacterized neuropeptide is an active ligand on these receptors and we named the peptide insect PTH (iPTH). The phenotypes from RNA interferences (RNAi) and the RNAseq data suggest that the iPTH system is involved in the maturation of the exoskeletal cuticle in the wings at the time of adult eclosion.

Results

Evolution and diversity of iPTH and iPTHR in insects

In our survey of iPTHs, the neuropeptide originally described as PXXXamide [16], we identified iPTHs in many species of insect and in other arthropod taxa, but lacking in Diptera and Lepidoptera genome sequences. In T. castaneum , the gene encoding iPTH was found in the linkage group LG6 (Fig 1A) [17, 18]. The conceptual translation contains a signal peptide at the N-terminus, followed by a dibasic cleavage site KR (Lys-Arg) and the conserved mature peptide ending with putative amidation and a dibasic cleavage signal (GRR) with further extension of the C-terminal sequence. Although this mature peptide has no sequence homology to vertebrate PTH (lower than 12% identities with PTH, PTHrP, and Tip39 of Homo sapiens ), highly conserved iPTHs were found in various insects and other taxa (Fig 1B and S1 Table). In the predicted mature peptides, N-terminal amino acids SD (Ser and Asp at positions 7 and 8 in the alignment Fig 1B) are highly conserved. The C-terminal P (Pro at the −4 position) for the original description of the peptide as PXXXamide was found to be mildly conserved in an expanded survey of the sequences. The C-terminal amidation and cleavage signal GXX (X = positively charged amino acids)[19] are strictly conserved in different species.

Conserved iPTH and iPTHR sequences.
Fig 1
(A) The iPTH gene structure and the conceptual translation in T. castaneum. (B) An alignment showing the consensus of iPTH in various insect species. The letters with black backgrounds are identical and gray are similar amino acids with the 50% majority rule. (C) The sequence logo of the high to moderate levels of conserved residues in the predicted mature peptides in various insect species.Conserved iPTH and iPTHR sequences.

The iPTHRs belong to the class B1 G protein-coupled receptors (GPCRs) and are considered counterparts of vertebrate PTHRs based on sequence similarity (Fig 2A; S1 Fig). The evolutionary processes of vertebrate PTHRs have previously been described as the events of multiple gene duplications and deletions [7, 8]. We also found support for gene duplication events of iPTHR in the insect lineages (Fig 2A; S1 Fig). The genome sequences of many insect species revealed two or more copies of the PTH receptor, although it is difficult to resolve whether the gene duplication was a single event or independent multiple events because of incongruent branching patterns of the tree constructed based on the amino acid sequences. In T. castaneum, the genome data revealed two iPTHRs at LG4 and LG3 for Tc-PTHR1 and R2 , respectively. There were no iPTHR homologs in Diptera and Lepidoptera genome sequences, consistent with the absence of iPTH sequence in those taxa (S1 Fig; S1 Table).

(A) Phylogenetic tree of the PTHRs of the metazoan species constructed by the Maximum likelihood method with the JTT model (see more details in S1 Fig for 1000 bootstrapping tests using MEGA 7.0 software). (B) Dose-response curves of the iPTH on the iPTHRs expressed in the CHO cells. The EC50 of Tc-iPTH to Tc-iPTHR1 and Tc-iPTHR2 were 63.29 nM and 413.36 nM, respectively.
Fig 2
(A) Phylogenetic tree of the PTHRs of the metazoan species constructed by the Maximum likelihood method with the JTT model (see more details in S1 Fig for 1000 bootstrapping tests using MEGA 7.0 software). (B) Dose-response curves of the iPTH on the iPTHRs expressed in the CHO cells. The EC50 of Tc-iPTH to Tc-iPTHR1 and Tc-iPTHR2 were 63.29 nM and 413.36 nM, respectively.

Functional expression of PTH receptors and iPTH activity assays

The Tribolium gene annotations (Tcas5.2) [18, 20] were used to obtain the open reading frames of the Tc-iPTHR1 (TC008110) and R2 (TC010267) by polymerase chain reaction (PCR) (GenBank with the accession numbers XM_008193931.2 and XM_964860.3, respectively) (S2A Fig). Conservation of seven cysteines (marked as stars on top of the alignment in S2B Fig) in the predicted extracellular loops including the N-terminal extracellular region implies a significant role of disulfide bridges for formation of the ligand binding pocket, as stabilization of extracellular structure by disulfide bonds in GPCR-B has previously shown [2124]. The C-terminal region showed large diversity in the sequences compared to the transmembrane regions. The iPTHRs were functionally expressed in a Chinese hamster ovarian (CHO) cell line to which components of aequorin as the responder to calcium mobilization in CHO cells had been added, which is an approach that has been described before [25, 26]. We found that the CHO cells transfected for expression of either receptor, Tc-iPTHR1 or Tc-iPTHR2, were activated by chemically synthesized Tc-iPTH. Tc-iPTHR1 is ~7× more sensitive (EC50 = 63 nM) than Tc-iPTHR2 (EC50 = 413 nM). The cells without the Tc-iPTHR constructs did not respond to the Tc-iPTH at 10 μM (Fig 2B).

Expression patterns of Tc-iPTH and Tc-iPTHRs

Transcript levels of the genes, Tc-iPTH, Tc-iPTHR1, and Tc-iPTHR2, were investigated to understand spatial and temporal expression patterns. Quantitative real time reverse transcriptase PCR (Q-PCR) revealed high levels of the transcripts in late larval and late adult stages for Tc-iPTHR1, the late adult stage for Tc-iPTHR2, and the late pupal stage for Tc-iPTH (Fig 3). The gut and central nervous system (CNS) had high transcript levels of all three genes in both the larva and adult stages (Fig 3D to 3I). Further investigations with separate foregut, midgut, and hindgut in larval, pupal, and adult stages found that the transcript levels were highest in the pupal foregut and midgut (S3A and S3B Fig). Interestingly, the Tc-iPTHR2 transcript was found to be moderately expressed in adult elytra and epidermis (Fig 3H). The Tc-iPTHR2 transcript was further confirmed by Q-PCR for the samples of pupal and adult elytra, and we found that adult elytra, but not pupal elytra, showed a moderate level of the Tc-iPTHR2 transcript (S3D Fig).

Temporal and spatial expression patterns of Tc-iPTH, Tc-iPTHR1, and Tc-iPTHR2 measured by Q-PCR.
Fig 3
(A-C) Expression patterns at different developmental stages. Different developmental stages are: EE, early embryos; LE, late embryos; EL, early larvae; LL, late larvae; EP, early pupae; LP, late pupae; EA, early adults; LA, late adults.(D-F) Expression patterns in various tissues from late larvae. The tissues are epidermis, gut, fat body and CNS (central nervous system). (G-I) Expression patterns in various tissues from late adults. The Tribolium ribosomal protein S3 (rps3) transcript was the internal control. The averages and standard errors are for three biological replicates. Different lowercase letters indicate statistically significant differences (p≤0.05).Temporal and spatial expression patterns of Tc-iPTH, Tc-iPTHR1, and Tc-iPTHR2 measured by Q-PCR.

Immunohistochemistry of Tc-iPTH revealed a number of neural cells and axonal projections in the central nervous systems of larval (Fig 4A to 4G, 4R and 4S), and pupal ganglia (Fig 4H to 4J, 4R and 4S). Visualization of cell bodies in the Tc-iPTH staining rarely occurred while the neural projections always appeared. Inconsistent immunohistochemical staining patterns of cell bodies depending on the samples were common for various neuropeptides in T. castaneum including Natalisin [25]. Combining the results from multiple samples (n>50 in each larvae and pharate adults), we concluded that there are at least two segmental ganglionic cells for Tc-iPTH expression, including the brain, subesophageal ganglion, and other segmental ganglia (Fig 4G to 4J and 4S). Two pairs of dorso longitudinal projections connected these varicosities (Fig 4D and 4E), while in some cases, transverse projections of neurons from ganglia were also observed (Fig 4B, 4C and 4S). The brain had extensive varicosities in the protocerebral and deutocerebral regions in both larva (Fig 4G and 4R) and pupa (Fig 4H and 4R).

Immunohistochemistry of Tc-iPTH in the central nervous system (A-J) and gut (K-Q) in larval (A-G, K-M), and pharate adult (H-J, N-Q) stages. The arrowheads (G-J) indicate positive neural cells for anti- iPTH, and the white arrows (K, M, O, Q) indicate positive neuronal projections for anti- iPTH. Bars are 200 μm. Diagrams of the brain (R), ganglion (S) showed the neuronal cells and axonal projections in the CNS. Diagrams of anterior part of midgut (T) showed stomatogastric nervous projections. The posterior part of midgut (U) showed the enteroendocrine cells.
Fig 4
Immunohistochemistry of Tc-iPTH in the central nervous system (A-J) and gut (K-Q) in larval (A-G, K-M), and pharate adult (H-J, N-Q) stages. The arrowheads (G-J) indicate positive neural cells for anti- iPTH, and the white arrows (K, M, O, Q) indicate positive neuronal projections for anti- iPTH. Bars are 200 μm. Diagrams of the brain (R), ganglion (S) showed the neuronal cells and axonal projections in the CNS. Diagrams of anterior part of midgut (T) showed stomatogastric nervous projections. The posterior part of midgut (U) showed the enteroendocrine cells.

Strong immunoreactivities were also found in the gut (Fig 4K to 4Q, 4T and 4U) including stomatogastric nervous projections and the puncta surrounding the projections on the anterior half of the midgut ran along the longitudinal axis (Fig 4K and 4T), which are similar to the bipolar enteric plexus cells in Manduca sexta [27]. Immunoreactive stomatogastric nerves and the puncta were also observed in the pharate adults (Fig 4N and 4Q). A separate neural projection was found in the anterior part of the ileum (Fig 4L, 4M and 4Q). In addition to the neural projection in innervation of the guts, enteroendocrine cells in the gut were also immunoreactive in both larva and pharate adults (Fig 4K, 4L, 4M, 4P and 4T). Large numbers of these cells were shown in the anterior and posterior part of the midgut (Fig 4T and 4U), while a small number of immunoreactive enteroendocrine cells were also observed in the middle part of the gut. Pharate adults also had positive reactions in very small enteroendocrine cells (<~3 μM in diameter) (Fig 4P). Preimmune serum used as the negative control did not show any specific staining patterns (S4A to S4G Fig).

Immunohistochemistry of the receptor Tc-iPTHR1, which has shown a strong activity by Tc-iPTH in the receptor-ligand activity assay, was carried out. Strong immunostaining was observed in the neural cells and axonal projections in the CNS in larval (Fig 5A to 5C), pupal (Fig 5D to 5F) and adult stage (Fig 5G to 5I). One or two pairs of the neural cells were consistently captured in the brain (Fig 5A and 5G) and thoracic ganglia (Fig 5B, 5D and 5H) in different individual samples. Immunoreactive projections also appeared in the thoracic ganglia (Fig 5B, 5E and 5H), as well as abdominal ganglia (Fig 5C, 5F and 5I). The immunoreactivity in the gut showed large stomatogastric neural cell bodies located on the anterior tip of midgut (Fig 5K) and their projections running longitudinally on the midgut in the anterior half (Fig 5J, 5L and 5M). The staining pattern is similar to the immunoreactivity of the ligand Tc-iPTH, but generally thinner than the staining for Tc-iPTH. There were no immunoreactive enteroendocrine cells. Negative controls of the immunohistochemistry of anti-Tc-iPTHR1 antibody staining were carried out using the preimmune serum, which did not show any staining in the CNS and gut (S5 Fig). Since both antibodies, anti-Tc-iPTH and anti-Tc-iPTHR1, were raised from rabbit, we did not perform the double staining.

Immunohistochemistry of Tc-iPTHR1 in the central nervous system (A-I) and gut (J-M) in larval (A-C, J, K), pupal (D-F, L) and adult (G-I, M) stages. The arrowheads indicate positive neural cells or neuronal projections for anti-iPTHR1. Bars (A-I) are 50 μm and bars (JLM) are 200 μm.
Fig 5
Immunohistochemistry of Tc-iPTHR1 in the central nervous system (A-I) and gut (J-M) in larval (A-C, J, K), pupal (D-F, L) and adult (G-I, M) stages. The arrowheads indicate positive neural cells or neuronal projections for anti-iPTHR1. Bars (A-I) are 50 μm and bars (JLM) are 200 μm.

RNAi of Tc-iPTHRs resulting in eclosion and wing defects

Both larval RNAi and pupal RNAi of Tc-iPTHR1 and Tc-iPTHR2 , which successfully suppressed the expression of their specific target genes (Fig 6A and 6D), resulted in partial defects in adult eclosion (Fig 6B and 6C). The eclosion deficiencies were associated with incomplete shedding of the pupal cuticle, failure in spreading the wings, and irregularly wrinkled elytra (Fig 6G; Fig 7). RNAi of Tc-iPTHR1 in the larval stage resulted in 81.3±8.4% of individuals with ecdysial arrest at the time of adult eclosion (Fig 6B; Fig 6G) and 14.3±5.8% showing adult morphological deficiency (Fig 6B; Fig 7A). Tc-iPTHR2 RNAi showed 43.8±10.5% individuals arrested at the time of eclosion (Fig 6B; Fig 6G), and 8.0±2.8% had adult morphological deficiency (Fig 6B; Fig 7A and 7B). RNAi with the dsRNA injections at the pupal stage resulted in the same phenotypes but with lower numbers of individuals affected (Fig 6C). The individuals with Tc-iPTHR1 RNAi laid significantly fewer eggs, 2 eggs compared to 7 eggs per day in the control (Fig 6E). The percentages of hatching eggs were also significantly reduced to 10 and 40% in RNAi of Tc-iPTHR1 and Tc-iPTHR2 , respectively, compared to ~80% hatching rate in the control (Fig 6F).

RNAi results in suppression of the transcript levels with wing and fecundity phenotypes.
Fig 6
Different treatment group are wild-type, IB (buffer injection), dsiPTH, dsiPTHR1 and dsiPTHR2. (A) Larval RNAi suppresses only the target transcripts. (B) The percentage with abnormal phenotypes in the beetles among all the animals after larval RNAi (C) and after pupal RNAi. (D) Pupal RNAi suppresses only the target transcripts. (E) Reduced egg numbers in RNAi of iPTHR1 and (F) hatching rate of the offspring after pupal RNAi. (G) RNAi phenotypes of the beetles arrested at the time of eclosion and abnormal wings. Both elytra and hindwings are shown.RNAi results in suppression of the transcript levels with wing and fecundity phenotypes.
Phenotypes after dsRNA treatments for Tc-iPTHRs.
Fig 7
The pictures are taken of the individuals who survived the dsRNA treatments but with some morphological abnormalities (adult morphological deficiency, AMD). (A) Survival after RNAi of iPTHR1 showed incomplete closure of the elytra and rough elytra surface and survival after RNAi of iPTHR2 often showed incomplete tanning of the elytra. (B) The elytra of wild-type adult and dsiPTHR2 survival. (C) Scanning electron microscopy (SEM) of elytra after RNAi with the 10-day-old adults.Phenotypes after dsRNA treatments for Tc-iPTHRs.

However, the RNAi of the ligand Tc-iPTH did not show any phenotype (Fig 6B, 6C, 6E and 6F). Immunohistochemistry indicated that the suppression of Tc-iPTH by the RNAi may not have been sufficient to induce a phenotype (S4H to S4K Fig).

Gene expression changes induced by the RNAi of Tc-iPTHR1 and Tc-iPTHR2

RNAseq of the individuals with the RNAi of Tc-iPTHR1 and Tc-iPTHR2 revealed large numbers of differentially expressed (DE) genes. Principal component analysis (PCA) of gene expression patterns in the control, dsiPTHR1, and dsiPTHR2 showed that the RNAi of Tc-iPTHR1 and Tc-iPTHR2 resulted in highly correlated patterns of gene expression changes compared to the control (Fig 8A; S6 Fig). Likewise, three-dimensional scatter plots also showed that DE genes in the dsiPTHR1 and dsiPTHR2 treatments have similar patterns (Fig 8B). When the DE genes (>4× up or down) were separately plotted for the log2 of the expression ratios, the correlation between dsiPTHR1 and dsiPTHR2 treatments was obvious (R2 = 0.78) (Fig 8C). Large portions of DE genes were common in response to the dsiPTHR1 and dsiPTHR2 treatments (110 and 242 genes for >4× up or down genes, respectively, Fig 8D). Manual examination followed by computational annotation and gene ontology term enrichment tests of the list found some noticeable enriched gene groups in the list of DE genes.

RNAseq showing the similar pattern of changes in the expression of genes in dsiPTHR1 and dsiPTHR2 compared to the control.
Fig 8
(A) Principal component analysis of gene expression patterns in control, dsiPTHR1, and dsiPTHR2. (B) Three-dimensional scatter plots analysis of differentially expressed genes in contrast with the control, dsiPTHR1 and dsiPTHR2. (C) Differentially expressed genes (>4× up or down) separately plotted for the log2 of the expression ratios. The expression changes in dsiPTHR1 and dsiPTHR2 share a similar pattern. (D) Venn diagrams show the common genes differentially expressed in dsiPTHR1 and dsiPTHR2 treatments.RNAseq showing the similar pattern of changes in the expression of genes in dsiPTHR1 and dsiPTHR2 compared to the control.

To verify the RNAseq-based DE data, we randomly selected 22 genes for each treatment, and Q-PCR was performed. The Q-PCR results were highly correlated with the RNAseq-based analyses by having R2 >0.96 in the linear regression (S7 Fig).

In the gene ontology (GO) term enrichment analyses for the DE genes, the largest portion encodes extracellular structural components including cuticular proteins; 68 are down- and 40 genes are up-regulated (Fig 8D; S1 Dataset and S2 Dataset). The 68 down-regulated genes were mainly cuticular proteins and included keratins and mucins, while the 40 up-regulated genes included 17 cuticle proteins and six different keratins and seven prismalin-14s. In the phylogeny of down-regulated cuticle proteins, it showed that most of the cuticle proteins belong to the CPR (Cuticular Protein with the Rebers and Riddiford Consensus) family, in which 10 of them belong to RR-1 cluster and 18 are in RR-2 cluster (S8 Fig).

In addition, genes for cuticular metabolism were also captured, in which chitinase 8 was up-regulated and six genes down-regulated (including chorion peroxidase-like, chitinase 5, chitinase 10, and knickkopf3). Three genes encoding cytochrome P-450 (Cyp450) were up-regulated, while five Cyp450 genes including Halloween gene spook (Cyp450 307A1) were significantly down-regulated. Four different antimicrobial peptides were up-regulated. Down-regulated genes were ecdysone-induced 74EF, E78, DHR3, HR4, E3, neuropeptide F receptor, and rickets, which codes the bursicon receptor. Six genes encoding sensory molecular components, gustatory receptor, odorant receptor, and chemosensory proteins, were also down-regulated.

RNAi of Tc-iPTHRs increased chitin content but did not affect 20E amount

Based on the wrinkled elytra phenotype and the large numbers of genes encoding cuticle protein in the DEs after the RNAi of Tc-iPTHR1 and Tc-iPTHR2, we measured the chitin content and 20-hydroxyecdysone (20E) in the control and in the dsiPTHR1 and dsiPTHR2 beetles. Significantly higher amounts of chitin, measured by the amount of hydrolyzed N-acetylglucosamine (GlcNAc), were found in the RNAi of Tc-iPTHR1 (Fig 9A and 9B). More specially, chitin contents in the isolated elytra showed RNAi of both Tc-iPTHR1 and Tc-iPTHR2 were significantly higher than that in the control (Fig 9C).

(A-B) Chitin content of the individuals in the LP (A) and LA (B) stage after dsRNA treatment. (C) Chitin content of the tissues of elytra, gut and epidermis in the adults after larval dsRNA treatment. The statistical results are for the treatments compared to the controls. *, p ≤ 0.05; **, p ≤ 0.01 and ***, p ≤ 0.001.
Fig 9
(A-B) Chitin content of the individuals in the LP (A) and LA (B) stage after dsRNA treatment. (C) Chitin content of the tissues of elytra, gut and epidermis in the adults after larval dsRNA treatment. The statistical results are for the treatments compared to the controls. *, p ≤ 0.05; **, p ≤ 0.01 and ***, p ≤ 0.001.

The 20E titer was measured in two different developmental stages, 3-day-old and 5-day-old pupa. No significant effect on the 20E titer after RNAi of Tc-iPTH, Tc-iPTHR1 or Tc-iPTHR2 was found (S9 Fig), using a20E antibody that has a relatively high specificity against 20E (10× higher affinity for 20E than that for ecdysone) [28].

Discussion

The initial purpose of this study was to identify the authentic ligand for the iPTHR, which is an ortholog of vertebrate PTHR. Similar phylogenetic distributions of the genes and the presence/absence of the neuropeptide candidates and iPTHR in the survey of different taxa allowed initial identification of iPTH. Activity of iPTH on the iPTHRs showed that iPTH does indeed acts a ligand for these receptors. RNAi of iPTHRs in T. castaneum suggests that iPTH signaling is required for cuticle maturation in the adult eclosion. It is interesting to find a function of iPTHR in exoskeletal development in T. castaneum while the major PTH system function in vertebrates includes skeleton development, although these findings could well be a consequence of analogous functional evolution.

The RNAi phenotypes for Tc-iPTHR1 and Tc-iPTHR2 together with the DE analyses after the RNAi provided insights into the functions of the iPTH system. The similarities of the RNAi phenotypes and the DE gene patterns after the RNAi of Tc-iPTHR1 and Tc-iPTHR2 indicate that although these two receptors may serve complementary functions, both are necessary as they are unable to compensate each other.

The mechanism behind the requirement of iPTHR for exoskeletal cuticle development is still enigmatic. Based on the DE data obtained 5 days after the injections of dsiPTHR (~prepupal stages), noticeable DE genes with biological significance for wing development are in two different categories: the genes for cuticle formation and signaling pathways. The major category of the DE genes consisted of genes involved in cuticle formation such as the proteins making up the structural components of the exoskeleton and chitin metabolism. Among the down-regulated genes, chitinase 5 and chitinase 10 were previously shown to affect degradation of old cuticle at the molting stage, ultimately leading to molting deficiency in T. castaneum [29, 30], while knickkopf 3 and tyrosine hydroxylases are also well characterized for their essential roles in adult cuticle formations and sclerotization [31, 32].

Among the down-regulated cuticle proteins, it has been reported that some of the CPRs are highly expressed in the elytra and hindwings of newly molted Tribolium adults (4 and 9 respectively) [33]. Additionally, three CPAP (Cuticular Proteins Analogous to Peritrophins) proteins identified here have been reported to be expressed only in cuticle-forming tissues, and at relatively high levels in the elytra and hindwings of the 2 week-old Tribolium adults [34], in which TC011140 and TC009263 knockdown can result in adult lethality and adults with rough elytra unable to completely cover the abdomen, respectively [35]. A similar phenotype was also observed after RNAi of Tc-iPTHR1 and Tc-iPTHR2 . Among down-regulated extracellular matrix proteins, keratin is a component of skin and epidermal appendages of vertebrates [36] and prismalin-14 [37] is a protein characterized as a novel matrix protein in pearl oyster.

In addition, up-regulated expressions of Chitinase 8, 9, 16, and 3-like were found after RNAi of Tc-iPTHRs. Three of them (Chitinase 8, 9, and 16) belong to the Group IV chitinases, in which Chitinase 8 contains only one chitin binding domain while Chitinase 9 and 16 lack this domain [38]. It was reported there were no detectable abnormalities after RNAi of these three chitinases individually [29]. These results suggested that they are not essential while the phenotypes of over/mis-expressions of these genes are potentially disruptive in T. castaneum development. Likewise, a number of cuticular proteins were also up-regulated. These data together with the large number of down-regulated cuticular proteins indicate that the expressions of different sets of cuticular proteins were disrupted in the RNAi of Tc-iPTHRs. The conclusions drawn from the DE data analysis were also supported by cuticle contents changing in the RNAi pupae and adults, especially in the abnormal adults’ elytra, as the direct causal factor of the phenotype during the eclosion.

The down-regulated genes involved in signaling pathways included the ecdysone metabolism genes spook (Cyp307A1) and Cyp303A1 which has an unknown substrate. Spook is required in the early step of ecdysteroid synthesis for ketodiol formation in D. melanogaster and in Schistocerca gregaria [39, 40]. Cyp303A1 in D. melanogaster is a gene required for bristle morphogenesis in cuticle metabolism [41] and its knockdown culminated in insects with a phenotype of defects in wing development at the adult eclosion [42]. Moreover, ecdysone-inducible genes [43, 44] 74EF, E78, DHR3, HR4, and E3 and Ftz-F1 were also down-regulated (S1 Dataset). However, when the individuals with the iPTHR RNAi were examined in the stage where 20E is reaching its peak (3 days after pupation), we were unable to detect any differences in 20E immunoreactivity compared to the control. The current study and the knowledge available on ecdysteroid functions are insufficient to draw a conclusion as to whether iPTH signaling is required for ecdysteroid metabolism including the biosynthesis or inactivation steps.

In addition, other signaling molecule receptors were also down-regulated: neuropeptide F receptor, rickets, and CG13579, an amine receptor. Among these, rickets is well characterized as the receptor for bursicon that is required for cuticle sclerotization and wing expansion at the time of eclosion [4549]. Involvement of neuropeptide F and biogenic amine in metamorphic development are as yet unknown.

The phenotype after RNAi of Tc-iPTHRs appeared at eclosion and early adult stage, while significant DE genes started to appear in the earlier prepupal stage. It is likely that the phenotype is a consequence of the culmination of accumulated physiological and developmental imbalances occurring in the earlier stages. In addition, although Q-PCR result showed a similar marked reduction of both Tc-iPTH and its receptors when the RNAi is applied, the inhibition may have not been sufficient for the phenotype in the RNAi of Tc-iPTH as the dsRNA of Tc-iPTH did not show any phenotype. We found that the RNAi for suppression of Tc-iPTH still displayed residual iPTH immunoreactivity, indicating that the Tc-iPTH RNAi did not provide an appropriate level of transcript knockdown and the leaking level of protein production was sufficient for the required levels of bioactivity. An alternative possibility is that there may be another ligand for PTH receptors that yet needs to be discovered in T. castaneum.

In summary, we identified a new insect peptidergic system, iPTH, which is the orthologous system to vertebrate PTH, showing that the PTH system is an ancestral signaling system dating back to the evolutionary time before the divergence of protostomes and deuterostomes, although it has been lost in Diptera and Lepidoptera. High expression at both mRNA and protein levels of the iPTH and its receptors in the gut and the CNS of Tribolium indicated that this neuropeptide may have important regulatory functions. We investigated this by RNAi and showed that iPTH system plays an essential role in cuticle maturation and fecundity in insects.

Materials and methods

Experimental insects and phylogenetic analysis

A laboratory colony of T. castaneum Georgia-1 line insects was used for all of these experiments and reared in whole wheat flour with 5% brewer yeast powder at 30°C and 40% relative humidity on a 14h:10h light cycle under standard conditions [50].

Blast searches of iPTH and iPTHR sequences were made in the databases of NCBI: nr, RefSeq genome, genome, Expressed sequence tag, Whole-genome shotgun contigs, Transcriptome shotgun assembly, and of i5k workspace. The search algorithm and parameters were BLOSUM60 or PAM30 with the expected threshold of 50, followed by manual screens.

Transmembrane domains of the iPTHR sequences were determined by the on-line software TMHMM (v2.0) (http://www.cbs.dtu.dk/services/TMHMM/) [51]. Multiple sequence alignments of the amino acid sequences were performed using ClustalW [52]. The alignments were viewed and edited in BOXSHADE software (http://www.ch.embnet.org/software/BOX_form.html). The phylogenetic trees were constructed using MEGA 7.0 software using the Maximum likelihood method of the JTT model with 1000 bootstrapping tests [53].

Q-PCR analyses for temporal and spatial expression patterns

For the stage-specific expression study, total RNAs were extracted from pools of multiple individuals in the following developmental stages: early embryo (EE, 1 day old, ~50 mg), late embryo (LE, 3 day old, ~50 mg), early larva (EL, 1 day old, ~50mg), late larva (LL, ~15 days old, 3 individuals), early pupa (EP, 1 day old, 3 individuals), late pupa (LP, 5 days old, 3 individuals), early adult (EA, 1 day old, 3 individuals), and late adult (LA, 10 days old, 3 individuals). For the tissue-specific expression study, total RNA of various tissues from approximately 100 LL and 100 LA was also extracted. The LL tissues included epidermis, gut, fat body, central nervous system (CNS), and in LA, besides the four tissues mentioned, accessory gland, testis and ovary were also collected. For the further expression investigations in the gut and epidermal system, separate foregut, midgut, hindgut, epidermis of head and thorax, epidermis of abdomen in the LL, LP, and LA stages were prepared, and in the LP and LA stage, the elytra were also collected. Reverse transcription was performed by using 1 μg total RNA. Q-PCR analyses were performed in MiniOpticon (Bio-Rad) with the SYBR green kit (TaKaRa, Japan). The relative mRNA levels were standardized to the control gene ribosomal protein S3 (rpS3, GenBank accession numbers CB335975), using the 2−ΔΔCT method [54]. The primers are listed in S2 Table.

RNA interference

Systemic RNAi was performed as previously described [55]. Primers containing Tc-iPTHR1 or Tc-iPTHR2 sequences and the T7 polymerase promoter at the 5’-end of both the forward primer and reverse primer were used (S2 Table). The regions that were selected to amplify dsDNA showed low sequence similarity between Tc-iPTHR1 and Tc-iPTHR2 to minimize the potential cross activities of the RNAi. The PCR product was used as the template for dsRNA synthesis with the HiScribe T7 Quick High Yield RNA Synthesis Kit (New England Biolabs, USA). 30 beetles in two different stages, early 6th instar stage and 1th day pupal stage, were used for injections of 200 nl of 2 μg/μl dsRNA per beetle. Negative controls consisted of non-injection (wild- type group, WT) or injection of an equal volume of buffer only (IB group). RNAi efficiency was examined by Q-PCR using a pool of three individuals on the 5th day after dsRNA injections. The beetles were daily observed and mortality and phenotype were recorded. Egg numbers and hatching rate were obtained from 14–16 day old beetles of pupal RNAi. At least three replications were carried out for each RNAi phenotype count.

GPCR assays

Tc-iPTHR1 and Tc-iPTHR2 cloning were based on TC008110 and TC010267 gene predictions (GenBank accession number XP_008192153.1 and XP_969953.1, respectively). The gene was cloned into pcDNA3.1+ and was used for heterologous expression in Chinese hamster ovarian (CHO) cells. Transient transfection was followed by an aquorine based functional assay as previously described [1, 26].

Immunohistochemistry and microscopy

Tissue samples, CNS and gut from late larva, pupa(pharate adult) and early adult were dissected in PBS (pH 7.4). The tissues were fixed for 1~2 h at room temperature in 4% paraformaldehyde in PBS. After extensive washing, the tissues were incubated in primary antibody (1:1000 for anti-Tc-iPTH or 1:2000 for anti-Tc-iPTHR1) for 20 h at room temperature followed by incubation in secondary antibody (1:300 for goat anti-rabbit IgG-Alexa 488) overnight. Images were acquired with a Zeiss LSM 700. Anti-Tc-iPTH and anti-Tc-iPTHR1 was separately raised in a rabbit after three immunizations with the CTMPVGGGRIDPGRIamide (C-terminus 14 amino acids of iPTH) or SNEYIKWTRNNYATRQCamide (N-terminus 17 amino acids of iPTHR1) conjugated to keyhole limpet hymocyanin (KLH). The final bleed was purified by using the antigenic peptide. Pre-immune serum was used as the negative control.

The elytra phenotypes of iPTHR1 and iPTHR2 dsRNA injections at the EP stage were obtained for 10 day-old adults. The samples were dehydrated in an ethanol gradient of 60, 70, 80, 90, and 100 for 5 min each. Samples were coated with platinum and observed with a KYKY10008B electron microscope (JSM-5610, Tokyo, Japan) [56].

RNA-seq analysis

Total RNA of 6 larvae was isolated at day 5 after dsRNA injections in the early 6th instar stage using a Trizol Reagent (TaKaRa, Japan) according to the manufacturer’s protocol. The mRNA was enriched by using oligo(dT) magnetic beads. After fragmentation into ~200 bp lengths, it was used for the first strand cDNA synthesis by random hexamers. Library construction and sequencing was performed by BGI-Tech with Illumina HiSeq2000 [57]. Raw sequence data were submitted to the Short Read Archive (SRA) database of NCBI under the accession numbers SRR8559987 (dsiPTHR1), SRR8559988 (dsiPTHR2) and SRR1176913 (control).

The raw data were analyzed after filtering the low-quality sequences. Sequences were aligned to the T. castaneum genome (http://www.beetlebase.org/) using SOAPaligner/soap2. The expression level of genes from the RNA sequencing was normalized by the RPKM method (Reads Per Kilobase of exon model per Million mapped reads). The number of total mapped reads was approximately 7 million for each library. GO terms for each T. castaneum genes were obtained using Blast2GO (version 2.3.5) (http://www.blast2go.org/). Only those with more than 4× differences were kept for the DE analyses.

20E titer and chitin content assay

For measuring 20E titer, a single individual was homogenized in 500 μl high grade 100% methanol by using a plastic pestle, and then the samples were placed at 4°C overnight with gentle rocking. The tubes were spun at 13,000 rpm and the supernatant was collected. One more round of methanol extraction of the pellet was applied. Pooled methanol extractions were concentrated in the Eppendorf Centrifugal Vacuum Concentrator 5301 (Hamburg, Germany). Enzyme immunoassay for 20E quantification was performed as was previously described [28].

To detect the effect of the RNAi of iPTHRs on chitin synthesis, the insects treated with RNAi in the EP stage were collected on the 5th day after pupation and the 10th day after eclosion. Epidermal tissues (elytra, epidermis of abdomen, epidermis of head and thorax, and gut) of the 10th day adults were dissected. Chitin content of 3 pooled whole body or tissues of 10 individuals were assayed based on the method previously described [58]. Briefly, chitin is decomposed into chitosan under the action of high temperature (130°C) and alkali (14 M KOH), and then the chitosan could be depolymerized to glucosamine residues under a mixture of 10% NaNO2 and 10% KHSO4. After addition of NH4SO3NH2 and freshly prepared MBTH (3-methyl-2-benzothiazolone hydrazone hydrochloride hydrate), the mixtures were incubated at 100°C for 5 min in a water bath. Absorbance at 650 nm was determined in the Synergy H1 microplate reader (BioTek Instruments Inc., USA). Chitin content was expressed as a glucosamine equivalent according to a standard curve constructed by using known concentrations of glucosamine (Sigma-Aldrich).

References

1 

MJ Aikins, DA Schooley, K Begum, M Detheux, RW Beeman, Y Park. . Vasopressin-like peptide and its receptor function in an indirect diuretic signaling pathway in the red flour beetle. Insect Biochem Mol Biol. 2008;38(7):, pp.740–8. , doi: 10.1016/j.ibmb.2008.04.006

2 

I Beets, T Janssen, E Meelkop, L Temmerman, N Suetens, S Rademakers, et al. Vasopressin/oxytocin-related signaling regulates gustatory associative learning in C. elegans. Science. 2012;338(6106):, pp.543–5. , doi: 10.1126/science.1226860

3 

E Van Sinay, O Mirabeau, G Depuydt, MB Van Hiel, K Peymen, J Watteyne, et al. Evolutionarily conserved TRH neuropeptide pathway regulates growth in Caenorhabditis elegans. Proc Natl Acad Sci U S A. 2017;114(20):, pp.E4065–E74. , doi: 10.1073/pnas.1617392114

4 

B Li, R Predel, S Neupert, F Hauser, Y Tanaka, G Cazzamali, et al. Genomics, transcriptomics, and peptidomics of neuropeptides and protein hormones in the red flour beetle Tribolium castaneum. Genome Res. 2008;18(1):, pp.113–22. , doi: 10.1101/gr.6714008

5 

B Lanske, AC Karaplis, K Lee, A Luz, A Vortkamp, A Pirro, et al. PTH/PTHrP receptor in early development and Indian hedgehog-regulated bone growth. Science. 1996;273(5275):, pp.663–6. , doi: 10.1126/science.273.5275.663

6 

AD White, F Fang, FG Jean-Alphonse, LJ Clark, HJ An, H Liu, et al. Ca(2+) allostery in PTH-receptor signaling. Proc Natl Acad Sci U S A. 2019;116(8):, pp.3294–9. , doi: 10.1073/pnas.1814670116

7 

JS On, BK Chow, LT Lee. . Evolution of parathyroid hormone receptor family and their ligands in vertebrate. Front Endocrinol (Lausanne). 2015;6:, pp.28.

8 

PL Pinheiro, JC Cardoso, DM Power, AV Canario. . Functional characterization and evolution of PTH/PTHrP receptors: insights from the chicken. BMC Evol Biol. 2012;12:, pp.110, doi: 10.1186/1471-2148-12-110

9 

P Suarez-Bregua, L Cal, C Canestro, J Rotllant. . PTH Reloaded: A New Evolutionary Perspective. Front Physiol. 2017;8:, pp.776, doi: 10.3389/fphys.2017.00776

10 

O Mirabeau, JS Joly. . Molecular evolution of peptidergic signaling systems in bilaterians. Proc Natl Acad Sci U S A. 2013;110(22):, pp.E2028–37. , doi: 10.1073/pnas.1219956110

11 

G Jekely. . Global view of the evolution and diversity of metazoan neuropeptide signaling. Proc Natl Acad Sci U S A. 2013;110(21):, pp.8702–7. , doi: 10.1073/pnas.1221833110

12 

MR Elphick, O Mirabeau, D Larhammar. . Evolution of neuropeptide signalling systems. J Exp Biol. 2018;221(Pt 3).

13 

CJ Li, M Chen, M Sang, X Liu, W Wu, B Li. . Comparative genomic analysis and evolution of family-B G protein-coupled receptors from six model insect species. Gene. 2013;519(1):, pp.1–12. , doi: 10.1016/j.gene.2013.01.061

14 

AJ Harmar. . Family-B G-protein-coupled receptors. Genome Biol. 2001;2(12):, pp.REVIEWS3013, doi: 10.1186/gb-2001-2-12-reviews3013

15 

Y Tanaka, Y Suetsugu, K Yamamoto, H Noda, T Shinoda. . Transcriptome analysis of neuropeptides and G-protein coupled receptors (GPCRs) for neuropeptides in the brown planthopper Nilaparvata lugens. Peptides. 2014;53:, pp.125–33. , doi: 10.1016/j.peptides.2013.07.027

16 

C Zatylny-Gaudin, V Cornet, A Leduc, B Zanuttini, E Corre, G Le Corguille, et al. Neuropeptidome of the Cephalopod Sepia officinalis: Identification, Tissue Mapping, and Expression Pattern of Neuropeptides and Neurohormones during Egg Laying. J Proteome Res. 2016;15(1):, pp.48–67. , doi: 10.1021/acs.jproteome.5b00463

17 

HS Kim, T Murphy, J Xia, D Caragea, Y Park, RW Beeman, et al. BeetleBase in 2010: revisions to provide comprehensive genomic information for Tribolium castaneum. Nucleic Acids Res. 2010;38:, pp.D437–D42. , doi: 10.1093/nar/gkp807

18 

S Richards, RA Gibbs, GM Weinstock, SJ Brown, R Denell, RW Beeman, et al. The genome of the model beetle and pest Tribolium castaneum. Nature. 2008;452(7190):, pp.949–55. , doi: 10.1038/nature06784

19 

JA Veenstra. . Mono- and dibasic proteolytic cleavage sites in insect neuroendocrine peptide precursors. Arch Insect Biochem. 2000;43(2):, pp.49–63.

20 

J Donitz, C Schmitt-Engel, D Grossmann, L Gerischer, M Tech, M Schoppmeier, et al. iBeetle-Base: a database for RNAi phenotypes in the red flour beetle Tribolium castaneum. Nucleic Acids Res. 2015;43(D1):, pp.D720–D5.

21 

A Wilmen, B Goke, R Goke. . The isolated N-terminal extracellular domain of the glucagon-like peptide-1 (GLP)-1 receptor has intrinsic binding activity. FEBS Lett. 1996;398(1):, pp.43–7. , doi: 10.1016/s0014-5793(96)01214-8

22 

P Gaudin, A Couvineau, JJ Maoret, C Rouyer-Fessard, M Laburthe. . Mutational analysis of cysteine residues within the extracellular domains of the human vasoactive intestinal peptide (VIP) 1 receptor identifies seven mutants that are defective in VIP binding. Biochem Biophys Res Commun. 1995;211(3):, pp.901–8. , doi: 10.1006/bbrc.1995.1897

23 

LJ Qi, AT Leung, Y Xiong, KA Marx, AB Abou-Samra. . Extracellular cysteines of the corticotropin-releasing factor receptor are critical for ligand interaction. Biochemistry. 1997;36(41):, pp.12442–8. , doi: 10.1021/bi970997r

24 

V Karageorgos, M Venihaki, S Sakellaris, M Pardalos, G Kontakis, MT Matsoukas, et al. Current understanding of the structure and function of family B GPCRs to design novel drugs. Hormones (Athens). 2018;17(1):, pp.45–59.

25 

HB Jiang, A Lkhagva, I Daubnerova, HS Chae, L Simo, SH Jung, et al. Natalisin, a tachykinin-like signaling system, regulates sexual activity and fecundity in insects. P Natl Acad Sci USA. 2013;110(37):, pp.E3526–E34.

26 

D Kim, L Simo, Y Park. . Molecular characterization of neuropeptide elevenin and two elevenin receptors, IsElevRl and IsElevR2, from the blacklegged tick, Ixodes scapularis. Insect Biochem Molec. 2018;101:, pp.66–75.

27 

PF Copenhaver, PH Taghert. Neurogenesis in the insect enteric nervous system: Generation of pre-migratory neurons from an epithelial placode Development. 1990;109:, pp.17–28.

28 

K Ruang-Rit, Y Park. . Endocrine system in supernumerary molting of the flour beetle, Tribolium freemani, under crowded conditions. Insect Biochem Mol Biol. 2018;101:, pp.76–84. , doi: 10.1016/j.ibmb.2018.08.002

29 

Q Zhu, Y Arakane, RW Beeman, KJ Kramer, S Muthukrishnan. . Functional specialization among insect chitinase family genes revealed by RNA interference. Proc Natl Acad Sci U S A. 2008;105(18):, pp.6650–5. , doi: 10.1073/pnas.0800739105

30 

Q Zhu, Y Arakane, RW Beeman, KJ Kramer, S Muthukrishnan. . Characterization of recombinant chitinase-like proteins of Drosophila melanogaster and Tribolium castaneum. Insect Biochem Mol Biol. 2008;38(4):, pp.467–77. , doi: 10.1016/j.ibmb.2007.06.011

31 

SS Chaudhari, MY Noh, B Moussian, CA Specht, KJ Kramer, RW Beeman, et al. Knickkopf and retroactive proteins are required for formation of laminar serosal procuticle during embryonic development of Tribolium castaneum. Insect Biochem Mol Biol. 2015;60:, pp.1–6. , doi: 10.1016/j.ibmb.2015.02.013

32 

MJ Gorman, Y Arakane. . Tyrosine hydroxylase is required for cuticle sclerotization and pigmentation in Tribolium castaneum. Insect Biochem Mol Biol. 2010;40(3):, pp.267–73. , doi: 10.1016/j.ibmb.2010.01.004

33 

NT Dittmer, Y Hiromasa, JM Tomich, N Lu, RW Beeman, KJ Kramer, et al. Proteomic and transcriptomic analyses of rigid and membranous cuticles and epidermis from the elytra and hindwings of the red flour beetle, Tribolium castaneum. J Proteome Res. 2012;11(1):, pp.269–78. , doi: 10.1021/pr2009803

34 

S Jasrapuria, Y Arakane, G Osman, KJ Kramer, RW Beeman, S Muthukrishnan. . Genes encoding proteins with peritrophin A-type chitin-binding domains in Tribolium castaneum are grouped into three distinct families based on phylogeny, expression and function. Insect Biochem Mol Biol. 2010;40(3):, pp.214–27. , doi: 10.1016/j.ibmb.2010.01.011

35 

S Jasrapuria, CA Specht, KJ Kramer, RW Beeman, S Muthukrishnan. . Gene Families of Cuticular Proteins Analogous to Peritrophins (CPAPs) in Tribolium castaneum Have Diverse Functions. Plos One. 2012;7(11).

36 

B Wang, W Yang, J McKittrick, MA Meyers. . Keratin: Structure, mechanical properties, occurrence in biological organisms, and efforts at bioinspiration. Prog Mater Sci. 2016;76:, pp.229–318.

37 

M Suzuki, E Murayama, H Inoue, N Ozaki, H Tohse, T Kogure, et al. Characterization of Prismalin-14, a novel matrix protein from the prismatic layer of the Japanese pearl oyster (Pinctada fucata). Biochem J. 2004;382(Pt 1):, pp.205–13. , doi: 10.1042/BJ20040319

38 

Q Zhu, Y Arakane, D Banerjee, RW Beeman, KJ Kramer, S Muthukrishnan. . Domain organization and phylogenetic analysis of the chitinase-like family of proteins in three species of insects. Insect Biochem Mol Biol. 2008;38(4):, pp.452–66. , doi: 10.1016/j.ibmb.2007.06.010

39 

E Marchal, L Badisco, H Verlinden, T Vandersmissen, S Van Soest, P Van Wielendaele, et al. Role of the Halloween genes, Spook and Phantom in ecdysteroidogenesis in the desert locust, Schistocerca gregaria. J Insect Physiol. 2011;57(9):, pp.1240–8. , doi: 10.1016/j.jinsphys.2011.05.009

40 

H Ono, KF Rewitz, T Shinoda, K Itoyama, A Petryk, R Rybczynski, et al. Spook and Spookier code for stage-specific components of the ecdysone biosynthetic pathway in Diptera. Dev Biol. 2006;298(2):, pp.555–70. , doi: 10.1016/j.ydbio.2006.07.023

41 

AT Willingham, T Keil. . A tissue specific cytochrome P450 required for the structure and function of Drosophila sensory organs. Mech Develop. 2004;121(10):, pp.1289–97.

42 

L Wu, Q Jia, X Zhang, X Zhang, S Liu, Y Park, et al. CYP303A1 has a conserved function in adult eclosion in Locusta migratoria and Drosophila melanogaster. Insect Biochem Mol Biol. 2019;113:, pp.103210, doi: 10.1016/j.ibmb.2019.103210

43 

P Hurban, CS Thummel. . Isolation and characterization of fifteen ecdysone-inducible Drosophila genes reveal unexpected complexities in ecdysone regulation. Mol Cell Biol. 1993;13(11):, pp.7101–11. , doi: 10.1128/mcb.13.11.7101

44 

AR Frand, S Russel, G Ruvkun. . Functional genomic analysis of C. elegans molting. PLoS Biol. 2005;3(10):, pp.e312, doi: 10.1371/journal.pbio.0030312

45 

Y Arakane, B Li, S Muthukrishnan, RW Beeman, KJ Kramer, Y Park. . Functional analysis of four neuropeptides, EH, ETH, CCAP and bursicon, and their receptors in adult ecdysis behavior of the red flour beetle, Tribolium castaneum. Mech Develop. 2008;125(11–12):, pp.984–95.

46 

JD Baker, JW Truman. . Mutations in the Drosophila glycoprotein hormone receptor, rickets, eliminate neuropeptide-induced tanning and selectively block a stereotyped behavioral program. J Exp Biol. 2002;205(17):, pp.2555–65.

47 

CP Costa, M Elias-Neto, T Falcon, RP Dallacqua, JR Martins, MMG Bitondi. . RNAi-Mediated Functional Analysis of Bursicon Genes Related to Adult Cuticle Formation and Tanning in the Honeybee, Apis mellifera. Plos One. 2016;11(12).

48 

FQ Diao, BH White. . A Novel Approach for Directing Transgene Expression in Drosophila: T2A-Gal4 In-Frame Fusion. Genetics. 2012;190(3):, pp.1139–U356. , doi: 10.1534/genetics.111.136291

49 

CW Luo, EM Dewey, S Sudo, J Ewer, SY Hsu, HW Honegger, et al. Bursicon, the insect cuticle-hardening hormone, is a heterodimeric cystine knot protein that activates G protein-coupled receptor LGR2. P Natl Acad Sci USA. 2005;102(8):, pp.2820–5.

50 

J Xie, XX Hu, MF Zhai, XJ Yu, XW Song, SS Gao, et al. Characterization and functional analysis of hsp18.3 gene in the red flour beetle, Tribolium castaneum. Insect Sci. 2019;26(2):, pp.263–73. , doi: 10.1111/1744-7917.12543

51 

A Krogh, B Larsson, G von Heijne, EL Sonnhammer. . Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001;305(3):, pp.567–80. , doi: 10.1006/jmbi.2000.4315

52 

JD Thompson, TJ Gibson, DG Higgins. . Multiple sequence alignment using ClustalW and ClustalX. Curr Protoc Bioinformatics. 2002;Chapter 2:Unit 2 3.

53 

S Kumar, G Stecher, K Tamura. . MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol Biol Evol. 2016;33(7):, pp.1870–4. , doi: 10.1093/molbev/msw054

54 

TD Schmittgen, KJ Livak. . Analyzing real-time PCR data by the comparative C-T method. Nat Protoc. 2008;3(6):, pp.1101–8. , doi: 10.1038/nprot.2008.73

55 

X Song, F Huang, J Liu, C Li, S Gao, W Wu, et al. Genome-wide DNA methylomes from discrete developmental stages reveal the predominance of non-CpG methylation in Tribolium castaneum. DNA Res. 2017;24(5):, pp.445–57. , doi: 10.1093/dnares/dsx016

56 

C Li, Y Lu, S Ma, P Lu, B Li, K Chen. . Crinkled employs wingless pathway for wing development in Tribolium castaneum. Arch Insect Biochem Physiol. 2018;99(2):, pp.e21496, doi: 10.1002/arch.21496

57 

S Gao, X Liu, J Liu, W Xiong, X Song, W Wu, et al. Identification and evolution of latrophilin receptor gene involved in Tribolium castaneum devolopment and female fecundity. Genesis. 2017;55(12).

58 

J Zhang, KY Zhu. . Characterization of a chitin synthase cDNA and its increased mRNA level associated with decreased chitin synthesis in Anopheles quadrimaculatus exposed to diflubenzuron. Insect Biochem Mol Biol. 2006;36(9):, pp.712–25. , doi: 10.1016/j.ibmb.2006.06.002


Transfer Alert

This paper was transferred from another journal. As a result, its full editorial history (including decision letters, peer reviews and author responses) may not be present.

8 Jan 2020

Dear Dr Li,

Thank you very much for submitting your Research Article entitled 'A new neuropeptide insect parathyroid hormone iPTH in the red flour beetle Tribolium castaneum' to PLOS Genetics. Your manuscript was fully evaluated at the editorial level and by independent peer reviewers. The reviewers appreciated the attention to an important problem, but raised some substantial concerns about the current manuscript. Based on the reviews, we will not be able to accept this version of the manuscript, but we would be willing to review again a much-revised version. We cannot, of course, promise publication at that time.

Should you decide to revise the manuscript for further consideration here, your revisions should address the specific points made by each reviewer. We will also require a detailed list of your responses to the review comments and a description of the changes you have made in the manuscript.

If you decide to revise the manuscript for further consideration at PLOS Genetics, please aim to resubmit within the next 60 days, unless it will take extra time to address the concerns of the reviewers, in which case we would appreciate an expected resubmission date by email to plosgenetics@plos.org.

If present, accompanying reviewer attachments are included with this email; please notify the journal office if any appear to be missing. They will also be available for download from the link below. You can use this link to log into the system when you are ready to submit a revised version, having first consulted our Submission Checklist.

To enhance the reproducibility of your results, we recommend that you deposit your laboratory protocols in protocols.io, where a protocol can be assigned its own identifier (DOI) such that it can be cited independently in the future. For instructions see our guidelines.

Please be aware that our data availability policy requires that all numerical data underlying graphs or summary statistics are included with the submission, and you will need to provide this upon resubmission if not already present. In addition, we do not permit the inclusion of phrases such as "data not shown" or "unpublished results" in manuscripts. All points should be backed up by data provided with the submission.

While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool.  PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email us at figures@plos.org.

PLOS has incorporated Similarity Check, powered by iThenticate, into its journal-wide submission system in order to screen submitted content for originality before publication. Each PLOS journal undertakes screening on a proportion of submitted articles. You will be contacted if needed following the screening process.

To resubmit, use the link below and 'Revise Submission' in the 'Submissions Needing Revision' folder.

[LINK]

We are sorry that we cannot be more positive about your manuscript at this stage. Please do not hesitate to contact us if you have any concerns or questions.

Yours sincerely,

Liliane Schoofs

Associate Editor

PLOS Genetics

Gregory Barsh

Editor-in-Chief

PLOS Genetics

Reviewer's Responses to Questions

Comments to the Authors:

Please note here if the review is uploaded as an attachment.

Reviewer #1: Comments:

Introduction

line [87] Add other references of relevant insect GPCR papers/reviews?

Results

[129] Change (Fig. 2A) to (Fig. 2A; Fig. S1). Figure S1 displays the gene duplication in insects more clearly.

[136] According to the supporting information file, Fig S1 only displays PTHRs and no iPTH sequences. Thus, there is either a reference to a Figure missing or a reference to literature.

[149-151] This is described in literature as characteristic of B-group GPCRs. Add this fact in your description/discussion, together with relevant references.

[166-167] Remove line (see the definition of early and late in the materials and methods), as it speaks for itself that the description would and should be in this section.

Discussion

[346-349] Agreed, this is a very interesting finding. However, how does this correlate with the neural projections of iPTH (Fig. 4) and the expression levels of both iPTH and both iPTHR (as described in Fig. 3 & 5)? Both receptors and iPTH itself seem to be mainly expressed in Gut and CNS (Fig 3.). There should be at least mention of this fact in your discussion, even if the effect would probably be indirect.

[354-355] You have not really proven that the receptors are unable to compensate each other, see comment line 460-462. Suggested experiment: RT-qPCR of both receptors in experimental tissues used of these experiments. As you already have the samples, this addition should be easily conducted. These results must be added in supplementary data. Or alternatively, show this fact explicitly in a figure separate from of your RNAseq data.

Materials and methods

[439] Define “multiple individuals” as this could describe any quantity, from 2 to infinity. The reader should be informed about the size of the biological pool by which one can determine the relevance of the described results.

[460-462] This potential has not been experimentally excluded, or it is not shown. One should. What kind of controls or additional experiments have been performed to confirm the data and to show that these are not due to off-target effects?

[466-467] This RNAi efficiency at RNA/protein level is not shown. Add to Supplementary data.

Reviewer #2: The review is uploaded as an attachment.

Reviewer #3: This paper describes the first identification of insect peptides that activate “orphan” receptors (iPTHRs) homologous to the vertebrate parathyroid hormone receptor (PTHR). On this basis, the authors designate these peptides “insect parathyroid hormones” (iPTHs). Interestingly, while vertebrate PTH regulates calcium homeostasis by modulating bone metabolism, iPTHs are implicated in regulation of exoskeleton maturation and fecundity. This is a major advance in assigning functional roles for iPTHRs and likely will further our understanding of ligand-receptor coevolution. The manuscript summarizes a great deal of work and is well presented. The bioinformatic analysis is excellent and RNAi knockdown creates phenotypes that involve integumentary deficits. Some questions remain should be addressed in the interests of clarifying the source(s) and targets of iPTHs.

Specific Comments

1. Designation of peptides as iPTHs. On one hand, it is understandable, since functions of iPTHs and vertebrate PTHs appear to be quite different. But why not simply refer to them according to species (e.g., TcPTH)? Perhaps the authors could explain the rationale for the designation of iPTH.

2. Sequences shown in Fig. E 1 fail to depict amidation at the C-terminus; it might be an idea to add a lower case “a” for amidation at the C-terminus of each sequence.

3. No mention of sequence similarity between iPTHs and vertebrate PTH. This reviewer found 47% similarity and 31% identity in a pairwise comparison with the consensus sequence shown in Fig. 1C. Addition of such information would complement analysis of sequence similarities of iPTHR and vertebrate PTHR.

4. Is there any evidence that iPTHs regulate hemolymph Ca levels, as is the case for vertebrate PTH?

5. Unlike the tissue source of vertebrate PTH, which is a well-defined gland behind the thyroid gland, iPTHs apparently originate from diverse sources (CNS, gut, etc). This is a matter of some uncertainty, however, since the antiserum used to locate iPTH-positive cells is not well characterized (see below). Hence, the possibility of cross-reactivity with other epitopes cannot be ruled out. Control experiments reported here are confined to negative results with pre-immune serum. This is not convincing. Additional controls could be: 1) determine if staining intensity decreases following RNAi of ligands and receptors. Although this has been shown to some extent in Fig. S4 for iPTH, convincing side-by-side comparisons with and without dsRNA are not provided. Another control could be pre-absorption of antiserum with ligand, following by staining, which should be negative. Likewise, no characterization of antisera raised against iPTHR is provided. Difficulties in specific immunohistochemical identification of GPCRs are well-known. Finally, it would be much more convincing if corroborating in situ hybridizations were included to demonstrate presence of the predicted transcripts for iPTHs and iPTHRs.

6. A second problem with the immunohistochemical data shown in Figs 4 and 5 is that the images are not described in enough detail. It is often difficult to discern which parts of the nervous system are depicted in Figs. They are simply referred to en masse; consequently, the reader cannot discern anything specific.

7. Any indication that iPTH is secreted by endocrine cells or glands, as is vertebrate PTH?

8. With regard to functionality of iPTH signaling, no mention is made in the abstract about reduction of fecundity following RNAi treatment.

9. Acronyms are not adequately defined in many instances; e.g. CPR family (line 325), DE analyses (line 512).

10. Abstract, line 39-40. Reference formatting is inconsistent and it is unusual to cite references in the abstract.

11. Lines 158-159. Fig 2B contains no information to support this statement of negative results with untransfected cells.

12. The manuscript requires further editing to correct English grammar and spelling.

**********

Have all data underlying the figures and results presented in the manuscript been provided?

Large-scale datasets should be made available via a public repository as described in the PLOS Geneticsdata availability policy, and numerical data that underlies graphs or summary statistics should be provided in spreadsheet form as supporting information.

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Yes

**********

PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: No

Reviewer #2: No

Reviewer #3: No

Submitted filename: review_Xie_plosgenet2019_rev_Mirabeau.docx

16 Feb 2020

Submitted filename: Response to Reviewers.docx

20 Mar 2020

* Please note while forming your response, if your article is accepted, you may have the opportunity to make the peer review history publicly available. The record will include editor decision letters (with reviews) and your responses to reviewer comments. If eligible, we will contact you to opt in or out. *

Dear Dr Li,

Thank you very much for submitting your Research Article entitled 'A new neuropeptide insect parathyroid hormone iPTH in the red flour beetle Tribolium castaneum' to PLOS Genetics. Your manuscript was fully evaluated at the editorial level and by independent peer reviewers. The reviewers appreciated the attention to an important topic but still identified some aspects of the manuscript that should be improved.

We therefore ask you to modify the manuscript according to the review recommendations before we can consider your manuscript for acceptance. Your revisions should address the specific points made by each reviewer.

In addition we ask that you:

1) Provide a detailed list of your responses to the review comments and a description of the changes you have made in the manuscript.

2) Upload a Striking Image with a corresponding caption to accompany your manuscript if one is available (either a new image or an existing one from within your manuscript). If this image is judged to be suitable, it may be featured on our website. Images should ideally be high resolution, eye-catching, single panel square images. For examples, please browse our archive. If your image is from someone other than yourself, please ensure that the artist has read and agreed to the terms and conditions of the Creative Commons Attribution License. Note: we cannot publish copyrighted images.

We hope to receive your revised manuscript within the next 30 days. If you anticipate any delay in its return, we would ask you to let us know the expected resubmission date by email to plosgenetics@plos.org.

If present, accompanying reviewer attachments should be included with this email; please notify the journal office if any appear to be missing. They will also be available for download from the link below. You can use this link to log into the system when you are ready to submit a revised version, having first consulted our Submission Checklist.

While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email us at figures@plos.org.

Please be aware that our data availability policy requires that all numerical data underlying graphs or summary statistics are included with the submission, and you will need to provide this upon resubmission if not already present. In addition, we do not permit the inclusion of phrases such as "data not shown" or "unpublished results" in manuscripts. All points should be backed up by data provided with the submission.

PLOS has incorporated Similarity Check, powered by iThenticate, into its journal-wide submission system in order to screen submitted content for originality before publication. Each PLOS journal undertakes screening on a proportion of submitted articles. You will be contacted if needed following the screening process.

To resubmit, you will need to go to the link below and 'Revise Submission' in the 'Submissions Needing Revision' folder.

[LINK]

Please let us know if you have any questions while making these revisions.

Yours sincerely,

Liliane Schoofs

Associate Editor

PLOS Genetics

Gregory Barsh

Editor-in-Chief

PLOS Genetics

Reviewer's Responses to Questions

Comments to the Authors:

Please note here if the review is uploaded as an attachment.

Reviewer #1: Several comments have been well addressed, but some are remaining:

Introduction

Line[97-98]: Described references are not added in the text, nor added in reference list. Please do so, check more carefully …

Results

Line[161] Add “,” between … binding pocket - and - as stabilization…

Discussion:

[441-445]: sentences are incoherent, revise wording.

More in depth: L442: Too vague: high expression of what? RNA, protein? Be more specific.

L443: Too vague: which additional functions? This does not add anything. Either describe which additional functions you imply, or don’t mention this.

[354-355]: Figure shows exactly what was asked for, but it should be shown to readers in order to make these claims. Add this figure either to Fig 7 or Fig 8 as an extra panel. Or if not shown, than the claim that the receptors are unable to compensate each other should be eliminated from your text.

Reviewer #2: I thank the authors for taking into account my remarks which I consider have been sufficiently addressed.

Some of the newly written sentences could be improved:

lines 58-59: Vertebrate parathyroid hormone (PTH) and its receptors have been extensively studied with RESPECT TO [the] THEIR function [of] IN bone remodeling and calcium metabolism.

lines 53-54: Moreover, RNAi of iPTHRs also led TO significant reductions in egg numbers and hatching rates after parental RNAi.

line 59-61: Insect parathyroid hormone receptors (iPTHRs) [were] HAVE BEEN previously described as counterparts of vertebrate PTHRs.

line 82-83: Functional studies of ancestral bilaterian neuropeptides have been successful by starting from the sequence similarities, e.g., vasopressin AND TRH [1-4]

line150: coincided with the [lack] ABSENCE of this specific

line 152: Maximum likelihood method with THE JTT model

line 220: neuronal cells and axonal projections in THE CNS

line 284: beetles arrested at the time of [ecolsion] ECLOSION

line 297: suppression of Tc-iPTH by the RNAi may [have not] NOT HAVE been sufficient to induce

lines 306-307: plots also showed THAT DE genes in the dsiPTHR1 and dsiPTHR2 treatments [to] have similar patterns

line 310: the correlation between dsiPTHR1 and dsiPTHR2 treatments [with a R-square showing 0.78] was obvious (R-SQUARED=0.78)

line 430: The CONCLUSIONS DRAWN FROM THE DE data ANALYSIS [together] were also supported by cuticle contents

line 460: [yet] STILL displayed [remaining] RESIDUAL iPTH immunoreactivity

line 472: THAT THIS NEUROPEPTIDE MAY HAVE additional functions.

line 805: "Two-dimensional scatter plots" instead of "Two dimensional scatter plots"

**********

Have all data underlying the figures and results presented in the manuscript been provided?

Large-scale datasets should be made available via a public repository as described in the PLOS Geneticsdata availability policy, and numerical data that underlies graphs or summary statistics should be provided in spreadsheet form as supporting information.

Reviewer #1: No: See comments to authors

Reviewer #2: Yes

**********

PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: No

Reviewer #2: No


23 Mar 2020

Submitted filename: Response to reviewers.docx

9 Apr 2020

Dear Dr Li,

We are pleased to inform you that your manuscript entitled "A new neuropeptide insect parathyroid hormone iPTH in the red flour beetle Tribolium castaneum" has been editorially accepted for publication in PLOS Genetics. Congratulations!

Before your submission can be formally accepted and sent to production you will need to complete our formatting changes, which you will receive in a follow up email. Please be aware that it may take several days for you to receive this email; during this time no action is required by you. Please note: the accept date on your published article will reflect the date of this provisional accept, but your manuscript will not be scheduled for publication until the required changes have been made.

Once your paper is formally accepted, an uncorrected proof of your manuscript will be published online ahead of the final version, unless you’ve already opted out via the online submission form. If, for any reason, you do not want an earlier version of your manuscript published online or are unsure if you have already indicated as such, please let the journal staff know immediately at plosgenetics@plos.org.

In the meantime, please log into Editorial Manager at https://www.editorialmanager.com/pgenetics/, click the "Update My Information" link at the top of the page, and update your user information to ensure an efficient production and billing process. Note that PLOS requires an ORCID iD for all corresponding authors. Therefore, please ensure that you have an ORCID iD and that it is validated in Editorial Manager. To do this, go to ‘Update my Information’ (in the upper left-hand corner of the main menu), and click on the Fetch/Validate link next to the ORCID field.  This will take you to the ORCID site and allow you to create a new iD or authenticate a pre-existing iD in Editorial Manager.

If you have a press-related query, or would like to know about one way to make your underlying data available (as you will be aware, this is required for publication), please see the end of this email. If your institution or institutions have a press office, please notify them about your upcoming article at this point, to enable them to help maximise its impact. Inform journal staff as soon as possible if you are preparing a press release for your article and need a publication date.

Thank you again for supporting open-access publishing; we are looking forward to publishing your work in PLOS Genetics!

Yours sincerely,

Liliane Schoofs

Associate Editor

PLOS Genetics

Gregory Barsh

Editor-in-Chief

PLOS Genetics

www.plosgenetics.org

Twitter: @PLOSGenetics

----------------------------------------------------

Comments from the reviewers (if applicable):

----------------------------------------------------

Data Deposition

If you have submitted a Research Article or Front Matter that has associated data that are not suitable for deposition in a subject-specific public repository (such as GenBank or ArrayExpress), one way to make that data available is to deposit it in the Dryad Digital Repository. As you may recall, we ask all authors to agree to make data available; this is one way to achieve that. A full list of recommended repositories can be found on our website.

The following link will take you to the Dryad record for your article, so you won't have to re‐enter its bibliographic information, and can upload your files directly: 

http://datadryad.org/submit?journalID=pgenetics&manu=PGENETICS-D-19-01925R2

More information about depositing data in Dryad is available at http://www.datadryad.org/depositing. If you experience any difficulties in submitting your data, please contact help@datadryad.org for support.

Additionally, please be aware that our data availability policy requires that all numerical data underlying display items are included with the submission, and you will need to provide this before we can formally accept your manuscript, if not already present.

----------------------------------------------------

Press Queries

If you or your institution will be preparing press materials for this manuscript, or if you need to know your paper's publication date for media purposes, please inform the journal staff as soon as possible so that your submission can be scheduled accordingly. Your manuscript will remain under a strict press embargo until the publication date and time. This means an early version of your manuscript will not be published ahead of your final version. PLOS Genetics may also choose to issue a press release for your article. If there's anything the journal should know or you'd like more information, please get in touch via plosgenetics@plos.org.


24 Apr 2020

PGENETICS-D-19-01925R2

A new neuropeptide insect parathyroid hormone iPTH in the red flour beetle Tribolium castaneum

Dear Dr Li,

We are pleased to inform you that your manuscript entitled "A new neuropeptide insect parathyroid hormone iPTH in the red flour beetle Tribolium castaneum" has been formally accepted for publication in PLOS Genetics! Your manuscript is now with our production department and you will be notified of the publication date in due course.

The corresponding author will soon be receiving a typeset proof for review, to ensure errors have not been introduced during production. Please review the PDF proof of your manuscript carefully, as this is the last chance to correct any errors. Please note that major changes, or those which affect the scientific understanding of the work, will likely cause delays to the publication date of your manuscript.

Soon after your final files are uploaded, unless you have opted out or your manuscript is a front-matter piece, the early version of your manuscript will be published online. The date of the early version will be your article's publication date. The final article will be published to the same URL, and all versions of the paper will be accessible to readers.

Thank you again for supporting PLOS Genetics and open-access publishing. We are looking forward to publishing your work!

With kind regards,

Matt Lyles

PLOS Genetics

On behalf of:

The PLOS Genetics Team

Carlyle House, Carlyle Road, Cambridge CB4 3DN | United Kingdom

plosgenetics@plos.org | +44 (0) 1223-442823

plosgenetics.org | Twitter: @PLOSGenetics

https://www.researchpad.co/tools/openurl?pubtype=article&doi=10.1371/journal.pgen.1008772&title=A new neuropeptide insect parathyroid hormone iPTH in the red flour beetle <i>Tribolium castaneum</i>&author=Jia Xie,Ming Sang,Xiaowen Song,Sisi Zhang,Donghun Kim,Jan A. Veenstra,Yoonseong Park,Bin Li,Liliane Schoofs,Gregory S. Barsh,Liliane Schoofs,Gregory S. Barsh,Liliane Schoofs,Gregory S. Barsh,Liliane Schoofs,Gregory S. Barsh,Liliane Schoofs,&keyword=&subject=Research Article,Biology and life sciences,Genetics,Epigenetics,RNA interference,Biology and life sciences,Genetics,Gene expression,RNA interference,Biology and life sciences,Genetics,Genetic interference,RNA interference,Biology and life sciences,Biochemistry,Nucleic acids,RNA,RNA interference,Biology and Life Sciences,Organisms,Eukaryota,Animals,Invertebrates,Arthropoda,Insects,Biology and Life Sciences,Biochemistry,Hormones,Parathyroid Hormone,Biology and Life Sciences,Genetics,Gene Expression,Research and Analysis Methods,Database and Informatics Methods,Bioinformatics,Sequence Analysis,Sequence Alignment,Biology and Life Sciences,Anatomy,Nervous System,Central Nervous System,Medicine and Health Sciences,Anatomy,Nervous System,Central Nervous System,Biology and Life Sciences,Organisms,Eukaryota,Animals,Vertebrates,Physical Sciences,Chemistry,Polymer Chemistry,Macromolecules,Polymers,Chitin,Physical Sciences,Materials Science,Materials,Polymers,Chitin,Physical Sciences,Chemistry,Polymer Chemistry,Polymers,Chitin,