Transcriptomic basis of sex loss in the pea aphid.

BACKGROUND: Transitions from sexual to asexual reproduction are common in eukaryotes, but the underlying mechanisms remain poorly known. The pea aphid-Acyrthosiphon pisum-exhibits reproductive polymorphism, with cyclical parthenogenetic and obligate parthenogenetic lineages, offering an opportunity to decipher the genetic basis of sex loss. Previous work on this species identified a single 840 kb region controlling reproductive polymorphism and carrying 32 genes. With the aim of identifying the gene(s) responsible for sex loss and the resulting consequences on the genetic programs controlling sexual or asexual embryogenesis, we compared the transcriptomic response to photoperiod shortening-the main sex-inducing cue-of a sexual and an obligate asexual lineage of the pea aphid, focusing on heads (where the photoperiodic cue is detected) and embryos (the final target of the cue). RESULTS: Our analyses revealed that four genes (one expressed in the head, and three in the embryos) of the region responded differently to photoperiod in the two lineages. We also found that the downstream genetic programs expressed during embryonic development of a future sexual female encompass ∼1600 genes, among which miRNAs, piRNAs and histone modification pathways are overrepresented. These genes mainly co-localize in two genomic regions enriched in transposable elements (TEs). CONCLUSIONS: Our results suggest that the causal polymorphism(s) in the 840 kb region somehow impair downstream epigenetic and post-transcriptional regulations in obligate asexual lineages, thereby sustaining asexual reproduction.

An integrated protocol for targeted mutagenesis with CRISPR-Cas9 system in the pea aphid.

CRISPR-Cas9 technology is a very efficient functional analysis tool and has been developed in several insects to edit their genome through injection of eggs with guide RNAs targeting coding sequences of genes of interest. However, its implementation in aphids is more challenging. Aphids are major pests of crops worldwide that alternate during their life cycle between clonality and sexual reproduction. The production of eggs after mating of sexual individuals is a single yearly event and is necessarily triggered by a photoperiod decrease. Fertilized eggs then experience an obligate 3-month diapause period before hatching as new clonal colonies. Taking into consideration these particularities, we developed in the pea aphid Acyrthosiphon pisum a step-by-step protocol of targeted mutagenesis based on the microinjection within fertilized eggs of CRISPR-Cas9 components designed for the editing of a cuticular protein gene (stylin-01). This protocol includes the following steps: i) the photoperiod-triggered induction of sexual morphs (2 months), ii) the mating and egg collection step (2 weeks), iii) egg microinjection and melanization, iv) the 3-month obligate diapause, v) the hatching of new lineages from injected eggs (2 weeks) and vi) the maintenance of stable lineages (2 weeks). Overall, this 7-month long procedure was applied to three different crosses in order to estimate the impact of the choice of the genetic combination on egg production dynamics by females as well as hatching rates after diapause. Mutation rates within eggs before diapause were estimated at 70-80%. The hatching rate of injected eggs following diapause ranged from 1 to 11% depending on the cross and finally a total of 17 stable lineages were obtained and maintained clonally. Out of these, 6 lineages were mutated at the defined sgRNAs target sites within stylin-01 coding sequence, either at the two alleles (2 lineages) or at one allele (4 lineages). The final germline transmission rate of the mutations was thus around 35%. Our protocol of an efficient targeted mutagenesis opens the avenue for functional studies through genome editing in aphids.

Comparison of gene repertoires and patterns of evolutionary rates in eight aphid species that differ by reproductive mode.

In theory, the loss of sexual reproduction is expected to result in the accumulation of deleterious mutations. In aphids, two main types of life cycle, cyclic and obligate parthenogenesis, represent respectively "sexual" and "asexual" reproductive modes. We used the complete pea aphid genome and previously published expressed sequence tags (ESTs) from two other aphid species. In addition, we obtained 100,000 new ESTs from five more species. The final set comprised four sexual and four asexual aphid species and served to test the influence of the reproductive mode on the evolutionary rates of genes. We reconstructed coding sequences from ESTs and annotated these genes, discovering a novel peptide gene family that appears to be among the most highly expressed transcripts from several aphid species. From 203 genes found to be 1:1 orthologs among the eight species considered, we established a species tree that partly conflicted with taxonomy (for Myzus ascalonicus). We then used this topology to evaluate the dynamics of evolutionary rates and mutation accumulation in the four sexual and four asexual taxa. No significant increase of the nonsynonymous to synonymous ratio or of nonsynonymous mutation numbers was found in any of the four branches for asexual taxa. We however found a significant increase of the synonymous rate in the branch leading to the asexual species Rhopalosiphum maidis, which could be due to a change in the mutation rate or to an increased number of generations implied by its change of life cycle.

