Identification and characterization of microRNAs in the screwworm flies Cochliomyia hominivorax and Cochliomyia macellaria (Diptera: Calliphoridae).

MicroRNAs (miRNAs) are small noncoding RNAs that modulate gene expression through post-transcriptional regulation. Here, we report the identification and characterization of miRNAs in two closely related screwworm flies with different feeding habits: Cochliomyia hominivorax and Cochliomyia macellaria. The New World screwworm, C. hominivorax, is an obligatory parasite of warm-blooded vertebrates, whereas the secondary screwworm, C. macellaria, is a free-living organism that feeds on decaying organic matter. Here, the small RNA transcriptomes of adults and third-instar larvae of both species were sequenced. A total of 110 evolutionarily conserved miRNAs were identified, and 10 putative precursor miRNAs (pre-miRNAs) were predicted. The relative expression of six selected miRNAs was further investigated, including miRNAs that are related to reproduction and neural processes in other insects. Mature miRNAs were also characterized across an evolutionary time scale, suggesting that the majority of them have been conserved since the emergence of the Arthropoda [540 million years ago (Ma)], Hexapoda (488 Ma) and Brachycera (195 Ma) lineages. This study is the first report of miRNAs for screwworm flies. We also performed a comparative analysis with the hereby predicted miRNAs from the sheep blowfly, Lucilia cuprina. The results presented may advance our understanding of parasitic habits within Calliphoridae and assist further functional studies in blowflies.

Detection and genetic diversity of a heliothine invader (Lepidoptera: Noctuidae) from north and northeast of Brazil.

The cotton bollworm, Helicoverpa armigera (Hübner), was recently introduced in Brazil. During the 2012-2013 harvest, producers reported reduced yields up to 35% on major crops. The economic losses reached US$ 1 billion only in western Bahia, triggering a phytosanitary crisis. The deficiencies in existing taxonomic keys to deal with the morphologically indistinct larvae of H. armigera and the native Helicoverpa zea (Boddie) constrained the detection of new incursions of this heliothine invader. This study explored the identity of heliothine larvae that were found infesting soybean- and corn-growing areas from Roraima state, northern Brazil, through sequences of the mitochondrial cytochrome c oxidase subunit I gene. The inter- and intraspecies sequence variations of DNA barcodes in H. armigera and H. zea were analyzed. The genetic diversity and population structure of the specimens from Roraima and two populations from Piauí and Bahia states, northeastern Brazil, were assessed by adding the cytochrome c oxidase subunit II gene to the analysis. Owing to the lack of studies on genetic introgression for the two species, the suitability of using three different nuclear genes to distinguish the two species was also investigated. The results showed strong evidence that the heliothine larvae from north and northeast of Brazil are conspecific with H. armigera, suggesting that this invasive moth has already crossed the Amazon basin. Surveys in the north of South America should start as soon as possible to monitor the entry or spread of this moth in the Caribbean, Central America, and the United States.

Structural characterization of the internal transcribed spacer 2 (ITS2) of the ribosomal DNA (rDNA) cluster in Calyptratae (Diptera: Schizophora) and its implications for molecular phylogenetic analyses.

The internal transcribed spacer 2 (ITS2) of the eukaryotic ribosomal DNA (rDNA) cluster plays an essential role in processing of the ribosomal RNA, which is primarily accomplished by the secondary structures acquired by the molecule after transcription. Two possible structural conformation models have been proposed for the ITS2 region, the "ring model" and the "hairpin model," and the former has been widely used in many molecular phylogenetic analyses incorporating structural information available to date. To evaluate the validity of this model, in vitro transcribed ITS2 molecules from species representing the three superfamilies of the Calyptratae clade (Diptera: Schizophora), namely Cochliomyia hominivorax, Musca domestica, and Glossina morsitans, were submitted to enzymatic digestion with single- and double-stranded specific nucleases (RNases I, A, T1, and V1). The resulting fragments were analyzed by capillary electrophoresis and digestion sites were mapped in the secondary structure models which were obtained by in silico prediction with further refinement by homology comparisons. The pattern of RNA fragments generated by these RNases show a high degree of correlation to most of the predicted helix-loop regions and structural motifs. Discrepancies to the models can be explained by alternative structural conformation dynamics (in M. domestica and G. morsitans) and by higher-order factors (such as tertiary interactions) that may stabilize thermodynamically unfavored structures (in C. hominivorax).

