A novel gene REPTOR2 activates the autophagic degradation of wing disc in pea aphid
Abstract
Wing dimorphism in insects is an evolutionarily adaptive trait to maximize insect fitness under various environments, by which the population could be balanced between dispersing and reproduction. Most studies concern the regulatory mechanisms underlying the stimulation of wing morph in aphids, but relatively little research addresses the molecular basis of wing loss. Here, we found that, while developing normally in winged-destined pea aphids, the wing disc in wingless-destined aphids degenerated 30-hr postbirth and that this degeneration was due to autophagy rather than apoptosis. Activation of autophagy in first instar nymphs reduced the proportion of winged aphids, and suppression of autophagy increased the proportion. REPTOR2, associated with TOR signaling pathway, was identified by RNA-seq as a differentially expressed gene between the two morphs with higher expression in the thorax of wingless-destined aphids. Further genetic analysis indicated that REPTOR2 could be a novel gene derived from a gene duplication event that occurred exclusively in pea aphids on autosome A1 but translocated to the sex chromosome. Knockdown of REPTOR2 reduced autophagy in the wing disc and increased the proportion of winged aphids. In agreement with REPTOR’s canonical negative regulatory role of TOR on autophagy, winged-destined aphids had higher TOR expression in the wing disc. Suppression of TOR activated autophagy of the wing disc and decreased the proportion of winged aphids, and vice versa. Co-suppression of TOR and REPTOR2 showed that dsREPTOR2 could mask the positive effect of dsTOR on autophagy, suggesting that REPTOR2 acted as a key regulator downstream of TOR in the signaling pathway. These results revealed that the TOR signaling pathway suppressed autophagic degradation of the wing disc in pea aphids by negatively regulating the expression of REPTOR2.
Editor's evaluation
This study addresses the genetic basis of an iconic example of developmental plasticity. The background and hypothesis are clear, and the experiments are well conducted. The paper shows that winglessness in aphids involves autophagy (rather than apoptosis) of the wing primordium, and that this is regulated by a novel gene, REPTOR2, which is a target of TOR signaling.
https://doi.org/10.7554/eLife.83023.sa0Introduction
Flight ability is a great evolutionary success for animals that favors seeking new habitats, mates, and food resources (Roff, 1990). The loss of such ability has occurred in numerous insects (Roff, 1990; Harrison, 1980; Roff and Fairbairn, 1991). The wingless morph has evolved to optimize resource allocation from dispersal to reproduction, which facilitates population expansion under suitable conditions (Zera and Denno, 1997; Mole and Zera, 1993). Alternation between the wing morphs in aphid offspring is strongly influenced by maternal experience (Braendle et al., 2006). Pea aphids are recognized as a compelling laboratory model to study regulatory mechanisms underlying wing dimorphism, most notably through maternal population density (herein refers to as maternal density). It has been shown that ecdysteroids, miR-9b-ABCG4-insulin signaling, and two laterally transferred viral genes (Apns-1 and Apns-2) are responsible for producing winged offspring in parthenogenetic females (Vellichirammal et al., 2017; Shang et al., 2020; Parker and Brisson, 2019). Nevertheless, the phylogenetic evidence suggests that winged morph in aphids is a default developmental pattern, while the wingless morph is secondarily derived, from their winged ancestor (Braendle et al., 2006; Miyazaki, 1987). Thus, how wingless morph is controlled under favorable conditions needs to be studied, as such information will help to elucidate the differentiation of wing dimorphism in aphids.
Tissue degradation is a common developmental process in insects executed by programmed cell death (PCD), such as autophagy and apoptosis (Xu et al., 2019). During metamorphosis, Drosophila salivary glands, midgut and other larval tissues are removed by PCD (Xu et al., 2019). For example, suppression of apoptosis or autophagy or both could delay the degeneration of salivary glands (Berry and Baehrecke, 2007). Suppressing autophagy, but not apoptosis, severely interferes with the removal of the midgut during metamorphosis (Denton et al., 2012). These findings suggest that tissue-specific degeneration required different types of PCD. Furthermore, histological examinations in pea aphids show that the second instar nymphs of the winged-destined morph obviously possess a pair of wing discs, whereas wingless-destined morph does not (Ishikawa et al., 2008; Ogawa and Miura, 2013). Since the wing bud and the full wing are developed from the wing disc, it would be interesting to know whether and how PCD controls the fate of the wing disc cells, ultimately resulting in the wingless morph.
The evolutionarily conserved nutrient-sensing TOR (the target of rapamycin) pathway has been implicated in regulation of fecundity, life span as well as tissues degradation of insects under different nutritional conditions (Xu et al., 2019; Katewa and Kapahi, 2011). TOR coordinately regulates the balance between growth and autophagy in response to cellular and physiological conditions as well as environmental stresses such as intracellular amino acid availability and accessible energy levels (Saxton and Sabatini, 2017). Shortage of dietary amino acid depletion or starvation, inhibits TOR activity, which is known as a crucial step for autophagy induction in eukaryotes (Shimobayashi and Hall, 2016). Downregulation of TOR signaling pathway seems to be necessary for autophagic degradation to occur in different insect tissues such as the salivary gland and midgut (Berry and Baehrecke, 2007; Denton et al., 2012). TOR suppresses autophagy by directly phosphorylating multiple residues of autophagy-related proteins or repressing the transcription of autophagy-related genes. TOR phosphorylates Atg13, causing reduced affinity between Atg1 and its binding proteins and subsequently repressed autophagy (Puente et al., 2016). Alternatively, TOR also phosphorylates Ser527 and Ser530 of REPTOR (repressed by TOR) to prevent its entry into nucleus where REPTOR targets autophagy-related genes and suppresses their expression (Tiebe et al., 2015). In this study, we identify a novel gene REPTOR2 in pea aphids by comparing the transcriptome of wingless-destined versus winged-destined aphids at 24-hr postbirth. Since REPTOR2 has a higher expression in thorax of wingless-destined aphids, we hypothesized that REPTOR2 activates autophagic degradation in the wing disc to produce the wingless morph. Likewise, winged-destined aphids have a stronger TOR activity to suppress autophagy in the wing disc to develop into full wings. Experimentally, we specifically determined: (i) the time point to initiate PCD in the wing disc, (ii) the PCD type that was responsible for wing disc degeneration, and (iii) the role of REPTOR2 and TOR signaling cascade in determination wing dimorphism.
Results
Effects of maternal density and duration on the proportion of winged offspring
To efficiently sample winged offspring, we examined how maternal density and duration affected the proportion of winged offspring. Compared to the adults individually housed, the adults housed in groups of 2, 4, and 8 in a petri dish for 8-hr led to higher proportions of winged offspring, respectively (Figure 1A). Since two adults were sufficient to produce over 80% winged offspring, we then tested different contact time periods on the proportion of winged morph (Figure 1B). Apparently, 4-hr was sufficient to induce >90% winged offspring, and 2-hr induced a moderate effect (approximately 50%) (Figure 1B). Furthermore, two adults with 4-hr contact duration continuously produced over 80% winged offspring for 3 days, followed by a small reduction afterward (Figure 1C). Furthermore, to eliminate a strain-specific induction, three other pea aphid strains were also tested for two adults and 4-hr duration; even the strain with the lowest proportion of winged offspring exceeded 40% (Figure 1—figure supplement 1).
-
Figure 1—source data 1
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig1-data1-v1.xlsx
-
Figure 1—source data 2
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig1-data2-v1.xlsx
-
Figure 1—source data 3
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig1-data3-v1.xlsx
-
Figure 1—source data 4
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig1-data4-v1.xlsx
Developmental differentiation in the wing and wing disc of winged versus wingless morphs
Previous studies indicated that the wing disc was degenerated at the second instar in wingless morph of the pea aphid (Ishikawa et al., 2008; Ogawa and Miura, 2013). We compared the wing morphology across different nymphal stages between the two wing morphs. Winged- and wingless-destined aphids could not be distinguished morphogenetically until the third instar (Figure 1D). The wing bud of winged-destined aphids became visible during the third and fourth instars, which later developed into full unfolded wing once adult emergence. By contrast, wingless-destined aphids did not possess obvious wing buds at the third and fourth instar stages. To further determine the time point of initiation of the wing disc degeneration, the first instar nymphs at 24-, 30-, and 36-hr postbirth were collected for histological analyses. The wing disc sizes were not significantly different between the two morphs 24-hr postbirth, but were then rapidly degenerated in wingless-destined morphs while continuously developed and expended in winged-destined morphs 30-hr postbirth (Figure 1E and F).
Autophagic rather than apoptotic degradation was responsible for degeneration of wing disc
Since wing disc degeneration initiated at 30-hr post birth, PCD occurred in wing disc at this time point was characterized for its type and process by transmission electron microscopy (TEM), immunofluorescence, TdT-mediated dUTP Nick-End Labeling (TUNEL), and quantitative polymerase chain reaction (qPCR). Typical autophagic features such as numerous cytoplasmic vacuoles, large autophagosomes with sequestered visible remnants of organelles, and membranous whorls were observed by TEM within the cells of the wing disc of the wingless-destined morph. By contrast, only mitophagy and a few cytoplasmic vacuoles were found in wing disc cells of winged-destined morph. Signs of autophagy were observed in the wing disc of 54% cells of the wingless-destined morph, but only 15% cells of winged-destined aphids (Figure 2A). Furthermore, cell shrinkage, plasma membrane blebbing, apoptotic body, chromatin condensation, nucleolus disorganization, and nuclear fragmentation, typical features of apoptosis, were hardly observed in wing discs of either winged- or wingless-destined morphs.
-
Figure 2—source data 1
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig2-data1-v1.xlsx
-
Figure 2—source data 2
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig2-data2-v1.xlsx
-
Figure 2—source data 3
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig2-data3-v1.xlsx
-
Figure 2—source data 4
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig2-data4-v1.xlsx
-
Figure 2—source data 5
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig2-data5-v1.xlsx
-
Figure 2—source data 6
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig2-data6-v1.xlsx
-
Figure 2—source data 7
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig2-data7-v1.xlsx
-
Figure 2—source data 8
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig2-data8-v1.xlsx
-
Figure 2—source data 9
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig2-data9-v1.xlsx
The immunofluorescence and TUNEL results showed that autophagy, instead of apoptosis, was activated in wing disc cells of the wingless-destined morph (Figure 2B and C). The autophagy-related genes ATG1a, ATG1b, ATG3a, ATG3b, ATG4, ATG10, ATG13, and ATG14 had higher expression in the wingless-destined morph (Figure 2D), while all apoptosis-related genes except for Decay2 were not significantly different between two wing morphs (Figure 2E). These results indicated that autophagy rather than apoptosis was responsible for degeneration of the wing disc in pea aphids, and thereby reduced the proportion of the winged morph. To verify this hypothesis, the group with two adults placed in a petri dish for 2-hr was applied to produce a moderate proportion (approximately 50%) of winged offspring (Figure 2—figure supplement 1). Newborn nymphs were fed with autophagy agonist rapamycin or antagonist 3-MA for 24-hr. As anticipated, rapamycin reduced winged proportion from 67% to 36% by enhancing autophagy in the wing disc of first instar nymphs (Figure 2F and H), whereas 3-MA attenuated autophagy in the wing disc and increased the winged aphids from 54% to 67% (Figure 2G and I).
A novel gene REPTOR2 in pea aphid positively regulated the wing disc autophagy
Of the 387 differentially expressed genes (DEGs), in total, the first instar nymphs of the wingless-destined morph had 220 genes expressed higher and 167 genes expressed lower than the winged-destined morph at 24-hr postbirth (Figure 3A). Strong evidence has suggested that autophagic cell death is regulated by TORC1 (Berry and Baehrecke, 2007; Denton et al., 2012), we therefore focused on genes in the TOR signaling pathway (Figure 3B) that had significantly higher expression in wingless-destined morph than the winged counterpart. Only REPTOR2 (repress by TOR 2, LOC103309809) and SLC38A9 (sodium-coupled neutral amino acid transporter 9, LOC100158916) fell into this category. SLC38A9 is a lysosomal membrane-resident protein involved in the amino acid transport and responsible for amino acids-induced TORC1 activation (Rebsamen et al., 2015). REPTOR has been recognized as a transcription factor that relays the TOR signaling cascade by activating the expression of autophagy-related genes (Tiebe et al., 2015).
-
Figure 3—source data 1
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig3-data1-v1.xlsx
-
Figure 3—source data 2
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig3-data2-v1.xlsx
-
Figure 3—source data 3
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig3-data3-v1.xlsx
-
Figure 3—source data 4
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig3-data4-v1.xlsx
-
Figure 3—source data 5
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig3-data5-v1.xlsx
-
Figure 3—source data 6
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig3-data6-v1.xlsx
-
Figure 3—source data 7
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig3-data7-v1.xlsx
-
Figure 3—source data 8
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig3-data8-v1.xlsx
-
Figure 3—source data 9
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig3-data9-v1.xlsx
When aligned with the Acyrthosiphon pisum genome, two REPTOR homologs, ApREPTOR1 (LOC100168197) and ApREPTOR2 (LOC103309809), were identified. Seven other aphid species contained only a single REPTOR copy, including Myzus cerasi, Myzus persicae, Melanaphis secchari, Aphis gossyppii, Rhopalosiphum padi, Rhopalosiphum maidis, and Eriosoma lanigerum. Both REPTOR genes in A. pisum consisted of a conserved basic region leucine zipper (BRLZ) domain (Figure 3C). An unrooted phylogenetic tree with bootstrap consensus was constructed across eight aphid species using coding sequence alignments (Figure 3D). MUMmer analyses were used to determine the similarity and homologous regions of ApREPTOR2 in four aphid species with chromosome-level genomes in Aphidbase. ApREPTOR1 and MpREPTOR had 9 exons, while ApREPTOR2 had 7 exons (Figure 3E). ApREPTOR1 shared high sequence similarity with MpREPTOR, RmREPTOR, and ApREPTOR2, but not with ElREPTOR. The variations in the flanking region of REPTOR indicated that ApREPTOR2 was quite different from ApREPTOR1 (Figure 3F). Chromosomal distribution results further showed that REPTOR1 was located in A1 chromosome while REPTOR2 was located in the X chromosome of A. pisum (Figure 3G). These results indicated that REPTOR1 in pea aphid may have experienced a gene duplication to generate a novel gene REPTOR2 which probably translocates from autosome A1 to the sex chromosome.
REPTOR2 had higher expression in the wingless-destined morph at 12- and 24-hr postbirth than that of the winged-destined morph (Figure 3H), and the highest expression occurred in the thorax (Figure 3I), suggesting that REPTOR2 was essential in activating autophagic degradation of the wing disc. To further support this notion, we fed newborn nymphs with dsREPTOR2-RNA and observed reduced REPTOR2 transcripts in the wing disc and whole body by 53% and 65%, respectively (Figure 3K and L), as well as decreased gene expressions of ATG13 and ATG14 (Figure 3J). Knockdown of REPTOR2 attenuated autophagy in the wing disc (Figure 3M), resulting in higher proportion of the winged morph, from 39% to 51% (Figure 3N).
The effect of silencing TOR on wing disc autophagy and the proportion of the winged morph
Suppression of TOR is known to induce autophagic cell death during tissue degeneration in insects (Berry and Baehrecke, 2007; Denton et al., 2012). As expected in the wing disc, the fluorescence in situ hybridization (FISH) data indicated that TOR had higher abundance in winged- than wingless-destined morphs (Figure 4A). Feeding newborn nymphs with dsTOR-RNA reduced the TOR transcripts in the whole body and the wing disc by 48% and 34%, respectively (Figure 4B and C), and also declined winged proportion from 69% to 52% (Figure 4D). By contrast, activation of TOR via pharmacological agonist MHY1485 suppressed autophagy in the wing disc (Figure 4E) and increased winged proportion from 37% to 55% (Figure 4F). Since availability of amino acids is an important factor regulating TORC1 activity, depletion of dietary amino acids was used in the current study to suppress the TOR signaling pathway in pea aphids. Suppression of TOR via depletion of amino acids by 50% enhanced autophagy in wing disc (Figure 4G), and decreased winged proportion from 49% to 15% (Figure 4H).
-
Figure 4—source data 1
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig4-data1-v1.xlsx
-
Figure 4—source data 2
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig4-data2-v1.xlsx
-
Figure 4—source data 3
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig4-data3-v1.xlsx
-
Figure 4—source data 4
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig4-data4-v1.xlsx
-
Figure 4—source data 5
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig4-data5-v1.xlsx
-
Figure 4—source data 6
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig4-data6-v1.xlsx
-
Figure 4—source data 7
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig4-data7-v1.xlsx
-
Figure 4—source data 8
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig4-data8-v1.xlsx
-
Figure 4—source data 9
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig4-data9-v1.xlsx
-
Figure 4—source data 10
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig4-data10-v1.xlsx
-
Figure 4—source data 11
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig4-data11-v1.xlsx
-
Figure 4—source data 12
- https://cdn.elifesciences.org/articles/83023/elife-83023-fig4-data12-v1.xlsx
Effects of co-suppression of TOR and REPTOR2 on wing disc autophagy and proportion of the winged morph
To establish the causal link between REPTOR2 and TOR, both genes were silenced via dsRNA feeding. Interestingly, aphids displayed lower wing disc autophagy and had higher winged proportion than dsTOR aphids (Figure 4I and J). Furthermore, dsTOR aphids showed increased REPTOR2 expression in the wing disc and the whole body by 1.97- and 1.75-fold, respectively (Figure 4K and L), suggesting that REPTOR2 was downstream of TOR in the signaling pathway.
Discussion
The wing disc consists of undifferentiated and proliferating cells designated to develop into a full insect wing during the nymphal stage (Neto-Silva et al., 2009). Fate of the wing disc is presumably determined by maternal signals triggered by environmental signals such as population density, physical contact, and the presence of natural enemies, as wing dimorphism in aphids is transgenerational (Braendle et al., 2006). It has been reported that ecdysteroids and miR-9b are the molecules that transduce maternal crowding signals to embryos, leading to a high proportion of winged offspring (Vellichirammal et al., 2017; Shang et al., 2020). However, few studies have focused on the role of the wing disc degradation at the onset of the nymphal stage in forming wing dimorphism. In this study, we demonstrated that a novel gene REPTOR2, highly expressed in the first instar nymphs of wingless-destined pea aphid morph 12-h postbirth, activated autophagy in the wing disc and resulted in its degeneration as early as 30-h postbirth. Different types of PCD have been reported to activate in tissue degeneration in many insect orders (Xu et al., 2019). Wing discs of soldier ant and some sex-dependent wingless morphs of lepidopterans are degenerated by apoptotic cell death (Fujiwara and Ogai, 2001; Niitsu, 2001; Sameshima et al., 2004; Rajakumar et al., 2012; Rajakumar et al., 2018). On the other hand, apoptosis did not seem to be necessary in wing disc degradation in aphids. Instead, typical features of autophagic cell death occurred at 30-h postbirth in wingless-destined morphs. Our results supported a model concerning aphid wing disc degeneration and development: Once experienced crowding, maternal aphids increased TOR expression in first instar nymphs. This was followed by the repression of REPTOR2 and autophagic degradation in the wing disc, thereby causing a high proportion of winged offspring. For the solitary maternal group on the other hand, TOR activity was arrested in the first instar offspring, and REPTOR2 and autophagy were subsequently activated in wing disc, leading to a high proportion of wingless offspring (Figure 5). Given the high expression in thorax of wingless-destined morph at early first instar nymphal stage, it was likely that REPTOR2 transduced the TOR signaling to activate autophagic degeneration of the wing disc in a tissue-specific manner.
The TOR signaling pathway has been reported to modulate developmental plasticity in some insects (Koyama et al., 2013). For instance, activation of TOR activity in honeybees is responsible for the developmental transition from facultatively sterile workers to reproductive queens, while repressing TOR activity could induce worker characters in queen-destined individuals (Michalak et al., 2007). Transcriptomic analysis in maternal crowding of pea aphids suggested that genes functionally related to TOR signaling pathway were differentially expressed in mother aphids (Parker et al., 2021). We therefore hypothesized that the canonical pathway of TOR-REPTOR-autophagy may be responsible for developmental plasticity of the wing disc in aphids, in accordance with lower TOR abundance in wingless-destined morph. Suppression of TOR in Drosophila upregulates the expression of autophagy-related genes via phosphorylation of Ser527 and Ser530 in REPTOR (Tiebe et al., 2015). Similarly, pea aphid REPTOR2 is predicted to contain a conserved BRLZ domain that is capable of binding to GTAAACAA, the binding motif of FOXO as well (Tiebe et al., 2015). FOXO is a transcription factor that relays the insulin signaling to insect wing dimorphism. In the planthopper, activating insulin signaling transforms short-winged to long-winged morph via repressing the FOXO activity (Xu et al., 2015). Activation of FOXO inhibits the development of the long-winged morph by suppressing vestigial (Vg) expression during the wing-morph decision-making stage (Zhang et al., 2021). It appears that FOXO and REPTOR2, respectively, relay the insulin and TOR signaling cascades to control wing dimorphism in planthoppers and aphids. Alternatively, the BRLZ domain is required to form REPTOR/REPTOR-BP homo/hetero-dimer complex which conveys negative regulation of TOR signaling on genes related to autophagy (Tiebe et al., 2015). Our results are consistent with the canonical pathway regulated by TOR in that either knockdown of REPTOR2 or activation of TOR could attenuate autophagy in the wing disc. These results suggest that REPTOR2 from pea aphids served as a transcription factor to exert a negative effect on TOR-regulated autophagy initiation.
Insulin signaling is known to determine the wing-morph switch in planthoppers (Xu et al., 2015) and is also responsible for production of winged offspring regulated by miR-9b-ABCG4 cascade in brown citrus aphid (Shang et al., 2020). Furthermore, insulin could enhance TOR signaling by releasing the negative feedback of Akt on the tuberous sclerosis proteins TSC1 (hamartin) and TSC2 (tuberin) (Blagosklonny et al., 2009). It has been shown that embryos derived from pea aphids subjected to crowding treatment have higher insulin signaling (Grantham et al., 2019). This may lead to stronger TOR signaling in the first instar offspring, which could suppress wing disc autophagy and increase the proportion of winged aphids.
Gene duplication is a major source for novel gene production (Kaessmann, 2010; Prince and Pickett, 2002), which has become an important mechanism regulating phenotypic plasticity (Chen et al., 2013), including manipulating wing polyphenism in insects. For instance, a duplicated fs3 (follistatin 3) determined the wingless morph formation for male pea aphids (Li et al., 2020). Another duplicated gene InR2 (insulin receptor 2) was highly expressed in the wing bud of the planthopper that served as a regulatory switch of short-winged morph formation (Xu et al., 2015). Similarly, different from other seven aphid species with single copies of REPTOR, pea aphids possess two REPTOR genes. Although it has lost 2 exons, ApREPTOR2 has higher similarity in coding sequence with ApREPTOR1, but shares little similarity in the flanking regions. This is in agreement with the chromosomal distribution results indicating that ApREPTOR1 is located in A1 chromosome while ApREPTOR2 is in X chromosome. ApREPTOR2 is likely a duplicate from ApREPTOR1, followed by translocating from A1 to X chromosome. This may also warrant the reproductive role of sexual female in pea aphids because most sexual females are monomorphic wingless, whereas males that contain just single X chromosome are dimorphic (Brisson, 2010).
In the aphid wing polyphenism, differential expression of diverged paralogous genes could represent an adaptive response to particular environmental challenges (Shigenobu et al., 2010; Brisson et al., 2010). Evidently, ApREPTOR2-regulated wing disc degradation is significant for the winged- versus wingless-destined aphids, whereas the expression of its paralog ApREPTOR1 was unaffected by wing morphs or tissue types (Figure 3—figure supplement 1A). Knockdown of ApREPTOR1 did not affect the expression of ApREPTOR2 and ATGs and the proportion of winged aphids (Figure 3—figure supplement 1B–E). These results indicated that ApREPTOR2 rather than the canonical REPTOR activated a tissue-specific degradation in the wing disc and not a systemic PCD. Apparently, more research is necessary to further elucidate how other aphid species that contain a single REPTOR regulate autophagic degradation in the wing disc.
Materials and methods
Aphid rearing
Request a detailed protocolThe pea aphid A. pisum (strain: Ningxia Red) was originally collected from Medicago sativa in Ningxia Province and had been reared in the laboratory for over 8 years. Nymphs from the same parthenogenetic pea aphid female were reared on Vicia faba at 18–20°C, with 60% relative humidity and a photoperiod of a 16 hr: 8 hr (light/dark) cycle. To eliminate the transgenerational effects on offspring morphs, females were maintained at low density (three per plant) on V. faba seedlings for more than three generations (Sutherland, 1969). Other three strains of A. pisum, including Ningxia Green, Gansu Red, and Gansu Green, were originally collected from M. sativa in Ningxia Province and Gansu Province, respectively, and were used to test the inductive effect of two-adult contacts on the proportion of winged offspring.
Induction of winged- and wingless-destined morphs
Request a detailed protocolTo effectively induce a high proportion of winged offspring, the maternal density and contacting duration were evaluated within petri dishes (35 mm in diameter) as described by Sutherland, 1969. Four groups of aphids consisting of 1, 2, 4, or 8 wingless adults, respectively, were placed in a petri dish for 8-hr. Later, the two wingless adult groups were placed in a petri dish for 2-, 4-, and 8-hr, respectively. After treatment, each adult was then transferred to a freshly detached V. faba leaf kept in petri dishes with 1% agar, and allowed to reproduce for 24-hr. After newborns reached the fourth instar, the proportion of winged offspring from each adult was calculated. Furthermore, oviposition of each adult from the two wingless adult groups with 4-hr contact was documented for 5 days. Offspring were transferred daily to detached V. faba leaves that were changed every 3 days. Winged or wingless of each offspring was counted until the fourth instar. To eliminate the strain-specific induction, aphid strains of Ningxia Green, Gansu Red, and Gansu Green were also tested using two wingless adults contacting for 4-hr.
Developmental time and wing morphology of winged- and wingless-destined morphs
Request a detailed protocolNewly born nymphs of winged- or wingless-destined morphs were raised individually on fresh, detached leaf of V. faba in petri dishes. The nymphal developmental time and the pre-reproductive period were recorded every 12-hr. The digital microscope (Keyence VHX-5000, Osaka, Japan) was used to determine the wing morphology in nymphs of each instar.
Histology
Request a detailed protocolWinged- and wingless-destined morphs at 24-, 30-, and 36-hr postbirth were collected and placed in 4% paraformaldehyde fixative overnight at 4°C, washed in phosphate-buffered saline for three times and dehydrated at 4°C for 24-hr in 30% sucrose solution. After freezing at –20 °C in Tissue Freezing Medium (Leica), aphids were sectioned at 10 μm using a Leica CM1950 freezing microtome, and stained with H&E (Beyotime, C0105) according to the standard protocol. Images of wing discs were taken using a Nikon light microscope, equipped with a DS-Fi1c camera (Nikon), and the images were generated using the NIS-Element D software (Nikon). The relative area of wing disc was analyzed using ImageJ software.
Transmission electron microscope
Request a detailed protocolSince wing disc degeneration initiated at 30-hr postbirth, winged- and wingless-destined morphs at this time point were collected and fixed with 2.5% (vol/vol) glutaraldehyde, and washed four times in PB (0.1 M, pH 7.4). Aphids were then fixed with 1% (wt/vol) osmium tetraoxide in PB for 2-hr at 4°C, and dehydrated through a graded ethanol series (30%, 50%, 70%, 80%, 90%, 100%, and 100%, 7 min each) into pure acetone (2×10 min). Samples were infiltrated in a graded mixture (3:1, 1:1, and 1:3) of acetone and SPI-PON812 resin 16.2 g SPI-PON812, 10 g dodecyl succinic anhydride, and 8.9 g N-methylolacrylamide, then in pure resin. Aphids were embedded in pure resin with 1.5% benzyl dimethyl amine and polymerized for 12-hr at 45°C, and 48-hr at 60°C. The ultrathin sections (70 nm thick) were sectioned with microtome (Leica EM UC6, Austria), and then double-stained by uranyl acetate and lead citrate, and imaged under transmission electron microscope (FEI Tecnai Spirit 120kV, USA). Ultrathin sections were used to count the cell number that had signs of autophagy (mitophagy, cytoplasmic vacuoles, large autophagosomes sequestering visible remnants of organelles, and membranous whorls) (Simonet et al., 2018). The percentage of autophagic cells was determined by randomly counted cells in the wing disc.
Immunofluorescence detection and TUNEL assays
Request a detailed protocolThe frozen sections of aphids with different treatments were rinsed three times in TBST (TBS with 0.05% Tween-20). TUNEL assays were carried out by One Step TUNEL Apoptosis Assay Kit (Beyotime, C1089) according to the manufacturer’s protocol. For immunofluorescence, frozen sections of aphids were blocked with SuperBlock T20 (Pierce) for 20 min. The samples were incubated with the primary antibody (anti-ATG8, 1:200) overnight at 4°C, rinsed three times in TBST, and incubated with the secondary antibody (1:500) at room temperature for 2-hr. A negative control was performed in each independent experiment. The samples were rinsed for three times in TBST, and mounted in Fluoroshield Mounting Medium with DAPI (Abcam). Sections were imaged using a Zeiss LSM710 confocal microscope (Zeiss, Germany). The polyclonal antibody rabbit anti-ATG8 was kindly provided by Professor Le Kang (Institute of Zoology, Chinese Academy of Sciences). Anti-rabbit IgG (H+L), F(ab')2 Fragment (Alexa Fluor 594 Conjugate, 8889S) was purchased from CST. Relative intensities of TUNEL and ATG8 were analyzed using the ImageJ software.
Fluorescence in situ hybridization
Request a detailed protocolFISH was performed with the technique modified by Kliot et al., 2014. After washing by TBST (TBS with 0.2% Triton-X) for 10 min, the frozen sections were rinsed for three times in hybridization buffer (20 mM Tris-HCl, pH 8.0, 0.9 M NaCl, 0.01% [wt/vol] sodium dodecyl sulfate, and 30% [vol/vol] formamide) for pre-hybrid (without the probe). The frozen sections of aphids were then hybridized overnight in hybridization buffer containing 10 pmol of the fluorescent RNA probe (conjugated with Cy3). Afterward, the frozen sections were rinsed for three times in TBST, and mounted in Fluoroshield Mounting Medium with DAPI (Abcam). Sections were imaged using a Zeiss LSM710 confocal microscope (Zeiss, Germany).
Quantitative PCR
Request a detailed protocolTotal RNA of aphid samples was extracted by TRIzol Reagent (Thermo Fisher) according to the manufacturer’s protocol, and cDNA was synthesized using the FastQuant RT Kit with gDNase (Tiangen). The qPCR reactions were performed with the PowerUp SYBR Green Master Mix (Applied Biosystems). RT-qPCR reactions were carried out on the QuantStudio 12 K Flex Real-Time PCR System (ABI, A25742) as follows: 30 s at 95°C; followed by 40 cycles of 10 s at 95°C, 30 s at 60°C; and finally one cycle of 15 s at 95°C, 60 s at 60°C, and 15 s at 95°C. The melting curves were used to determine the specificity of the PCR products. The housekeeping gene actin (ACYPI000064) was used as the internal qPCR standard to analyze gene expression (Simonet et al., 2018). The relative level of each target gene was standardized by comparing the copy numbers of target mRNA with copy numbers of actin. Data were analyzed by the 2−ΔΔCT relative quantification method. Specific primers for each gene were designed from the aphid sequences using Primer Premier 6 software. The availability of each pair of primers had been tested in the preliminary experiments. The oligonucleotides are listed in Supplementary file 1.
We tested for enrichment of expression of REPTOR1 or REPTOR2 in aphid heads, thorax, and abdomen using qPCR. Aphid dissection was performed with a technique modified by Chen et al., 2015. Winged- and wingless-destined aphids at 24-hr postbirth were collected in 2 ml RNAlater (Thermo Fisher Scientific) and kept at 4℃ overnight, and 100 μl 0.1% PBST was added to help RNAlater penetration. Head, thorax, and abdomen from around 80 aphids were dissected in cold RNAlater. RNA extraction and qPCR of dissected tissues were carried out as above.
Pharmacological experiments and RNA interference via feeding
Request a detailed protocolTwenty newborn aphids were placed in cylindrical containers, 20 mm diameter × 20 mm height. Each container was provided with a 100 μl artificial diet (Haribal and Jander, 2015). The autophagy inhibitor 3-MA (Selleck, S2767) was dissolved in double distilled water and incorporated into diet at 10 μM. The autophagy agonist rapamycin (Selleck, S1039) was dissolved in ethanol and mixed in diet at a concentration of 10 μM, and an equivalent amount of ethanol was used as a feeding control. The TOR agonist MHY1485 (Selleck, S7811) was dissolved in DMSO and mixed in diet at a concentration of 10 μM, and an equivalent amount of DMSO was used as a feeding control. Diet with 50% amino acids (each amino acid was supplied at 50% of standard formulation) was used to repress TOR activity, and a normal diet was used as a control. All chemicals were administered via 24-hr feeding.
The T7 RiboMAX Express RNAi System (Promega, P1700) was used to synthesize double-stranded RNA (dsRNA) according to the manufacturer’s protocol. dsRNAs, that is, dsREPTOR1, dsREPTOR2, dsTOR (5 μl, 10 μg/μl), or dsGFP control were added to 95 μl artificial diet to feed aphids. The interference efficiencies were measured by qPCR.
Co-silencing of TOR and REPTOR2
Request a detailed protocolNewborn nymphs were fed with a dsRNA mixture of dsTOR+dsREPTOR2, and dsTOR+dsGFP served as the control (50 μg for each dsRNA). Nymphs at 30-hr postbirth were collected for immunofluorescence detection. In addition, each winged or wingless aphid was counted until the fourth instar.
RNA-seq and data analysis
Request a detailed protocolTo determine DEGs in winged- and wingless-destined morphs, the whole body of 20 aphids at 24-hr postbirth was collected for further RNA extraction. RNA purity and quantification were evaluated using the NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific). RNA integrity was assessed using the Agilent 2100 Bioanalyzer (Agilent Technologies). The filtered data were obtained by removing low-quality reads from the raw data with HTseq-count software (Anders et al., 2014). After quality control, paired-end clean reads were aligned to the reference genome downloaded from NCBI (ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/005/508/785/GCF_005508785.1_pea_aphid_22Mar2018_4r6ur/GCF_005508785.1). FPKM of each gene was calculated using Cufflinks (v2.1.1). Differential expression analysis was performed using the DESeq (2012) R package (Anders and Huber, 2012). log2 |fold change|>0.58, p value<0.05 was set as the threshold for significantly differential expression. Sequence data have been deposited in the National Center for Biotechnology Information’s Sequence Read Archive, https://www.ncbi.nlm.nih.gov/sra (accession no. PRJNA857589).
Analysis of REPTOR
Request a detailed protocolTwo REPTOR genes were identified in the A. pisum reference genome (GCF_005508785.1): LOC100168197 (REPTOR1) and LOC103309809 (REPTOR2).
The REPTOR homologs in the M. cerasi, M. persicae, M. secchari, A. gossyppii, R. padi, R. maidis, and E. lanigerum were identified by BLASTp using pea aphid REPTOR1. Coding sequences of homologous genes of REPTOR, including LOC100168197 (REPTOR1, A. pisum), LOC103309809 (REPTOR2, A. pisum), Mca07599.t1 (M. cerasi), g26487.t1 (M. persicae), LOC112600553 (M. secchari), LOC114124537 (A. gossyppii), g1773.t1 (R. padi), LOC113558611 (R. maidis), and jg23767.t1 (E. lanigerum), were obtained from NCBI and AphidBase. Sequences were aligned with MUSCLE v3.8.1551 (Edgar, 2004). The phylogenetic analysis was performed in IQ-TREE v2.1.4-beta with HKY+F+G4 model (Minh et al., 2020) and the phylogenetic tree was visualized in FigTree v1.4.4 (http://tree.bio.ed.ac.uk/software/figtree/). Genomic sequences of orthologous regions of REPTOR were obtained from M. persicae (v2.0), R. maidis (v1.0), and E. lanigerum (v1.0), which have chromosome-level genome assemblies available at Aphidbase. Sequences were compared with NCBI blast v2.2.26 and homologous sequences were linked by the gray color when sequences matched to positive strand and by the red color when sequences matched to negative strand. Dot plots were generated using MUMmer v3.23. MapGene2Chrom v2.1 was used to analyze the chromosomal distribution of two pea aphid REPTOR genes (http://mg2c.iask.in/mg2c_v2.1/). Conserved domains of REPTOR were identified by SMART (http://smart.embl-heidelberg.de/). The amino acid sequences of REPTOR homologs from A. pisum, D. melanogaster, and H. sapiens were downloaded from GenBank and aligned with ApREPTOR2 using MUSCLE v3.8.1551.
Statistical analysis
Request a detailed protocolStatistical analyses were performed with SPSS software (Chicago, IL) and GraphPad software. All data were checked for normality by the Shapiro-Wilk test. Student’s t tests were used to analyze the relative area of the wing disc, percentage of autophagic cells, the number of autophagosomes per wing disc, intensity ratios of ATG8/DAPI, TUNEL/ DAPI, REPTOR2/DAPI, and TOR/ DAPI, percentage of winged aphids, and quantification of gene expression for two-group comparisons. ANOVA was used to analyze the effects of maternal density and contacting duration on the proportion of winged offspring, and proportion of winged offspring in each day for 4-hr contacting treatment with two female adults.
Data availability
Sequencing data have been deposited in SRA under accession ID PRJNA857589.
-
NCBI BioProjectID PRJNA857589. Wing dimorphism in Acyrthosiphon pisum.
References
-
SoftwareDifferential expression of RNA-seq data at the gene level – the deseq packageDESeq Package.
-
Aphid wing dimorphisms: Linking environmental and genetic control of trait variationPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 365:605–616.https://doi.org/10.1098/rstb.2009.0255
-
New genes as drivers of phenotypic evolutionNature Reviews. Genetics 14:645–660.https://doi.org/10.1038/nrg3521
-
Relationship between growth arrest and autophagy in midgut programmed cell death in DrosophilaCell Death and Differentiation 19:1299–1307.https://doi.org/10.1038/cdd.2012.43
-
Muscle: Multiple sequence alignment with high accuracy and high throughputNucleic Acids Research 32:1792–1797.https://doi.org/10.1093/nar/gkh340
-
Ecdysteroid-induced programmed cell death and cell proliferation during pupal wing development of the silkworm, Bombyx moriDevelopment Genes and Evolution 211:118–123.https://doi.org/10.1007/s004270000131
-
Stable isotope studies reveal pathways for the incorporation of non-essential amino acids in acyrthosiphon Pisum (pea aphids)The Journal of Experimental Biology 218:3797–3806.https://doi.org/10.1242/jeb.129189
-
Dispersal polymorphisms in insectsAnnual Review of Ecology and Systematics 11:95–118.https://doi.org/10.1146/annurev.es.11.110180.000523
-
Origins, evolution, and phenotypic impact of new genesGenome Research 20:1313–1326.https://doi.org/10.1101/gr.101386.109
-
Role of TOR signaling in aging and related biological processes in Drosophila melanogasterExperimental Gerontology 46:382–390.https://doi.org/10.1016/j.exger.2010.11.036
-
Adaptation to nicotine in the facultative tobacco-feeding hemipteran bemisia tabaciPest Management Science 70:1595–1603.https://doi.org/10.1002/ps.3739
-
IQ-TREE 2: New models and efficient methods for phylogenetic inference in the genomic eraMolecular Biology and Evolution 37:1530–1534.https://doi.org/10.1093/molbev/msaa015
-
BookForms and morphs of aphidsIn: Minks AK, Harrewijin P, editors. Aphids, their biology, natural enemies, and control. Elsevier. pp. 27–50.
-
Mechanisms of growth and homeostasis in the Drosophila wingAnnual Review of Cell and Developmental Biology 25:197–220.https://doi.org/10.1146/annurev.cellbio.24.110707.175242
-
Wing degeneration due to apoptosis in the female of the winter moth nyssiodes lefuarius (lepidoptera, geometridae)Entomological Science 4:1–7.
-
A laterally transferred viral gene modifies aphid wing plasticityCurrent Biology 29:2098–2103.https://doi.org/10.1016/j.cub.2019.05.041
-
Wing plasticity and associated gene expression varies across the pea aphid biotype complexEvolution; International Journal of Organic Evolution 75:1143–1149.https://doi.org/10.1111/evo.14174
-
Splitting pairs: the diverging fates of duplicated genesNature Reviews. Genetics 3:827–837.https://doi.org/10.1038/nrg928
-
Nutrient-Regulated phosphorylation of Atg13 inhibits starvation-induced autophagyThe Journal of Biological Chemistry 291:6026–6035.https://doi.org/10.1074/jbc.M115.689646
-
The evolution of flightlessness in insectsEcological Monographs 60:389–421.https://doi.org/10.2307/1943013
-
Wing dimorphisms and the evolution of migratory polymorphisms among the InsectaAmerican Zoologist 31:243–251.https://doi.org/10.1093/icb/31.1.243
-
Multiple amino acid sensing inputs to mTORC1Cell Research 26:7–20.https://doi.org/10.1038/cr.2015.146
-
The role of crowding in the production of winged forms by two strains of the pea aphid, acyrthosiphon PisumJournal of Insect Physiology 15:1385–1410.https://doi.org/10.1016/0022-1910(69)90199-1
-
Ecdysone controlled cell and tissue deletionCell Death and Differentiation 27:1–14.https://doi.org/10.1038/s41418-019-0456-9
-
Physiology and ecology of dispersal polymorphism in insectsAnnual Review of Entomology 42:207–230.https://doi.org/10.1146/annurev.ento.42.1.207
Article and author information
Author details
Funding
Ministry of Science and Technology of the People's Republic of China (the National Key R&D Program of China (no. 2022YFD1400800))
- Yucheng Sun
Chinese Academy of Sciences (the Strategic Priority Research Program of the Chinese Academy of Sciences (no. XDPB16))
- Yucheng Sun
National Natural Science Foundation of China (the National Natural Science Foundation of China (no. 32250002))
- Yucheng Sun
National Natural Science Foundation of China (the National Natural Science Foundation of China (no. 31970453))
- Yiyang Yuan
National Natural Science Foundation of China (the National Natural Science Foundation of China (no. 31870394))
- Huijuan Guo
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors thank our colleagues Prof. Jinfeng Chen for his support in genetic analysis and Prof. Chenzhu Wang for his advises in project coordination and manuscript writing. The authors appreciate TEM sample preparing and analysis support from Can Peng, Center for Biological Imaging, Institute of Biophysics, CAS. This project was supported by the National Key R&D Program of China (No. 2022YFD1400800), the Strategic Priority Research Program of the Chinese Academy of Sciences (No. XDPB16), and the National Natural Science Foundation of China (Nos. 32250002, 31970453, and 31870394).
Copyright
© 2023, Yuan et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 985
- views
-
- 279
- downloads
-
- 7
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
Cells called alveolar myofibroblasts, which have a central role in the development of the lung after birth, receive an orchestrated input from a range of different signaling pathways.
-
- Developmental Biology
We have generated single cell transcriptomic atlases of vomeronasal organs (VNO) from juvenile and adult mice. Combined with spatial molecular imaging, we uncover a distinct, previously unidentified class of cells that express the vomeronasal receptors (VRs) and a population of canonical olfactory sensory neurons in the VNO. High-resolution trajectory and cluster analyses reveal the lineage relationship, spatial distribution of cell types, and a putative cascade of molecular events that specify the V1r, V2r, and OR lineages from a common stem cell population. The expression of vomeronasal and olfactory receptors follow power law distributions, but there is high variability in average expression levels between individual receptor and cell types. Substantial co-expression is found between receptors across clades, from different classes, and between olfactory and VRs, with nearly half from pairs located on the same chromosome. Interestingly, the expression of V2r, but not V1r, genes is associated with various transcription factors, suggesting distinct mechanisms of receptor choice associated with the two cell types. We identify association between transcription factors, surface axon guidance molecules, and individual VRs, thereby uncovering a molecular code that guides the specification of the vomeronasal circuitry. Our study provides a wealth of data on the development and organization of the accessory olfactory system at both cellular and molecular levels to enable a deeper understanding of vomeronasal system function.