Sympatric ecological speciation meets pyrosequencing: sampling the transcriptome of the apple maggot Rhagoletis pomonella

The role of sexual isolation during rapid ecological divergence: Evidence for a new dimension of isolation in Rhagoletis pomonella

The pace of divergence and likelihood of speciation often depends on how and when different types of reproductive barriers evolve. Questions remain about how reproductive isolation evolves after initial divergence. We tested for the presence of sexual isolation (reduced mating between populations due to divergent mating preferences and traits) in Rhagoletis pomonella flies, a model system for incipient ecological speciation. We measured the strength of sexual isolation between two very recently diverged (~170 generations) sympatric populations, adapted to different host fruits (hawthorn and apple). We found that flies from both populations were more likely to mate within than between populations. Thus, sexual isolation may play an important role in reducing gene flow allowed by early-acting ecological barriers. We also tested how warmer temperatures predicted under climate change could alter sexual isolation and found that sexual isolation was markedly asymmetric under warmer temperatures - apple males and hawthorn females mated randomly while apple females and hawthorn males mated more within populations than between. Our findings provide a window into the early speciation process and the role of sexual isolation after initial ecological divergence, in addition to examining how environmental conditions could shape the likelihood of further divergence.

Role of Genes in Regulating Host Plants Expansion in Tephritid Fruit Flies (Diptera) and Potential for RNAi-Based Control

Host plant expansion is an important survival strategy for tephritids as they expand their range. Successful host expansion requires tephritids to adapt to the chemical and nonchemical properties of a novel host fruit, such as fruit color, phenology, and phytochemicals. These plant properties trigger a series of processes in tephritids, with each process having its own genetic basis, which means that various genes are involved in regulating host plant expansion by tephritids. This review summarizes current knowledge on the categories and roles of genes involved in host plant expansion in several important tephritid species, including genes related to chemoreception (olfactory and gustation), vision, digestion, detoxification, development, ribosomal and energy metabolism. Chemoreception- and detoxification- and digestion-related genes are stimulated by volatile chemicals and secondary chemicals of different hosts, respectively, which are involved in the regulation of nervous signal transduction that triggers behavioral, physical, and chemical responses to the novel host fruit. Vision-, nerve-, and development-related genes and metabolism-associated genes are activated in response to nonchemical stimuli from different hosts, such as color and phenology, to regulate a comprehensive adaptation of the extending host for tephritids. The chemical and nonchemical signals of hosts activate ribosomal and energy-related genes that result in the basic regulation of many processes of host expansion, including detoxification and development. These genes do not regulate novel host use individually, but multiple genes regulate multilevel adaptation to novel host fruits via multiple mechanisms. These genes may also be potential target genes for RNAi-based control of tephritid pests.

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.

Tephritid Fruit Fly Semiochemicals: Current Knowledge and Future Perspectives

The Dipteran family Tephritidae (true fruit flies) comprises more than 5000 species classified in 500 genera distributed worldwide. Tephritidae include devastating agricultural pests and highly invasive species whose spread is currently facilitated by globalization, international trade and human mobility. The ability to identify and exploit a wide range of host plants for oviposition, as well as effective and diversified reproductive strategies, are among the key features supporting tephritid biological success. Intraspecific communication involves the exchange of a complex set of sensory cues that are species- and sex-specific. Chemical signals, which are standing out in tephritid communication, comprise long-distance pheromones emitted by one or both sexes, cuticular hydrocarbons with limited volatility deposited on the surrounding substrate or on the insect body regulating medium- to short-distance communication, and host-marking compounds deposited on the fruit after oviposition. In this review, the current knowledge on tephritid chemical communication was analysed with a special emphasis on fruit fly pest species belonging to the Anastrepha, Bactrocera, Ceratitis, and Rhagoletis genera. The multidisciplinary approaches adopted for characterising tephritid semiochemicals, and the real-world applications and challenges for Integrated Pest Management (IPM) and biological control strategies are critically discussed. Future perspectives for targeted research on fruit fly chemical communication are highlighted.

A First Glimpse of the Mexican Fruit Fly Anastrepha ludens (Diptera: Tephritidae) Antenna Morphology and Proteome in Response to a Proteinaceous Attractant

Anastrepha ludens is a key pest of mangoes and citrus from Texas to Costa Rica but the mechanisms of odorant perception in this species are poorly understood. Detection of volatiles in insects occurs mainly in the antenna, where molecules penetrate sensillum pores and link to soluble proteins in the hemolymph until reaching specific odor receptors that trigger signal transduction and lead to behavioral responses. Scrutinizing the molecular foundation of odorant perception in A. ludens is necessary to improve biorational management strategies against this pest. After exposing adults of three maturity stages to a proteinaceous attractant, we studied antennal morphology and comparative proteomic profiles using nano-LC-MS/MS with tandem mass tags combined with synchronous precursor selection (SPS)-MS3. Antennas from newly emerged flies exhibited dense agglomerations of olfactory sensory neurons. We discovered 4618 unique proteins in the antennas of A. ludens and identified some associated with odor signaling, including odorant-binding and calcium signaling related proteins, the odorant receptor co-receptor (Orco), and putative odorant-degrading enzymes. Antennas of sexually immature flies exhibited the most upregulation of odor perception proteins compared to mature flies exposed to the attractant. This is the first report where critical molecular players are linked to the odor perception mechanism of A. ludens.

Gene Expression and Diet Breadth in Plant-Feeding Insects: Summarizing Trends

Transcriptomic studies lend insights into the role of transcriptional plasticity in adaptation and specialization. Recently, there has been growing interest in understanding the relationship between variation in herbivorous insect gene expression and the evolution of diet breadth. We review the studies that have emerged on insect gene expression and host plant use, and outline the questions and approaches in the field. Many candidate genes underlying herbivory and specialization have been identified, and a few key studies demonstrate increased transcriptional plasticity associated with generalist compared with specialist species. Addressing the roles that transcriptional variation plays in insect diet breadth will have important implications for our understanding of the evolution of specialization and the genetic and environmental factors that govern insect-plant interactions.

Use of DNA Markers for Grape Phylloxera Population and Evolutionary Genetics: From RAPDs to SSRs and Beyond

Grape phylloxera (Daktulosphaira vitifoliae Fitch) is a major pest of cultivated grapevines (Vitis spp.), occurring in virtually all viticultural regions around the world. Different grape phylloxera strains can be found at varying levels on leaves and roots on both own-rooted plants and in plants grafted onto partially resistant rootstocks. Considering its relevance for the adequate management of the pest in infested vineyards, the analysis of its genetic diversity has received considerable attention from the scientific community in the last decades. Here, we review 25 years of DNA-based molecular markers applied to the analysis of the genetic structure and the reproductive mode of grape phylloxera in its native range and in different introduced regions. The use given to RAPD, AFLP, mtDNA sequencing and microsatellite (SSR) genetic markers for the analysis of grape phylloxera diversity is discussed, and an overview of the main findings obtained after their application to different populations collected in diverse regions all around the world is shown. Lastly, we explore how recent advancements in molecular biology and in modern high throughput genotyping technologies may be applied to better understand grape phylloxera natural diversity at a molecular level.

