Ana Carolina M. Junqueira

Structure and evolution of the mitochondrial genomes of Haematobia irritans and Stomoxys calcitrans : The Muscidae (Diptera: Calyptratae) perspective

We present the first two mitochondrial genomes of Muscidae dipterans for the species Haematobia irritans (the horn fly) and Stomoxys calcitrans (the stable fly). Typical insect mtDNA features are described, such as a high A+T content (79.1% and 78.9%, respectively), the preference for A+T-rich codons, and the evidence of a non-optimal codon usage. The strong A+T enrichment partially masks another nucleotide content bias maintained by A+C mutation pressure in these Muscidae mtDNAs. The analysis of this data provides a model of metazoans tRNA anticodon evolution, based on the selection hypothesis of anticodon versatility. H. irritans mitochondrial genome (16078 bp) is structurally similar to the hypothetical ancestral mitochondrial genome of arthropods and its control region (A+ T-rich region in insects) organization is consistent with the structure described for Brachycera dipterans. On the other hand, the mitochondrial genome of S. calcitrans is approximately 2kb longer (18 kb), characterized by the presence of approximately 550 bp tandem repeats in the control region, and an extra copy of trnI remarkably similar to a duplicated element of blowflies mtDNA. Putative sequence elements, involved in the regulation of transcription and replication of the mtDNA, were reliably identified in S. calcitrans control region despite the 0.8-1.5 kb gap uncovered from this genome. The use of amino acid and nucleotide sequences of concatenated mitochondrial protein-coding genes (PCGs) in phylogenetic reconstructions of Diptera does not support the monophyly of Muscomorpha, as well as the monophyly of Acalyptratae. Within the Calyptratae group, the inclusion of Muscidae (Muscoidea) as a sister group of Calliphoridae (Oestroidea) implies in a potential conflict concerning the monophyly of the superfamily Oestroidea.

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.

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

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

Molecular phylogeny of the Drosophila tripunctata and closely related species groups (Diptera: Drosophilidae)

We suggest a new phylogenetic hypothesis for the tripunctata radiation based on sequences of mitochondrial genes. Phylogenetic trees were reconstructed by parsimony, maximum likelihood and Bayesian methods. We performed tests for hypotheses of monophyly for taxonomic groups and other specific hypotheses. Results reject the monophyly for the tripunctata group whereas monophyly is not rejected for the tripunctata radiation and other specific groups within the radiation. Although most of the basal nodes were unresolved we were able to identify four clusters within the tripunctata radiation. These results suggest the collection of additional data before a proper taxonomic revision could be proposed.

The Mitochondrial Control Region of Blowflies (Diptera: Calliphoridae): A Hot Spot for Mitochondrial Genome Rearrangements

Abstract The family Calliphoridae consists of myiasis-causing flies, including species of economic, forensic, and medical importance. In this study, the complete control regions (CRs) of mitochondrial DNA from 15 calliphorid species were sequenced and structurally characterized. The CRs had a high content of adenines (A) and thymines (T) and varied in length from 854 to 2,018 bp, showing intraspecific variations in sequence and length. Two major domains were identified: the conserved domain containing conserved sequence blocks and cis-regulatory structures that may be related to the transcription and the origin of replication of mitochondrial DNA, and the variable domain, containing high sequence and length variation. Within the variable domain, duplication of the tRNAIle gene, previously reported for three Chrysomya species, was identified in two more species of this genus and in two species of two other genera. The structural characterization shows the plasticity of the mitochondrial genome in dipterans. The organizational similarities of the duplicated region found in different species and the possible origin of the duplicated genes are discussed.

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.

Large-scale mitogenomics enables insights into Schizophora (Diptera) radiation and population diversity.

True flies are insects of the order Diptera and encompass one of the most diverse groups of animals on Earth. Within dipterans, Schizophora represents a recent radiation of insects that was used as a model to develop a pipeline for generating complete mitogenomes using various sequencing platforms and strategies. 91 mitogenomes from 32 different species were sequenced and assembled with high fidelity, using amplicon, whole genome shotgun or single molecule sequencing approaches. Based on the novel mitogenomes, we estimate the origin of Schizophora within the Cretaceous-Paleogene (K-Pg) boundary, about 68.3 Ma. Detailed analyses of the blowfly family (Calliphoridae) place its origin at 22 Ma, concomitant with the radiation of grazing mammals. The emergence of ectoparasitism within calliphorids was dated 6.95 Ma for the screwworm fly and 2.3 Ma for the Australian sheep blowfly. Varying population histories were observed for the blowfly Chrysomya megacephala and the housefly Musca domestica samples in our dataset. Whereas blowflies (n = 50) appear to have undergone selective sweeps and/or severe bottlenecks in the New World, houseflies (n = 14) display variation among populations from different zoogeographical zones and low levels of gene flow. The reported high-throughput mitogenomics approach for insects enables new insights into schizophoran diversity and population history of flies.

The microbiomes of blowflies and houseflies as bacterial transmission reservoirs

Blowflies and houseflies are mechanical vectors inhabiting synanthropic environments around the world. They feed and breed in fecal and decaying organic matter, but the microbiome they harbour and transport is largely uncharacterized. We sampled 116 individual houseflies and blowflies from varying habitats on three continents and subjected them to high-coverage, whole-genome shotgun sequencing. This allowed for genomic and metagenomic analyses of the host-associated microbiome at the species level. Both fly host species segregate based on principal coordinate analysis of their microbial communities, but they also show an overlapping core microbiome. Legs and wings displayed the largest microbial diversity and were shown to be an important route for microbial dispersion. The environmental sequencing approach presented here detected a stochastic distribution of human pathogens, such as Helicobacter pylori, thereby demonstrating the potential of flies as proxies for environmental and public health surveillance.

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.

Complete genome sequence of bacillus altitudinis type strain SGAir0031 isolated from tropical air collected in Singapore

ABSTRACT Bacillus altitudinis strain SGAir0031 (Firmicutes) was isolated from tropical air samples collected in Singapore. Its genome was assembled using short reads and single-molecule real-time sequencing, comprising one chromosome with 3.81 Mb and one plasmid with 32 kb. The genome consists of 3,820 protein-coding genes, 81 tRNAs, and 24 rRNAs.