ORIGINAL RESEARCH article

Front. Physiol., 20 May 2020

Sec. Invertebrate Physiology

Volume 11 - 2020 | https://doi.org/10.3389/fphys.2020.00424

Molecular Disruption of Ion Transport Peptide Receptor Results in Impaired Water Homeostasis and Developmental Defects in Bombyx mori

  • 1. College of Forestry, Northeast Forestry University, Harbin, China

  • 2. CAS Center for Excellence in Biotic Interactions, University of Chinese Academy of Sciences, Beijing, China

  • 3. Key Laboratory of Insect Developmental and Evolutionary Biology, CAS Center for Excellence in Molecular Plant Sciences, Shanghai Institute of Plant Physiology and Ecology, Chinese Academy of Sciences, Shanghai, China

Abstract

Insect ion transport peptides (ITPs) are important regulators of many physiological processes and they exert their functions by interacting with their receptors (ITPRs). In the current study, we comprehensively investigated the physiological functions of ITPR in the lepidopteran model insect, the silkworm (Bombyx mori), using the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein-9 nuclease (Cas9) genome editing technique. Mutations in silkworm ITPR (BNGR-A2) resulted in a prolongnation of the larval stage by 3.5-day as well as failure in wing expansion of moths. The BNGR-A2 mutation accelerated food transition throughout the digestive tract, which is 1.55-fold that of wild type (WT) insects. Excretion was 1.56-fold of WT insects during the larval stage, resulting in the loss of body water content. Loss of BNGR-A2 function induced significant upregulation of nitric oxide synthase (NOS) enzyme activity and nitric oxide (NO) content, as well as downstream Ca2+/NO/cGMP signaling pathways. Key genes in insulin and ecdysone signaling pathways were also affected by BNGR-A2 disruption. Our data show that ITPR plays key roles in regulating insect water homeostasis and development.

Introduction

Ion transport peptides (ITPs) and its alternatively spliced homologous ITP-like (ITPL) products belong to the crustacean hyperglycemic hormone (CHH) family of peptides. They are homologs of the crustacean hyperglycemic hormones of molting-inhibiting hormone (MIH) and gonad/vitellogenesis-inhibiting hormone (GIH) (Audsley et al., 2006; Webster et al., 2012). In crustaceans, CHH family peptides play roles in energy metabolism, molting, reproduction, immune defense, and osmotic regulation, as well as in homeostatic regulation of stress responses (Sonobe et al., 2001; Webster et al., 2012). In insects, ITP and ITPL with Cl, Na+, and Ca2+ transport, as well as fluid reabsorption functions, have been identified (Audsley et al., 1992; Drexler et al., 2007; Dircksen et al., 2008; Webster et al., 2012; Nagai et al., 2014; Yu et al., 2016). In Drosophila melanogaster, ITP plays a critical role in development and also participates in the control of locomotor rhythms (Johard et al., 2009; Hermann-Luibl et al., 2014; Gáliková et al., 2018), and the mutation of ITP is embryonically lethal (Park et al., 2004). Gáliková et al. (2018) work identified master regulatory roles of ITP in water homeostasis of Drosophila; ITP levels increase under desiccation stress and protect the fly from water loss by increasing thirst, reducing excretion rate, and promoting ingestion of water instead of food. Several studies have shown that ITPs function during ecdysis in Manduca sexta (Drexler et al., 2007). The insect ITPL also participates in ovarian maturation in Tribolium castaneum (Begum et al., 2009) and regulates wing expansion and cuticle melanism in Nilaparvata lugens (Yu et al., 2016).

Insect ITP and ITPL act through their receptors (ITPRs) to function in many physiological processes. However, no ITPRs have been identified in the model insect D. melanogaster thus far (Gáliková et al., 2018; Nässel and Zandawala, 2019). In the lepidopteran model insect, Bombyx mori, it has been reported that two Bombyx neuropeptide G protein-coupled receptors (BNGRs), BNGR-A2 and BNGR-A34, are native receptors for ITP (Nagai et al., 2014). The ITP receptor, BNGR-A34, is conserved across insect species, whereas BNGR-A2 may have a species-specific role in lepidopterans (Nagai et al., 2014). Although a positive correlation between response of the receptor and bioactivity of the ligand has been verified in vitro (Nagai et al., 2014), the physiological functions of both ITP/ITPR and their receptors in insects in vivo remain to be elucidated.

In this study, we investigated the physiological functions of BNGR-A2 using the CRISPR/Cas9 genome editing technique. Loss of BNGR-A2 function extended the larval developmental stage and caused failure of wing expansion in moths. The BNGR-A2 mutation also accelerated food transition throughout the digestive tract as well as excretion. Our data therefore reveals for the first time that ITPR plays critical roles in regulating water homeostasis and developement in B. mori.

Materials and Methods

Silkworm Strains

A multivoltine silkworm strain, Nistari, was used in all experiments. Larvae were reared with fresh mulberry leaves at 25°C and 75% relative humidity (RH) (Tan et al., 2013).

Plasmid Construction and Germline Transformation

The activator line of pBac [IE1-EGFP-Nos-Cas9] (Nos-Cas9) was used to construct the transgenic CRISPR/Cas9 system. Cas9 was driven by the germ-cell-specific promoter, Nos, as described previously (Xu et al., 2017b). The effector line pBac [IE1-DsRed × 2-U6-2 × BNGR-A2-sgRNAs] (BNGR-A2-sgRNAs) was under the control of the silkworm small nuclear RNA promoter U6. The plasmids targeting BNGR-A2 were constructed through a series of cloning steps (Xu et al., 2017a); primer sequences are listed in Supplementary Table S1.

For silkworm germline transformation, G0 Nistari embryos within 8 h of oviposition were separately injected with a plasmid mixture of Nos-Cas9 or BNGR-A2-sgRNAs. Hatched larvae were reared to the adult stage and G0 moths were sib-mated or backcrossed with wild-type (WT) moths. Using a fluorescence microscope (Nikon AZ100), G1 offspring were screened for the marker gene during the embryonic stage. The Nos-Cas9 line and the BNGR-A2-sgRNA line were crossed to produce heterozygous F1 progeny (mutant was obtained by indicator of two-color fluorescence), which was then used in subsequent experiments.

Genomic DNA Extraction and Mutagenesis Analysis

