Transcriptome sequencing, and rapid development and application of SNP markers for the legume pod borer Maruca vitrata (Lepidoptera: Crambidae)

De Novo Mining and Validating Novel Microsatellite Markers to Assess Genetic Diversity in Maruca vitrata (F.), a Legume Pod Borer

Maruca vitrata (Fabricius) is an invasive insect pest capable of causing enormous economic losses to a broad spectrum of leguminous crops. Microsatellites are valuable molecular markers for population genetic studies; however, an inadequate number of M. vitrata microsatellite loci are available to carry out population association studies. Thus, we utilized this insect's public domain databases for mining expressed sequence tags (EST)-derived microsatellite markers. In total, 234 microsatellite markers were identified from 10053 unigenes. We discovered that trinucleotide repeats were the most predominant microsatellite motifs (61.53%), followed by dinucleotide repeats (23.50%) and tetranucleotide repeats (14.95%). Based on the analysis, twenty-five markers were selected for validation in M. vitrata populations collected from various regions of India. The number of alleles (Na), observed heterozygosity (Ho), and expected heterozygosity (He) ranged from 2 to 5; 0.00 to 0.80; and 0.10 to 0.69, respectively. The polymorphic loci showed polymorphism information content (PIC), ranging from 0.09 to 0.72. Based on the genetic distance matrix, the unrooted neighbor-joining dendrogram differentiated the selected populations into two discrete groups. The SSR markers developed and validated in this study will be helpful in population-level investigations of M. vitrata to understand the gene flow, demography, dispersal patterns, biotype differentiation, and host dynamics.

Analysis of Genetic Diversity of Banana Weevils (Cosmopolites sordidus) (Coleoptera: Curculionidae) Using Transcriptome-Derived Simple Sequence Repeat Markers

The banana weevil, Cosmopolites sordidus (Germar) (Coleoptera: Curculionidae) is an economically important insect pest of bananas. It causes up to 100% yield losses and substantial lifespan reduction in bananas. Advances in genomics, proteomics, and sequencing technologies have provided powerful pathways to genotyping disastrous pests such as C. sordidus. However, such technologies are often not available to the majority of rural subtropical African banana growers and pest control managers. This study was therefore motivated by the need to create cheap and easily accessible C. sordidus genotyping methods that could be deployed by banana pest control managers to the benefit of C. sordidus control programs in the tropics where such advanced technologies are not readily accessible. We used an in-house C. sordidus transcriptome from the an-ongoing study from which we mined an array of simple sequence repeat (SSR) markers. Of these, six highly polymorphic transcriptome-derived SSR markers were used to successfully genotype within and among banana weevil population genetic diversity of 12 C. sordidus populations collected from four banana-growing agro-ecological zones (AEZs) in Uganda. The developed transcriptome-derived SSR markers can be used by researchers in population genetics for characterization of the C. sordidus and identification of new genes that are linked to traits of particular interest. The significant genetic diversity revealed in C. sordidus provides pertinent information for integrated pest management strategies.

Natural Pest Regulation and Its Compatibility with Other Crop Protection Practices in Smallholder Bean Farming Systems

Simple Summary

Bean production by smallholder farmers in sub-Saharan Africa is frequently constrained by insect pests, two of the most serious being Maruca vitrata and Aphis fabae. For many bean farmers, the options available to control these pests are limited. A few can access synthetic insecticides, but these have negative consequences for their health and the environment. Natural pest regulation (NPR) offers environmentally benign approaches for smallholders to manage bean pests. For example, here, we focus on biological control whereby beneficial organisms predate or parasitize the pests. Field studies show this is a feasible strategy for controlling M. vitrata and A. fabae. In particular, we highlight how compatible biological control is with other NPR options, such as the use of biopesticides (including plant extracts), resistant varieties, and cultural control. We recommend that smallholder farmers consider biological control alongside other NPR strategies for reducing the populations of A. fabae and M. vitrata in the common bean, increasing the yields and reducing the negative impacts of the synthetic pesticides.

Abstract

Common bean (Phaseolus vulgaris) production and storage are limited by numerous constraints. Insect pests are often the most destructive. However, resource-constrained smallholders in sub-Saharan Africa (SSA) often do little to manage pests. Where farmers do use a control strategy, it typically relies on chemical pesticides, which have adverse effects on the wildlife, crop pollinators, natural enemies, mammals, and the development of resistance by pests. Nature-based solutions —in particular, using biological control agents with sustainable approaches that include biopesticides, resistant varieties, and cultural tools—are alternatives to chemical control. However, significant barriers to their adoption in SSA include a lack of field data and knowledge on the natural enemies of pests, safety, efficacy, the spectrum of activities, the availability and costs of biopesticides, the lack of sources of resistance for different cultivars, and spatial and temporal inconsistencies for cultural methods. Here, we critically review the control options for bean pests, particularly the black bean aphid (Aphis fabae) and pod borers (Maruca vitrata). We identified natural pest regulation as the option with the greatest potential for this farming system. We recommend that farmers adapt to using biological control due to its compatibility with other sustainable approaches, such as cultural tools, resistant varieties, and biopesticides for effective management, especially in SSA.

