Large-scale gene discovery in the pea aphid Acyrthosiphon pisum (Hemiptera)

Selenoprofiles: profile-based scanning of eukaryotic genome sequences for selenoprotein genes

Motivation: Selenoproteins are a group of proteins that contain selenocysteine (Sec), a rare amino acid inserted co-translationally into the protein chain. The Sec codon is UGA, which is normally a stop codon. In selenoproteins, UGA is recoded to Sec in presence of specific features on selenoprotein gene transcripts. Due to the dual role of the UGA codon, selenoprotein prediction and annotation are difficult tasks, and even known selenoproteins are often misannotated in genome databases.

Results: We present an homology-based in silico method to scan genomes for members of the known eukaryotic selenoprotein families: selenoprofiles. The core of the method is a set of manually curated highly reliable multiple sequence alignments of selenoprotein families, which are used as queries to scan genomic sequences. Results of the scan are processed through a number of steps, to produce highly accurate predictions of selenoprotein genes with little or no human intervention. Selenoprofiles is a valuable tool for bioinformatic characterization of eukaryotic selenoproteomes, and can complement genome annotation pipelines.

Availability and Implementation: Selenoprofiles is a python-built pipeline that internally runs psitblastn, exonerate, genewise, SECISearch and a number of custom-made scripts and programs. The program is available at http://big.crg.cat/services/selenoprofiles. The predictions presented in this article are available through DAS at http://genome.crg.cat:9000/das/Selenoprofiles_ensembl.

Contact:marco.mariotti@crg.es

Supplementary information:Supplementary data are available at Bioinformatics online.

The gut transcriptome of a gall midge, Mayetiola destructor

The Hessian fly, Mayetiola destructor, is a serious pest of wheat and an experimental organism for the study of gall midge-plant interactions. In addition to food digestion and detoxification, the gut of Hessian fly larvae is also an important interface for insect-host interactions. Analysis of the genes expressed in the Hessian fly larval gut will enhance our understanding of the overall gut physiology and may also lead to the identification of critical molecules for Hessian fly-host plant interactions. Over 10,000 Expressed Sequence Tags (ESTs) were generated and assembled into 2007 clusters. The most striking feature of the Hessian fly larval transcriptome is the existence of a large number of transcripts coding for so-called small secretory proteins (SSP) with amino acids less than 250. Eleven of the 30 largest clusters were SSP transcripts with the largest cluster containing 11.3% of total ESTs. Transcripts coding for diverse digestive enzymes and detoxification proteins were also identified. Putative digestive enzymes included trypsins, chymotrypsins, cysteine proteases, aspartic protease, endo-oligopeptidase, aminopeptidases, carboxypeptidases, and alpha-amylases. Putative detoxification proteins included cytochrome P450s, glutathione S-transferases, peroxidases, ferritins, a catalase, peroxiredoxins, and others. This study represents the first global analysis of gut transcripts from a gall midge. The identification of a large number of transcripts coding for SSPs, digestive enzymes, detoxification proteins in the Hessian fly larval gut provides a foundation for future studies on the functions of these genes.

Combining next-generation sequencing strategies for rapid molecular resource development from an invasive aphid species, Aphis glycines

Background

Aphids are one of the most important insect taxa in terms of ecology, evolutionary biology, genetics and genomics, and interactions with endosymbionts. Additionally, many aphids are serious pest species of agricultural and horticultural plants. Recent genetic and genomic research has expanded molecular resources for many aphid species, including the whole genome sequencing of the pea aphid, Acrythosiphon pisum. However, the invasive soybean aphid, Aphis glycines, lacks in any significant molecular resources.

Methodology/Principal Findings

Two next-generation sequencing technologies (Roche-454 and Illumina GA-II) were used in a combined approach to develop both transcriptomic and genomic resources, including expressed genes and molecular markers. Over 278 million bp were sequenced among the two methods, resulting in 19,293 transcripts and 56,688 genomic sequences. From this data set, 635 SNPs and 1,382 microsatellite markers were identified. For each sequencing method, different soybean aphid biotypes were used which revealed potential biotype specific markers. In addition, we uncovered 39,822 bp of sequence that were related to the obligatory endosymbiont, Buchnera aphidicola, as well as sequences that suggest the presence of Hamiltonella defensa, a facultative endosymbiont.

Conclusions and Significance

Molecular resources for an invasive, non-model aphid species were generated. Additionally, the power of next-generation sequencing to uncover endosymbionts was demonstrated. The resources presented here will complement ongoing molecular studies within the Aphididae, including the pea aphid whole genome, lead to better understanding of aphid adaptation and evolution, and help provide novel targets for soybean aphid control.

Comparative analysis of the Acyrthosiphon pisum genome and expressed sequence tag-based gene sets from other aphid species

To study gene repertoires and their evolution within aphids, we compared the complete genome sequence of Acyrthosiphon pisum (reference gene set) and expressed sequence tag (EST) data from three other species: Myzus persicae, Aphis gossypii and Toxoptera citricida. We assembled ESTs, predicted coding sequences, and identified potential pairs of orthologues (reciprocical best hits) with A. pisum. Pairwise comparisons show that a fraction of the genes evolve fast (high ratio of non-synonymous to synonymous rates), including many genes shared by aphids but with no hit in Uniprot. A detailed phylogenetic study for four fast-evolving genes (C002, JHAMT, Apo and GH) shows that rate accelerations are often associated with duplication events. We also compare compositional patterns between the two tribes of aphids, Aphidini and Macrosiphini.

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.

