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

RNAi of Neuropeptide CCHamide-1 and Its Receptor Indicates Role in Feeding Behavior in the Pea Aphid, Acyrthosiphon pisum

Neuropeptide CCHamide-1 (abbreviated as CCHa1) is a recently discovered peptide that is present in many arthropods and is the ligand of the CCHa1R, a member of the G protein-coupled receptors (GPCRs) superfamily, which plays a regulatory role in diverse physiological processes such as feeding, circadian rhythm, insulin production, lipid metabolism, growth, and reproduction. However, the function of this gene in aphids is still unknown. Here, we characterized and determined the potential role of CCHa1/CCHa1R signaling in the pea aphid, Acyrthosiphon pisum, which is a notorious pest in agriculture. The docking analysis revealed that the CCHa1 peptide binds to its receptor CCHa1R through specific amino acid residues, which are critical for maintaining the structural and functional integrity of the peptide-receptor complex. Quantitative real-time reverse transcription-PCR (qRT-PCR) revealed the expression levels of CCHa1/CCHa1R transcripts in different development stages and different tissues, indicating that the CCHa1 expression was high in the first nymphal instar compared to the upcoming nymphal instars and adults, and was predominantly high in the brain. The CCHa1/CCHa1R transcript levels were significantly upregulated in starved aphids compared to fed aphids. Moreover, RNAi knockdown by the injection of dsRNA-CCHa1 and dsRNA-CCHa1R significantly reduced the corresponding expression of the target gene and reduced their food intake in adult aphids, as revealed by the electrical penetration graph results. CCHa1/CCHa1R-silencing also reduced the reproduction, but not the survival, in A. pisum. Our data demonstrated that CCHa1/CCHa1R play a role in the regulation of feeding in A. pisum, suggesting a role of the CCHa1 signaling pathway in the aphids relating to their nutritional status.

The ATP-binding cassette transporter subfamily G member 4 mediates cuticular hydrocarbon transport to regulate drought tolerance in Acyrthosiphon pisum

ABC transporters are a highly conserved membrane protein class that promote the transport of substances across membranes. Under drought conditions, insects primarily regulate the content of cuticular hydrocarbon (CHC) to retain water and prevent evaporative loss. Involvement of ABC transporter protein G (ABCG) subfamily genes in insect CHC transport has been relatively understudied. In this study, we demonstrated that ABCG4 gene in Acyrthosiphon pisum (ApABCG4) is involved in CHC transport and affects drought tolerance by regulating CHC accumulation. ApABCG4 is strongly expressed in the abdominal cuticle and embryonic stages of A. pisum. Effective silencing of ApABCG4 was achieved using RNAi, and the silencing duration was analyzed. ApABCG4 silencing resulted in a significant decrease in the total and component contents of the CHC and cuticular waxy coatings of A. pisum. Nevertheless, the internal hydrocarbon content remained unchanged. The lack of cuticular hydrocarbons significantly reduced the drought tolerance of A. pisum, shortening its survival time under drought stress. Drought stress caused significant upregulation of ApABCG4. Molecular docking showed that ApABCG4 has a high binding affinity for nine n-alkanes of CHC through electrostatic interactions. These results indicate that ApABCG4 is a novel RNAi target with key applications in aphid biological control.

Antioxidant Enzyme, Transcriptomic, and Metabolomic Changes in Lily (Lilium spp.) Leaves Induced by Aphis gossypii Glover

Cotton aphids (Aphis gossypii Glover) cause harm by feeding on phloem sap and spreading plant viruses to lily. Understanding the mechanisms by which aphids infest lily plants is crucial for effective aphid management and control. In this study, we investigated the activity of antioxidants, integrated nontargeted metabolomes and transcriptomes of lilies infested by cotton aphids to explore the changes in lily leaves. Overall, the results indicated that the catalase (CAT) activity in the leaves of the lily plants was greater than that in the leaves of the control plants. A comprehensive identification of 604 substances was conducted in the leaves. Furthermore, the differentially abundant metabolite analysis revealed the enrichment of phenylalanine metabolism and α-linolenic acid metabolism. Moreover, 3574 differentially expressed genes (DEGs), whose expression tended to increase, were linked to glutathione metabolism and phenylpropanoid biosynthesis. In addition, the integrated analysis revealed that the defensive response of lily leaves to aphids is manifested through antioxidant reactions, phenylpropane and flavonoid biosynthesis, and α-linolenic acid metabolism. Finally, the key metabolites were CAT, glutathione, coumaric acid, and jasmonic acid, along with the key genes chalcone synthase (CHS), phenylalanine ammonia-lyase (PAL), and 12-oxo-phytodienoic acid reductase (OPR). Accordingly, the findings of this research elucidate the molecular and metabolic reactions of A. gossypii in lily plants, offering valuable insights for developing aphid resistance strategies in lily farming.

A new gene finding tool GeneMark-ETP significantly improves the accuracy of automatic annotation of large eukaryotic genomes

Large-scale genomic initiatives, such as the Earth BioGenome Project, require efficient methods for eukaryotic genome annotation. Here we present an automatic gene finder, GeneMark-ETP, integrating genomic-, transcriptomic- and protein-derived evidence that has been developed with a focus on large plant and animal genomes. GeneMark-ETP first identifies genomic loci where extrinsic data is sufficient for making gene predictions with 'high confidence'. The genes situated in the genomic space between the high confidence genes are predicted in the next stage. The set of high confidence genes serves as an initial training set for the statistical model. Further on, the model parameters are iteratively updated in the rounds of gene prediction and parameter re-estimation. Upon reaching convergence, GeneMark-ETP makes the final predictions and delivers the whole complement of predicted genes. GeneMark-ETP outperformed gene finders using a single type of extrinsic evidence. Comparisons with gene finders utilizing both transcript- and protein-derived extrinsic evidence, MAKER2, and TSEBRA, demonstrated that GeneMark-ETP delivered state-of-the-art gene prediction accuracy with the margin of outperforming existing approaches increasing in its applications to larger and more complex eukaryotic genomes.

Systematic assessment of long-read RNA-seq methods for transcript identification and quantification

The Long-read RNA-Seq Genome Annotation Assessment Project (LRGASP) Consortium was formed to evaluate the effectiveness of long-read approaches for transcriptome analysis. The consortium generated over 427 million long-read sequences from cDNA and direct RNA datasets, encompassing human, mouse, and manatee species, using different protocols and sequencing platforms. These data were utilized by developers to address challenges in transcript isoform detection and quantification, as well as de novo transcript isoform identification. The study revealed that libraries with longer, more accurate sequences produce more accurate transcripts than those with increased read depth, whereas greater read depth improved quantification accuracy. In well-annotated genomes, tools based on reference sequences demonstrated the best performance. When aiming to detect rare and novel transcripts or when using reference-free approaches, incorporating additional orthogonal data and replicate samples are advised. This collaborative study offers a benchmark for current practices and provides direction for future method development in transcriptome analysis.

Biochemical probing of phloem sap defensive traits in Brassica juncea-B. fruticulosa introgression lines following Lipaphis erysimi infestation

