The transposable element-rich genome of the cereal pest Sitophilus oryzae.

BACKGROUND: The rice weevil Sitophilus oryzae is one of the most important agricultural pests, causing extensive damage to cereal in fields and to stored grains. S. oryzae has an intracellular symbiotic relationship (endosymbiosis) with the Gram-negative bacterium Sodalis pierantonius and is a valuable model to decipher host-symbiont molecular interactions. RESULTS: We sequenced the Sitophilus oryzae genome using a combination of short and long reads to produce the best assembly for a Curculionidae species to date. We show that S. oryzae has undergone successive bursts of transposable element (TE) amplification, representing 72% of the genome. In addition, we show that many TE families are transcriptionally active, and changes in their expression are associated with insect endosymbiotic state. S. oryzae has undergone a high gene expansion rate, when compared to other beetles. Reconstruction of host-symbiont metabolic networks revealed that, despite its recent association with cereal weevils (30 kyear), S. pierantonius relies on the host for several amino acids and nucleotides to survive and to produce vitamins and essential amino acids required for insect development and cuticle biosynthesis. CONCLUSIONS: Here we present the genome of an agricultural pest beetle, which may act as a foundation for pest control. In addition, S. oryzae may be a useful model for endosymbiosis, and studying TE evolution and regulation, along with the impact of TEs on eukaryotic genomes.

Systemic control of nodule formation by plant nitrogen demand requires autoregulation-dependent and independent mechanisms.

In legumes interacting with rhizobia, the formation of symbiotic organs involved in the acquisition of atmospheric nitrogen gas (N2) is dependent on the plant nitrogen (N) demand. We used Medicago truncatula plants cultivated in split-root systems to discriminate between responses to local and systemic N signaling. We evidenced a strong control of nodule formation by systemic N signaling but obtained no clear evidence of a local control by mineral nitrogen. Systemic signaling of the plant N demand controls numerous transcripts involved in root transcriptome reprogramming associated with early rhizobia interaction and nodule formation. SUPER NUMERIC NODULES (SUNN) has an important role in this control, but we found that major systemic N signaling responses remained active in the sunn mutant. Genes involved in the activation of nitrogen fixation are regulated by systemic N signaling in the mutant, explaining why its hypernodulation phenotype is not associated with higher nitrogen fixation of the whole plant. We show that the control of transcriptome reprogramming of nodule formation by systemic N signaling requires other pathway(s) that parallel the SUNN/CLE (CLAVATA3/EMBRYO SURROUNDING REGION-LIKE PEPTIDES) pathway.

Host-specific competitiveness to form nodules in Rhizobium leguminosarum symbiovar viciae.

Fabeae legumes such as pea and faba bean form symbiotic nodules with a large diversity of soil Rhizobium leguminosarum symbiovar viciae (Rlv) bacteria. However, bacteria competitive to form root nodules (CFN) are generally not the most efficient to fix dinitrogen, resulting in a decrease in legume crop yields. Here, we investigate differential selection by host plants on the diversity of Rlv. A large collection of Rlv was collected by nodule trapping with pea and faba bean from soils at five European sites. Representative genomes were sequenced. In parallel, diversity and abundance of Rlv were estimated directly in these soils using metabarcoding. The CFN of isolates was measured with both legume hosts. Pea/faba bean CFN were associated to Rlv genomic regions. Variations of bacterial pea and/or faba bean CFN explained the differential abundance of Rlv genotypes in pea and faba bean nodules. No evidence was found for genetic association between CFN and variations in the core genome, but variations in specific regions of the nod locus, as well as in other plasmid loci, were associated with differences in CFN. These findings shed light on the genetic control of CFN in Rlv and emphasise the importance of host plants in controlling Rhizobium diversity.

Responses of mature symbiotic nodules to the whole-plant systemic nitrogen signaling.