The BNGR-A2 mutant genomic DNA was extracted using a nucleic acid isolation system DP323 (TianGen, Beijing, China), and the DNA fragments, including the designed sgRNA targeting site, were amplified using EasyTaq® DNA Polymerase (TransGen Biotech, Beijing, China). The PCR conditions were as follows: 94°C for 5 min, 35 cycles at 94°C for 30 s, 60°C for 30 s, and 72°C for 1 min, followed by a final extension period of 72°C for 10 min. The PCR products were cloned into pMD19-T-Simple (Takara, Dalian, China) and subsequently sequenced. The primers that were designed to detect mutagenesis in targeted sites are listed in Supplementary Table S1.

Bioassays

The cumulative numbers of fecal pellets from 20 fifth instar silkworms (0 day of fifth instar larvae, L5D0) were counted in the breeding box every hour. Droppings were counted in three replicates. The statistical significance was analyzed by two-way ANOVA, and the number of fresh fecal masses that were deposited within 2 h (the given period) was used.

To measure water content, L5D0 silkworm larvae were dehydrated at 80°C until a constant weight. The weight of 10 larvae was recorded before and after dehyration using a Mettler MT5 analytical microbalance (Columbus, OH, United States). The difference in water content between the fresh and dry weights was calculated and is expressed as the percent of the fresh body weight. Three replicates (each consisting of 10 silkworms) were tested per WT and mutant silkworm.

Ecdysteroid Titer and NOS Activity Determination

Haemolymph was extracted from transgenic or WT larvae (L5D3). Total ecdysteroids in the haemolymph samples were extracted as previously described (Mirth et al., 2005). Concentrations of total ecdysteroids were quantified using a 20-Hydroxyecdysone EIA kit (Cayman Chemicals, Ann Arbor, MI, United States) according to the manufacturer’s instructions.

The nitric oxide synthase (NOS) kit (Nanjing Jiancheng Bioengineering Institute, Nanjing, China) was used to detect the total NOS activity in silkworm hindgut. The total nitric oxide (NO) was determined using a Nitric Oxide Assay Kit (Beyotime, Shanghai, China). Protein concentrations were measured using the Bradford Protein Assay Kit (Beyotime, Shanghai, China). Each measurement was from at least four biological replicates including two males and two females. Data were normalized to the protein concentration.

RNA Isolation, cDNA Synthesis, and Quantitative Real-Time PCR (qRT-PCR) Analysis

Total RNA was extracted from the fat body and hindgut of L5D3 animals using Trizol reagent (Invitrogen, Carlsbad, CA, United States) in accordance with the manufacturer’s instructions. Then, 0.5 μg of RNA was reverse-transcribed for cDNA using TransScript® RT/RI Enzyme Mix (TransGen Biotech, Beijing, China). qRT-PCR was performed using a SYBR Green PCR Kit (Toyobo, Osaka, Japan) and an ABI Stepone plus real-time PCR system (Applied Biosystems, United States). Each 20 μL reaction contained 10 μL of SYBR Green real-time PCR master mix (Toyobo, Osaka, Japan), 0.5 μM of each forward and reverse primer, and 2 μL of cDNA template (equivalent to 100 ng of total RNA). Primers used for amplification are listed in Supplementary Table S1. Thermal cycling parameters were 94°C for 5 min followed by 40 cycles at 94°C for 30 s and 60°C for 1 min. For each sample, a melting curve was generated at the end of each run to allow assessment of product purity. Expression levels were calculated according to the 2–ΔΔCT method (Pfaffl et al., 2002), and relative expression levels were calculated by dividing mutant transcription levels by levels recorded in wild type, WT insects.

RNA-Seq Analysis

Hindguts from WT and mutant L5D3 animals were separately collected and used to construct a transcriptome library. Each type of sample was collected from four animals (two males and two females). The cDNA libraries were sequenced using an Illumina BGISEQ-500 platform (BGI, Wuhan, China). Total RNA was isolated as described previously, and cDNA libraries were prepared according to Illumina’s protocol. The raw data were filtered using Tophat software, corrected, and mapped to the silkworm genome sequence. The total reading was normalized by a plurality of normalization factors. Transcript levels were calculated using standard readings per petaflop of mapping readings. The difference between the mutant and the control groups is represented by the P value and fold change (ΔBNGR-A2/WT). Differentially expressed genes with a significance level of P < 0.001 and fold change (ΔBNGR-A2/WT) > 2 were enriched in each comparison.

Statistical Analyses

Two-tailed Student’s t-test and one-way ANOVA were used to analyze measurement variables. Asterisk symbols reflect significant differences (P < 0.05, ∗∗P < 0.01, ∗∗∗P < 0.001, ****P < 0.0001). Error bars represent SEM. Measurement variable data was analyzed using Excel Workbook (version 2016) and data related to nominal variables were analyzed by GraphPad Prism (version 7.0).

Results

Generation of BNGR-A2 Targeted Mutations Using the Transgenic CRISPR/Cas9 System

A binary transgenic CRISPR/Cas9 system (Li et al., 2015) was used to generate somatic BNGR-A2 mutants. This system included two transgenic silkworm lines. One line expressed Cas9 under the control of Nos promoter and the other line expressed two sequence-specific sgRNAs under the control of U6 promoter to target BNGR-A2 exon 1 (Figures 1A,B). Transgenic animals carrying nos-Cas9 or U6-BNGR-A2–sgRNA were fully viable and fertile, indicating that the accumulation of neither Cas9 nor BNGR-A2-sgRNA had deleterious effects on silkworm physiology. The offspring crossed by Nos-Cas9 and U6-BNGR-A2-sgRNA lines specifically expressed the active Cas9-sgRNA complex in the germ line. The mutations were located at the targeted genomic loci of the BNGR-A2 gene (Figure 1B). To assess mutation efficiency, five newly hatched silkworms with both red and green fluoresence were randomly selected to extract the genomic DNA. The DNA samples were run through PCR amplification. PCR-amplified fragments were subcloned and sequenced. The results showed that various deletion mutations located at two sgRNA targets and deletions ranged from 215 to 349 bp, indicating that successful mutations were produced using the transgenic CRISPR/Cas9 system. Moreover, the BNGR-A2 gene was successfully disrupted (Figure 1C).

FIGURE 1

Water Content and Excretion Rate of BNGR-A2 Mutant Larvae