The lack of resistance to Lipaphis erysimi in cultivated Brassicas makes caused this pest highly devastating resulting in significant loss of rapeseed-mustard productivity in India. B. fruticulosa, a wild crucifer is known as an excellent source of resistance to L. erysimi. Therefore, we planned to assess defense associated biochemical alterations and molecular components of B. juncea-B. fruticulosa ILs to mustard aphid. Phenotypic assessment of ILs on the basis of aphid population per plant (APP) categorized genotypes into resistant (7.15-18.50 APP), moderately susceptible (42.29-53.33 APP) and susceptible (70.00-77.07 APP) genotypes. Mustard aphid infested minimally B. fruticulosa (0.80 APP) among tested genotypes. The maximum increase in catalase (CAT) activity was determined in B. fruticulosa and resistant ILs after 48 h (2.03 and 1.76-fold, respectively) and one week (2.98 and 1.79-fold, respectively) of mustard aphid infestation. The strong induction of CAT2 transcripts (19.25-fold) and CAT activity (5.88-fold) along with low aphid count in resistant IL, Ad4-64 (13.85 APP) suggested the pivotal role of CAT in resistance to mustard aphid. Guaiacol peroxidase (GPX) was significantly decreased following pest infestation at both infestation stages. The ascorbate content was highest in resistant IL, ADV-6RD (2.14-fold) after one week of aphid infestation. H2O2 content rapidly increased in B. juncea-B. fruticulosa derived lines after 48 h of aphid infestation. The negative and significant association between APP and CAT (- 0.56** and - 0.48*, respectively), glutathione (- 0.43* and - 0.40*, respectively), H2O2 (- 0.57** and - 0.43*, respectively) at both 48 h and one week infestation stages signified their role in deterring mustard aphid infestation. The positive and significant association between total sugars (0.33* at 7 DPI), reducing sugars (0.33* at 7 DPI), sucrose (0.36** at 48 h) and APP indicated that higher the sugars content, higher will be mustard aphid infestation in B. juncea derived ILs. The information being generated and key candidates (CAT2, ascorbate and H2O2) being identified may help in effective deployment of B. fruticulosa resistance in mustard breeding.

Lipophorin receptor regulates the cuticular hydrocarbon accumulation and adult fecundity of the pea aphid Acyrthosiphon pisum

Cuticular hydrocarbons form a barrier that protects terrestrial insects from water loss via the epicuticle. Lipophorin loads and transports lipids, including hydrocarbons, from one tissue to another. In some insects, the lipophorin receptor (LpR), which binds to lipophorin and accepts its lipid cargo, is essential for female fecundity because it mediates the incorporation of lipophorin by developing oocytes. However, it is unclear whether LpR is involved in the accumulation of cuticular hydrocarbons and its precise role in aphid reproduction remains unknown. We herein present the results of our molecular characterization, phylogenetic analysis, and functional annotation of the pea aphid (Acyrthosiphon pisum) LpR gene (ApLpR). This gene was transcribed throughout the A. pisum life cycle, but especially during the embryonic stage and in the abdominal cuticle. Furthermore, we optimized the RHA interference (RNAi) parameters by determining the ideal dose and duration for gene silencing in the pea aphid. We observed that the RNAi-based ApLpR suppression significantly decreased the internal and cuticular hydrocarbon contents as well as adult fecundity. Additionally, a deficiency in cuticular hydrocarbons increased the susceptibility of aphids to desiccation stress, with decreased survival rates under simulated drought conditions. Moreover, ApLpR expression levels significantly increased in response to the desiccation treatment. These results confirm that ApLpR is involved in transporting hydrocarbons and protecting aphids from desiccation stress. Furthermore, this gene is vital for aphid reproduction. Therefore, the ApLpR gene of A. pisum may be a novel RNAi target relevant for insect pest management.

Nitrogen Supply and Host-Plant Genotype Modulate the Transcriptomic Profile of Plasmodiophora brassicae

Nitrogen fertilization can affect the susceptibility of Brassica napus to the telluric pathogen Plasmodiophora brassicae. Our previous works highlighted that the influence of nitrogen can strongly vary regarding plant cultivar/pathogen strain combinations, but the underlying mechanisms are unknown. The present work aims to explore how nitrogen supply can affect the molecular physiology of P. brassicae through its life epidemiological cycle. A time-course transcriptome experiment was conducted to study the interaction, under two conditions of nitrogen supply, between isolate eH and two B. napus genotypes (Yudal and HD-018), harboring (or not harboring) low nitrogen-conditional resistance toward this isolate (respectively). P. brassicae transcriptional patterns were modulated by nitrogen supply, these modulations being dependent on both host-plant genotype and kinetic time. Functional analysis allowed the identification of P. brassicae genes expressed during the secondary phase of infection, which may play a role in the reduction of Yudal disease symptoms in low-nitrogen conditions. Candidate genes included pathogenicity-related genes ("NUDIX," "carboxypeptidase," and "NEP-proteins") and genes associated to obligate biotrophic functions of P. brassicae. This work illustrates the importance of considering pathogen's physiological responses to get a better understanding of the influence of abiotic factors on clubroot resistance/susceptibility.

pH effect on strain-specific transcriptomes of the take-all fungus

The soilborne fungus Gaeumannomyces tritici (G. tritici) causes the take-all disease on wheat roots. Ambient pH has been shown to be critical in different steps of G. tritici life cycle such as survival in bulk soil, saprophytic growth, and pathogenicity on plants. There are however intra-specific variations and we previously found two types of G. tritici strains that grow preferentially either at acidic pH or at neutral/alkaline pH; gene expression involved in pH-signal transduction pathway and pathogenesis was differentially regulated in two strains representative of these types. To go deeper in the description of the genetic pathways and the understanding of this adaptative mechanism, transcriptome sequencing was achieved on two strains (PG6 and PG38) which displayed opposite growth profiles in two pH conditions (acidic and neutral). PG6, growing better at acidic pH, overexpressed in this condition genes related to cell proliferation. In contrast, PG38, which grew better at neutral pH, overexpressed in this condition genes involved in fatty acids and amino acid metabolisms, and genes potentially related to pathogenesis. This strain also expressed stress resistance mechanisms at both pH, to assert a convenient growth under various ambient pH conditions. These differences in metabolic pathway expression between strains at different pH might buffer the effect of field or soil variation in wheat fields, and explain the success of the pathogen.

Identification and Characterization of Neuropeptides and Their G Protein-Coupled Receptors (GPCRs) in the Cowpea Aphid Aphis craccivora

