Anterior development in the parthenogenetic and viviparous form of the pea aphid, Acyrthosiphon pisum: hunchback and orthodenticle expression

Germline specification and axis determination in viviparous and oviparous pea aphids: conserved and divergent features

Aphids are hemimetabolous insects that undergo incomplete metamorphosis without pupation. The annual life cycle of most aphids includes both an asexual (viviparous) and a sexual (oviparous) phase. Sexual reproduction only occurs once per year and is followed by many generations of asexual reproduction, during which aphids propagate exponentially with telescopic development. Here, we discuss the potential links between viviparous embryogenesis and derived developmental features in the pea aphid Acyrthosiphon pisum, particularly focusing on germline specification and axis determination, both of which are key events of early development in insects. We also discuss potential evolutionary paths through which both viviparous and oviparous females might have come to utilize maternal germ plasm to drive germline specification. This developmental strategy, as defined by germline markers, has not been reported in other hemimetabolous insects. In viviparous females, furthermore, we discuss whether molecules that in other insects characterize germ plasm, like Vasa, also participate in posterior determination and how the anterior localization of the hunchback orthologue Ap-hb establishes the anterior-posterior axis. We propose that the linked chain of developing oocytes and embryos within each ovariole and the special morphology of early embryos might have driven the formation of evolutionary novelties in germline specification and axis determination in the viviparous aphids. Moreover, based upon the finding that the endosymbiont Buchnera aphidicola is closely associated with germ cells throughout embryogenesis, we propose presumptive roles for B. aphidicola in aphid development, discussing how it might regulate germline migration in both reproductive modes of pea aphids. In summary, we expect that this review will shed light on viviparous as well as oviparous development in aphids.

Reproductive constraint in the social aphid Ceratovacuna japonica: Sterility regulation in the soldier caste of a viviparous insect

Differentiation of the non-reproductive caste is a unique feature of eusocial insects. Apoptosis in oocytes plays a major role in constraining the reproductivity of the eusocial insects including bees, ants, and termites. However, the regulation of reproductive constraint in non-reproductives of primitively eusocial insects other than hymenopterans and blattodeans is almost unknown. Here, we investigated the soldier sterility in a hemipteran insect, the social aphid Ceratovacuna japonica. We compared the gonads of soldiers, that are completely sterile, with those of reproductives in their viviparous development. We found that soldiers possess a pair of ovaries and the same number of germaria as reproductives, but soldiers' ovarioles were small and lacking gastrulating embryos. Unlike in most model social insects, the staining of cleaved Caspase-3 showed apoptosis in the maternal nutritive cells, rather in the oocyte, of soldier ovaries. In addition, the ubiquitous C. japonica vasa1 and piwi2a expression indicates the developmental failure of embryos in soldier ovaries. The absence of posterior nos1, an insect posterior determinant, indicates deficient posterior patterning in soldier ovarioles. Our findings suggest a different mode of reproductive constraint, which regulates both oogenesis and embryogenesis in a viviparous insect ovary. This is the first report of the reproductive constraint in a viviparous social insect at the molecular level.

Perspectives on divergence of early developmental regulatory pathways: Insight from the evolution of echinoderm double negative gate

Evolution of gene regulatory networks (GRN) that orchestrate the highly coordinated course of development, is made possible by the network's robust nature for incorporating change without detrimental developmental outcome. It can be considered that the upstream network regulating early development, has immense influence over succeeding pathways thus may be less subjected to evolutionary modification. However, recent studies show incorporation of novel genes in such early developmental pathways such as the echinoderm pmar1 as evidence for drastic change occurring high in the GRN hierarchy. Here we discuss the mechanisms that underlie divergence of early developmental pathways utilizing promising insights from the evolution of echinoderm early mesoderm specification pathway of Pmar1-HesC double negative gate found solely in the euechinoid sea urchin lineage, as well as examples from other groups such as Spiralia and Drosophila.

The evolution of the gene regulatory networks patterning the Drosophila Blastoderm

The Drosophila blastoderm gene regulatory network is one of the best studied networks in biology. It is composed of a series of tiered sub-networks that act sequentially to generate a primary segmental pattern. Many of these sub-networks have been studied in other arthropods, allowing us to reconstruct how each of them evolved over the transition from the arthropod ancestor to the situation seen in Drosophila today. I trace the evolution of each of these networks, showing how some of them have been modified significantly in Drosophila relative to the ancestral state while others are largely conserved across evolutionary timescales. I compare the putative ancestral arthropod segmentation network with that found in Drosophila and discuss how and why it has been modified throughout evolution, and to what extent this modification is unusual.

