Exploring Life History Choices: Using Temperature and Substrate Type as Interacting Factors for Blowfly Larval and Female Preferences

Blowflies (Diptera: Calliphoridae) present a wide range of larval lifestyles, typically classified as obligate parasitism, facultative parasitism, and complete sapro-necrophagy. Several parasitic species, both obligate and facultative, are considered to be of sanitary and economic importance, as their larvae can cause myiasis (maggot infestation in live tissue). However, it is noteworthy that the adult female plays a decisive role as she chooses the oviposition site, and, therefore, largely determines the feeding habit and developmental conditions of the larvae. In this study, two protocols are proposed to test larval feeding preference and female oviposition site preference considering two interacting factors: meat substrate type and temperature. The setups presented here allowed to test Lucilia cuprina larvae and gravid females in a four-choice assay with two temperatures (33 ± 2 °C and 25 ± 2 °C) and two types of meat substrates (fresh meat supplemented with blood and 5-day-old rotten meat). Larvae or gravid females can choose to burrow or lay their eggs, respectively, in either of the following: rotten meat at 25 °C (simulating a necrophagous species condition), fresh meat supplemented with blood at 33 °C (simulating a parasitic species condition), and two controls, rotten meat at 33 °C, or fresh meat supplemented with blood at 25 °C. The preference is assessed by counting the number of larvae or eggs laid in each option for each replicate. Comparing the observed results to a random distribution allowed for the estimation of the statistical significance of the preference. The results indicated that L. cuprina larvae have a strong preference for the rotten substrate at 25 °C. Conversely, oviposition-site preference by females was more varied for the meat type. This methodology can be adapted to test the preference of other insect species of similar size. Other questions can also be explored by using alternative conditions.

Biocatalytic potential of Pseudolycoriella CAZymes (Sciaroidea, Diptera) in degrading plant and fungal cell wall polysaccharides

Soil biota has a crucial impact on soil ecology, global climate changes, and effective crop management and studying the diverse ecological roles of dipteran larvae deepens the understanding of soil food webs. A multi-omics study of Pseudolycoriella hygida comb. nov. (Diptera: Sciaroidea: Sciaridae) aimed to characterize carbohydrate-active enzymes (CAZymes) for litter degradation in this species. Manual curation of 17,881 predicted proteins in the Psl. hygida genome identified 137 secreted CAZymes, of which 33 are present in the saliva proteome, and broadly confirmed by saliva CAZyme catalytic profiling against plant cell wall polysaccharides and pNP-glycosyl substrates. Comparisons with two other sciarid species and the outgroup Lucilia cuprina (Diptera: Calliphoridae) identified 42 CAZyme families defining a sciarid CAZyme profile. The litter-degrading potential of sciarids corroborates their significant role as decomposers, yields insights to the evolution of insect feeding habits, and highlights the importance of insects as a source of biotechnologically relevant enzymes.

Characterization of the mitochondrial genomes of Bradysia hygida, Phytosciara flavipes and Trichosia splendens (Diptera: Sciaridae) and novel insights on the control region of sciarid mitogenomes

Sciarids, also called "fungus gnats" are small, almost entirely dark-coloured insects. Sciarid larvae feed on different substrates and can infest agricultural crops and mushroom nurseries, causing economic losses. Of the 2174 Diptera mitogenome sequences currently available in GenBank, only eight are from the Sciaridae family, none of which are complete circular molecules. Here we describe the mitogenome sequences of three sciarid species: Phytosciara flavipes, Trichosia splendens and Bradysia hygida and provide novel insights on the control region of sciarid mitogenomes. The assembled mitogenomes range from 16,062 bp in P. flavipes to 17,095 bp in B. hygida. All 13 protein coding genes, 22 tRNAs and 2 rRNAs characteristic of insect mitogenomes were identified, but the sequence of the control region could not be determined. Experimental results suggest that the B. hygida control region is about 21 kb long resulting in a 37 kb long mitogenome which constitutes the largest insect mitochondrial genome described so far. Phylogenetic analysis using all Bibionomorpha mitogenome sequences available in GenBank strongly supports the Sciaridae monophyly and led to the identification of species and subfamily specific gene rearrangements. Our study extends the knowledge of this large and diverse insect family that includes agricultural pest species.