Neuropeptides are the most abundant and diverse signal molecules in insects. They act as neurohormones and neuromodulators to regulate the physiology and behavior of insects. The majority of neuropeptides initiate downstream signaling pathways through binding to G protein-coupled receptors (GPCRs) on the cell surface. In this study, RNA-seq technology and bioinformatics were used to search for genes encoding neuropeptides and their GPCRs in the cowpea aphid Aphis craccivora. And the expression of these genes at different developmental stages of A. craccivora was analyzed by quantitative real-time PCR (qRT-PCR). A total of 40 candidate genes encoding neuropeptide precursors were identified from the transcriptome data, which is roughly equivalent to the number of neuropeptide genes that have been reported in other insects. On this basis, software analysis combined with homologous prediction estimated that there could be more than 60 mature neuropeptides with biological activity. In addition, 46 neuropeptide GPCRs were obtained, of which 40 belong to rhodopsin-like receptors (A-family GPCRs), including 21 families of neuropeptide receptors and 7 orphan receptors, and 6 belong to secretin-like receptors (B-family GPCRs), including receptors for diuretic hormone 31, diuretic hormone 44 and pigment-dispersing factor (PDF). Compared with holometabolous insects such as Drosophila melanogaster, the coding genes for sulfakinin, corazonin, arginine vasopressin-like peptide (AVLP), and trissin and the corresponding receptors were not found in A. craccivora. It is speculated that A. craccivora likely lacks the above neuropeptide signaling pathways, which is consistent with Acyrthosiphon pisum and that the loss of these pathways may be a common feature of aphids. In addition, expression profiling revealed neuropeptide genes and their GPCR genes that are differentially expressed at different developmental stages and in different wing morphs. This study will help to deepen our understanding of the neuropeptide signaling systems in aphids, thus laying the foundation for the development of new methods for aphid control targeting these signaling systems.

(E)-β-farnesene synthase gene affects aphid behavior in transgenic Medicago sativa

BACKGROUND

Pea aphid (Acyrthosiphon pisum Harris) is one of the major pests in alfalfa crops, causing significant yield losses. (E)-β-farnesene (EβF), an alarm pheromone released by pea aphid, is generic to many species of aphids, and is used to minimize potential danger from predators and parasitoids by avoiding the source of the pheromone.

RESULTS

In this study, EβF synthase gene was constructed into a plant expression vector, and overexpressed in alfalfa (Medicago sativa L.), with expression among transgenic lines confirmed by qRT-PCR. Subcellular localization analysis showed that EβF synthase gene was expressed in the plasma membrane and nucleus of the leaf. GC/MS of extraction from transgenic alfalfa indicated emission of EβF ranging from 5.92 to 13.09 ng day^-1 g^-1 fresh tissue. Behavior assays in Y-olfactometers demonstrated that transgenic alfalfa expressing AaEβF gene could repel pea aphids, with aphids taking a significantly longer time to select a transgenic line compared with the control line (P < 0.01).

CONCLUSION

We have demonstrated a potentially valuable strategy of using EβF synthase genes for aphid control in alfalfa. © 2018 Society of Chemical Industry.

Solid-phase microextraction-based cuticular hydrocarbon profiling for intraspecific delimitation in Acyrthosiphon pisum

Insect cuticular hydrocarbons (CHCs) play critical roles in reducing water loss and chemical communication. Species-specific CHC profiles have been used increasingly as an excellent character for species classification. However, considerably less is known about their potential for population delimitation within species. The aims of this study were to develop a solid-phase microextraction (SPME)-based CHC collection method and to investigate whether CHC profiles could serve as potential chemotaxonomic tools for intraspecific delimitation in Acyrthosiphon pisum. Optimization of fibers for SPME sampling revealed that 7 μm polydimethylsiloxane (PDMS) demonstrated the most efficient adsorption of CHCs among five different tested fibers. SPME sampling showed good reproducibility with repeated collections of CHCs from a single aphid. Validation of SPME was performed by comparing CHC profiles with those from conventional hexane extractions. The two methods showed no qualitative differences in CHCs, although SPME appeared to extract relatively fewer short-chained CHCs. While CHC profiles of a given population differed among developmental stages, wing dimorphism types, and host plants, wingless adult aphids showed very low variance in relative proportions of individual CHC components. Reproducibility of CHC profiles was explored further to classify wingless adult morphs of A. pisum from five different geographic regions that showed no variation in mitochondrial COI gene sequences. Our results demonstrate that CHC profiles are useful in intraspecific delimitation in the field of insect chemotaxonomy.

Cytochrome P450 gene, CYP4G51, modulates hydrocarbon production in the pea aphid, Acyrthosiphon pisum

Terrestrial insects deposit a layer of hydrocarbons (HCs) as waterproofing agents on their epicuticle. The insect-specific CYP4G genes, subfamily members of P450, have been found in all insects with sequenced genomes to date. They are critical for HC biosynthesis in Drosophila; however, their functional roles in other insects including the piercing-sucking hemipterous aphids remain unknown. In this study, we presented the molecular characterization and a functional study of the CYP4G51 gene in the pea aphid, Acyrthosiphon pisum (Harris). CYP4G51 transcript was detectable across the whole life cycle of A. pisum, and was prominently expressed in the aphid head and abdominal cuticle. Up-regulation of CYP4G51 under desiccation stress was more significant in the third instar nymphs compared with the adults. Also, up-regulation of CYP4G51 was observed when the aphids fed on an artificial diet compared with those fed on the broad bean plant, and was positively correlated with a high level of cuticular HCs (CHCs). RNAi knockdown of CYP4G51 significantly reduced its expression and caused reductions in both internal and external HCs. A deficiency in CHCs resulted in aphids being more susceptible to desiccation, with increased mortality under desiccation stress. The current results confirm that CYP4G51 modulates HC biosynthesis to protect aphids from desiccation. Moreover, our data also indicate that saturated and straight-chain HCs play a major role in cuticular waterproofing in the pea aphid. A. pisum CYP4G51 could be considered as a novel RNAi target in the field of insect pest management.

Molecular characterization of two isoforms of a farnesyl pyrophosphate synthase gene in wheat and their roles in sesquiterpene synthesis and inducible defence against aphid infestation

Aphids are important pests of wheat (Triticum aestivum) that affect crop production globally. Herbivore-induced emission of sesquiterpenes can repel pests, and farnesyl pyrophosphate synthase (FPS) is a key enzyme involved in sesquiterpene biosynthesis. However, fps orthologues in wheat and their functional roles in sesquiterpene synthesis and defence against aphid infestation are unknown. Here, two fps isoforms, Tafps1 and Tafps2, were identified in wheat. Quantitative real-time polymerase chain reaction (qRT-PCR) and in vitro catalytic activity analyses were conducted to investigate expression patterns and activity. Heterologous expression of these isoforms in Arabidopsis thaliana, virus-induced gene silencing (VIGS) in wheat and aphid behavioural assays were performed to understand the functional roles of these two isoforms. We demonstrated that Tafps1 and Tafps2 played different roles in induced responses to aphid infestation and in sesquiterpene synthesis. Heterologous expression in A. thaliana resulted in repulsion of the peach aphid (Myzus persicae). Wheat plants with these two isoforms transiently silenced were significantly attractive to grain aphid (Sitobion avenae). Our results provide new insights into induced defence against aphid herbivory in wheat, in particular, the different roles of the two Tafps isoforms in both sesquiterpene biosynthesis and defence against aphid infestation.

Gene expression analysis of parthenogenetic embryonic development of the pea aphid, Acyrthosiphon pisum, suggests that aphid parthenogenesis evolved from meiotic oogenesis