Previous studies showed that ITP has a positive effect on the water balance in Drosophila (Gáliková et al., 2018). However, whether ITPR affects the water balance in insects is unclear. In order to investigate whether BNGR-A2 is involved in water homeostasis in B. mori, the water contents of WT and mutant 5th instar larvae were determined. The results showed that the water contents of ΔBNGR-A2 animals was (87.34 ± 0.03)%, which was significantly lower than that of WT animals (88.76 ± 0.05)%; P < 0.0001), indicating that the loss of BNGR-A2 function affected the water homeostasis regulation in B. mori (Figure 2A).

FIGURE 2

In insects, the primary urete enters the hindgut, and water that is not reabsorbed by the posterior intestinal epithelium is usually excreted in the same way as feces (Gäde, 2004). We investigated whether BNGR-A2 knockout affects excretion. The excretion rate of the 5th instar larvae was calculated by counting the number of fecal particles produced every 2 h. The results showed that each mutant larvae produced 29 fecal particles per 2 h period, which was significantly higher than in WT animals (18.67 fecal particles per 2 h, P < 0.05). Loss of BNGR-A2 function accelerated food transition throughout the digestive tract (Figure 2B), resulting in more frequent defecation events.

Phenotypic Defects in BNGR-A2 Mutants

The BNGR-A2 mutant has no deleterious phenotype during embryonic development. From the end of the 2nd larval instar, mutant larvae exhibited a significant delay in growth rate relative to WT animals (Figures 3A,B). Compared with WT animals, the development of mutant larvae was delayed by half a day during the 3rd larval stage. The weight of mutant larvae was 57.41% of WT animals at L3D2. During the 4th instar larval stage, the weight of the mutant larvae was 55.58% that of WT animals, and development was delayed by 2.5 days. The weight of the 5th instar mutant larvae was 73.84% of the WT animals at L5D0, and the weight of the L5D4 mutant was similar to that of the L5D2 WT animals (Figures 3A,B). The half pupation time (PT50) of the mutant silkworm was delayed by 3 days relative to WT silkworms (Figure 3A). Pupal weight did not show significant differences between mutant and WT animals. After eclosion, the fore and hindwings of mutant moths became crumpled (Figure 3C).

FIGURE 3

Ca2+/NO/cGMP Signaling in BNGR-A2 Mutants

To explore the molecular mechanisms underlying defective phenotypes induced by BNGR-A2 mutations, we perfomed RNA-seq analysis by using hindgut tissues between WT and mutant larvae. A total of 1,288 differentially expressed genes (612 up-regulated and 676 down-regulated) were identified (Figure 4A and Supplementary Figure 1SA). The up- and down-regulated genes were allocated into three categories according to the Gene Ontology (GO) terms: molecular function, cellular component, and biological process (Supplementary Figure 1SB). Among GO categories, the “catalytic activity” in “molecular function,” “membrane” and “membrane part” in “cellular component,” and “metabolic process” in “biological process” were the most abundant protein classes. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that water-transport-related pathways, such as the calcium signaling pathway, cGMP-PKG signaling pathway and cAMP signaling pathway, were annotated (Table 1 and Supplementary Figure 1SC).

FIGURE 4

TABLE 1

Signaling pathwayFunction annotationP valueGene IDlog2 (ΔBNGR-A2/WT)Nr
GPCR0NGR-A243.36NP_001127722.1| 3.9e-225| neuropeptide receptor A24 [Bombyx mori]
CALM0LOC1017428902.94XP_004927631.1| 4.9e-135| neo-calmodulin [Bombyx mori]
GPCR3.24E-21NGR-A171.66NP_001127715.1| 7.0e-204| neuropeptide receptor A17 [Bombyx mori]
NOS5.22E-26NOS21.44XP_012545275.1| 0.0e+00| nitric oxide synthase-like protein isoform X1 [Bombyx mori]
GPCR3.45E-07NGR-A321.17NP_001127748.1| 1.2e-215| neuropeptide receptor A32 [Bombyx mori]
GPCR4.10E-131LOC101736629–1.04XP 022130216.1| 0.0e+00| cytosolic carboxypeptidase 4-like isoform X6 [Pieris rapae]
CalciumROC6.28E-16BGI novel G000814–1.23
signaling pathwaySERCA2.54E-13LOC101735353–1.28XP_021207909.1| 6.5e-220| uncharacterized protein LOC101735353 [Bombyx mori]
CALM7.40E-08LOC101741700–1.28XP_004922728.1| 1.3e-163| sarcoplasmic calcium-binding protein [Bombyx mori]
SPHK1.76E-90LOC101742466–1.31XP_004931675.2| 0.0e+00| sphingosine kinase 1 [Bombyx mori]
GPCR1.82E-05LOC101737969–1.55XP_012551041.1| 3.2e-262| beta-3 adrenergic receptor [Bombyx mori]
IP3 3K2.91E-192LOC101739177–1.83XP_012547652.1| 1.2e-260| uncharacterized protein LOC101739177 [Bombyx mori]
PLC β0.000261562LOC101743946–2.50XP_021206545.1| 0.0e+00| uncharacterized protein LOC101743946 isoform X1 [Bombyx mori]
GPCR1.07E-22LOC101743833–2.94XP_021204873.1| 0.0e+00| cytosolic carboxypeptidase 2 [Bombyx mori]
cGMP-PKGCAM0LOC1017428902.94XP_004927631.1| 4.9e-135| neo-calmodulin [Bombyx mori]
signaling pathwayATPase1.03E-108LOC1017460121.21XP_004929151.1| 1.3e-193| sodium/potassium-transporting ATPase subunit beta-2 [Bombyx mori]
INS3.57E-08BmILP1.15XP_012548888.1| 1.1e-66| insulin-like peptide isoform X1 [Bombyx mori]
ROCK6.10E-07BGI novel G000435–1.01
ROCK0.000248524BGI novel G000569–1.09
CREB6.19E-07BGI novel G000948–1.19
SERCA2.54E-13LOC101735353–1.28XP_021207909.1| 6.5e-220| uncharacterized protein LOC101735353 [Bombyx mori]
CAM7.40E-08LOC101741700–1.28XP_004922728.1| 1.3e-163| sarcoplasmic calcium-binding protein [Bombyx mori]
ROCK9.43E-05BGI novel G000712–1.32
ROCK2.33E-07LOC105843016–1.43XP_011555716.1| 1.9e-112| PREDICTED: uncharacterized protein LOC105386785 [Plutella xylostella]
SRF7.66E-32LOC101744587–1.62XP_012552250.1| 8.2e-88| serum response factor homolog [Bombyx mori]
GATA-46.67E-33LOC101737855–2.36XP_004924467.1| 0.0e+00| transcription factor GATA-6 isoform X1 [Bombyx mori]
PLC β0.000261562LOC101743946–2.50XP_021206545.1| 0.0e+00| uncharacterized protein LOC101743946 isoform X1 [Bombyx mori]
CNG2.99E-44LOC101738542–3.71XP_004931049.1| 0.0e+00| cyclic nucleotide-gated cation channel beta-1 isoform X2 [Bombyx mori]
cAMPCAM0LOC1017428902.94XP_004927631.1| 4.9e-135| neo-calmodulin [Bombyx mori]
signaling pathwayPI3K2.29E-05BGI novel G0001782.57
GPCR3.88E-09Ser-42.32NP_001037502.1| 1.1e-251| 5-hydroxytryptamine receptor [Bombyx mori]
ACO0.000104462LOC1017466362.12XP_021206441.1| 0.0e+00| LOW QUALITY PROTEIN: glutamate receptor ionotropic, delta-2 [Bombyx mori]
ATP1.03E-108LOC1017460121.21XP_004929151.1| 1.3e-193| sodium/potassium-transporting ATPase subunit beta-2 [Bombyx mori]
ROCK6.10E-07BGI novel G000435–1.01
GPCR4.10E-131LOC101736629–1.04XP 022130216.1| 0.0e+00| cytosolic carboxypeptidase 4-like isoform X6 [Pieris rapae]
ROCK0.000248524BGI novel G000569–1.09
GPCR1.33E-12LOC101745640–1.10XP_004934259.1| 5.7e-94| piggyBac transposable element-derived protein 4-like, partial [Bombyx mori]
NMDAE6.28E-16BGI novel G000814–1.23
CAM7.40E-08LOC101741700–1.28XP_004922728.1| 1.3e-163| sarcoplasmic calcium-binding protein [Bombyx mori]
ROCK9.43E-05BGI novel G000712–1.32
GPCR0.000327112LOC105842734–1.35XP_004923418.1| 4.7e-85| piggyBac transposable element-derived protein 4 [Bombyx mori]
ROCK2.33E-07LOC105843016–1.43XP_011555716.1| 1.9e-112| PREDICTED: uncharacterized protein LOC105386785 [Plutella xylostella]
GPCR2.01E-38LOC101739439–1.81XP_004923418.1| 0.0e+00| piggyBac transposable element-derived protein 4 [Bombyx mori]
GPCR1.07E-22LOC101743833–2.94XP_021204873.1| 0.0e+00| cytosolic carboxypeptidase 2 [Bombyx mori]
CNGC2.99E-44LOC101738542–3.71XP_004931049.1| 0.0e+00| cyclic nucleotide-gated cation channel beta-1 isoform X2 [Bombyx mori]