Evidence of positive selection on six spider developmental genes

Spiders constitute more than 49,000 described species distributed all over the world, and all ecological environments. Their order, Araneae, is defined by a set of characteristics with no parallel among their arachnid counterparts (e.g., spinnerets, silk glands, chelicerae that inoculate venom, among others). Changes in developmental pathways often underlie the evolution of morphological synapomorphies, and as such spiders are a promising model to study the role of developmental genes in the origin of evolutionary novelties. With that in mind, we investigated changes in the evolutionary regime of a set of six developmental genes, using spiders as our model. The genes were mainly chosen for their roles in spinneret ontogeny, yet they are pleiotropic, and it is likely that the origins of other unique morphological phenotypes are also linked to changes in their sequences. Our results indicate no great differences in the selective pressures on those genes when comparing spiders to other arachnids, but a few site-specific positive selection evidence were found in the Araneae lineage. These findings lead us to new insights on spider evolution that are to be further tested.

Genomic analysis on Brazilian strains of Anaplasma marginale

Anaplasma marginale is a vector-borne pathogen that causes a disease known as anaplasmosis. No sequenced genomes of Brazilian strains are yet available. The aim of this work was to compare whole genomes of Brazilian strains of A. marginale (Palmeira and Jaboticabal) with genomes of strains from other regions (USA and Australia strains). Genome sequencing of Brazilian strains was performed by means of next-generation sequencing. Reads were mapped using the genome of the Florida strain of A. marginale as a reference sequence. Single nucleotide polymorphisms (SNPs) and insertions/deletions (INDELs) were identified. The data showed that two Brazilian strains grouped together in one particular clade, which grouped in a larger American group together with North American strains. Moreover, some important differences in surface proteins between the two Brazilian isolates can be discerned. These results shed light on the evolutionary history of A. marginale and provide the first genome information on South American isolates. Assessing the genome sequences of strains from different regions is essential for increasing knowledge of the pan-genome of this bacteria.

Evolution of coding sequence and gene expression of blowflies and botflies with contrasting feeding habits

The Oestroidea superfamily is characterized by the diversity of feeding preferences among closely-related species; these flies are saprophagous, obligate parasites, or facultative parasites. We used gene expression and coding sequence data from five species (Cochliomyia hominivorax, Chrysomya megacephala, Lucilia cuprina, Dermatobia hominis, and Oestrus ovis) to identify underlying genetic differences involved in the diverse lifestyles. We tested whether 1287 orthologs have different expression and evolutionary constraints under different scenarios. We found two up-regulated genes; one in species causing cutaneous myiasis that is involved in iron transportation/metabolization (ferritin), and another in species causing traumatic myiasis that responds to reduced oxygen levels (anoxia up-regulated-like). Our evolutionary analysis showed a similar result. In the Co. hominivorax branch, we found one gene with the same function as ferritin that may be evolving under positive selection, spook. This is the first step towards understanding origins and evolution of parasitic strategy diversity in Oestroidea.

Traditional versus 3' RNA-seq in a non-model species

One limitation of the widely used RNA-seq method is that long transcripts are represented by more reads than shorter transcripts, resulting in a biased estimation of expression levels. The 3′ RNA-seq method, which yields only one sequence per transcript, bypasses this limitation. Here, RNA was extracted from two samples, in which we expected to find differentially expressed genes. Each was processed by both traditional and 3′ RNA-seq protocols. Both methods yielded similar differentially expressed genes and estimated expression levels in a comparable way, confirming they both represent valid tools for RNA-seq analysis. Notably, however, we identified more differentially expressed transcripts with the 3′ RNA-seq method, suggesting a greater power to detect expression variation using this method. Hence, when little genomic information is available for the species studied, the standard RNA-seq presents a better cost-benefit compromise, whereas for model species, the 3′ RNA-seq method might more accurately detect differential expression.