Aphids exhibit a form of phenotypic plasticity, called polyphenism, in which genetically identical females reproduce sexually during one part of the life cycle and asexually (via parthenogenesis) during the remainder of the life cycle. The molecular basis for aphid parthenogenesis is unknown. Cytological observations of aphid parthenogenesis suggest that asexual oogenesis evolved either through a modification of meiosis or from a mitotic process. As a test of these alternatives, we assessed the expression levels and expression patterns of canonical meiotic recombination and germline genes in the sexual and asexual ovaries of the pea aphid, Acyrthosiphon pisum. We observed expression of all meiosis genes in similar patterns in asexual and sexual ovaries, with the exception that some genes encoding Argonaute-family members were not expressed in sexual ovaries. In addition, we observed that asexual aphid tissues accumulated unspliced transcripts of Spo11, whereas sexual aphid tissues accumulated primarily spliced transcripts. In situ hybridization revealed Spo11 transcript in sexual germ cells and undetectable levels of Spo11 transcript in asexual germ cells. We also found that an obligately asexual strain of pea aphid produced little spliced Spo11 transcript. Together, these results suggest that parthenogenetic oogenesis evolved from a meiosis-like, and not a mitosis-like, process and that the aphid reproductive polyphenism may involve a modification of Spo11 gene activity.

Engineering plants for aphid resistance: current status and future perspectives

The current status of development of transgenic plants for improved aphid resistance, and the pros and cons of different strategies are reviewed and future perspectives are proposed. Aphids are major agricultural pests that cause significant yield losses of crop plants each year. Excessive dependence on insecticides for aphid control is undesirable because of the development of insecticide resistance, the potential negative effects on non-target organisms and environmental pollution. Transgenic plants engineered for resistance to aphids via a non-toxic mode of action could be an efficient alternative strategy. In this review, the distribution of major aphid species and their damages on crop plants, the so far isolated aphid-resistance genes and their applications in developments of transgenic plants for improved aphid resistance, and the pros and cons of these strategies are reviewed and future perspectives are proposed. Although the transgenic plants developed through expressing aphid-resistant genes, manipulating plant secondary metabolism and plant-mediated RNAi strategy have been demonstrated to confer improved aphid resistance to some degree. So far, no aphid-resistant transgenic crop plants have ever been commercialized. This commentary is intended to be a helpful insight into the generation and future commercialization of aphid-resistant transgenic crops in a global context.

Advances in Human Biology: Combining Genetics and Molecular Biophysics to Pave the Way for Personalized Diagnostics and Medicine

Advances in several biology-oriented initiatives such as genome sequencing and structural genomics, along with the progress made through traditional biological and biochemical research, have opened up a unique opportunity to better understand the molecular effects of human diseases. Human DNA can vary significantly from person to person and determines an individual’s physical characteristics and their susceptibility to diseases. Armed with an individual’s DNA sequence, researchers and physicians can check for defects known to be associated with certain diseases by utilizing various databases. However, for unclassified DNA mutations or in order to reveal molecular mechanism behind the effects, the mutations have to be mapped onto the corresponding networks and macromolecular structures and then analyzed to reveal their effect on the wild type properties of biological processes involved. Predicting the effect of DNA mutations on individual’s health is typically referred to as personalized or companion diagnostics. Furthermore, once the molecular mechanism of the mutations is revealed, the patient should be given drugs which are the most appropriate for the individual genome, referred to as pharmacogenomics. Altogether, the shift in focus in medicine towards more genomic-oriented practices is the foundation of personalized medicine. The progress made in these rapidly developing fields is outlined.

SIBIS: a Bayesian model for inconsistent protein sequence estimation

MOTIVATION

The prediction of protein coding genes is a major challenge that depends on the quality of genome sequencing, the accuracy of the model used to elucidate the exonic structure of the genes and the complexity of the gene splicing process leading to different protein variants. As a consequence, today's protein databases contain a huge amount of inconsistency, due to both natural variants and sequence prediction errors.

RESULTS

We have developed a new method, called SIBIS, to detect such inconsistencies based on the evolutionary information in multiple sequence alignments. A Bayesian framework, combined with Dirichlet mixture models, is used to estimate the probability of observing specific amino acids and to detect inconsistent or erroneous sequence segments. We evaluated the performance of SIBIS on a reference set of protein sequences with experimentally validated errors and showed that the sensitivity is significantly higher than previous methods, with only a small loss of specificity. We also assessed a large set of human sequences from the UniProt database and found evidence of inconsistency in 48% of the previously uncharacterized sequences. We conclude that the integration of quality control methods like SIBIS in automatic analysis pipelines will be critical for the robust inference of structural, functional and phylogenetic information from these sequences.

AVAILABILITY AND IMPLEMENTATION

Source code, implemented in C on a linux system, and the datasets of protein sequences are freely available for download at http://www.lbgi.fr/∼julie/SIBIS.

Plant-generated artificial small RNAs mediated aphid resistance

BACKGROUND

RNA silencing is an important mechanism for regulation of endogenous gene expression and defense against genomic intruders in plants. This natural defense system was adopted to generate virus-resistant plants even before the mechanism of RNA silencing was unveiled. With the clarification of that mechanism, transgenic antiviral plants were developed that expressed artificial virus-specific hairpin RNAs (hpRNAs) or microRNAs (amiRNAs) in host plants. Previous works also showed that plant-mediated RNA silencing technology could be a practical method for constructing insect-resistant plants by expressing hpRNAs targeting essential genes of insects.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we chose aphid Myzus persicae of order Hemiptera as a target insect. To screen for aphid genes vulnerable to attack by plant-mediated RNA silencing to establish plant aphid resistance, we selected nine genes of M. persicae as silencing targets, and constructed their hpRNA-expressing vectors. For the acetylcholinesterase 2 coding gene (MpAChE2), two amiRNA-expressing vectors were also constructed. The vectors were transformed into tobacco plants (Nicotiana tabacum cv. Xanti). Insect challenge assays showed that most of the transgenic plants gained aphid resistance, among which those expressing hpRNAs targeting V-type proton ATPase subunit E-like (V-ATPaseE) or tubulin folding cofactor D (TBCD) genes displayed stronger aphicidal activity. The transgenic plants expressing amiRNAs targeting two different sites in the MpAChE2 gene exhibited better aphid resistance than the plants expressing MpAChE2-specific hpRNA.

CONCLUSIONS/SIGNIFICANCE

Our results indicated that plant-mediated insect-RNA silencing might be an effective way to develop plants resistant to insects with piercing-sucking mouthparts, and both the selection of vulnerable target genes and the biogenetic type of the small RNAs were crucial for the effectiveness of aphid control. The expression of insect-specific amiRNA is a promising and preferable approach to engineer plants resistant to aphids and, possibly, to other plant-infesting insects.

Aphids evolved novel secreted proteins for symbiosis with bacterial endosymbiont