In symbiotic root nodules of legumes, terminally differentiated rhizobia fix atmospheric N2 producing an NH4+ influx that is assimilated by the plant. The plant, in return, provides photosynthates that fuel the symbiotic nitrogen acquisition. Mechanisms responsible for the adjustment of the symbiotic capacity to the plant N demand remain poorly understood. We have investigated the role of systemic signaling of whole-plant N demand on the mature N2-fixing nodules of the model symbiotic association Medicago truncatula/Sinorhizobium using split-root systems. The whole-plant N-satiety signaling rapidly triggers reductions of both N2 fixation and allocation of sugars to the nodule. These responses are associated with the induction of nodule senescence and the activation of plant defenses against microbes, as well as variations in sugars transport and nodule metabolism. The whole-plant N-deficit responses mirror these changes: a rapid increase of sucrose allocation in response to N-deficit is associated with a stimulation of nodule functioning and development resulting in nodule expansion in the long term. Physiological, transcriptomic, and metabolomic data together provide evidence for strong integration of symbiotic nodules into whole-plant nitrogen demand by systemic signaling and suggest roles for sugar allocation and hormones in the signaling mechanisms.

Coupling clearing and hybridization chain reaction approaches to investigate gene expression in organs inside intact insect heads.

Detecting RNA molecules within their natural environment inside intact arthropods has long been challenging, particularly in small organisms covered by a tanned and pigmented cuticle. Here, we have developed a methodology that enables high-resolution analysis of the spatial distribution of transcripts of interest without having to dissect tiny organs or tissues, thereby preserving their integrity. We have combined an in situ amplification approach based on hybridization chain reaction, which enhances the signal-to-noise ratio, and a clearing approach that allows the visualization of inner organs beneath the cuticle. We have implemented this methodology for the first time in Hemiptera, mapping two salivary aphid (Acyrthosiphon pisum) transcripts, the effector c002 and the salivary sheath protein SHP. With a multiplex approach, we could simultaneously detect different mRNAs in mounted pea aphid head-thorax samples and show that they were distributed in distinct secretory cells of salivary glands. RESEARCH HIGHLIGHTS: Combining hybridisation chain reaction and clearing allows the detection of transcripts in intact aphids heads. The transcripts of the two salivary proteins c002 and SHP are compartmentalized in distinct secretory cells of the principal glands.

Tyrosine pathway regulation is host-mediated in the pea aphid symbiosis during late embryonic and early larval development.

BACKGROUND: Nutritional symbioses play a central role in insects' adaptation to specialized diets and in their evolutionary success. The obligatory symbiosis between the pea aphid, Acyrthosiphon pisum, and the bacterium, Buchnera aphidicola, is no exception as it enables this important agricultural pest insect to develop on a diet exclusively based on plant phloem sap. The symbiotic bacteria provide the host with essential amino acids lacking in its diet but necessary for the rapid embryonic growth seen in the parthenogenetic viviparous reproduction of aphids. The aphid furnishes, in exchange, non-essential amino acids and other important metabolites. Understanding the regulations acting on this integrated metabolic system during the development of this insect is essential in elucidating aphid biology. RESULTS: We used a microarray-based approach to analyse gene expression in the late embryonic and the early larval stages of the pea aphid, characterizing, for the first time, the transcriptional profiles in these developmental phases. Our analyses allowed us to identify key genes in the phenylalanine, tyrosine and dopamine pathways and we identified ACYPI004243, one of the four genes encoding for the aspartate transaminase (E.C. 2.6.1.1), as specifically regulated during development. Indeed, the tyrosine biosynthetic pathway is crucial for the symbiotic metabolism as it is shared between the two partners, all the precursors being produced by B. aphidicola. Our microarray data are supported by HPLC amino acid analyses demonstrating an accumulation of tyrosine at the same developmental stages, with an up-regulation of the tyrosine biosynthetic genes. Tyrosine is also essential for the synthesis of cuticular proteins and it is an important precursor for cuticle maturation: together with the up-regulation of tyrosine biosynthesis, we observed an up-regulation of cuticular genes expression. We were also able to identify some amino acid transporter genes which are essential for the switch over to the late embryonic stages in pea aphid development. CONCLUSIONS: Our data show that, in the development of A. pisum, a specific host gene set regulates the biosynthetic pathways of amino acids, demonstrating how the regulation of gene expression enables an insect to control the production of metabolites crucial for its own development and symbiotic metabolism.

Localized osmotic stress activates systemic responses to N limitation in Medicago truncatula-Sinorhizobium symbiotic plants.