Germline expression of the hunchback orthologues in the asexual viviparous aphids: a conserved feature within the Aphididae

In animals, differentiation of germline from soma usually takes place during embryogenesis. Genes and their products that are preferentially expressed in the embryonic germ cells are regarded as candidates for maintaining germline fate or promoting germline identity. In Drosophila, for example, the protein encoded by the germline gene vasa is specifically restricted to the germ cells, while products of the gap gene hunchback (hb), a somatic gene, are preferentially expressed in the neuroblasts. In this study, we report the expression of both messenger RNA and protein encoded by Aphb, an hb orthologue in the asexual viviparous pea aphid Acyrthosiphon pisum, in germ cells as well as in neuroblasts. We infer that expression of Aphb messenger RNA in the germ cells during the formation of germaria is required for the anterior localization of Aphb in the protruding oocytes. Germarial expression and anterior localization of ApKrüppel was also identified but, unlike Aphb, its expression was not detected in the migrating germ cells. Very similar patterns of hb expression were also identified in the green peach aphid Myzus persicae, suggesting that germline expression of hb is conserved within the Aphididae. To date, this pattern of hb germline expression has not been reported in other insects.

A feed-forward relay integrates the regulatory activities of Bicoid and Orthodenticle via sequential binding to suboptimal sites

Datta et al. define three major classes of enhancers that are differentially sensitive to binding and transcriptional activation by Bicoid (Bcd) and Orthodenticle (Otd). The specific activities of enhancers in each class are mediated by DNA motif variants preferentially bound by Bcd or Otd and the presence or absence of sites for cofactors that interact with these proteins.

The K50 (lysine at amino acid position 50) homeodomain (HD) protein Orthodenticle (Otd) is critical for anterior patterning and brain and eye development in most metazoans. In Drosophila melanogaster, another K50HD protein, Bicoid (Bcd), has evolved to replace Otd's ancestral function in embryo patterning. Bcd is distributed as a long-range maternal gradient and activates transcription of a large number of target genes, including otd. Otd and Bcd bind similar DNA sequences in vitro, but how their transcriptional activities are integrated to pattern anterior regions of the embryo is unknown. Here we define three major classes of enhancers that are differentially sensitive to binding and transcriptional activation by Bcd and Otd. Class 1 enhancers are initially activated by Bcd, and activation is transferred to Otd via a feed-forward relay (FFR) that involves sequential binding of the two proteins to the same DNA motif. Class 2 enhancers are activated by Bcd and maintained by an Otd-independent mechanism. Class 3 enhancers are never bound by Bcd, but Otd binds and activates them in a second wave of zygotic transcription. The specific activities of enhancers in each class are mediated by DNA motif variants preferentially bound by Bcd or Otd and the presence or absence of sites for cofactors that interact with these proteins. Our results define specific patterning roles for Bcd and Otd and provide mechanisms for coordinating the precise timing of gene expression patterns during embryonic development.

Amino acid transporters implicated in endocytosis of Buchnera during symbiont transmission in the pea aphid

Background

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.

Results

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.

Conclusions

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-016-0061-7) contains supplementary material, which is available to authorized users.

Diversity of molecules and mechanisms in establishing insect anterior-posterior polarity

Anterior-posterior (AP) patterning is an essential process that requires the generation of large amounts of positional information to properly specify many distinct cell fates along the long axis of the insect embryo. While the general molecular basis of this process has long been known in the fly Drosophila, detailed understanding of this process is still emerging in other insect species. What is now clear is that this process in extremely labile, and distinct AP patterning programs can exist even within a single species. This review presents recent progress on this topic in an attempt to synthesize the disparate data and provide an outlook on the future of the field.

Reliable protocols for whole-mount fluorescent in situ hybridization (FISH) in the pea aphid Acyrthosiphon pisum: a comprehensive survey and analysis