Aphids evolved novel cells, called bacteriocytes, that differentiate specifically to harbour the obligatory mutualistic endosymbiotic bacteria Buchnera aphidicola. The genome of the host aphid Acyrthosiphon pisum contains many orphan genes that display no similarity with genes found in other sequenced organisms, prompting us to hypothesize that some of these orphan genes are related to lineage-specific traits, such as symbiosis. We conducted deep sequencing of bacteriocytes mRNA followed by whole mount in situ hybridizations of over-represented transcripts encoding aphid-specific orphan proteins. We identified a novel class of genes that encode small proteins with signal peptides, which are often cysteine-rich, that are over-represented in bacteriocytes. These genes are first expressed at a developmental time point coincident with the incorporation of symbionts strictly in the cells that contribute to the bacteriocyte and this bacteriocyte-specific expression is maintained throughout the aphid's life. The expression pattern suggests that recently evolved secretion proteins act within bacteriocytes, perhaps to mediate the symbiosis with beneficial bacterial partners, which is reminiscent of the evolution of novel cysteine-rich secreted proteins of leguminous plants that regulate nitrogen-fixing endosymbionts.

Evidence of a conserved functional role for DNA methylation in termites

Many organisms are capable of developing distinct phenotypes from the same genotype. This developmental plasticity is particularly prevalent in insects, which can produce alternate adaptive forms in response to distinct environmental cues. The ability to develop divergent phenotypes from the same genotype often relies on epigenetic information, which affects gene function and is transmitted through cell divisions. One of the most important epigenetic marks, DNA methylation, has been lost in several insect lineages, yet its taxonomic distribution and functional conservation remain uninvestigated in many taxa. In the present study, we demonstrate that the signature of high levels of DNA methylation exists in the expressed genes of two termites, Reticulitermes flavipes and Coptotermes formosanus. Further, we show that DNA methylation is preferentially targeted to genes with ubiquitous expression among morphs. Functional associations of DNA methylation are also similar to those observed in other invertebrate taxa with functional DNA methylation systems. Finally, we demonstrate an association between DNA methylation and the long-term evolutionary conservation of genes. Overall, our findings strongly suggest DNA methylation is present at particularly high levels in termites and may play similar roles to those found in other insects.

Transcriptome of Dickeya dadantii infecting Acyrthosiphon pisum reveals a strong defense against antimicrobial peptides

The plant pathogenic bacterium Dickeya dadantii has recently been shown to be able to kill the aphid Acyrthosiphon pisum. While the factors required to cause plant disease are now well characterized, those required for insect pathogeny remain mostly unknown. To identify these factors, we analyzed the transcriptome of the bacteria isolated from infected aphids. More than 150 genes were upregulated and 300 downregulated more than 5-fold at 3 days post infection. No homologue to known toxin genes could be identified in the upregulated genes. The upregulated genes reflect the response of the bacteria to the conditions encountered inside aphids. While only a few genes involved in the response to oxidative stress were induced, a strong defense against antimicrobial peptides (AMP) was induced. Expression of a great number of efflux proteins and transporters was increased. Besides the genes involved in LPS modification by addition of 4-aminoarabinose (the arnBCADTEF operon) and phosphoethanolamine (pmrC, eptB) usually induced in Gram negative bacteria in response to AMPs, dltBAC and pbpG genes, which confer Gram positive bacteria resistance to AMPs by adding alanine to teichoic acids, were also induced. Both types of modification confer D. dadantii resistance to the AMP polymyxin. A. pisum harbors symbiotic bacteria and it is thought that it has a very limited immune system to maintain these populations and do not synthesize AMPs. The arnB mutant was less pathogenic to A. pisum, which suggests that, in contrast to what has been supposed, aphids do synthesize AMP.

Evolution of the lamin protein family: what introns can tell

Lamins are the major components of the nuclear lamina, a filamentous layer found at the interphase between chromatin and the inner nuclear membrane. The lamina supports the nuclear envelope and provides anchorage sites for chromatin. Lamins and their associated proteins are required for most nuclear activities, mitosis, and for linking the nucleoskeleton to the network of cytoskeletal filaments. Mutations in lamins and their associated proteins give rise to a wide range of diseases, collectively called laminopathies. This review focuses on the evolution of the lamin protein family. Evolution from basal metazoans to man will be described on the basis of protein sequence comparisons and analyses of their gene structure. Lamins are the founding members of the family of intermediate filament proteins. How genes encoding cytoplasmic IF proteins could have arisen from the archetypal lamin gene progenitor, can be inferred from a comparison of the respective gene structures. The lamin/IF protein family seems to be restricted to the metazoans. In general, invertebrate genomes harbor only a single lamin gene encoding a B-type lamin. The archetypal lamin gene structure found in basal metazoans is conserved up to the vertebrate lineage. The completely different structure of lamin genes in Caenorhabditis and Drosophila are exceptions rather than the rule within their systematic groups. However, variation in the length of the coiled-coil forming central domain might be more common than previously anticipated. The increase in the number of lamin genes in vertebrates can be explained by two rounds of genome duplication. The origin of lamin A by exon shuffling might explain the processing of prelamin A to the mature non-isoprenylated form of lamin A. By alternative splicing the number of vertebrate lamin proteins has increased even further. Lamin C, an alternative splice form of the LMNA gene, is restricted to mammals. Amphibians and mammals express germline-specific lamins that differ in their protein structure from that of somatic lamins. Evidence is provided that there exist lamin-like proteins outside the metazoan lineage.

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.

Genome scans reveal candidate regions involved in the adaptation to host plant in the pea aphid complex

A major goal in evolutionary biology is to uncover the genetic basis of adaptation. Divergent selection exerted on ecological traits may result in adaptive population differentiation and reproductive isolation and affect differentially the level of genetic divergence along the genome. Genome-wide scan of large sets of individuals from multiple populations is a powerful approach to identify loci or genomic regions under ecologically divergent selection. Here, we focused on the pea aphid, a species complex of divergent host races, to explore the organization of the genomic divergence associated with host plant adaptation and ecological speciation. We analysed 390 microsatellite markers located at variable distances from predicted genes in replicate samples of sympatric populations of the pea aphid collected on alfalfa, red clover and pea, which correspond to three common host-adapted races reported in this species complex. Using a method that accounts for the hierarchical structure of our data set, we found a set of 11 outlier loci that show higher genetic differentiation between host races than expected under the null hypothesis of neutral evolution. Two of the outliers are close to olfactory receptor genes and three other nearby genes encoding salivary proteins. The remaining outliers are located in regions with genes of unknown functions, or which functions are unlikely to be involved in interactions with the host plant. This study reveals genetic signatures of divergent selection across the genome and provides an inventory of candidate genes responsible for plant specialization in the pea aphid, thereby setting the stage for future functional studies.

Localizing F(ST) outliers on a QTL map reveals evidence for large genomic regions of reduced gene exchange during speciation-with-gene-flow