In mature symbiotic root nodules, differentiated rhizobia fix atmospheric dinitrogen and provide ammonium to fulfill the plant nitrogen (N) demand. The plant enables this process by providing photosynthates to the nodules. The symbiosis is adjusted to the whole plant N demand thanks to systemic N signaling controlling nodule development. Symbiotic plants under N deficit stimulate nodule expansion and activate nodule senescence under N satiety. Besides, nodules are highly sensitive to drought. Here, we used split-root systems to characterize the systemic responses of symbiotic plants to a localized osmotic stress. We showed that polyéthylène glycol (PEG) application rapidly inhibited the symbiotic dinitrogen fixation activity of nodules locally exposed to the treatment, resulting to the N limitation of the plant supplied exclusively by symbiotic dinitrogen fixation. The localized PEG treatment triggered systemic signaling stimulating nodule development in the distant untreated roots. This response was associated with an enhancement of the sucrose allocation. Our analyses showed that transcriptomic reprogramming associated with PEG and N deficit systemic signaling(s) shared many targets transcripts. Altogether, our study suggests that systemic N signaling is a component of the adaptation of the symbiotic plant to the local variations of its edaphic environment.

Multimodal dynamic response of the Buchnera aphidicola pLeu plasmid to variations in leucine demand of its host, the pea aphid Acyrthosiphon pisum.

Aphids, important agricultural pests, can grow and reproduce thanks to their intimate symbiosis with the γ-proteobacterium Buchnera aphidicola that furnishes them with essential amino acids lacking in their phloem sap diet. To study how B. aphidicola, with its reduced genome containing very few transcriptional regulators, responds to variations in the metabolic requirements of its host, we concentrated on the leucine metabolic pathway. We show that leucine is a limiting factor for aphid growth and it displays a stimulatory feeding effect. Our metabolic analyses demonstrate that symbiotic aphids are able to respond to leucine starvation or excess by modulating the neosynthesis of this amino acid. At a molecular level, this response involves an early important transcriptional regulation (after 12 h of treatment) followed by a moderate change in the pLeu plasmid copy number. Both responses are no longer apparent after 7 days of treatment. These experimental data are discussed in the light of a re-annotation of the pLeu plasmid regulatory elements. Taken together, our data show that the response of B. aphidicola to the leucine demand of its host is multimodal and dynamically regulated, providing new insights concerning the genetic regulation capabilities of this bacterium in relation to its symbiotic functions.

Analyzing Multifactorial RNA-Seq Experiments with DicoExpress.

The proper use of statistical modeling in NGS data analysis requires an advanced level of expertise. There has recently been a growing consensus on using generalized linear models for differential analysis of RNA-Seq data and the advantage of mixture models to perform co-expression analysis. To offer a managed setting to use these modeling approaches, we developed DiCoExpress that provides a standardized R pipeline to perform an RNA-Seq analysis. Without any particular knowledge in statistics or R programming, beginners can perform a complete RNA-Seq analysis from quality controls to co-expression through differential analysis based on contrasts inside a generalized linear model. An enrichment analysis is proposed both on the lists of differentially expressed genes, and the co-expressed gene clusters. This video tutorial is conceived as a step-by-step protocol to help users take full advantage of DiCoExpress and its potential in empowering the biological interpretation of an RNA-Seq experiment.

Identification and characterization of Vietnamese coffee bacterial endophytes displaying in vitro antifungal and nematicidal activities.

The endophytic bacteria were isolated from coffee roots and seeds in Vietnam and identified with 16S rDNA sequencing as belonging to the Actinobacteria, Firmicutes and Proteobacteria phyla with the Nocardia, Bacillus and Burkholderia as dominant genera, respectively. Out of the thirty genera recovered from Coffea canephora and Coffea liberica, twelve were reported for the first time in endophytic association with coffee including members of the genera Brachybacterium, Caballeronia, Kitasatospora, Lechevalieria, Leifsonia, Luteibacter, Lysinibacillus, Mycolicibacterium, Nakamurella, Paracoccus, Sinomonas and Sphingobium. A total of eighty bacterial endophytes were characterized in vitro for several plant growth promoting and biocontrol traits including: the phosphate solubilization, the indolic compounds, siderophores, HCN, esterase, lipase, gelatinase and chitinase production. A subset of fifty selected bacteria were tested for their potential as biocontrol agents with in vitro confrontations with the fungal pathogen Fusarium oxysporum as well as the coffee parasitic nematodes Radopholus duriophilus and Pratylenchus coffeae. The three most efficient isolates on F. oxysporum belonging to the Bacillus, Burkholderia, and Streptomyces genera displayed a growth inhibition rate higher than 40%. Finally, five isolates from the Bacillus genus were able to lead to 100% of mortality in 24 h on both R. duriophilus and P. coffeae.