RNA in situ hybridization (ISH), including chromogenic ISH (CISH) and fluorescent ISH (FISH), has become a powerful tool for revealing the spatial distribution of gene transcripts in model organisms. Previously, we developed a robust protocol for whole-mount RNA CISH in the pea aphid Acyrthosiphon pisum, an emerging insect genomic model. In order to improve the resolving capacity of gene detection, we comprehensively surveyed current protocols of whole-mount RNA-FISH and developed protocols that allow, using confocal microscopy, clearer visualization of target messenger RNAs (mRNAs) - including those subcellularly localized and those with spatially overlapping expression. We find that Fast dye-based substrate fluorescence (SF), tyramide signal amplification (TSA), and TSA Plus all enable identifying gene expression thanks to multiplex amplification of fluorescent signals. By contrast, methods of direct fluorescence (DF) do not allow visualizing signals. Detection of a single gene target was achieved with SF and TSA Plus for most mRNAs, whereas TSA only allowed visualization of abundant transcripts such as Apvas1 and Appiwi2 in the germ cells. For detection of multiple gene targets using double FISH, we recommend: (i) TSA/TSA, rather than TSA Plus/TSA Plus for colocalized mRNAs abundantly expressed in germ cells, as proteinase K treatment can be omitted; and (ii) SF/TSA Plus for other gene targets such as Apen1 and Apen2 as inactivation of enzyme conjugates is not required. SF/SF is not ideal for double FISH experiments due to signal blurring. Based on these new conditions for RNA-FISH, we have obtained a better understanding of germline specification and embryonic segmentation in the pea aphid. We anticipate that the RNA-FISH protocols for the pea aphid may also be used for other aphids and possibly other insect species, thus expanding the range of species from which useful insights into development and evolution may be obtained.

Plant-mediated RNAi of a gap gene-enhanced tobacco tolerance against the Myzus persicae

Plant-mediated RNAi has been developed as a powerful weapon in the fight against agricultural insect pests. The gap gene hunchback (hb) is of crucial importance in insect axial patterning and knockdown of hb is deforming and lethal to the next generation. The peach potato aphid, Myzus persicae (Sulzer), has many host plants and can be found throughout the world. To investigate the effect of plant-mediated RNAi on control of this insect, the hb gene in M. persicae was cloned, plant RNAi vector was constructed, and transgenic tobacco expressing Mphb dsRNA was developed. Transgenic tobacco had a different integration pattern of the transgene. Bioassays were performed by applying neonate aphids to homozygous transgenic plants in the T2 generation. Results revealed that continuous feeding of transgenic diet reduced Mphb mRNA level in the fed aphids and inhibited insect reproduction, indicating successful knockdown of the target gene in M. persicae by plant-mediated RNAi.

Hunchback is required for abdominal identity suppression and germband growth in the parthenogenetic embryogenesis of the pea aphid, Acyrthosiphon pisum

Aphid, a short germband insect, displays an embryogenesis different from that of long germband insect species. Furthermore, the development of its parthenogenetic and viviparous embryo is different from that of the embryo resulting from sexual reproduction. To better understand the genetic regulation of this type of embryogenesis, the functions of hunchback in asexual Acyrthosiphon pisum were investigated by parental RNAi. Microinjection of Aphb double-stranded RNA yielded several defective phenotypes. Quantitative real-time PCR analysis revealed that these defects resulted from reduction of Aphb mRNA level in injected aphids. All these results suggested that the hb gene in parthenogenetic and viviparous Acyrthosiphon pisum was involved in abdominal identity suppression and germband growth as its homologue does in sexual insects.

Noncanonical expression of caudal during early embryogenesis in the pea aphid Acyrthosiphon pisum: maternal cad-driven posterior development is not conserved

Previously we identified anterior localization of hunchback (Aphb) mRNA in oocytes and early embryos of the parthenogenetic and viviparous pea aphid Acyrthosiphon pisum, suggesting that the breaking of anterior asymmetry in the oocytes leads to the formation of the anterior axis in embryos. In order to study posterior development in the asexual pea aphid, we cloned and analysed the developmental expression of caudal (Apcad), a posterior gene highly conserved in many animal phyla. We found that transcripts of Apcad were not detected in germaria, oocytes and embryos prior to the formation of the blastoderm in the asexual (viviparous) pea aphid. This unusual expression pattern differs from that of the existing insect models, including long- and short-germ insects, where maternal cad mRNA is passed to the early embryos and forms a posterior-anterior gradient. The first detectable Apcad expression occurred in the newly formed primordial germ cells and their adjacent blastodermal cells during late blastulation. From gastrulation onward, and as in other insects, Apcad mRNA is restricted to the posteriormost region of the germ band. Similarly, in the sexual (oviparous) oocytes we were able to identify anterior localization of Aphb mRNA but posterior localization of Apcad was not detected. This suggests that cad-driven posterior development is not conserved during early embryogenesis in asexual and sexual pea aphids.