Populations that maintain phenotypic divergence in sympatry typically show a mosaic pattern of genomic divergence, requiring a corresponding mosaic of genomic isolation (reduced gene flow). However, mechanisms that could produce the genomic isolation required for divergence-with-gene-flow have barely been explored, apart from the traditional localized effects of selection and reduced recombination near centromeres or inversions. By localizing F(ST) outliers from a genome scan of wild pea aphid host races on a Quantitative Trait Locus (QTL) map of key traits, we test the hypothesis that between-population recombination and gene exchange are reduced over large 'divergence hitchhiking' (DH) regions. As expected under divergence hitchhiking, our map confirms that QTL and divergent markers cluster together in multiple large genomic regions. Under divergence hitchhiking, the nonoutlier markers within these regions should show signs of reduced gene exchange relative to nonoutlier markers in genomic regions where ongoing gene flow is expected. We use this predicted difference among nonoutliers to perform a critical test of divergence hitchhiking. Results show that nonoutlier markers within clusters of F(ST) outliers and QTL resolve the genetic population structure of the two host races nearly as well as the outliers themselves, while nonoutliers outside DH regions reveal no population structure, as expected if they experience more gene flow. These results provide clear evidence for divergence hitchhiking, a mechanism that may dramatically facilitate the process of speciation-with-gene-flow. They also show the power of integrating genome scans with genetic analyses of the phenotypic traits involved in local adaptation and population divergence.

Host generated siRNAs attenuate expression of serine protease gene in Myzus persicae

Background

Sap sucking hemipteran aphids damage diverse crop species. Although delivery of ds-RNA or siRNA through microinjection/feeding has been demonstrated, the efficacy of host-mediated delivery of aphid-specific dsRNA in developing aphid resistance has been far from being elucidated.

Methodology/Principal Findings

Transgenic Arabidopsis expressing ds-RNA of Myzus persicae serine protease (MySP) was developed that triggered the generation of corresponding siRNAs amenable for delivery to the feeding aphids. M. persicae when fed on the transgenic plants for different time intervals under controlled growth conditions resulted in a significant attenuation of the expression of MySP and a commensurate decline in gut protease activity. Although the survivability of these aphids was not affected, there was a noticeable decline in their fecundity resulting in a significant reduction in parthenogenetic population.

Conclusions/Significance

The study highlighted the feasibility of developing host based RNAi-mediated resistance against hemipteran pest aphids.

Metabolic engineering of plant-derived (E)-β-farnesene synthase genes for a novel type of aphid-resistant genetically modified crop plants

Aphids are major agricultural pests that cause significant yield losses of crop plants each year. Excessive dependence on insecticides for long-term aphid control is undesirable because of the development of insecticide resistance, the potential negative effects on non-target organisms and environmental pollution. Transgenic crops engineered for resistance to aphids via a non-toxic mode of action could be an efficient alternative strategy. (E)-β-Farnesene (EβF) synthases catalyze the formation of EβF, which for many pest aphids is the main component of the alarm pheromone involved in the chemical communication within these species. EβF can also be synthesized by certain plants but is then normally contaminated with inhibitory compounds. Engineering of crop plants capable of synthesizing and emitting EβF could cause repulsion of aphids and also the attraction of natural enemies that use EβF as a foraging cue, thus minimizing aphid infestation. In this review, the effects of aphids on host plants, plants' defenses against aphid herbivory and the recruitment of natural enemies for aphid control in an agricultural setting are briefly introduced. Furthermore, the plant-derived EβF synthase genes cloned to date along with their potential roles in generating novel aphid resistance via genetically modified approaches are discussed.

Isolation and characterization of polymorphic microsatellite markers in Tetranychus urticae and cross amplification in other Tetranychidae and Phytoseiidae species of economic importance

Tetranychus urticae Koch is a cosmopolitan phytophagous mite considered as the most polyphagous species among spider mites. Population genetic studies using molecular markers such as microsatellites have proven to be extremely informative to address questions about population structure, phylogeography and host preferences. The aim of this study was to increase the available molecular tools to gain insight into the genetic structure of T. urticae populations of citrus orchards, which might help in their management. Five microsatellite DNA libraries were developed using probes with the motifs CT, CTT, GT and CAC following the FIASCO protocol. Positive clones, those that included the insert with the microsatellite, were detected using the PIMA-PCR technique. Combinations of primers were designed on 22 out of 32 new microsatellites loci and their polymorphism was tested in four populations sampled along the eastern coast of Spain. Eleven successful amplifications were obtained. Cross amplification was tested in the tetranychids Aphlonobia histricina, Eutetranychus banksi, E. orientalis, Oligonychus perseae, Panonychus citri, Tetranychus evansi, T. okinawanus and T. turkestani, and the phytoseiids Amblyseius swirskii, A. cucumeris, A. andersoni, Euseius stipulatus, Neoseiulus barkeri, N. californicus, Phytoseiulus persimilis and Typhlodromus phialatus. Eight successful cross amplifications were obtained.

Comparative proteomics reveals a significant bias toward alternative protein isoforms with conserved structure and function

Advances in high-throughput mass spectrometry are making proteomics an increasingly important tool in genome annotation projects. Peptides detected in mass spectrometry experiments can be used to validate gene models and verify the translation of putative coding sequences (CDSs). Here, we have identified peptides that cover 35% of the genes annotated by the GENCODE consortium for the human genome as part of a comprehensive analysis of experimental spectra from two large publicly available mass spectrometry databases. We detected the translation to protein of “novel” and “putative” protein-coding transcripts as well as transcripts annotated as pseudogenes and nonsense-mediated decay targets.

We provide a detailed overview of the population of alternatively spliced protein isoforms that are detectable by peptide identification methods. We found that 150 genes expressed multiple alternative protein isoforms. This constitutes the largest set of reliably confirmed alternatively spliced proteins yet discovered. Three groups of genes were highly overrepresented. We detected alternative isoforms for 10 of the 25 possible heterogeneous nuclear ribonucleoproteins, proteins with a key role in the splicing process. Alternative isoforms generated from interchangeable homologous exons and from short indels were also significantly enriched, both in human experiments and in parallel analyses of mouse and Drosophila proteomics experiments. Our results show that a surprisingly high proportion (almost 25%) of the detected alternative isoforms are only subtly different from their constitutive counterparts.

Many of the alternative splicing events that give rise to these alternative isoforms are conserved in mouse. It was striking that very few of these conserved splicing events broke Pfam functional domains or would damage globular protein structures. This evidence of a strong bias toward subtle differences in CDS and likely conserved cellular function and structure is remarkable and strongly suggests that the translation of alternative transcripts may be subject to selective constraints.

Aphids: a model for polyphenism and epigenetics

Environmental conditions can alter the form, function, and behavior of organisms over short and long timescales, and even over generations. Aphid females respond to specific environmental cues by transmitting signals that have the effect of altering the development of their offspring. These epigenetic phenomena have positioned aphids as a model for the study of phenotypic plasticity. The molecular basis for this epigenetic inheritance in aphids and how this type of inheritance system could have evolved are still unanswered questions. With the availability of the pea aphid genome sequence, new genomics technologies, and ongoing genomics projects in aphids, these questions can now be addressed. Here, we review epigenetic phenomena in aphids and recent progress toward elucidating the molecular basis of epigenetics in aphids. The discovery of a functional DNA methylation system, functional small RNA system, and expanded set of chromatin modifying genes provides a platform for analyzing these pathways in the context of aphid plasticity. With these tools and further research, aphids are an emerging model system for studying the molecular epigenetics of polyphenisms.

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.

Acyrthosiphon pisum AQP2: a multifunctional insect aquaglyceroporin