Juvenile hormone titre and related gene expression during the change of reproductive modes in the pea aphid.

Most aphids show reproductive polyphenism, i.e. they alternate their reproductive modes from parthenogenesis to sexual reproduction in response to short photoperiods. Although juvenile hormone (JH) has been considered a likely candidate for regulating the transition from asexual to sexual reproduction after photoperiod sensing, there are few studies investigating the direct relationship between JH titres and the reproductive-mode change. In addition, the sequencing of the pea aphid genome has allowed identification of the genes involved in the JH pathway, which in turn allows us to examine their expression levels in relation to the reproductive-mode change. Using liquid chromatography-mass spectrometry in the pea aphid, JHIII titre was shown to be lower in aphids producing sexual morphs under short-day conditions than in aphids producing parthenogenetic morphs under long-day conditions. The expression levels of genes upstream and downstream of JH action were quantified by real-time quantitative reverse-transcription-PCR across the reproductive-mode change. The expression level of JH esterase, which is responsible for JH degradation, was significantly higher in aphids reared under short-day conditions. This suggests that the upregulation of the JH degradation pathway may be responsible for the lower JHIII titre in aphids exposed to short-days, leading to the production of sexual morphs.

AphidBase: a centralized bioinformatic resource for annotation of the pea aphid genome.

AphidBase is a centralized bioinformatic resource that was developed to facilitate community annotation of the pea aphid genome by the International Aphid Genomics Consortium (IAGC). The AphidBase Information System designed to organize and distribute genomic data and annotations for a large international community was constructed using open source software tools from the Generic Model Organism Database (GMOD). The system includes Apollo and GBrowse utilities as well as a wiki, blast search capabilities and a full text search engine. AphidBase strongly supported community cooperation and coordination in the curation of gene models during community annotation of the pea aphid genome. AphidBase can be accessed at http://www.aphidbase.com.

Neuropeptide and neurohormone precursors in the pea aphid, Acyrthosiphon pisum.

Aphids respond to environmental changes by developing alternative phenotypes with differing reproductive modes. Parthenogenetic reproduction occurs in spring and summer, whereas decreasing day lengths in autumn provoke the production of sexual forms. Changing environmental signals are relayed by brain neuroendocrine signals to the ovarioles. We combined bioinformatic analyses with brain peptidomics and cDNA analyses to establish a catalogue of pea aphid neuropeptides and neurohormones. 42 genes encoding neuropeptides and neurohormones were identified, of which several were supported by expressed sequence tags and/or peptide mass analyses. Interesting features of the pea aphid peptidome are the absence of genes coding for corazonin, vasopressin and sulfakinin and the presence of 10 different genes coding insulin related peptides, one of which appears to be very abundantly expressed.

A functional DNA methylation system in the pea aphid, Acyrthosiphon pisum.

Methylation of cytosine is one of the main epigenetic mechanisms involved in controlling gene expression. Here we show that the pea aphid (Acyrthosiphon pisum) genome possesses homologues to all the DNA methyltransferases found in vertebrates, and that 0.69% (+/-0.25%) of all cytosines are methylated. Identified methylation sites are predominantly restricted to the coding sequence of genes at CpG sites. We identify twelve methylated genes, including genes that interact with juvenile hormone, a key endocrine signal in insects. Bioinformatic prediction using CpG ratios for all predicted genes suggest that a large proportion of genes are methylated within the pea aphid.

Transcriptomic and proteomic analyses of seasonal photoperiodism in the pea aphid.

BACKGROUND: Aphid adaptation to harsh winter conditions is illustrated by an alternation of their reproductive mode. Aphids detect photoperiod shortening by sensing the length of the night and switch from viviparous parthenogenesis in spring and summer, to oviparous sexual reproduction in autumn. The photoperiodic signal is transduced from the head to the reproductive tract to change the fate of the future oocytes from mitotic diploid embryogenesis to haploid formation of gametes. This process takes place in three consecutive generations due to viviparous parthenogenesis. To understand the molecular basis of the switch in the reproductive mode, transcriptomic and proteomic approaches were used to detect significantly regulated transcripts and polypeptides in the heads of the pea aphid Acyrthosiphon pisum. RESULTS: The transcriptomic profiles of the heads of the first generation were slightly affected by photoperiod shortening. This suggests that trans-generation signalling between the grand-mothers and the viviparous embryos they contain is not essential. By analogy, many of the genes and some of the proteins regulated in the heads of the second generation are implicated in visual functions, photoreception and cuticle structure. The modification of the cuticle could be accompanied by a down-regulation of the N-beta-alanyldopamine pathway and desclerotization. In Drosophila, modification of the insulin pathway could cause a decrease of juvenile hormones in short-day reared aphids. CONCLUSION: This work led to the construction of hypotheses for photoperiodic regulation of the switch of the reproductive mode in aphids.