Annotated DEGs of Ca2+/NO/cAMP signaling pathway.

Among these differentially expressed genes (DEGs), several classes of genes and inorganic salt ion transporters, ATP-sensitive inward rectifier potassium channels, uridine nucleosides, potassium voltage-gated channel proteins, calcium/calmodulin-dependent proteins, solute carriers, and organic anions were identified. A large of DEGs were associated with transporters, aquaporins, and water transport genes were heavily annotated, and a class of important signaling molecular genes were annotated. Using qRT-PCR, several types of genes were confirmed in two comparisons (Figure 4A and Table 2), showing elevated transcriptional levels in the hindgut of mutants. The expression patterns of these genes, including inorganic salt ion transporters, calcium/calmodulin-dependent proteins, kynureninase, aquaporins and NOS, were confirmed in mutants and WT comparisons (Figure 4B and Table 2).

TABLE 2

Other gene IDlog2 (ΔBNGR-A2/WT)FDRP-valueNr
aqp1.3191.46E-321.26E-33NP_001106228.1| 2.0e-146| aquaporin [Bombyx mori]
Kin2.4242.23E-153.19E-16NP_001124363.1| 1.1e-195| leucokinin precursor [Bombyx mori]
Kynu6.64000NP_001136399.1| 8.9e-254| kynureninase [Bombyx mori]
LOC101735517–1.5414.32E-117.54E-12XP_012547671.1| 0.0e+00| potassium voltage-gated channel protein eag isoform X3 [Bombyx mori]
LOC1017358713.8103.03E-1727.47E-174XP_012546046.1| 0.0e+00| sodium-independent sulfate anion transporter [Bombyx mori]
LOC101736269–1.0222.58E-421.91E-43XP_004926573.1| 1.9e-274| synaptic vesicle glycoprotein 2C isoform X1 [Bombyx mori]
LOC1017367632.3565.14E-167.17E-17XP_012546620.1| 5.9e-287| putative transporter SVOPL [Bombyx mori]
LOC1017369012.1587.24E-151.06E-15XP_004931037.1| 1.5e-182| carbohydrate sulfotransferase 11 [Bombyx mori]
LOC101736988–2.9635.36E-503.51E-51XP_004929740.2| 0.0e+00| bumetanide-sensitive sodium-(potassium)-chloride cotransporter isoform X1 [Bombyx mori]
LOC1017375792.6703.21E-1271.04E-128XP_012545206.1| 1.1e-287| putative inorganic phosphate cotransporter [Bombyx mori]
LOC1017378701.8532.06E-621.15E-63XP_004929455.1| 2.0e-285| putative transporter SVOPL isoform X1 [Bombyx mori]
LOC1017382229.4883.08E-1221.03E-123XP_004930039.1| 9.8e-305| monocarboxylate transporter 9 isoform X1 [Bombyx mori]
LOC101738542–3.7094.13E-432.99E-44XP_004931049.1| 0.0e+00| cyclic nucleotide-gated cation channel beta-1 isoform X2 [Bombyx mori]
LOC1017386241.8571.50E-063.64E-07XP_004930042.1| 7.0e-151| trypsin, alkaline C [Bombyx mori]
LOC1017387038.7785.04E-622.85E-63XP_004925899.1| 7.9e-308| putative inorganic phosphate cotransporter [Bombyx mori]
LOC1017393826.6722.17E-162.97E-17XP_004934350.2| 0.0e+00| oxysterol-binding protein-related protein 1-like [Bombyx mori]
LOC101739498–3.14900XP_004925267.1| 4.6e-294| synaptic vesicle glycoprotein 2A [Bombyx mori]
LOC101739574–1.0292.17E-381.70E-39XP_021209181.1| 1.5e-239| adenosine 3’-phospho 5’-phosphosulfate transporter 1 isoform X2 [Bombyx mori]
LOC10173972010.6631.41E-2082.96E-210XP_021209245.1| 6.8e-298| putative inorganic phosphate cotransporter [Bombyx mori]
LOC1017397341.8172.74E-057.50E-06XP_021204132.1| 3.4e-245| glutamate receptor 2 isoform X1 [Bombyx mori]
LOC1017398351.2031.02E-241.07E-25XP_012548689.2| 3.7e-310| sodium/hydrogen exchanger 9B2 isoform X2 [Bombyx mori]
LOC1017399984.4502.17E-272.12E-28XP_004924083.1| 3.6e-222| ATP-sensitive inward rectifier potassium channel 12 [Bombyx mori]
LOC1017401531.2418.00E-052.32E-05XP_012548956.1| 0.0e+00| cyclic nucleotide-gated cation channel subunit A [Bombyx mori]
LOC1017403184.5555.66E-177.57E-18XP_012545113.2| 2.1e-181| synaptic vesicle glycoprotein 2B [Bombyx mori]
LOC1017408801.3453.43E-303.14E-31XP_012549174.2| 0.0e+00| potassium voltage-gated channel subfamily KQT member 1 isoform X3 [Bombyx mori]
LOC101741058–1.15000XP_004925076.1| 0.0e+00| ATP-dependent Clp protease ATP-binding subunit clpX-like, mitochondrial [Bombyx mori]
LOC1017411311.3721.11E-072.42E-08XP_004934299.1| 2.8e-82| voltage-dependent L-type calcium channel subunit beta-2-like, partial [Bombyx mori]
LOC101741156–1.3351.90E-957.67E-97XP_004926607.1| 2.2e-293| synaptic vesicle 2-related protein isoform X1 [Bombyx mori]
LOC101741524–1.21400XP_004926609.1| 9.4e-297| synaptic vesicle glycoprotein 2C [Bombyx mori]
LOC101741700–1.2793.27E-077.40E-08XP_004922728.1| 1.3e-163| sarcoplasmic calcium-binding protein [Bombyx mori]
LOC101741853–1.4765.92E-081.28E-08XP_004923628.1| 0.0e+00| potassium channel subfamily K member 9 [Bombyx mori]
LOC101742903–9.21600XP_004931132.1| 2.9e-298| synaptic vesicle glycoprotein 2B isoform X2 [Bombyx mori]
LOC1017434531.0761.85E-054.99E-06XP_012545879.1| 0.0e+00| solute carrier organic anion transporter family member 4A1 isoform X1 [Bombyx mori]
LOC101743804–24.54E-051.28E-05XP_004932640.2| 0.0e+00| glutamate receptor ionotropic, kainate 2 [Bombyx mori]
LOC1017441191.9241.97E-212.28E-22XP_012545589.1| 9.2e-306| putative inorganic phosphate cotransporter [Bombyx mori]
LOC1017441322.5851.93E-103.49E-11NP_001299580.1| 0.0e+00| potassium voltage-gated channel subfamily H member 6-like [Bombyx mori]
LOC1017444274.26000XP_004930921.1| 9.6e-294| monocarboxylate transporter 12 [Bombyx mori]
LOC1017445224.3926.42E-474.41E-48XP_012549165.1| 1.1e-262| monocarboxylate transporter 1 isoform X1 [Bombyx mori]
LOC1017445681.3132.53E-1357.85E-137XP_004929500.1| 2.9e-205| calcium/calmodulin-dependent protein kinase type 1 [Bombyx mori]
LOC101744605–1.1411.53E-251.57E-26XP_004925299.1| 1.4e-287| synaptic vesicle glycoprotein 2B [Bombyx mori]
LOC101744652–1.2644.80E-1571.30E-158XP_012546196.1| 4.5e-299| synaptic vesicle glycoprotein 2A [Bombyx mori]
LOC101744748–2.4291.46E-063.54E-07XP_012545105.1| 1.3e-148| synaptic vesicle glycoprotein 2C-like [Bombyx mori]
LOC1017448571.5461.67E-1165.80E-118XP_004929389.1| 1.9e-296| sodium/hydrogen exchanger 9B2 isoform X1 [Bombyx mori]
LOC1017451377.0995.44E-1751.32E-176XP_004927807.1| 7.2e-259| monocarboxylate transporter 7 [Bombyx mori]
LOC1017454301.4105.33E-138.45E-14XP_012547869.1| 0.0e+00| regulating synaptic membrane exocytosis protein 2 isoform X1 [Bombyx mori]
LOC101745500–1.1830.0003556050.000111834XP_004933911.1| 4.4e-82| heat shock 70 kDa protein, partial [Bombyx mori]
LOC1017456591.8487.47E-111.32E-11XP_012544306.1| 1.1e-246| potassium voltage-gated channel protein Shaw [Bombyx mori]
LOC1017458231.3591.90E-142.83E-15XP_004927809.1| 1.3e-252| proton-coupled amino acid transporter-like protein pathetic [Bombyx mori]
LOC1017458611.5634.86E-147.34E-15XP_004924747.1| 8.7e-277| proton-coupled folate transporter [Bombyx mori]
LOC1017458632.0071.09E-309.81E-32XP_004926413.2| 0.0e+00| sodium-independent sulfate anion transporter isoform X1 [Bombyx mori]
LOC1017460121.2152.77E-1071.03E-108XP_004929151.1| 1.3e-193| sodium/potassium-transporting ATPase subunit beta-2 [Bombyx mori]
LOC1017460433.9371.26E-241.33E-25XP_004925372.1| 4.0e-272| putative inorganic phosphate cotransporter [Bombyx mori]
LOC1017460941.8151.56E-737.68E-75XP_004924179.1| 2.5e-232| 5-hydroxytryptamine receptor 1A isoform X1 [Bombyx mori]
LOC101746144–4.0232.64E-2584.57E-260XP_004927040.1| 0.0e+00| sodium-dependent nutrient amino acid transporter 1-like [Bombyx mori]
LOC1017462261.4131.30E-381.00E-39NP_001296484.1| 0.0e+00| transient receptor potential channel pyrexia [Bombyx mori]
LOC1017466362.1150.0003333380.000104462XP_021206441.1| 0.0e+00| LOW QUALITY PROTEIN: glutamate receptor ionotropic, delta-2 [Bombyx mori]
LOC1017466514.4591.00E-091.90E-10XP_004923028.1| 5.4e-188| probable uridine nucleosidase 2 [Bombyx mori]
LOC1017468873.7169.27E-062.42E-06XP_004923030.1| 1.3e-188| probable uridine nucleosidase 2 isoform X2 [Bombyx mori]
LOC101747040–1.3222.92E-115.04E-12XP_004934379.2| 5.0e-69| proton-coupled amino acid transporter-like protein pathetic [Bombyx mori]
LOC1017471473.3066.38E-266.44E-27XP_004932959.1| 0.0e+00| sodium-coupled monocarboxylate transporter 2 [Bombyx mori]
LOC101747179–3.1032.10E-2294.05E-231XP_012547428.1| 1.6e-283| proton-coupled amino acid transporter-like protein pathetic [Bombyx mori]
LOC1058412431.4761.09E-171.41E-18XP_012548721.1| 3.1e-175| ATP-binding cassette sub-family G member 1-like [Bombyx mori]
LOC105841453–2.19500XP_012544828.1| 1.3e-263| proton-coupled amino acid transporter 1 [Bombyx mori]
LOC1058414946.6870.0002240416.86E-05XP_012545113.2| 6.9e-53| synaptic vesicle glycoprotein 2B [Bombyx mori]
LOC105841866–1.1171.75E-122.85E-13XP_021204603.1| 7.2e-157| glutamate receptor ionotropic, kainate 2-like [Bombyx mori]
Pi-tl1.2411.94E-222.18E-23NP_001036907.1| 7.7e-205| inorganic phosphate transporter 1 [Bombyx mori]
Hsp23.7–3.1151.81E-1674.64E-169NP_001036942.1| 4.0e-115| heat shock protein hsp23.7 precursor [Bombyx mori]
Kmo5.3731.49E-1015.80E-103NP_001106135.1| 6.6e-271| kynurenine 3-monooxygenase [Bombyx mori]
Ser-42.3221.86E-083.88E-09NP_001037502.1| 1.1e-251| 5-hydroxytryptamine receptor [Bombyx mori]
Adh5.3811.95E-093.79E-10NP_001037610.1| 1.9e-147| putative alcohol dehydrogenase [Bombyx mori]