Expression Patterns, Molecular Characterization, and Response to Host Stress of CYP Genes from Phenacoccus solenopsis (Hemiptera: Pseudococcidae)

The quarantine insect pest Phenacoccus solenopsis (Hemiptera: Pseudococcidae) has a broad host range and is distributed worldwide. Each year, P. solenopsis causes significant crop losses. The detoxification of various xenobiotic compounds involves the cytochrome P450 monooxygenase (CYP) superfamily of enzymes. However, the functions of CYPs in P. solenopsis are poorly understood. In the present study, P. solenopsis was reared from the egg to the adult stage on three host plants: Tomato, cotton, and hibiscus. Thirty-seven P. solenopsis CYP genes were identified and their phylogenetic relationships were analyzed. Eleven CYP genes (PsCYP4NT1, PsCYP4G219, PsCYP6PZ1, PsCYP6PZ5, PsCYP301B1, PsCYP302A1, PsCYP305A22, PsCYP315A1, PsCYP353F1, PsCYP3634A1, and PsCYP3635A2) were selected for quantitative real-time PCR analysis. The results demonstrated marked differences in CYP expression levels in P. solenopsis grown on different host plants. The results will aid the molecular characterization of CYPs and will increase our understanding of CYP expression patterns in P. solenopsis during development and growth on different hosts.

Transcriptome-Based Analysis Reveals a Crucial Role of BxGPCR17454 in Low Temperature Response of Pine Wood Nematode (Bursaphelenchus xylophilus)

Background: The causal agent of pine wilt disease is the pine wood nematode (PWN) (Bursaphelenchus xylophilus), whose ability to adapt different ecological niches is a crucial determinant of their invasion to colder regions. To discover the molecular mechanism of low temperature response mechanism, we attempted to study the molecular response patterns under low temperature from B. xylophilus with a comprehensive RNA sequencing analysis and validated the differentially expressed genes (DEGs) with quantitative real-time polymerase chain reaction (qRT-PCR). Bioinformatic software was utilized to isolate and identify the low-temperature-related BxGPCR genes. Transcript abundance of six low-temperature-related BxGPCR genes and function of one of the BxGPCR genes are studied by qRT-PCR and RNA interference. Results: The results showed that we detected 432 DEGs through RNA sequencing between low-temperature-treated and ambient-temperature-treated groups nematodes. The transcript level of 6 low-temperature-related BxGPCR genes increased at low temperature. And, the survival rates of BxGPCR17454 silenced B. xylophilus revealed a significant decrease at low temperature. Conclusion: in conclusion, this transcriptome-based study revealed a crucial role of BxGPCR17454 in low temperature response process of pine wood nematode. These discoveries would assist the development of management and methods for efficient control of this devastating pine tree pest.

Understanding the regulation of overwintering diapause molecular mechanisms in Culex pipiens pallens through comparative proteomics

To reveal overwintering dormancy (diapause) mechanisms of Culex pipiens pallens (L.), global protein expression differences at three separate time points represent nondiapause, diapause preparation and overwintering diapause phases of Cx. pipiens pallens were compared using iTRAQ. Cx. pipiens pallens females accumulate more lipid droplets during diapause preparation and overwintering diapause maintenance than during the nondiapause phase. A total of 1030 proteins were identified, among which 1020 were quantified and compared. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes (KEGG), Domain and Clusters of Orthologous Groups (COG) analyses revealed key groups of proteins, pathways and domains differentially regulated during diapause preparation and overwintering diapause maintenance phases in this mosquito, including major shifts in energy production and conversion, fatty acid metabolism, the citrate (TCA) cycle, and the cytoskeletal reorganization pathway. Our results provide novel insight into the molecular bases of diapause in mosquitoes and corroborate previously reported diapause-associated features in invertebrates. More interestingly, the phototransduction pathway exists in Cx. pipiens pallens, in particular, actin, rather than other proteins, appears to have substantial role in diapause regulation. In addition, the differential changes in calmodulin protein expression in each stage implicate its important regulatory role of the Cx. pipiens pallens biological clock. Finally, 24 proteins were selected for verification of differential expression using a parallel reaction monitoring strategy. The findings of this study provide a unique opportunity to explore the molecular modifications underlying diapause in mosquitoes and might therefore enable the future design and development of novel genetic tools for improving management strategies in mosquitoes.

Divergence in larval diapause induction between the rice and water-oat populations of the striped stem borer, Chilo suppressalis (Walker) (Lepidoptera: Crambidae)

Differences in diapause traits can result in the seasonal reproductive isolation of host plant-associated insect populations and thereby facilitate the population divergence. The striped stem borer, Chilo suppressalis, has two host plant-associated populations: rice population and water-oat population. Several studies have found evidence that seasonal reproductive isolation between these populations is at least partially due to interpopulation differences in diapause. However, there still lack unambiguous evidence comparing characteristics of diapause induction for both populations. We compared the photoperiodic response and the age of peak photoperiod sensitivity of these populations and used RNA-Seq to compare the molecular response of diapause induction between populations. The photoperiodic response of the two populations differed at 25 °C; the critical night length of larvae from the rice population was 11 h and 20 min, whereas no obvious critical night length was in those from the water-oat population. In rice population, larvae were most sensitive to photoperiod at 9-12 days of age, whereas in water-oat population, larvae were the most sensitive to photoperiod at 9-10 days of age. The RNA-Seq results indicated that there were several differences in the molecular response of diapause induction and small overlap in differentially expressed genes (DEGs) between populations. Furthermore, GO analysis indicated that both rice and water-oat population's DEGs were significantly enriched in heme and iron binding. Besides, water-oat population's DEGs were significantly enriched in metabolizing nutrients but rice population's DEGs do not. Thus, our results described differences in diapause induction between rice and water-oat populations of C. suppressalis which could affect the timing of diapause and thereby contribute to the seasonal reproductive isolation of these host plant-associated populations. In conclusion, this work suggests that difference in diapause induction could promote the population divergence in insects associated with different host plants.

Molecular characterization and evolution of a chemosensory receptor gene family in three notorious rice planthoppers, Nilaparvata lugens, Sogatella furcifera and Laodelphax striatellus, based on genome and transcriptome analyses

BACKGROUND

The white-backed planthopper (WBPH) Sogatella furcifera, the brown planthopper (BPH) Nilaparvata lugens, and the small brown planthopper (SBPH) Laodelphax striatellus (Hemiptera: Delphacidae) are rice pests that damage rice plants by sap-sucking and by transmitting viruses. Host-seeking behavior involves chemosensory receptor genes that include odorant receptors (ORs), ionotropic receptors (IRs) and gustatory receptors (GRs).

RESULTS