Transcriptomic Analysis of Aedes aegypti in Response to Mosquitocidal Bacillus thuringiensis LLP29 Toxin

Globally, Aedes aegypti is one of the most dangerous mosquitoes that plays a crucial role as a vector for human diseases, such as yellow fever, dengue, and chikungunya. To identify (1) transcriptomic basis of midgut (2) key genes that are involved in the toxicity process by a comparative transcriptomic analysis between the control and Bacillus thuringiensis (Bt) toxin (LLP29 proteins)-treated groups. Next-generation sequencing technology was used to sequence the midgut transcriptome of A. aegypti. A total of 17130 unigenes, including 574 new unigenes, were identified containing 16358 (95.49%) unigenes that were functionally annotated. According to differentially expressed gene (DEG) analysis, 557 DEGs were annotated, including 226 upregulated and 231 downregulated unigenes in the Bt toxin-treated group. A total of 442 DEGs were functionally annotated; among these, 33 were specific to multidrug resistance, 6 were immune-system-related (Lectin, Defensin, Lysozyme), 28 were related to putative proteases, 7 were lipase-related, 8 were related to phosphatases, and 30 were related to other transporters. In addition, the relative expression of 28 DEGs was further confirmed through quantitative real time polymerase chain reaction. The results provide a transcriptomic basis for the identification and functional authentication of DEGs in A. aegypti.

Transcriptome analysis in the beet webworm, Spoladea recurvalis (Lepidoptera: Crambidae)

The beet webworm, Spoladea recurvalis Fabricius, is a destructive pest on vegetable crops in tropics and subtropics; its main host plant is amaranth. It has become imperative to develop non-chemical methods to control S. recurvalis on amaranth. However, the lack of molecular information about this species has hindered the development of novel pest management strategies. In this study, high-throughput RNA sequencing covering de novo sequence assemblies, functional annotation of transcripts, gene function classification and enrichment was performed on S. recurvalis. Illumina sequencing generated a total of 120 435 transcript contigs ranging from 201 to 22 729 bases with a mean length of 688 bases. The assembled transcripts were subjected to Basic Local Alignment Search Tool-X (BLASTX) to obtain the annotations against non-redundant, Swiss-Prot, Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) protein databases. A subset of 58 225 transcript sequences returned hits from known proteins in the National Center for Biotechnology Information database, and the majority of the transcript sequences had the highest number of hits for Danaus plexippus (50.43%). A total of 1217 Gene Ontology-level 3 annotations were assigned to 51 805 transcripts, while 39 650 transcripts were predicted as functional protein-coding genes in the COG database and 20 037 transcripts were enriched to KEGG pathways. We identified 40 putative genes related to pheromone production and reception in S. recurvalis, with the expression of one gene between 0.29 and 1141.79 fragments per kilo base per million (FPKM) reads. The transcriptome sequence of S. recurvalis is a first step toward offering a comprehensive genomic resource which would enable better understanding of molecular mechanisms to enable development of effective pest management practices for this species.

Quantitative proteomic analysis of the fall armyworm saliva

Lepidopteran larvae secrete saliva on plant tissues during feeding. Components in the saliva may aid in food digestion, whereas other components are recognized by plants as cues to elicit defense responses. Despite the ecological and economical importance of these plant-feeding insects, knowledge of their saliva composition is limited to a few species. In this study, we identified the salivary proteins of larvae of the fall armyworm (FAW), Spodoptera frugiperda; determined qualitative and quantitative differences in the salivary proteome of the two host races-corn and rice strains-of this insect; and identified changes in total protein concentration and relative protein abundance in the saliva of FAW larvae associated with different host plants. Quantitative proteomic analyses were performed using labeling with isobaric tags for relative and absolute quantification followed by liquid chromatography-tandem mass spectrometry. In total, 98 proteins were identified (>99% confidence) in the FAW saliva. These proteins were further categorized into five functional groups: proteins potentially involved in (1) plant defense regulation, (2) herbivore offense, (3) insect immunity, (4) detoxification, (5) digestion, and (6) other functions. Moreover, there were differences in the salivary proteome between the FAW strains that were identified by label-free proteomic analyses. Thirteen differentially identified proteins were present in each strain. There were also differences in the relative abundance of eleven salivary proteins between the two FAW host strains as well as differences within each strain associated with different diets. The total salivary protein concentration was also different for the two strains reared on different host plants. Based on these results, we conclude that the FAW saliva contains a complex mixture of proteins involved in different functions that are specific for each strain and its composition can change plastically in response to diet type.