DiCoExpress: a tool to process multifactorial RNAseq experiments from quality controls to co-expression analysis through differential analysis based on contrasts inside GLM models.

BACKGROUND: RNAseq is nowadays the method of choice for transcriptome analysis. In the last decades, a high number of statistical methods, and associated bioinformatics tools, for RNAseq analysis were developed. More recently, statistical studies realised neutral comparison studies using benchmark datasets, shedding light on the most appropriate approaches for RNAseq data analysis. RESULTS: DiCoExpress is a script-based tool implemented in R that includes methods chosen based on their performance in neutral comparisons studies. DiCoExpress uses pre-existing R packages including FactoMineR, edgeR and coseq, to perform quality control, differential, and co-expression analysis of RNAseq data. Users can perform the full analysis, providing a mapped read expression data file and a file containing the information on the experimental design. Following the quality control step, the user can move on to the differential expression analysis performed using generalized linear models thanks to the automated contrast writing function. A co-expression analysis is implemented using the coseq package. Lists of differentially expressed genes and identified co-expression clusters are automatically analyzed for enrichment of annotations provided by the user. We used DiCoExpress to analyze a publicly available RNAseq dataset on the transcriptional response of Brassica napus L. to silicon treatment in plant roots and mature leaves. This dataset, including two biological factors and three replicates for each condition, allowed us to demonstrate in a tutorial all the features of DiCoExpress. CONCLUSIONS: DiCoExpress is an R script-based tool allowing users to perform a full RNAseq analysis from quality controls to co-expression analysis through differential analysis based on contrasts inside generalized linear models. DiCoExpress focuses on the statistical modelling of gene expression according to the experimental design and facilitates the data analysis leading the biological interpretation of the results.

GenoSNP: a variational Bayes within-sample SNP genotyping algorithm that does not require a reference population.

UNLABELLED: Current genotyping algorithms typically call genotypes by clustering allele-specific intensity data on a single nucleotide polymorphism (SNP) by SNP basis. This approach assumes the availability of a large number of control samples that have been sampled on the same array and platform. We have developed a SNP genotyping algorithm for the Illumina Infinium SNP genotyping assay that is entirely within-sample and does not require the need for a population of control samples nor parameters derived from such a population. Our algorithm exhibits high concordance with current methods and >99% call accuracy on HapMap samples. The ability to call genotypes using only within-sample information makes the method computationally light and practical for studies involving small sample sizes and provides a valuable independent quality control metric for other population-based approaches. AVAILABILITY: http://www.stats.ox.ac.uk/~giannoul/GenoSNP/.

hsa-miR-210 Is induced by hypoxia and is an independent prognostic factor in breast cancer.

PURPOSE: MicroRNA (miRNA) expression alterations have been described in cancer. Many cancers are characterized by areas of hypoxia, enhanced hypoxia-inducible factor (HIF) levels, and increased expression of hypoxically regulated genes, all of which correlate with patient outcome. We examined hypoxia-induced miRNA expression changes to identify markers of survival in breast cancer. EXPERIMENTAL DESIGN: We used microarrays to analyze miRNA expression changes induced by hypoxia in MCF7 breast cancer cell lines and validated results by quantitative-PCR (Q-PCR). Small interfering RNA against HIF-1alpha and HIF-2alpha, and RCC4 cells transfected with the von Hippel-Lindau (VHL) protein were used to investigate HIF dependency of miRNA expression. miRNA Q-PCR assays were done on 219 early breast cancer samples with long-term follow-up. Correlation of expression with clinical variables was done using Pearson and Spearman's rank tests, univariate, and Cox multivariate analysis. RESULTS: hsa-miR-210 induction was the most significant change under hypoxia by microarray analysis (3.4-fold, P < 0.001). hsa-miR-210 expression changes were validated by Q-PCR and detected in other cancer cell lines. Using small interfering RNAs and RCC4 cells transfected with VHL, we showed that the regulation by hypoxia of hsa-miR-210 was mediated by the HIF-1alpha/VHL transcriptional system but not HIF-2alpha. hsa-miR-210 expression levels in breast cancer samples correlated directly with a hypoxia score based on the expression of 99 genes. hsa-miR-210 expression levels showed an inverse correlation with disease-free and overall survival, significant in both univariate and multivariate analyses. CONCLUSIONS: We show that hsa-miR-210 overexpression is induced by hypoxia in a HIF-1alpha- and VHL-dependent fashion and its expression levels in breast cancer samples are an independent prognostic factor.