We used genome and transcriptome data to identify 141 ORs, 28 GRs and 25 IRs in BPH; 135 ORs, 18 GRs and 16 IRs in WBPH; and 37 ORs, 14 GRs and 6 IRs in SBPH. A phylogenetic analysis identified several specific OR clades of rice planthoppers, the results indicating that these OR members might be used to respond to specific host volatiles. OR co-receptor (Orco) is the most conserved and essential OR gene among these species and RNA interference (RNAi) can decrease their mRNA expression level to <50%. RNAi knockdown rice planthoppers were anosmia and were unable to seek or locate rice plants in behavioral tests.

CONCLUSION

The results demonstrate the importance of the planthopper Orco genes in locating rice plants. This information may aid in the development of RNAi-based transgenic rice and other pest management technologies. © 2018 Society of Chemical Industry.

Current source density mapping of antennal sensory selectivity reveals conserved olfactory systems between tephritids and Drosophila

Ecological specialization of insects involves the functional and morphological reshaping of olfactory systems. Little is known about the degree to which insect sensitivity to odorant compounds is conserved between genera, tribes, or families. Here we compared the olfactory systems of six tephritid fruit fly species spanning two tribes and the distantly related Drosophila melanogaster at molecular, functional, and morphological levels. Olfaction in these flies is mediated by a set of olfactory receptors (ORs) expressed in different functional classes of neurons located in distinct antennal regions. We performed a phylogenetic analysis that revealed both family-specific OR genes and putative orthologous OR genes between tephritids and Drosophila. With respect to function, we then used a current source density (CSD) analysis to map activity across antennae. Functional maps mirrored the intrinsic structure of antennae observed with scanning electron microscopy. Together, the results revealed partial conservation of the olfactory systems between tephritids and Drosophila. We also demonstrate that the mapping of olfactory responses is necessary to decipher antennal sensory selectivity to olfactory compounds. CSD analysis can be easily applied to map antennae of other species and therefore enables the rapid deriving of olfactory maps and the reconstructing of the target organisms' history of evolution.

Comparative transcriptome analysis reveals host-associated differentiation in Chilo suppressalis (Lepidoptera: Crambidae)

The striped stem borer, Chilo suppressalis Walker (Lepidoptera: Crambidae), is one of the most serious rice pests. Besides attacking rice, it also feeds on an economically important vegetable crop, water-oat Zizania latifolia. The species feeding on water-oat has higher growth and survival rate than those on rice, suggesting their success in adaptation to the new host plant. However, little is known about the molecular mechanisms of host plant adaptation. Here we investigated the midgut transcriptome responses of C. suppressalis larvae reared on rice and water-oat. A total of 1,633 differentially expressed genes were identified, with a greater number up-regulated on the more delicious new host. The up-regulation of most digestive and detoxification-related genes may be the result of adaptation to the changes in nutritional requirements and toxic chemicals during host shift. In contrast, down-regulation of ribosomal genes may be related to their better development performance when feeding on the new host. In conclusion, our results suggest that transcriptional regulation of genes related to digestion, detoxification and ribosome may play an important role in adaptation of C. suppressalis to a new host plant.

Odorant-binding proteins expression patterns in recently diverged species of Anastrepha fruit flies

We studied two species of closely related South American fruit flies, Anastrepha fraterculus and Anastrepha obliqua which, despite being able to interbreed, still show some ecological and reproductive differences. Because part of these differences, such as host and mate preferences, may be related to olfactory perception, we focused our investigation on the differential expression of Odorant-binding protein (OBP) gene family, which participate in initial steps of the olfactory signal transduction cascade. We investigated patterns of expression of eight OBP genes by qPCR in male and female head tissues of both species. The expression patterns of these OBPs suggest that some OBP genes are more likely involved with the location of food resources, while others seem to be associated with mate and pheromone perception. Furthermore, the expression patterns obtained at different reproductive stages indicate that OBP expression levels changed significantly after mating in males and females of both species. All eight OBP genes analyzed here showed significant levels of differential expression between A. fraterculus and A. obliqua, suggesting that they may hold important roles in their olfactory perception differences, and consequently, may potentially be involved in their differentiation.

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.

Tephritid Integrative Taxonomy: Where We Are Now, with a Focus on the Resolution of Three Tropical Fruit Fly Species Complexes

Accurate species delimitation underpins good taxonomy. Formalization of integrative taxonomy in the past decade has provided a framework for using multidisciplinary data to make species delimitation hypotheses more rigorous. We address the current state of integrative taxonomy by using as a case study an international project targeted at resolving three important tephritid species complexes: Bactrocera dorsalis complex, Anastrepha fraterculus complex, and Ceratitis FAR (C. fasciventris, C. anonae, C. rosa) complex. The integrative taxonomic approach has helped deliver significant advances in resolving these complexes: It has been used to identify some taxa as belonging to the same biological species as well as to confirm hidden cryptic diversity under a single taxonomic name. Nevertheless, the general application of integrative taxonomy has not been without issue, revealing challenges that must be considered when undertaking an integrative taxonomy project. Scrutiny of this international case study provides a unique opportunity to document lessons learned for the benefit of not only tephritid taxonomists, but also the wider taxonomic community.

Molecular evolution of Odorant-binding proteins gene family in two closely related Anastrepha fruit flies

BACKGROUND

Odorant-binding proteins (OBPs) are of great importance for survival and reproduction since they participate in initial steps of the olfactory signal transduction cascade, solubilizing and transporting chemical signals to the olfactory receptors. A comparative analysis of OBPs between closely related species may help explain how these genes evolve and are maintained under natural selection and how differences in these proteins can affect olfactory responses. We studied OBP genes in the closely related species Anastrepha fraterculus and A. obliqua, which have different host preferences, using data from RNA-seq cDNA libraries of head and reproductive tissues from male and female adults, aiming to understand the speciation process occurred between them.

RESULTS

We identified 23 different OBP sequences from Anastrepha fraterculus and 24 from A. obliqua, which correspond to 20 Drosophila melanogaster OBP genes. Phylogenetic analysis separated Anastrepha OBPs sequences in four branches that represent four subfamilies: classic, minus-C, plus-C and dimer. Both species showed five plus-C members, which is the biggest number found in tephritids until now. We found evidence of positive selection in four genes and at least one duplication event that preceded the speciation of these two species. Inferences on tertiary structures of putative proteins from these genes revealed that at least one positively selected change involves the binding cavity (the odorant binding region) in the plus-C OBP50a.

CONCLUSIONS

A. fraterculus and A. obliqua have a bigger OBP repertoire than the other tephritids studied, though the total number of Anastrepha OBPs may be larger, since we studied only a limited number of tissues. The contrast of these closely related species reveals that there are several amino acid changes between the homologous genes, which might be related to their host preferences. The plus-C OBP that has one amino acid under positive selection located in the binding cavity may be under a selection pressure to recognize and bind a new odorant. The other positively selected sites found may be involved in important structural and functional changes, especially ones in which site-specific changes would radically change amino acid properties.