De novo transcriptome analysis and microsatellite marker development for population genetic study of a serious insect pest, Rhopalosiphum padi (L.) (Hemiptera: Aphididae)

The bird cherry-oat aphid, Rhopalosiphum padi (L.), is one of the most abundant aphid pests of cereals and has a global distribution. Next-generation sequencing (NGS) is a rapid and efficient method for developing molecular markers. However, transcriptomic and genomic resources of R. padi have not been investigated. In this study, we used transcriptome information obtained by RNA-Seq to develop polymorphic microsatellites for investigating population genetics in this species. The transcriptome of R. padi was sequenced on an Illumina HiSeq 2000 platform. A total of 114.4 million raw reads with a GC content of 40.03% was generated. The raw reads were cleaned and assembled into 29,467 unigenes with an N50 length of 1,580 bp. Using several public databases, 82.47% of these unigenes were annotated. Of the annotated unigenes, 8,022 were assigned to COG pathways, 9,895 were assigned to GO pathways, and 14,586 were mapped to 257 KEGG pathways. A total of 7,936 potential microsatellites were identified in 5,564 unigenes, 60 of which were selected randomly and amplified using specific primer pairs. Fourteen loci were found to be polymorphic in the four R. padi populations. The transcriptomic data presented herein will facilitate gene discovery, gene analyses, and development of molecular markers for future studies of R. padi and other closely related aphid species.

Genomics of Lepidoptera saliva reveals function in herbivory

Lepidoptera herbivores deposit copious amounts of saliva when feeding. Their saliva is produced by the paired mandibular and labial glands and evidence indicates that it may play an important role in allowing an herbivore to establish on its host plant. Genomic studies of Lepidoptera saliva are beginning to reveal the role of saliva in herbivory. Molecules involved in digestion, detoxification, immunity, defense against plant secondary chemicals, chemoreception and so on have been identified using high throughput genomic tools. These genomic tools have also revealed changes that occur in Lepidoptera saliva when caterpillars feed on different host plants. However, there are other factors either biotic or abiotic (e.g., larval stage, larval health, temperature, water stress, etc.) that might also affect its composition. Though further functional and ecological studies are still necessary to fully understand the role of Lepidoptera saliva on herbivory, here we review current trends.

Transcriptome Analysis and Identification of Major Detoxification Gene Families and Insecticide Targets in Grapholita Molesta (Busck) (Lepidoptera: Tortricidae)

The oriental fruit moth, Grapholita molesta (Busck) (Lepidoptera: Tortricidae), is an important pest of most stone and pome fruits and causes serious damage to the fruit industry worldwide. This insect pest has been primarily controlled through the application of insecticides; as a result, G. molesta has developed resistance to many different types of insecticides. To identify detoxification genes, we have, de novo, sequenced the transcriptome of G. molesta (Lepidoptera: Tortricidae) and yielded 58,970 unigenes of which 26,985 unigenes matched to known proteins. In total, 2,040 simple sequence repeats have been identified. The comprehensive transcriptome data set has permitted us to identify members of important gene families related to detoxification in G. molesta, including 77 unigenes of putative cytochrome P450s, 28 of glutathione S-transferases, 46 of Carboxylesterases, and 31 of insecticide targets. Orthologs of some of these unigenes have shown to play a pivotal role in insecticide resistance in other insect species and those unigenes likely have similar functions in G. molesta.

Sex pheromone recognition and characterization of three pheromone-binding proteins in the legume pod borer, Maruca vitrata Fabricius (Lepidoptera: Crambidae)

Pheromone-binding proteins (PBPs) are essential for the filtering, binding and transporting of sex pheromones across sensillum lymph to membrane-associated pheromone receptors of moths. In this study, three novel PBP genes were expressed in Escherichia coli to examine their involvement in the sex pheromone perception of Maruca vitrata. Fluorescence binding experiments indicated that MvitPBP1-3 had strong binding affinities with four sex pheromones. Moreover, molecular docking results demonstrated that six amino acid residues of three MvitPBPs were involved in the binding of the sex pheromones. These results suggested that MvitPBP1-3 might play critical roles in the perception of female sex pheromones. Additionally, the binding capacity of MvitPBP3 with the host-plant floral volatiles was high and was similar to that of MvitGOBP2. Furthermore, sequence alignment and docking analysis showed that both MvitGOBP2 and MvitPBP3 possessed an identical key binding site (arginine, R130/R140) and a similar protein pocket structure around the binding cavity. Therefore, we hypothesized that MvitPBP3 and MvitGOBP2 might have synergistic roles in binding different volatile ligands. In combination, the use of synthetic sex pheromones and floral volatiles from host-plant may be used in the exploration for more efficient monitoring and integrated management strategies for the legume pod borer in the field.