Molecular phylogenetics of Oestroidea (Diptera: Calyptratae) with emphasis on Calliphoridae: insights into the inter-familial relationships and additional evidence for paraphyly among blowflies.

The superfamily Oestroidea, comprising ∼15,000 species, is a large and ecologically diverse clade within the order Diptera. Among its six commonly recognized families, Calliphoridae seems to be crucial for understanding evolutionary relationships in the group, as it is recognized as a controversial paraphyletic grouping. To further investigate this matter, the ITS2, 28S, COI and 16S regions were used to infer phylogenetic relationships in Oestroidea with maximum-parsimony (MP), maximum-likelihood (ML) and Bayesian inference (BI) methods. For the BI analyses, a deep evaluation of different data partitioning strategies was conducted, including consideration of structural conformation (ITS2 and 16S) and codon position (COI) information. Results suggest the existence of two main clades in Oestroidea: (Tachinidae+Mesembrinellinae) and (Rhiniinae, (Sarcophagidae+Calliphoridae sensu stricto)). Oestridae was recovered as sister group of the remaining Oestroidea in the MP trees while it was placed closer to the (Rhiniinae+Sarcophagidae+Calliphoridae sensu stricto) group in the ML and BI trees. A paraphyletic Calliphoridae was recovered, confirming the exclusion of Rhiniinae, a clade recently promoted to family status and therefore already excluded. Mesembrinellinae could also be considered a distinct group apart from Calliphoridae, although further studies are required. Consideration of structural and codon position information led to a significant increase in the log-likelihoods of the analyses, which were accompanied by small changes in the inferred topologies, branch lengths and posterior probability support values. However, as model complexity increases, so does uncertainty across the estimated parameters, including tree topologies, and phylogenies inferred under very parameter-rich models may be less reliable even when possessing higher log-likelihoods.

Exposure to the COVID-19 pandemic environment and generosity.

We report data from an online experiment which allows us to study how generosity changed over a 6-day period during the initial explosive growth of the COVID-19 pandemic in Andalusia, Spain, while the country was under a strict lockdown. Participants (n = 969) could donate a fraction of a €100 prize to an unknown charity. Our data are particularly rich in the age distribution and we complement them with daily public information about COVID-19-related deaths, infections and hospital admissions. We find correlational evidence that donations decreased in the period under study, particularly among older individuals. Our analysis of the mechanisms behind the detected decrease in generosity suggests that expectations about others' behaviour, perceived mortality risk and (alarming) information play a key-but independent-role for behavioural adaptation. These results indicate that social behaviour is quickly adjusted in response to the pandemic environment, possibly reflecting some form of selective prosociality.

Evaluation of the internal transcribed spacer 2 (ITS2) as a molecular marker for phylogenetic inference using sequence and secondary structure information in blow flies (Diptera: Calliphoridae).

The internal transcribed spacer 2 (ITS2) is a small non-coding region located inside the nuclear ribosomal DNA cluster. ITS2 sequence variability is thought to be appropriate to differentiate species and for phylogenetic reconstructions analyses, which can be further improved if structural information is considered. We evaluated the potential of ITS2 as a molecular marker for phylogenetic inference in Calliphoridae (Diptera: Brachycera) using a broad range of inference methods and different substitution models, accounting or not for structural information. Sequence analyses revealed a hierarchically organized pattern of sequence variation and a small level of nucleotide substitution saturation. Intragenomic variation due to small sequence repeats was found mainly in the most variable domain (IV), but it has no significant impact on the phylogenetic signal at the species level. Inferred secondary structures revealed that GC pairs are more frequently found flanking bulges and loops regions in more conserved domains, thus ensuring structure stability. In the phylogenetic analyses, the use of substitution models accounting for structural information significantly improves phylogenetic inference in both neighbour-joining and Bayesian analyses, although the former provides limited resolution for dealing with highly divergent sequences. For Bayesian analyses, a significant improvement in likelihood was observed when considering structure information, although with small changes in topology and overall support, probably reflecting better evolutionary rates estimates. Based on these findings, ITS2 is a suitable molecular marker for phylogenetic analyses in Calliphoridae, at both species and generic level.

Wing morphometry as a tool for correct identification of primary and secondary New World screwworm fly.