BdorOBP83a-2 Mediates Responses of the Oriental Fruit Fly to Semiochemicals

The oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae), is one of the most destructive pests throughout tropical and subtropical regions in Asia. This insect displays remarkable changes during different developmental phases in olfactory behavior between sexually immature and mated adults. The olfactory behavioral changes provide clues to examine physiological and molecular bases of olfactory perception in this insect. We comparatively analyzed behavioral and neuronal responses of B. dorsalis adults to attractant semiochemicals, and the expression profiles of antenna chemosensory genes. We found that some odorant-binding proteins (OBPs) were upregulated in mated adults in association with their behavioral and neuronal responses. Ligand-binding assays further showed that one of OBP83a orthologs, BdorOBP83a-2, binds with high affinity to attractant semiochemicals. Functional analyses confirmed that the reduction in BdorOBP83a-2 transcript abundance led to a decrease in neuronal and behavioral responses to selected attractants. This study suggests that BdorOBP83a-2 mediates behavioral responses to attractant semiochemicals and could be a potential efficient target for pest control.

Transcriptome Analysis in Domesticated Species: Challenges and Strategies

Domesticated species occupy a special place in the human world due to their economic and cultural value. In the era of genomic research, domesticated species provide unique advantages for investigation of diseases and complex phenotypes. RNA sequencing, or RNA-seq, has recently emerged as a new approach for studying transcriptional activity of the whole genome, changing the focus from individual genes to gene networks. RNA-seq analysis in domesticated species may complement genome-wide association studies of complex traits with economic importance or direct relevance to biomedical research. However, RNA-seq studies are more challenging in domesticated species than in model organisms. These challenges are at least in part associated with the lack of quality genome assemblies for some domesticated species and the absence of genome assemblies for others. In this review, we discuss strategies for analyzing RNA-seq data, focusing particularly on questions and examples relevant to domesticated species.

Identification and preliminary characterization of chemosensory perception-associated proteins in the melon fly Bactrocera cucurbitae using RNA-seq

An investigation into proteins involved in chemosensory perception in the melon fly, Bactrocera cucurbitae (Diptera: Tephritidae) is described here using a newly generated transcriptome dataset. The melon fly is a major agricultural pest, widely distributed in the Asia-Pacific region and some parts of Africa. For this study, a transcriptome dataset was generated using RNA extracted from 4-day-old adult specimens of the melon fly. The dataset was assembled and annotated via Gene Ontology (GO) analysis. Based on this and similarity searches to data from other species, a number of protein sequences putatively involved in chemosensory reception were identified and characterized in the melon fly. This included the highly conserved "Orco" along with a number of other less conserved odorant binding protein sequences. In addition, several sequences representing putative ionotropic and gustatory receptors were also identified. This study provides a foundation for future functional studies of chemosensory proteins in the melon fly and for making more detailed comparisons to other species. In the long term, this will ultimately help in the development of improved tools for pest management.

Divergence of the diapause transcriptome in apple maggot flies: winter regulation and post-winter transcriptional repression

Duration of dormancy regulates seasonal timing in many organisms and may be modulated by day length and temperature. Though photoperiodic modulation has been well studied, temperature modulation of dormancy has received less attention. Here, we leverage genetic variation in diapause in the apple maggot fly, Rhagoletis pomonella, to test whether gene expression during winter or following spring warming regulates diapause duration. We used RNAseq to compare transcript abundance during and after simulated winter between an apple-infesting population and a hawthorn-infesting population where the apple population ends pupal diapause earlier than the hawthorn-infesting population. Marked differences in transcription between the two populations during winter suggests that the ‘early’ apple population is developmentally advanced compared to the ‘late’ hawthorn population prior to spring warming, with transcripts participating in growth and developmental processes relatively up-regulated in apple pupae during the winter cold period. Thus, regulatory differences during winter ultimately drive phenological differences that manifest themselves in the following summer. Expression and polymorphism analysis identify candidate genes in the Wnt and insulin signaling pathways that contribute to population differences in seasonality. Both populations remained in diapause and displayed a pattern of up- and then down-regulation (or vice versa) of growth-related transcripts following warming, consistent with transcriptional repression. The ability to repress growth stimulated by permissive temperatures is likely critical to avoid mismatched phenology and excessive metabolic demand. Compared to diapause studies in other insects, our results suggest some overlap in candidate genes/pathways, though the timing and direction of changes in transcription are likely species-specific.

Molecular Characterization and Differential Expression of an Olfactory Receptor Gene Family in the White-Backed Planthopper Sogatella furcifera Based on Transcriptome Analysis

The white-backed planthopper, Sogatella furcifera, a notorious rice pest in Asia, employs host plant volatiles as cues for host location. In insects, odor detection is mediated by two types of olfactory receptors: odorant receptors (ORs) and ionotropic receptors (IRs). In this study, we identified 63 SfurORs and 14 SfurIRs in S. furcifera based on sequences obtained from the head transcriptome and bioinformatics analysis. The motif-pattern of 130 hemiptera ORs indicated an apparent differentiation in this order. Phylogenetic trees of the ORs and IRs were constructed using neighbor-joining estimates. Most of the ORs had orthologous genes, but a specific OR clade was identified in S. furcifera, which suggests that these ORs may have specific olfactory functions in this species. Our results provide a basis for further investigations of how S. furcifera coordinates its olfactory receptor genes with its plant hosts, thereby providing a foundation for novel pest management approaches based on these genes.

Discovery of Chemosensory Genes in the Oriental Fruit Fly, Bactrocera dorsalis

The oriental fruit fly, Bactrocera dorsalis, is a devastating fruit fly pest in tropical and sub-tropical countries. Like other insects, this fly uses its chemosensory system to efficiently interact with its environment. However, our understanding of the molecular components comprising B. dorsalis chemosensory system is limited. Using next generation sequencing technologies, we sequenced the transcriptome of four B. dorsalis developmental stages: egg, larva, pupa and adult chemosensory tissues. A total of 31 candidate odorant binding proteins (OBPs), 4 candidate chemosensory proteins (CSPs), 23 candidate odorant receptors (ORs), 11 candidate ionotropic receptors (IRs), 6 candidate gustatory receptors (GRs) and 3 candidate sensory neuron membrane proteins (SNMPs) were identified. The tissue distributions of the OBP and CSP transcripts were determined by RT-PCR and a subset of nine genes were further characterized. The predicted proteins from these genes shared high sequence similarity to Drosophila melanogaster pheromone binding protein related proteins (PBPRPs). Interestingly, one OBP (BdorOBP19c) was exclusively expressed in the sex pheromone glands of mature females. RT-PCR was also used to compare the expression of the candidate genes in the antennae of male and female B. dorsalis adults. These antennae-enriched OBPs, CSPs, ORs, IRs and SNMPs could play a role in the detection of pheromones and general odorants and thus could be useful target genes for the integrated pest management of B. dorsalis and other agricultural pests.

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.