Transcriptome Analysis and Discovery of Genes Relevant to Development in Bradysia odoriphaga at Three Developmental Stages

Bradysia odoriphaga (Diptera: Sciaridae) is the most important pest of Chinese chive (Allium tuberosum) in Asia; however, the molecular genetics are poorly understood. To explore the molecular biological mechanism of development, Illumina sequencing and de novo assembly were performed in the third-instar, fourth-instar, and pupal B. odoriphaga. The study resulted in 16.2 Gb of clean data and 47,578 unigenes (≥125 bp) contained in 7,632,430 contigs, 46.21% of which were annotated from non-redundant protein (NR), Gene Ontology (GO), Clusters of Orthologous Groups (COG), Eukaryotic Orthologous Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. It was found that 19.67% of unigenes matched the homologous species mainly, including Aedes aegypti, Culex quinquefasciatus, Ceratitis capitata, and Anopheles gambiae. According to differentially expressed gene (DEG) analysis, 143, 490, and 309 DEGs were annotated as involved in the developmental process in the GO database respectively, in the comparisons of third-instar and fourth-instar larvae, third-instar larvae and pupae, and fourth-instar larvae and pupae. Twenty-five genes were closely related to these processes, including developmental process, reproduction process, and reproductive organs development and programmed cell death (PCD). The information of unigenes assembled in B. odoriphaga through transcriptome and DEG analyses could provide a detailed genetic basis and regulated information for elaborating the developmental mechanism from the larval, pre-pupal to pupal stages of B. odoriphaga.

Generation of a Transcriptome in a Model Lepidopteran Pest, Heliothis virescens, Using Multiple Sequencing Strategies for Profiling Midgut Gene Expression

Heliothine pests such as the tobacco budworm, Heliothis virescens (F.), pose a significant threat to production of a variety of crops and ornamental plants and are models for developmental and physiological studies. The efforts to develop new control measures for H. virescens, as well as its use as a relevant biological model, are hampered by a lack of molecular resources. The present work demonstrates the utility of next-generation sequencing technologies for rapid molecular resource generation from this species for which lacks a sequenced genome. In order to amass a de novo transcriptome for this moth, transcript sequences generated from Illumina, Roche 454, and Sanger sequencing platforms were merged into a single de novo transcriptome assembly. This pooling strategy allowed a thorough sampling of transcripts produced under diverse environmental conditions, developmental stages, tissues, and infections with entomopathogens used for biological control, to provide the most complete transcriptome to date for this species. Over 138 million reads from the three platforms were assembled into the final set of 63,648 contigs. Of these, 29,978 had significant BLAST scores indicating orthologous relationships to transcripts of other insect species, with the top-hit species being the monarch butterfly (Danaus plexippus) and silkworm (Bombyx mori). Among identified H. virescens orthologs were immune effectors, signal transduction pathways, olfactory receptors, hormone biosynthetic pathways, peptide hormones and their receptors, digestive enzymes, and insecticide resistance enzymes. As an example, we demonstrate the utility of this transcriptomic resource to study gene expression profiling of larval midguts and detect transcripts of putative Bacillus thuringiensis (Bt) Cry toxin receptors. The substantial molecular resources described in this study will facilitate development of H. virescens as a relevant biological model for functional genomics and for new biological experimentation needed to develop efficient control efforts for this and related Noctuid pest moths.

Novel insights into the insect trancriptome response to a natural DNA virus

Background

Little is known about invertebrate responses to DNA viruses. Here, we infect a commercially important pest moth species Plodia interpunctella with its naturally infecting DNA virus. We sequenced, assembled and annotated the complete transcriptome of the moth, and a partial transcriptome of the virus. We then tested for differential gene expression between moths that were exposed to the virus and controls.

Results

We found 51 genes that were differentially expressed in moths exposed to a DNA baculovirus compared to controls. Gene set enrichment analysis revealed that cuticle proteins were significantly overrepresented in this group of genes. Interestingly, 6 of the 7 differentially expressed cuticle proteins were downregulated, suggesting that baculoviruses are able to manipulate its host’s response. In fact, an additional 29 of the 51 genes were also downregulated in exposed compared with control animals, including a gram-negative binding protein. In contrast, genes involved in transposable element movement were upregulated after infection.

Conclusions

We present the first experiment to measure genome-wide gene expression in an insect after infection with a natural DNA virus. Our results indicate that cuticle proteins might be key genes underpinning the response to DNA viruses. Furthermore, the large proportion of genes that were downregulated after viral exposure suggests that this virus is actively manipulating the insect immune response. Finally, it appears that transposable element activity might increase during viral invasion. Combined, these results provide much needed host candidate genes that respond to DNA viral invaders.

Transcriptome profiling of the testis reveals genes involved in spermatogenesis and marker discovery in the oriental fruit fly, Bactrocera dorsalis