Evolution of genes involved in feeding preference and metabolic processes in Calliphoridae (Diptera: Calyptratae)

BACKGROUND

The genotype-phenotype interactions among traits governing feeding preference are of fundamental importance to behavioral genetics and evolutionary biology. The genetic basis of behavioral traits has been explored in different taxa using different approaches. However, the complex nature of the genetic mechanisms undergirding behavior is poorly understood. Here, we present an evolutionary study of candidate genes related to parasitism in Calliphoridae (Diptera: Calyptratae). Closely related species in this family exhibit distinct larval feeding habits, most notably necro-saprophagy and obligate parasitism.

METHODS

To understand the genetic and molecular bases underlying these habits, expression levels of eight candidate genes for feeding behavior-Cyp6g2, foraging, glutamate dehydrogenase, Jonah65aiv, Malvolio, PGRP-SC2, RPS6-p70-protein kinase, and smooth-were measured in four species using qPCR. Moreover we used expression values and sequence information to reconstruct the relationship among species and the d_N/d_S rate to infer possible sites under selection.

RESULTS

For most candidate genes, no statistically significant differences were observed, indicating a high degree of conservation in expression. However, Malvolio was differentially expressed between habits. Evolutionary analyses based on transcript levels and nucleotide sequences of Malvolio coding region suggest that transcript levels were correlated to feeding habit preferences among species, although deviations under a strictly neutral model were also observed in statistical tests.

DISCUSSION

Malvolio was the only gene demonstrating a possible connection to feeding habit. Differences in gene expression may be involved in (or be a result of) the genetic regulation of Calliphoridae feeding habit. Our results are the first steps towards understanding the genetic basis and evolution of feeding behavior in Calliphoridae using a functional approach.

Clustering of water bodies in unpolluted and polluted environments based on Escherichia coli phylogroup abundance using a simple interaction database

Different types of water bodies, including lakes, streams, and coastal marine waters, are often susceptible to fecal contamination from a range of point and nonpoint sources, and have been evaluated using fecal indicator microorganisms. The most commonly used fecal indicator is Escherichia coli, but traditional cultivation methods do not allow discrimination of the source of pollution. The use of triplex PCR offers an approach that is fast and inexpensive, and here enabled the identification of phylogroups. The phylogenetic distribution of E. coli subgroups isolated from water samples revealed higher frequencies of subgroups A1 and B23 in rivers impacted by human pollution sources, while subgroups D1 and D2 were associated with pristine sites, and subgroup B1 with domesticated animal sources, suggesting their use as a first screening for pollution source identification. A simple classification is also proposed based on phylogenetic subgroup distribution using the w-clique metric, enabling differentiation of polluted and unpolluted sites.

Selection and validation of reference genes for functional studies in the Calliphoridae family

The genera Cochliomyia and Chrysomya contain both obligate and saprophagous flies, which allows the comparison of different feeding habits between closely related species. Among the different strategies for comparing these habits is the use of qPCR to investigate the expression levels of candidate genes involved in feeding behavior. To ensure an accurate measure of the levels of gene expression, it is necessary to normalize the amount of the target gene with the amount of a reference gene having a stable expression across the compared species. Since there is no universal gene that can be used as a reference in functional studies, candidate genes for qPCR data normalization were selected and validated in three Calliphoridae (Diptera) species, Cochliomyia hominivorax Coquerel, Cochliomyia macellaria Fabricius, and Chrysomya albiceps Wiedemann . The expression stability of six genes ( Actin, Gapdh, Rp49, Rps17, α -tubulin, and GstD1) was evaluated among species within the same life stage and between life stages within each species. The expression levels of Actin, Gapdh, and Rp49 were the most stable among the selected genes. These genes can be used as reliable reference genes for functional studies in Calliphoridae using similar experimental settings.

Expression profiling of Drosophila mitochondrial genes via deep mRNA sequencing

