Sexual and asexual oogenesis require the expression of unique and shared sets of genes in the insect Acyrthosiphon pisum
- Aurore Gallot1,
- Shuji Shigenobu2,
- Tomomi Hashiyama2,
- Stéphanie Jaubert-Possamai†1Email author and
- Denis Tagu†1Email author
© Gallot et al; licensee BioMed Central Ltd. 2012
Received: 21 September 2011
Accepted: 15 February 2012
Published: 15 February 2012
Although sexual reproduction is dominant within eukaryotes, asexual reproduction is widespread and has evolved independently as a derived trait in almost all major taxa. How asexuality evolved in sexual organisms is unclear. Aphids, such as Acyrthosiphon pisum, alternate between asexual and sexual reproductive means, as the production of parthenogenetic viviparous females or sexual oviparous females and males varies in response to seasonal photoperiodism. Consequently, sexual and asexual development in aphids can be analyzed simultaneously in genetically identical individuals.
We compared the transcriptomes of aphid embryos in the stages of development during which the trajectory of oogenesis is determined for producing sexual or asexual gametes. This study design aimed at identifying genes involved in the onset of the divergent mechanisms that result in the sexual or asexual phenotype. We detected 33 genes that were differentially transcribed in sexual and asexual embryos. Functional annotation by gene ontology (GO) showed a biological signature of oogenesis, cell cycle regulation, epigenetic regulation and RNA maturation. In situ hybridizations demonstrated that 16 of the differentially-transcribed genes were specifically expressed in germ cells and/or oocytes of asexual and/or sexual ovaries, and therefore may contribute to aphid oogenesis. We categorized these 16 genes by their transcription patterns in the two types of ovaries; they were: i) expressed during sexual and asexual oogenesis; ii) expressed during sexual and asexual oogenesis but with different localizations; or iii) expressed only during sexual or asexual oogenesis.
Our results show that asexual and sexual oogenesis in aphids share common genetic programs but diverge by adapting specificities in their respective gene expression profiles in germ cells and oocytes.
Sexual reproduction involves two main events: meiosis and fertilization, and creates new genotypes by shuffling allelic combinations. Although the predominance of sexual reproduction in eukaryotes supports this innovation as a successful reproduction strategy, asexuality has evolved independently multiple times from sexual ancestors in almost all major taxa [1–4], such as in stick insects  and Ranunculus plants . How asexuality has evolved in sexual organisms is unclear. In aphids, asexuality was acquired once about 250 million years ago by a common sexual ancestor . Most aphid species alternate between sexual reproduction and asexual parthenogenetic reproduction according to seasonal variations. In spring and summer, aphids reproduce asexually by parthenogenesis and produce clonal parthenogenetic female progeny by viviparity. The autumnal shortening of the photoperiod induces the concentration of juvenile hormone (JH) to decrease in the aphid haemolymph., and particular form of parthenogenetic female called the sexuparae are produced. Sexuparae females produce sexual females and males that subsequently mate to produce overwintering eggs. Although parthenogenetic viviparous females and sexual oviparous females exhibit major differences in morphology and behavior, they share the same genome. This phenomenon, called reproductive polyphenism, is an example of aphid phenotypic plasticity .
The cellular and cytogenetic bases of reproductive polyphenism have been described for several aphid species [7, 8]. In aphids, three generations are represented within one viviparous female: the mature embryos developing inside the maternal abdomen carry the first developmental stage of the third generation. This phenomenon is known as the 'telescoping of generations' . The embryonic developments of asexual and sexual females are similar until the formation of undifferentiated germaria, which takes place at the end of germ band retraction (stage 18) . The 32 germ cells within the germaria enter into early meiotic prophase (zygotene) and a partial association between bivalents occurs . After this stage, the embryonic development and maturation of germ cells diverge between sexual and asexual females. Ancestral sexual development occurs under short photoperiod: the basal cells of germaria progress through meiosis; the sexual oocyte displays chromosome pairing, synapsis and meiotic recombinations. Meiosis is blocked in first meiotic prophase until the egg is fertilized and laid. Under long photoperiod, sexual oogenesis is inhibited . The asexual oocyte ovulates in the posterior part of the germaria and forms a new follicle where it skips the first meiotic division and undergoes a single maturation division resulting in a discarded polar body and a diploid clonal oocyte. Synchronous mitotic divisions follow immediately without fertilization to initiate early embryogenesis, which continues within the parthenogenetic female abdomen and progresses through a series of 20 developmental stages .
We identified genes regulated in sexual and asexual aphid oogenesis by combining transcriptomic analyses and in situ hybridizations. In aphids, sexual and asexual reproduction can be analyzed in individuals possessing identical genomes. Sequencing of the LSR1 pea aphid genome  and the recent development of new genomic resources offer an unique opportunity to understand the molecular bases of sexual and asexual reproduction. Our results demonstrate that sexual and asexual reproduction pathways require unique and shared genetic programs.
Synchronization of sexual and asexual embryos
Differentially expressed transcripts
Description of the 33 differentially-expressed transcripts in sexual and asexual embryos
oocyte dorsal/ventral axis formation, egg activation
oocyte axis specification, mRNA polyadenylation, germ cell development
mitosis, oogenesis, ovarian follicle cell migration
female germline ring canal formation, karyosome formation, oogenesis
female meiosis chromosome segregation, germ cell development, oogenesis
female germline ring canal formation, karyosome formation, oogenesis
putative kelch-like repeat
female germline ring canal formation, karyosome formation, oogenesis
termination of RNA polymerase II transcription, female meiosis
RNA polyadenylation, histone mRNA catabolic process, mRNA processing
poly(A)+ mRNA export from nucleus
RNA metabolic process, nuclear-transcribed mRNA poly(A) tail shortening
gene silencing by miRNA, production of siRNA involved in RNA interference
chromatin assembly or disassembly, nucleosome assembly
chromatin assembly or disassembly, nucleosome assembly
histone methylation, regulation of gene expression, epigenetic
DNA repair, cell cycle, methylated histone residues binding
Cell cycle (3)
mitotic cell cycle, embryonic
cell division, exit from mitosis, condensed chromosome kinetochore
cytokinesis, intracellular distribution of mitochondria
cotranslational protein targeting to membrane
lipid particle organization, regulation of lipid storage
RNA dependant DNA Polymerase/Reverse Transcriptase
fat body development, gonad development, mesoderm development
No similarity (10)
Based on gene annotation in other organisms, 4 functional groups were defined (Table 1). Seven transcripts were related to oogenesis (e.g. nanos-1, orb, kelch); five transcripts were related to post-transcriptional regulation, such as mRNA polyadenylation (e.g. gld2) or mRNA silencing via RNA interference (cbp20); four transcripts were related to epigenetic regulation, such as histone modifications (suv4-20, uhrf1) or histone variants (H2B.3, H1); and three transcripts were assigned to cell cycle regulation (e.g. cyclin J, clasp1). Fourteen transcripts were not assigned to any functional category: among them, ten showed no similarity with characterized transcripts.
Our results displayed multiple genes related to oogesis (7/33). This oogenesis signature suggests that our experimental design, using artificially-induced asexual aphids, successfully targeted genes defining asexual or sexual oogenesis within the germaria or the first ovulated oocytes of the developing embryo.
Localization of regulated transcripts in sexual and asexual ovaries
We compared the localization of these differentially-expressed transcripts in the ovaries of untreated asexual and sexual aphids using whole-mount in situ hybridization (WISH), a technique that is not suitable for transcript quantification. Specific antisense riboprobes were produced for 25 transcripts: five produced no signal (Additional file 1), four showed ubiquitous distribution (Additional file 2) and 16 transcripts were specifically expressed in germaria (containing germ and nurse cells) and/or in oocytes. These results on untreated aphids were consistent with the microarray data from acetone- or kinoprene-treated aphids: most of the candidate genes are expressed in the ovarioles. These 16 transcripts were grouped based on their patterns of expression in the two types of ovaries.
By coupling quantitative transcriptomic and qualitative WISH approaches, we identified 33 transcripts associated with the onset of divergent transcripts expression leading to the establishment of the sexual or asexual phenotype in aphids. The involvement of these candidate genes in determining the developmental direction of oogenesis is supported by: i) a strong signature of oogenesis among the regulated transcripts (7/33); and ii) the oocyte- and/or germarial-specific localizations of 16 of the transcripts. These data illustrate both shared and unique regulatory patterns of transcript expression for the two modes of reproduction.
First, we identified transcript expression profiles shared by sexual and asexual oogenesis. These transcripts are expressed in sexual and asexual germaria and/or oocytes but display differences in their expression levels and/or in their localization. In various organisms, gld2, bicC, pop2 and orb are all involved in regulating the poly(A) tail length of maternal mRNAs. This process is a common mechanism of translational control essential for progression of meiosis and axis patterning during oocyte maturation [12–14]. The regulation of poly(A) tail length of maternal mRNAs results from a balance between concomitant deadenylation (translation repression) and polyadenylation (translation activation) . In Drosophila melanogaster, gld2, bicC and orb have been described as part of the cytoplasmic polyadenylation complex that polyadenylates mRNAs and activates their translation [13, 16]. Gld2-dependent mRNA polyadenylation and translational activation have been shown to be essential for oogenic meiosis . On the other hand, mRNA deadenylation by the major deadenylation complex in Drosophila, CCR4-NOT, which includes the deadenylases CCR4, POP2, and four NOT proteins , promotes negative regulation of target mRNAs. Deadenylation of specific maternal mRNAs could be involved in the precise temporal and spatial mRNA localization necessary for axial patterning during oogenesis . Expression of Api-gld2, Api-bicC, Api-pop2 and Api-orb in sexual and asexual germaria and/or in oocytes suggests that the balance between cytoplasmic polyadenylation and deadenylation in oocytes acts in the maturation of haploid as well as diploid oocytes in the pea aphid as it has been shown for other sexual organisms. Moreover, expression of the oocyte-specific linkers histoneH1 and uhrf1 in asexual and sexual oocytes suggests that epigenetic regulation is involved in pea aphid oogenesis. The linker histone H1 interacts with DNA to set-up and maintain the organization of the chromatin between the nucleosomes. Histone H1 proteins form a complex family of related proteins that have distinct species, tissue and developmental specificities and are involved in the epigenetic control of gene expression. The replacement of a somatic linker H1 by an oocyte specific one (H1oo) during oogenesis has been described in a wide range of species (reviewed in [18, 19]) and may promote a loosening in chromatin structure that appears to correlate with the initiation of meiosis . The presence of Api-histoneH1 in germaria and oocytes suggests that this oocyte-specific variant of H1 might be involved in both sexual and asexual oogenesis in aphid. uhrf1 encodes an ubiquitin-like protein that is involved in two epigenetic silencing pathways through its regulation of DNA and histone methylation [20, 21]. Although, uhrf1 has not been implicated in any oogenesis process to date, our results suggest a role for this protein in asexual and sexual oogenesis in aphids. Altogether, our results suggest that epigenetic regulatory mechanisms and the balance between cytoplasmic mRNA polyadenylation and deadenylation may have roles in regulating haploid and diploid oocyte maturation. Moreover, as orb transcript is the earliest to be differentially expressed during the developmental time course, we propose that Orb might play a major role in determining oocyte fate.
Second, we identified one gene, lsd1, that is specifically expressed in pea aphid sexual germaria and oocytes. The lsd1 transcript was not significantly detected either in asexual embryos according to microarray data (Additional file 3) or in asexual ovaries according to WISH. lsd1 encodes a protein involved in the regulation of lipid droplet storage in eukaryotic cells . Lipid droplets, like mature yolk bodies, are the major energy storage for oocytes and developing embryos in Drosophila. Our results suggest that lsd1 is required only for sexual oogenesis during which it may be involved in yolk accumulation. The absence of lsd1 transcripts in asexual embryos is consistent with the absence of yolk in parthenogenetic oocytes . Thus, identification of lsd1 expression as a specific feature of sexual oogenesis validates the accuracy of our approach. Despite the existence of genes known to be involved in meiotic recombination in the pea aphid genome , we identified no such genes as specific markers of sexual oocytes. However, germline expression of meiotic genes such as spo11 was shown to be similar in sexual and asexual aphids. Expression of this gene seems to be modulated between asexual and sexual aphids by alternative splicing .
Finally, we identified three transcripts that were specific to asexual oogenesis. No significant transcription signal has been observed for these genes either in sexual embryos using microarrays (Additional file 3) or in sexual ovaries using WISH. That the polypeptide encoded by ACYPI39770 does not bear any similarity to a protein or protein domain of any other organism, suggests that ACYPI39770 is either an orphan gene or a gene with an unusual accelerated evolution rate. In Drosophila, cyclin J and lodestar are involved in sexual oogenesis and early embryogenesis [26, 27]. In aphids, their expression profiles suggest that they may haveroles in asexual oogenesis only. The absence of lodestar and cyclin J transcripts from aphid sexual oocytes and embryos was particularly unexpected. Three lodestar paralogs were identified within the pea aphid genome (ACYPI000753, ACYPI000590 and ACYPI006971). Although the three copies were represented on the microarray by copy-specific oligonucleotides, only ACYPI000753 was differentially expressed in sexual and asexual embryos. Localization of ACYPI000753 transcripts by in situ hybridization showed that its expression was specific to asexual oocytes. However, although the sequence of the riboprobe used to locate ACYPI000753 transcripts was absent from the ACYPI006971 sequence, it shared 87% nucleotide identity with the unregulated paralog ACYPI000590. Therefore, cross hybridization between the ACYPI000753 riboprobe and ACYPI000590 transcripts cannot be excluded. Nevertheless, we detected no ACYPI000753 expression in sexual oogenesis. This unexpected result might reflect the functional specialization of at least one of the lodestar paralogs in asexual oogenesis.
Among the 16 transcripts found to be specifically localized in germaria and/or oocytes , lsd1 was the only one for which the expression level was found to be higher in sexual ovaries than in asexual ovaries according to the microarray results. This disproportion is consistent with the striking over-representation of transcripts that were upregulated in asexual embryos (29/33), as evidenced by the microarray analysis. Sexual oogenesis and asexual oogenesis display major inherent differences. Within embryos, the sexual oocyte nucleus remains blocked in meiosis prophase 1 whereas the asexual oocyte undergoes a modified but complete meiotic division immediately after ovulation . This modified meiosis involves specific changes within the oocyte; self-organized asters are formed, that then recruit centriole precursors and pericentriolar material . This process may be accompanied by the activation of microtubule stabilization factors or by the inactivation of microtubule destabilization factors specific to asexual oogenesis . The predominance of genes that are upregulated in asexual embryos may be partially explained by the production of molecular components that are specifically required for asexual oogenesis but not for sexual oogenesis.
In this study, we identified a number of putative embryonic target genes of the signaling cascade initiated by photoperiodic cues and transduced by JH. These genes are likely to be among those that direct sexual or asexual differentiation of the developing embryos. Our analysis of their expression patterns represents the first time the molecular bases of aphid sexual and asexual oogenesis have been addressed. Altogether, our results showed that aphid asexual oogenesis requires the same genes as sexual oogenesis but that these two distinctive differentiation processes display specific gene expression programs.
Aphid rearing, dissection and extraction
Clonal descendants from the sequenced Acyrthosiphon pisum LSR1 clone were reared on Vicia fabae at 18°C and at low density (5 individuals per plant) to prevent the production of winged morphs. Parthenogenesis was maintained under long photoperiod (16 h). In our conditions, the shortening of the photoperiod (12 h) requires three generations to induce the production of sexual morphs. The first asexual generation was reared under long photoperiod and transferred to short photoperiod (12 h) after the third instar moult. Adult asexual virginoparae females produced the second asexual generation, named sexuparae. Adult sexuparae, reared in short photoperiods, produced the third generation of sexual females oviparous, males, and asexual viviparous females.
Aphid sexuparae were synchronised at the fourth instar moult, during a 6 hour window. 24 h after, 400 ng of kinoprene (Sigma Aldrich) diluted in 50 μL of acetone were ectopically applied on sexuparae abdomens of 100 females. The progeny of treated aphids was collected daily for 5 days. When they reached adulthood, collected aphids were dissected and their reproductive type was determined according to the type of ovaries they contained within their abdomen: sexual ovaries that are only made of haploid eggs, ambiphasic ovaries that contain a mixture of eggs and embryos, and asexual ovaries that are only constituted of embryos at different stages of development. Ectopic applications of 50 μL of acetone were performed as a negative control.
Sampling, RNA isolation and microarray hybridization
Synchronized sexuparae were randomly divided into 2 batches and treated with kinoprene or acetone (as a control) 24 h after fourth instar moult (Additional file 4). 25 treated sexuparae were collected 24, 48, and 72 hours after kinoprene or acetone application. For each condition, the 5 most developed embryos were isolated from each of the 25 treated sexuparae by dissection, pooled together, frozen into liquid nitrogen and stored at -80°C. These most developed embryos within treated sexuparae collected 24, 48, and 72 hours after treatment correspond to the developemental stage 18, 19 or 20 respectively. This procedure was repeated 5 times to generate as much independent biological replicates. Total RNAs were isolated from each sample by using the RNeasy Mini kit (Qiagen) according to manufacturer's instructions. RNA quality was checked on Bioanalyser (Agilent) and quantified on Nanodrop (Thermo scientific). For each sample, 20 μg of total RNAs were sent to the NimbleGen expression array platform (Roche). Double stranded cDNA synthesis and Cy3 end-labelling were performed by NimbleGen.
Microarrays design and analysis
Custom microarrays were constructed on NimbleGen (Roche) 385K 4-plex (4 × 72,000 probes). 24,011 transcripts were represented by 3 60-mers oligonucleotides probes (MIAMExpress, http://www.ebi.ac.uk/miamexpress/, microarray: INRA-BF2I_A.pisum_Nimblegen-ACYPI_4x72k_v1; Array express accession: A-MEXP-1999). Hybridization and scanning were performed by NimbleGen, providing the final raw data file (experiment: Aphid_embryo, ArrayExpress accession: E-MEXP-3481). The limma 2.16 package  in R 2.9.2  was used for statistical analyses. For each transcript the average of the 3 probe signals was considered. Median normalization between microarrays was performed. A design matrix incorporating effects of treatment K (kinoprene) or A (acetone) and developing stage (18, 19, 20) was constructed. Variance was adjusted for Bayesian fitting of the model. Differential expression was determined with a two-way ANOVA considering the effect of time and treatment between 3 contrasts: 'A18 versus K18', 'A19 versus K19' and 'A20 versus K20'. Significance of differential expression was assigned with a 10% false discovery rate (FDR) . Annotation of differentially expressed transcripts was curated using pea aphid RNA-seq data (published in SRA of NCBI), the software CAP3  and Blastx on NCBI non-redundant protein sequence data bank. Functional annotations were performed using the Gene Ontology (http://www.geneontology.org/).
Riboprobe synthesis for in situ hybridization
Templates for synthesis of riboprobes were obtained from full-length cDNAs collection (ACYPI003103, ACYPI010052) or amplified by RT-PCR and cloned (Additional file 5). Total RNAs were extracted from parthenogenetic virginoparae females with RNeasy plant kit (Qiagen). DNA contaminations were removed by a treatment with RQ1 RNase-free DNAse (Promega). First strand cDNAs were produced from 1 μg of total RNAs by using random primer 9 (New England BioLabs) and SuperScript® III Reverse Transcriptase (Invitrogen) following the supplier's instructions. cDNAs were used as a matrix for PCR amplification with specific primers (Additional file 5). Amplified fragments were cloned into the StrataClone PCR Cloning Vector pSC-A-amp/kan (StrataClone) or pENTR Directional TOPO (Invitrogen) and sequenced (Genoscreen). Linear PCR products were amplified from cloned sequence with universal primers M13 and used as a matrix for synthesis of sense and antisense riboprobes by using digoxigenin-labelled dNTPs and the appropriate RNA polymerase (T3/T7/SP6) supplied in the DIG RNA labelling kit (Roche). Remaining DNA was removed with RQ1 RNase-free DNAse treatment (Promega) and riboprobes were purified with the RNeasy mini kit (Qiagen). Riboprobe quality and quantity were checked on Nanodrop (Thermo scientific).
Whole mount in situ hybridization and microscopy
Whole mount in situ hybridization was performed on ovaries of sexual or asexual nymphs. Ovaries were dissected and fixed in 4% paraformaldehyde in 1× PBS at room temperature (RT) for 30 min. Ovaries were washed in 50% methanol for 30 min, then dehydrated and stored in methanol at -20°C. They were rehydrated in graded methanol/PTw (1× PBS, 0.2% Tween-20) solutions (70%; 50%; 30%, 10 min each), post fixated in 4% paraformaldehyde in 1× PBS for 20 min and washed 3 times for 5 min in PTw. Ovaries were washed for 45 min in 1% SDS, 0.5% Tween-20, 50 mM Tris-HCl (pH 7,5), 1 mM EDTA (pH 8), 150 mM NaCl and 5 times for 5 min with PTw. Pre-hybridization was performed 10 min in hyb-wash solution (5× SSC, 50% formamide, 0.1% Tween-20) and 1 hour at 65°C in hyb solution (0.3% SDS, 5× SSC, 50% formamide, 100 μg/mL heparin, 0.1% Tween-20, 100 μg/mL yeast RNA, 10 mM DTT). Hybridization was performed overnight at 65°C with 500 ng/mL of denatured sense or antisense RNA probe diluted in hyb solution. Nonspecific hybridizations were washed off twice at 65°C in hyb-wash solution for 30 min, 3 times in 50% hyb-wash solution/PTw 30 min and 3 times in PTw for 10 min at RT. Ovaries were incubated in blocking solution (0.2% BSA in PTw) for 5 min and 1 h at RT, and bound digoxigenin labelled probes were detected overnight at 4°C with anti-DIG-alkaline phosphatase (AP) Fab fragments (Roche) diluted 1:2000 in blocking solution. Ovaries were washed 4 times 20 min and twice for 1 h in blocking solution at RT. After 3 washes in AP reaction buffer (100 mM Tris (pH 9,5), 100 mM NaCl, 5 mM MgCl2, 0.2% Tween-20), signal was revealed with 4 μl NitroBlue Tetrazolium/5-Bromo-4-Chloro-3-Indolyl Phosphate (NBT/BCIP) Stock Solution (Roche)/ml AP reaction buffer. Finally, ovaries were rinsed at least 3 times for 5 min in PTw, dehydrated 5 min in methanol, and rinsed twice again 5 min in PTw, before mounting in 70% glycerol in PBS. Samples were photographed with a microscope Nikon 90i connected to a Nikon type DS-Ri1 camera or with Olympus BX61 connected to Nikon DS-Fi1 camera.
For the fluorescence protocol, the same procedure was followed until hybridization step. Then TSA™ Biotin System (Perkin Elmer, Waltham, USA) was used for detection as follows. After overnight hybridization of the probes, nonspecific probes were washed off twice at 65°C in hyb-wash solution for 30 min, once in 50% hyb-wash solution/TNT buffer (0.1 M Tris-HCl (pH 7.5), 0.15 M NaCl, 0.05% Tween-20) for 30 min at RT, and 3 times in TNT buffer for 10 min at RT. Ovaries were blocked with the TSA™ blocking reagent, for 30 min at RT. Bound digoxigenin labelled probes were detected with HRP conjugated anti-DIG (Boehringer-Mannheim) (1:250) diluted in TNT buffer at 4°C for overnight incubation. Ovaries were washed four times in TNT buffer for 15 min. The TSA™ Biotin System Amplification was used by incubation 20 min into a diluted solution of the Biotinyl Tyramide (Amplification Reagent) in 1× Amplification Diluent (1:50). HRP catalyzes the formation of TSA free radicals, which form covalent bounds to tyrosine residues proximal to HRP. They were then washed 3 times 15 min in TNT buffer, and 4 times 20 min in PTw. Ovaries were incubated 2 h into Alexa 594-conjugated Streptavidin (Invitrogen) diluted in TNT buffer (1:500) and rinsed into PTw before performing a nuclear stain with TO-PRO3® (1:1000 in PTw) for 30 min. Finally, ovaries were washed 3 times in PTw before mounting. Images were acquired and processed using an Olympus FV1000 confocal microscope.
Phalloidin and propidium iodide staining
Oviparous ovaries were dissected and fixed as previously described. Samples were washed 3 times in PTw, followed by PTw for 30 min and 2% Normal Goat Serum (NGS)/BPTw for 60 min. Samples were then incubated overnight at 4°C with Alexa Fluor 488 phalloidin (Invitrogen) at a concentration of 1:200 to visualize F-actin. The next day samples were washed 4 times in PTw, and placed in 20 μg/mL RNaseA (QIAGEN)/PTw for 3 h at RT. Samples were washed twice in PTw, and incubated with 10 μg/mL propidium iodide (Invitrogen) in 20 μg/mL RNaseA/PTw for 60 min at RT to visualize DNA. Samples were washed 4 times in PTw and mounted in Vectashield (Vector Laboratories). Images were acquired and processed using an Olympus FV1000 confocal microscope.
AG, SJP and DT: INRA, UMR 1099 BiO3P, Biologie des Organismes et des Populations appliquée à la Protection des Plantes, 35653 Le Rheu cedex, France
SS and TH: Okazaki Institute for Integrative Bioscience, National Institute for Basic Biology, National Institutes of Natural Sciences, Higashiyama, Myodaiji, Okazaki, 444-8787, Japan
Hubert Charles, Stefano Colella (INRA INSA UMR BF2I, France) and Grégory Guernec (INRA UMR SCRIBE, France) are acknowledged for microarrays design and analyses. We acknowledge Sylvie Hudaverdian (INRA UMR BiO3P, France), Chantal Cauty and Jérôme Burgeon (INRA UMR SCRIBE, France) for technical support. Jim Hardie (Imperial College of London, UK) provided advices on the use of JH analogs. Owain Edwards (CSIRO, Australia), Jean-René Huynh (CNRS Institut Curie, France), Luc Paillard (CNRS, UMR 6061, France), Gaël Le Trionnaire and Jean-Christophe Simon (INRA UMR BiO3P, France) provided valuable comments on the manuscript. This work was funded by ANR programs 'Aphicible' and 'Holocentrisme', by the Japan Society for the Promotion of Sciences and by french "Ministère de l'enseignement supérieur et de la recherche".
- Lost sex, the evolutionary biology of parthenogenesis. Edited by: Schön I, Martens K, van Dijk P. 2009, Springer
- Simon JC, Delmotte F, Rispe C, Crease T: Phylogenetic relationships between parthenogens and their sexual relatives: the possible routes to parthenogenesis in animals. Biol J Linn Soc. 2003, 79: 151-163. 10.1046/j.1095-8312.2003.00175.x.View Article
- Schwander T, Crespi BJ: Multiple direct transitions from sexual reproduction to apomistic parthenogenesis in tinema stick insects. Evolution. 2009, 63: 84-103. 10.1111/j.1558-5646.2008.00524.x.View ArticlePubMed
- Hörandl E: Geographical Parthenogenesis: Opportunities for Asexuality. Lost sex, the evolutionary biology of parthenogenesis. Edited by: Schön I, Martens, K., and van Dijk, P. 2009, Springer, 161-186.
- Dixon AFG: Aphid ecology. 1998, Chapman & Hall, 2
- Le Trionnaire G, Hardie J, Jaubert-Possamai S, Simon JC, Tagu D: Shifting from clonal to sexual reproduction in aphids: physiological and developmental aspects. Biol Cell. 2008, 100: 441-451. 10.1042/BC20070135.View ArticlePubMed
- Blackman RL: Reproduction, cytogenetics and development. Aphids: their biology, natural enemies and control. Edited by: Minks AKaH, A.P. 1987, Amsterdam: Elsevier, 2A: 163-195.
- 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 Part B. 2003, 295B: 59-81. 10.1002/jez.b.3.View Article
- Stevens NM: A study of the germ cells of Aphis rosae and Aphis oenotherae. Journal Experimental Zoology. 1904, 313-333.
- International Aphid Genomics Consortium: Genome sequence of the pea aphid Acyrthosiphon pisum. PLoS Biol. 2010, 8: e1000313-10.1371/journal.pbio.1000313.View Article
- Corbitt TS, Hardie J: Juvenile-hormone effects on polymorphism in the pea aphid, Acyrthosiphon pisum. Entomol Exp Appl. 1985, 38: 131-135. 10.1111/j.1570-7458.1985.tb03509.x.View Article
- Chicoine J, Benoit P, Gamberi C, Paliouras M, Simonelig M, Lasko P: Bicaudal-C recruits CCR4-NOT deadenylase to target mRNAs and regulates oogenesis, cytoskeletal organization, and its own expression. Dev Cell. 2007, 13: 691-704. 10.1016/j.devcel.2007.10.002.View ArticlePubMed
- Benoit P, Papin C, Kwak JE, Wickens M, Simonelig M: PAP- and GLD-2-type poly(A) polymerases are required sequentially in cytoplasmic polyadenylation and oogenesis in Drosophila. Development. 2008, 135: 1969-1979. 10.1242/dev.021444.View ArticlePubMed
- Eckmann CR, Crittenden SL, Suh N, Kimble J: GLD-3 and control of the mitosis/meiosis decision in the germline of Caenorhabditis elegans. Genetics. 2004, 168: 147-160. 10.1534/genetics.104.029264.PubMed CentralView ArticlePubMed
- Kim JH, Richter JD: Opposing polymerase-deadenylase activities regulate cytoplasmic polyadenylation. Mol Cell. 2006, 24: 173-183. 10.1016/j.molcel.2006.08.016.View ArticlePubMed
- Castagnetti S, Ephrussi A: Orb and a long poly(A) tall are required for efficient oskar translation at the posterior pole of the Drosophila oocyte. Development. 2003, 130: 835-843. 10.1242/dev.00309.View ArticlePubMed
- Temme C, Zhang LB, Kremmer E, Ihling C, Chartier A, Sinz A, Simonelig M, Wahle E: Subunits of the Drosophila CCR4-NOT complex and their roles in mRNA deadenylation. RNA-Publ RNA Soc. 2010, 16: 1356-1370.View Article
- Izzo A, Kamieniarz K, Schneider R: The histone H1 family: specific members, specific functions?. Biol Chem. 2008, 389: 333-343. 10.1515/BC.2008.037.View ArticlePubMed
- Godde JS, Ura K: Dynamic alterations of linker histone variants during development. Int J Dev Biol. 2009, 53: 215-229. 10.1387/ijdb.082644jg.View ArticlePubMed
- Qin WH, Leonhardt H, Spada F: Usp7 and Uhrf1 Control Ubiquitination and Stability of the Maintenance DNA Methyltransferase Dnmt1. J Cell Biochem. 2011, 112: 439-444. 10.1002/jcb.22998.View ArticlePubMed
- Hashimoto H, Horton JR, Zhang X, Cheng XD: UHRF1, a modular multi-domain protein, regulates replication-coupled crosstalk between DNA methylation and histone modifications. Epigenetics. 2009, 4: 8-14. 10.4161/epi.4.1.7370.PubMed CentralView ArticlePubMed
- Beller M, Bulankina AV, Hsiao HH, Urlaub H, Jackle H, Kuhnlein RP: PERILIPIN-Dependent Control of Lipid Droplet Structure and Fat Storage in Drosophila. Cell Metab. 2010, 12: 521-532. 10.1016/j.cmet.2010.10.001.View ArticlePubMed
- Cermelli S, Guo Y, Gross SP, Welte MA: The lipid-droplet proteome reveals that droplets are a protein-storage depot. Curr Biol. 2006, 16: 1783-1795. 10.1016/j.cub.2006.07.062.View ArticlePubMed
- Srinivasan DG, Fenton B, Jaubert-Possamai S, Jaouannet M: Analysis of meiosis and cell cycle genes of the facultatively asexual pea aphid, Acyrthosiphon pisum (Hemiptera: Aphididae). Insect Mol Biol. 2010, 19: 229-239.View ArticlePubMed
- Srinivasan D, Año L, Fenton B, Jaubert-Possamai S, Stern DL: Molecular mechanism of facultative parthenogenesis in the pea aphid, Acyrthosiphon pisum. 2009, 8th International symposium on aphids; 8-12 June 2009; Catania (Italy)
- Kolonin MG, Finley RL: A role for cyclin J in the rapid nuclear division cycles of early Drosophila embryogenesis. Dev Biol. 2000, 227: 661-672. 10.1006/dbio.2000.9916.View ArticlePubMed
- Szalontai T, Gaspar I, Belecz I, Kerekes I, Erdelyi M, Boros I, Szabad J: Horka(D), a Chromosome Instability-Causing Mutation in Drosophila, Is a Dominant-Negative Allele of lodestar. Genetics. 2009, 181: 367-377.PubMed CentralView ArticlePubMed
- Riparbelli MG, Tagu D, Bonhomme J, Callaini G: Aster self-organization at meiosis: a conserved mechanism in insect parthenogenesis?. Dev Biol. 2005, 278: 220-230. 10.1016/j.ydbio.2004.11.009.View ArticlePubMed
- Smyth GK: Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Statistical Applications in Genetics and Molecular Biology. 2004, 1-26.
- R Development Core Team: R: a language and environment for statistical computing. 2009, Vienna, Austria: R foundation for statistical computing
- Benjamini Y, Hochberg Y: Controlling the false discovery rate - A practical and powerful approach ot multiple testing. J R Stat Soc Ser B-Methodol. 1995, 57: 289-300.
- Huang XQ, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9: 868-877. 10.1101/gr.9.9.868.PubMed CentralView ArticlePubMed
- Shigenobu S, Richards S, Cree AG, Morioka M, Fukatsu T, Kudo T, Miyagishima S, Gibbs RA, Stern DL, Nakabachi A: A full-length cDNA resource for the pea aphid, Acyrthosiphon pisum. Insect Mol Biol. 2010, 19: 23-31.PubMed CentralView ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.