Aphid wing dimorphisms: linking environmental and genetic control of trait variation

Both genetic and environmental factors underlie phenotypic variation. While research at the interface of evolutionary and developmental biology has made excellent advances in understanding the contribution of genes to morphology, less well understood is the manner in which environmental cues are incorporated during development to influence the phenotype. Also virtually unexplored is how evolutionary transitions between environmental and genetic control of trait variation are achieved. Here, I review investigations into molecular mechanisms underlying phenotypic plasticity in the aphid wing dimorphism system. Among aphids, some species alternate between environmentally sensitive (polyphenic) and genetic (polymorphic) control of wing morph determination in their life cycle. Therefore, a traditional molecular genetic approach into understanding the genetically controlled polymorphism may provide a unique avenue into not only understanding the molecular basis of polyphenic variation in this group, but also the opportunity to compare and contrast the mechanistic basis of environmental and genetic control of similar dimorphisms.

A review of FMRFamide- and RFamide-like peptides in metazoa

Neuropeptides are a diverse class of signalling molecules that are widely employed as neurotransmitters and neuromodulators in animals, both invertebrate and vertebrate. However, despite their fundamental importance to animal physiology and behaviour, they are much less well understood than the small molecule neurotransmitters. The neuropeptides are classified into families according to similarities in their peptide sequence; and on this basis, the FMRFamide and RFamide-like peptides, first discovered in molluscs, are an example of a family that is conserved throughout the animal phyla. In this review, the literature on these neuropeptides has been consolidated with a particular emphasis on allowing a comparison between data sets in phyla as diverse as coelenterates and mammals. The intention is that this focus on the structure and functional aspects of FMRFamide and RFamide-like neuropeptides will inform understanding of conserved principles and distinct properties of signalling across the animal phyla.

Immunity and other defenses in pea aphids, Acyrthosiphon pisum

BACKGROUND

Recent genomic analyses of arthropod defense mechanisms suggest conservation of key elements underlying responses to pathogens, parasites and stresses. At the center of pathogen-induced immune responses are signaling pathways triggered by the recognition of fungal, bacterial and viral signatures. These pathways result in the production of response molecules, such as antimicrobial peptides and lysozymes, which degrade or destroy invaders. Using the recently sequenced genome of the pea aphid (Acyrthosiphon pisum), we conducted the first extensive annotation of the immune and stress gene repertoire of a hemipterous insect, which is phylogenetically distantly related to previously characterized insects models.

RESULTS

Strikingly, pea aphids appear to be missing genes present in insect genomes characterized to date and thought critical for recognition, signaling and killing of microbes. In line with results of gene annotation, experimental analyses designed to characterize immune response through the isolation of RNA transcripts and proteins from immune-challenged pea aphids uncovered few immune-related products. Gene expression studies, however, indicated some expression of immune and stress-related genes.

CONCLUSIONS

The absence of genes suspected to be essential for the insect immune response suggests that the traditional view of insect immunity may not be as broadly applicable as once thought. The limitations of the aphid immune system may be representative of a broad range of insects, or may be aphid specific. We suggest that several aspects of the aphid life style, such as their association with microbial symbionts, could facilitate survival without strong immune protection.

Meeting report: a workshop on Best Practices in Genome Annotation

Efforts to annotate the genomes of a wide variety of model organisms are currently carried out by sequencing centers, model organism databases and academic/institutional laboratories around the world. Different annotation methods and tools have been developed over time to meet the needs of biologists faced with the task of annotating biological data. While standardized methods are essential for consistent curation within each annotation group, methods and tools can differ between groups, especially when the groups are curating different organisms. Biocurators from several institutes met at the Third International Biocuration Conference in Berlin, Germany, April 2009 and hosted the ‘Best Practices in Genome Annotation: Inference from Evidence’ workshop to share their strategies, pipelines, standards and tools. This article documents the material presented in the workshop.

Caterpillars did not evolve from onychophorans by hybridogenesis

The evolution and loss of distinctive larval forms in animal life cycles have produced complex patterns of similarity and difference among life-history stages and major animal lineages. One example of this similarity is the morphological forms of Onychophora (velvet worms) and the caterpillar-like larvae of some insects. Williamson [(2009) Proc Natl Acad Sci USA 106:15786-15790] has made the astonishing and unfounded claim that the ancestors of the velvet worms directly gave rise to insect caterpillars via hybridization and that evidence of this ancient "larval transfer" could be found in comparisons among the genomes of extant onychophorans, insects with larvae, and insects without larvae. Williamson has made a series of predictions arising from his hypothesis and urged genomicists to test them. Here, we use data already in the literature to show these predictions to be false. Hybridogenesis between distantly related animals does not explain patterns of morphological and life-history evolution in general, and the genes and genomes of animals provide strong evidence against hybridization or larval transfer between a velvet worm and an insect in particular.

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-β-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.

Large gene family expansions and adaptive evolution for odorant and gustatory receptors in the pea aphid, Acyrthosiphon pisum