Differentially expressed genes associated with water transport in △BNGR-A2 and WT animals.

In D. melanogaster tubules, as in vertebrate epithelia, convergence of distinct but interacting signal transduction pathways (Ca2+/NO/cGMP) stimulates calcium entry into principal cells via several classes of calcium channels, resulting in increased fluid transport (Choi et al., 1995; Davies, 2000). In this study, we examined whether the NOS pathway could be influenced after BNGR-A2 mutation. In mutant animals, the transcriptional level of NOS2 was significantly up-regulated by 143% (Figure 4C), NOS enzyme activity increased by 15.39% (Figure 4D), and NO concentration increased by 6.85% (Figure 4E), as compared with WT animals. These results supported that the NOS pathway was accelerated in mutant animals.

Insulin and Ecdysone Signaling Pathways in BNGR-A2 Mutants

To further investigate the molecular mechanisms underlying phenotypic defects in the growth of BNGR-A2 animals, transcription levels of key genes in the insulin and ecdysone signaling pathways were investigated in L5D3 larval fat body of WT and BNGR-A2 mutants. Relative mRNA expressions of key genes in insulin signaling pathways including PI3K, AKT, S6K, and 4EBP were down-regulated, suggesting that knockout of BNGR-A2 affected the insulin signaling pathway (Figure 5A). Lower ecdysteroid titers and downregulation of 20E-related key genes, including ECRA, E75A, E75B and HR3, suggested that the 20E signaling pathway was affected after BNGR-A2 disruption (Figures 5B,C).