Annotation of the recently sequenced genome of the pea aphid (Acyrthosiphon pisum) identified a gene ApAQP2 (ACYPI009194, Gene ID: 100168499) with homology to the Major Intrinsic Protein/aquaporin superfamily of membrane channel proteins. Phylogenetic analysis suggests that ApAQP2 is a member of an insect-specific clade of this superfamily. Homology model structures of ApAQP2 showed a novel array of amino acids comprising the substrate selectivity-determining "aromatic/arginine" region of the putative transport pore. Subsequent characterization of the transport properties of ApAQP2 upon expression in Xenopus oocytes supports an unusual substrate selectivity profile. Water permeability analyses show that the ApAQP2 protein exhibits a robust mercury-insensitive aquaporin activity. However unlike the water-specific ApAQP1 protein, ApAQP2 forms a multifunctional transport channel that shows a wide permeability profile to a range of linear polyols, including the potentially biologically relevant substrates glycerol, mannitol and sorbitol. Gene expression analysis indicates that ApAQP2 is highly expressed in the insect bacteriocytes (cells bearing the symbiotic bacteria Buchnera) and the fat body. Overall the results demonstrate that ApAQP2 is a novel insect aquaglyceroporin which may be involved in water and polyol transport in support of the Buchnera symbiosis and aphid osmoregulation.

The biocontrol bacterium Pseudomonas fluorescens Pf29Arp strain affects the pathogenesis-related gene expression of the take-all fungus Gaeumannomyces graminis var. tritici on wheat roots

The main effects of antagonistic rhizobacteria on plant pathogenic fungi are antibiosis, fungistasis or an indirect constraint through the induction of a plant defence response. To explore different biocontrol mechanisms, an in vitro confrontation assay was conducted with the rhizobacterium Pseudomonas fluorescens Pf29Arp as a biocontrol agent of the fungus Gaeumannomyces graminis var. tritici (Ggt) on wheat roots. In parallel with the assessment of disease extension, together with the bacterial and fungal root colonization rates, the transcript levels of candidate fungal pathogenicity and plant-induced genes were monitored during the 10-day infection process. The bacterial inoculation of wheat roots with the Pf29Arp strain reduced the development of Ggt-induced disease expressed as attack frequency and necrosis length. The growth rates of Ggt and Pf29Arp, monitored through quantitative polymerase chain reaction of DNA amounts with a part of the Ggt 18S rDNA gene and a specific Pf29Arp strain detection probe, respectively, increased throughout the interactions. Bacterial antagonism and colonization had no significant effect on root colonization by Ggt. The expression of fungal and plant genes was quantified in planta by quantitative reverse transcription-polymerase chain reaction during the interactions thanks to the design of specific primers and an innovative universal reference system. During the early stages of the tripartite interaction, several of the fungal genes assayed were down-regulated by Pf29Arp, including two laccases, a β-1,3-exoglucanase and a mitogen-activated protein kinase. The plant host glutathione-S-transferase gene was induced by Ggt alone and up-regulated by Pf29Arp bacteria in interaction with the pathogen. We conclude that Pf29Arp antagonism acts through the alteration of fungal pathogenesis and probably through the activation of host defences.

The origins, evolution, and functional potential of alternative splicing in vertebrates

Alternative splicing (AS) has the potential to greatly expand the functional repertoire of mammalian transcriptomes. However, few variant transcripts have been characterized functionally, making it difficult to assess the contribution of AS to the generation of phenotypic complexity and to study the evolution of splicing patterns. We have compared the AS of 309 protein-coding genes in the human ENCODE pilot regions against their mouse orthologs in unprecedented detail, utilizing traditional transcriptomic and RNAseq data. The conservation status of every transcript has been investigated, and each functionally categorized as coding (separated into coding sequence [CDS] or nonsense-mediated decay [NMD] linked) or noncoding. In total, 36.7% of human and 19.3% of mouse coding transcripts are species specific, and we observe a 3.6 times excess of human NMD transcripts compared with mouse; in contrast to previous studies, the majority of species-specific AS is unlinked to transposable elements. We observe one conserved CDS variant and one conserved NMD variant per 2.3 and 11.4 genes, respectively. Subsequently, we identify and characterize equivalent AS patterns for 22.9% of these CDS or NMD-linked events in nonmammalian vertebrate genomes, and our data indicate that functional NMD-linked AS is more widespread and ancient than previously thought. Furthermore, although we observe an association between conserved AS and elevated sequence conservation, as previously reported, we emphasize that 30% of conserved AS exons display sequence conservation below the average score for constitutive exons. In conclusion, we demonstrate the value of detailed comparative annotation in generating a comprehensive set of AS transcripts, increasing our understanding of AS evolution in vertebrates. Our data supports a model whereby the acquisition of functional AS has occurred throughout vertebrate evolution and is considered alongside amino acid change as a key mechanism in gene evolution.

Genomics of environmentally induced phenotypes in 2 extremely plastic arthropods

Understanding how genes and the environment interact to shape phenotypes is of fundamental importance for resolving important issues in adaptive evolution. Yet, for most model species with mature genetics and accessible genomic resources, we know little about the natural environmental factors that shape their evolution. By contrast, animal species with deeply understood ecologies and well characterized responses to environmental cues are rarely subjects of genomic investigations. Here, we preview advances in genomics in aphids and waterfleas that may help transform research on the regulatory mechanisms of phenotypic plasticity. This insect and crustacean duo has the capacity to produce extremely divergent phenotypes in response to environmental stimuli. Sexual fate and reproductive mode are condition-dependent in both groups, which are also capable of altering morphology, physiology and behavior in response to biotic and abiotic cues. Recently, the genome sequences for the pea aphid Acyrthosiphon pisum and the waterflea Daphnia pulex were described by their respective research communities. We propose that an integrative study of genome biology focused on the condition-dependent transcriptional basis of their shared plastic traits and specialized mode of reproduction will provide broad insight into adaptive plasticity and genome by environment interactions. We highlight recent advances in understanding the genome regulation of alternative phenotypes and environmental cue processing, and we propose future research avenues to discover gene networks and epigenetic mechanisms underlying phenotypic plasticity.

Analysis of the workers head transcriptome of the Asian subterranean termite, Coptotermes gestroi

The lower termite, Coptotermes gestroi (Isoptera: Rhinotermitidae), is originally from Southeast Asia and has become a pest in Brazil. The main goal of this study was to survey C. gestroi transcriptome composition. To accomplish this, we sequenced and analyzed 3003 expressed sequence tags (ESTs) isolated from libraries of worker heads. After assembly, 695 uniESTs were obtained from which 349 have similarity with known sequences. Comparison with insect genomes demonstrated similarity, primarily with genes from Apis mellifera (28%), Tribolium castaneum (28%) and Aedes aegypti (10%). Notably, we identified two endogenous cellulases in the sequences, which may be of interest for biotechnological applications. The results presented in this work represent the first genomic study of the Asian subterranean termite, Coptotermes gestroi.

Responses of the pea aphid transcriptome to infection by facultative symbionts