Cochliomyia hominivorax and Cochliomyia macellaria are endemic Neotropical Calliphoridae species. The former causes severe myiasis in hosts while the latter is Sarcosaprophagous, but commonly found as a second invader in wounds. Due to the morphological similarity between them and the potential losses that C. hominivorax represents for cattle breeders, the rapid and correct identification of these two species is very important. In addition to a correct identification of these species, a good knowledge of C. hominivorax biology can be helpful for designing control programs. We applied geometric morphometric methods to assess wing differences between C. hominivorax and C. macellaria and conduct a preliminary analysis of wing morphological variation in C. hominivorax populations. Canonical variate analysis, using wing shape data, correctly classified 100% of the individuals analyzed according to sex and species. This result demonstrates that wing morphometry is a simple and reliable method for identifying C. hominivorax and C. macellaria samples and can be used to monitor C. hominivorax. Both species show sexual dimorphism, but in C. hominivorax it is magnified. We suggest that this may reflect different histories of selection pressures operating on males and females. Significant differences in wing size and shape were obtained among C. hominivorax populations, with little correlation with latitude. This result suggests that wing variation is also a good morphological marker for studying population variation in C. hominivorax.

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.

Complex patterns of genetic variability in populations of the New World screwworm fly revealed by mitochondrial DNA markers.

Cochliomyia hominivorax (Coquerel) (Diptera: Calliphoridae), the New World screwworm fly, is an important agent of traumatic myiasis, which is endemic in the Neotropical region and which has great economic impact on the development of the livestock industry. International efforts have been aimed at designing programmes to control and eradicate this species from endemic areas. Thorough knowledge of the population genetics of an insect pest is a fundamental component to ensuring the success of a pest management strategy because it enables the determination of an appropriate geographic scale for carrying out effective treatments. This study undertook an analysis of mtDNA polymerase chain reaction-restricted fragment length polymorphism (PCR-RFLP) in 34 populations of C. hominivorax from 10 countries, encompassing almost all the current distribution of the species. Results showed high levels of mitochondrial DNA variability (pi= 2.9%) and a complex pattern of population genetic structure for this species. Significant population structure (Phi st= 0.5234) and low variability were found in Caribbean populations, suggesting that, in general, islands constitute independent evolutionary entities connected by restricted gene flow. By contrast, high variability and low, but significant, differentiation was found among mainland populations (Phi st= 0.0483), which could not be attributed to geographic distance. Several processes may be acting to maintain the observed patterns, with different implications for establishing control programmes.

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.

Photographic map of the polytene chromosomes of Cochliomyia hominivorax.