FIGURE 5

Discussion

ITP belongs to the class of CHH family neuropeptides and has been functionally characterized in many insect species (Drexler et al., 2007; Dircksen et al., 2008; Webster et al., 2012; Nagai et al., 2014; Yu et al., 2016). However, little is known regarding ITP signaling pathways in insects, and no insect ITPRs have been identified thus far, except in B. mori (Nagai et al., 2014). Here, we comprehensively investigated physiological functions of BNGR-A2 in B. mori by using CRISPR/Cas9-mediated gene disruption. The NOS enzyme activity, NO content, and calmodulin transcription levels were significantly higher when BNGR-A2 was knocked out. The downstream Ca2+/NO/cGMP signaling pathways were activated, resulting in the loss of water in the body. The decrease in the water content and concurrent increase in the amount of excretion in the mutant animals resulted in insufficient pressure in the adult body after emergence and ultimately led to the failure of wing expansion (Figure 6). BNGR-A2 deficiency affected silkworm growth and development by inhibiting ecdysone and insulin signaling. These results reveal the important role BNGR-A2 plays in insect larval growth rate, wing expansion, and water homeostasis (Figure 6).

FIGURE 6

In the insect excretory system, the Malpighian tubules and hindgut are involved in water preservation and absorption, which aid in maintaining water balance, allowing the organism to tolerate external hypertonicity or dryness (Meredith et al., 1996). This process is strictly regulated by neuropeptides, including diuretic hormone 31 (DH31) (Coast et al., 2001), diuretic hormone 44 (DH44), leucokinin (Cannell et al., 2016), and capa (Terhzaz et al., 2015), as well as antidiuretic hormones such as ITP. In insects, the diuretic peptide acts as a scavenging hormone, promoting the formation of primary urine in malpighian tubules and inducing rapid circulation of hemolymph and frequent hypoperistalsis of the digestive tract. Diuretic peptides promote the entry of raw urine into the hindgut, aid in moving dry food residue through the intestine (Gäde, 2004). The antidiuretic hormone acts on the hindgut to promote water reabsorption and plays an important role in insect water retention (Gäde, 2004). In previous studies, the NO signaling pathway was reported to be involved in the preservation of water in Diptera (Davies, 2000). NO is produced from the amino acid L-arginine by the enzyme action of nitric oxide synthase (NOS). NOS is encoded by a large multigene family, and occurs several isoforms. NOS isoforms have been classified into either calcium-dependent or calcium-independent enzymes, and belong to constitutive (ecNOS and ncNOS) or inducible NOS (iNOS) (Davies, 2000). Liquid transport was found to be controlled by NOS in Drosophila Malpighian tubules (Davies, 2000). In the late stage of Plasmodium berghei infection of A. stephensi in Diptera, NOS activity was significantly increased, likely due to the increase of hemolymph clearance (Dimopoulos et al., 1998). The NOS isoform in mosquitoes also contains Ca2+/calmodulin binding sites, indicating that NOS activation is associated with a regulatory calcium signaling process (Davies, 2000). In B. mori, two types of NO synthase were found in the malpighian tube. During the gut purging phase of B. mori larvae, here considered the excretion phase, the activity of NOS in malpighian tubule was significantly increased (Choi et al., 1995). We also found NOS activity and NO content to be significantly higher in the mutant B. mori than in WT animals. Based on hindgut transcriptome and qPCR data, downstream Ca2+/NO/cGMP signaling pathways were activated and 46 DEGs were annotated in the pathways named above when BNGR-A2 was knocked out (Table 2). The transcription levels of the DEGs in mutant B. mori, including those encoding aquaporins, ion cotransporters, calmodulin (CaM), and NOS, were significantly higher than those of the WT animals. This is the first work to show that knockout of BNGR-A2 in silkworms caused a (1.42 ± 0.035)% reduction in body water. Previous research showed that in humans, 2% dehydration leads to diminished mental activity and cognitive ability (Grandjean and Grandjean, 2007). Our results revealed that the Ca2+/NO/cGMP signaling pathway was involved in the transportation of water in B. mori, indicating that BNGR-A2 affects the resorption of the hindgut water.