Sex versus parthenogenesis: a transcriptomic approach of photoperiod response in the model aphid Acyrthosiphon pisum (Hemiptera: Aphididae).

Most aphids develop a cyclic parthenogenesis life-cycle. After several generations of viviparous parthenogenetic females, it follows a single annual generation of sexual individuals, usually in autumn, that mate and lay the sexual eggs. Shortening of photoperiod at the end of the summer is a key factor inducing the sexual response. With the survey here reported we aimed at identifying a collection of candidate genes to participate at some point in the cascade of events that lead to the sexual phenotypes. Following a suppression subtractive hybridization methodology (SSH) on the model aphid Acyrthosiphon pisum, we built and characterised two reciprocal cDNA libraries (SDU and SDD) enriched respectively in genes up-regulated or down-regulated by short photoperiod conditions that lead to the sexual response in this aphid species. A total of 557 ESTs were obtained altogether representing 223 non-overlapping contigs. 29% of these were new sequences not present in previous aphid EST libraries. BLAST searches allowed putative identification of about 54% of the contigs present in both libraries. Relative quantification of expression through real-time quantitative PCR demonstrated the differential expression in relation with the photoperiod of 6 genes (3 up-regulated and 3 down-regulated by shortening the day length). Among these, expression of a tubulin gene, two cuticular proteins and a yet unidentified sequence along the day-night cycle was further investigated. Implications for current studies on gene regulation of the dichotomy sex vs. parthenogenesis in aphids are discussed.

Seasonal photoperiodism regulates the expression of cuticular and signalling protein genes in the pea aphid.

Seasonal photoperiodism in aphids is responsible for the spectacular switch from asexual to sexual reproduction. However, little is known on the molecular and physiological mechanisms involved in reproductive mode shift through the action of day length. Earlier works showed that aphid head, but not eyes, directly perceives the photoperiodic signal through the cuticle. In order to identify genes regulating the photoperiodic response, a 3321 cDNA microarray developed for the pea aphid, Acyrthosiphon pisum was used to compare RNA populations extracted from heads of short- and long-day reared aphids. Microarray analyses revealed that 59 different transcripts were significantly regulated, among which a majority encoded cuticular proteins and several encoded proteins involved in cellular signalling or signal transduction. These results were confirmed by quantitative RT-PCR experiments on two cuticular and three signalling protein genes. Complementary experiments eliminated moulting and circadian rhythms as putative confounding effects. Quantitative RT-PCR performed at additional developmental stages demonstrated the regulation of expression of cuticular and signalling protein genes during the whole process of photoperiod shortening. This suggests that photoperiodic changes could affect cuticle structure and cell to cell communication in the head of aphids in relation with the switch of reproductive modes.

T-DNA transfer from Agrobacterium tumefaciens to the ectomycorrhizal fungus Pisolithus microcarpus.

The model ectomycorrhizal fungus Pisolithus microcarpus isolate 441 was transformed by using Agrobacterium tumefaciens LBA1100 and AGL-1. The selection marker was the Shble gene of Streptoallotecius hidustanus, conferring resistance to phleomycin, under the control of the gpd gene promoter and terminator of Schizophyllum commune. Transformation resulted in phleomycin resistant clones which were confirmed by PCR to contain the resistance cassette. A. tumefaciens-mediated gene transfer would allow the development of RNA interference technology in P. microcarpus.

T-DNA transfer from Agrobacterium tumefaciens to the ectomycorrhizal fungus Pisolithus microcarpus / Transferencia de T-DNA de Agrobacterium tumefaciens al hongo ectomicorrícico Pisolithus microcarpus

The model ectomycorrhizal fungus Pisolithus microcarpus isolate 441 was transformed by using Agrobacterium tumefaciens LBA1100 and AGL-1. The selection marker was the Shble gene of Streptoallotecius hidustanus, conferring resistance to phleomycin, under the control of the gpd gene promoter and terminator of Schizophyllum commune. Transformation resulted in phleomycin resistant clones which were confirmed by PCR to contain the resistance cassette. A. tumefaciens-mediated gene transfer would allow the development of RNA interference technology in P. microcarpus.