The testis is a highly specialized tissue that plays a vital role in ensuring fertility by producing spermatozoa, which are transferred to the female during mating. Spermatogenesis is a complex process, resulting in the production of mature sperm, and involves significant structural and biochemical changes in the seminiferous epithelium of the adult testis. The identification of genes involved in spermatogenesis of Bactrocera dorsalis (Hendel) is critical for a better understanding of its reproductive development. In this study, we constructed a cDNA library of testes from male B. dorsalis adults at different ages, and performed de novo transcriptome sequencing to produce a comprehensive transcript data set, using Illumina sequencing technology. The analysis yielded 52 016 732 clean reads, including a total of 4.65 Gb of nucleotides. These reads were assembled into 47 677 contigs (average 443 bp) and then clustered into 30 516 unigenes (average 756 bp). Based on BLAST hits with known proteins in different databases, 20 921 unigenes were annotated with a cut-off E-value of 10(-5). The transcriptome sequences were further annotated using the Clusters of Orthologous Groups, Gene Orthology and the Kyoto Encyclopedia of Genes and Genomes databases. Functional genes involved in spermatogenesis were analysed, including cell cycle proteins, metalloproteins, actin, and ubiquitin and antihyperthermia proteins. Several testis-specific genes were also identified. The transcripts database will help us to understand the molecular mechanisms underlying spermatogenesis in B. dorsalis. Furthermore, 2913 simple sequence repeats and 151 431 single nucleotide polymorphisms were identified, which will be useful for investigating the genetic diversity of B. dorsalis in the future.

High-throughput transcriptome sequencing and preliminary functional analysis in four Neotropical tree species

Background

The Amazonian rainforest is predicted to suffer from ongoing environmental changes. Despite the need to evaluate the impact of such changes on tree genetic diversity, we almost entirely lack genomic resources.

Results

In this study, we analysed the transcriptome of four tropical tree species (Carapa guianensis, Eperua falcata, Symphonia globulifera and Virola michelii) with contrasting ecological features, belonging to four widespread botanical families (respectively Meliaceae, Fabaceae, Clusiaceae and Myristicaceae). We sequenced cDNA libraries from three organs (leaves, stems, and roots) using 454 pyrosequencing. We have developed an R and bioperl-based bioinformatic procedure for de novo assembly, gene functional annotation and marker discovery. Mismatch identification takes into account single-base quality values as well as the likelihood of false variants as a function of contig depth and number of sequenced chromosomes. Between 17103 (for Symphonia globulifera) and 23390 (for Eperua falcata) contigs were assembled. Organs varied in the numbers of unigenes they apparently express, with higher number in roots. Patterns of gene expression were similar across species, with metabolism of aromatic compounds standing out as an overrepresented gene function. Transcripts corresponding to several gene functions were found to be over- or underrepresented in each organ. We identified between 4434 (for Symphonia globulifera) and 9076 (for Virola surinamensis) well-supported mismatches. The resulting overall mismatch density was comprised between 0.89 (S. globulifera) and 1.05 (V. surinamensis) mismatches/100 bp in variation-containing contigs.

Conclusion

The relative representation of gene functions in the four transcriptomes suggests that secondary metabolism may be particularly important in tropical trees. The differential representation of transcripts among tissues suggests differential gene expression, which opens the way to functional studies in these non-model, ecologically important species. We found substantial amounts of mismatches in the four species. These newly identified putative variants are a first step towards acquiring much needed genomic resources for tropical tree species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-238) contains supplementary material, which is available to authorized users.

Comparative analysis of two phenologically divergent populations of the pine processionary moth (Thaumetopoea pityocampa) by de novo transcriptome sequencing

The pine processionary moth Thaumetopoea pityocampa is a Mediterranean lepidopteran defoliator that experiences a rapid range expansion towards higher latitudes and altitudes due to the current climate warming. Its phenology - the time of sexual reproduction - is certainly a key trait for the local adaptation of the processionary moth to climatic conditions. Moreover, an exceptional case of allochronic differentiation was discovered ca. 15 years ago in this species. A population with a shifted phenology (the summer population, SP) co-exists near Leiria, Portugal, with a population following the classical cycle (the winter population, WP). The existence of this population is an outstanding opportunity to decipher the genetic bases of phenology. No genomic resources were so far available for T. pityocampa. We developed a high-throughput sequencing approach to build a first reference transcriptome, and to proceed with comparative analyses of the sympatric SP and WP. We pooled RNA extracted from whole individuals of various developmental stages, and performed a transcriptome characterisation for both populations combining Roche 454-FLX and traditional Sanger data. The obtained sequences were clustered into ca. 12,000 transcripts corresponding to 9265 unigenes. The mean transcript coverage was 21.9 reads per bp. Almost 70% of the de novo assembled transcripts displayed significant similarity to previously published proteins and around 50% of the transcripts contained a full-length coding region. Comparative analyses of the population transcriptomes allowed to investigate genes specifically expressed in one of the studied populations only, and to identify the most divergent homologous SP/WP transcripts. The most divergent pairs of transcripts did not correspond to obvious phenology-related candidate genes, and 43% could not be functionally annotated. This study provides the first comprehensive genome-wide resource for the target species T. pityocampa. Many of the assembled genes are orthologs of published Lepidoptera genes, which allows carrying out gene-specific re-sequencing. Data mining has allowed the identification of SNP loci that will be useful for population genomic approaches and genome-wide scans of population differentiation to identify signatures of selection.