BNGR-A2 was also found to play a critical role in the development of B. mori. The mutant larvae had significantly lower body weights than WT larvae at the same stage. The mutants showed slower preadult development, and larval development was extended by 3 days. Although development was prolonged, the pupal body weight showed no significant difference between mutant and WT B. mori. The size of adults was primarily determined by two factors: duration of growth and growth rate. The ecdysone and juvenile hormones synergistically regulate insect developmental metamorphosis and determine the duration of insect growth. The growth and death signals, such as insulin, nutrition, and cell contact inhibition, as well as their conduction pathways, control cell division, growth, differentiation, death, and ultimately determine the insect’s growth rate (Colombani et al., 2005; OrmeSally and Leevers, 2005). We discovered that the ecdysone titer is lower and the ecdysone pathway genes were down-regulated in the Δ BNGR-A2. The key genes of the insulin signaling pathway in the mutant were down-regulated, and the results indicated that the mutation of BNGR-A2 disrupted the expression of other genes in the insulin signaling pathway, affecting the nutritional signaling pathway of the silkworm. This ultimately led to changes in the duration of development of mutant silkworms. However, the regulatory interaction of BNGR-A2 between water homeostasis and ecdysone, as well as the insulin signaling pathways, needs further investigation.

Statements

Data availability statement

All datasets generated for this study are included in the article/Supplementary Material.

Author contributions

LS and AT designed the research and wrote the manuscript. ZZ, RZ, YY, and FY performed the experiments and analyzed the data. AT revised the manuscript. All authors listed have approved the manuscript for publication.

Funding

This work was supported by grants from the National Science Foundation of China (31925007 and U1738110).

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.2020.00424/full#supplementary-material

FIGURE S1

Transcriptomic profiles influenced by BNGR-A2-loss in silkworm hindguts by RNA-seq. (A) M-versus-A plot, MA plot diagram of DEGs in the transcriptome. (B) Gene Ontology (GO) term of DEGs in the transcriptome. (C) Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations of differentially expressed genes (DEGs) in the transcriptome.

TABLE S1

List of primers used in this study.