Gaining insight into the mechanisms of chemoreception in aphids is of primary importance for both integrative studies on the evolution of host plant specialization and applied research in pest control management because aphids rely on their sense of smell and taste to locate and assess their host plants. We made use of the recent genome sequence of the pea aphid, Acyrthosiphon pisum, to address the molecular characterization and evolution of key molecular components of chemoreception: the odorant (Or) and gustatory (Gr) receptor genes. We identified 79 Or and 77 Gr genes in the pea aphid genome and showed that most of them are aphid-specific genes that have undergone recent and rapid expansion in the genome. By addressing selection within sets of paralogous Or and Gr expansions, for the first time in an insect species, we show that the most recently duplicated loci have evolved under positive selection, which might be related to the high degree of ecological specialization of this species. Although more functional studies are still needed for insect chemoreceptors, we provide evidence that Grs and Ors have different sets of positively selected sites, suggesting the possibility that these two gene families might have different binding pockets and bind structurally distinct classes of ligand. The pea aphid is the most basal insect species with a completely sequenced genome to date. The identification of chemoreceptor genes in this species is a key step toward further exploring insect comparative genetics, the genomics of ecological specialization and speciation, and new pest control strategies.

Dissimilar molecular defense responses are elicited in Triticum aestivum after infestation by different Diuraphis noxia biotypes

For sustainable solutions to the problem of insect infestation, the study of molecular plant-insect interactions is integral to resistance breeding strategies. This also holds true in the case of wheat (Triticum aestivum), where the Russian wheat aphid (Diuraphis noxia, Kurdjumov, RWA) is responsible for significant crop losses in most major wheat producing countries around the world. Our study is focused on gaining a greater understanding of the resistance mechanisms activated by the RWA resistance gene Dn7 by comparing responses following infestation with three different aphid biotypes (RWA-SA, RWA-US1 and RWA-US2). This consisted of analyzing the resistant wheat line 94M370 (containing Dn7) and its susceptible counterpart (Gamtoos) on a transcriptional level with complementary DNA-amplified fragment length polymorphisms (cDNA-AFLPs) using 17 primer combinations, as well as quantitative reverse transcription polymerase chain reaction (qRT-PCR) of 10 differentially expressed transcripts. The results of this expression profile analysis suggest that Dn7 activates similar responses against the two US aphid biotypes, which differ noticeably from the response following infestation with a South African aphid biotype. This is consistent with recent research showing limited molecular variations between the two US aphid biotypes (approximately 0.12%), compared with a distinctly different South African biotype. We therefore conclude that Dn7 recognizes and interacts in a highly specific manner with different aphid's putative eliciting agents, which in turn activates specific defense pathways unique to that interaction.

WebGMAP: a web service for mapping and aligning cDNA sequences to genomes