Mitochondria play an essential role in several cellular processes. Nevertheless, very little is known about patterns of gene expression of genes encoded by the mitochondrial DNA (mtDNA). In this study, we used next-generation sequencing (NGS) for transcription profiling of genes encoded in the mitochondrial genome of Drosophila melanogaster and D. pseudoobscura. The analysis of males and females in both species indicated that the expression pattern was conserved between the two species, but differed significantly between both sexes. Interestingly, mRNA levels were not only different among genes encoded by separate transcription units, but also showed significant differences among genes located in the same transcription unit. Hence, mRNA abundance of genes encoded by mtDNA seems to be heavily modulated by post-transcriptional regulation. Finally, we also identified several transcripts with a noncanonical structure, suggesting that processing of mitochondrial transcripts may be more complex than previously assumed.

Molecular characterization of esterase E3 gene associated with organophosphorus insecticide resistance in the New World screwworm fly, Cochliomyia hominivorax

The New World screwworm, Cochliomyia hominivorax (Coquerel) (Diptera: Calliphoridae), is one of the most important myiasis-causing flies in South America. It is responsible for severe economic losses to livestock producers, mainly because it causes mortality in newborn calves and reductions in the quality of leather and in the production of milk and meat. The economic losses caused by myiasis, along with those caused by other internal and external parasites, are the main factors limiting meat production. In Brazil, C. hominivorax has been controlled by applying insecticides, particularly organophosphate (OP)-based compounds. However, the improper and continuous use of these chemicals can lead to the selection of OP-resistant strains. This, associated with the fast development of OP resistance in other myiasis-causing flies, shows the importance of investigating resistance in C. hominivorax. Based on the findings of previous studies, the objective of the current work was to isolate and sequence the E3 gene in C. hominivorax. Mutations at the positions (Gly137 and Trp251) responsible for conferring OP resistance in Lucilia cuprina and Musca domestica L. (Muscidae) were identified in C. hominivorax. In addition, the orthologous region in C. hominivorax contained motifs that are highly conserved among carboxyl/cholinesterases and contribute to the catalytic mechanism of the active site. The characterization of this gene in natural populations of New World screwworm can be an important tool for monitoring resistance to insecticides throughout its current geographic distribution. This will provide information for the selection and implementation of more effective pest management programmes.

Population genetics of New World screwworm from the Caribbean: insights from microsatellite data

Diseases affecting livestock can have a significant impact on animal productivity and on trade of live animals, meat and other animal products, which, consequently, affects the overall process of economic development. The New World screwworm, Cochliomyia hominivorax (Coquerel) (Diptera: Calliphoridae), is an important parasitic insect pest in Neotropical regions. This species has been successfully eradicated from North and most of Central America by the sterile insect technique, but continues to affect the development of the livestock sector in most Caribbean economies. Here, we provide some insight into the patterns of genetic variation and structure and gene flow of C. hominivorax populations from the Caribbean. Analysis of populations from 10 geographical sites in four islands revealed a moderate genetic variability within the populations. Surprisingly, a high population differentiation was found even in intra-island comparisons between populations. This observation can reflect either highly structured populations resulting from a lack of gene flow or a source-sink dynamic. Our study also suggests that New World screwworm populations can recover very rapidly from population contractions. This is valuable information that should be required prior to any investment in large-scale efforts aiming at controlling this pest.

PanGEA: identification of allele specific gene expression using the 454 technology

Background

Next generation sequencing technologies hold great potential for many biological questions. While mainly used for genomic sequencing, they are also very promising for gene expression profiling. Sequencing of cDNA does not only provide an estimate of the absolute expression level, it can also be used for the identification of allele specific gene expression.

Results

We developed PanGEA, a tool which enables a fast and user-friendly analysis of allele specific gene expression using the 454 technology. PanGEA allows mapping of 454-ESTs to genes or whole genomes, displaying gene expression profiles, identification of SNPs and the quantification of allele specific gene expression. The intuitive GUI of PanGEA facilitates a flexible and interactive analysis of the data. PanGEA additionally implements a modification of the Smith-Waterman algorithm which deals with incorrect estimates of homopolymer length as occuring in the 454 technology