Sugarcane giant borer transcriptome analysis and identification of genes related to digestion

Sugarcane is a widely cultivated plant that serves primarily as a source of sugar and ethanol. Its annual yield can be significantly reduced by the action of several insect pests including the sugarcane giant borer (Telchin licus licus), a lepidopteran that presents a long life cycle and which efforts to control it using pesticides have been inefficient. Although its economical relevance, only a few DNA sequences are available for this species in the GenBank. Pyrosequencing technology was used to investigate the transcriptome of several developmental stages of the insect. To maximize transcript diversity, a pool of total RNA was extracted from whole body insects and used to construct a normalized cDNA database. Sequencing produced over 650,000 reads, which were de novo assembled to generate a reference library of 23,824 contigs. After quality score and annotation, 43% of the contigs had at least one BLAST hit against the NCBI non-redundant database, and 40% showed similarities with the lepidopteran Bombyx mori. In a further analysis, we conducted a comparison with Manduca sexta midgut sequences to identify transcripts of genes involved in digestion. Of these transcripts, many presented an expansion or depletion in gene number, compared to B. mori genome. From the sugarcane giant borer (SGB) transcriptome, a number of aminopeptidase N (APN) cDNAs were characterized based on homology to those reported as Cry toxin receptors. This is the first report that provides a large-scale EST database for the species. Transcriptome analysis will certainly be useful to identify novel developmental genes, to better understand the insect's biology and to guide the development of new strategies for insect-pest control.

De novo assembly and analysis of crow lungs transcriptome

The jungle crow (Corvus macrorhynchos) belongs to the order Passeriformes of bird species and is important for avian ecological and evolutionary genetics studies. However, there is limited information on the transcriptome data of this species. In the present study, we report the characterization of the lung transcriptome of the jungle crow using GS FLX Titanium XLR70. Altogether, 1,510,303 high-quality sequence reads with 581,198,230 bases was de novo assembled into 22,169 isotigs (isotig represents an individual transcript) and 784,009 singletons. Using these isotigs and 581,681 length-filtered (greater than 300 bp) singletons, 20,010 unique protein-coding genes were identified by BLASTx comparison against a nonredundant (nr) protein sequence database. Comparative analysis revealed that 46,604 (70.29%) and 51,642 (72.48%) of the assembled transcripts have significant similarity to zebra finch and chicken RefSeq proteins, respectively. As determined by GO annotation and KEGG pathway mapping, functional annotation of the unigenes recovered diverse biological functions and processes. Transcripts putatively involved in the immune response were identified. Furthermore, 20,599 single nucleotide polymorphisms (SNPs) and 7525 simple sequence repeats (SSRs) were retrieved from the assembled transcript database. This resource should lay an important base for future ecological, evolutionary, and conservation genetic studies on this species and in other related species.

How functional genomics will impact fruit fly pest control: the example of the Mediterranean fruit fly, Ceratitis capitata

The highly invasive agricultural insect pest Ceratitis capitata (Diptera: Tephritidae) is the most thoroughly studied tephritid fruit fly at the genetic and molecular levels. It has become a model for the analysis of fruit fly invasions and for the development of area-wide integrated pest management (AW-IPM) programmes based on the environmentally-friendly Sterile Insect Technique (SIT). Extensive transcriptome resources and the recently released genome sequence are making it possible to unravel several aspects of the medfly reproductive biology and behaviour, opening new opportunities for comparative genomics and barcoding for species identification. New genes, promotors and regulatory sequences are becoming available for the development/improvement of highly competitive sexing strains, for the monitoring of sterile males released in the field and for determining the mating status of wild females. The tools developed in this species have been transferred to other tephritids that are also the subject of SIT programmes.

An experimental and bioinformatics protocol for RNA-Seq analyses of photoperiodic diapause in the Asian tiger mosquito, Aedes albopictus

Photoperiodic diapause is an important adaptation that allows individuals to escape harsh seasonal environments via a series of physiological changes, most notably developmental arrest and reduced metabolism. Global gene expression profiling via RNA-Seq can provide important insights into the transcriptional mechanisms of photoperiodic diapause. The Asian tiger mosquito, Aedes albopictus, is an outstanding organism for studying the transcriptional bases of diapause due to its ease of rearing, easily induced diapause, and the genomic resources available. This manuscript presents a general experimental workflow for identifying diapause-induced transcriptional differences in A. albopictus. Rearing techniques, conditions necessary to induce diapause and non-diapause development, methods to estimate percent diapause in a population, and RNA extraction and integrity assessment for mosquitoes are documented. A workflow to process RNA-Seq data from Illumina sequencers culminates in a list of differentially expressed genes. The representative results demonstrate that this protocol can be used to effectively identify genes differentially regulated at the transcriptional level in A. albopictus due to photoperiodic differences. With modest adjustments, this workflow can be readily adapted to study the transcriptional bases of diapause or other important life history traits in other mosquitoes.

Transcriptome characterization analysis of Bactrocera minax and new insights into its pupal diapause development with gene expression analysis

Bactrocera minax is a major citrus pest distributed in China, Bhutan and India. The long pupal diapause duration of this fly is a major bottleneck for artificial rearing and underlying mechanisms remain unknown. Genetic information on B. minax transcriptome and gene expression profiles are needed to understand its pupal diapause. High-throughput RNA-seq technology was used to characterize the B. minax transcriptome and to identify differentially expressed genes during pupal diapause development. A total number of 52,519,948 reads were generated and assembled into 47,217 unigenes. 26,843 unigenes matched to proteins in the NCBI database using the BLAST search. Four digital gene expression (DGE) libraries were constructed for pupae at early diapause, late diapause, post-diapause and diapause terminated developmental status. 4,355 unigenes showing the differences expressed across four libraries revealed major shifts in cellular functions of cell proliferation, protein processing and export, metabolism and stress response in pupal diapause. When diapause was terminated by 20-hydroxyecdysone (20E), many genes involved in ribosome and metabolism were differentially expressed which may mediate diapause transition. The gene sets involved in protein and energy metabolisms varied throughout early-, late- and post-diapause. A total of 15 genes were selected to verify the DGE results through quantitative real-time PCR (qRT-PCR); qRT-PCR expression levels strongly correlated with the DGE data. The results provided the extensive sequence resources available for B. minax and increased our knowledge on its pupal diapause development and they shed new light on the possible mechanisms involved in pupal diapause in this species.

Identification of pheromone components and their binding affinity to the odorant binding protein CcapOBP83a-2 of the Mediterranean fruit fly, Ceratitis capitata