The pea aphid (Acyrthosiphon pisum) genome encodes two divergent early developmental programs

The pea aphid (Acyrthosiphon pisum) can reproduce either sexually or asexually (parthenogenetically), giving rise, in each case, to almost identical adults. These two modes of reproduction are accompanied by differences in ovarian morphology and the developmental environment of the offspring, with sexual forms producing eggs that are laid, whereas asexual development occurs within the mother. Here we examine the effect each mode of reproduction has on the expression of key maternal and axis patterning genes; orthodenticle (otd), hunchback (hb), caudal (cad) and nanos (nos). We show that three of these genes (Ap-hb, Ap-otd and Ap-cad) are expressed differently between the sexually and asexually produced oocytes and embryos of the pea aphid. We also show, using immunohistochemistry and cytoskeletal inhibitors, that Ap-hb RNA is localized differently between sexually and asexually produced oocytes, and that this is likely due to differences in the 3' untranslated regions of the RNA. Furthermore, Ap-hb and Ap-otd have extensive expression domains in early sexually produced embryos, but are not expressed at equivalent stages in asexually produced embryos. These differences in expression likely correspond with substantial changes in the gene regulatory networks controlling early development in the pea aphid. These data imply that in the evolution of parthenogenesis a new program has evolved to control the development of asexually produced embryos, whilst retaining the existing, sexual, developmental program. The patterns of modification of these developmental processes mirror the changes that we see in developmental processes between species, in that early acting pathways in development are less constrained, and evolve faster, than later ones. We suggest that the evolution of the novel asexual development pathway in aphids is not a simple modification of an ancestral system, but the evolution of two very different developmental mechanisms occurring within a single species.

Canonical terminal patterning is an evolutionary novelty

Patterning of the terminal regions of the Drosophila embryo is achieved by an exquisitely regulated signal that passes between the follicle cells of the ovary, and the developing embryo. This pathway, however, is missing or modified in other insects. Here we trace the evolution of this pathway by examining the origins and expression of its components. The three core components of this pathway: trunk, torso and torso-like have different evolutionary histories and have been assembled step-wise to form the canonical terminal patterning pathway of Drosophila and Tribolium. Trunk, torso and a gene unrelated to terminal patterning, prothoraciotrophic hormone (PTTH), show an intimately linked evolutionary history, with every holometabolous insect, except the honeybee, possessing both PTTH and torso genes. Trunk is more restricted in its phylogenetic distribution, present only in the Diptera and Tribolium and, surprisingly, in the chelicerate Ixodes scapularis, raising the possibility that trunk and torso evolved earlier than previously thought. In Drosophila torso-like restricts the activation of the terminal patterning pathway to the poles of the embryo. Torso-like evolved in the pan-crustacean lineage, but based on expression of components of the canonical terminal patterning system in the hemimetabolous insect Acyrthosiphon pisum and the holometabolous insect Apis mellifera, we find that the canonical terminal-patterning system is not active in these insects. We therefore propose that the ancestral function of torso-like is unrelated to terminal patterning and that torso-like has become co-opted into terminal patterning in the lineage leading to Coleoptera and Diptera. We also show that this co-option has not resulted in changes to the molecular function of this protein. Torso-like from the pea aphid, honeybee and Drosophila, despite being expressed in different patterns, are functionally equivalent. We propose that co-option of torso-like into restricting the activity of trunk and torso facilitated the final step in the evolution of this pathway; the capture of transcriptional control of target genes such as tailless and huckebein by this complex and novel patterning pathway.

Expansion of genes encoding piRNA-associated argonaute proteins in the pea aphid: diversification of expression profiles in different plastic morphs

Piwi-interacting RNAs (piRNAs) are known to regulate transposon activity in germ cells of several animal models that propagate sexually. However, the role of piRNAs during asexual reproduction remains almost unknown. Aphids that can alternate sexual and asexual reproduction cycles in response to seasonal changes of photoperiod provide a unique opportunity to study piRNAs and the piRNA pathway in both reproductive modes. Taking advantage of the recently sequenced genome of the pea aphid Acyrthosiphon pisum, we found an unusually large lineage-specific expansion of genes encoding the Piwi sub-clade of Argonaute proteins. In situ hybridisation showed differential expressions between the duplicated piwi copies: while Api-piwi2 and Api-piwi6 are "specialised" in germ cells their most closely related copy, respectively Api-piwi5 and Api-piwi3, are expressed in the somatic cells. The differential expression was also identified in duplicated ago3: Api-ago3a in germ cells and Api-ago3b in somatic cells. Moreover, analyses of expression profiles of the expanded piwi and ago3 genes by semi-quantitative RT-PCR showed that expressions varied according to the reproductive types. These specific expression patterns suggest that expanded aphid piwi and ago3 genes have distinct roles in asexual and sexual reproduction.