QuantiSNP: an Objective Bayes Hidden-Markov Model to detect and accurately map copy number variation using SNP genotyping data.

Array-based technologies have been used to detect chromosomal copy number changes (aneuploidies) in the human genome. Recent studies identified numerous copy number variants (CNV) and some are common polymorphisms that may contribute to disease susceptibility. We developed, and experimentally validated, a novel computational framework (QuantiSNP) for detecting regions of copy number variation from BeadArray SNP genotyping data using an Objective Bayes Hidden-Markov Model (OB-HMM). Objective Bayes measures are used to set certain hyperparameters in the priors using a novel re-sampling framework to calibrate the model to a fixed Type I (false positive) error rate. Other parameters are set via maximum marginal likelihood to prior training data of known structure. QuantiSNP provides probabilistic quantification of state classifications and significantly improves the accuracy of segmental aneuploidy identification and mapping, relative to existing analytical tools (Beadstudio, Illumina), as demonstrated by validation of breakpoint boundaries. QuantiSNP identified both novel and validated CNVs. QuantiSNP was developed using BeadArray SNP data but it can be adapted to other platforms and we believe that the OB-HMM framework has widespread applicability in genomic research. In conclusion, QuantiSNP is a novel algorithm for high-resolution CNV/aneuploidy detection with application to clinical genetics, cancer and disease association studies.

Matrix-assisted laser desorption/ionisation, time-of-flight mass spectrometry in genomics research.

The beginning of this millennium has seen dramatic advances in genomic research. Milestones such as the complete sequencing of the human genome and of many other species were achieved and complemented by the systematic discovery of variation at the single nucleotide (SNP) and whole segment (copy number polymorphism) level. Currently most genomics research efforts are concentrated on the production of whole genome functional annotations, as well as on mapping the epigenome by identifying the methylation status of CpGs, mainly in CpG islands, in different tissues. These recent advances have a major impact on the way genetic research is conducted and have accelerated the discovery of genetic factors contributing to disease. Technology was the critical driving force behind genomics projects: both the combination of Sanger sequencing with high-throughput capillary electrophoresis and the rapid advances in microarray technologies were keys to success. MALDI-TOF MS-based genome analysis represents a relative newcomer in this field. Can it establish itself as a long-term contributor to genetics research, or is it only suitable for niche areas and for laboratories with a passion for mass spectrometry? In this review, we will highlight the potential of MALDI-TOF MS-based tools for resequencing and for epigenetics research applications, as well as for classical complex genetic studies, allele quantification, and quantitative gene expression analysis. We will also identify the current limitations of this approach and attempt to place it in the context of other genome analysis technologies.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome.

BACKGROUND: The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae. RESULTS: The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding. CONCLUSIONS: With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Bacteriocyte Reprogramming to Cope With Nutritional Stress in a Phloem Sap Feeding Hemipteran, the Pea Aphid Acyrthosiphon pisum.