Assessment of the geographic origins of pinewood nematode isolates via single nucleotide polymorphism in effector genes

The pinewood nematode, Bursaphelenchus xylophilus, is native to North America but it only causes damaging pine wilt disease in those regions of the world where it has been introduced. The accurate detection of the species and its dispersal routes are thus essential to define effective control measures. The main goals of this study were to analyse the genetic diversity among B. xylophilus isolates from different geographic locations and identify single nucleotide polymorphism (SNPs) markers for geographic origin, through a comparative transcriptomic approach. The transcriptomes of seven B. xylophilus isolates, from Continental Portugal (4), China (1), Japan (1) and USA (1), were sequenced in the next generation platform Roche 454. Analysis of effector gene transcripts revealed inter-isolate nucleotide diversity that was validated by Sanger sequencing in the genomic DNA of the seven isolates and eight additional isolates from different geographic locations: Madeira Island (2), China (1), USA (1), Japan (2) and South Korea (2). The analysis identified 136 polymorphic positions in 10 effector transcripts. Pairwise comparison of the 136 SNPs through Neighbor-Joining and the Maximum Likelihood methods and 5-mer frequency analysis with the alignment-independent bilinear multivariate modelling approach correlated the SNPs with the isolates geographic origin. Furthermore, the SNP analysis indicated a closer proximity of the Portuguese isolates to the Korean and Chinese isolates than to the Japanese or American isolates. Each geographic cluster carried exclusive alleles that can be used as SNP markers for B. xylophilus isolate identification.

Comparative analyses of two Geraniaceae transcriptomes using next-generation sequencing

BACKGROUND

Organelle genomes of Geraniaceae exhibit several unusual evolutionary phenomena compared to other angiosperm families including accelerated nucleotide substitution rates, widespread gene loss, reduced RNA editing, and extensive genomic rearrangements. Since most organelle-encoded proteins function in multi-subunit complexes that also contain nuclear-encoded proteins, it is likely that the atypical organellar phenomena affect the evolution of nuclear genes encoding organellar proteins. To begin to unravel the complex co-evolutionary interplay between organellar and nuclear genomes in this family, we sequenced nuclear transcriptomes of two species, Geranium maderense and Pelargonium x hortorum.

RESULTS

Normalized cDNA libraries of G. maderense and P. x hortorum were used for transcriptome sequencing. Five assemblers (MIRA, Newbler, SOAPdenovo, SOAPdenovo-trans [SOAPtrans], Trinity) and two next-generation technologies (454 and Illumina) were compared to determine the optimal transcriptome sequencing approach. Trinity provided the highest quality assembly of Illumina data with the deepest transcriptome coverage. An analysis to determine the amount of sequencing needed for de novo assembly revealed diminishing returns of coverage and quality with data sets larger than sixty million Illumina paired end reads for both species. The G. maderense and P. x hortorum transcriptomes contained fewer transcripts encoding the PLS subclass of PPR proteins relative to other angiosperms, consistent with reduced mitochondrial RNA editing activity in Geraniaceae. In addition, transcripts for all six plastid targeted sigma factors were identified in both transcriptomes, suggesting that one of the highly divergent rpoA-like ORFs in the P. x hortorum plastid genome is functional.

CONCLUSIONS

The findings support the use of the Illumina platform and assemblers optimized for transcriptome assembly, such as Trinity or SOAPtrans, to generate high-quality de novo transcriptomes with broad coverage. In addition, results indicated no major improvements in breadth of coverage with data sets larger than six billion nucleotides or when sampling RNA from four tissue types rather than from a single tissue. Finally, this work demonstrates the power of cross-compartmental genomic analyses to deepen our understanding of the correlated evolution of the nuclear, plastid, and mitochondrial genomes in plants.

Development of reference transcriptomes for the major field insect pests of cowpea: a toolbox for insect pest management approaches in west Africa