El hongo ectomicorrícico modelo Pisolithus microcarpus aislamiento 441 fue transformado utilizando Agrobacterium tumefaciens LBA 1100 y AGL-1. El marcador de selección fue el gen Shble de Streptoallotecius hidustanus, el cual confiere resistencia a fleomicina, bajo el control del promotor y terminador del gen gpd de Schizophyllum commune. La transformación resultó en clones resistentes a fleomicina comprobándose por PCR la presencia del transgen. La transferencia génica mediada por Agrobacterium podría permitir el desarrollo de la tecnología de interferencia por ARN en P. microcarpus.

Genetic analysis of phenotypic variation for ectomycorrhiza formation in an interspecific F1 poplar full-sib family.

A plant's capability to develop ectomycorrhizal symbiosis is under the control of both genetic and environmental factors. In order to determine the roles played by these different factors, we have performed a quantitative genetic analysis of the ability of poplar trees to form ectomycorrhizas. Quantitative genetics were applied to an interspecific family of poplar for which the two parental genetic maps had already been described, and for which data analyses concerning fungal aggressors were obtained. Quantitative trait loci (QTL) related to ectomycorrhiza formation were identified and located in the genetic maps of the two parents. One QTL was located at a linkage group of the genetic map of Populus trichocarpa showing a high concentration of several QTL involved in the pathogenic interaction with the fungus Melampsora larici-populina, the causal agent of leaf rust.

Annotated expressed sequence tags for studies of the regulation of reproductive modes in aphids.

The damaging effect of aphids to crops is largely determined by the spectacular rate of increase of populational expansion due to their parthenogenetic generations. Despite this, the molecular processes triggering the transition between the parthenogenetic and sexual phases between their annual life cycle have received little attention. Here, we describe a collection of genes from the cereal aphid Rhopalosiphum padi expressed during the switch from parthenogenetic to sexual reproduction. After cDNA cloning and sequencing, 726 expressed sequence tags (EST) were annotated. The R. padi EST collection contained a substantial number (139) of bacterial endosymbiont sequences. The majority of R. padi cDNAs encoded either unknown proteins (56%) or housekeeping polypeptides (38%). The large proportion of sequences without similarities in the databases is related to both their small size and their high GC content, corresponding probably to the presence of 5'-unstranslated regions. Fifteen genes involved in developmental and differentiation events were identified by similarity to known genes. Some of these may be useful candidates for markers of the early steps of sexual differentiation.

Hydrophobins in ectomycorrhizas: heterologous transcription of the Pisolithus HydPt-1 gene in yeast and Hebeloma cylindrosporum.

Hydrophobins are fungal cell wall proteins involved in aggregation of hyphae. Upon the development of the ectomycorrhizal symbiosis between tree roots and fungal hyphae, the transcripts of hydrophobin genes markedly accumulated. As the precise role of these proteins in symbiosis is not yet known, we develop heterologous expression system of the Pisolithus hydrophobin HYDPt-1. This gene has been introduced in Saccharomyces cerevisiae and in the ectomycorrhizal basidiomycete Hebeloma cylindrosporum. Introns were required for hydPt-1 transcript accumulation in the basidiomycete H. cylindrosporum. Heterologous transcript accumulation did not alter the phenotype of either species. The lack of altered phenotype resulted from the absence of HYDPt-1 polypeptide accumulation in transformed strains.

Isolation and characterization of a new peroxiredoxin from poplar sieve tubes that uses either glutaredoxin or thioredoxin as a proton donor.

A sequence coding for a peroxiredoxin (Prx) was isolated from a xylem/phloem cDNA library from Populus trichocarpa and subsequently inserted into an expression plasmid yielding the construction pET-Prx. The recombinant protein was produced in Escherichia coli cells and purified to homogeneity with a high yield. The poplar Prx is composed of 162 residues, a property that makes it the shortest plant Prx sequence isolated so far. It was shown that the protein is monomeric and possesses two conserved cysteines (Cys). The Prx degrades hydrogen peroxide and alkyl hydroperoxides in the presence of an exogenous proton donor that can be either thioredoxin or glutaredoxin (Grx). Based on this finding, we propose that the poplar protein represents a new type of Prx that differs from the so-called 2-Cys and 1-Cys Prx, a suggestion supported by the existence of natural fusion sequences constituted of a Prx motif coupled to a Grx motif. The protein was shown to be highly expressed in sieve tubes where thioredoxin h and Grx are also major proteins.

Cloning and expression analysis of a new hydrophobin cDNA from the ectomycorrhizal basidiomycete Pisolithus.