Nutritional symbioses play a central role in the ability of insects to thrive on unbalanced diets and in ensuring their evolutionary success. A genomic model for nutritional symbiosis comprises the hemipteran Acyrthosiphon pisum, and the gamma-3-proteobacterium, Buchnera aphidicola, with genomes encoding highly integrated metabolic pathways. A. pisum feeds exclusively on plant phloem sap, a nutritionally unbalanced diet highly variable in composition, thus raising the question of how this symbiotic system responds to nutritional stress. We addressed this by combining transcriptomic, phenotypic and life history trait analyses to determine the organismal impact of deprivation of tyrosine and phenylalanine. These two aromatic amino acids are essential for aphid development, are synthesized in a metabolic pathway for which the aphid host and the endosymbiont are interdependent, and their concentration can be highly variable in plant phloem sap. We found that this nutritional challenge does not have major phenotypic effects on the pea aphid, except for a limited weight reduction and a 2-day delay in onset of nymph laying. Transcriptomic analyses through aphid development showed a prominent response in bacteriocytes (the core symbiotic tissue which houses the symbionts), but not in gut, thus highlighting the role of bacteriocytes as major modulators of this homeostasis. This response does not involve a direct regulation of tyrosine and phenylalanine biosynthetic pathway and transporter genes. Instead, we observed an extensive transcriptional reprogramming of the bacteriocyte with a rapid down-regulation of genes encoding sugar transporters and genes required for sugar metabolism. Consistently, we observed continued overexpression of the A. pisum homolog of RRAD, a small GTPase implicated in repressing aerobic glycolysis. In addition, we found increased transcription of genes involved in proliferation, cell size control and signaling. We experimentally confirmed the significance of these gene expression changes detecting an increase in bacteriocyte number and cell size in vivo under tyrosine and phenylalanine depletion. Our results support a central role of bacteriocytes in the aphid response to amino acid deprivation: their transcriptional and cellular responses fine-tune host physiology providing the host insect with an effective way to cope with the challenges posed by the variability in composition of phloem sap.

Rapid transcriptional plasticity of duplicated gene clusters enables a clonally reproducing aphid to colonise diverse plant species.

BACKGROUND: The prevailing paradigm of host-parasite evolution is that arms races lead to increasing specialisation via genetic adaptation. Insect herbivores are no exception and the majority have evolved to colonise a small number of closely related host species. Remarkably, the green peach aphid, Myzus persicae, colonises plant species across 40 families and single M. persicae clonal lineages can colonise distantly related plants. This remarkable ability makes M. persicae a highly destructive pest of many important crop species. RESULTS: To investigate the exceptional phenotypic plasticity of M. persicae, we sequenced the M. persicae genome and assessed how one clonal lineage responds to host plant species of different families. We show that genetically identical individuals are able to colonise distantly related host species through the differential regulation of genes belonging to aphid-expanded gene families. Multigene clusters collectively upregulate in single aphids within two days upon host switch. Furthermore, we demonstrate the functional significance of this rapid transcriptional change using RNA interference (RNAi)-mediated knock-down of genes belonging to the cathepsin B gene family. Knock-down of cathepsin B genes reduced aphid fitness, but only on the host that induced upregulation of these genes. CONCLUSIONS: Previous research has focused on the role of genetic adaptation of parasites to their hosts. Here we show that the generalist aphid pest M. persicae is able to colonise diverse host plant species in the absence of genetic specialisation. This is achieved through rapid transcriptional plasticity of genes that have duplicated during aphid evolution.

ArthropodaCyc: a CycADS powered collection of BioCyc databases to analyse and compare metabolism of arthropods.

Arthropods interact with humans at different levels with highly beneficial roles (e.g. as pollinators), as well as with a negative impact for example as vectors of human or animal diseases, or as agricultural pests. Several arthropod genomes are available at present and many others will be sequenced in the near future in the context of the i5K initiative, offering opportunities for reconstructing, modelling and comparing their metabolic networks. In-depth analysis of these genomic data through metabolism reconstruction is expected to contribute to a better understanding of the biology of arthropods, thereby allowing the development of new strategies to control harmful species. In this context, we present here ArthropodaCyc, a dedicated BioCyc collection of databases using the Cyc annotation database system (CycADS), allowing researchers to perform reliable metabolism comparisons of fully sequenced arthropods genomes. Since the annotation quality is a key factor when performing such global genome comparisons, all proteins from the genomes included in the ArthropodaCyc database were re-annotated using several annotation tools and orthology information. All functional/domain annotation results and their sources were integrated in the databases for user access. Currently, ArthropodaCyc offers a centralized repository of metabolic pathways, protein sequence domains, Gene Ontology annotations as well as evolutionary information for 28 arthropod species. Such database collection allows metabolism analysis both with integrated tools and through extraction of data in formats suitable for systems biology studies.Database URL: http://arthropodacyc.cycadsys.org/.