- State Key Laboratory of Crop Stress Biology for Arid Areas, Northwest A&F University, Yangling, China
In aphids there is a fecundity-dispersal trade-off between wingless and winged morphs. Recent research on the molecular mechanism of wing morphs associated with dispersal reveals that insulin receptors in the insulin signaling (IS) pathway regulate alternation of wing morphs in planthoppers. However, little is known about whether genes in the IS pathway are involved in developmental regulation in aphid nymphs with different wing morphs. In this study, we show that expression of the insulin-related peptide 5 gene (Apirp5) affects biochemical composition and embryo development of wingless pea aphids, Acyrthosiphon pisum. After comparing expression levels of major genes in the IS pathway between third instar winged and wingless nymphs, we found that Apirp5 showed higher expression in head and thorax in the wingless nymphs than in the winged nymphs. Although microinjection treatment affects physical performance in aphids, nymphs with RNA interference of Apirp5 had less weight, smaller embryos, and higher carbohydrate and protein contents compared to the control group. Comparison between winged and wingless nymphs showed a similar trend. These results indicate that Apirp5 is involved in embryo development and metabolic regulation in wing dimorphic pea aphid.
Introduction
Fecundity-dispersal trade-offs strongly drive the life-history evolution of insects with wing polyphenism. Because of their short lifespan, allocation of limited nutritional resources differs for insects that engage in optimal reproduction or long-distance migration. Active dispersal can decrease the fecundity of insects (Bonte et al., 2012), whereas reproduction may limit migration during a lifecycle. The trade-off between fecundity and dispersal is best studied in species with the trait of wing dimorphism, including pea aphid (Acyrthosiphon pisum; Braendle et al., 2006; Brisson, 2010), sand cricket (Gryllus firmus; Mole and Zera, 1993), and planthopper (Prokelisia dolus; Denno et al., 1989). Typically, the flight-capable morph (winged or long-wing) has fully developed wings and muscle apparatus and undeveloped ovaries, but in contrast, ovaries grow rapidly and early in the flightless morph (wingless or short-wing), and these differences can be identified during nymphal stages. The choice to reproduce or disperse is usually governed by environmental conditions and/or genetic backgrounds. The molecular mechanism governing the wing morph has been revealed for the migratory brown planthopper (Nilaparvata lugens); two insulin receptors in the insulin signaling (IS) pathway play a regulatory role in governing wing morph determination (Xu et al., 2015). However, very little is known about the molecular mechanism of developmental and physiological changes in nymphs with wing polyphenism.
Photoperiod could affect the reproductive mode of aphids. Sexual reproduction occurs when the days grow shorter in the late summer, and eggs are produced by mated female aphids for overwintering. Under long day and warm conditions, developing embryos can be found in the ovarioles of nymphs and adult parthenogenetic viviparous aphids. The aphid also shows wing polyphenism, in which asexual aphids with an identical genetic background can produce winged or wingless progeny under long day conditions (Sutherland, 1969). Aphids with different wing morphs show a fecundity-dispersal trade-off and distinct growth rate in the same tissues (Dixon et al., 1993; Ishikawa et al., 2008; Ishikawa and Miura, 2009). Newly molted wingless adult aphids have highly developed ovaries for rapid reproduction, and winged adult aphids possess indirect flight muscles and wings for long distance migration (Ishikawa and Miura, 2009). Wing dimorphism is evident during the aphid nymphal phase. However, during the first and second instar stages, it is impossible to distinguish winged nymphs from wingless nymphs by wing apparatus morphology (Ishikawa et al., 2008). In the third instar, nymphs destined to be winged aphids have visible wing primordia on mesothorax, whereas no wing primordia are present on wingless nymphs (Ishikawa et al., 2008; Ogawa and Miura, 2013). Embryos in the ovary of wingless nymphs develop faster than those of winged aphids during the entire nymphal stage (Dixon and Howard, 1986; Ishikawa and Miura, 2009). Such morphological and physiological differences clearly separate the development of winged aphids from wingless aphids (Braendle et al., 2006; Simpson et al., 2011).
The insulin and insulin growth factor signaling pathways play diverse roles in regulation of growth, metabolism, reproduction, behavior, and aging in invertebrates and vertebrates (Garofalo, 2002; Taniguchi et al., 2006; Antonova et al., 2012). The insulin signaling (IS) pathway is activated by insulin-like peptides (ILPs, also called insulin-related peptides, IRPs) through the insulin receptor (InR), which activates a sequence of downstream proteins, such as the insulin receptor substrate (IRS). Signaling in this pathway can be downregulated by the insulin degrading enzyme (IDE), which degrades insulin intracellularly (Shen et al., 2006). Evidence shows that the IS pathway regulates alternative reproductive phenotypes in insects by controlling resource allocation during tissue development (Emlen et al., 2006; Wheeler et al., 2006; de Azevedo and Hartfelder, 2008). Amilp2 dsRNA treatment results in a reduction in ovary size in the honey bee (Wang et al., 2013), and ovariole number is decreased by irs dsRNA feeding (Wolschin et al., 2011). In another example, the FOXO gene, a downstream nexus in the IS pathway, directly regulates the size of copulatory organs in male dung beetle Onthophagus nigriventris, and is affected by nutritional polyphenism (Snell-Rood and Moczek, 2012).
The pea aphid A. pisum (Harris) is an important crop pest and a study model with wingless and winged morphs (Brisson and Stern, 2006; Figures S1A,B) that reflect a fecundity-dispersal trade-off. External environmental signals, including crowding, poor-quality diet, and predators, stimulate adult female aphids to produce nymphs that undergo wing differentiation and other related changes in organs and tissues (Sutherland, 1969; Brisson, 2010; Ogawa and Miura, 2013). How these signals are transduced internally in nymphs is not known. Based on the role of the IS pathway in regulating wing dimorphism in other insects, we first examined the conservation of ILPs and signaling proteins, as revealed in the pea aphid genome (Consortium, 2010). Previous studies identified 10 insulin-related peptide genes (Apirps), two insulin receptor genes (Apinrs), one insulin receptor substrate gene (Apirs), and two insulin degrading enzyme genes (Apides) in pea aphid (Consortium, 2010; Huybrechts et al., 2010; Antonova et al., 2012). The ApIRPs are divided into three groups based on cleavage sites and cysteine spacing (Huybrechts et al., 2010). The first group includes ApIRP1 to ApIRP4 with expected cysteine spacing and canonical prohormone convertase cleavage sites. The second group (ApIRP5, ApIRP6, and ApIRP7) is likely cleaved by furin. Notably, Apirp5 transcripts are highly represented in EST sequences, suggestive of high expression and an important role in rapid development of aphids (Huybrechts et al., 2010). The third group, including ApIRP8 to ApIRP10, has a different cysteine spacing from the first group (Huybrechts et al., 2010). We hypothesized that genes in the IS pathway would show different expression patterns, and that some of the genes were involved in regulation of development and biochemical components in aphid nymphs with or without wing primordia. In this study, we focused on the third instar pea aphid, because it is the earliest phase when winged nymphs can be identified by visual inspection. After comparing the transcript levels of genes in the IS pathway between third instar winged and wingless aphids, we selected Apirp5 for RNA interference (RNAi). Analysis of morphological and physiological changes of pea aphids after Apirp5 RNAi treatment revealed functions of Apirp5 on regulation of embryo development and biochemical components in wingless nymphs.
Methods
Pea Aphids
A. pisum used in this study was collected from Yunnan Province, China, and maintained on potted broad bean (Vicia faba L.) under long-day (LD) conditions (14 h light: 10 h dark, 20 ± 2°C).
Candidate Genes in the Insulin Signaling Pathway
Sequences of Apirps, Apinrs, Apirs, and Apides were obtained from Huybrechts et al. (2010) and from the National Center for Biotechnology Information (NCBI) via Blast searches of the pea aphid genome (www.ncbi.nlm.nih.gov/BLAST/) with homologous proteins sequences from Drosophila melanogaster and A. mellifera. The protein domains of predicted Apinr, Apirs, and Apide genes were analyzed and identified by searching the Self-Monitoring, Analysis and Reporting Technology database (http://smart.embl-heidelberg.de/). Multiple sequence alignments were performed by using MEGA 6.0 (Tamura et al., 2013).
Experimental Treatments
Winged aphids were induced under a high-density (HD) condition (30 wingless adult aphids per seedling) as described by Ishikawa et al. (2008) and Sutherland (1969). Thirty wingless adult aphids were caged on a broad bean seedling that was ~2 cm in height. After 24 h, adult aphids were removed, and the offspring were reared on the broad beans until they developed to third instar when winged nymphs show visible wing primordia (Ishikawa et al., 2008). The induction rate of winged aphids was 47.8 ± 5.3% (n = 836). In order to minimize the effects of development rate of winged and wingless nymphs on gene expression pattern and biochemical contents, we weighed and dissected third instar nymphs on the fifth day after birth because each of the first three instars in both winged and wingless aphids lasted for 2 days, and winged and wingless nymphs could not be visually identified until the third instar. Five days after removing adult aphids, all third instar nymphs were individually weighed, and 30–40 third instar aphids with or without wing primordia were dissected in 0.1 M phosphate-buffered saline (PBS, pH 7.4). The antenna, head without antenna, thorax, midgut, and abdomen without midgut were separately collected. For the treatments above, each had four independent replicates. To avoid the known effects of HD conditions on promoting the wing polyphenism (MacKay and Wellington, 1976), female asexual aphids were reared at low density (LD; three aphids per seedling) for longer than three generations before experiments.
RNA Extraction and cDNA Synthesis
Total RNA was extracted using RNAiso Plus reagent (TaKaRa, Dalian, China) following the manufacturer's protocol. The RNA was quantified using a NanoDrop® 2000c (Thermo Fisher Scientific, Middletown, VA, USA) at 260 nm. To generate the first strand cDNA, 1 μg of total RNA was reverse-transcribed in 20 μl volume using the PrimeScript RT reagent Kit With gDNA Eraser (TaKaRa, Dalian, China).
Real-Time Quantitative PCR
The reactions were performed with iQ5 real-time cycler (Bio-Rad, Hercules, CA, USA) and SYBR® Premix Ex Taq™ II (TaKaRa, Dalian, China). The primers to each target transcript were listed in Table 1. The qRT-PCR amplified fragment of Apirp5 did not overlap with Apirp5 dsRNA (Figure 1). Running parameters were 94°C for 5 min, followed by 40 cycles at 94°C for 30 s, 57°C for 30 s, and 72°C for 1 min. Considering the high homology present in several insulin genes, the real-time PCR products were extracted from 1% agarose gel and purified using a gel extraction kit (Axygen Scientific, Union City, CA, USA) for cloning into pGEM-T Easy Vector System (Promega, Madison, WI, USA) and submitted for sequencing (Life technology Inc., Shanghai, China). To ensure that only a single product was amplified, we carried out a melting curve analysis for each qRT-PCR reaction. Mean PCR efficiency values for each gene were calculated from individual amplification by using LinRegPCR software (Ramakers et al., 2003), which also provided the baseline and cycle threshold (ct). Relative gene expression data were analyzed following the Pfaffl method (Pfaffl, 2001) with normalization to the A. pisum ribosomal protein (Rp) L7 transcript (Nakabachi et al., 2005). Non-expression was defined as a relative quantity lower than 0.0001.
Table 1. List of accession numbers for genes encoding insulin signal components in the pea aphid, and primers for real-time quantitative PCR and dsRNA synthesis.
Figure 1. Comparison of Apirp5 mRNA sequences. Two transcript variants of Apirp5 (line 1 and line 2) share the same predicted ORF (box of line 3), and there is no overlap between the products of Apirp5 qRT-PCR and Apirp5 dsRNA (line 4 and line 5). In line 1 and line 2, the black lines show the same sequences in the two transcript variants, and gray solid line and gray dotted line show different sequences in the 5′UTR. In lines 4 and 5, dotted lines show the products of Apirp5 qRT-PCR and dsRNA, and short lines are primers of two products.
Double Stranded RNA (dsRNA) Synthesis and Delivery
Specific primers designed for dsApirp5 were used to amplify target sequence from cDNA derived from adult aphid heads (Table 1). Cycling conditions were 95°C for 5 min, followed by 30 cycles at 95°C for 30 s, 55°C for 30 s, 72°C for 30 s, and then 72°C for 5 min. The PCR product was extracted from 1% agarose gel and purified using a gel extraction kit (Axygen Scientific, Union City, CA, USA). The product was cloned (pGEM-T Easy Vector System, Promega, Madison, WI, USA) and submitted for sequencing (Life Technology Inc., Shanghai, China). To synthesize dsApirp5, two separate PCR reactions were performed using a gene specific primer and a primer with terminal 5′ T7 promoter site (Table 1). Amplification strategies were 95°C for 1 min, followed by initial 10 cycles at 95°C for 30 s, 55°C for 30 s, 72°C for 30 s, and then 25 cycles at 95°C for 30 s, 60°C for 30 s, 72°C for 30 s, and then 72°C for 5 min. Purified PCR products were used as templates for dsRNA synthesis using the T7 RiboMAX™ Express RNAi System (Promega, Madison, WI, USA). The control dsRNA was prepared using gfp gene with primers listed in Table 1.
About 80–100 second instar pea aphids were maintained under a LD condition for 24 h, and dsRNA (200 nl, 8 mg/ml) was injected into each aphid with an MM33 links microinjector (Märzhäuser, Wetzlar, Germany). The needles were prepared by a P-97 Flaming/Brown Micropipette Puller (Sutter Co, Novato, CA, USA) using Borosilicate glass capillaries. The precise injection volume was controlled by an MM33 micromanipulator (Märzhäuser, Wetzlar, Germany). Each injected aphid was reared on a leaf disc kept on a thin layer water-agar gel (0.9% agar) in a 24-well plate. Two days later, aphids molting into third instar nymphs were individually weighed using Cubis® Ultramicro Balance with 0.001 mg readability (Sartorius, Göttingen, Germany), and then 22 treated and 22 control aphids were frozen in liquid nitrogen, and kept at –80°C before use, and 15 treated and 15 control aphids were dissected. The treatments were repeated four times. The aphids that died during experiments were removed.
Size of Embryos
The embryo size was measured from four groups of third instar nymphs: aphids with or without wing primordia produced under the HD condition, and aphids injected with dsApirp5 or dsgfp. The aphids were dissected under a dissecting microscope, and the embryos surrounded by a very thin layer of ovaries were photographed through an EVOS FL microscope (Life Technology Inc., Gaithersburg, MD, USA). The length of the five largest embryos in the ovary of each aphid was measured using ImageJ and analyzed (Abràmoff et al., 2004; Ishikawa and Miura, 2009).
Protein, Glycogen, and Soluble Carbohydrate Assay
Protein, glycogen, and soluble carbohydrate were measured from third instar winged and wingless nymphs reared under HD condition, and also from third instar nymphs subjected to RNAi (dsgfp injected and dsApirp5 injected).
Eleven third instar nymphs with or without wing primordia were randomly collected from plants under HD inducing conditions for metabolic analysis. Third instar nymphs from dsRNA treatments were collected. Each whole aphid was measured individually as reported (Foray et al., 2012). After homogenizing in 180 μl lysis buffer solution (100 mM KH2PO4, 1 mM dithiothreitol, and 1 mM ethylenediaminetetraacetic acid, pH 7.4) at room temperature, 50 μl was removed to a new 1.5 ml microtube for protein assay following the method of Bradford (Kruger, 2009). The protein content was measured using a NanoDrop® 2000c (Thermo Fisher Scientific, Middletown, VA, USA) at 595 nm with bovine serum albumin as the standard. The rest was added with 37.5 μl of 20% Na2SO4, 112.5 μl of lysis buffer solution, 300 μl of chloroform, and 600 μl of methanol. After vigorous vortexing, the sample was centrifuged at 180 g for 15 min at 4°C to remove glycogen from the supernatant. One hundred and fifty microliters of supernatant were transferred into a new microtube for soluble carbohydrate analyses. The pellet was kept for glycogen content determination. All glycogen and soluble carbohydrate contents were determined by the anthrone method (Van Handel, 1965) and the absorbance was measured at 625 nm with D-glucose as the standard. The treatment was repeated three times.
Statistical Analyses
Two-way ANOVAs were used to determine whether the levels of gene expression were significantly different across the five tissues or parts of aphid, and whether this variation was affected by wing morph of nymphs. Then, multiple comparisons were corrected by Sidak test. Prior to comparing other values using Student's t-test, homogeneity of variance and normality assumptions (Shapiro-Wilk normality test) of the raw data were tested. Significance was ascribed at p < 0.05. All analyses were performed, and the data were graphed in GraphPad Prism (GraphPad Prism version 6.01 for Windows, GraphPad Software, San Diego, CA, USA). No statistical methods were used to predetermine sample size. Figures and schematic illustrations were prepared using Adobe Illustrator CS5 (Adobe Systems, San Jose, CA, USA).
Results
Comparison of Sequence Features for Selected Genes in the is Pathway
By comparing the products of real-time PCR with the NCBI database, we confirmed that the fragments amplified for all selected genes of A. pisum clone in this study were the same as sequences retrieved from pea aphid genome. In all retrieved Apirps, only Apirp5 was expressed in two transcript variants (Apirp5 transcript variant x1 and x2) which shared the same ORF (open reading frame; sequence information of Apirp5 is shown in Figure 1). Accession numbers of Apirp7, 9, and 10 mRNAs were not retrieved from NCBI, so the predicted mRNA sequences were spliced based on the pea aphid genome and their amino acid sequences (Table 1).
Two Apinr sequences were retrieved from the genome of A. pisum (Table 1), and one ApIRS and two ApIDEs were retrieved based on the amino acid sequences of D. melanogaster (Table 1). Apinr1, Apirs, and Apide2 were transcribed to three, four and two transcript variants, respectively, and then were translated into different protein isoforms (Table 1 and Figure S5).
The functional domains of deduced ApInRs, ApIRS, and ApIDEs of A. pisum were aligned with those of other model insects (Figures S2–S4).
Expression Profiles of Genes in the is Pathway in Third Instar Nymphs with Different Wing Morphs
The third instar is the earliest phase in which winged nymphs could be identified, so the expressions of insulin related genes in the IS pathway were compared between the third instar nymphs with different wing morphs. Aside from the signal of amplification of Apirp2, 3, 8, 9, and 10, other transcripts were detected.
All expression differences of Apirps between the two wing-morph nymphs showed the same trend: transcript expression was body-part specific and higher in wingless nymphs than in winged nymphs (Figures 2A–E). The expression of Apirp1 was identified in antenna and head, although there was a significant difference in the heads between the winged and the wingless nymphs (Figure 2A). The transcripts of Apirp4 existed in head, thorax and abdomen of the two wing-morph nymphs, and the expression level of Apirp4 in the thorax of the wingless nymphs increased significantly compared with that of the winged nymphs (Figure 2B). The expression of Apirp5 was detected in all tissues (Figure 2C). However, Apirp5 transcript levels were significantly higher in wingless vs. winged nymphs in both head and thorax (Figure 2C). The transcripts of Apirp6 were also detected in all tissues, and Apirp7 existed only in the thorax (Figures 2D,E). The expression levels of Apirp6 and Apirp7 also showed significant increase in the thorax of the wingless nymphs (Figures 2D,E). Apinrs, Apirs, and Apide2 transcripts were detected in all body parts (Figures 2F–H,J), but a higher level of Apide1 transcription was evident in midgut of winged nymphs (Figure 2I).
Figure 2. Expression levels of Apirps, Apinrs, Apirs, and Apides in different tissues of winged (black column) and wingless (gray column) third instar pea aphid. The abundance of Apirps (A–E), Apinrs (F,G), Apirs (H), Apides (I,J) transcripts were compared among different tissues and between winged and wingless third instar aphids. Numbers under key are P-values resulting from two-way ANOVAs analyzing the effects of body parts (P), wing morphs (M), and parts × wing morphs (P × M). If the P-value of P is < 0.05, it means variation of target gene expression levels was affected by wing morph type. If the P-value of M is < 0.05, it means that the expression level of target gene show significant difference between two wing morphs considering all body parts. If the P-value of P × M is < 0.05, it means that the body parts and two morphs significantly affected difference comparison between the two wing morphs in each body part. Means and standard errors are shown from four independent biological replicates. Relative gene expression data were normalized to Rpl7, and multiple comparisons were corrected using Sidak test: *p < 0.05; ***p < 0.001; ****p < 0.0001.
Body Weight, Ovaries or Embryos, and Nutrient Content of Wingless and Winged Nymphs, Apirp5 RNAi Wingless Nymphs and Controls
Differences in the phenotype of insects is associated with differentiation of organs or tissues and allocation of nutrient content. Third instar nymphs with wing primordia weighed less than wingless nymphs (Figure 3A). Embryos in the ovary of the wingless nymphs were significantly larger than those of the winged nymphs (Figure 3B and Figures S6A,B). Carbohydrate and protein contents of winged nymphs were significantly higher than those of wingless nymphs (Figures 4A,C), whereas glycogen content was similar between the wingless and winged nymphs (Figure 4B).
Figure 3. Comparison of body weights and embryo lengths among different groups of third instar pea aphid nymphs and expression levels of Apirp5 between wingless nymphs treated with dsgfp and dsApirp5. (A) Body weights of wingless nymphs is greater than that of winged nymphs (n = 70 wingless, 92 winged), as is embryo lengths (B) (n = 14 wingless, 14 winged). (C) Apirp5 expression is reduced in dsApirp5 treated wingless nymphs relative to that of dsgfp treated nymphs from four independent biological repeats. Relative gene expression data were normalized to Rpl7. (D) Body weight is reduced in dsApirp5 treated wingless nymphs (n = 32 dsgfp, 32 dsApirp5), as is embryo lengths (E) (n = 25 dsgfp, 29 dsApirp5), relative to that of dsgfp treated nymphs. Means and standard errors are shown. Data analyzed with Student's t-tests: *p < 0.05; **p < 0.01; ***p < 0.001.
Figure 4. Contents of carbohydrate (A,D), glycogen (B,E) and protein (C,F) in third instar nymphs. The comparison of carbohydrate (A) (n = 31 wingless, 28 winged), glycogen (B) (n = 32 wingless, 32 winged) and protein (C) (n = 32 wingless, 32 winged) contents between winged and wingless nymphs is shown in the top row. In the second row, the carbohydrate (D) (n = 12 dsgfp, 13 dsApirp5), glycogen (E) (n = 42 dsgfp, 28 dsApirp5) and proteins (F) (n = 29 dsgfp, 42 dsApirp5) are compared between dsgfp and dsApirp5 treatments in wingless nymphs. Means and standard errors are shown. Data analyzed with Student's t-tests: *p < 0.05; **p < 0.01.
According to previous reports, the large majority of Apirp5 ESTs in aphid libraries is likely involved in stimulating fast growth (Huybrechts et al., 2010). And in both of head and thorax, Apirp5 showed differential expression between third instar nymphs with and without wing primordia (Figure 2C), so the functions of Apirp5 on phenotypic plasticity of wing dimorphism were tested. Three days after dsApirp5 was injected into hemolymph of the second instar nymphs, expression levels of Apirp5 were reduced significantly compared with dsgfp controls (Figure 3C). Compared with dsgfp control, the dsApirp5 group showed lighter body weight and smaller embryo size (Figures 3D,E and Figures S6C,D). However, the formation of wing structures was not observed in dsApirp5-treated and control nymphs.
Wingless nymphs treated with dsApirp5 also had significantly higher carbohydrate and protein contents than the dsgfp control nymphs (Figures 4D,F), but glycogen content was similar between the two group nymphs (Figure 4E).
Discussion
The IS pathway directly or/and indirectly affects polyphenism in insects (Wolschin et al., 2011; Snell-Rood and Moczek, 2012; Wang et al., 2013; Xu et al., 2015). For wing dimorphic insects, nymphs with or without wing primoridia often show different development status and physiological performance. In this study, combined with previous reports, we presented the expression of selected genes in the IS pathway, and identified one strong candidate gene—Apirp5 (Figure 5). This report is the first study linking wing polyphenism with the IS pathway in aphids. Our results indicated that the IS pathway might have effects of regulation on embryo development and certain biochemical substance during wing differentiation in nymphs with wing polyphenism.
Figure 5. Summary diagram of embryo development and biochemical components between the winged and wingless third instar nymphs (A), and dsApirp5 and dsgfp treated third instar nymphs (B). The knockdown experiment showed that Apirp5 was involved in embryo development and regulation of carbohydrate and protein in wingless nymphs.
Of the 10 Apirps identified, Apirp2, 3, 8, 9, and 10 were not detected in the third instar nymphs. These undetected genes were probably expressed at low-levels in the third instar pea aphids or in special tissues or cells. The expression profiles of the other five Apirps varied dramatically among different body parts, indicating that the expression of Apirps might be associated with specific tissues. For instance, Apirp1 was only detected in aphid's head (Figure 2A), which contains median neurosecretory cell (MNC) that secretes many important ILPs in other insects, whereas Apirp7 transcript was only found in thorax (Figure 2E). However, Apirp5 and 6 transcripts were detected in whole body (Figures 2C,D), implying that Apirp5 and Apirp6 might serve universal roles in the third instar nymphs of A. pisum. In the brown planthopper (N. lugens), Nilp1 and Nilp3 are highly expressed in head, fat bodies, legs, wing buds, and cuticle in the fourth instar nymphs (Xu et al., 2015). Therefore, different irp or ilp genes could be expressed in the same or different tissues of these hemipteran insects. Similar findings are reported for Drosophila (Ikeya et al., 2002; Broughton et al., 2005; Veenstra and Sellami, 2008) and Aedes aegypti (Riehle et al., 2006). Although the MNC in the medial region of brain is the main source of irp or ilp expression in many insects, the fat body is another important source (Antonova et al., 2012). Unlike winged nymphs, there are no muscle tissues inside the thorax of wingless nymphs. Fat bodies filling in the thorax cavity may play a role in increasing expression levels of Apirp4, 5, 6, and 7. Because the IS pathway has been shown to affect juvenile hormone (JH) synthesis in other insects (Tatar et al., 2001; Tu et al., 2005), JH can inhibit wing development in aphids (Ishikawa et al., 2013). The higher expression levels of Apirps in thorax of wingless nymph than winged nymph might promote higher JH in wingless nymphs. It seemed that the decrease in expression levels of Apirps was associated with the winged phenotype of A. pisum.
Based on gene expression pattern, we found Apirp5 showed significant expression differences between two wing-morph nymphs in both head and thorax (Figure 2C), which are the core areas for secreting hormone that is associating with polyphenism in some insects (Hartfelder and Emlen, 2012). Besides Apirp5 and 7, expression levels of the other Apirps with significant differences between two morphs were low (Figures 2A,B,D,E). Because two pesudogenes are very similar to Apirp7 (Huybrechts et al., 2010), transcripts from pesudogenes probably affected quantification of transcripts of Apirp7. Besides, due to previous reports that Apirp5 with high expression abundance is likely a serious candidate for rapid growth in many aphids (Huybrechts et al., 2010), Apirp5 was chosen to knockdown by RNA interference.
High reproduction rate has previously been observed in wingless aphids (van Emden and Harrington, 2007). Our results showed that the embryos in the ovaries of wingless nymphs were larger than those of winged nymphs, and similar results are also found in other aphid species (Ishikawa and Miura, 2009). Meanwhile we found that the bodies of the third instar wingless nymphs were heavier than the third instar winged nymphs. This was consistent with previous findings that wingless aphids have a shorter development time than winged aphids (Dixon and Howard, 1986), and winged aphids are smaller than wingless aphids (Zera and Denno, 1997). The results of RNAi experiments suggested Apirp5 played a pivotal role in body weight and embryo development of A. pisum. The reduction of Apirp5 transcripts resulted in weight loss and smaller embryo development of wingless nymphs, which were features of winged nymphs. In A. mellifera, ovariole numbers connected with queen and worker development decrease when Amilp2 RNAi is conducted (Wang et al., 2013). In Drosophila, inhibiting expression of dilp2, 3, and 5 also causes development delay and serious weight loss (Rulifson et al., 2002). For another dipteran insect, A. aegypti, Aaegilp3 is shown to control egg production (Brown et al., 2008). Therefore, influences of Apirp5 on weights and embryo size of pea aphid indicated that this gene was likely involved in regulating embryo development during wing differentiation.
Further, we also found that the carbohydrate and protein contents were higher in winged nymphs than those in wingless nymphs. Such differences in nutritional metabolism between the two wing-morph nymphs could also be the result of IS pathway regulation and might be associated with their different morphological and physiological attributes. For instance, winged aphids possess indirect flight muscle and need fuel to provide energy for flight (Zera and Denno, 1997; Ishikawa et al., 2008). Therefore, high abundances of protein and trehalose, the main component of soluble carbohydrate in haemolymph (Wyatt and Kalf, 1957), are essential for flight muscle development and flight, respectively (Beenakkers et al., 1984; Marden, 2000). Our results showed that the aphid nymphs with reduced expression of Apirp5 presented an increase of soluble carbohydrate and protein contents of whole body. Compared with dsgfp controls, after dsApirp5 treatments, these changing trends of soluble carbohydrate and protein contents in wingless nymphs were similar to that in winged nymphs. This suggested that the changing of nutrient metabolism regulated by Apirp5 was likely to meet the demand of dispersal of A. pisum. Previous studies demonstrate that knockdown and knockout of dilp2 induce the increase of trehalose in whole body of Drosophila (Broughton et al., 2008; Grönke et al., 2010). This implied that the function of irp or ilp on regulation of soluble carbohydrates was likely to be conserved in insects. Since a swelling mesothorax was not observed after dsApirp5 injection, it was not clear whether an increase of protein content caused by Apirp5 down-regulation was connected with flight muscle construction. Future histological evidence in the mesothorax would explain whether increased proteins were used to construct flight muscles or not.
The dsRNA treated nymphs showed lower body weight, lower soluble carbohydrate, and higher glycogen than untreated nymphs (Figures 3A,D and Figures 4A,B,D,E). Because of the soft body of aphids, a stab wound made by the capillary resulted in the loss of hemolymph (Altincicek et al., 2008), as suggested by body weight loss. Therefore, it may be reasonable that content of soluble carbohydrate per body weight is decreased while that of glycogen increased, because the main soluble carbohydrate is trehalose in hemolymph and glycogen is contained in the non-hemolymph tissues including fat boy and muscle. The protein concentration, which was similar between injected and non-injected animals (Figures 4C,F), is probably similar in the hemolymph and other tissues.
Although dsApirp5 treatments affected body weights, embryo sizes and body composition, the decrease of Apirp5 expression did not promote the development of wing apparatus in the nymphs. There were some possible explanations. First, the expression of Apirp5 was associated with the development of embryos in the ovaries of aphid nymphs during wing differentiation, but Apirp5 probably did not affect the formation of wing apparatus. The second explanation could be the time point of microinjection. The process of wing differentiation starts in the second instar (Ogawa and Miura, 2013), but it was difficult to inject dsRNA into the first instar nymphs (up to 100% mortality within 48 h in our preliminary analysis). The effects of redundancy among irps on lifespan are found in Drosophila (Broughton et al., 2008) and thus may be another reason that only changing Apirp5 expression levels could not affect wing morph of aphid nymphs, because Apirp1 in head and Apirp4, 6, and 7 in thorax showed significant differences between the two wing morphs. Lastly, low RNAi efficiency has been observed in aphids (Jaubert-Possamai et al., 2007; Christiaens et al., 2014). Our RNAi test also showed the similar result, although we used high concentrations of dsRNA for treatments. Future research will focus on analyzing the compensation and redundancy of Apirps using more efficient gene knockdown methods. Based on existing results, we could not rule out a role for Apirp5 in wing development.
Other genes in the IS pathway are involved in the regulation of traits associated with trade-offs in other insects. There is only one insulin receptor in Diptera (Fernandez et al., 1995; Ruan et al., 1995), but two insulin receptor homologs are identified in Coleoptera, Hymenoptera, and Hemiptera (Antonova et al., 2012). That raised a question of whether the two insulin receptors had different functions or redundant functions. Although two insulin receptors have been proven to regulate wing polyphenism in the brown planthopper (Xu et al., 2015), in our study, Apinr1 only showed a significant difference in expression in the midgut between wingless and winged nymphs; therefore, it was hard to establish a link between Apinrs and wing polyphenism. Amirs RNAi tests on A. mellifer have shown a strong link between the IRS and caste fate (Mutti et al., 2011; Wolschin et al., 2011). We found that there were no significant differences in expression levels of Apirs in any of the tissues between winged and wingless nymphs in pea aphid. This suggested that Apirs might not be the main factor in the IS pathway for regulation of wing differentiation. For Apides, we found that expression of Apide1 was higher in winged nymphs than in wingless. There is no evidence that IDEs degrade insect ILPs (Antonova et al., 2012), but knockout of dide increases the content of trehalose in Drosophila (Tsuda et al., 2010). The connection between sugar metabolism in pea aphid and Apides is worth pursuing in future study.
In this study, we tried to investigate the relationship between the IS pathway and wing polyphenism of pea aphid in the nymphal phase by comparing gene expression levels in the IS pathway and RNAi for one Apirp. Our results indicated that Apirp5 played an important role in regulation of metabolism and development of embryos during wing differentiation in pea aphid. These results gave further recognition on the roles of IS pathway on phenotypic plasticity in aphids.
Author Contributions
SG designed and conducted this study, collected data, analyzed the data, and wrote the paper; MZ helped in all phrases of this study; Prof. TL designed and supervised the whole study, and wrote the paper.
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.
Acknowledgments
We are grateful to Mark Brown and Patricia Pietrantonio for constructive and valuable comments and editing on the earlier versions of this manuscript. We also thank Ron Salzman for revising manuscript. This work was supported by the National Key Basic Research Program of China (973 Program) (No. 2013CB127600 and No. 2012CB114105).
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fphys.2016.00031
Abbreviations
ILP, insulin-like peptide; IRP, insulin-related peptide; InR, insulin receptor; IRS, insulin receptor substrate; IDE, insulin degrading enzyme; HD, high density; LD, low density; IS, insulin signaling; ilps, insulin-like peptide genes; irps, insulin-related peptide genes; inrs, insulin receptor genes; irs, insulin receptor substrate gene; ides, insulin degrading enzyme genes.
References
Abràmoff, M. D., Magalhães, P. J., and Ram, S. J. (2004). Image processing with Image. J. Biophotonics Int. 11, 36–43. Available online at: http://dspace.library.uu.nl/handle/1874/204900
Altincicek, B., Gross, J., and Vilcinskas, A. (2008). Wounding-mediated gene expression and accelerated viviparous reproduction of the pea aphid Acyrthosiphon pisum. Insect Mol. Biol. 17, 711–716. doi: 10.1111/j.1365-2583.2008.00835.x
Antonova, Y., Arik, A. J., Moore, W., Riehle, M., Brown, M. R., and Gilbert, L. (2012). “Insulin-like peptides: structure, signaling, and function,” in Insect Endocrinology, ed L. I. Gilbert (Amsterdam: Academic Press; Elsevier), 63–92.
Beenakkers, A. T., Van der Horst, D., and Van Marrewijk, W. (1984). Insect flight muscle metabolism. Insect Biochem. 14, 243–260. doi: 10.1016/0020-1790(84)90057-X
Bonte, D., Van Dyck, H., Bullock, J. M., Coulon, A., Delgado, M., Gibbs, M., et al. (2012). Costs of dispersal. Biol. Rev. 87, 290–312. doi: 10.1111/j.1469-185X.2011.00201.x
Braendle, C., Davis, G. K., Brisson, J. A., and Stern, D. L. (2006). Wing dimorphism in aphids. Heredity 97, 192–199. doi: 10.1038/sj.hdy.6800863
Brisson, J. A. (2010). Aphid wing dimorphisms: linking environmental and genetic control of trait variation. Phil. Trans. R. Soc. B Biol. Sci. 365, 605–616. doi: 10.1098/rstb.2009.0255
Brisson, J. A., and Stern, D. L. (2006). The pea aphid, Acyrthosiphon pisum: an emerging genomic model system for ecological, developmental and evolutionary studies. Bioessays 28, 747–755. doi: 10.1002/bies.20436
Broughton, S., Alic, N., Slack, C., Bass, T., Ikeya, T., Vinti, G., et al. (2008). Reduction of DILP2 in Drosophila triages a metabolic phenotype from lifespan revealing redundancy and compensation among DILPs. PLoS ONE 3:e3721. doi: 10.1371/journal.pone.0003721
Broughton, S. J., Piper, M. D., Ikeya, T., Bass, T. M., Jacobson, J., Driege, Y., et al. (2005). Longer lifespan, altered metabolism, and stress resistance in Drosophila from ablation of cells making insulin-like ligands. Proc. Natl. Acad. Sci. U.S.A. 102, 3105–3110. doi: 10.1073/pnas.0405775102
Brown, M. R., Clark, K. D., Gulia, M., Zhao, Z., Garczynski, S. F., Crim, J. W., et al. (2008). An insulin-like peptide regulates egg maturation and metabolism in the mosquito Aedes aegypti. Proc. Natl. Acad. Sci. 105, 5716–5721. doi: 10.1073/pnas.0800478105
Christiaens, O., Swevers, L., and Smagghe, G. (2014). DsRNA degradation in the pea aphid (Acyrthosiphon pisum) associated with lack of response in RNAi feeding and injection assay. Peptides 53, 307–314. doi: 10.1016/j.peptides.2013.12.014
Consortium, I. A. G. (2010). Genome sequence of the pea aphid Acyrthosiphon pisum. PLoS Biol. 8:e1000313. doi: 10.1371/journal.pbio.1000313
de Azevedo, S. V., and Hartfelder, K. (2008). The insulin signaling pathway in honey bee (Apis mellifera) caste development—differential expression of insulin-like peptides and insulin receptors in queen and worker larvae. J. Insect Physiol. 54, 1064–1071. doi: 10.1016/j.jinsphys.2008.04.009
Denno, R. F., Olmstead, K. L., and McCloud, E. S. (1989). Reproductive cost of flight capability: a comparison of life history traits in wing dimorphic planthoppers. Ecol. Entomol. 14, 31–44. doi: 10.1111/j.1365-2311.1989.tb00751.x
Dixon, A., Horth, S., and Kindlmann, P. (1993). Migration in insects: cost and strategies. J. Anim. Ecol. 62, 182–190. doi: 10.2307/5492
Dixon, A., and Howard, M. (1986). “Dispersal in aphids, a problem in resource allocation,” in Insect Flight, ed W. Danthanarayana (Heidelberg: Springer Press), 145–151.
Emlen, D., Szafran, Q., Corley, L., and Dworkin, I. (2006). Insulin signaling and limb-patterning: candidate pathways for the origin and evolutionary diversification of beetle ‘horns’. Heredity 97, 179–191. doi: 10.1038/sj.hdy.6800868
Fernandez, R., Tabarini, D., Azpiazu, N., Frasch, M., and Schlessinger, J. (1995). The Drosophila insulin receptor homolog: a gene essential for embryonic development encodes two receptor isoforms with different signaling potential. EMBO J. 14, 3373.
Foray, V., Pelisson, P.-F., Bel-Venner, M.-C., Desouhant, E., Venner, S., Giron, D., et al. (2012). A handbook for uncovering the complete energetic budget in insects: the van Handel's method (1985) revisited. Physiol. Entomol. 37, 295–302. doi: 10.1111/j.1365-3032.2012.00831.x
Garofalo, R. S. (2002). Genetic analysis of insulin signaling in Drosophila. Trends Endocrinol. Metab. 13, 156–162. doi: 10.1016/S1043-2760(01)00548-3
Grönke, S., Clarke, D.-F., Broughton, S., Andrews, T. D., and Partridge, L. (2010). Molecular evolution and functional characterization of Drosophila insulin-like peptides. PLoS Genet. 6:e1000857. doi: 10.1371/journal.pgen.1000857
Hartfelder, K., and Emlen, D. (2012). “Endocrine control of insect polyphenism,” in Insect Endocrinology, ed L. I. Gilbert (Amsterdam: Academic Press (Elsevier)), 484–522.
Huybrechts, J., Bonhomme, J., Minoli, S., Prunier-Leterme, N., Dombrovsky, A., Abdel-Latief, M., et al. (2010). Neuropeptide and neurohormone precursors in the pea aphid, Acyrthosiphon pisum. Insect Mol. Biol. 19, 87–95. doi: 10.1111/j.1365-2583.2009.00951.x
Ikeya, T., Galic, M., Belawat, P., Nairz, K., and Hafen, E. (2002). Nutrient-dependent expression of insulin-like peptides from neuroendocrine cells in the CNS contributes to growth regulation in Drosophila. Curr. Biol. 12, 1293–1300. doi: 10.1016/S0960-9822(02)01043-6
Ishikawa, A., Gotoh, H., Abe, T., and Miura, T. (2013). Juvenile hormone titer and wing-morph differentiation in the vetch aphid Megoura crassicauda. J. Insect Physiol. 59, 444–449. doi: 10.1016/j.jinsphys.2013.02.004
Ishikawa, A., Hongo, S., and Miura, T. (2008). Morphological and histological examination of polyphenic wing formation in the pea aphid Acyrthosiphon pisum (Hemiptera, Hexapoda). Zoomorphology 127, 121–133. doi: 10.1007/s00435-008-0057-5
Ishikawa, A., and Miura, T. (2009). Differential regulations of wing and ovarian development and heterochronic changes of embryogenesis between morphs in wing polyphenism of the vetch aphid. Evol. Dev. 11, 680–688. doi: 10.1111/j.1525-142X.2009.00375.x
Jaubert-Possamai, S., Le Trionnaire, G., Bonhomme, J., Christophides, G. K., Rispe, C., and Tagu, D. (2007). Gene knockdown by RNAi in the pea aphid Acyrthosiphon pisum. BMC Biotechnol. 7:63. doi: 10.1186/1472-6750-7-63
Kruger, N. J. (2009). “The Bradford method for protein quantitation,” in The Protein Protocols Handbook, ed J. M. Walker (Heidelberg: Springer Press), 17–24.
MacKay, P. A., and Wellington, W. (1976). Maternal age as a source of variation in the ability of an aphid to produce dispersing forms. Res. Popul. Ecol. 18, 195–209. doi: 10.1007/BF02754093
Marden, J. H. (2000). Variability in the size, composition, and function of insect flight muscles. Annu. Rev. Physiol. 62, 157–178. doi: 10.1146/annurev.physiol.62.1.157
Mole, S., and Zera, A. J. (1993). Differential allocation of resources underlies the dispersal-reproduction trade-off in the wing-dimorphic cricket, Gryllus rubens. Oecologia 93, 121–127. doi: 10.1007/BF00321201
Mutti, N. S., Dolezal, A. G., Wolschin, F., Mutti, J. S., Gill, K. S., and Amdam, G. V. (2011). IRS and TOR nutrient-signaling pathways act via juvenile hormone to influence honey bee caste fate. J. Exp. Biol. 214, 3977–3984. doi: 10.1242/jeb.061499
Nakabachi, A., Shigenobu, S., Sakazume, N., Shiraki, T., Hayashizaki, Y., Carninci, P., et al. (2005). Transcriptome analysis of the aphid bacteriocyte, the symbiotic host cell that harbors an endocellular mutualistic bacterium, Buchnera. Proc. Natl. Acad. Sci. U.S.A. 102, 5477–5482. doi: 10.1073/pnas.0409034102
Ogawa, K., and Miura, T. (2013). Two developmental switch points for the wing polymorphisms in the pea aphid Acyrthosiphon pisum. Evodevo 4, 30. doi: 10.1186/2041-9139-4-30
Pfaffl, M. W. (2001). A new mathematical model for relative quantification in real-time RT–PCR. Nucleic Acids Res. 29, e45. doi: 10.1093/nar/29.9.e45
Ramakers, C., Ruijter, J. M., Deprez, R. H. L., and Moorman, A. F. (2003). Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci. Lett. 339, 62–66. doi: 10.1016/S0304-3940(02)01423-4
Riehle, M. A., Fan, Y., Cao, C., and Brown, M. R. (2006). Molecular characterization of insulin-like peptides in the yellow fever mosquito, Aedes aegypti: expression, cellular localization, and phylogeny. Peptides 27, 2547–2560. doi: 10.1016/j.peptides.2006.07.016
Ruan, Y., Chen, C., Cao, Y., and Garofalo, R. S. (1995). The Drosophila insulin receptor contains a novel carboxyl-terminal extension likely to play an important role in signal transduction. J. Biol. Chem. 270, 4236–4243. doi: 10.1074/jbc.270.9.4236
Rulifson, E. J., Kim, S. K., and Nusse, R. (2002). Ablation of insulin-producing neurons in flies: growth and diabetic phenotypes. Science 296, 1118–1120. doi: 10.1126/science.1070058
Shen, Y., Joachimiak, A., Rosner, M. R., and Tang, W.-J. (2006). Structures of human insulin-degrading enzyme reveal a new substrate recognition mechanism. Nature 443, 870–874. doi: 10.1038/nature05143
Simpson, S. J., Sword, G. A., and Lo, N. (2011). Polyphenism in insects. Curr. Biol. 21, R738–R749. doi: 10.1016/j.cub.2011.06.006
Snell-Rood, E. C., and Moczek, A. P. (2012). Insulin signaling as a mechanism underlying developmental plasticity: the role of FOXO in a nutritional polyphenism. PLoS ONE 7:e34857. doi: 10.1371/journal.pone.0034857
Sutherland, O. (1969). The role of crowding in the production of winged forms by two strains of the pea aphid, Acyrthosiphon pisum. J. Insect Physiol. 15, 1385–1410. doi: 10.1016/0022-1910(69)90199-1
Tamura, K., Stecher, G., Peterson, D., Filipski, A., and Kumar, S. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 30, 2725–2729. doi: 10.1093/molbev/mst197
Taniguchi, C. M., Emanuelli, B., and Kahn, C. R. (2006). Critical nodes in signalling pathways: insights into insulin action. Nat. Rev. Mol. Cell Biol. 7, 85–96. doi: 10.1038/nrm1837
Tatar, M., Kopelman, A., Epstein, D., Tu, M.-P., Yin, C.-M., and Garofalo, R. (2001). A mutant Drosophila insulin receptor homolog that extends life-span and impairs neuroendocrine function. Science 292, 107–110. doi: 10.1126/science.1057987
Tsuda, M., Kobayashi, T., Matsuo, T., and Aigaki, T. (2010). Insulin-degrading enzyme antagonizes insulin-dependent tissue growth and Aβ-induced neurotoxicity in Drosophila. FEBS Lett. 584, 2916–2920. doi: 10.1016/j.febslet.2010.05.010
Tu, M.-P., Yin, C.-M., and Tatar, M. (2005). Mutations in insulin signaling pathway alter juvenile hormone synthesis in Drosophila melanogaster. Gen. Comp. Endocrinol. 142, 347–356. doi: 10.1016/j.ygcen.2005.02.009
Van Handel, E. (1965). Microseparation of glycogen, sugars, and lipids. Anal. Biochem. 11, 266–271. doi: 10.1016/0003-2697(65)90014-X
Veenstra, J. A., and Sellami, A. (2008). Regulatory peptides in fruit fly midgut. Cell Tissue Res. 334, 499–516. doi: 10.1007/s00441-008-0708-3
Wang, Y., Azevedo, S. V., Hartfelder, K., and Amdam, G. V. (2013). Insulin-like peptides (AmILP1 and AmILP2) differentially affect female caste development in the honey bee (Apis mellifera L.). J. Exp. Biol. 216, 4347–4357. doi: 10.1242/jeb.085779
Wheeler, D., Buck, N., and Evans, J. D. (2006). Expression of insulin pathway genes during the period of caste determination in the honey bee, Apis mellifera. Insect Mol. Biol. 15, 597–602. doi: 10.1111/j.1365-2583.2006.00681.x
Wolschin, F., Mutti, N. S., and Amdam, G. V. (2011). Insulin receptor substrate influences female caste development in honeybees. Biol. Lett. 7, 112–115. doi: 10.1098/rsbl.2010.0463
Wyatt, G., and Kalf, G. (1957). The chemistry of insect hemolymph II. Trehalose and other carbohydrates. J. Gen. Physiol. 40, 833–847. doi: 10.1085/jgp.40.6.833
Xu, H.-J., Xue, J., Lu, B., Zhang, X.-C., Zhuo, J.-C., He, S.-F., et al. (2015). Two insulin receptors determine alternative wing morphs in planthoppers. Nature 519, 464–467. doi: 10.1038/nature14286
Keywords: aphid, insulin-like peptide, wing polyphenism, wing differentiation, biochemical components
Citation: Guo S-S, Zhang M and Liu T-X (2016) Insulin-Related Peptide 5 is Involved in Regulating Embryo Development and Biochemical Composition in Pea Aphid with Wing Polyphenism. Front. Physiol. 7:31. doi: 10.3389/fphys.2016.00031
Received: 14 October 2015; Accepted: 22 January 2016;
Published: 09 February 2016.
Edited by:
Senthil-Nathan Sengottayan, Manonmaniam Sundaranar University, IndiaReviewed by:
Akira Mizoguchi, Nagoya University, JapanKrzysztof Jagla, Institut National de la Santé et de la Recherche Médicale, France
Copyright © 2016 Guo, Zhang and Liu. 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) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tong-Xian Liu, dHhsaXVAbndzdWFmLmVkdS5jbg==
†These authors have contributed equally to this work.