Conclusion

To our knowledge, PanGEA is the first tool which facilitates the identification of allele specific gene expression. PanGEA is distributed under the Mozilla Public License and available at:

Gene expression profiling by massively parallel sequencing

Massively parallel sequencing holds great promise for expression profiling, as it combines the high throughput of SAGE with the accuracy of EST sequencing. Nevertheless, until now only very limited information had been available on the suitability of the current technology to meet the requirements. Here, we evaluate the potential of 454 sequencing technology for expression profiling using Drosophila melanogaster. We show that short (< approximately 80 bp) and long (> approximately 300-400 bp) cDNA fragments are under-represented in 454 sequence reads. Nevertheless, sequencing of 3' cDNA fragments generated by nebulization could be used to overcome the length bias of the 454 sequencing technology. Gene expression measurements generated by restriction analysis and nebulization for fragments within the 80- to 300-bp range showed correlations similar to those reported for replicated microarray experiments (0.83-0.91); 97% of the cDNA fragments could be unambiguously mapped to the genomic DNA, demonstrating the advantage of longer sequence reads. Our analyses suggest that the 454 technology has a large potential for expression profiling, and the high mapping accuracy indicates that it should be possible to compare expression profiles across species.

A survey of mutations in the Cochliomyia hominivorax (Diptera: Calliphoridae) esterase E3 gene associated with organophosphate resistance and the molecular identification of mutant alleles

Cochliomyia hominivorax (Calliphoridae) is one of the most important myiasis-causing flies and is responsible for severe economic losses to the livestock industry throughout the Neotropical region. In Brazil, C. hominivorax has been controlled mainly with organophosphate (OP) insecticides, although the inappropriate use of these chemicals can result in the selection of resistant flies. Changes in carboxylesterase activity have been associated with OP insecticides in some arthopodan species. In this work, we isolated and characterized part of the E3 gene in C. hominivorax (ChalphaE7), which contained the same substitutions responsible for the acquisition of OP hydrolase activity in Lucilia cuprina (Calliphoridae). Digestion of the polymerase chain reaction products with a restriction enzyme that specifically recognized the mutation site unambiguously differentiated wild and mutated esterase alleles. The PCR-RFLP assay therefore provided a fast, reliable DNA-based method for identifying C. hominivorax individuals with a mutation in the esterase gene. Further bioassays to determine the association of this mutation with OP resistance in C. hominivorax should allow the development of more effective strategies for managing this species.

The mitochondrial genome of the blowfly Chrysomya chloropyga (Diptera: Calliphoridae)

In view of the medical, sanitary and forensic importance of Chrysomya species, a knowledge of their nucleotide sequences would be useful for the molecular characterization of this genus, and would help in designing primers and in improving the molecular identification of Calliphoridae species. In this work, the mitochondrial genome of the blowfly Chrysomya chloropyga (Diptera: Calliphoridae) was completely sequenced. The entire mitochondrial DNA (mtDNA) molecule was 15,837 bp long and was sequenced using the shotgun approach. The overall nucleotide composition was heavily biased towards As and Ts, which accounted for 76.7% of the whole genome. The cox1 gene had a serine as the start codon, while incomplete termination codons mediated by tRNA signals were found for cox2, nd4 and nd5. The C. chloropyga genes were in the same order and orientation as the mitochondrial genome of other dipteran species, except for the occurrence of a 123 bp region that included a complete duplication of tRNA(Ile) and a partial duplication of tRNA(Gln) genes. C. chloropyga is the first species of Diptera with 23 tRNA genes instead of the usual 22 already described. A phylogenetic analysis showed a split of Brachycera into Calyptratae and Acalyptratae subdivisions. The complete sequence of C. chloropyga mtDNA described here will be a useful source of sequence information for general molecular and evolutionary studies in Diptera.