Amino acid transporters implicated in endocytosis of Buchnera during symbiont transmission in the pea aphid
© The Author(s) 2016
Received: 25 June 2016
Accepted: 4 November 2016
Published: 21 November 2016
Many insects host their obligate, maternally transmitted symbiotic bacteria in specialized cells called bacteriocytes. One of the best-studied insect nutritional endosymbioses is that of the aphid and its endosymbiont, Buchnera aphidicola. Aphids and Buchnera are metabolically and developmentally integrated, but the molecular mechanisms underlying Buchnera transmission and coordination with aphid development remain largely unknown. Previous work using electron microscopy to study aphid asexual embryogenesis has revealed that Buchnera transmission involves exocytosis from a maternal bacteriocyte followed by endocytotic uptake by a blastula. While the importance of exo- and endocytic cellular processes for symbiont transmission is clear, the molecular mechanisms that regulate these processes are not known. Here, we shed light on the molecular mechanisms that regulate Buchnera transmission and developmental integration.
We present the developmental atlas of ACYPI000536 and ACYPI008904 mRNAs during asexual embryogenesis in the pea aphid, Acyrthosiphon pisum. Immediately before Buchnera invasion, transcripts of both genes were detected by whole-mount in situ hybridization in the posterior syncytial nuclei of late blastula embryos. Following Buchnera invasion, expression of both genes was identified in the region occupied by Buchnera throughout embryogenesis. Notably during Buchnera migration, expression of both genes was not concomitant with the entirety of the bacterial mass but rather expression colocalized with Buchnera in the anterior region of the bacterial mass. In addition, we found that ACYPI000536 was expressed in nuclei at the leading edge of the bacterial mass, joining the bacterial mass in subsequent developmental stages. Finally, quantitative reverse transcription real-time PCR suggested that early in development both transcripts were maternally provisioned to embryos.
We venture that ACYPI000536 and ACYPI008904 function as nutrient sensors at the site of symbiont invasion to facilitate TOR-pathway-mediated endocytosis of Buchnera by the aphid blastula. Our data support earlier reports of bacteriocyte determination involving a two-step recruitment process but suggest that the second wave of recruitment occurs earlier than previously described. Finally, our work highlights that bacteriocyte-enriched amino acid transporter paralogs have additionally been retained to play novel developmental roles in both symbiont recruitment and bacteriome development.
KeywordsAmino acid transporter Aphid development Coevolution Host/symbiont developmental integration Symbiosis
A compelling challenge in symbiosis research concerns identifying and characterizing the mechanisms of host/symbiont metabolic and developmental integration . Recently, we began to address this challenge by studying host/symbiont developmental integration in the pea aphid Acyrthosiphon pisum and its obligate endosymbiont Buchnera aphidicola . Like all aphids, A. pisum feeds exclusively on plant phloem sap, a diet that is rich in sugar but limited in provision of essential amino acids and vitamins [3–5]. To overcome the nutritional deficiency of their diet, aphids anciently acquired B. aphidicola, an endosymbiont that metabolically compensates the nutritional shortcomings of their diet [1, 6, 7]. Buchnera is transmitted from generation to generation through transovarial transmission , and remarkably aphids cured of Buchnera infection fail to reproduce [9–12].
Aphids are cyclical parthenogens, capable of both sexual and asexual modes of reproduction. In asexual aphids, one ovariole contains a series of embryos at various developmental stages . Parthenogenetic embryogenesis has been described by Will  and Miura et al. . The parthenogenetic embryo originates from a maturation division of the oocyte by skipping the reduction division of meiosis I, producing a polar body and an oocyte nucleus [8, 15]. The oocyte nucleus is considered to be the incipient parthenogenetic embryo that divides synchronously five times in a common cytoplasm to form the syncytial embryo with 32 syncytial nuclei. Cellularization of the syncytial nuclei then generates the blastula embryo. Primordial germ cells that are specified in the posterior region of the blastula separate the remaining uncellularized syncytial nuclei into two locations: (1) the central syncytium (in the center of the blastoderm) and (2) the posterior syncytium (at the posterior end of the blastoderm) .
Developmental integration of Buchnera into host aphids starts with transmission of a small population of Buchnera from maternal bacteriocytes, the specialized aphid cells accommodating Buchnera, to a blastula . Transmission of Buchnera involves exocytosis from maternal bacteriocytes followed by endocytosis by the developing blastula. During exocytosis, individual Buchnera lose their host-derived symbiosomal membrane (a perisymbiont membrane that encloses individual Buchnera cells) and migrate naked across the extracellular space between the maternal bacteriocyte and the posterior pole of the blastula to reacquire their symbiosomal membrane by endocytotic passage into the blastoderm (stage 7 of development) [8, 16]. Following transmission, the Buchnera population, together with aphid uncellularized syncytial nuclei, are subdivided into uninucleate bacteriocyte cells . After this, bacteriocyte cells proliferate , and maturation of the bacteriome proceeds . Aphid bacteriomes, the organ of symbiosis, comprise two cell types: bacteriocytes that house Buchnera and intercalating sheath cells that do not house Buchnera [18, 19]. Localization of the transcription factor Distal-less (Dll) reveals that bacteriocytes are recruited through two steps. In the first step, bacteriocyte cells originate from the syncytial nuclei that localize prior to Buchnera invasion in the posterior syncytium of the blastoderm embryo . In the second step, bacteriocyte cells originate from approximately 40–60 Dll-expressing cells that migrate across the germband and intercalate between the original bacteriocytes during germband extension . In contrast to the original proposal by Braendle et al.  that bacteriocytes are recruited in two steps, Koga et al.  hypothesized, on the basis of the smaller nuclei of the second group of cells and their localization between the original bacteriocyte population, that the second group of cells may in fact be sheath cells and not bacteriocytes. In late embryogenesis, aphid embryos undergo katatrepsis when the head is repositioned from the posterior to the anterior pole of the egg chamber resulting in inversion of the body axis; during this maturation, bacteriocyte cells are positioned to the dorsal side of the embryo where they later form the mature bacteriome [2, 8, 16].
Work to date in asexually reproducing A. pisum has demonstrated that specification of bacteriocytes occurs independently of symbiont invasion from mother to progeny  and that Buchnera transmission involves exo- and endocytotic cellular processes . Bacteriocytes are widely found in sap-feeding and blood-feeding insect lineages, and yet little is known about the molecular network underlying bacteriocyte morphogenesis or the mechanisms that regulate transmission and developmental integration of symbionts.
Here, by studying the expression of two aphid amino acid transporter mRNAs during asexual embryogenesis, we gain insight into the molecular mechanisms that regulate transmission and developmental integration of Buchnera into aphid hosts. Amino acid transporters are the membrane proteins responsible for cellular nutrient exchange and nutrient sensing . As nutrient sensors, amino acid transporters act as transceptors: dual-function molecules that both transport amino acids and function as receptors sensing extracellular amino acid availability and modulating intracellular nutrient-responsive signaling pathways [20–23]. Amino acid transporters have been implicated in a broad diversity of cellular and developmental functions that include control of cell growth [24–27], cell proliferation [28–30], regulation of translational capacity , lysosome biogenesis , and coordination of endocytosis/phagocytosis [26, 33]. The A. pisum genome encodes 40 putative amino acid polyamine organocation superfamily amino acid transporters, eight of which are expressed in bacteriocyte cells . So far, bacteriocyte-enriched amino acid transporters have been implicated in regulating host and symbiont mediated amino acid biosynthesis in bacteriocytes , and comparative work has demonstrated that retention of duplicated amino acid transporters can result from selection to function in symbiosis [36–38]. Within A. pisum bacteriocytes, the two amino acid transporters that are most highly expressed are the eukaryotic-specific amino acid/auxin permease (AAAP) family (TC# 2.A.18) member, ACYPI000536, and the amino acid/polyamine/organocation (APC) family (TC# 2.A.3) member, ACYPI008904 . ACYPI000536 and ACYPI008904 are gene identifiers that have been assigned in the ACY rthosiphon PI sum reference genome available at AphidBase . To enhance readability in this paper, we will refer to ACYPI000536 as AAAP-536 and ACYPI008904 as APC-8904. Using whole-mount in situ hybridization and quantitative PCR, we present the developmental atlas of AAAP-536 and APC-8904 mRNAs during asexual embryogenesis. In interpreting the developmental atlas of AAAP-536 and APC-8904 expression, we propose that these amino acid transporters play novel roles in endocytosis of Buchnera by aphid blastula and in bacteriocyte development during asexual aphid embryogenesis.
Acyrthosiphon pisum line LSR1  and a Taiwanese line NTU  were maintained on Vicia faba and incubated at 20 °C under a 16-h light/8-h dark cycle. Oocytes and embryos for whole-mount in situ hybridization were dissected from wingless adult females in phosphate-buffered saline (PBS; 10 mM phosphate buffer, 154 mM NaCl, pH 7.4; Sigma-Aldrich). Oogenesis and embryogenesis staging was according to Miura et al.  and germ cell locations according to Chang et al. .
Cloning of AAAP-536 and APC-8904 gene fragments
Primer pairs for amplifying AAAP-536 and APC-8904 gene fragments were designed via MacVector version 7.2.2 (Accelrys Inc.) according to the gene sequences in AphidBase 2.1 . Primers designed to AAAP-536 coding sequence were further checked using Primer-BLAST against the pea aphid genome Acyr_2.0 assembly, confirming that each primer hit only its isoform targets (amplicon located at 1719–2387 of XM_008185850.2 and 1150–1818 of NM_001246304.1). Due to the high similarity of APC-8904 nucleotide sequence to that of its paralogs: 82.56% coding sequence similarity to ACYPI008323 and 76.43% coding sequence similarity to ACYPI002633, we designed APC-8904 gene-specific primers for hybridization to the 3′ untranslated region (UTR) of APC-8904. Nucleotide similarity in the 3′ UTR of APC-8904 and its paralogs is much lower such that similarity between APC-8904 and ACYPI008323 is 66.32% across the 193 base pairs (bp) of alignable 3′ UTR sequence and that between APC-8904 and ACYPI002633 is 37.99% across the 766 bp of alignable 3′ UTR sequence. For AAAP-536, 669 bp of open reading frame (ORF) was amplified from a plasmid containing the gene’s full-length cDNA  with forward primer 5′-TGGGATGCTAAGTATTCACTTCC-3′ and reverse primer 5′-TTCAATGGTAGAATTTATAGTGT-3′. For APC-8904, 566 bp of the 3′UTR was amplified from whole insect LSR1 cDNA  with forward primer 5′-TCCTCCCCTATGTTTCCACG-3′ and reverse primer 5′-TTCTCAGCGAAGACACACCGTC-3′. Amplification of gene fragments was performed starting with an initial denaturation at 94 °C for 30 s, followed by 30 cycles at 94 °C for 10 s, 50 °C for 30 s, 72 °C for 1 min and a final extension at 72 °C for 5 min. Both amplicons were subcloned into pGEM-T vector (Promega), and the insert sequences and orientations were confirmed by DNA sequencing before RNA probe synthesis.
Whole-mount in situ hybridization and microscopy
Hybridizations in situ with digoxigenin (DIG)-labeled RNA probes to whole-mount ovaries were performed as described by Chang et al. . The antisense and sense DIG-labeled RNA probes were transcribed with SP6 or T7 RNA polymerase (New England Biolabs). Transcription was carried out from 1 μg linearized plasmid DNA in 20 μl reaction mixtures containing 2 μl 10 × RNA polymerase reaction buffer (400 mM Tris–HCl, 60 mM MgCl2, 20 mM spermidine, 10 mM dithiothreitol, pH 7.9; New England Biolab), 2 μl 10 × DIG RNA labeling mix (Roche), 0.5 μl RNAse inhibitor (Invitrogen), and 2 μl RNA polymerase (New England Biolabs). The reaction mixtures were incubated for 4 h at 37 °C, and the products were verified by size on a 1.5% agarose gel containing 2.2 M formaldehyde. After adding sodium acetate to 0.3 M, nucleic acids were precipitated with ethanol until further use.
Dissected ovaries from A. pisum line LSR1 containing developing oocytes and embryos were fixed in 3.8% formaldehyde (VWR) in PBS at 4 °C overnight. Each experiment included antisense riboprobe (experimental group) and sense riboprobe (control group) treatments, and the working concentration for each riboprobe was 2.0 ng/μl. Other steps followed the protocol of Chang et al. . Probe hybridizations for AAAP-536 and APC-8904 were carried out at 68 and 60 °C, respectively, and specimens were incubated overnight at 4 °C with anti-DIG Fab fragments conjugated to alkaline phosphatase (Roche), diluted 1:500 in 1× blocking reagent (Roche). Nitroblue tetrazolium (NBT)/5-bromo-4-chloro-3-indolyl phosphate (BCIP) (Roche) was applied to detect the hybridized probes. We terminated the substrate reaction concurrently for experimental (antisense riboprobe) and control (sense riboprobe) samples based on the relative intensity of signal in the region of embryos occupied by Buchnera in the antisense riboprobe (experimental) and sense riboprobe (control) treatments. Following signal development, samples were counterstained with 2 μg/ml of 4′,6-diamidino-2-phenylindole (DAPI) (Sigma-Aldrich) at room temperature for 1 h and mounted in 70% glycerol (Sigma-Aldrich) in PBS at 4 °C overnight. Photographs of NBT/BCIP staining were acquired with a Zeiss Axiovert 200 microscope connected to a Zeiss AxioCam ICc1 camera. Photographs of DAPI staining were taken with the Leica SP5 laser scanning confocal microscope in the University of Miami, Department of Biology Microscopy Core Facility. The experiment was performed 13 times for AAAP-536 and 6 times for APC-8904. Samples collected from each developmental stage were counted to report on the repeatability of our results; count data are available in Additional file 1: Table S1 for AAAP-536 and Additional file 2: Table S2 for APC-8904.
Total RNA for each biological replicate was prepared from asexual adult female ovarioles (equivalent to 600–700 ovarioles) from A. pisum line NTU. Each ovariole contains germaria, oocytes and developing embryos at various stages of development enveloped within egg chambers. Dissection was performed under a Leica EZ4 Stereo Microscope with non-removable 16× eyepieces at 40× magnification on a 12-well watch glass with one individual adult dissected into each well. To avoid contamination, ovaries were separated from other somatic tissues and the somatic tissues were removed using a 200-μl pipetman (Eppendorf). Ovaries were washed several times with ice-cold PBS until no maternal bacteriocytes or other somatic tissue remained in the well. Washes were performed with great care to avoid damaging the ovary structure. Under 56× magnification on the Leica Stereo Microscope, using a black enameled pin (diameter 0.25 mm; Entomoravia), the egg chambers were separated on the follicular stalk between embryos that were stage 5 or younger and stage 6 or older. Importantly, at developmental stage 6, embryos are clearly identifiable morphologically by the appearance of their primordial germ cells and enlarged follicle cells . Two sets of embryos were collected: set A that included germaria and stage 1 to stage 5 embryos and set B that included stage 6 to stage 20 embryos. Total RNA was extracted from embryo sets using a RNeasy Mini kit (Qiagen). To ensure that total RNA was free of DNA, a DNase I digestion was performed twice during RNA purification with the RNase-Free DNase Set (Qiagen). The first treatment was an on-column digestion, and DNase was removed in subsequent wash steps. The second treatment was performed after RNA purification by an in-solution digestion prior to RNA cleanup. Reaction setup followed the manufacturer’s instructions (Qiagen). Total RNA was quantified using an Eppendorf Biophotometer spectrophotometer. Two hundred nanograms of total RNA was reverse-transcribed into cDNA with QuantiTect Reverse Transcription kit (Qiagen) in standard 20 μl reaction mixture. All experiments were performed with three independent biological replicates; RNA for each biological replicate was prepared using an individual set of dissected embryos from the NTU A. pisum line.
Primer pairs used to detect gene expression as well as amplicon sizes and amplification efficiencies are given in Additional file 3: Table S3. Ten microliter PCR mixture contained 1× QuantiFast SYBR Green PCR Master Mix (Qiagen), 200 nM of forward primer, 200 nM of reverse primer, and cDNA derived from 9.2 ng of total RNA from set A or set B tissue samples. Each sample was amplified using triplicate technical replicates. Each plate included no reverse transcription controls and no template controls for each primer set. Samples were run on a Rotor-Gene Q real-time PCR cycler (Qiagen). The thermal cycling conditions were started by activating the enzyme reaction with 5 min at 95 °C followed by 40 cycles of 95 °C for 15 s, 52 °C for 15 s, and 60 °C for 20 s. A final melt-curve step was included post-PCR to confirm the absence of any non-specific amplification using Rotor-Gene Q series software (version 1.7). We used the programs geNorm [44, 45], Normfinder , and BestKeeper  to assess the stability of five control genes (GAPDH , EF1α , RPL7 , RPL32 , and βTUB ) between the samples of set A and set B. Based on the gene stability analyses, GAPDH, EF1α, and RPL32 were selected as endogenous controls. Gene expression was compared to embryos before stage 5 using 2−∆∆CT methodology  and normalized to the geometric mean of GAPDH, EF1α, and RPL32 .
Dissected ovaries from A. pisum line NTU containing developing oocytes and embryos were fixed in 4% paraformaldehyde (Sigma-Aldrich) in PBS at room temperature for 20 min. Following fixation, ovaries were washed with 0.2% of Triton-X 100 (Sigma-Aldrich) in PBS three times, 10 min for each wash. Ovaries were then incubated with 2 μg/ml of DAPI (Sigma-Aldrich) at room temperature for 1 h and mounted in 70% glycerol (Sigma-Aldrich) in PBS at 4 °C overnight. Photographs of DAPI staining were taken with the Leica SP5 laser scanning confocal microscope at National Taiwan University.
The developmental atlas of AAAP-536 and APC-8904 expression with reference to the symbiotic bacterial mass
Unlike AAAP-536, we observed ubiquitous background signals of APC-8904 in the germarium, the newly segregated oocyte (Fig. 1g), and mature divided oocyte (Fig. 1h). In syncytial embryos and early blastoderm embryos, transcripts of APC-8904 were randomly distributed (Fig. 1i–k, indicated with arrowheads). Similar to AAAP-536, in late stage 6 blastoderm embryos, we observed preferential localization of APC-8904 transcripts in the posterior syncytium (Fig. 1l); importantly, these preferential posterior syncytium APC-8904 signals were not found in the sense control group (Fig. 1m, n). Similar to AAAP-536, we observed transcripts of APC-8904 in the region occupied by the bacterial mass: during Buchnera invasion (Fig. 2e–g), bacteriocyte formation (stages 10–14) (Figs. 2h, 3f–i), and bacteriome maturation (stages 15–20) (Figs. 3j, 4f–i).
To summarize, we highlight the similarities and differences in expression of AAAP-536 and APC-8904 mRNAs. With respect to similarities, (1) preferential expression of both transporter genes was first recognized in nuclei of the posterior syncytium of the blastula (developmental stage 6 in Fig. 1f, l); (2) transcripts of both genes did not always colocalize with all Buchnera cells (Fig. 2b, c, f, g, h, indicated with arrows), but rather both appeared at the Buchnera migration front such that we observed colocalization of both transcripts and Buchnera in the anterior region of the bacterial mass (stages 8–9); and (3) signals for both transporters were detected during Buchnera invasion, bacteriocyte formation, and bacteriome maturation (Figs. 2, 3, 4). In contrast, where AAAP-536 transcripts were detected in the periphery of nuclei outside the bacterial region during Buchnera invasion (Fig. 2b, c, indicated with arrowheads), transcripts of APC-8904 were not (Fig. 2f, g).
AAAP-536 and APC-8904 were expressed very early in development
Developmental integration of Buchnera into aphid hosts
Aphid work that targets elucidating the mechanisms of host/endosymbiont developmental integration has shown that: (1) specification of bacteriocytes involves a two-step determination process , (2) cellular mechanisms of Buchnera transmission from maternal bacteriocytes to embryos involve exocytotic and endocytotic cellular processes , and (3) maternal and embryonic bacteriomes are not equivalent [2, 54, 55]. Here, we extend work elucidating the mechanisms of host/endosymbiont developmental integration in aphids by describing the developmental atlas of AAAP-536 and APC-8904 mRNAs with respect to bacteriocyte specification and Buchnera transmission. In this study, we demonstrated that AAAP-536 and APC-8904 are expressed at the site of Buchnera invasion ahead of Buchnera transmission from mother to progeny (Fig. 1f, l, indicated with arrowheads). Following Buchnera invasion, we found that AAAP-536 and APC-8904 expression was not concomitant with the entirety of the bacterial mass, but rather expression of AAAP-536 and APC-8904 overlaps with only the anterior part of the bacterial mass in association with the syncytial nuclei (Fig. 2b, c, f–h). In addition, we typically observed signal from AAAP-536 in the periphery of nuclei outside the bacterial mass (Fig. 2b, c, indicated with arrowheads); remarkably these nuclei localize to the areas that the bacterial mass moves to in the next developmental stage (Figs. 2b–d, 3a, b). Our results are consistent with the observations of Braendle et al.  and the hypothesis of Miura et al.  that the developmental origin of bacteriocytes involves a two-step determination process, the first from nuclei in the posterior syncytium and the second from the nuclei in the central syncytium. In contrast to Braendle et al. , who observed development of the second set of bacteriocytes at stage 13, our data suggest that recruitment of nuclei to the second set of bacteriocytes occurs as early as stage 8 of development (Fig. 2b, c). Alternatively, it is possible that AAAP-536 and Dll, the marker used by Braendle et al. , are tracking the recruitment of two different sets of cells. By synthesizing the work of Braendle et al.  and Koga et al.  with the work we present here, we propose that the nuclei recruited at stage 8, defined by AAAP-536 expression, are destined to form bacteriocyte cells, where the nuclei recruited at stages 13–15 are destined to form sheath cells (cells that are not occupied by Buchnera with nuclei much smaller than those of bacteriocytes [16, 17]).
While the clear and discrete expression of AAAP-536 and APC-8904 mRNA at the site of Buchnera invasion prior to Buchnera invasion suggests that both genes are involved in Buchnera recruitment to the blastula, the equally clear and discrete expression of AAAP-536 and APC-8904 in the posterior syncytium in association with the bacterial mass throughout embryogenesis suggests that both genes are involved additionally in bacteriocyte development and possibly play roles in regulating the transovarial transmission of Buchnera.
Amino acid transporters AAAP-536 and APC-8904 do more than move amino acids across membranes
A handful of Drosophila orthologs of AAAP-536 and APC-8904 have been functionally characterized. On the basis of our observations of AAAP-536 and APC-8904 expression and the functional characterization of Drosophila orthologs, we hypothesize that AAAP-536 and APC-8904 function not to move amino acids across membranes but rather serve as effectors of the endocytotic transfer of Buchnera from maternal bacteriocytes to embryos and as nutrient sensors. AAAP-536, the most highly expressed bacteriocyte-expressed amino acid transporter, is orthologous to Drosophila CG8785 , a gene whose functionally characterized paralogs include pathetic  and polyphemus (CG12943) . Encoding a proton-assisted transporter (PAT) related to the mammalian solute carrier 36 (SLC36) amino acid transporters, pathetic is a well-known mediator of cell growth [24, 27, 31, 56]. Localized to the cell surface and late endosomes and lysosomes (LELs) of Drosophila fat body, pathetic senses extracellular amino acids and activates the target of rapamycin (TOR) signaling pathway [24, 29, 56]. The other functionally characterized CG8785 paralog, polyphemus, encodes a glutamate transporter that functions to maintain immune-related function by modulating phagocytosis of microbial pathogens and controlling bacterial growth during infection . The second most highly expressed bacteriocyte-expressed amino acid transporter, APC-8904, is orthologous to Drosophila slimfast [21, 34] and closely related to the slimfast paralogs minidiscs  and genderblind . Drosophila slimfast and minidiscs functions are associated with nutrition and tissue growth involving the TOR signaling pathway [20, 21, 30, 58]. Encoding a glutamate transporter, genderblind functions in the nervous system with mutations in genderblind yielding a homosexual behavior phenotype [25, 57]. In addition to controlling glutamatergic synapse strength, genderblind functions with polyphemus to affect phagocytosis of microbial pathogens . Significantly, taken together, functionally characterized Drosophila transporters related to AAAP-536 and APC-8904 are nutrient sensors that function as upstream mediators that activate the TOR signaling pathway . Remarkably, activation of TOR signaling stimulates bulk endocytic uptake via clathrin-mediated endocytosis [26, 56]. In light of this genetic evidence from Drosophila and the expression of AAAP-536 and APC-8904 in the posterior syncytium, we hypothesize that AAAP-536 and APC-8904 gene products coordinate with the TOR signaling pathway to acquire Buchnera endosymbionts through an endocytic process that results in reacquisition of the symbiosomal membrane at the time of blastula invasion.
An evolutionary perspective on bacteriocyte formation in insects
In contrast to the little that is known at a molecular level about the developmental specification of host bacteriocytes, other aspects of host/symbiont biology are becoming clear. It is clear that host/symbiont genome coevolution involves host/symbiont metabolic collaboration within metabolic pathways for amino acid [1, 6, 7, 54, 59–61] and vitamin biosynthesis . It is also clear that while host/symbiont genome coevolution is constrained by host genome gene content [63–66], it can also involve the acquisition of new genes by lateral gene transfer [55, 63–65, 67, 68], and duplication of existing genes [34, 36–38, 69, 70]. While characterizing patterns of host/symbiont genome coevolution has been facilitated greatly by the recent accumulation of genome data from host/symbiont pairs, it has also been facilitated by studying bacteriocyte transcriptomes [36, 54, 59, 63, 71]. Bacteriocyte cells that house heritable endosymbionts are common among sap-feeding insects . Substantial insights are likely to be gained by systematically studying similarities and differences in the genes and gene families expressed in bacteriocytes within a comparative framework. Further, which gene sets are conserved in their expression across a diversity of sap-feeding insects through development will likely facilitate determination of the functions that modulate development of the bacteriocyte.
The molecular mechanisms of bacteriocyte development were recently documented in the seed bug, Nysius plebeius . Bacteriocyte cells of N. plebeius differentiate under regulation of homeobox transcription factors that include Ultrabithorax (Ubx), abdominal-A (abd-A), and Antennapedia (Antp). In gene knockdown experiments, suppression of Ubx results in the disappearance of bacteriocyte cells and failure of symbionts to localize, while suppression of abd-A alters the spatial organization of bacteriocyte cells, and suppression of Antp results in mislocalization of bacteriomes . Notably, homeobox transcription factors that include Ubx/Abd-A, Distal-less (Dll), and Engrailed (En) are also expressed in aphid bacteriocytes . While it is known that Dll is involved in aphid bacteriocyte determination , the way that Dll and other homeobox transcription factors determine bacteriocyte formation and bacteriome organization in aphids remains an open question.
We are motivated to understand the metabolic and developmental integration of hosts and symbionts in order to understand host/symbiont coevolution. We selected the two amino acid transporters that are most highly expressed in A. pisum bacteriocytes and applied whole-mount in situ hybridization to analyze their developmental expression patterns. While we conjecture that AAAP-536 and APC-8904 gene products are integral to bacteriocyte development, we think that the most intriguing implication of this work concerns our inference that AAAP-536 and APC-8904 gene products modulate TOR signaling. Interpreting our data with knowledge that blastula acquire Buchnera by endocytosis  leads us to venture that AAAP-536 and APC-8904 stimulate the endocytotic process by coordinating with the TOR signaling pathway to faithfully facilitate the vertical transmission of Buchnera from mother to embryo across tens of millions of years of evolution. Notable, within an evolutionary framework, is the fact that the orthologs of AAAP-536 and APC-8904 show the same patterns of bacteriocyte-enriched expression in the green peach aphid, Myzus persicae . Thus, we speculate that Mper-AAP11 and Mper-APC08  are similarly expressed during asexual embryogenesis in M. persicae. This study highlights that duplicated amino acid transporters have potentially been retained in the aphid genome to play novel developmental roles in symbiont recruitment and genesis of the aphid bacteriome.
amino acid/auxin permease family
- abd-A :
- Antp :
amino acid/polyamine/organocation family
- βTUB :
- EF1α :
elongation factor 1-alpha
- GAPDH :
late endosomes and lysosomes
open reading frame
proton-dependent amino acid transporter
polymerase chain reaction
- RPL7 :
ribosomal protein L7
- RPL32 :
ribosomal protein L32
mammalian solute carrier 36
target of rapamycin
- Ubx :
ACCW conceived of the experiments. ACCW, HLL, and CC designed the experiments. HLL performed the experiments. HLL, CC, and ACCW contributed to data analysis and preparation of figures for publication. HLL and ACCW drafted the manuscript. All authors contributed to preparation of the final manuscript. All authors read and approved the final manuscript.
This work benefitted from useful discussions with Daniel Price, Rebecca Duncan, Honglin Feng, Athula Wikramanayake, and Yi-min Hsiao. We thank James Baker of the Confocal Microscopy Core facility in the Department of Biology at the University of Miami, Yi-Chun Chuang of Technology Commons, College of Life Science, National Taiwan University, together with Lingyu Wang and Wei Wu for technical support. Financial support for the project came from Ministry of Science and Technology Awards 102-2917-I-002-006, 104-2811-B-002-134 and 105-2811-B-002-135 (to HLL), and 104-2313-B-002-022-MY3 (to CC and to HLL) and National Science Foundation award IOS-1121847 and IOS-1354154 (to ACCW).
The confocal and in situ localization work reported here was completed by HLL while she was a fellow of the Taiwanese Ministry of Science and Technology Graduate Student Study Abroad Program in the laboratory of ACCW at the University of Miami. HLL completed her Ph.D. in the laboratory of CC studying developmental plasticity of germline development in the pea aphid. CC is a developmental biologist who is studying oogenesis and embryogenesis in aphids and other insects. ACCW is an evolutionary biologist who is studying the coevolution of plant sap-feeding insects and their obligate intracellular bacterial symbionts. Recent work in Wilson’s group focuses on functional characterization of amino acid transporters in sternorrhyncha insects.
The authors declare that they have no competing interests.
Availability of data and materials
Materials described in the manuscript, including all relevant raw data, will be freely available to any scientist wishing to use them for non-commercial purposes.
This work was supported by the following Grants from the Ministry of Science and Technology (102-2917-I-002-006, 104-2811-B-002-134, 105-2811-B-002-135 and 104-2313-B-002-022-MY3) and the National Science Foundation (IOS-1121847 and IOS-1354154). The funding bodies had no role in the design of the study, or the collection, analysis, and interpretation of data.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Shigenobu S, Wilson ACC. Genomic revelations of a mutualism: the pea aphid and its obligate bacterial symbiont. Cell Mol Life Sci. 2011;68:1297–309.View ArticlePubMedPubMed CentralGoogle Scholar
- Lu HL, Price DRG, Wikramanayake A, Chang CC, Wilson ACC. Ontogenetic differences in localization of glutamine transporter ApGLNT1 in the pea aphid demonstrate that mechanisms of host/symbiont integration are not similar in the maternal versus embryonic bacteriome. EvoDevo. 2016;7:1.View ArticlePubMedPubMed CentralGoogle Scholar
- Douglas AE. Phloem-sap feeding by animals: problems and solutions. J Exp Bot. 2006;57:747–54.View ArticlePubMedGoogle Scholar
- Karley AJ, Douglas AE, Parker WE. Amino acid composition and nutritional quality of potato leaf phloem sap for aphids. J Exp Biol. 2002;205:3009–18.PubMedGoogle Scholar
- Sandstrom J, Pettersson J. Amino acid composition of phloem sap and the relation to intraspecific variation in pea aphid (Acyrthosiphon pisum) performance. J Insect Physiol. 1994;40:947–55.View ArticleGoogle Scholar
- Russell CW, Poliakov A, Haribal M, Jander G, van Wijk KJ, Douglas AE. Matching the supply of bacterial nutrients to the nutritional demand of the animal host. Proc Biol Sci. 2014;281:20141163.View ArticlePubMedPubMed CentralGoogle Scholar
- Wilson ACC, Ashton PD, Calevro F, Charles H, Colella S, Febvay G, et al. Genomic insight into the amino acid relations of the pea aphid, Acyrthosiphon pisum, with its symbiotic bacterium Buchnera aphidicola. Insect Mol Biol. 2010;19(Suppl 2):249–58.View ArticlePubMedGoogle Scholar
- Miura T, Braendle C, Shingleton A, Sisk G, Kambhampati S, Stern DL. A comparison of parthenogenetic and sexual embryogenesis of the pea aphid Acyrthosiphon pisum (Hemiptera: Aphidoidea). J Exp Zool B Mol Dev Evol. 2003;295:59–81.View ArticlePubMedGoogle Scholar
- Caillaud CM, Rahbe Y. Aposymbiosis in a cereal aphid: reproductive failure and influence on plant utilization. Ecol Entomol. 1999;24:111–4.View ArticleGoogle Scholar
- Douglas AE. Reproductive diapause and the bacterial symbiosis in the sycamore aphid Drepanosiphum platanoidis. Ecol Entomol. 2000;25:256–61.View ArticleGoogle Scholar
- Ishikawa H, Yamaji M. Symbionin, an aphid endosymbiont-specific protein-I production of insects deficient in symbiont. Insect Biochem. 1985;15:155–63.View ArticleGoogle Scholar
- Jayaraj S, Ehrhardt P, Schmutterer H. The effect of certain antibiotics on reproduction of the black bean aphid, Aphis fabae Scop. Ann Appl Biol. 1967;59:13–21.View ArticlePubMedGoogle Scholar
- Blackman RL. Reproduction, cytogenetics and development. In: Minks AK, Harrewijn P, editors. Aphids: their biology, natural enemies and control, vol. A. Amsterdam: Elsevier; 1987. p. 163–95.Google Scholar
- Will L. Entwicklungsgeschichte der Viviparen Aphiden. Jena: Gustav Fischer; 1888.Google Scholar
- Blackman RL. Early development of parthenogenetic egg in 3 species of aphids (Homoptera Aphididae). Int J Insect Morphol Embryol. 1978;7:33–44.View ArticleGoogle Scholar
- Koga R, Meng XY, Tsuchida T, Fukatsu T. Cellular mechanism for selective vertical transmission of an obligate insect symbiont at the bacteriocyte-embryo interface. Proc Natl Acad Sci U S A. 2012;109:E1230–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Braendle C, Miura T, Bickel R, Shingleton AW, Kambhampati S, Stern DL. Developmental origin and evolution of bacteriocytes in the aphid-Buchnera symbiosis. PLoS Biol. 2003;1:E21.View ArticlePubMedPubMed CentralGoogle Scholar
- Fukatsu T, Nikoh N, Kawai R, Koga R. The secondary endosymbiotic bacterium of the pea aphid Acyrthosiphon pisum (Insecta: homoptera). Appl Environ Microbiol. 2000;66:2748–58.View ArticlePubMedPubMed CentralGoogle Scholar
- Griffiths GW, Beck SD. Intracellular symbiotes of the pea aphid, Acyrthosiphon pisum. J Insect Physiol. 1973;19:75–84.View ArticleGoogle Scholar
- Hundal HS, Taylor PM. Amino acid transceptors: gate keepers of nutrient exchange and regulators of nutrient signaling. Am J Physiol Endocrinol Metab. 2009;296:E603–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Colombani J, Raisin S, Pantalacci S, Radimerski T, Montagne J, Leopold P. A nutrient sensor mechanism controls Drosophila growth. Cell. 2003;114:739–49.View ArticlePubMedGoogle Scholar
- Hyde R, Taylor PM, Hundal HS. Amino acid transporters: roles in amino acid sensing and signalling in animal cells. Biochem J. 2003;373:1–18.View ArticlePubMedPubMed CentralGoogle Scholar
- Taylor PM. Role of amino acid transporters in amino acid sensing. Am J Clin Nutr. 2014;99:223S–30S.View ArticlePubMedGoogle Scholar
- Goberdhan DCI, Meredith D, Richard Boyd CA, Wilson C. PAT-related amino acid transporters regulate growth via a novel mechanism that does not require bulk transport of amino acids. Development. 2005;132:2365–75.View ArticlePubMedGoogle Scholar
- Grosjean Y, Grillet M, Augustin H, Ferveur JF, Featherstone DE. A glial amino-acid transporter controls synapse strength and courtship in Drosophila. Nat Neurosci. 2008;11:54–61.View ArticlePubMedGoogle Scholar
- Hennig KM, Colombani J, Neufeld TP. TOR coordinates bulk and targeted endocytosis in the Drosophila melanogaster fat body to regulate cell growth. J Cell Biol. 2006;173:963–74.View ArticlePubMedPubMed CentralGoogle Scholar
- Lin WY, Williams C, Yan C, Koledachkina T, Luedke K, Dalton J, et al. The SLC36 transporter Pathetic is required for extreme dendrite growth in Drosophila sensory neurons. Genes Dev. 2015;29:1120–35.View ArticlePubMedPubMed CentralGoogle Scholar
- Armstrong AR, Laws KM, Drummond-Barbosa D. Adipocyte amino acid sensing controls adult germline stem cell number via the amino acid response pathway and independently of target of rapamycin signaling in Drosophila. Development. 2014;141:4479–88.View ArticlePubMedPubMed CentralGoogle Scholar
- Heublein S, Kazi S, Ogmundsdottir MH, Attwood EV, Kala S, Richard Boyd CA, et al. Proton-assisted amino-acid transporters are conserved regulators of proliferation and amino-acid-dependent mTORC1 activation. Oncogene. 2010;29:4068–79.View ArticlePubMedPubMed CentralGoogle Scholar
- Martin JF, Hersperger E, Simcox A, Shearn A. minidiscs encodes a putative amino acid transporter subunit required non-autonomously for imaginal cell proliferation. Mech Dev. 2000;92:155–67.View ArticlePubMedGoogle Scholar
- Lin WY, Williams CR, Yan C, Parrish JZ. Functions of the SLC36 transporter Pathetic in growth control. Fly. 2015;9:99–106.View ArticlePubMedGoogle Scholar
- Boll M, Daniel H, Gasnier B. The SLC36 family: proton-coupled transporters for the absorption of selected amino acids from extracellular and intracellular proteolysis. Pflugers Arch. 2004;447:776–9.View ArticlePubMedGoogle Scholar
- Gonzalez EA, Garg A, Tang J, Nazario-Toole AE, Wu LP. A glutamate-dependent redox system in blood cells is integral for phagocytosis in Drosophila melanogaster. Curr Biol. 2013;23:2319–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Price DRG, Duncan RP, Shigenobu S, Wilson ACC. Genome expansion and differential expression of amino acid transporters at the aphid/Buchnera symbiotic interface. Mol Biol Evol. 2011;28:3113–26.View ArticlePubMedGoogle Scholar
- Price DRG, Feng H, Baker JD, Bavan S, Luetje CW, Wilson ACC. Aphid amino acid transporter regulates glutamine supply to intracellular bacterial symbionts. Proc Natl Acad Sci U S A. 2014;111:320–5.View ArticlePubMedGoogle Scholar
- Duncan RP, Feng H, Nguyen DM, Wilson ACC. Gene family expansions in aphids maintained by endosymbiotic and non-symbiotic traits. Genome Biol Evol. 2016;8:753–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Dahan RA, Duncan RP, Wilson ACC, Davalos LM. Amino acid transporter expansions associated with the evolution of obligate endosymbiosis in sap-feeding insects (Hemiptera: sternorrhyncha). BMC Evol Biol. 2015;15:52.View ArticlePubMedPubMed CentralGoogle Scholar
- Duncan RP, Husnik F, Van Leuven JT, Gilbert DG, Davalos LM, McCutcheon JP, Wilson ACC. Dynamic recruitment of amino acid transporters to the insect/symbiont interface. Mol Ecol. 2014;23:1608–23.View ArticlePubMedGoogle Scholar
- AphidBase: an aphid genomics database. INRA, Rennes. 2016. http://bipaa.genouest.org/is/aphidbase/. Accessed 28 Oct 2016.
- The International Aphid Genomics Consortium. Genome sequence of the pea aphid Acyrthosiphon pisum. PLoS Biol. 2010;8:e1000313.View ArticlePubMed CentralGoogle Scholar
- Chang CC, Lee WC, Cook CE, Lin GW, Chang T. Germ-plasm specification and germline development in the parthenogenetic pea aphid Acyrthosiphon pisum: Vasa and Nanos as markers. Int J Dev Biol. 2006;50:413–21.View ArticlePubMedGoogle Scholar
- Chang CC, Lin GW, Cook CE, Horng SB, Lee HJ, Huang TY. Apvasa marks germ-cell migration in the parthenogenetic pea aphid Acyrthosiphon pisum (Hemiptera: Aphidoidea). Dev Genes Evol. 2007;217:275–87.View ArticlePubMedGoogle Scholar
- Chang CC, Huang TY, Shih CL, Lin GW, Chang TP, Chiu H, Chang WC. Whole-mount identification of gene transcripts in aphids: protocols and evaluation of probe accessibility. Arch Insect Biochem Physiol. 2008;68:186–96.View ArticlePubMedGoogle Scholar
- Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J. qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007;8:R19.View ArticlePubMedPubMed CentralGoogle Scholar
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3:RESEARCH0034.
- Andersen CL, Jensen JL, Orntoft TF. Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004;64:5245–50.View ArticlePubMedGoogle Scholar
- Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper–Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004;26:509–15.View ArticlePubMedGoogle Scholar
- Ishikawa A, Ogawa K, Gotoh H, Walsh TK, Tagu D, Brisson JA, et al. Juvenile hormone titre and related gene expression during the change of reproductive modes in the pea aphid. Insect Mol Biol. 2012;21:49–60.View ArticlePubMedGoogle Scholar
- Whyard S, Singh AD, Wong S. Ingested double-stranded RNAs can act as species-specific insecticides. Insect Biochem Mol Biol. 2009;39:824–32.View ArticlePubMedGoogle Scholar
- Shakesby AJ, Wallace IS, Isaacs HV, Pritchard J, Roberts DM, Douglas AE. A water-specific aquaporin involved in aphid osmoregulation. Insect Biochem Mol Biol. 2009;39:1–10.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25:402–8.View ArticlePubMedGoogle Scholar
- Chang CC, Hsiao YM, Huang TY, Cook CE, Shigenobu S, Chang TH. Noncanonical expression of caudal during early embryogenesis in the pea aphid Acyrthosiphon pisum: maternal cad-driven posterior development is not conserved. Insect Mol Biol. 2013;22:442–55.View ArticlePubMedGoogle Scholar
- Huang TY, Cook CE, Davis GK, Shigenobu S, Chen RPY, Chang CC. Anterior development in the parthenogenetic and viviparous form of the pea aphid, Acyrthosiphon pisum: hunchback and orthodenticle expression. Insect Mol Biol. 2010;19(Suppl 2):75–85.View ArticlePubMedGoogle Scholar
- Hansen AK, Degnan PH. Widespread expression of conserved small RNAs in small symbiont genomes. ISME J. 2014;8:2490–502.View ArticlePubMedPubMed CentralGoogle Scholar
- Nakabachi A, Ishida K, Hongoh Y, Ohkuma M, Miyagishima SY. Aphid gene of bacterial origin encodes a protein transported to an obligate endosymbiont. Curr Biol. 2014;24:R640–1.View ArticlePubMedGoogle Scholar
- Ogmundsdottir MH, Heublein S, Kazi S, Reynolds B, Visvalingam SM, Shaw MK, Goberdhan DCI. Proton-assisted amino acid transporter PAT1 complexes with Rag GTPases and activates TORC1 on late endosomal and lysosomal membranes. PLoS ONE. 2012;7:e36616.View ArticlePubMedPubMed CentralGoogle Scholar
- Augustin H, Grosjean Y, Chen K, Sheng Q, Featherstone DE. Nonvesicular release of glutamate by glial xCT transporters suppresses glutamate receptor clustering in vivo. J Neurosci. 2007;27:111–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Reynolds B, Roversi P, Laynes R, Kazi S, Richard Boyd CA, Goberdhan DCI. Drosophila expresses a CD98 transporter with an evolutionarily conserved structure and amino acid-transport properties. Biochem J. 2009;420:363–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Hansen AK, Moran NA. Aphid genome expression reveals host-symbiont cooperation in the production of amino acids. Proc Natl Acad Sci U S A. 2011;108:2849–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Rabatel A, Febvay G, Gaget K, Duport G, Baa-Puyoulet P, Sapountzis P, et al. Tyrosine pathway regulation is host-mediated in the pea aphid symbiosis during late embryonic and early larval development. BMC Genom. 2013;14:235.View ArticleGoogle Scholar
- Russell CW, Bouvaine S, Newell PD, Douglas AE. Shared metabolic pathways in a coevolved insect-bacterial symbiosis. Appl Environ Microbiol. 2013;79:6117–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Price DRG, Wilson ACC. A substrate ambiguous enzyme facilitates genome reduction in an intracellular symbiont. BMC Biol. 2014;12:110.View ArticlePubMedPubMed CentralGoogle Scholar
- Husnik F, Nikoh N, Koga R, Ross L, Duncan RP, Fujie M, et al. Horizontal gene transfer from diverse bacteria to an insect genome enables a tripartite nested mealybug symbiosis. Cell. 2013;153:1567–78.View ArticlePubMedGoogle Scholar
- Moran NA, Jarvik T. Lateral transfer of genes from fungi underlies carotenoid production in aphids. Science. 2010;328:624–7.View ArticlePubMedGoogle Scholar
- Sloan DB, Nakabachi A, Richards S, Qu J, Murali SC, Gibbs RA, Moran NA. Parallel histories of horizontal gene transfer facilitated extreme reduction of endosymbiont genomes in sap-feeding insects. Mol Biol Evol. 2014;31:857–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Wilson ACC, Duncan RP. Signatures of host/symbiont genome coevolution in insect nutritional endosymbioses. Proc Natl Acad Sci U S A. 2015;112:10255–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Nikoh N, McCutcheon JP, Kudo T, Miyagishima SY, Moran NA, Nakabachi A. Bacterial genes in the aphid genome: absence of functional gene transfer from Buchnera to its host. PLoS Genet. 2010;6:e1000827.View ArticlePubMedPubMed CentralGoogle Scholar
- Nikoh N, Nakabachi A. Aphids acquired symbiotic genes via lateral gene transfer. BMC Biol. 2009;7:12.View ArticlePubMedPubMed CentralGoogle Scholar
- Duncan RP, Nathanson L, Wilson ACC. Novel male-biased expression in paralogs of the aphid slimfast nutrient amino acid transporter expansion. BMC Evol Biol. 2011;11:253.View ArticlePubMedPubMed CentralGoogle Scholar
- Shigenobu S, Stern DL. Aphids evolved novel secreted proteins for symbiosis with bacterial endosymbiont. Proc Biol Sci. 2013;280:20121952.View ArticlePubMedPubMed CentralGoogle Scholar
- Nakabachi A, Shigenobu S, Sakazume N, Shiraki T, Hayashizaki Y, Carninci P, et al. Transcriptome analysis of the aphid bacteriocyte, the symbiotic host cell that harbors an endocellular mutualistic bacterium, Buchnera. Proc Natl Acad Sci U S A. 2005;102:5477–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Douglas AE. Mycetocyte symbiosis in insects. Biol Rev Camb Philos Soc. 1989;64:409–34.View ArticlePubMedGoogle Scholar
- Matsuura Y, Kikuchi Y, Miura T, Fukatsu T. Ultrabithorax is essential for bacteriocyte development. Proc Natl Acad Sci U S A. 2015;112:9376–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Chang CC, Huang TY, Cook CE, Lin GW, Shih CL, Chen RPY. Developmental expression of Apnanos during oogenesis and embryogenesis in the parthenogenetic pea aphid Acyrthosiphon pisum. Int J Dev Biol. 2009;53:169–76.View ArticlePubMedGoogle Scholar