References

  • 1

    AudsleyN.McintoshC.PhillipsJ. E. (1992). Actions of Ion-transport peptide from locust corpus cardiacum on several hindgut transport processes.J. Exp. Biol.173275288.

  • 2

    AudsleyN.MeredithJ.PhillipsJ. E. (2006). Haemolymph levels of Schistocerca gregaria ion transport peptide and ion transport-like peptide.Physiol. Entomol.31154163. 10.1111/j.1365-3032.2006.00500.x

  • 3

    BegumK.LiB.BeemanR. W.ParkY. (2009). Functions of ion transport peptide and ion transport peptide-like in the red flour beetle Tribolium castaneum.Insect Biochem. Mol.Biol.39717725. 10.1016/j.ibmb.2009.08.005

  • 4

    CannellE.DornanA. J.HalbergK. A.TerhzazS.DowJ. A. T.DaviesS. A. (2016). The corticotropin-releasing factor- like diuretic hormone 44 (DH44) and kinin neuropeptides modulate desiccation and starvation tolerance in Drosophila melanogaster.Peptides8096107. 10.1016/j.peptides.2016.02.004

  • 5

    ChoiS. K.ChoiH. K.Kadono-OkudaK.TaniaiK.KatoY.YamamotoM.et al (1995). Occurrence of novel types of nitric oxide synthase in the silkworm, Bombyx mori.Biochem. Biophys. Res. Commun.207452459. 10.1006/bbrc.1995.1209

  • 6

    CoastG. M.WebsterS. G.ScheggK. M.TobeS. S.SchooleyD. A. (2001). The Drosophila melanogaster homologue of an insect calcitonin-like diuretic peptide stimulates V-ATPase activity in fruit fly Malpighian tubules.J. Exp. Biol.20417951804.

  • 7

    ColombaniJ.BianchiniL.LayalleS.PondevilleE.Dauphin-VillemantC.AntoniewskiC.et al (2005). Antagonistic actions of ecdysone and insulins determine final size in Drosophila.Science310667670. 10.1126/science.1119432

  • 8

    DaviesS. A. (2000). Nitric oxide signalling in insects.Insect Biochem. Mol. Biol.3011231138. 10.1016/s0965-1748(00)00118-1

  • 9

    DimopoulosG.SeeleyD.WolfA.KafatosF. C. (1998). Malaria infection of the mosquito Anopheles gambiae activates immune-responsive genes during critical transition stages of the parasite life cycle.EMBO J.1761156123. 10.1093/emboj/17.21.6115

  • 10

    DircksenH.TesfaiL. K.AlbusC.NässelD. R. (2008). Ion transport peptide splice forms in central and peripheral neurons throughout postembryogenesis of Drosophila melanogaster.J. Comp. Neurol.5092341. 10.1002/cne.21715

  • 11

    DrexlerA. L.HarrisC. C.delaPenaM. G.Asuncion-UchiM.ChungS.WebsterS.et al (2007). Molecular characterization and cell-specific expression of an ion transport peptide in the tobacco hornworm, Manduca sexta.Cell Tissue Res.329391408. 10.1007/s00441-007-0391-9

  • 12

    GädeG. (2004). Regulation of intermediary metabolism and water balance of insects by neuropeptides.Annu. Rev. Entomol.4993113.

  • 13

    GálikováM.DircksenH.NässelD. R. (2018). The thirsty fly: ion transport peptide (ITP) is a novel endocrine regulator of water homeostasis in Drosophila.PLoS Genet.14:e1007618. 10.1371/journal.pgen.1007618

  • 14

    GrandjeanA. C.GrandjeanN. R. (2007). Dehydration and cognitive performance.J. Am. Coll. Nutr.26549S544S.

  • 15

    Hermann-LuiblC.YoshiiT.SenthilanP. R.DircksenH.Helfrich-FörsterC. (2014). The ion transport peptide is a new functional clock neuropeptide in the fruit fly Drosophila melanogaster.J. Neurosci.3495229536. 10.1523/JNEUROSCI.0111-14.2014

  • 16

    JohardH. A. D.YoishiiT.DircksenH.CusumanoP.RouyerF.Helfrich-FörsterC.et al (2009). Peptidergic clock neurons in Drosophila: ion transport peptide and short neuropeptide F in subsets of dorsal and ventral lateral neurons.J. Comp. Neurol.5165973. 10.1002/cne.22099

  • 17

    LiZ. Q.YouL.ZengB. S.LingL.XuJ.ChenX.et al (2015). Ectopic expression of ecdysone oxidase impairs tissue degeneration in Bombyx mori.Proc. Biol. Sci.282:20150513. 10.1098/rspb.2015.0513

  • 18

    MeredithJ.RingM.MacinsA.MarschallJ.ChengN. N.TheilmannD.et al (1996). Locust ion transport peptide (ITP): primary structure, cDNA and expression in a baculovirus system.J. Exp. Biol.19910531061.

  • 19

    MirthC.TrumanJ. W.RiddifordL. M. (2005). The role of the prothoracic gland in determining critical weight for metamorphosis in Drosophila melanogaster.Curr. Biol.1517961807. 10.1016/j.cub.2005.09.017

  • 20

    NagaiC.Mabashi-AsazumaH.NagasawaH.NagataS. (2014). Identification and characterization of receptors for ion transport peptide (ITP) and ITP-like (ITPL) in the silkworm Bombyx mori.J. Biol. Chem.2893216632177. 10.1074/jbc.M114.590646

  • 21

    NässelD. R.ZandawalaM. (2019). Recent advances in neuropeptide signaling in Drosophila, from genes to physiology and behavior.Prog. Neurobiol.179:101607. 10.1016/j.pneurobio.2019.02.003

  • 22

    OrmeSallyM. H.LeeversJ. (2005). Flies on steroids: the interplay between ecdysone and insulin signaling.Cell Metab.2277278. 10.1016/j.cmet.2005.10.005

  • 23

    ParkY.KimH.LiD.AdamsM. (2004). “A Novel function of ion transport peptide in Drosophila ecdysis,” in Proceedings of the Program and Abstracts 45th Annual Drosophila Research Conference, Washington, DC.

  • 24

    PfafflM. W.HorganG. W.DempfleL. (2002). Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR.Nucleic Acids Res.30:e36.

  • 25

    SonobeH.NishimuraT.SonobeM.NakatsujiT.YanagiharaR.KawakamiT.et al (2001). The molt-inhibiting hormone in the American crayfish Procambarus clarkii. its chemical and biological activity.Gen. Comp. Endocrinol.121196204. 10.1006/gcen.2000.7586

  • 26

    TanA.FuG.JinL.GuoQ.LiZ.NiuB.et al (2013). Transgene-based, female-specific lethality system for genetic sexing of the silkworm, Bombyx mori.Proc. Natl. Acad. Sci. U.S.A.11067666770. 10.1073/pnas.1221700110

  • 27

    TerhzazS.TeetsN. M.CabreroP.HendersonL.RitchieM. G.NachmanR. J.et al (2015). Insect capa neuropeptides impact desiccation and cold tolerance.Proc. Natl. Acad. Sci. U.S.A.11228822887. 10.1073/pnas.1501518112

  • 28

    WebsterS. G.KellerR.DircksenH. (2012). The CHH-superfamily of multifunctional peptide hormones controlling crustacean metabolism, osmoregulation, moulting, and reproduction.Gen. Comp. Endocr.175217233. 10.1016/j.ygcen.2011.11.035

  • 29

    XuJ.ChenS. Q.ZengB. S.JamesA. A.TanA. J.HuangY. P. (2017a). Bombyx mori P-element somatic inhibitor (BmPSI) is a key auxiliary factor for silkworm male sex determination.PLoS Genet.13:e1006576. 10.1371/journal.pgen.1006576

  • 30

    XuJ.ZhanS.ChenS. Q.ZengB. S.LiZ. Q.JamesA. A.et al (2017b). Sexually dimorphic traits in the silkworm, Bombyx mori, are regulated by doublesex.Insect Biochem. Mol. Biol.804251. 10.1016/j.ibmb.2016.11.005

  • 31

    YuB.LiD. T.WangS. L.XuH. J.BaoY. Y.ZhangC. X. (2016). Ion transport peptide (ITP) regulates wing expansion and cuticle melanism in the brown planthopper, Nilaparvata lugens.Insect Mol. Biol.25778787. 10.1111/imb.12262

Summary

Keywords

Bombyx mori, ion transport peptide receptor, CRISPR/Cas9, water homeostasis, development

Citation

Sun L, Zhang Z, Zhang R, Yu Y, Yang F and Tan A (2020) Molecular Disruption of Ion Transport Peptide Receptor Results in Impaired Water Homeostasis and Developmental Defects in Bombyx mori. Front. Physiol. 11:424. doi: 10.3389/fphys.2020.00424

Received

13 December 2019

Accepted

07 April 2020

Published

20 May 2020

Volume

11 - 2020

Edited by

Bin Tang, Hangzhou Normal University, China

Reviewed by

Shinji Nagata, The University of Tokyo, Japan; Keping Chen, Jiangsu University, China

Updates

Copyright

*Correspondence: Anjiang Tan,

This article was submitted to Invertebrate Physiology, a section of the journal Frontiers in Physiology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics