Silencing the conserved small nuclear ribonucleoprotein SmD1 target gene alters susceptibility to root-knot nematodes in plants.

Root-knot nematodes (RKNs) are among the most damaging pests of agricultural crops. Meloidogyne is an extremely polyphagous genus of nematodes that can infect thousands of plant species. A few genes for resistance (R-genes) to RKN suitable for use in crop breeding have been identified, but virulent strains and species of RKN have emerged that render these R-genes ineffective. Secretion of RKN effectors targeting plant functions mediates the reprogramming of root cells into specialized feeding cells, the giant cells, essential for RKN development and reproduction. Conserved targets among plant species define the more relevant strategies for controlling nematode infection. The EFFECTOR18 (EFF18) protein from M. incognita interacts with the spliceosomal small nuclear ribonucleoprotein D1 (SmD1) in Arabidopsis (Arabidopsis thaliana), disrupting its function in alternative splicing regulation and modulating the giant cell transcriptome. We show here that EFF18 is a conserved RKN-specific effector that targets this conserved spliceosomal SmD1 protein in Solanaceae. This interaction modulates alternative splicing events produced by tomato (Solanum lycopersicum) in response to M. incognita infection. The alteration of SmD1 expression by virus-induced gene silencing in Solanaceae affects giant cell formation and nematode development. Thus, our work defines a promising conserved SmD1 target gene to develop broad resistance for the control of Meloidogyne spp. in plants.

AUXIN RESPONSIVE FACTOR8 regulates development of the feeding site induced by root-knot nematodes in tomato.

Root-knot nematodes (RKN) from the genus Meloidogyne induce the dedifferentiation of root vascular cells into giant multinucleate feeding cells. These feeding cells result from an extensive reprogramming of gene expression, and auxin is known to be a key player in their development. However, little is known about how the auxin signal is transmitted during giant cell development. Integrative analyses combining transcriptome and small non-coding RNA datasets with the specific sequencing of cleaved transcripts identified genes targeted by miRNAs in tomato (Solanum lycopersicum) galls. The two auxin-responsive transcription factors ARF8A and ARF8B, and their miRNA167 regulators, were identified as robust gene-miRNA pair candidates to be involved in the tomato response to M. incognita. Spatiotemporal expression analysis using promoter-ß-glucuronidase (GUS) fusions showed the up-regulation of ARF8A and ARF8B in RKN-induced feeding cells and surrounding cells. The generation and phenotyping of CRISPR (clustered regularly interspaced palindromic repeats) mutants demonstrated the role of ARF8A and ARF8B in giant cell development and allowed the characterization of their downstream regulated genes.

Copper microRNAs modulate the formation of giant feeding cells induced by the root knot nematode Meloidogyne incognita in Arabidopsis thaliana.

Root-knot nematodes (RKNs) are root endoparasites that induce the dedifferentiation of a few root cells and the reprogramming of their gene expression to generate giant hypermetabolic feeding cells. We identified two microRNA families, miR408 and miR398, as upregulated in Arabidopsis thaliana and Solanum lycopersicum roots infected by RKNs. In plants, the expression of these two conserved microRNA families is known to be activated by the SPL7 transcription factor in response to copper starvation. By combining functional approaches, we deciphered the network involving these microRNAs, their regulator and their targets. MIR408 expression was located within nematode-induced feeding cells like its regulator SPL7 and was regulated by copper. Moreover, infection assays with mir408 and spl7 knockout mutants or lines expressing targets rendered resistant to cleavage by miR398 demonstrated the essential role of the SPL7/MIR408/MIR398 module in the formation of giant feeding cells. Our findings reveal how perturbation of plant copper homeostasis, via the SPL7/MIR408/MIR398 module, modulates the development of nematode-induced feeding cells.

Characterization of siRNAs clusters in Arabidopsis thaliana galls induced by the root-knot nematode Meloidogyne incognita.

BACKGROUND: Root-knot nematodes (RKN), genus Meloidogyne, are plant parasitic worms that have the ability to transform root vascular cylinder cells into hypertrophied, multinucleate and metabolically over-active feeding cells. Redifferentiation into feeding cells is the result of a massive transcriptional reprogramming of root cells targeted by RKN. Since RKN are able to induce similar feeding cells in roots of thousands of plant species, these worms are thought to manipulate essential and conserved plant molecular pathways. RESULTS: Small non-coding RNAs of uninfected roots and infected root galls induced by M. incognita from Arabidopsis thaliana were sequenced by high throughput sequencing. SiRNA populations were analysed by using the Shortstack algorithm. We identified siRNA clusters that are differentially expressed in infected roots and evidenced an over-representation of the 23-24 nt siRNAs in infected tissue. This size corresponds to heterochromatic siRNAs (hc-siRNAs) which are known to regulate expression of transposons and genes at the transcriptional level, mainly by inducing DNA methylation. CONCLUSIONS: Correlation of siRNA clusters expression profile with transcriptomic data identified several protein coding genes that are candidates to be regulated by siRNAs at the transcriptional level by RNA directed DNA methylation (RdDM) pathway either directly or indirectly via silencing of neighbouring transposable elements.

Characterization of microRNAs from Arabidopsis galls highlights a role for miR159 in the plant response to the root-knot nematode Meloidogyne incognita.

Root knot nematodes (RKN) are root parasites that induce the genetic reprogramming of vascular cells into giant feeding cells and the development of root galls. MicroRNAs (miRNAs) regulate gene expression during development and plant responses to various stresses. Disruption of post-transcriptional gene silencing in Arabidopsis ago1 or ago2 mutants decrease the infection rate of RKN suggesting a role for this mechanism in the plant-nematode interaction. By sequencing small RNAs from uninfected Arabidopsis roots and from galls 7 and 14 d post infection with Meloidogyne incognita, we identified 24 miRNAs differentially expressed in gall as putative regulators of gall development. Moreover, strong activity within galls was detected for five miRNA promoters. Analyses of nematode development in an Arabidopsis miR159abc mutant had a lower susceptibility to RKN, suggesting a role for the miR159 family in the plant response to M. incognita. Localization of mature miR159 within the giant and surrounding cells suggested a role in giant cell and gall. Finally, overexpression of miR159 in galls at 14 d post inoculation was associated with the repression of the miR159 target MYB33 which expression is restricted to the early stages of infection. Overall, these results implicate the miR159 in plant responses to RKN.

Differentially expressed small RNAs in Arabidopsis galls formed by Meloidogyne javanica: a functional role for miR390 and its TAS3-derived tasiRNAs.

Root-knot nematodes (RKNs) induce inside the vascular cylinder the giant cells (GCs) embedded in the galls. The distinctive gene repression in early-developing GCs could be facilitated by small RNAs (sRNA) such as miRNAs, and/or epigenetic mechanisms mediated by 24nt-sRNAs, rasiRNAs and 21-22nt-sRNAs. Therefore, the sRNA-population together with the role of the miR390/TAS3/ARFs module were studied during early gall/GC formation. Three sRNA libraries from 3-d-post-inoculation (dpi) galls induced by Meloidogyne javanica in Arabidopsis and three from uninfected root segments were sequenced following Illumina-Solexa technology. pMIR390a::GUS and pTAS3::GUS lines were assayed for nematode-dependent promoter activation. A sensor line indicative of TAS3-derived tasiRNAs binding to the ARF3 sequence (pARF3:ARF3-GUS) together with a tasiRNA-resistant ARF3 line (pARF3:ARF3m-GUS) were used for functional analysis. The sRNA population showed significant differences between galls and controls, with high validation rate and correspondence with their target expression: 21-nt sRNAs corresponding mainly to miRNAs were downregulated, whilst 24-nt-sRNAs from the rasiRNA family were mostly upregulated in galls. The promoters of MIR390a and TAS3, active in galls, and the pARF3:ARF3-GUS line, indicated a role of TAS3-derived-tasiRNAs in galls. The regulatory module miR390/TAS3 is necessary for proper gall formation possibly through auxin-responsive factors, and the abundance of 24-nt sRNAs (mostly rasiRNAs) constitutes a gall hallmark.

Evolutionary study of duplications of the miRNA machinery in aphids associated with striking rate acceleration and changes in expression profiles.

BACKGROUND: The sequencing of the genome of the pea aphid Acyrthosiphon pisum revealed an unusual expansion of the miRNA machinery, with two argonaute-1, two dicer-1 and four pasha gene copies. In this report, we have undertaken a deeper evolutionary analysis of the phylogenetic timing of these gene duplications and of the associated selective pressures by sequencing the two copies of ago-1 and dcr-1 in different aphid species of the subfamily Aphidinae. We have also carried out an analysis of the expression of both copies of ago-1 and dcr-1 by semi-quantitative PCR in different morphs of the pea aphid life cycle. RESULTS: The analysis has shown that the duplication of ago-1 occurred in an ancestor of the subfamily Aphidinae while the duplication of dcr-1 appears to be more recent. Besides, it has confirmed a pattern of one conserved copy and one accelerated copy for both genes, and has revealed the action of positive selection on several regions of the fast-evolving ago-1b. On the other hand, the semi-quantitative PCR experiments have revealed a differential expression of these genes between the morphs of the parthenogenetic and the sexual phases of Acyrthosiphon pisum. CONCLUSIONS: The discovery of these gene duplications in the miRNA machinery of aphids opens new perspectives of research about the regulation of gene expression in these insects. Accelerated evolution, positive selection and differential expression affecting some of the copies of these genes suggests the possibility of a neofunctionalization of these duplicates, which might play a role in the display of the striking phenotypic plasticity of aphids.

Sexual and asexual oogenesis require the expression of unique and shared sets of genes in the insect Acyrthosiphon pisum.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

Expansion of the miRNA pathway in the hemipteran insect Acyrthosiphon pisum.

The pathways that allow short noncoding RNAs such as the microRNAs (miRNAs) to mediate gene regulation and control critical cellular and developmental processes involve a limited number of key protein components. These proteins are the Dicer-like RNases, double-stranded RNA (dsRNA)-binding proteins, and the Argonaute (AGO) proteins that process stem-loop hairpin transcripts of endogenous genes to generate miRNAs or long dsRNA precursors (either exogenous or endogenous). Comparative genomics studies of metazoans have shown the pathways to be highly conserved overall; the major difference observed is that the vertebrate pathways overlap in sharing a single Dicer (DCR) and AGO proteins, whereas those of insects appear to be parallel, with distinct Dicers and AGOs required for each pathway. The genome of the pea aphid is the first available for a hemipteran insect and discloses an unexpected expansion of the miRNA pathway. It has two copies of the miRNA-specific dicr-1 and ago1 genes and four copies of pasha a cofactor of drosha involved in miRNA biosynthesis. For three of these expansions, we showed that one copy of the genes diverged rapidly and in one case (ago1b) shows signs of positive selection. These expansions occurred concomitantly within a brief evolutionary period. The pea aphid, which reproduces by viviparous parthenogenesis, is able to produce several adapted phenotypes from one single genotype. We show by reverse transcriptase-polymerase chain reaction that all the duplicated copies of the miRNA machinery genes are expressed in the different morphs. Investigating the function of these novel genes offers an exciting new challenge in aphid biology.

The Meloidogyne incognita Nuclear Effector MiEFF1 Interacts With Arabidopsis Cytosolic Glyceraldehyde-3-Phosphate Dehydrogenases to Promote Parasitism.

Root-knot nematodes are obligate endoparasites that maintain a biotrophic relationship with their hosts over a period of several weeks. They induce the differentiation of root cells into specialized multinucleate hypertrophied feeding cells known as giant cells. Nematode effectors synthesized in the esophageal glands and injected into the plant tissue through the syringe-like stylet play a key role in giant cell ontogenesis. The Meloidogyne incognita MiEFF1 is one of the rare effectors of phytopathogenic nematodes to have been located in vivo in feeding cells. This effector specifically targets the giant cell nuclei. We investigated the Arabidopsis functions modulated by this effector, by using a yeast two-hybrid approach to identify its host targets. We characterized a universal stress protein (USP) and cytosolic glyceraldehyde-3-phosphate dehydrogenases (GAPCs) as the targets of MiEFF1. We validated the interaction of MiEFF1 with these host targets in the plant cell nucleus, by bimolecular fluorescence complementation (BiFC). A functional analysis with Arabidopsis GUS reporter lines and knockout mutant lines showed that GAPCs were induced in giant cells and that their non-metabolic functions were required for root-knot nematode infection. These susceptibility factors are potentially interesting targets for the development of new root-knot nematode control strategies.

Bioinformatic prediction, deep sequencing of microRNAs and expression analysis during phenotypic plasticity in the pea aphid, Acyrthosiphon pisum.

BACKGROUND: Post-transcriptional regulation in eukaryotes can be operated through microRNA (miRNAs) mediated gene silencing. MiRNAs are small (18-25 nucleotides) non-coding RNAs that play crucial role in regulation of gene expression in eukaryotes. In insects, miRNAs have been shown to be involved in multiple mechanisms such as embryonic development, tissue differentiation, metamorphosis or circadian rhythm. Insect miRNAs have been identified in different species belonging to five orders: Coleoptera, Diptera, Hymenoptera, Lepidoptera and Orthoptera. RESULTS: We developed high throughput Solexa sequencing and bioinformatic analyses of the genome of the pea aphid Acyrthosiphon pisum in order to identify the first miRNAs from a hemipteran insect. By combining these methods we identified 149 miRNAs including 55 conserved and 94 new miRNAs. Moreover, we investigated the regulation of these miRNAs in different alternative morphs of the pea aphid by analysing the expression of miRNAs across the switch of reproduction mode. Pea aphid microRNA sequences have been posted to miRBase: http://microrna.sanger.ac.uk/sequences/. CONCLUSIONS: Our study has identified candidates as putative regulators involved in reproductive polyphenism in aphids and opens new avenues for further functional analyses.

Transcriptomic analysis of intestinal genes following acquisition of pea enation mosaic virus by the pea aphid Acyrthosiphon pisum.

Viruses in the family Luteoviridae are strictly transmitted by aphids in a non-propagative, circulative and persistent mode. Virions ingested by aphids successively cross the gut and the accessory salivary gland epithelia before being released, together with saliva, into the plant vasculature. Virion transport through aphid cells occurs by a transcytosis mechanism. This study conducted a transcriptomic analysis of intestinal genes of the pea aphid Acyrthosiphon pisum following uptake of pea enation mosaic virus. Among the 7166 transcripts analysed, 128 were significantly regulated (105 genes downregulated and 23 upregulated). Of these genes, 5 % were involved in intracellular trafficking, endocytosis and signal transduction, three important steps in the internalization and transport of virions. The limited levels of downregulation (maximum of 3.45-fold) and upregulation (maximum of 1.37-fold) suggest that the virus hijacks a constitutive endocytosis-exocytosis mechanism without heavily perturbing cell metabolism. Although limited to about 20 % of the pea aphid genes, this work represents the first large-scale analysis of aphid gene regulation following virus acquisition. A better knowledge of this virus-vector interaction will be possible only when tools representing the complete genomic capacity of the aphid become available.

Shifting from clonal to sexual reproduction in aphids: physiological and developmental aspects.

Developmental biology is one of the fastest growing and fascinating research fields in life sciences. Among the wide range of embryonic development, a fundamental difference exists between organisms with sexual or asexual development. Aphids are unusual organisms which display alternative pathways of sexual and asexual development, the orientation of the pathway being determined by environmental conditions. These insects offer an adapted system in which to study developmental plasticity, because a side-by-side comparison of sexual and asexual development can be made in individuals with the same genotype. In this review, we describe the developmental mechanisms that have evolved in aphids for alternative sexual and asexual reproduction. In particular, we discuss how environmental cues orientate the reproductive mode of aphids from signal perception to endocrine regulation, and propose a comparative analysis of sexual and asexual gametogenesis and embryogenesis, which has been possible due to the development of molecular methods. As a result of the recent development of genomic resources in aphids, we expect these species will permit major advances in the study of the genomic basis underlying the choice of developmental fate and multiple reproduction strategies.

MicroRNAs, New Players in the Plant-Nematode Interaction.

Plant-parasitic root-knot and cyst nematodes are microscopic worms that cause severe damage to crops and induce major agricultural losses worldwide. These parasites penetrate into host roots and induce the formation of specialized feeding structures, which supply the resources required for nematode development. Root-knot nematodes induce the redifferentiation of five to seven root cells into giant multinucleate feeding cells, whereas cyst nematodes induce the formation of a multinucleate syncytium by targeting a single root cell. Transcriptomic analyses have shown that the induction of these feeding cells by nematodes involves an extensive reprogramming of gene expression within the targeted root cells. MicroRNAs are small noncoding RNAs that act as key regulators of gene expression in eukaryotes by inducing the posttranscriptional silencing of protein coding genes, including many genes encoding transcription factors. A number of microRNAs (miRNAs) displaying changes in expression in root cells in response to nematode infection have recently been identified in various plant species. Modules consisting of miRNAs and the transcription factors they target were recently shown to be required for correct feeding site formation. Examples include miR396 and GRF in soybean syncytia and miR159 and MYB33 in Arabidopsis giant cells. Moreover, some conserved miRNA/target modules seem to have similar functions in feeding site formation in different plant species. These miRNAs may be master regulators of the reprogramming of expression occurring during feeding site formation. This review summarizes current knowledge about the role of these plant miRNAs in plant-nematode interactions.

In situ Hybridization (ISH) in Preparasitic and Parasitic Stages of the Plant-parasitic Nematode Meloidogyne spp.

The spatio-temporal expression pattern of a gene provides important indications to better understand its biological function. In situ hybridization (ISH) uses a labeled complementary single-stranded RNA or DNA probe to localize gene transcripts in a whole organism, a whole organ or a section of tissue. We adapted the ISH technique to the plant parasite Meloidogyne spp. (root-knot nematode) to visualize RNAs both in free-living preparasitic juveniles and in parasitic stages settled in the plant tissues. We describe each step of the probe synthesis, digoxigenin (DIG) labeling, nematode extraction from plant tissue, and ISH procedure.

Gene knockdown by RNAi in the pea aphid Acyrthosiphon pisum.

BACKGROUND: RNA interference (RNAi) is a powerful method to inhibit gene expression in a sequence specific manner. RESULTS: Here, we described the development of RNAi by micro-injection of double-stranded RNA (dsRNA) in the pea aphid Acyrthosiphon pisum. Injection of dsRNA into whole aphid body induced the silencing of two marker genes with different expression patterns: the ubiquitously expressed Ap-crt genes encoding a calreticulin and the gut specific Ap-cath-L gene encoding a cathepsin-L. Time-course analysis of the silencing showed similar temporal patterns for both genes: inhibition started at 1 day after injection, reached its maximum at 5 days and stopped at 7 days. A comparable 40% decrease of gene expression was observed for Ap-crt and Ap-cath-L. CONCLUSION: The pea aphid is the first Hemipteran insect for which genome sequence will be available soon. The gene knockdown technique developed in this study will be an essential post-genomic tool for further investigations in aphidology.

Gall-forming root-knot nematodes hijack key plant cellular functions to induce multinucleate and hypertrophied feeding cells.

Among plant-parasitic nematodes, the root-knot nematodes (RKNs) of the Meloidogyne spp. are the most economically important genus. RKN are root parasitic worms able to infect nearly all crop species and have a wide geographic distribution. During infection, RKNs establish and maintain an intimate relationship with the host plant. This includes the creation of a specialized nutritional structure composed of multinucleate and hypertrophied giant cells, which result from the redifferentiation of vascular root cells. Giant cells constitute the sole source of nutrients for the nematode and are essential for growth and reproduction. Hyperplasia of surrounding root cells leads to the formation of the gall or root-knot, an easily recognized symptom of plant infection by RKNs. Secreted effectors produced in nematode salivary glands and injected into plant cells through a specialized feeding structure called the stylet play a critical role in the formation of giant cells. Here, we describe the complex interactions between RKNs and their host plants. We highlight progress in understanding host plant responses, focusing on how RKNs manipulate key plant processes and functions, including cell cycle, defence, hormones, cellular scaffold, metabolism and transport.

Neonicotinoid binding, toxicity and expression of nicotinic acetylcholine receptor subunits in the aphid Acyrthosiphon pisum.

Neonicotinoid insecticides act on nicotinic acetylcholine receptor and are particularly effective against sucking pests. They are widely used in crops protection to fight against aphids, which cause severe damage. In the present study we evaluated the susceptibility of the pea aphid Acyrthosiphon pisum to the commonly used neonicotinoid insecticides imidacloprid (IMI), thiamethoxam (TMX) and clothianidin (CLT). Binding studies on aphid membrane preparations revealed the existence of high and low-affinity binding sites for [3H]-IMI (Kd of 0.16 ± 0.04 nM and 41.7 ± 5.9 nM) and for the nicotinic antagonist [125I]-α-bungarotoxin (Kd of 0.008 ± 0.002 nM and 1.135 ± 0.213 nM). Competitive binding experiments demonstrated that TMX displayed a higher affinity than IMI for [125I]-α-bungarotoxin binding sites while CLT affinity was similar for both [125I]-α-bungarotoxin and [3H]-IMI binding sites. Interestingly, toxicological studies revealed that at 48 h, IMI (LC50 = 0.038 µg/ml) and TMX (LC50 = 0.034 µg/ml) were more toxic than CLT (LC50 = 0.118 µg/ml). The effect of TMX could be associated to its metabolite CLT as demonstrated by HPLC/MS analysis. In addition, we found that aphid larvae treated either with IMI, TMX or CLT showed a strong variation of nAChR subunit expression. Using semi-quantitative PCR experiments, we detected for all insecticides an increase of Apisumα10 and Apisumß1 expressions levels, whereas Apisumß2 expression decreased. Moreover, some other receptor subunits seemed to be differently regulated according to the insecticide used. Finally, we also demonstrated that nAChR subunit expression differed during pea aphid development. Altogether these results highlight species specificity that should be taken into account in pest management strategies.

New insight into the RNA interference response against cathepsin-L gene in the pea aphid, Acyrthosiphon pisum: molting or gut phenotypes specifically induced by injection or feeding treatments.

RNA interference (RNAi) has been widely and successfully used for gene inactivation in insects, including aphids, where dsRNA administration can be performed either by feeding or microinjection. However, several aspects related to the aphid response to RNAi, as well as the influence of the administration method on tissue response, or the mixed success to observe phenotypes specific to the gene targeted, are still unclear in this insect group. In the present study, we made the first direct comparison of two administration methods (injection or feeding) for delivery of dsRNA targeting the cathepsin-L gene in the pea aphid, Acyrthosiphon pisum. In order to maximize the possibility of discovering specific phenotypes, the effect of the treatment was analyzed in single individual aphids at the level of five body compartments: the bacteriocytes, the gut, the embryonic chains, the head and the remaining body carcass. Our analysis revealed that gene expression knockdown effect in each single body compartment was dependent on the administration method used, and allowed us to discover new functions for the cathepsin-L gene in aphids. Injection of cathepsin-L dsRNA was much more effective on carcass and head, inducing body morphology alterations, and suggesting a novel role of this gene in the molting of these insects. Administration by feeding provoked cathepsin-L knockdown in the gut and specific gut epithelial cell alteration, therefore allowing a better characterization of tissue specific role of this gene in aphids.

Transcriptomic profiling of the reproductive mode switch in the pea aphid in response to natural autumnal photoperiod.

Aphids are among the rare organisms that can change their reproductive mode across their life cycle. During spring and summer they reproduce clonally and efficiently by parthenogenesis. At the end of summer aphids perceive the shortening of day length which triggers the production of sexual individuals - males and oviparous females - that will mate and lay overwintering cold-resistant eggs. Recent large scale transcriptomic studies allowed the discovery of transcripts and functions such as nervous and hormonal signaling involved in the early steps of detection and transduction of the photoperiodic signal. Nevertheless these experiments were performed under controlled conditions when the photoperiod was the only varying parameter. To characterize the response of aphids under natural conditions, aphids were reared outdoor both in summer and autumn and material was collected to compare their transcriptomic profile using a cDNA microarray containing around 7000 transcripts. Statistical analyses revealed that close to 5% of these transcripts (367) were differentially expressed at two developmental stages of the process in response to the autumnal environmental conditions. Functional classification of regulated transcripts confirmed the putative contribution of the neuro-endocrine system in the process. Furthermore, these experiments revealed the regulation of transcripts involved in juvenile hormone synthesis and signaling pathway, confirming the key role played by these molecules in the reproductive mode switch. Aphids placed under outdoor conditions were confronted to a range of abiotic factors such as temperature fluctuations which was confirmed by the differential expression of an important proportion of heat shock protein transcripts between the two seasons. Finally, this original approach completed the understanding of genetic programs involved in aphid phenotypic plasticity.