The genomes of thousands of organisms are being sequenced, often with accompanying sequences of cDNAs or ESTs. One of the great challenges in bioinformatics is to make these genomic sequences and genome annotations accessible in a user-friendly manner to general biologists to address interesting biological questions. We have created an open-access web service called WebGMAP (http://www.bioinfolab.org/software/webgmap) that seamlessly integrates cDNA-genome alignment tools, such as GMAP, with easy-to-use data visualization and mining tools. This web service is intended to facilitate community efforts in improving genome annotation, determining accurate gene structures and their variations, and exploring important biological processes such as alternative splicing and alternative polyadenylation. For routine sequence analysis, WebGMAP provides a web-based sequence viewer with many useful functions, including nucleotide positioning, six-frame translations, sequence reverse complementation, and imperfect motif detection and alignment. WebGMAP also provides users with the ability to sort, filter and search for individual cDNA sequences and cDNA-genome alignments. Our EST-Genome-Browser can display annotated gene structures and cDNA-genome alignments at scales from 100 to 50 000 nt. With its ability to highlight base differences between query cDNAs and the genome, our EST-Genome-Browser allows biologists to discover potential point or insertion-deletion variations from cDNA-genome alignments.

A water-specific aquaporin involved in aphid osmoregulation

The osmotic pressure of plant phloem sap is generally higher than that of insect body fluids. Water cycling from the distal to proximal regions of the gut is believed to contribute to the osmoregulation of aphids and other phloem-feeding insects, with the high flux of water mediated by a membrane-associated aquaporin. A putative aquaporin referred to as ApAQP1 was identified by RT-PCR of RNA isolated from the guts of pea aphids Acyrthosiphon pisum. The ApAQP1 protein has a predicted molecular mass 28.94kDa. Molecular modeling suggests that ApAQP1 has the general aquaporin topology and possesses the conserved pore properties of water-specific aquaporins. When expressed in Xenopus oocytes, ApAQP1 showed the hallmarks of aquaporin-mediated water transport, including an 18-fold increase in the osmotic water permeability of the oolemma, a reduced activation energy, and inhibition of elevated water transport activity by Hg ions. The ApAQP1 transcript was localised to the stomach and distal intestine, and RNAi-mediated knockdown of its expression resulted in elevated osmotic pressure of the haemolymph. Taken together, these data suggest that ApAQP1 contributes to the molecular basis of water cycling in the aphid gut.

Wounding-mediated gene expression and accelerated viviparous reproduction of the pea aphid Acyrthosiphon pisum

Most insects mount a potent antimicrobial defence upon contact with microbes or microbe-associated pattern molecules. Using a combined set of methods for analysis of insect innate immunity, we report here that piercing of the pea aphid Acyrthosiphon pisum with a bacteria-contaminated needle elicits lysozyme-like activity in the haemolymph but no detectable activities against live bacteria. Confirming these results, we found no homologues of known antimicrobial peptides in our cDNA library generated by using the suppression subtractive hybridization method or in over 90,000 public expressed sequence tag (EST) sequences, but lysozyme genes have recently been described in A. pisum. Interestingly, we discovered that production of viviparous offspring was significantly accelerated upon wounding. Therefore, we postulate that aphids may increase terminal reproductive investment and limit antibacterial defence in response to a threat to their survival.

Whole-mount identification of gene transcripts in aphids: protocols and evaluation of probe accessibility

In situ hybridization has become a powerful tool for detecting the temporal and spatial distribution of gene transcripts in prokaryotes and eukaryotes. We report an efficient protocol for whole-mount identification of the expression of mRNAs in the parthenogenetic pea aphid Acyrthosiphon pisum, an emerging model organism with a growing accumulation of genome sequencing data. In addition to steps common for most animal in situ hybridization protocols, we describe processing methods specific to aphids, the accessibility of antisense riboprobes of different lengths in whole-mounted aphids, and signal intensity versus probe lengths. To find substrate combinations that clearly contrast single and double in situ signals in A. pisum, we tested our protocols using riboprobes constructed from two conserved germline markers, Apvasa and Apnanos, and examined colocalized signals in the germaria and developing oocytes. Finally, we propose conditions for stringent permeabilization that may be applied to tissues deep within the aphid embryo.

ProSOM: core promoter prediction based on unsupervised clustering of DNA physical profiles

MOTIVATION

More and more genomes are being sequenced, and to keep up with the pace of sequencing projects, automated annotation techniques are required. One of the most challenging problems in genome annotation is the identification of the core promoter. Because the identification of the transcription initiation region is such a challenging problem, it is not yet a common practice to integrate transcription start site prediction in genome annotation projects. Nevertheless, better core promoter prediction can improve genome annotation and can be used to guide experimental work.

RESULTS

Comparing the average structural profile based on base stacking energy of transcribed, promoter and intergenic sequences demonstrates that the core promoter has unique features that cannot be found in other sequences. We show that unsupervised clustering by using self-organizing maps can clearly distinguish between the structural profiles of promoter sequences and other genomic sequences. An implementation of this promoter prediction program, called ProSOM, is available and has been compared with the state-of-the-art. We propose an objective, accurate and biologically sound validation scheme for core promoter predictors. ProSOM performs at least as well as the software currently available, but our technique is more balanced in terms of the number of predicted sites and the number of false predictions, resulting in a better all-round performance. Additional tests on the ENCODE regions of the human genome show that 98% of all predictions made by ProSOM can be associated with transcriptionally active regions, which demonstrates the high precision.

AVAILABILITY

Predictions for the human genome, the validation datasets and the program (ProSOM) are available upon request.

Early progress in aphid genomics and consequences for plant-aphid interactions studies

Aphids occupy a niche comprising two conceptual realms: a micron-scale feeding site beneath the plant surface, in which a syringe-like appendage mediates chemical exchange with a specific plant cell type; and the larger realm of a metazoan with sensory organs, a nervous system, and behavior, all responsive to the condition of the host plant and the broader environment. The biology that connects these realms is not well understood, but new details are emerging with the help of genomic tools. The power of these tools is set to increase substantially now that the first genome of an aphid is being sequenced and annotated. This has been possible because a community of aphid researchers focused their efforts to develop and share genomic resources through an international consortium. This complete genome sequence, along with other resources, should permit major advances in understanding the complex and peculiar biological traits responsible for aphids' evolutionary success and their damaging effects on agriculture. This review highlights early progress in the application of aphid genomics and identifies key issues of plant-aphid interactions likely to benefit as molecular tools are further developed. Use of this new knowledge could make significant contributions to crop protection against these and other phloem-feeding insects.

Identification of aphid salivary proteins: a proteomic investigation of Myzus persicae

The role of insect saliva in the first contact between an insect and a plant is crucial during feeding. Some elicitors, particularly in insect regurgitants, have been identified as inducing plant defence reactions. Here, we focused on the salivary proteome of the green peach aphid, Myzus persicae. Proteins were either directly in-solution digested or were separated by 2D SDS-PAGE before trypsin digestion. Resulting peptides were then identified by mass spectrometry coupled with database investigations. A homemade database was constituted of expressed sequence tags from the pea aphid Acyrtosiphon pisum and M. persicae. The databases were used to identify proteins related to M. persicae with a nonsequenced genome. This procedure enabled us to discover glucose oxidase, glucose dehydrogenase, NADH dehydrogenase, alpha-glucosidase and alpha-amylase in M. persicae saliva. The presence of these enzymes is discussed in terms of plant-aphid interactions.

Annotated ESTs from various tissues of the brown planthopper Nilaparvata lugens: a genomic resource for studying agricultural pests

Background

The brown planthopper (BPH), Nilaparvata lugens (Hemiptera, Delphacidae), is a serious insect pests of rice plants. Major means of BPH control are application of agricultural chemicals and cultivation of BPH resistant rice varieties. Nevertheless, BPH strains that are resistant to agricultural chemicals have developed, and BPH strains have appeared that are virulent against the resistant rice varieties. Expressed sequence tag (EST) analysis and related applications are useful to elucidate the mechanisms of resistance and virulence and to reveal physiological aspects of this non-model insect, with its poorly understood genetic background.

Results

More than 37,000 high-quality ESTs, excluding sequences of mitochondrial genome, microbial genomes, and rDNA, have been produced from 18 libraries of various BPH tissues and stages. About 10,200 clusters have been made from whole EST sequences, with average EST size of 627 bp. Among the top ten most abundantly expressed genes, three are unique and show no homology in BLAST searches. The actin gene was highly expressed in BPH, especially in the thorax. Tissue-specifically expressed genes were extracted based on the expression frequency among the libraries. An EST database is available at our web site.

Conclusion

The EST library will provide useful information for transcriptional analyses, proteomic analyses, and gene functional analyses of BPH. Moreover, specific genes for hemimetabolous insects will be identified. The microarray fabricated based on the EST information will be useful for finding genes related to agricultural and biological problems related to this pest.

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.

Analysis of heterochromatic epigenetic markers in the holocentric chromosomes of the aphid Acyrthosiphon pisum

Monomethylated-K9 H3 histones (Me9H3) and heterochromatin protein 1 (HP1) are reported as heterochromatin markers in several eukaryotes possessing monocentric chromosomes. In order to confirm that these epigenetic markers are evolutionarily conserved, we sequenced the HP1 cDNA and verified the distribution of Me9H3 histones and HP1 in the holocentric chromosomes of the aphid Acyrthosiphon pisum. Sequencing indicates that A. pisum HP1 cDNA (called ApHP1) is 1623 bp long, including a 170 bp long 5'UTR and a 688 bp long 3'UTR. The ApHP1 protein consists of 254 amino acidic residues, has a predicted molecular mass of 28 kDa and a net negative charge. At the structural level, it shows an N terminal chromo domain and a chromo shadow domain at the C terminus linked by a short hinge region. At the cytogenetic level, ApHP1 is located exclusively in the heterochromatic regions of the chromosomes. The same heterochromatic regions were labelled after immuno-staining with antibodies against Me9H3 histones, confirming that Hp1 and Me9H3 co-localize at heterochromatic chromosomal areas. Surprisingly, aphid heterochromatin lacks DNA methylation and methylated cytosine residues were mainly spread at euchromatic regions. Finally, the absence of DNA methylation is observed also in aphid rDNA genes that have been repeatedly described as mosaic of methylated and unmethylated units in vertebrates.

Genomic resources for Myzus persicae: EST sequencing, SNP identification, and microarray design

Background

The green peach aphid, Myzus persicae (Sulzer), is a world-wide insect pest capable of infesting more than 40 plant families, including many crop species. However, despite the significant damage inflicted by M. persicae in agricultural systems through direct feeding damage and by its ability to transmit plant viruses, limited genomic information is available for this species.

Results

Sequencing of 16 M. persicae cDNA libraries generated 26,669 expressed sequence tags (ESTs). Aphids for library construction were raised on Arabidopsis thaliana, Nicotiana benthamiana, Brassica oleracea, B. napus, and Physalis floridana (with and without Potato leafroll virus infection). The M. persicae cDNA libraries include ones made from sexual and asexual whole aphids, guts, heads, and salivary glands. In silico comparison of cDNA libraries identified aphid genes with tissue-specific expression patterns, and gene expression that is induced by feeding on Nicotiana benthamiana. Furthermore, 2423 genes that are novel to science and potentially aphid-specific were identified. Comparison of cDNA data from three aphid lineages identified single nucleotide polymorphisms that can be used as genetic markers and, in some cases, may represent functional differences in the protein products. In particular, non-conservative amino acid substitutions in a highly expressed gut protease may be of adaptive significance for M. persicae feeding on different host plants. The Agilent eArray platform was used to design an M. persicae oligonucleotide microarray representing over 10,000 unique genes.

Conclusion

New genomic resources have been developed for M. persicae, an agriculturally important insect pest. These include previously unknown sequence data, a collection of expressed genes, molecular markers, and a DNA microarray that can be used to study aphid gene expression. These resources will help elucidate the adaptations that allow M. persicae to develop compatible interactions with its host plants, complementing ongoing work illuminating plant molecular responses to phloem-feeding insects.

Analysis of genetic variation within clonal lineages of grape phylloxera (Daktulosphaira vitifoliae Fitch) using AFLP fingerprinting and DNA sequencing

Two AFLP fingerprinting methods were employed to estimate the potential of AFLP fingerprints for the detection of genetic diversity within single founder lineages of grape phylloxera (Daktulosphaira vitifoliae Fitch). Eight clonal lineages, reared under controlled conditions in a greenhouse and reproducing asexually throughout a minimum of 15 generations, were monitored and mutations were scored as polymorphisms between the founder individual and individuals of succeeding generations. Genetic variation was detected within all lineages, from early generations on. Six to 15 polymorphic loci (from a total of 141 loci) were detected within the lineages, making up 4.3% of the total amount of genetic variation. The presence of contaminating extra-genomic sequences (e.g., viral material, bacteria, or ingested chloroplast DNA) was excluded as a source of intraclonal variation. Sequencing of 37 selected polymorphic bands confirmed their origin in mostly noncoding regions of the grape phylloxera genome. AFLP techniques were revealed to be powerful for the identification of reproducible banding patterns within clonal lineages.

Initial sequence and comparative analysis of the cat genome

The genome sequence (1.9-fold coverage) of an inbred Abyssinian domestic cat was assembled, mapped, and annotated with a comparative approach that involved cross-reference to annotated genome assemblies of six mammals (human, chimpanzee, mouse, rat, dog, and cow). The results resolved chromosomal positions for 663,480 contigs, 20,285 putative feline gene orthologs, and 133,499 conserved sequence blocks (CSBs). Additional annotated features include repetitive elements, endogenous retroviral sequences, nuclear mitochondrial (numt) sequences, micro-RNAs, and evolutionary breakpoints that suggest historic balancing of translocation and inversion incidences in distinct mammalian lineages. Large numbers of single nucleotide polymorphisms (SNPs), deletion insertion polymorphisms (DIPs), and short tandem repeats (STRs), suitable for linkage or association studies were characterized in the context of long stretches of chromosome homozygosity. In spite of the light coverage capturing approximately 65% of euchromatin sequence from the cat genome, these comparative insights shed new light on the tempo and mode of gene/genome evolution in mammals, promise several research applications for the cat, and also illustrate that a comparative approach using more deeply covered mammals provides an informative, preliminary annotation of a light (1.9-fold) coverage mammal genome sequence.

Strong heterogeneity in nucleotidic composition and codon bias in the pea aphid (Acyrthosiphon pisum) shown by EST-based coding genome reconstruction

The aim of this study was to analyze patterns of nucleotidic composition and codon usage in the pea aphid genome (Acyrthosiphon pisum). A collection of 60,000 expressed sequence tags (ESTs) in the pea aphid has been used to automatically reconstruct 5809 coding sequences (CDSs), based on similarity with known proteins and on coding style recognition. Reconstructions were manually checked for ribosomal proteins, leading to tentatively reconstruct the nea-complete set of this category. Pea aphid coding sequences showed a shift toward AT (especially at the third codon position) compared to drosophila homologues. Genes with a putative high level of expression (ribosomal and other genes with high EST support) remained more GC3-rich and had a distinct codon usage from bulk sequences: they exhibited a preference for C-ending codons and CGT (for arginine), which thus appeared optimal for translation. However, the discrimination was not as strong as in drosophila, suggesting a reduced degree of translational selection. The space of variation in codon usage for A. pisum appeared to be larger than in drosophila, with a substantial fraction of genes that remained GC3-rich. Some of those (in particular some structural proteins) also showed high levels of codon bias and a very strong preference for C-ending codons, which could be explained either by strong translational selection or by other mechanisms. Finally, genomic traces were analyzed to build 206 fragments containing a full CDS, which allowed studying the correlations between GC contents of coding and those of noncoding (flanking and introns) sequences.

Cross-species transferability of microsatellite markers from six aphid (Hemiptera: Aphididae) species and their use for evaluating biotypic diversity in two cereal aphids

The abundance and distribution of microsatellites, or simple sequence repeats (SSRs) were explored in the expressed sequence tag (EST) and genomic sequences of the pea aphid, Acyrthosiphon pisum (Harris), and the green peach aphid, Myzus persicae (Sulzer). A total of 108 newly developed, together with 40 published, SSR markers were investigated for their cross-species transferability among six aphid species. Genetic diversity among six greenbug, Schizaphis graminum (Rondani) and two Russian wheat aphid, Diuraphis noxia (Kurdjumov) biotypes was further examined with 67 transferable SSRs. It was found that the pea aphid genome is abundant in SSRs with a unique frequency and distribution of SSR motifs. Cross-species transferability of EST-derived SSRs is dependent on phylogenetic closeness between SSR donor and target species, but is higher than that of genomic SSRs. Neighbor-joining analysis of SSR data revealed host-adapted genetic divergence as well as regional differentiation of greenbug biotypes. The two Russian wheat aphid biotypes are genetically as diverse as the greenbug ones although it was introduced into the USA only 20 years ago. This is the first report of large-scale development of SSR markers in aphids, which are expected to have wide applications in aphid genetic, ecological and evolutionary studies.

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.

Molecular profiling of proteolytic and lectin transcripts in Homalodisca vitripennis (Hemiptera: Auchenorrhyncha: Cicadellidae) feeding on sunflower and cowpea

Homalodisca vitripennis Germar 1821 (Hemiptera: Cicadellidae) [Takiya et al. (2006) Ann Entomol Soc Am 99:648-655; syn. H. coagulata (Say)] salivary gland and gut EST libraries were used to isolate cDNA fragments of the gene transcripts encoding for cathepsin L, asparaginyl endopeptidase, cathepsin B, metalloendopeptidase, cathepsin D, multicatalytic endopeptidase, and a sugar-binding C-type lectin. Transcript levels were evaluated in immature and adult H. vitripennis feeding on sunflower (Helianthus annuus) or cowpea (Vigna unguiculata). Northern blot hybridization results showed that expression of most of the transcripts were similar for all developmental stages and feeding on the two diets examined. However, the expression of the transcript for asparaginyl endopeptidase was less expressed in sunflower-fed adult females compared to sunflower-fed immatures. Also, the expression of the C-type lectin transcript was up-regulated in adults compared to immatures when fed on either diet. Documenting both the presence and variation of transcript expression involved in putative digestive proteolysis in this xylem-feeding leafhopper is noteworthy and aids efforts to design specific diet formulations for mass production of hosts and parasitoids to be used as effective biological control measures.

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.

Promoting transcriptome diversity

Although the number of protein-encoding human genes is more limited than many had estimated, the human transcript repertoire is much more diverse than anticipated. In part, transcript diversity is generated through the use of alternative promoters and alternate splicing. In addition, based on discoveries using technologies such as full-length cDNA libraries and whole genome tiling microarrays, it is now likely that non-protein-encoding transcripts comprise a substantial fraction of the human RNA population. Much attention is currently focused on understanding the role of alternative promoters in generating transcript diversity, both for non-protein-encoding (ncRNAs) and protein-encoding RNAs.

Common genome-wide patterns of transcript accumulation underlying the wing polyphenism and polymorphism in the pea aphid (Acyrthosiphon pisum)

The pea aphid, Acyrthosiphon pisum, exhibits several environmentally cued polyphenisms, in which discrete, alternative phenotypes are produced. At low-density, parthenogenetic females produce unwinged female progeny, but at high-density females produce progeny that develop with wings. These alternative phenotypes represent a solution to the competing demands of dispersal and reproduction. Males also develop as either winged or unwinged, but these alternatives are determined by a genetic polymorphism. Winged and unwinged males are morphologically less distinct from each other than winged and unwinged females, possibly because males experience fewer trade-offs between dispersal and reproduction. To assess whether shared physiological differences mirror the shared morphological differences that characterize the wing polyphenism and polymorphism, we used a cDNA microarray representing an estimated 10% of the coding genome (1734 genes) to examine differential transcript accumulation between winged and unwinged females and males. We identified several transcripts that differentially accumulate between winged and unwinged morphs in both sexes, the majority of which are involved in energy production. Unexpectedly, the extent of differential transcript accumulation between winged and unwinged morphs was greater for adult males than for adult females. Together, these results suggest not only that similar physiological differences underlie the polyphenism and polymorphism, but that male morphs, like females, are subject to trade-offs between reproduction and dispersal that are reflected in levels of transcript accumulation and possibly genome-wide patterns of gene regulation. These data also provide a baseline for future studies of the molecular and physiological basis of life-history trade-offs.

Proteomic profiling of aphid Macrosiphum euphorbiae responses to host-plant-mediated stress induced by defoliation and water deficit

Abiotic and biotic host-plant stress, such as desiccation and herbivory, may strongly affect sap-sucking insects such as aphids via changes in plant chemicals of insect nutritional or plant defensive value. Here, we examined (i) water deprivation and (ii) defoliation by the beetle Leptinotarsa decemlineata as stresses indirectly affecting the aphid Macrosiphum euphorbiae via its host plant Solanum tuberosum. For plant-induced stress, aphids were reared on healthy vs. continuously stressed potato for 14 days (no watering; defoliation maintained at approximately 40%). Aphid performance under stress was correlated with metabolic responses monitored by profiling of the aphid proteome. M. euphorbiae was strongly affected by water stress, as adult survival, total aphid number and biomass were reduced by 67%, 64%, and 79%, respectively. Aphids performed normally on defoliated potato, indicating that they were unaffected or able to compensate any stress induced by plant defoliation. Stressed aphid proteomes revealed 419-453 protein spots, including 27 that were modulated specifically or jointly under each kind of host-plant stress. Reduced aphid fitness on water-stressed plants mostly correlated with modulation of proteins involved in energy metabolism, apparently to conserve energy in order to prioritize survival. Despite normal performance, several aphid proteins that are known to be implicated in cell communication were modulated on defoliated plants, possibly suggesting modified aphid behaviour. The GroEL protein (or symbionin) of the endosymbiont Buchnera aphidicola was predominant under all conditions in M. euphorbiae. Its expression level was not significantly affected by aphid host-plant stresses, which is consistent with the high priority of symbiosis in stressed aphids.

Characterization of a membrane-bound aminopeptidase purified from Acyrthosiphon pisum midgut cells. A major binding site for toxic mannose lectins

A single membrane-bound aminopeptidase N (APN) occurs in the pea aphid (Acyrthosiphon pisum Harris) midgut, with a pH optimum of 7.0, pI of 8.1 and molecular mass of 130 kDa. This enzyme accounts for more than 15.6% of the total gut proteins. After being solubilized in detergent, APN was purified to homogeneity. The enzyme is a glycoprotein rich in mannose residues, which binds the entomotoxic lectins of the concanavalin family. The internal sequence of APN is homologous with a conservative domain in APNs, and degenerated primers of highly conserved APN motifs were used to screen a gut cDNA library. The complete sequence of APN has standard residues involved in zinc co-ordination and catalysis and a glycosyl-phosphatidylinositol anchor, as in APNs from Lepidoptera. APN has a broad specificity towards N-terminal amino acids, but does not hydrolyze acidic aminoacyl-peptides, thus resembling the mammalian enzyme (EC 3.4.11.2). The kcat/Km ratios for different di-, tri-, tetra-, and penta-peptides suggest a preference for tripeptides, and that subsites S1, S2' and S3' are pockets able to bind bulky aminoacyl residues. Bestatin and amastatin bound APN in a rapidly reversible mode, with Ki values of 1.8 microM and 0.6 microM, respectively. EDTA inactivates this APN (k(obs) 0.14 M(-1) x s(-1), reaction order of 0.44) at a rate that is reduced by competitive inhibitors. In addition to oligopeptide digestion, APN is proposed to be associated with amino-acid-absorption processes which, in contrast with aminopeptidase activity, may be hampered on lectin binding.

AphidBase: a database for aphid genomic resources

AphidBase aims to (i) store recently acquired genomic resources on aphids and (ii) compare them to other insect resources as functional annotation tools. For that, the Drosophila melanogaster genome has been loaded in the database using the GMOD open source software for a comparison with the 17 069 pea aphid unique transcripts (contigs) and the 13 639 gene transcripts of the Anopheles gambiae. Links to FlyBase and A.gambiae Entrez databases allow a rapid characterization of the putative functions of the aphid sequences. Text mining of the D.melanogaster literature was performed to construct a network of co-cited gene or protein names, which should facilitate functional annotation of aphid homolog sequences. AphidBase represents one of the first genomic databases for a hemipteran insect.

AVAILABILITY

http://w3.rennes.inra.fr/AphidBase.

Codon usage bias and tRNA over-expression in Buchnera aphidicola after aromatic amino acid nutritional stress on its host Acyrthosiphon pisum

Codon usage bias and relative abundances of tRNA isoacceptors were analysed in the obligate intracellular symbiotic bacterium, Buchnera aphidicola from the aphid Acyrthosiphon pisum, using a dedicated 35mer oligonucleotide microarray. Buchnera is archetypal of organisms living with minimal metabolic requirements and presents a reduced genome with high-evolutionary rate. Codonusage in Buchnera has been overcome by the high mutational bias towards AT bases. However, several lines of evidence for codon usage selection are given here. A significant correlation was found between tRNA relative abundances and codon composition of Buchnera genes. A significant codon usage bias was found for the choice of rare codons in Buchnera: C-ending codons are preferred in highly expressed genes, whereas G-ending codons are avoided. This bias is not explained by GC skew in the bacteria and might correspond to a selection for perfect matching between codon-anticodon pairs for some essential amino acids in Buchnera proteins. Nutritional stress applied to the aphid host induced a significant overexpression of most of the tRNA isoacceptors in bacteria. Although, molecular regulation of the tRNA operons in Buchnera was not investigated, a correlation between relative expression levels and organization in transcription unit was found in the genome of Buchnera.

A dual-genome microarray for the pea aphid, Acyrthosiphon pisum, and its obligate bacterial symbiont, Buchnera aphidicola

Background

The best studied insect-symbiont system is that of aphids and their primary bacterial endosymbiont Buchnera aphidicola. Buchnera inhabits specialized host cells called bacteriocytes, provides nutrients to the aphid and has co-speciated with its aphid hosts for the past 150 million years. We have used a single microarray to examine gene expression in the pea aphid, Acyrthosiphon pisum, and its resident Buchnera. Very little is known of gene expression in aphids, few studies have examined gene expression in Buchnera, and no study has examined simultaneously the expression profiles of a host and its symbiont. Expression profiling of aphids, in studies such as this, will be critical for assigning newly discovered A. pisum genes to functional roles. In particular, because aphids possess many genes that are absent from Drosophila and other holometabolous insect taxa, aphid genome annotation efforts cannot rely entirely on homology to the best-studied insect systems. Development of this dual-genome array represents a first attempt to characterize gene expression in this emerging model system.

Results

We chose to examine heat shock response because it has been well characterized both in Buchnera and in other insect species. Our results from the Buchnera of A. pisum show responses for the same gene set as an earlier study of heat shock response in Buchnera for the host aphid Schizaphis graminum. Additionally, analyses of aphid transcripts showed the expected response for homologs of known heat shock genes as well as responses for several genes with unknown functional roles.

Conclusion

We examined gene expression under heat shock of an insect and its bacterial symbiont in a single assay using a dual-genome microarray. Further, our results indicate that microarrays are a useful tool for inferring functional roles of genes in A. pisum and other insects and suggest that the pea aphid genome may contain many gene paralogs that are differentially regulated.

RNAi knockdown of a salivary transcript leading to lethality in the pea aphid, Acyrthosiphon pisum

Abstract Injection of siRNA (small interfering RNA) into parthenogenetic adult pea aphids (Acyrthosiphon pisum) is shown here to lead to depletion of a target salivary gland transcript. The siRNA was generated from double stranded RNA that covered most of the open reading frame of the transcript, which we have called Coo2. The Coo2 transcript level decreases dramatically over a 3-day period after injection of siRNA. With a lag of 1 to 2 days, the siCoo2-RNA injected insects died, on average 8 days before the death of control insects injected with siRNA for green fluorescent protein. It appears, therefore, that siRNA injections into adults will be a useful tool in studying the roles of individual transcripts in aphid salivary glands and suggests that siCoo2-RNA injections can be a useful positive control in such studies.