Front. Physiol. Frontiers in Physiology Front. Physiol. 1664-042X Frontiers Media S.A. 10.3389/fphys.2018.00057 Physiology Original Research Identification and Functional Analysis of the First Aquaporin from Striped Stem Borer, Chilo suppressalis Lu Ming-Xing 1 2 Pan Dan-Dan 1 2 Xu Jing 1 2 Liu Yang 3 Wang Gui-Rong 3 Du Yu-Zhou 1 2 * 1College of Horticulture and Plant Protection and Institute of Applied Entomology, Yangzhou University, Yangzhou, China 2Jiangsu Key Laboratory of Crop Genetics and Physiology, Co-Innovation Center for Modern Production Technology of Grain Crops, Yangzhou University, Yangzhou, China 3State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing, China

Edited by: Fernando Ariel Genta, Fundação Oswaldo Cruz (Fiocruz), Brazil

Reviewed by: Wen-Jia Yang, Guiyang University, China; Yifan Zhai, Shandong Academy of Agricultural Sciences, China

*Correspondence: Yu-Zhou Du yzdu@yzu.edu.cn

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

 05 02 2018 2018 9 57 03 10 2017 17 01 2018 Copyright © 2018 Lu, Pan, Xu, Liu, Wang and Du. 2018 Lu, Pan, Xu, Liu, Wang and Du

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Aquaporins are integral membrane proteins some of which form high capacity water-selective channels, promoting water permeation across cell membranes. In this study, we isolated the aquaporin transcript (CsDrip1) of Chilo suppressalis, one of the important rice pests. CsDrip1 included two variants, CsDrip1_v1 and CsDrip1_v2. Although CsDrip1_v2 sequence (>409 bp) was longer than CsDrip1_v1, they possessed the same open reading frame (ORF). Protein structure and topology of CsDrip1 was analyzed using a predicted model, and the results demonstrated the conserved properties of insect water-specific aquaporins, including 6 transmembrane domains, 2 NPA motifs, ar/R constriction region (Phe69, His194, Ser203, and Arg209) and the C-terminal peptide sequence ending in “SYDF.” Our data revealed that the Xenopus oocytes expressing CsDrip1 indicated CsDrip1 could transport water instead of glycerol, trehalose and urea. Further, the transcript of CsDrip1 expressed ubiquitously but differentially in different tissues or organs and developmental stages of C. suppressalis. CsDrip1 mRNA exhibited the highest level of expression within hindgut and the third instar larvae. Regardless of pupae and adults, there were significantly different expression levels of CsDrip1 gene between male and female. Different from at low temperature, the transcript of CsDrip1 in larvae exposed to high temperature was increased significantly. Moreover, the mRNA levels of CsDrip1 in the third instar larvae, the fifth instar larvae, pupae (male and female), and adults (male and female) under different humidities were investigated. However, the mRNA levels of CsDrip1 of only female and male adults were changed remarkably. In conclusions, CsDrip1 plays important roles in maintaining water homeostasis in this important rice pest.

 Drosophila integral protein (Drip) Chilo suppressalis structure expression functional assay 31371937 31401733 National Natural Science Foundation of China10.13039/501100001809 Introduction

Obviously, water is one of the most fundamental molecules for all living organisms. Scientists have confirmed that in addition to simple diffusion, there were high capacity water-selective channels to account for the high water permeability in certain tissues of the animals (Preston et al., 1992; Shakesby et al., 2009). Aquaporins (AQPs), often known as water channels, are integral membrane proteins that regulate the flow of water driven by osmotic gradients (Campbell et al., 2008; Benga, 2009). AQPs exist as tetramers in the cell membrane with each monomer functioning as a water channel (Agre et al., 1987; van Hoek et al., 1991). These AQPs are present in both prokaryotes and eukaryotes species and play an important role in the water transport system (Heymann and Engel, 1999; Agre, 2006; Campbell et al., 2008). Some AQPs can also transport a number of small solutes, typically glycerol or urea (Rojek et al., 2008). Current research on AQPs primarily focuses on plants and vertebrates (Verkman, 2005; Campbell et al., 2008; Maurel et al., 2008). To date, 13 AQP types have been isolated from mammals, named AQP0-AQP12 (Yasui, 2004; Ishibashi, 2006), However, only a few AQP genes were identified and characterized from insects, among which just four species come from Lepidoptera, Bombyx mori, Grapholita molesta, Spodoptera litura, and Ectropis obliqua (Dow et al., 1995; LeCahérec et al., 1996a,b; Echevarria et al., 2001; Yanochko and Yool, 2002; Duchesne et al., 2003; Kaufmann et al., 2005; Kikawada et al., 2008; Kambara et al., 2009; Kataoka et al., 2009a,b; Shakesby et al., 2009; Goto et al., 2011; Herraiz et al., 2011; Mathew et al., 2011; Philip et al., 2011; Azuma et al., 2012; Liu et al., 2013; Fabrick et al., 2014; Ibanez et al., 2014; Li et al., 2016; Van Ekert et al., 2016; Liu and Piermarini, 2017). The phylogenetic analyses of insect AQPs had revealed the existence of five major subfamilies, including the Drip, Prip, Bib, Eglps, and AQP12L, and members of the Drip and Prip subfamilies typically were water selective channels (Chawn and Nicolson, 2004; Kambara et al., 2009; Herraiz et al., 2011; Mathew et al., 2011; Drake et al., 2015; Finn and Cerda, 2015; Van Ekert et al., 2016). Meantime, according to the database of genome and transcriptome, the only one AQP in Chilo suppressalis, two AQPs in Bombyx mori, three AQPs in Danaus plexippus and Manduca sexta, and five AQPs in Manduca sexta was identified (http://www.insect-genome.com/data/detail.php?id=7) (Yin et al., 2014), and basing on the homology of AQPs, the above species all could possess the water-selective AQPs. Therefore, Drip and Prip are very important to maintain the balance of water in insect.

Otherwise, some studies suggested that AQPs played an important role in the physiological functions of insect (Goto et al., 2011; Benoit et al., 2014b; Fabrick et al., 2014; Drake et al., 2015). For example, freeze tolerance of insects was related to the ability to remove water from cells by AQPs (Philip et al., 2011). The removal of water from the cells could suspend metabolic processes or avoid damaging ice crystal formation (Spring et al., 2009). Freeze tolerance also needed to accumulate glycerol in the cells, a role admirably suited to the aquaglyceroporins (Spring et al., 2009). Down regulation of AQPs in Aedes aegypti enhanced mosquito desiccation resistance (Drake et al., 2015). The female tsetse flies of Glossina morsitans morsitans have been studied to elucidate the role of AQPs in heat tolerance (Benoit et al., 2014a,b).

The striped stem borer, Chilo suppressalis (Walker) (Insecta: Lepidoptera: Pyralidae), an important rice pest widely distributed in Asia, has caused significant damage to rice crops in China, especially to hybrid rice varieties in recent years. In the district of Yangzhou (32.23°N, 119.26°E), Jiangsu province, China, C. suppressalis has two complete and a partial third generation each year, and the larvae enter facultative diapause in winter (Lu et al., 2013). According to our studies, in March of 2010, field-collected larvae could survive at −21oC (Lu et al., 2012). Overwintering larvae of C. suppressalis could acquire freeze tolerance (Tsumuki and Hirai, 2007). Some studies demonstrated AQPs might play very important role in the cold tolerance of C. suppressalis (Izumi et al., 2006, 2007). However, they did not further study any aquaporin gene. It is well-known that C. suppressalis need very high humidity condition in the life cycle (Shang et al., 1979; Lu et al., 2014), and AQPs may help to maintain water homeostasis.

Thus, in this paper we firstly described the characteristics of C. suppressalis Drip1 (CsDrip1), and assayed the abundance of CsDrip1 in different tissues or organs and developmental stages. Secondly, in order to understand the CsDrip1 regulation under different humidities, the CsDrip1 mRNA levels of different developmental stages and sexes of C. suppressalis were investigated under different humidities. Moreover, in order to understand the relationship between CsDrip1 regulation and temperature, we also studied the expression patterns of CsDrip1 mRNA under various temperatures. Last but not the least, to further demonstrate the functions of CsDrip1, functional oocyte swelling assays were executed by water and three kinds of solutes. These studies will help us understand the role of AQPs in C. suppressalis, and also may provide insights in developing strategies for the control of this pest.

 Materials and methods Insects

The population of C. suppressalis was collected from the suburb of Yangzhou (32.39°N, 119.42°E). The rice stem borers were reared in an environmental chamber at 28 ± 1°C, 16:8 (L: D) and RH = 70 ± 5% (Shang et al., 1979).

 Cloning and RACEs

Total RNA of the fifth instar larva was extracted by the SV Total RNA Isolation System (Promega Z3100) combined with DNase digestion to eliminate DNA contamination. Total cDNA was synthesized by oligo(dT)18 primer (TaKaRa). Degenerate primers DP-F and DP- R (Table 1) were used to amplify the partial segments of AQP. The full-length cDNA of the CsDrip1 gene was determined using 5′- and 3′-RACE (SMART RACE, Clontech). The primers used (VA and VB) were shown in Table 1. The full length sequence of CsDrip1 was confirmed by the template of RACE 5′ cDNA.

Primers used in this study.

Gene Primer name Primer sequences (5′ → 3′) Tm (°C) Length (bp)
RT-PCR
CsDrip1 F CACATCAAYCCMGCBGTCAC 431
R GGNCCCRCCCARTAMACCCA
RACE PCR
CsDrip1 3′ CCACGAAGACAGCACCAACCGCAAGA 68.0 666
variant A 5′ CCTTCATGCGGCTAACGTACTCCT 68.0 812
CsDrip1 3′ CCTCCATCACCCGTGCTCGTTTT 68.0 757
variant B 5′ TGGCAGCAGCAGTGGGCTCGTTG 68.0 1093
FUNCTIONAL ASSAYS
CsDrip1 F TCAACTAGTGCCACCATGAAAACGGATTACGCTGT 774
R TCAGCGGCCGCTTAGAAGTCATAGGAGCCGC
QPCR
CsDrip1 F GTGGATGAGATGCAGAAAAGGA 59.4 120
R AAGCGATGTCAGCACAAAGGT
H3 F CTGCACCAAGCACTGGTGGA 56.0 184
R TAGCGGCGGACTGGAAACG
EF1 F AAAATGGACTCGACTGAACCCC 56.6 137
R TCTCCGTGCCAACCAGAAATA
18S F GTGATGGGACGAGTGCTTTTATT 62.5 258
R GCTGCCTTCCTTGGATGTGG
ACTIN F AAAGAAACAGCAAAAGTCGGGG 56.0 243
R GTTCAATGGAGGTTCGGTAAGTAAA
TUB F GAGGGCATGGACGAGATGGA 60.4 178
R ACGACGGTACGAGTATGACGGG
UBI F TCACCGACAGCAAACCAGACT 60.2 219
R GGAAGAAAACACCCCCCTCATATA

The qPCR primers used in this study were validated.

 Sample preparation

The rice stem borers were reared successively to the third generation in the seedlings. Then, the egg masses, the first, second, third, fourth, fifth instar larvae, pupae (male and female), and 1-day adults (male and female) were randomly selected for the experiment. Each experiment was repeated four times. The fifth instar deep diapausing larvae of C. suppressalis collected from rice stubble were anesthetized on ice before dissection. Heads (HE), epidermis (EP), fat body (FB), foregut (FG), midgut (MG), hindgut (HG), Malpighian tubules (MT) and haemocytes (HC) were collected and rinsed with 0.9% NaCl. The samples were frozen immediately in liquid nitrogen and stored at −70°C until the experiment.

 Temperature stress treatment

The larvae used in experiments were all 5th instars of similar body size and were assigned randomly to each experimental group. Each group of larvae was confined individually in glass tubes (relative humidity, 90 ± 10%) and exposed to a given temperature (including −11, −9, −3, 0, 27, 36, 39, and 42°C) for 2 h in a constant temperature incubator (DC-3010, Jiangnan equipment). The larvae were recovered at 27 ± 1°C for 2 h, and the survived larvae were frozen in liquid nitrogen and then stored at −70°C.

 Humidity treatment

The third instar larvae, fifth instar larvae, pupae (male and female) and adults (male and female) were treated in humidity chamber (HTC-100, SANTN, Shanghai, China) at 27°C, respectively under four different relative humidity levels (RH) (25, 50, 75, and 95%) for 24 h, and additionally the third instar larvae were treated for 12 h. Each treatment was replicated 30 times, and each replicate consisted of one insect. The humidities chosen were based on a prior pilot experiment. At the same time, survival in each treatment was assessed, and the survived larvae were treated as above.

 Functional oocyte swelling assays

The vector construction followed previously reported protocol (Chang et al., 2015). Firstly, the entire coding region of CsDrip1 was amplified with a high-fidelity polymerase (PrimeSTAR®, HS DNA polymerase. TaKaRa, Tokyo, Japan) using primers with Kozak sequence and restriction enzyme cutting sites (SpeI and NotI). And the PCR products were digested with SpeI and NotI, and subcloned into pT7Ts vector (Invitrogen, Carlsbad, California, USA), and then plasmids were fully linearized with SmaI. The cRNAs of CsDrip1 were synthesized in vitro using mMESSAGE mMACHINE T7 kit (Ambion, Austin, TX, USA). Purified cRNAs were resuspended in nuclease-free water at a concentration of 0.2 μg/μl and stored at −80°C.

Unfertilized stage V and VI oocytes of Xenopus were defolliculated with 2 mg/ml collagenase I (GIBCO, Carlsbad, CA) in washing buffer (96 mM NaCl, 2 mM KCl, 5 mM MgCl2, and 5 mM HEPES [pH 7.6]) for about 1 h at room temperature (26°C). After being cultured overnight at 18°C, oocytes were microinjected with 27.6 nl CsDrip1 cRNA (5.52 ng) and 27.6 nl nuclease-free water as control. After injection, oocytes were incubated for 3 days at 18°C in 1X Ringer's solution (96 mM NaCl, 2 mM KCl, 5 mM MgCl2, 0.8 mM CaCl2, and 5 mM HEPES [pH = 7.6]) supplemented with 5% dialysed horse serum, 50 mg/ml tetracycline, 100 mg/ml streptomycin and 550 mg/ml sodium pyruvate.

Osmotic water permeability (Pf) was measured as previously described (Kataoka et al., 2009a,b). Oocytes was transferred to a 3-fold dilution of 1X Ringer's solution with distilled water and images were acquired of the oocyte silhouette every 15 s through a CCD camera DP-72 (Olympus, Tokyo, Japan) attached to a Olympus SZX16 stereomicroscope up to 5 min. The osmotic water permeability (Pf) was calculated as in previous reports (Zhang et al., 1990; Preston et al., 1992; Kataoka et al., 2009a,b) by the following equation: Pf = V0×d(V/V0)/dt/[S × Vw×(Osmin-Osmout)] where V0 is the oocyte initial volume (V0 = 9 × 10−4 cm3), S is the oocyte surface area (S = 0.045 cm2), Vw is the molecular volume of the water (Vw = 18 cm3/mol) and Osmin is 202 mmol·kg−1 and Osmout is 59 mmol·kg−1. Relative oocyte volume (V/V0) was calculated from the relative area (A/A0) in the focal plane, V/V0 = (A/A0)3/2.

Oocytes were transferred in an isotonic solution which containing 140 mM of solutes (glycerol, trehalose or urea) for solutes transport assays. To maintain the osmotic equilibrium, the increase in oocyte volume corresponds to the water influx accompanying the solute uptake. The volume changes were recorded for 5 min in the same way as described above. Apparent solute permeability was calculated from the equation: Psol = [d(V/V0)/dt] × (V0/S) (LeCahérec et al., 1996a,b; Duchesne et al., 2003). Water and solute permeabilities were performed at least for nine different Xenopus oocytes.

 Quantitative real-time PCR (qPCR) analysis

Total RNA was extracted by the SV Total RNA Isolation System (Promega Z3100), followed by DNase treatment to eliminate DNA contamination. The integrity of the RNA in all samples was verified by comparing the ribosomal RNA bands in ethidium bromide-stained gels. RNA sample purity was estimated using spectrophotometric measurements at 260 and 280 nm (Eppendorf BioPhotometer plus). Total RNA (500 ng) was reverse-transcribed into first-strand cDNA using the SuperScript II reverse transcription (RT)-PCR kit (Bio-Rad). The volume of reaction mixture was 20 μl. Each reaction mixture contained 10 μl of iTaq Universal SYBR Green supermix (2x) (Bio-Rad), 1 μl of each of gene specific primers (Table 1), 2 μl of cDNA templates, and 6 μl ddH2O. Reactions were carried out on a CFX-96 real-time PCR system (Bio-Rad). The amplification efficiencies of the target and reference genes were similar in this study. Therefore, the quantity of CsDrip1 mRNA was calculated using the 2−ΔΔCt method (Nolan et al., 2006; Schmittgen and Livak, 2008; Bustin et al., 2009). Relative expression levels of CsDrip1 in different tissue or organs were normalized with histone 3 (H3), elongation factor 1 (EF1) for different developmental stages and temperature stress (Xu et al., 2017). And previous tests of stability of the reference gene demonstrated 18S rRNA for the third larvae, ACTIN for the fifth larvae, TUB for the male pupae, UBI for the female pupae and male adults, and EF1 for the female adults under different humidity was most suitable respectively, and corresponding reference genes were selected to normalize (data not shown). Following qPCR, the homogeneity of the PCR products was confirmed by melting curve analysis, which was read 5 s per 0.5°C, increment from 65°C to 95°C. Every treatment included four biological replicates, and every repeat was run in triplicate.

 Bioinformatic analysis

The open reading frames (ORFs) were identified with the aid of the ORF Finder software (http://www.ncbi.nlm.nih.gov/gorf/gorf.html). The deduced amino acid sequences were aligned using ClustalX software. Sequence analysis tools of the ExPASy Molecular Biology Server of Swiss Institute of Bioinformatics, including Translate, Compute pI/MW, Blast and TMHMM, were used to analyze the deduced CsDrip1 protein sequence. Phosphorylation and kinases sites were predicted using NetPhos 2.0 and NetPhosK 1.0, respectively (Blom et al., 1999, 2004). Amino acid sequences were used to estimate phylogeny with the neighbor-joining, minimum evolution, maximum likelihood and maximum parsimony methods. Phylogenetic trees were constructed with 1000 bootstrap replicates using MEGA version 7.0 (Kumar et al., 2016).

 Computational molecular modeling

Homology models were generated using Protein Homology/analogy recognition engine V 2.0 (http://www.sbg.bio.ic.ac.uk/~phyre2/html) (Kelley and Sternberg, 2009). Briefly, the CsDrip1 amino acid sequence was aligned by the Phyre2, and the best model of bovine AQP1 X-ray A* structure (PDB ID: 1J4N) was used for modeling analyses (Sui et al., 2001). The Chimera Tool was used to visualize the three-dimensional coordinates for the atoms of the model (Pettersen et al., 2004).

 Statistical analysis

The data were tested for normality using the Shapiro-Wilk's test. Homogeneity of variances among different groups was evaluated by the Levene test. All the data was log-transformed when necessary. Then, the significance of differences between treatments was identified with either a Tukey's test (Homogeneity of variances) or Dunnett's T3 test (Non-homogeneous) for multiple comparisons. In the experiments with two groups, significant differences were determined by Student's t-test. The data were analyzed using SPSS16.0 software, and denoted as means ± SE (standard error) (Pallant, 2005).

 Results Isolation, cloning, sequencing, and structure of <italic>Cs</italic>Drip1

Degenerate primers based on conserved regions from several insect aquaporins were used to amplify a 431 bp partial fragment from C. suppressalis cDNA. The cloned fragment was sequenced, and a BLAST analysis of its deduced amino acid fragment revealed apparent sequence homology with insect Drip1 (data not shown). The 1,416 and 1,825 bp full-length of variant A and variant B, including the UTRs was obtained, respectively through 5′and 3′ RACE (GenBank accession no. JQ011314 and JQ011315). Two Drip1 types of C. suppressalis (CsDrip1_v1 and CsDrip1_v2) were obtained in the present study. Variant A and variant B possessed the same open reading frame (ORF) of 258 amino acids with a predicted molecular weight of 26.9 kDa and an isoelectric point of 6.5. The protein sequence that we referred to as CsDrip1 possessed the hallmarks of the aquaporin family, including the classical “NPA” boxes (residue 89–91; residue 206–208) and 6 putative transmembrane regions (Figures 1A,B). According to the hydropathy analysis, CsDrip1 contained cytoplasmic N- and C-terminal domains and the C-terminus ended with “SYDF”. In the amino acid sequence, seven potential phosphorylation sites (Tyr5, Ser15, Ser16, Ser176, Ser210, Ser250, and Ser253) and three potential protein kinase C-specific sites (Thr3, Ser21, Ser255) were also identified (Figures 1A,B). To investigate the potential structure-function relationship of CsDrip1, we generated its homology model with Phyre using bovine aquaporin-1 (1J4N) as a template (Sui et al., 2001). The structure of CsDrip1 model was very similar to that of aquaporin-1 (confidence, 100% statistically and identity, 36% based on sequence alignment). CsDrip1 models showed two tandem structural repeats, each consisting of three transmembrane helices (TM1-3 and TM4-6) and a short α-helix in loops B and E, each containing an NPA motif predicted to line one side of the pore (Figure 1C), which was called the “aquaporin fold” (Murata et al., 2000). Conserved NPAs formed the canonical structure in the center of the pore (Fu et al., 2000), and the structural model provided the defining force that orients water as it passed through the midpoint of the channel (Hoa et al., 2009). Residues that comprised the Ar/R constriction (Phe69, His194, Ser203, and Arg209) were found in CsDrip1 and predicted to establish water selectivity (Figures 1C,D).

Sequence analysis, predicted topology and structure of CsDrip1. (A,B) Nucleotide and deduced amino acid sequence of CsDrip1 variant A (JQ011314) and variant B (JQ011315). Six transmembrane domains (TMs) are also shown by the red underline and C-terminal “SYDF” sequence found in insect AQPs with dark underline. The NPA boxes are shown with gray highlight. Orange indicates amino acids residues constituting the ar/R constriction region. Blue dots and green dots, respectively are potential phosphorylation sites and protein kinase C-specific sites. (C) Homology modeling of the CsDrip1 is compared with bovine aquaporin-1 (PDB ID: 1J4N) as template using Phyre2. Intracellular (B,D) and extracellular (A,C,E) loops and transmembrane helices are shown (TM1 in blue, TM2 in red, TM3 in yellow, TM4 in green, TM5 in purple, TM6 in orange). (D) The structure depicted is from the extracellular side of the membrane. The classical NPA motifs are shown in gray sphere, and Ar/R selectivity residues regions (Phe-69, His-194, Ser-203, and Arg-209) are shown in gray stick.

 Phylogenetic analysis

We used CLUSTALX and MEGA 7.0 phylogenetic analysis to compare CsDrip1 with other insect aquaporins. Because neighbor-joining, minimum evolution, maximum likelihood, and maximum parsimony methods constructed the similar phylogenetic tree, Figure 2 only illustrated the phylogenetic tree constructed by the neighbor joining method. The result exhibited the phylogenetic tree include two clusters: Drip1 and Prip. And CsDrip1 (C. suppressalis Drip1) was most closely related to BmDrip1 (B. mori Drip1) and DpDrip1 (Danaus plexippus Drip1), to which they were 79% identical at the amino acid level. Our phylogenetic tree showed BmDrip1, DpDrip1, and CsDrip1 belonged to the same group, which was water-specific Drip subfamily. All the Lepidoptera except BmPrip (B. mori Prip) fell into the well-supported cluster, which was consistent with a previous publication about BmPrip classification in the Prip subfamily (Azuma et al., 2012).

Neighbour-joining phylogenetic tree showing the Chilo suppressalis Drip1 (JQ011314 and JQ011315) with respect to other insect AQPs. BmDrip1 (Bombyx mori, AB178640), CmDrip1 (Cydia molesta, AB469882), AaDrip1 (Aedes aegypti, AF218314), PhDrip1 (Pediculus humanus, DS235154), TcDrip1 (Tribolium castaneum, CM000285), LlDrip1 (Lutzomyia longipalpis, EU124628), AgDrip1 (Anopheles gambiae, AAAB01008984), Cf Drip1 (Coptotermes formosanus, AB433197), HieDrip1 (Haematobia irritans exigua, U51638), CvDrip1 (Cicadella viridis, X97159), DpDrip1 (Drosophila pseudoobscura, CM000071), BtDrip1 (Bemisia tabaci, EU127479), GmDrip1 (Glossina morsitans, EZ422826), CqDrip1 (Culex quinquefasciatus, DS232465), ApDrip1 (Acyrthosiphon pisum, EU925136), D. plexippusDrip1 (Danaus plexippus, EHJ75085), EoDrip1 (Ectropis oblique, AMO02270.1), SlDrip1 (Spodoptera litura, AGT50942.1), ObDrip1 (Operophtera brumata, KOB73271.1), PrDrip1 (Phormia regina, BAM26200.1), MdDrip1 (Musca domestica, AIB09141.1), LhDrip1 (Lygus hesperus, AHI85743.1), EsPrip (Eurosta solidaginis, FJ489680), BgPrip (Blattella germanica, FR744897), Cf Prip (Camponotus floridanus, GL435707), BaPrip (Belgica Antarctica, AB602341), PvPrip (Polypedilum vanderplanki, AB281619), and BmPrip (Bombyx mori, AB458833). Numbers on the branches are the bootstrap values obtained from 1,000 replicates (only bootstrap values >50 are shown).

 Functional assay

In order to confirm further that CsDrip1 was a water-selective aquaporin, the permeability of CsDrip1 to water, glycerol, trehalose and urea transport was performed using the Xenopus oocyte expression system, respectively. Compared to control oocytes, oocytes expressing CsDrip1 showed an 11-fold increase in the osmotic water permeability coefficient (Pf) (N = 10) (Figure 3A), and as a result of continuous uptake, even some tested oocytes broke. For glycerol, trehalose and urea uptake in X. laevis oocytes expressing CsDrip1, no significant (p > 0.05) uptakes were observed in oocytes expressing CsDrip1, and these oocytes showed a slight shrinkage. For example, the Pf of CsDrip1 oocytes in glycerol, trehalose and urea was 5.031e−6 ± 5.766e−7 cm/s (N = 9), 5.670e−6 ± 6.0758e−7 cm/s (N = 10), and 5.029e−6 ± 3.542e−7 cm/s (N = 11), respectively (Figures 3B–D). These results clearly indicated that CsDrip1 was a specific water-selective aquaporin.

Functional assay of Chilo suppressalis aquaporin (CsDrip1). (A) Water permeability of CsDrip1. (B) Glycerol permeability of CsDrip1. (C) Trehalose permeability of CsDrip1. (D) Urea permeability of CsDrip1. The graphs of oocytes assayed in the experiments were indicated above each bar. Error bars represented SEM (n ≥ 9). Student's t-test was performed to find means that were significantly (p < 0.05) different from the water-injected control oocytes (indicated with asterisk).

 Expression of <italic>Cs</italic>DRIP1 in tissues or organs of <italic>C. suppressalis</italic> larvae

Real-time PCR verified that CsDrip1 mRNA was expressed in head, epidermis, fat body, foregut, midgut, hindgut, and Malpighian tubules [F(6, 14) = 4.441, P = 0.010], but there were not expression level in haemocytes. The highest level of CsDrip1 mRNA was observed within hindgut, which was 19.23-fold higher than that in the foregut. And CsDrip1 mRNA in Malpighian tubules also exhibited the second high abundance. Interestingly, the head was one of the organs that expressed a high level of CsDrip1 mRNA (Figure 4).

Relative mRNA expression levels of the CsDrip1 gene in different tissues (organs) of the Chilo suppressalis. HG, hindgut; MG, midgut; FG, foregut; HC, haemocytes; MT, Malpighian tubules; EP, epidermis; FB, fat body and HE, Heads. Values are denoted as the mean ± SE. Significant differences are indicated as P < 0.05.

 Expression of <italic>Cs</italic>DRIP1 in developmental stages of <italic>C. suppressalis</italic>

We investigated the expression patterns of CsDrip1 transcripts over the life cycle of C. suppressalis, including eggs, larvae (1st, 2nd, 3rd, 4th, and 5th instars larvae), pupae (male and female) and adults (male and female). The results demonstrated that the highest mRNA level of CsDrip1 was observed in the third instar larvae, which was 117.89-fold higher than that of male pupae, which followed by the first instar larvae, male adults, and egg mass. CsDrip1 transcript of the male pupae was the least. In addition, CsDrip1 mRNA exhibited significantly different level in male pupae and male adults than those of female pupae and female adults, respectively [F(9, 19) = 23.954, P < 0.001] (Figure 5).

Relative mRNA expression levels of the CsDrip1 gene in different developmental stages of the Chilo suppressalis. E, eggs; L1, the first instar larvae; L2, the second instar larvae; L3, the third instar larvae; L4, the fourth instar larvae; L5, the fifth instar larvae; MP, male pupae; FP, female pupae; MA, male adults; and FA, female adults. Values are denoted as the mean ± SE. Significant differences are indicated as P < 0.05.

 Expression of <italic>Cs</italic>Drip1 under various temperatures

Although the survival rate of larvae of C. suppressalis only reached to 55.56% at −11°C and100% at 42°C, respectively (Lu et al., 2014), CsDrip1 mRNA in larvae exposed to low temperatures for 2 h was not different significantly [F(7, 16) = 39.957, P = 0.560] (Figure 6). However, CsDrip1 displayed a different expressional pattern under heat stress. For example, CsDrip1 mRNA reached maximum at 36°C, which was 3.75-fold of that at control temperature (27°C) [F(7, 16) = 39.957, P < 0.001]. Subsequently, the abundance of CsDrip1 mRNA in larvae decreased by the elevated temperature (Figure 6).

Relative expression levels of CsDrip1 gene in the Chilo suppressalis under temperature stress. The larvae were exposed for −11°C, −9°C, −3°C, 0°C, 27°C, 36°C, 39°C, or 42°C for 2 h. Values are denoted as the mean ± SE. Significant differences are indicated as P < 0.05.

 Expression of <italic>Cs</italic>Drip1 in developmental stages of <italic>C. suppressalis</italic> under different humidities

According to prior experiments, C. suppressalis could tolerate these humidity treatments. The expression levels of CsDrip1 mRNA of the third instar larvae, fifth instar larvae, pupae (male and female) and adults (male and female) of C. suppressalis under different humidities have been determined. The results demonstrated that after exposure of the third instar larvae, fifth instar larvae, pupae (male and female) to 25–95% RH, CsDrip1 mRNA levels were not significantly different from each other [F(3, 12) = 1.309, P = 0.317; F(3, 8) = 1.285, P = 0.324; F(3, 8) = 3.453, P = 0.072; F(3, 8) = 0.964, P = 0.455] (Figures 7A–D). However, the CsDrip1 mRNA of female and male adults exposed to different humidities was changed significantly, and the CsDrip1 mRNA of female adults was up-regulated remarkably, which was contrary to male adults [F(3, 8) = 4.272, P = 0.045; F(3, 10) = 3.816, P = 0.047] (Figures 7E,F).

Relative expression levels of CsDrip1 gene in the Chilo suppressalis under humidity treatments. (A) The third instar larvae treated by 25, 50, 75, and 95% for 12 h, respectively. (B) The fifth instar larvae treated by 25, 50, 75, and 95% for 24 h, respectively. (C) The female pupae treated by 25, 50, 75, and 95% for 24 h, respectively. (D) The male pupae treated by 25, 50, 75, and 95% for 24 h, respectively. (E) The female adults treated by 25, 50, 75, and 95% for 24 h, respectively. (F) The male adults treated by 25, 50, 75, and 95% for 24 h, respectively. Values are denoted as the mean ± SE. Significant differences are indicated as P < 0.05.

 Discussion

Water-selective aquaporins are integral membrane proteins belonging to a large family of water channel proteins that assist in rapid movement of water across cellular membranes. This report was the first extensive analysis of Drip1 and its functions in C. suppressalis, an important rice pest of Lepidoptera. We isolated and characterized the first AQP from the C. suppressalis (CsDrip1). Deduced amino acid sequence of CsDrip1 resembled the features of other insect Drip1 (Ibanez et al., 2014). For example, two inverted hemi-helices on loops B and E that project opposing NPA motifs, which regulate the conductance of water (Wree et al., 2011). However, some insect AQPs lacked the first NPA motif, lost the ability to transport water but possessed the capacity of the glycerol transporter (Liu et al., 2016). Additionally, it had seven potential phosphorylation sites (Tyr5, Ser15, Ser16, Ser176, Ser210, Ser250, and Ser253), and three potential protein kinase C (PKC) specific sites (Thr3, Ser21, and Ser255) which somewhat contrasted to those of B. mori Drip1 (BmDrip1) that had only one PKC specific sites (Thr6) and one tyrosine kinase phosphorylation site (Tyr250) (Kataoka et al., 2009a). Although CsDrip1 variant A and variant B possessed the same ORF, the full-length sequence of variant A was 409 bp shorter than variant B (Figures 1A,B). Similar phenomenon was found in Prip of Anopheles gambiae, which included two variants, but two variants encoded different amino acids (Tsujimoto et al., 2013). It was also found that in Belgica antarctica there were three types of Prip (variants A-C) derived from the same gene, which was suggested to be due to alternative splicing (Goto et al., 2011). Hydrophobicity and structural prediction indicated that CsDrip1 possessed the conserved feature of water-specific AQPs. The pore of vertebrate AQPs is restricted by four residues (Phe58, His182, Cys191, and Arg197) that comprise the Ar/R constriction site (Sui et al., 2001; Horsefield et al., 2008; Hoa et al., 2009). Although the Ar/R constriction residues are generally conserved, Cys191 in vertebrate AQP is replaced by serine (Ser203 in TM5) in CsDrip1 (Figures 1C,D). Further comparative analysis of the vertebrate and insect AQPs revealed that either alanine or serine in all known insects Drip1 was substituted for Cys191 in vertebrate AQP1, and all Drip1 from Lepidoptera had a serine residue.

Aquaporins are most highly expressed in tissues where water movement is frequent and/or physiologically important. CsDrip1 transcripts showed specific expression patterns in various tissues (Figure 4). Water homeostasis in insects is achieved by a two-part system composed of Malpighian tubules and hindgut, and Malpighian tubules are the primary excretory and osmoregulatory organ in insects, analogous to the vertebrate renal tubules (Chawn and Nicolson, 2004). CsDrip1 mRNA was expressed abundantly in hindgut and Malpighian tubules, which was similar to the results observed in several other insects (Kambara et al., 2009; Kataoka et al., 2009a,b; Drake et al., 2010, 2015; Goto et al., 2011; Liu et al., 2011; Nagae et al., 2013; Fabrick et al., 2014). Fat body was one of the most cold-resistant tissues in the C. suppressalis because glycerol was accumulated in the fat body of cold acclimated and diapause larvae (Izumi et al., 2006, 2007). And our results revealed that CsDrip1 mRNA in fat body was high (Figure 4). However, low abundance or absence of Drip1 in the fat body had also been observed in several insect species (Kambara et al., 2009; Liu et al., 2011, 2013; Philip et al., 2011). Fat body was the site of glycerol biosynthesis (Kukal et al., 1988), which likely requires other kind of aquaporins, aquaglyceroporins for glycerol transportation as suggested in B. mori and G. molesta (Kataoka et al., 2009a,b). Therefore, CsDrip1 in the diapausing larvae of C. suppressalis couldn't contribute its cold hardiness. In summary, an unexpected diversity of AQPs was found in insect, and different AQPs contributed various functions. For example, there were six AQP genes in Ae. Aegypti. However, AQP 1 and AQP 2 were the strict water channels, which raised the of water permeability of midgut and Malpighian tubules. And AQP 5 demonstrated significant solute permeability for trehalose, which was important for insect temperature and dehydration tolerance (Drake et al., 2010, 2015; Van Ekert et al., 2016).

The levels of CsDrip1 transcript are high at early stages of development, and the expression of CsDrip1 was quite dynamic throughout development (Figure 5). CsDrip1 mRNA was highest in the third larvae of C. suppressalis, which just entered into the guzzled stage. The eggs of insect generally require high moist environment condition for hatching. High abundance of CsDrip1 mRNA in eggs of C. suppressalis was propitious to maintain the water balance of the embryonic development. However, BtDrip1 transcripts were most highly expressed in 2nd instar nymphs and least present in eggs (Mathew et al., 2011). AgDrip1 mRNA was higher in female adults than that in male adults (Liu et al., 2011), whereas CsDrip1 mRNA in male adults was significantly higher than in female adults. In addition, CsDrip1 mRNA was also greater in adults than that in pupae (Figure 5). We speculated that CsDrip1 expression might be related to reproduction.

When exposed to extreme temperatures, insects may respond in different ways: they can adopt behaviors to avoid or escape extreme temperatures, or they may regulate various proteins in response to adverse temperatures (Tursman et al., 1994; Rinehart and Denlinger, 2000; Tyshenko et al., 2005; Huang et al., 2007; Cui et al., 2010). In E. solidaginis by increasing the number of AQPs, cells likely improved their ability to rapidly redistribute water, better protecting themselves against the build-up of osmotic pressures across the membrane during freezing in winter (Philip et al., 2008, 2011). In Megaphorura arctica, up- or down-regulation of AQPs were contributed to exploit cryoprotective dehydration to enhance its cold tolerance (Clark et al., 2009). The AQPs of B. Antarctica also play very important roles in its freeze tolerance (Yi et al., 2011). And the studies suggested that the aquaporin might contribute to the cold tolerance of C. suppressalis (Izumi et al., 2006, 2007). The CsDrip1 might play the important role in the cold hardiness and diapause initiation of C. suppressalis (Lu et al., 2013). However, the CsDrip1 transcript in our studies was not up-regulated significantly under low temperature stress, but the CsDrip1 transcript increased significantly under high temperature stress (Figure 6). At the same time, our results confirm that CsDrip1 is a strict water channel. Therefore, C. suppressalis could utilize the up-regulation of CsDrip1 to exchange water in order to escape heat stress, but survival at low temperature they might largely depend on the ability of cells to accumulate the cryoprotectants among cellular compartments. Maybe, other AQPs existing in C. suppressalis could be coordinated to resist to the low stress.

Although the Xenopus oocyte expression system clearly demonstrated that CsDrip1 allowed water, but not glycerol, trehalose or urea, to pass through the cell, CsDrip1 mRNA of the third instar larvae, fifth instar larvae, and pupae (male and female) of C. suppressalis after exposure to 25–95% RH was not significantly regulated. However, CsDrip1 mRNA in female adults was greatly induced while expression of CsDrip1 mRNA in male adults was suppressed (Figure 7). Our data are inconsistent with that from a previous study in C. pipiens, which found that Prip mRNA level was significantly down-regulated in response to a low relative humidity (Liu and Piermarini, 2017). In the ovary of B. antarctica under water stress, no significant differences were observed in the levels of BaPrip mRNA (Goto et al., 2011). In response to dehydration, expression of PvDrip1 of Polypedilum vanderplanki larvae was greatly induced (Kikawada et al., 2008). Therefore, when confronted with water stress, different kinds of insects could possess different strategies by regulating different AQPs.

It is widely known that most insects possess multiple AQP genes. At least seven putative AQPs in A. gambia, six in Ae. Aegypti, eight in Drosophila, five in Lygus hesperus, eight in Bemisia tabaci, and seven in Triboleum castaneum had been identified (Adams et al., 2000; Holt et al., 2002; Drake et al., 2010; Fabrick et al., 2014; Van Ekert et al., 2016). And three AQPs in B. mori and two in G. molesta from Lepidoptera also had been isolated (Kataoka et al., 2009a,b; Azuma et al., 2012). The fat body of C. suppressalis could contain aquaglyceroporins related to the diapause and cold hardness (Izumi et al., 2006, 2007). Therefore, it's significant to identify and further analyse the other AQPs in C. suppressalis. We expect that the research of CsAQPs will be key in revealing the plot underlying mechanism of development, temperature tolerance, diapause and distribution of C. suppressalis in the future.

 Author contributions

Y-ZD and G-RW conceived and designed the experiments. M-XL, YL, D-DP, and JX preformed the experiments. M-XL and Y-ZD analyzed the data and wrote the manuscript.

Conflict of interest statement

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.

We would like to thank Jian-Yong Li for general advice and review; Xiao-Tian Tang, Wei-Fei Yan, and Qiu-Yu Li for insect feeding.

 References Adams M. D. Celniker S. E. Holt R. A. Evans C. A. Gocayne J. D. Amanatides P. G. . (2000). The genome sequence of Drosophila melanogaster. Science 287, 2185–2195. 10.1126/science.287.5461.218510731132 Agre P. (2006). The aquaporin water channels. Proc. Am. Thorac. Soc. 3, 5–13. 10.1513/pats.200510-109JH16493146 Agre P. Saboori A. M. Asimos A. Smith B. L. (1987). Purification and partial characterization of the Mr 30, 000 integral membrane protein associated with the erythrocyte Rh(D) antigen. J. Biol. Chem. 262, 17497–17503. 3121599 Azuma M. Nagae T. Maruyama M. Kataoka N. Miyake S. (2012). Two water-specific aquaporins at the apical and basal plasma membranes of insect epithelia: molecular basis for water recycling through the cryptonephric recta complex of Lepidopteran larvae. J. Insect Physiol. 58, 523–533. 10.1016/j.jinsphys.2012.01.00722285686 Benga G. (2009). Water channel proteins (later called aquaporins) and relatives: past, present, and future. IUBMB Life 61, 112–133. 10.1002/iub.15619165894 Benoit J. B. Hansen I. A. Attardo G. M. Mireji P. O. Bargul J. L. Drake L. L. . (2014a). Aquaporins are critical for provision of water during lactation and intrauterine progeny hydration to maintain tsetse fly reproductive success. PLoS Neg. Trop. Dis. 8, e2517. 10.1371/journal.pntd.0002517 Benoit J. B. Hansen I. A. Szuter E. M. Drake L. L. Burnett D. L. Attardo G. M.. (2014b). Emerging roles of aquaporins in relation to the physiology of blood-feeding arthropods. J. Comp. Physiol. B 184, 811–825. 10.1007/s00360-014-0836-x24942313 Blom N. Gammeltoft S. Brunak S. (1999). Sequence and structure-based prediction of eukaryotic protein phosphorylation sites. J. Mol. Biol. 294, 1351–1362. 10.1006/jmbi.1999.331010600390 Blom N. Sicheritz-Ponten T. Gupta R. Gammeltoft S. Brunak S. (2004). Prediction of post-translational glycosylation and phosphorylation of proteins from the amino acid sequence. Proteomics 4, 1633–1649. 10.1002/pmic.20030077115174133 Bustin S. A. Benes V. Garson J. A. Hellemans J. Huggett J. Kubista M. . (2009). The MIQE guidelines: minimum information for publication of quantitative Real-Time PCR experiments. Clin. Chem. 55, 611–622. 10.1373/clinchem.2008.11279719246619 Campbell E. M. Ball A. Hoppler S. Bowman A. S. (2008). Invertebrate aquaporins: a review. J. Comp. Physiol. B 178, 935–955. 10.1007/s00360-008-0288-218594835 Chang H. Liu Y. Yang T. Pelosi P. Dong S. Wang G. (2015). Pheromone binding proteins enhance the sensitivity of olfactory receptors to sex pheromones in Chilo suppressalis. Sci. Rep. 27, 5. 10.1038/srep13093 Chawn S. L. Nicolson S. W. (2004). Insect Physiological Ecology. New York, NY: Oxford University Press. Clark M. S. Thorne M. A. S. Purać J. Burns G. Hillyard G. Popović Ž. D. . (2009). Surviving the cold: molecular analyses of insect cryoprotective dehydration in the Arctic springtail Megaphorura arctica (Tullberg). BMC Genomics 10:328. 10.1186/1471-2164-10-32819622137 Cui Y. D. Du Y. Z. Lu M. X. Qiang C. K. (2010). Cloning of the heat shock protein 60 gene from the stem borer, Chilo suppressalis, and analysis of expression characteristics under heat stress. J. Insect Sci. 10, 1–13. 10.1673/031.010.1000120673188 Dow J. A. T. Kelly D. C. Davies S. A. Maddrell S. H. P. Brown D. (1995). A novel member of the major intrinsic protein family in Drosophila-aquaporins involved in insect Malpighian (Renal) tubule fluid secretion. J. Physiol. 489, 110–111. Drake L. L. Boudko D. Y. Marinotti O. Carpenter V. K. Dawe A. L. Hansen I. A. (2010). The aquaporin gene family of the yellow fever mosquito, Aedes aegypti. PLoS ONE 5:e15578. 10.1371/journal.pone.001557821249121 Drake L. L. Rodriguez S. D. Hansen I. A. (2015). Functional characterization of aquaporins and aquaglyceroporins of the yellow fever mosquito, Aedes aegypti. Sci. Rep. 5:7795. 10.1038/srep0779525589229 Duchesne L. Hubert J. F. Verbavatz J. M. Thomas D. Pietrantonio P. V. (2003). Mosquito (Aedes aegypti) aquaporin, present in tracheolar cells, transports water, not glycerol, and forms orthogonal arrays in Xenopus oocyte membranes. Eur. J. Biochem. 270, 422–429. 10.1046/j.1432-1033.2003.03389.x Echevarria M. Ramirez-Lorca R. Hernández C. S. Gutiérriez A. Méndez-Ferrer S. Toledo-Aral J. J. . (2001). Identification of a new water channel (Rp-MIP) in the Malpighian tubules of the insect Rhodnius prolixus. Pflvgers Arch. 442, 27–34. 10.1007/s00424000049411374065 Fabrick J. A. Pei J. Hull J. J. Yool A. J. (2014). Molecular and functional characterization of multiple aquaporin water channel proteins from the western tarnished plant bug, Lygus hesperus. Insect Biochem. Mol. Biol. 45, 125–140. 10.1016/j.ibmb.2013.12.00224333473 Fu D. Libson A. Miercke L. J. W. Weitzman C. Nollert P. Krucinski J. . (2000). Structure of a glycerol-conducting channel and the basis for its selectivity. Science 290, 481–486. 10.1126/science.290.5491.48111039922 Finn R. N. Cerda J. (2015). Evolution and functional diversity of aquaporins. Biol. Bull. 229:6e23. 10.1086/BBLv229n1p626338866 Goto S. G. Philip B. N. Teets N. M. Kawarasaki Y. Lee R. E. Jr. Denlinger D. L. (2011). Functional characterization of an aquaporin in the Antarctic midge Belgica Antarctica. J. Insect Physiol. 57, 1106–1114. 10.1016/j.jinsphys.2011.03.02321497603 Herraiz A. Chauvigné F. Cerdà J. Bellés X. Piulachs M. D. (2011). Identification and functional characterization of an ovarian aquaporin from the cockroach Blattella germanica L. (Dictyoptera, Blattellidae). J. Exp. Bio. 214, 3895–3903. 10.1242/jeb.06700921993792 Heymann J. B. Engel A. (1999). Aquaporins: phylogeny, structure, and physiology of water channels. News Physiol. Sci. 14, 187–193. 10.1152/physiologyonline.1999.14.5.18711390849 Hoa J. D. Yeh R. Sandstromb A. Chorny L. Harries W. E. C. Robbins R. A. . (2009). Crystal structure of human aquaporin 4 at 1.8 Å and its mechanism of conductance. Proc. Natl. Acad. Sci. U.S.A. 106, 7437–7442. 10.1073/pnas.0902725106 Holt R. A. Subramanian G. M. Halpern A. Sutton G. G. Charlab R. Nusskern D. R. . (2002). The genome sequence of the malaria mosquito Anopheles gambiae. Science 298, 129–149. 10.1126/science.107618112364791 Horsefield R. Nordén K. Fellert M. Backmark A. Törnroth-Horsefield S. Terwisscha van Scheltinga A. C. . (2008). High-resolution x-ray structure of human aquaporin 5. Proc. Natl. Acad. Sci. U.S.A. 105, 13327–13332. 10.1073/pnas.080146610518768791 Huang L. H. Chen B. Kang L. (2007). Impact of mild temperature hardening on thermotolerance, fecundity, and Hsp gene expression in Liriomyza huidobrensis. J. Insect Physiol. 53, 1199–1205. 10.1016/j.jinsphys.2007.06.01117651748 Ibanez F. Hancock J. Tamborindeguy C. (2014). Identification and expression analysis of aquaporins in the potato psyllid, Bactericera cockerelli. PLoS ONE 9:e111745. 10.1371/journal.pone.011174525354208 Ishibashi K. (2006). Aquaporin subfamily with unusual NPA boxes. Bioch. Biophy. Acta 1758, 989–993. 10.1016/j.bbamem.2006.02.02416579962 Izumi Y. Sonoda S. Tsumuki H. (2007). Effects of diapause and cold-acclimation on the avoidance of freezing injury in fat body tissue of the rice stem borer, Chilo suppressalis Walker. J. Insect Physiol. 53, 685–690. 10.1016/j.jinsphys.2007.04.00517543330 Izumi Y. Sonoda S. Yoshida H. (2006). Role of membrane transport of warter and glycerol in the freeze tolerance of the rice stem borer, Chilo suppressalis Walker (Lepidoptera: Pyralidae). J. Insect Physiol. 52, 215–220. 10.1016/j.jinsphys.2005.11.001 Kambara K. Takematsu Y. Azuma M. Kobayashi J. (2009). cDNA cloning of aquaporin gene expressed in the digestive tract of the Formosan subterranean termite, Coptotermes formosanus Shiraki (Isoptera: Rhinotermitidae). Appl. Entomol. Zool. 44, 315–321. 10.1303/aez.2009.315 Kataoka N. Miyake S. Azuma M. (2009a). Aquaporin and aquaglyceroporin in silkworms, differently expressed in the hindgut and midgut of Bombyx mori. Insect Mol. Biol. 18, 303–314. 10.1111/j.1365-2583.2009.00871.x19366375 Kataoka N. Miyake S. Azuma M. (2009b). Molecular characterization of aquaporin and aquaglyceroporin in the alimentary canal of Grapholita molesta (the oriental fruit moth) -comparison with Bombyx mori aquaporins. J. Insect Biotechnol. Sericol. 78, 81–90. 10.11416/jibs.78.2_81 Kaufmann N. Mathai J. C. Hill W. G. Dow J. A. T. Zeidel M. L. Brodsky J. L. (2005). Developmental expression and biophysical characterization of a Drosophila melanogaster aquaporin. Am. J. Physiol. Cell Phys. 289, C397–C407. 10.1152/ajpcell.00612.200415800049 Kelley L. A. Sternberg M. J. E. (2009). Protein structure prediction on the web: a case study using the phyre server. Nat. Protoc. 4, 363–371. 10.1038/nprot.2009.219247286 Kikawada T. Saito A. Kanamori Y. Fujita M. Snigórska K. Watanabe M. Okuda T. . (2008). Dehydration-inducible changes in expression of two aquaporins in the sleeping Chironomid, Polypedilum vanderplanki. Bioch. Biophy Acta 1778, 514–520. 10.1016/j.bbamem.2007.11.00918082130 Kukal O. Serianni A. S. Duman J. G. (1988). Glycerol metabolism in a freezetolerant arctic insect: an in vivo 13C NMR study. J. Comp. Physiol. B 158, 175–183. 10.1007/BF010758313170824 Kumar S. Stecher G. Tamura K. (2016). MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. 10.1093/molbev/msw05427004904 LeCahérec F. Bron P. Verbavatz J. M. Garret A. Morel G. Cavalier A. . (1996a). Incorporation of proteins into (Xenopus) oocytes by proteolipo-some microinjection: functional characterization of a novel aquaporin. J. Cell Sci. 109, 1285–1295. LeCahérec F. Deschamps S. Delamarche C. Pellerin I. Bonnec G. Guillam M. T. . (1996b). Molecular cloning and characterization of an insect aquaporin-functional comparison with aquaporin 1. Eur. J. Biochem. 241, 707–715. 10.1111/j.1432-1033.1996.00707.x Li L. D. Wang D. F. Liu F. J. Li H. L. Zhang H. Wu G. Y. (2016). cDNA cloning, preparation of polyclonal antibody and subcellular localization of aquaporin 1 (AQP1) in Ectropis obliqua (Lepidoptera: Geoqmetridae). Acta Entomol. Sin. 59, 382–391. 10.16380/j.kexb.2016.04.003 Liu H. Y. Shu B. S. Jiang C. L. Li L. D. Zhong G. H. (2013). Molecular cloning, characterization and expression analysis of aquaporin1 (AQP1) gene in Spodoptera litura (Lepidoptera: Noctuidae). Acta Entomol. Sin. 56, 339–349. Liu K. Tsujimoto H. Cha S. J. Agre P. Rasgon J. L. (2011). Aquaporin water channel AgAQP1 in the malaria vector mosquito Anopheles gambiae during blood feeding and humidity adaptation. Proc. Natl. Acad. Sci. U.S.A. 108, 6062–6066. 10.1073/pnas.110262910821444767 Liu K. Tsujimoto H. Huang Y. Z. Rasgon J. L. Agre P. (2016). Aquaglyceroporin function in the alaria mosquito Anopheles gambiae. Bio. Cell 108, 294–305. 10.1111/boc.20160003027406921 Liu Y. Piermarini P. M. (2017). Molecular expression of aquaporin mRNAs in the northern house mosquito, Culex pipiens. J. Insect Physiol. 96, 35–44. 10.1016/j.jinsphys.2016.10.011 Lu M. X. Du Y. Z. Liu Z. X. Hua J. Liu P. Y. Li J. Y. (2013). Diapause, signal and molecular characteristics of overwintering Chilo suppressalis. Sci. Rep. 3:3211. 10.1038/srep0321124226906 Lu M. X. Liu Z. X. Cui Y. D. Du Y. Z. (2014). Expression patterns of three heat shock proteins in Chilo suppressalis (Lepidoptera: Pyralidae). Ann. Entomol. Soc. Am. 107, 667–673. 10.1603/AN13092 Lu M. X. Liu Z. X. Wang X. Du Y. Z. (2012). Seasonal cold tolerance of Chilo suppressalis (Walker) in Yangzhou, China. Ann. Entomol. Soc. Am. 105, 479–483. 10.1603/AN11171 Mathew L. G. Campbell E. M. Yool A. J. Fabrick J. (2011). A identification and characterization of functional aquaporin water channel protein from alimentary tract of whitefly, Bemisia tabaci. Insect Biochem. Mol. Biol. 41, 178–190. 10.1016/j.ibmb.2010.12.00221146609 Maurel C. Verdoucq L. Luu D. T. Santoni V. (2008). Plant aquaporins: membrane channels with multiple integrated functions. Annu. Rev. Plant Biol. 59, 595–624. 10.1146/annurev.arplant.59.032607.09273418444909 Murata K. Mitsuoka K. Hirai T. Walz T. Agre P. Heymann J. B. . (2000). Structural determinants of water permeation through aquaporin-1. Nature 407, 599–605. 10.1038/3503651911034202 Nagae T. Miyake S. Kosaki S. Azuma M. (2013). Identification and characterisation of a functional aquaporin water channel (Anomala cuprea DRIP) in a coleopteran insect. J. Exp. Biol. 216, 2564–2572. 10.1242/jeb.08338623531819 Nolan T. Hand R. E. Bustin S. A. (2006). Quantification of mRNA using real-time RT-PCR. Nat. Ptotoc, 1, 1559–1552. 10.1038/nprot.2006.23617406449 Pallant J. (2005). SPSS Survival Manual: A Step by Step Guide to Data Analysis using SPSS for Windows (Version 12). Maidenhead: Open University Press. Pettersen E. F. Goddard T. D. Huang C. C. Couch G. S. Greenblatt D. M. Meng E. C. . (2004). UCSF Chimera–a visualization system for exploratory research and analysis. J. Comput. Chem. 25, 1605–1612. 10.1002/jcc.2008415264254 Philip B. N. Kiss A. J. Lee R. E. (2011). The protective role of aquaporins in the freeze-tolerant insect Eurosta solidaginis: functional characterization and tissue abundance of EsAQP1. J. Exp. Bio. 214, 848–857. 10.1242/jeb.05127621307072 Philip B. N. Yi S. X. Elnitsky M. A. Lee R. E. (2008). Aquaporins play a role in desiccation and freeze tolerance in larvae of the goldenrod gall fly, Eurosta solidaginis. J. Exp. Bio. 211, 1114–1119. 10.1242/jeb.01675818344486 Preston G. M. Carroll T. P. Guggino W. B. Agre P. (1992). Appearance of water channels in Xenopus oocytes expressing red cell CHIP28 protein. Science 256, 385–387. 10.1126/science.256.5055.3851373524 Rinehart J. P. Denlinger D. L. (2000). Heat-shock protein 90 is down-regulated during pupal diapause in the flesh fly, Sarcophaga crassipalpis, but remains responsive to thermal stress. Insect Mol. Biol. 9, 641–645. 10.1046/j.1365-2583.2000.00230.x11122473 Rojek A. Praetorius J. Frokiaer J. Nielsen S. Fenton R. A. (2008). A current view of the mammalian aquaglyceroporins. Annu. Rev. Physiol. 70, 301–327. 10.1146/annurev.physiol.70.113006.10045217961083 Schmittgen T. D. Livak K. J. (2008). Analyzing real-time PCR data by the comparative CT method. Nat. Protoc. 3, 1101–1108. 10.1038/nprot.2008.73 Shakesby A. J. Wallace I. S. Isaacs H. V. Pritchard J. Roberts D. M. Douglas A. E. (2009). A water-specific aquaporin involved in aphid osmoregulation. Insect Biochem. Mol. Biol. 39, 1–10. 10.1016/j.ibmb.2008.08.00818983920 Shang Z. Z. Wang Y. S. Zou Y. H. (1979). Study on rearing method of rice stem borer Chilo suppressalis Walker. Acta Entomol. Sin. 2, 164–167. Spring J. H. Robichaux S. R. Hamlin J. A. (2009). The role of aquaporins in excretion in insects. J. Exp. Bio. 212, 358–362. 10.1242/jeb.02479419151210 Sui H. X. Han B. G. Lee J. K. Walian P. Jap K. B. (2001). Structural basis of water-specific transport through the AQP1 water channel. Nature 414, 872–878. 10.1038/414872a11780053 Tsujimoto H. Liu K. Linser P. L. Agre P. Lasgon J. L. (2013). Organ-specific splice variants of aquaporin water channel AgAQP1 in the malaria vector Anopheles gambiae. PLoS ONE 8:e75888. 10.1371/journal.pone.007588824066188 Tsumuki H. Hirai M. (2007). Effects of photoperiod and temperature on endogenous ice nucleus production in larvae of the rice stem borer, Chilo suppressalis Walker (Lepidoptera: Pyralidae). Appl. Entomol. Zool. 42, 305–308. 10.1303/aez.2007.305 Tursman D. Duman J. G. Knight C. A. (1994). Freeze tolerance adaptations in the centipede, Lithobius forficatus. J. Exp. Zool. 268, 347–353. 10.1002/jez.1402680503 Tyshenko M. G. Doucet D. Walker V. K. (2005). Analysis of antifreeze proteins within spruce budworm sister species. Insect Mol. Biol. 14, 319–326. 10.1111/j.1365-2583.2005.00562.x15926901 Van Ekert E. Chauvigné F. Finn R. N. Mathew L. G. Hull J. J. Cerd,à J. . (2016). Molecular and functional characterization of Bemisia tabaci aquaporins reveals the water channel diversity of hemipteran insects. Insect Bio. Mol. Biol. 77, 39–51. 10.1016/j.ibmb.2016.07.01027491441 van Hoek A. N. Hom M. L. Luthjens L. H. de Jong M. D. Dempster J. A. van Os C. H. (1991). Functional unit of 30 kDa for proximal tubule water channels as revealed by radiation inactivation. J. Biol. Chem. 266, 16633–16635. 1885592 Verkman A. S. (2005). More than just water channels: unexpected cellular roles of aquaporins. J. Cell Sci. 118, 3225–3232. 10.1242/jcs.0251916079275 Wree D. Wu B. Zeuthen T. Beitz E.. (2011). Requirement for asparagine in the aquaporin NPA sequence signature motifs for cation exclusion. FEBS J. 278, 740–748. 10.1111/j.1742-4658.2010.07993.x21205205 Xu J. Lu M. X. Cui Y. D. Du Y. Z. (2017). Selection and evaluation of reference genes for expression analysis using qRT-PCR in Chilo suppressalis (Lepidoptera: Pyralidae). J. Econ. Entomol. 110, 683–691. 10.1093/jee/tow29728115499 Yanochko G. M. Yool A. J. (2002). Regulated cationic channel function in Xenopus oocytes expressing Drosophila big brain. Int. J. Neurosci 22, 2530–2540. 11923418 Yasui M. (2004). Molecular mechanisms and drug development in aquaporin water channel diseases: structure and function of aquaporins. J. Pharm. Sci. 96, 260–263. 10.1254/jphs.FMJ04004X415557734 Yi S. X. Benoit J. B. Elnitsky M. A. Kaufmann N. Brodsky J. L. Zeidel M. L. . (2011). Function and immune-localization of aquaporins in the Antarctic midge Belgica antarctica. J. Insect Physiol. 57, 1096–1105. 10.1016/j.jinsphys.2011.02.006 Yin C. L. Liu Y. Liu J. D. Xiao H. M. Huang S. Q. Lin Y. J. . (2014). ChiloDB: a genomic and transcriptome database for an important rice insect pest Chilo suppressalis. Database 15, 92–108. 10.1093/database/bau065 Zhang R. B. Logee K. A. Verkman A. S. (1990). Expression of mRNA coding for kidney and red cell water channels in Xenopus oocytes. J. Biol. Chem. 265, 15375–15378. 2394728

Funding. This work was supported by the National Natural Science Foundation of China (31371937, 31401733) and the State Key Laboratory for Biology of Plant Diseases and Insect Pests (SKLOF201406).