Serratia symbiotica is a facultative symbiont of pea aphids (Acyrthosiphon pisum) that provides tolerance to heat stress. Although the phenotypic effects of facultative symbionts upon hosts have been studied in some detail, little is known about the molecular and genomic basis of these interactions. Previous studies show a large impact of S. symbiotica upon the aphid metabolome. Whole-genome transcriptional profiling and next-generation sequencing demonstrated expression of 94% of RefSeq genes from the pea aphid genome, providing the largest dataset to date on aphid gene expression. However, only 28 genes showed changes in expression with S. symbiotica infection, and these changes were of small magnitude. No expression differences in genes involved in innate immunity in other insects were observed. Therefore, the large metabolic impact of S. symbiotica is most likely a result of metabolism of the symbiont itself, or of post-transcriptional modification of host gene expression. Although S. symbiotica has a major influence on its host's metabolome and resistance to heat, it induces little change in gene expression in its host.

Shotgun proteomics aids discovery of novel protein-coding genes, alternative splicing, and "resurrected" pseudogenes in the mouse genome

Recent advances in proteomic mass spectrometry (MS) offer the chance to marry high-throughput peptide sequencing to transcript models, allowing the validation, refinement, and identification of new protein-coding loci. We present a novel pipeline that integrates highly sensitive and statistically robust peptide spectrum matching with genome-wide protein-coding predictions to perform large-scale gene validation and discovery in the mouse genome for the first time. In searching an excess of 10 million spectra, we have been able to validate 32%, 17%, and 7% of all protein-coding genes, exons, and splice boundaries, respectively. Moreover, we present strong evidence for the identification of multiple alternatively spliced translations from 53 genes and have uncovered 10 entirely novel protein-coding genes, which are not covered in any mouse annotation data sources. One such novel protein-coding gene is a fusion protein that spans the Ins2 and Igf2 loci to produce a transcript encoding the insulin II and the insulin-like growth factor 2-derived peptides. We also report nine processed pseudogenes that have unique peptide hits, demonstrating, for the first time, that they are not just transcribed but are translated and are therefore resurrected into new coding loci. This work not only highlights an important utility for MS data in genome annotation but also provides unique insights into the gene structure and propagation in the mouse genome. All these data have been subsequently used to improve the publicly available mouse annotation available in both the Vega and Ensembl genome browsers (http://vega.sanger.ac.uk).

RNAcode: robust discrimination of coding and noncoding regions in comparative sequence data

With the availability of genome-wide transcription data and massive comparative sequencing, the discrimination of coding from noncoding RNAs and the assessment of coding potential in evolutionarily conserved regions arose as a core analysis task. Here we present RNAcode, a program to detect coding regions in multiple sequence alignments that is optimized for emerging applications not covered by current protein gene-finding software. Our algorithm combines information from nucleotide substitution and gap patterns in a unified framework and also deals with real-life issues such as alignment and sequencing errors. It uses an explicit statistical model with no machine learning component and can therefore be applied "out of the box," without any training, to data from all domains of life. We describe the RNAcode method and apply it in combination with mass spectrometry experiments to predict and confirm seven novel short peptides in Escherichia coli and to analyze the coding potential of RNAs previously annotated as "noncoding." RNAcode is open source software and available for all major platforms at http://wash.github.com/rnacode.

Genomic revelations of a mutualism: the pea aphid and its obligate bacterial symbiont

The symbiosis of the pea aphid Acyrthosphion pisum with the bacterium Buchnera aphidicola APS represents the best-studied insect obligate symbiosis. Here we present a refined picture of this symbiosis by linking pre-genomic observations to new genomic data that includes the complete genomes of the eukaryotic and prokaryotic symbiotic partners. In doing so, we address four issues central to understanding the patterns and processes operating at the A. pisum/Buchnera APS interface. These four issues include: (1) lateral gene transfer, (2) host immunity, (3) symbiotic metabolism, and (4) regulation.

Identification, cloning, and expression of a GHF9 cellulase from Tribolium castaneum (Coleoptera: Tenebrionidae)

The availability of sequenced insect genomes has allowed for discovery and functional characterization of novel genes and proteins. We report use of the Tribolium castaneum (Herbst) (red flour beetle) genome to identify, clone, express, and characterize a novel endo-β-1,4-glucanase we named TcEG1 (T. castaneum endoglucanase 1). Sequence analysis of a full-length TcEG1 cDNA clone (1356bp) revealed sequence homology to enzymes in glycosyl hydrolase family 9 (GHF9), and verified presence of a change (Gly for Ser) in the conserved catalytic domain for GHF9 cellulases. This TcEG1 cDNA clone was predicted to encode a 49.5kDa protein with a calculated pI of 5.39. Heterologous expression of TcEG1 in Drosophila S2 cell cultures resulted in secretion of a 51-kDa protein, as determined by Western blotting. The expressed protein was used to characterize TcEG1 enzymatic activity against two cellulose substrates to determine its specificity and stability. Our data support that TcEG1 as a novel endo-β-1,4-glucanase, the first functional characterization of a cellulase enzyme derived from an insect genome with potential applications in the biofuel industry due to its high relative activity at alkaline pH.

Discovering conserved insect microRNAs from expressed sequence tags

MicroRNAs (miRNA) participate in regulating diverse biological pathways by translational repression in animals. They have attracted increasing attention recently. However, little work has been done on the miRNA genes in agriculturally important pests. Because the transcripts of most miRNA genes are the products of type-II RNA polymerase, pri-miRNA has a poly(A) tail and appears in expressed sequence tags (EST). We developed a computational pipeline to identify miRNA genes from insect ESTs. First, 980,697 ESTs from 63 insects were collected and used to search the nr database. The ESTs which did not share significant similarities with any known protein-coding genes were treated as non-coding ESTs. Next, known mature miRNAs were used to align with non-coding ESTs. The ESTs which contain the sequence of mature miRNA were treated as candidate ESTs. Finally, putative precursors were extracted flanking the mature miRNA region in candidate ESTs and evaluated by the Triplet-SVM algorithm. As a result, 86 miRNAs from 30 insect species were found based on a strict criterion while 330 miRNAs from 51 species were found based on a loose criterion. Evolution analysis indicated that mir-467, mir-297 and mir-466 were the highest conserved miRNA families in insects. To confirm the reliability of putative insect miRNAs, the expression profile of nine predicted miRNAs in Locusta migratoria was investigated. Eight miRNAs were successfully detected by RT-PCR. Most miRNAs were expressed ubiquitously at all examined tissues and developmental stages whereas Lmi-mir-509 was specifically expressed in the thorax of the 2nd, 4th and 5th instars and adult locust. In all, our work reported an efficient computational strategy for predicting miRNA genes from insect ESTs and presented tens of miRNAs in diverse insect species which are expected to participate in many important physiological processes.

Exclusive sequences of different genomes

We studied the distribution of 1-7 bp words in a dataset that includes 139 complete eukaryotic genomes, 33 masked eukaryotic genomes and coding regions from 35 genomes. We tested different statistical models to determine over- and under-represented words. The method described by Karlin et al. has the strongest predictive power compared to other methods. Using this method we identified over- and under-represented words consistent within a large array of taxonomic groups. Some of those words have not yet been described as exclusive. For example, CGCG is over-represented in CG-deficient organisms. We also describe exceptions for widely known exclusive words, such as CG and TA.