Hydrophobins are fungal cell wall proteins which play a crucial role in cell adhesion and aggregative processes. We have identified a new hydrophobin cDNA (hydPt-3) in the symbiotic mycelium of Pisolithus tinctorius (putative P. albus) during the formation of ectomycorrhizae around eucalypt roots. This sequence is highly divergent from two other previously identified Pisolithus symbiosis-regulated hydrophobins, hydPt-1 and hydPt-2. Also, expression analyses demonstrated that hydPt-3 is up-regulated during the formation of ectomycorrhizae. In contrast to phytopathogenic fungi, changes in glucose or ammonium concentrations in the growth medium did not influence the accumulation of any Pisolithus hydrophobin mRNAs. This suggests that other factors act as regulators of hydrophobin gene expression in ectomycorrhizae.

Immunolocalization of hydrophobin HYDPt-1 from the ectomycorrhizal basidiomycete Pisolithus tinctorius during colonization of Eucalyptus globulus roots.

• The immunolocalization of one of the hydrophobins of Pisolithustinctorius (HYDPt-1) is reported. Hydrophobin proteins play key roles in adhesion and aggregation of fungal hyphae, and it is already known that formation of ectomycorrhizas on eucalypt roots enhances the accumulation of hydrophobin mRNAs in the mycelium of Pisolithus tinctorius. • Purification of SDS-insoluble proteins from the mycelium of P. tinctorius showed the presence of a 13 kDa polypeptide with properties of class I hydrophobin. • Polyconal antibodies were raised against a recombinant HYDPt-1 polypeptide, and these were used for immunofluorescence-coupled transmission electron microscopy. • HYDPt-1 is a cell wall protein located at the surface of the hyphae with no preferential accumulation in the fungal cells of the different tissues of the ectomycorrhiza (i.e. extraradical hyphae, mantle or Hartig net).

The Roles of Auxins and Cytokinins in Mycorrhizal Symbioses.

Most land plant species that have been examined exist naturally with a higher fungus living in and around their roots in a symbiotic partnership called a mycorrhiza. Several types of mycorrhizal symbiosis exist, defined by the host/partner combination and the morphology of the symbiotic structures. The arbuscular mycorrhiza (AM) is ancient and may have co-evolved with land plants. Emerging results from gene expression studies have suggested that subsets of AM genes were co-opted during the evolution of other biotrophic symbioses. Here we compare the roles of phytohormones in AM symbiosis and ectomycorrhizas (EC), a more recent symbiosis. To date, there is little evidence of physiologic overlap between the two symbioses with respect to phytohormone involvement. Research on AM has shown that cytokinin (CK) accumulation is specifically enhanced by symbiosis throughout the plant. We propose a pathway of events linking enhanced CK to development of the AM. Additional and proposed involvement of other phytohormones are also described. The role of auxin in EC symbiosis and recent research advances on the topic are reviewed. We have reflected the literature bias in reporting individual growth regulator effects. However, we consider that gradients and ratios of these molecules are more likely to be the causal agents of morphologic changes resulting from fungal associations. We expect that once the individual roles of these compounds are explained, the subtleties of their function will be more clearly addressed.

A novel class of ectomycorrhiza-regulated cell wall polypeptides in Pisolithus tinctorius.

Development of the ectomycorrhizal symbiosis leads to the aggregation of fungal hyphae to form the mantle. To identify cell surface proteins involved in this developmental step, changes in the biosynthesis of fungal cell wall proteins were examined in Eucalyptus globulus-Pisolithus tinctorius ectomycorrhizas by two-dimensional polyacrylamide gel electrophoresis. Enhanced synthesis of several immunologically related fungal 31- and 32-kDa polypeptides, so-called symbiosis-regulated acidic polypeptides (SRAPs), was observed. Peptide sequences of SRAP32d were obtained after trypsin digestion. These peptides were found in the predicted sequence of six closely related fungal cDNAs coding for ectomycorrhiza up-regulated transcripts. The PtSRAP32 cDNAs represented about 10% of the differentially expressed cDNAs in ectomycorrhiza and are predicted to encode alanine-rich proteins of 28.2 kDa. There are no sequence homologies between SRAPs and previously identified proteins, but they contain the Arg-Gly-Asp (RGD) motif found in cell-adhesion proteins. SRAPs were observed on the hyphal surface by immunoelectron microscopy. They were also found in the host cell wall when P. tinctorius attached to the root surface. RNA blot analysis showed that the steady-state level of PtSRAP32 transcripts exhibited a drastic up-regulation when fungal hyphae form the mantle. These results suggest that SRAPs may form part of a cell-cell adhesion system needed for aggregation of hyphae in ectomycorrhizas.