The Mediterranean fruit fly (or medfly), Ceratitis capitata (Wiedemann; Diptera: Tephritidae), is a serious pest of agriculture worldwide, displaying a very wide larval host range with more than 250 different species of fruit and vegetables. Olfaction plays a key role in the invasive potential of this species. Unfortunately, the pheromone communication system of the medfly is complex and still not well established. In this study, we report the isolation of chemicals emitted by sexually mature individuals during the "calling" period and the electrophysiological responses that these compounds elicit on the antennae of male and female flies. Fifteen compounds with electrophysiological activity were isolated and identified in male emissions by gas chromatography coupled to electroantennography (GC-EAG). Within the group of 15 identified compounds, 11 elicited a response in antennae of both sexes, whilst 4 elicited a response only in female antennae. The binding affinity of these compounds, plus 4 additional compounds known to be behaviourally active from other studies, was measured using C. capitata OBP, CcapOBP83a-2. This OBP has a high homology to Drosophila melanogaster OBPs OS-E and OS-F, which are associated with trichoid sensilla and co-expressed with the well-studied Drosophila pheromone binding protein LUSH. The results provide evidence of involvement of CcapOBP83a-2 in the medfly's odorant perception and its wider specificity for (E,E)-α-farnesene, one of the five major compounds in medfly male pheromone emission. This represents the first step in the clarification of the C. capitata and pheromone reception pathway, and a starting point for further studies aimed towards the creation of new powerful attractants or repellents applicable in the actual control strategies.

A pipeline for the de novo assembly of the Themira biloba (Sepsidae: Diptera) transcriptome using a multiple k-mer length approach

BACKGROUND

The Sepsidae family of flies is a model for investigating how sexual selection shapes courtship and sexual dimorphism in a comparative framework. However, like many non-model systems, there are few molecular resources available. Large-scale sequencing and assembly have not been performed in any sepsid, and the lack of a closely related genome makes investigation of gene expression challenging. Our goal was to develop an automated pipeline for de novo transcriptome assembly, and to use that pipeline to assemble and analyze the transcriptome of the sepsid Themira biloba.

RESULTS

Our bioinformatics pipeline uses cloud computing services to assemble and analyze the transcriptome with off-site data management, processing, and backup. It uses a multiple k-mer length approach combined with a second meta-assembly to extend transcripts and recover more bases of transcript sequences than standard single k-mer assembly. We used 454 sequencing to generate 1.48 million reads from cDNA generated from embryo, larva, and pupae of T. biloba and assembled a transcriptome consisting of 24,495 contigs. Annotation identified 16,705 transcripts, including those involved in embryogenesis and limb patterning. We assembled transcriptomes from an additional three non-model organisms to demonstrate that our pipeline assembled a higher-quality transcriptome than single k-mer approaches across multiple species.

CONCLUSIONS

The pipeline we have developed for assembly and analysis increases contig length, recovers unique transcripts, and assembles more base pairs than other methods through the use of a meta-assembly. The T. biloba transcriptome is a critical resource for performing large-scale RNA-Seq investigations of gene expression patterns, and is the first transcriptome sequenced in this Dipteran family.

Whole genome sequencing of the black grouse (Tetrao tetrix): reference guided assembly suggests faster-Z and MHC evolution

BACKGROUND

The different regions of a genome do not evolve at the same rate. For example, comparative genomic studies have suggested that the sex chromosomes and the regions harbouring the immune defence genes in the Major Histocompatability Complex (MHC) may evolve faster than other genomic regions. The advent of the next generation sequencing technologies has made it possible to study which genomic regions are evolutionary liable to change and which are static, as well as enabling an increasing number of genome studies of non-model species. However, de novo sequencing of the whole genome of an organism remains non-trivial. In this study, we present the draft genome of the black grouse, which was developed using a reference-guided assembly strategy.

RESULTS

We generated 133 Gbp of sequence data from one black grouse individual by the SOLiD platform and used a combination of de novo assembly and chicken reference genome mapping to assemble the reads into 4572 scaffolds with a total length of 1022 Mb. The draft genome well covers the main chicken chromosomes 1 ~ 28 and Z which have a total length of 1001 Mb. The draft genome is fragmented, but has a good coverage of the homologous chicken genes. Especially, 33.0% of the coding regions of the homologous genes have more than 90% proportion of their sequences covered. In addition, we identified ~1 M SNPs from the genome and identified 106 genomic regions which had a high nucleotide divergence between black grouse and chicken or between black grouse and turkey.

CONCLUSIONS

Our results support the hypothesis that the chromosome X (Z) evolves faster than the autosomes and our data are consistent with the MHC regions being more liable to change than the genome average. Our study demonstrates how a moderate sequencing effort can be combined with existing genome references to generate a draft genome for a non-model species.

Sniffing out chemosensory genes from the Mediterranean fruit fly, Ceratitis capitata

The Mediterranean fruit fly, Ceratitis capitata (medfly), is an extremely invasive agricultural pest due to its extremely wide host range and its ability to adapt to a broad range of climatic conditions and habitats. Chemosensory behaviour plays an important role in many crucial stages in the life of this insect, such as the detection of pheromone cues during mate pursuit and odorants during host plant localisation. Thus, the analysis of the chemosensory gene repertoire is an important step for the interpretation of the biology of this species and consequently its invasive potential. Moreover, these genes may represent ideal targets for the development of novel, effective control methods and pest population monitoring systems. Expressed sequence tag libraries from C. capitata adult heads, embryos, male accessory glands and testes were screened for sequences encoding putative odorant binding proteins (OBPs). A total of seventeen putative OBP transcripts were identified, corresponding to 13 Classic, three Minus-C and one Plus-C subfamily OBPs. The tissue distributions of the OBP transcripts were assessed by RT-PCR and a subset of five genes with predicted proteins sharing high sequence similarities and close phylogenetic affinities to Drosophila melanogaster pheromone binding protein related proteins (PBPRPs) were characterised in greater detail. Real Time quantitative PCR was used to assess the effects of maturation, mating and time of day on the transcript abundances of the putative PBPRP genes in the principal olfactory organs, the antennae, in males and females. The results of the present study have facilitated the annotation of OBP genes in the recently released medfly genome sequence and represent a significant contribution to the characterisation of the medfly chemosensory repertoire. The identification of these medfly OBPs/PBPRPs permitted evolutionary and functional comparisons with homologous sequences from other tephritids of the genera Bactrocera and Rhagoletis.

Ecological genomics of host shifts in Drosophila mojavensis

Advances in next-generation sequencing technologies have liberated our dependency on model laboratory species for answering genomic and transcriptomic level questions. These new techniques have dramatically expanded our breadth of study organisms and have allowed the analysis of species from diverse ecological environments. One such species is the cactophilic Drosophila mojavensis that inhabits the deserts of western North America. These insects feed and develop in the necrotic cacti, feeding largely on the microflora of the necrotic plant tissues. Drosophila mojavensis is composed of four geographically and ecologically separated populations. Each population (Baja California peninsula, mainland Sonoran Desert, Mojave Desert and Santa Catalina Island) utilizes the necrotic tissues of distinct cactus species. The differences in the nutritional and chemical composition of the necroses include a set of toxic compounds to which resident population must adapt. These ecological differences have facilitated many of the life history, behavior, physiological and genetic differences between the cactus host populations. Genomic resources have allowed investigators to examine the genomic and transcriptional level changes associated with the local adaptation of the four D. mojavensis populations, thereby providing further understanding of the genetic mechanism of adaptation and its role in the divergence of ecologically distinct populations.