Cochliomyia hominivorax (Coquerel) (Diptera: 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. A polytene chromosome map is an invaluable tool for the genetic analysis and manipulation of any species because it allows the integration of physical and genetic maps. Cochliomyia hominivorax has a diploid number of 12 chromosomes (2n = 12): five pairs of autosomes and one pair of sex chromosomes (XX/XY), which do not polytenize. We created a new photomap of the polytene chromosomes of C. hominivorax describing its five autosomes (chromosomes 2-6). Pupal trichogen cells, which have chromosomes with a high degree of polytenization, were used to elaborate this map. The photomap was made by comparing 20 different nuclei and choosing, for each chromosome segment, the region with the highest resolution. Thus, we present a new photomap of the five autosomes of this species, with a total resolution of 1450 bands.

Investigation of mutations associated with pyrethroid resistance in populations of the New World Screwworm fly, Cochliomyia hominivorax (Diptera: Calliphoridae).

Cochliomyia hominivorax larvae are known for their parasitic habit in living vertebrates, causing considerable economic losses to livestock industry. This ectoparasite has been controlled mainly by applying insecticides, but this method usually results in the selection of resistant individuals. The resistance mechanism known as knockdown resistance (kdr) is a generic term for amino acid substitutions in the sodium channel associated with pyrethroid resistance, and substitutions in residue 251 of the carboxylesterase E3 have been associated with organophosphate and pyrethroid hydrolysis. We looked for L1014F kdr and W251S mutations in the sodium channel and E3 genes, respectively, in C. hominivorax populations. Ten populations obtained from its current distribution were investigated using the polymerase chain reaction-restriction fragment length polymorphism technique. No mutant individuals were found for the kdr mutation. However, the W251S mutation was found in all populations investigated. Results from a bioassay with cypermethrin (a pyrethroid) indicated that the survival at the lowest concentration (Fisher exact test, P = 0.0003) and an intermediate concentration (P = 0.024) were associated with the W251S mutation. This correlation was not significant at the highest concentration tested (P = 0.221). We suggest that the W251S mutation in the C. hominivorax carboxylesterase E3 is also associated with pyrethroid hydrolysis. This information may contribute to the elaboration of improved management programs for this ectoparasite.

Registered Replication Report: Rand, Greene, and Nowak (2012).

In an anonymous 4-person economic game, participants contributed more money to a common project (i.e., cooperated) when required to decide quickly than when forced to delay their decision (Rand, Greene & Nowak, 2012), a pattern consistent with the social heuristics hypothesis proposed by Rand and colleagues. The results of studies using time pressure have been mixed, with some replication attempts observing similar patterns (e.g., Rand et al., 2014) and others observing null effects (e.g., Tinghög et al., 2013; Verkoeijen & Bouwmeester, 2014). This Registered Replication Report (RRR) assessed the size and variability of the effect of time pressure on cooperative decisions by combining 21 separate, preregistered replications of the critical conditions from Study 7 of the original article (Rand et al., 2012). The primary planned analysis used data from all participants who were randomly assigned to conditions and who met the protocol inclusion criteria (an intent-to-treat approach that included the 65.9% of participants in the time-pressure condition and 7.5% in the forced-delay condition who did not adhere to the time constraints), and we observed a difference in contributions of -0.37 percentage points compared with an 8.6 percentage point difference calculated from the original data. Analyzing the data as the original article did, including data only for participants who complied with the time constraints, the RRR observed a 10.37 percentage point difference in contributions compared with a 15.31 percentage point difference in the original study. In combination, the results of the intent-to-treat analysis and the compliant-only analysis are consistent with the presence of selection biases and the absence of a causal effect of time pressure on cooperation.

Improving access to the control region and tRNA gene clusters of dipteran mitochondrial DNA.

Despite the successful use of universal primers for amplifying insect mtDNA, specific regions remain difficult to recover and demand the use of taxon-specific primers. In this work, we describe a new set of primers for efficiently amplifying and sequencing the mtDNA control region and three tRNA gene clusters of dipterans of medical and veterinary importance, including Muscidae, Calliphoridae, and Oestridae species. These new primers were useful for investigating the nucleotide information and the structural organization of dipteran mtDNA.

Primary myiasis in dog caused by Phaenicia eximia (Diptera:Calliphoridae) and preliminary mitochondrial DNA analysis of the species in Brazil.

A case of primary myiasis in a dog caused by Phaenicia eximia (Robineau-Desvoidy) in Brazil is presented. A young and healthy female dog, Canis familiaris, approximately 10 d old and still under maternal care, was found to have several eggs and 1st instars larvae in its abdomen and urogenital regions. Samples were collected in Campinas, São Paulo, and transferred to the laboratory for rearing and identification. A comparative analysis of the mitochondrial DNA (mtDNA) with 12 restriction enzymes in 2 sampled populations of P. eximia collected in different hosts (live dog and bovine carcass) and in the same locality revealed that 4, EcoRI, EcoRV, HaeIII, and MspI were suitable for detecting mtDNA markers in the 2 populations.

New evidence on the systematic and phylogenetic position of Parides burchellanus (Lepidoptera: Papilionidae).

Parides burchellanus is considered a rare and threatened swallowtail species restricted to central Brazil. It shows considerable morphological similarity to Parides panthonus jaguarae, with which it shares both geographical range and larval host plant. At present, P. burchellanus and P. panthonus are believed to be distinct species, based on minor differences in male genitalia. In this study, the phylogenetic and systematic position of P. burchellanus in relation to three subspecies of P. panthonus (P. p. jaguarae, P. p. lysimachus and P. p. aglaope) was evaluated using molecular evidence: the complete sequence of the mtDNA genes COI and COII and of the nuclear gene EF-1α (c. 3300 bp). In addition, the informativeness of the 'barcode' region next to the 5' end of COI (c. 650 bp) was evaluated for delimiting these taxa. Individual analysis by neighbour-joining, using Kimura 2-parameter distance model, and by maximum parsimony showed that P. p. jaguarae + P. p. lysimachus + P. p. aglaope + P. burchellanus form a strongly supported monophyletic clade, and all molecular regions consistently recovered P. p. jaguarae and P. burchellanus as sister species. The genetic divergence among the subspecies of P. panthonus and P. burchellanus is equivalent to the divergence among conspecifics of other species of Parides, and smaller than the interspecific divergence among different sister species of this genus. The results support the proposal that P. p. jaguarae and P. burchellanus are likely to be synonymous, and suggest that P. burchellanus can be considered conspecific with P. panthonus. The taxonomic classification of P. burchellanus should be revised on the basis of the molecular data.

Genetic variability in mitochondrial DNA of the screwworm, Cochliomyia hominivorax (Diptera: Calliphoridae), from Brazil.

Restriction fragment length polymorphism (RFLP) analysis of mitochondrial DNA (mtDNA) was used to examine genetic variation and population structure of screwworm flies in four populations from São Paulo State, Brazil. The total DNA of 405 individuals was digested with 15 restriction endonucleases and probed with five cloned HindIII fragments representing the entire mitochondrial genome of Cochliomyia hominivorax. The survey revealed that four enzymes (HaeIII, HindIII, MspI, and PvuII) were suitable to detect mtDNA variation among all populations. Based on the fragment patterns obtained for these four enzymes, a total of 15 haplotypes in combination was detected. Heteroplasmic individuals for the PvuII pattern were obtained in one of the populations. The estimated average for nucleotide sequence divergence (delta) was 0.92%. The cladogram of the geographical distribution among the observed haplotypes suggests that the sampled screwworms probably belong to a single evolutionary lineage with populations interconnected by reduced gene flow.

Evolution and structural organisation of mitochondrial DNA control region of myiasis-causing flies.

This study reports the molecular characterization of the mtDNA control region (called the A+T-rich region in insects) of five dipteran species which cause myiasis: Cochliomyia hominivorax Coquerel, Cochliomyia macellaria Fabricius, Chrysomya megacephala Fabricius, Lucilia eximia Wiedemann (Diptera: Calliphoridae) and Dermatobia hominis Linnaeus Jr (Diptera: Oestridae). The control region in these species varies in length from 1000 to 1600 bp. Two structural domains with specific evolutionary patterns were identified. These were (1) conserved sequence blocks containing primary sequence motifs, including dinucleotide pyrimidine-purine series and long T-stretches, located at the 5' end adjacent to the tRNA(Ile) gene and (2) a hypervariable domain at the 3' end characterized by increased nucleotide divergence and size variation. A high frequency of A<-->T transversions at nucleotide substitution level indicated directional mutation pressure. The phylogenetic usefulness of the insect control region is discussed.

Methods for the recovery of mitochondrial DNA sequences from museum specimens of myiasis-causing flies.

Mitochondrial DNA (mtDNA) sequences from eight species of myiasis-causing flies, stored for up to 50 years, were amplified successfully. Universal primers were used to amplify six specific regions from total genomic DNA, including five mtDNA genes. The comparison of phenol/chloroform, DNAzol and Chelex techniques for DNA extraction showed that the DNAzol reagent was the most efficient in retrieving DNA from museum specimens, although the Chelex extraction procedure is currently the most frequently reported method. Comparison of the universal primer sequences with the homologous sequences of Cochliomyia hominivorax Coquerel and Chrysomya putoria Wiedemann (Diptera: Calliphoridae) revealed mismatches that could contribute to the low recovery of a short sequence from subunit II of cytochrome oxidase. The ability to characterize mtDNA markers from museum specimens should be useful in comparative studies of contemporary samples and should help in elucidating species introduction, colonization and dispersal.

Intraspecific mitochondrial DNA variation in the Colorado potato beetle, Leptinotarsa decemlineata (Coleoptera: Chrysomelidae).

Restriction endonuclease analyses of mitochondrial DNA (mtDNA) were used to examine genetic variability and population structure in Leptinotarsa decemlineata (Say). A group of three enzymes, EcoRI, HpaI, and PstI, was used to reveal polymorphism both within and among some of the 10 populations tested, yielding 16 haplotypes in combination. The frequencies of these 16 haplotypes differed significantly across geographic regions, indicating some partitioning of mtDNA haplotypes. Estimates of mtDNA sequence divergence (delta) between haplotypes ranged from 0.016 to 0.135%, suggesting local differentiation of mtDNA in some populations. Analysis of these data suggests that Texas was colonized by more than one mtDNA lineage, most likely originating in Mexico. We hypothesize that a larger founder size for the initial introductions or high levels of variability in the parent population at the edge of the CPB expanding range led to the initial partitioning of haplotypes observed in samples from Texas.