Cowpea is a widely cultivated and major nutritional source of protein for many people that live in West Africa. Annual yields and longevity of grain storage is greatly reduced by feeding damage caused by a complex of insect pests that include the pod sucking bugs, Anoplocnemis curvipes Fabricius (Hemiptera: Coreidae) and Clavigralla tomentosicollis Stål (Hemiptera: Coreidae); as well as phloem-feeding cowpea aphids, Aphis craccivora Koch (Hemiptera: Aphididae) and flower thrips, Megalurothrips sjostedti Trybom (Thysanoptera: Thripidae). Efforts to control these pests remain a challenge and there is a need to understand the structure and movement of these pest populations in order to facilitate the development of integrated pest management strategies (IPM). Molecular tools have the potential to help facilitate a better understanding of pest populations. Towards this goal, we used 454 pyrosequencing technology to generate 319,126, 176,262, 320,722 and 227,882 raw reads from A. curvipes, A. craccivora, C. tomentosicollis and M. sjostedti, respectively. The reads were de novo assembled into 11,687, 7,647, 10,652 and 7,348 transcripts for A. curvipes, A. craccivora, C. tomentosicollis and M. sjostedti, respectively. Functional annotation of the resulting transcripts identified genes putatively involved in insecticide resistance, pathogen defense and immunity. Additionally, sequences that matched the primary aphid endosymbiont, Buchnera aphidicola, were identified among A. craccivora transcripts. Furthermore, 742, 97, 607 and 180 single nucleotide polymorphisms (SNPs) were respectively predicted among A. curvipes, A. craccivora, C. tomentosicollis and M. sjostedti transcripts, and will likely be valuable tools for future molecular genetic marker development. These results demonstrate that Roche 454-based transcriptome sequencing could be useful for the development of genomic resources for cowpea pest insects in West Africa.

Application of Cydia pomonella expressed sequence tags: Identification and expression of three general odorant binding proteins in codling moth

The codling moth, Cydia pomonella, is one of the most important pests of pome fruits in the world, yet the molecular genetics and the physiology of this insect remain poorly understood. A combined assembly of 8 341 expressed sequence tags was generated from Roche 454 GS-FLX sequencing of eight tissue-specific cDNA libraries. Putative chemosensory proteins (12) and odorant binding proteins (OBPs) (18) were annotated, which included three putative general OBP (GOBP), one more than typically reported for other Lepidoptera. To further characterize CpomGOBPs, we cloned cDNA copies of their transcripts and determined their expression patterns in various tissues. Cloning and sequencing of the 698 nt transcript for CpomGOBP1 resulted in the prediction of a 163 amino acid coding region, and subsequent RT-PCR indicated that the transcripts were mainly expressed in antennae and mouthparts. The 1 289 nt (160 amino acid) CpomGOBP2 and the novel 702 nt (169 amino acid) CpomGOBP3 transcripts are mainly expressed in antennae, mouthparts, and female abdomen tips. These results indicate that next generation sequencing is useful for the identification of novel transcripts of interest, and that codling moth expresses a transcript encoding for a new member of the GOBP subfamily.

Frequency of hybridization between Ostrinia nubilalis E-and Z-pheromone races in regions of sympatry within the United States

Female European corn borer, Ostrinia nubilalis, produce and males respond to sex pheromone blends with either E- or Z-Δ11-tetradecenyl acetate as the major component. E- and Z-race populations are sympatric in the Eastern United States, Southeastern Canada, and the Mediterranean region of Europe. The E- and Z-pheromone races of O. nubilalis are models for incipient species formation, but hybridization frequencies within natural populations remain obscure due to lack of a high-throughput phenotyping method. Lassance et al. previously identified a pheromone gland-expressed fatty-acyl reductase gene (pgfar) that controls the ratio of Δ11-tetradecenyl acetate stereoisomers. We identified three single nucleotide polymorphism (SNP) markers within pgfar that are differentially fixed between E- and Z-race females, and that are ≥98.2% correlated with female pheromone ratios measured by gas chromatography. Genotypic data from locations in the United States demonstrated that pgfar-z alleles were fixed within historically allopatric Z-pheromone race populations in the Midwest, and that hybrid frequency ranged from 0.00 to 0.42 within 11 sympatric sites where the two races co-occur in the Eastern United States (mean hybridization frequency or heterozygosity (H O) = 0.226 ± 0.279). Estimates of hybridization between the E- and Z-races are important for understanding the dynamics involved in maintaining race integrity, and are consistent with previous estimates of low levels of genetic divergence between E- and Z-races and the presence of weak prezygotic mating barriers. This work describes the development of new single nucleotide polymorphism (SNP) markers within the pheromone gland expressed fatty acyl reductase (pgfar) gene of Ostrinia nubilalis. These SNPs were shown to segregate based upon female pheromone production, and thus provide the first description of an assay for genetic determination of O. nubilalis pheromone strain from field-collected samples. These assays were applied to estimate hybridization within field populations, and represent valuable tools for future population genetic studies of this species.

Molecular genetics and genomics generate new insights into invertebrate pest invasions