No boundaries: genomes, organisms, and ecological interactions responsible for divergence and reproductive isolation

Revealing the genetic basis of traits that cause reproductive isolation, particularly premating or sexual isolation, usually involves the same challenges as most attempts at genotype-phenotype mapping and so requires knowledge of how these traits are expressed in different individuals, populations, and environments, particularly under natural conditions. Genetic dissection of speciation phenotypes thus requires understanding of the internal and external contexts in which underlying genetic elements are expressed. Gene expression is a product of complex interacting factors internal and external to the organism including developmental programs, the genetic background including nuclear-cytotype interactions, epistatic relationships, interactions among individuals or social effects, stochasticity, and prevailing variation in ecological conditions. Understanding of genomic divergence associated with reproductive isolation will be facilitated by functional expression analysis of annotated genomes in organisms with well-studied evolutionary histories, phylogenetic affinities, and known patterns of ecological variation throughout their life cycles. I review progress and prospects for understanding the pervasive role of host plant use on genetic and phenotypic expression of reproductive isolating mechanisms in cactophilic Drosophila mojavensis and suggest how this system can be used as a model for revealing the genetic basis for species formation in organisms where speciation phenotypes are under the joint influences of genetic and environmental factors.

A large-scale gene discovery for the red palm weevil Rhynchophorus ferrugineus (Coleoptera: Curculionidae)

The red palm weevil (RPW; Rhynchophorus ferrugineus) is a devastating pest of palms, prevalent in the Middle East as well as many other regions of the world. Here, we report a large-scale de novo complementary DNA (cDNA) sequencing effort that acquired ∼5 million reads and assembled them into 26 765 contigs from 12 libraries made from samples of different RPW developmental stages based on the Roche/454 GS FLX platform. We annotated these contigs based on the publically available known insect genes and the Tribolium castaneum genome assembly. We find that over 80% of coding sequences (CDS) from the RPW contigs have high-identity homologs to known proteins with complete CDS. Gene expression analysis shows that the pupa and larval stages have the highest and lowest expression levels, respectively. In addition, we also identified more than 60 000 single nucleotide polymorphisms and 1 200 simple sequence repeat markers. This study provides the first large-scale cDNA dataset for RPW, a much-needed resource for future molecular studies.

Microsatellite analyses across three diverse vertebrate transcriptomes (Acipenser fulvescens, Ambystoma tigrinum, and Dipodomys spectabilis)

Historically, many population genetics studies have utilized microsatellite markers sampled at random from the genome and presumed to be selectively neutral. Recent studies, however, have shown that microsatellites can occur in transcribed regions, where they are more likely to be under selection. In this study, we mined microsatellites from transcriptomes generated by 454-pyrosequencing for three vertebrate species: lake sturgeon (Acipenser fulvescens), tiger salamander (Ambystoma tigrinum), and kangaroo rat (Dipodomys spectabilis). We evaluated (i) the occurrence of microsatellites across species; (ii) whether particular gene ontology terms were over-represented in genes that contained microsatellites; (iii) whether repeat motifs were located in untranslated regions or coding sequences of genes; and (iv) in silico polymorphism. Microsatellites were less common in tiger salamanders than in either lake sturgeon or kangaroo rats. Across libraries, trinucleotides were found more frequently than any other motif type, presumably because they do not cause frameshift mutations. By evaluating variation across reads assembled to a given contig, we were able to identify repeat motifs likely to be polymorphic. Our study represents one of the first comparative data sets on the distribution of vertebrate microsatellites within expressed genes. Our results reinforce the idea that microsatellites do not always occur in noncoding DNA, but commonly occur in expressed genes.

Genomic divergence during speciation driven by adaptation to altitude

Even though Darwin's "On the Origin of Species" implied selection being the main driver of species formation, the role of natural selection in speciation remains poorly understood. In particular, it remains unclear how selection at a few genes can lead to genomewide divergence and the formation of distinct species. We used a particularly attractive clear-cut case of recent plant ecological speciation to investigate the demography and genomic bases of species formation driven by adaptation to contrasting conditions. High-altitude Senecio aethnensis and low-altitude S. chrysanthemifolius live at the extremes of a mountain slope on Mt. Etna, Sicily, and form a hybrid zone at intermediate altitudes but remain morphologically distinct. Genetic differentiation of these species was analyzed at the DNA polymorphism and gene expression levels by high-throughput sequencing of transcriptomes from multiple individuals. Out of ≈ 18,000 genes analyzed, only a small number (90) displayed differential expression between the two species. These genes showed significantly elevated species differentiation (FST and Dxy), consistent with diversifying selection acting on these genes. Genomewide genetic differentiation of the species is surprisingly low (FST = 0.19), while ≈ 200 genes showed significantly higher (false discovery rate < 1%; mean outlier FST > 0.6) interspecific differentiation and evidence for local adaptation. Diversifying selection at only a handful of loci may be enough for the formation and maintenance of taxonomically well-defined species, despite ongoing gene flow. This provides an explanation of why many closely related species (in plants, in particular) remain phenotypically and ecologically distinct despite ongoing hybridization, a question that has long puzzled naturalists and geneticists alike.

Is ecological speciation a major trend in aphids? Insights from a molecular phylogeny of the conifer-feeding genus Cinara

Introduction

In the past decade ecological speciation has been recognized as having an important role in the diversification of plant-feeding insects. Aphids are host-specialised phytophagous insects that mate on their host plants and, as such, they are prone to experience reproductive isolation linked with host plant association that could ultimately lead to species formation. The generality of such a scenario remains to be tested through macroevolutionary studies. To explore the prevalence of host-driven speciation in the diversification of the aphid genus Cinara and to investigate alternative modes of speciation, we reconstructed a phylogeny of this genus based on mitochondrial, nuclear and Buchnera aphidicola DNA sequence fragments and applied a DNA-based method of species delimitation. Using a recent software (PhyloType), we explored evolutionary transitions in host-plant genera, feeding sites and geographic distributions in the diversification of Cinara and investigated how transitions in these characters have accompanied speciation events.

Results

The diversification of Cinara has been constrained by host fidelity to conifer genera sometimes followed by sequential colonization onto different host species and by feeding-site specialisation. Nevertheless, our analyses suggest that, at the most, only half of the speciation events were accompanied by ecological niche shifts. The contribution of geographical isolation in the speciation process is clearly apparent in the occurrence of species from two continents in the same clades in relatively terminal positions in our phylogeny. Furthermore, in agreement with predictions from scenarios in which geographic isolation accounts for speciation events, geographic overlap between species increased significantly with time elapsed since their separation.