Diversity in insect axis formation: two orthodenticle genes and hunchback act in anterior patterning and influence dorsoventral organization in the honeybee (Apis mellifera)

Axis formation is a key step in development, but studies indicate that genes involved in insect axis formation are relatively fast evolving. Orthodenticle genes have conserved roles, often with hunchback, in maternal anterior patterning in several insect species. We show that two orthodenticle genes, otd1 and otd2, and hunchback act as maternal anterior patterning genes in the honeybee (Apis mellifera) but, unlike other insects, act to pattern the majority of the anteroposterior axis. These genes regulate the expression domains of anterior, central and posterior gap genes and may directly regulate the anterior gap gene giant. We show otd1 and hunchback also influence dorsoventral patterning by regulating zerknült (zen) as they do in Tribolium, but that zen does not regulate the expression of honeybee gap genes. This suggests that interactions between anteroposterior and dorsal-ventral patterning are ancestral in holometabolous insects. Honeybee axis formation, and the function of the conserved anterior patterning gene orthodenticle, displays unique characters that indicate that, even when conserved genes pattern the axis, their regulatory interactions differ within orders of insects, consistent with relatively fast evolution in axis formation pathways.

Evolution of a genomic regulatory domain: the role of gene co-option and gene duplication in the Enhancer of split complex

The Drosophila Enhancer of split complex [E(spl)-C] is a remarkable complex of genes many of which are effectors or modulators of Notch signaling. The complex contains different classes of genes including four bearded genes and seven basic helix-loop-helix (bHLH) genes. We examined the evolution of this unusual complex by identifying bearded and bHLH genes in the genome sequences of Arthropods. We find that a four-gene E(spl)-C, containing three bHLH genes and one bearded gene, is an ancient component of the genomes of Crustacea and Insects. The complex is well conserved in insects but is highly modified in Drosophila, where two of the ancestral genes of the complex are missing, and the remaining two have been duplicated multiple times. Through examining the expression of E(spl)-C genes in honeybees, aphids, and Drosophila, we determined that the complex ancestrally had a role in Notch signaling. The expression patterns of genes found inserted into the complex in some insects, or that of ancestral E(spl)-C genes that have moved out of the complex, imply that the E(spl)-C is a genomic domain regulated as a whole by Notch signaling. We hypothesize that the E(spl)-C is a Notch-regulated genomic domain conserved in Arthropod genomes for around 420 million years. We discuss the consequence of this conserved domain for the recruitment of novel genes into the Notch signaling cascade.

Genome sequence of the pea aphid Acyrthosiphon pisum

Aphids are important agricultural pests and also biological models for studies of insect-plant interactions, symbiosis, virus vectoring, and the developmental causes of extreme phenotypic plasticity. Here we present the 464 Mb draft genome assembly of the pea aphid Acyrthosiphon pisum. This first published whole genome sequence of a basal hemimetabolous insect provides an outgroup to the multiple published genomes of holometabolous insects. Pea aphids are host-plant specialists, they can reproduce both sexually and asexually, and they have coevolved with an obligate bacterial symbiont. Here we highlight findings from whole genome analysis that may be related to these unusual biological features. These findings include discovery of extensive gene duplication in more than 2000 gene families as well as loss of evolutionarily conserved genes. Gene family expansions relative to other published genomes include genes involved in chromatin modification, miRNA synthesis, and sugar transport. Gene losses include genes central to the IMD immune pathway, selenoprotein utilization, purine salvage, and the entire urea cycle. The pea aphid genome reveals that only a limited number of genes have been acquired from bacteria; thus the reduced gene count of Buchnera does not reflect gene transfer to the host genome. The inventory of metabolic genes in the pea aphid genome suggests that there is extensive metabolite exchange between the aphid and Buchnera, including sharing of amino acid biosynthesis between the aphid and Buchnera. The pea aphid genome provides a foundation for post-genomic studies of fundamental biological questions and applied agricultural problems.