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.

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.

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.

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.

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.

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.

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.

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.

Molecular basis of resistance to organophosphate insecticides in the New World screw-worm fly.

BACKGROUND: The emergence of insecticide resistance is a fast-paced example of the evolutionary process of natural selection. In this study, we investigated the molecular basis of resistance in the myiasis-causing fly Cochliomyia hominivorax (Diptera: Calliphoridae) to dimethyl-organophosphate (OP) insecticides. METHODS: By sequencing the RNA from surviving larvae treated with dimethyl-OP (resistant condition) and non-treated larvae (control condition), we identified genes displaying condition-specific polymorphisms, as well as those differentially expressed. RESULTS: Both analyses revealed that resistant individuals have altered expression and allele-specific expression of genes involved in proteolysis (specifically serine-endopeptidase), olfactory perception and cuticle metabolism, among others. We also confirmed that resistant individuals carry almost invariably the Trp251Ser mutation in the esterase E3, known to confer OP and Pyrethroid resistance. Interestingly, genes involved in metabolic and detoxifying processes (notably cytochrome P450s) were found under-expressed in resistant individuals. An exception to this were esterases, which were found up-regulated. CONCLUSIONS: These observations suggest that reduced penetration and aversion to dimethyl-OP contaminated food may be important complementary strategies of resistant individuals. The specific genes and processes found are an important starting point for future functional studies. Their role in insecticide resistance merits consideration to better the current pest management strategies.

Microsatellite markers for population genetic studies of the blowfly Chrysomya putoria (Diptera: Calliphoridae).

The investigation of the genetic variation and population structure of Chrysomya species is of great interest for both basic and applied research. However, very limited genetic information is available for this genus across its geographical distribution. Here, we describe 12 polymorphic microsatellite loci isolated from Chrysomya putoria with expected heterozygosities ranging from 0.1402-0.8312. These markers are of potential applied interest for forensic entomologists and for the characterisation of the genetic structure of C. putoria from recently colonised regions, with great promise for understanding the colonisation dynamics and spread of the genus Chrysomya in the New World.

Detoxification mechanisms involved in ivermectin resistance in the cattle tick, Rhipicephalus (Boophilus) microplus.

The cattle tick Rhipicephalus microplus is one of the most important ectoparasites with great sanitary and economic impact for cattle rearing worldwide. Ivermectin is commonly used to control tick populations, but its use over the last 30 years has led to the development of resistant populations of R. microplus, and a concomitant loss of efficacy. In this context, we aimed to determine the metabolic mechanisms that contribute to ivermectin resistance in a resistant strain of this species. We performed lethal time bioassays with inhibitors of detoxifying enzymes and xenobiotic transporters (four detoxification pathways) using two strains of ticks: a susceptible strain, Mozo, and a resistant strain, Juarez. We used four inhibitors to test the involvement of different families of proteins responsible for detoxification of ivermectin, namely cytochrome P450, esterases, glutathione-S-transferase, and ATP Binding Cassette Transporters. We calculated the synergistic factor for each inhibitor and strain. To different degrees, all tested inhibitors altered the mortality rates in the strain Juarez, indicating that multiple mechanisms are responsible for the resistant phenotype. Detoxification mechanisms mediated by ABC transporters were observed to be the most important. Esterases, glutathione-S-transferases, and cytochrome-oxidases played less important roles in detoxification.

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.

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.

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 dN/dS 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.

Genomic analysis on Brazilian strains of Anaplasma marginale / Análise genômica de cepas Brasileiras de Anaplasma marginale

Abstract 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.

Resumo Anaplasma marginale é um patógeno transmitido por vetores que causam uma doença conhecida como anaplasmose. Até a presente data, não há genomas sequenciados de cepas brasileiras. O objetivo deste estudo foi comparar o genoma completo das cepas brasileiras de A. marginale (Palmeira e Jaboticabal) com os genomas de cepas de outras regiões (cepas dos EUA e Austrália). As sequências dos genomas das cepas brasileiras foram obtidas mediante sequenciamento de nova geração. As "reads" foram mapeadas usando-se como referência o genoma de A. marginale da cepa Florida. Foram identificados polimorfismos de nucleotídeo único (SNPs) e analisadas inserções/deleções (INDELs). As duas linhagens brasileiras se agruparam em um clado particular que, por sua vez, agrupou-se em um grupo maior junto com as linhagens norte-americanas. Além disso, foram identificadas diferenças significativas nas proteínas de superfície entre os dois isolados brasileiros. Esses resultados lançam luz sobre a história evolutiva de A. marginale e fornecem as primeiras informações de genomas de isolados sul-americanos. Avaliar as sequências de genomas de cepas de diferentes regiões é essencial para aumentar o conhecimento do pan-genoma dessa bactéria.

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.

Development of polymorphic microsatellite markers for the human botfly, Dermatobia hominis (Diptera: Oestridae).

In this report, we describe the development of 17 polymorphic microsatellite markers for the human botfly, Dermatobia hominis, an obligatory parasite of mammals of great veterinary importance in Latin America. The number of alleles ranged from 5 to 21 per locus, with a mean of 12.2 alleles per locus. The expected heterozygosity ranged from 0.2571 to 0.9206 and from 0.2984 to 0.9291 in two populations from Brazil. These markers should provide a high resolution tool for assessment of the fine-scale genetic structure of natural populations of the human botfly.

Characterization of polymorphic microsatellite markers for the blowfly Chrysomya albiceps (Diptera: Calliphoridae).

Chrysomya albiceps is a blowfly of great medical, sanitary and forensic importance widely distributed in the Afrotropical, southern Palaearctic, northern Oriental regions and, recently, in Central and South Americas. Here, we report the characterization of 13 polymorphic microsatellite markers for C. albiceps. The number of alleles ranged from three to 13 alleles with expected heterozygosities ranging from 0.4668 to 0.8408. These markers will be extremely useful for investigating many important aspects of this species such as population structure, dispersal and colonization dynamics.

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.