Conclusions

The history of Cinara offers a different perspective on the mode of speciation of aphids than that provided by classic models such as the pea aphid. In this genus of aphids, the role of climate and landscape history has probably been as important as host-plant specialisation in having shaped present-day diversity.

Ionotropic receptors (IRs): chemosensory ionotropic glutamate receptors in Drosophila and beyond

Ionotropic Receptors (IRs) are a recently characterized family of olfactory receptors in the fruit fly, Drosophila melanogaster. IRs are not related to insect Odorant Receptors (ORs), but rather have evolved from ionotropic glutamate receptors (iGluRs), a conserved family of synaptic ligand-gated ion channels. Here, we review the expression and function of IRs in Drosophila, highlighting similarities and differences with iGluRs. We also briefly describe the organization of the neuronal circuits in which IRs function, comparing and contrasting them with the sensory pathways expressing ORs. Finally, we summarize the bioinformatic identification and initial characterization of IRs in other species, which imply an evolutionarily conserved role for these receptors in chemosensation in insects and other protostomes.

Assessing De Novo transcriptome assembly metrics for consistency and utility

BACKGROUND

Transcriptome sequencing and assembly represent a great resource for the study of non-model species, and many metrics have been used to evaluate and compare these assemblies. Unfortunately, it is still unclear which of these metrics accurately reflect assembly quality.

RESULTS

We simulated sequencing transcripts of Drosophila melanogaster. By assembling these simulated reads using both a "perfect" and a modern transcriptome assembler while varying read length and sequencing depth, we evaluated quality metrics to determine whether they 1) revealed perfect assemblies to be of higher quality, and 2) revealed perfect assemblies to be more complete as data quantity increased.Several commonly used metrics were not consistent with these expectations, including average contig coverage and length, though they became consistent when singletons were included in the analysis. We found several annotation-based metrics to be consistent and informative, including contig reciprocal best hit count and contig unique annotation count. Finally, we evaluated a number of novel metrics such as reverse annotation count, contig collapse factor, and the ortholog hit ratio, discovering that each assess assembly quality in unique ways.

CONCLUSIONS

Although much attention has been given to transcriptome assembly, little research has focused on determining how best to evaluate assemblies, particularly in light of the variety of options available for read length and sequencing depth. Our results provide an important review of these metrics and give researchers tools to produce the highest quality transcriptome assemblies.

Assessing De Novo transcriptome assembly metrics for consistency and utility

Background

Transcriptome sequencing and assembly represent a great resource for the study of non-model species, and many metrics have been used to evaluate and compare these assemblies. Unfortunately, it is still unclear which of these metrics accurately reflect assembly quality.

Results

We simulated sequencing transcripts of Drosophila melanogaster. By assembling these simulated reads using both a “perfect” and a modern transcriptome assembler while varying read length and sequencing depth, we evaluated quality metrics to determine whether they 1) revealed perfect assemblies to be of higher quality, and 2) revealed perfect assemblies to be more complete as data quantity increased.

Several commonly used metrics were not consistent with these expectations, including average contig coverage and length, though they became consistent when singletons were included in the analysis. We found several annotation-based metrics to be consistent and informative, including contig reciprocal best hit count and contig unique annotation count. Finally, we evaluated a number of novel metrics such as reverse annotation count, contig collapse factor, and the ortholog hit ratio, discovering that each assess assembly quality in unique ways.

Conclusions

Although much attention has been given to transcriptome assembly, little research has focused on determining how best to evaluate assemblies, particularly in light of the variety of options available for read length and sequencing depth. Our results provide an important review of these metrics and give researchers tools to produce the highest quality transcriptome assemblies.

An EST database of the Caribbean fruit fly, Anastrepha suspensa (Diptera: Tephritidae)

Invasive tephritid fruit flies are a great threat to agriculture worldwide and warrant serious pest control measures. Molecular strategies that promote embryonic lethality in these agricultural pests are limited by the small amount of nucleotide sequence data available for tephritids. To increase the dataset for sequence mining, we generated an EST database by 454 sequencing of the caribfly, Anastrepha suspensa, a model tephritid pest. This database yielded 95,803 assembled sequences with 24% identified as independent transcripts. The percentage of caribfly sequences with hits to the closely related tephritid, Rhagoletis pomonella, transcriptome was higher (28%) than to Drosophila proteins/genes (18%) in NCBI. The database contained genes specifically expressed in embryos, genes involved in the cell death, sex-determination, and RNAi pathways, and transposable elements and microsatellites. This study significantly expands the nucleotide data available for caribflies and will be a valuable resource for gene isolation and genomic studies in tephritid insects.

Patterns of transcriptome divergence in the male accessory gland of two closely related species of field crickets

One of the central questions in evolutionary genetics is how much of the genome is involved in the early stages of divergence between populations, causing them to be reproductively isolated. In this article, we investigate genomic differentiation in a pair of closely related field crickets (Gryllus firmus and G. pennsylvanicus). These two species are the result of allopatric divergence and now interact along an extensive hybrid zone in eastern North America. Genes encoding seminal fluid proteins (SFPs) are often divergent between species, and it has been hypothesized that these proteins may play a key role in the origin and maintenance of reproductive isolation between diverging lineages. Hence, we chose to scan the accessory gland transcriptome to enable direct comparisons of differentiation for genes known to encode SFPs with differentiation in a much larger set of genes expressed in the same tissue. We have characterized differences in allele frequency between two populations for >6000 SNPs and >26,000 contigs. About 10% of all SNPs showed nearly fixed differences between the two species. Genes encoding SFPs did not have significantly elevated numbers of fixed SNPs per contig, nor did they seem to show larger differences than expected in their average allele frequencies. The distribution of allele frequency differences across the transcriptome is distinctly bimodal, but the relatively high proportion of fixed SNPs does not necessarily imply "ancient" divergence between these two lineages. Further studies of linkage disequilibrium and introgression across the hybrid zone are needed to direct our attention to those genome regions that are important for reproductive isolation.

Development of microsatellite markers using pyrosequencing in Galium trifidum (Rubiaceae), a rare species in Central Europe

We identify a large number of microsatellites from Galium trfidum, a plant species considered rare and endangered in Central and Western Europe. Using a combination of a total enriched genomic library and small-scale 454 pyrosequencing, we determined 9755 contigs with a length of 100 to 6192 bp. Within this dataset, we identified 153 SSR motifs in 144 contigs. Here, we tested 14 microsatellite loci in 2 populations of G. trifidum. The number of alleles and expected heterozygosity were 1–8 (mean 3.2) and 0.00–0.876 (0.549 on average), respectively. The markers described in this study will be useful for evaluating genetic diversity within and between populations, and gene flow between G. trifidum populations. These markers could also be applied to investigate the biological aspects of G. trifidum, such as the population dynamics and clonal structure, and to develop effective conservation programs for the Central European populations of this species.