Invertebrate pest invasions and outbreaks are associated with high social, economic, and ecological costs, and their significance will intensify with an increasing pressure on agricultural productivity as a result of human population growth and climate change. New molecular genetic and genomic techniques are available and accessible, but have been grossly underutilized in studies of invertebrate pest invasions, despite that they are useful tools for applied pest management and for understanding fundamental features of pest invasions including pest population demographics and adaptation of pests to novel and/or changing environments. Here, we review current applications of molecular genetics and genomics in the study of invertebrate pest invasions and outbreaks, and we highlight shortcomings from the current body of research. We then discuss recent conceptual and methodological advances in the areas of molecular genetics/genomics and data analysis, and we highlight how these advances will further our understanding of the demographic, ecological, and evolutionary features of invertebrate pest invasions. We are now well equipped to use molecular data to understand invertebrate dispersal and adaptation, and this knowledge has valuable applications in agriculture at a time when these are critically required.

De novo transcriptomic resources for two sibling species of moths: Ostrinia nubilalis and O. scapulalis

Background

This study aimed at enhancing the transcriptomic resources for two sibling species of moths, Ostrinia scapulalis (Adzuki bean borer) and Ostrinia nubilalis (European corn borer), as a foundation for future researches on their divergence history. Previous works on these species had shown that their genetic divergence was low, while they were reproductively isolated in natura and specialized on different host plants. Comparative genomic resources will help facilitate the understanding of the mechanisms involved in this isolation and adaptation to the host plants. Despite their fundamental interest, these species still lack the genomic resources to thoroughly identify candidate genes for functions of interest. We present here a high throughput sequencing and de novo transcriptome assembly for these two sibling species in line with this objective of comparative genomics.

Results

Based on 322,504 and 307,622 reads of 454 sequencing for O. scapulalis and O. nubilalis respectively, we reconstructed 11,231 and 10,773 transcripts, of which 40% were functionally annotated by BLAST analyzes. We determined the level of completeness of both assemblies as well as the recovery level of published Ostrinia genomic resources. Gene ontology (GO) of common and species-specific de novo transcripts did not reveal GO terms significantly enriched in one or the other species. By applying stringent homology searches on transcripts common to O. scapulalis and O. nubilalis, we identified a set of homologous transcripts, with a mean nucleotide identity value of 98.1%. In this set, the most divergent transcripts revealed candidate genes involved in developmental, sensorial and pathogen defense processes.

Conclusions

This data greatly increases the genomic resources of Ostrinia species and constitute a solid skeleton for future comparative analyzes of expression or diversity, despite we show that the transcriptomes for both species have not been assembled at full completion. In addition, we provide a set of homologous transcripts together with their annotation as a source of candidate genes for comparative analyzes.

Modeling evolution of resistance by Maruca vitrata (Lepidoptera: Crambidae) to transgenic insecticidal cowpea in Africa

We created a detailed model of the Maruca vitrata (F.) and cowpea [Vigna unguiculata (L.) Walp] system to study the possible evolution of resistance by the insect to transgenic insecticidal cowpea, which is under development. We focused on population dynamics and genetics in a region of west Africa. We simulated single-toxin and pyramided (two-toxin) cowpea and emphasized conservative, worst-case scenarios in our analysis. The results indicate that as long as a pyramided, transgenic cowpea can be developed, seed saving by farmers and reliance on natural refuge are not major problems for resistance management. Furthermore, it is possible that one or both toxins in the pyramid may not need to be high dose for evolution to be delayed significantly (>20 yr or 80 generations for resistance to become a concern if transgenic cowpea is deployed in areas where M. vitrata is endemic). If efforts are made to deploy transgenic cowpea only into the regions where M. vitrata is not endemic, then there is little to no concern with resistance emerging in the M. vitrata population.

Transcriptome analysis of the citrus red mite, Panonychus citri, and its gene expression by exposure to insecticide/acaricide

The citrus red mite, Panonychus citri, is known for its ability rapidly to evolve resistance to insecticides/acaricides and to adapt to hosts that produce toxins. In this study, we constructed an unprecedented four gigabase pair transcriptome of P. citri, which was assembled into 64 149 unique transcripts, the functions of which were annotated by five public databases. A total of 116 unique transcripts were identified as representatives of potential involvement in the detoxification of xenobiotics. Genes recorded to encoding insecticide/acaricide target proteins were also obtained from the P. citri transcriptome. In order to explore novel candidate genes potentially involved in the pesticide detoxification of P. citri, we also constructed digital gene expression libraries of short-term transcriptome responses of P. citri to pesticides, which resulted in the identification of 120 unique transcripts potentially associated with insecticide/acaricide detoxification. Our study will facilitate molecular research on pesticide resistance in citrus red mites, as well as in other phytophagous mites.