Species-specific mitochondrial gene rearrangements in biting midges and vector species identification

Comparative mitogenomic analysis of subterranean and surface amphipods (Crustacea, Amphipoda) with special reference to the family Crangonyctidae

Mitochondrial genomes play important roles in studying genome evolution, phylogenetic analyses, and species identification. Amphipods (Class Malacostraca, Order Amphipoda) are one of the most ecologically diverse crustacean groups occurring in a diverse array of aquatic and terrestrial environments globally, from freshwater streams and lakes to groundwater aquifers and the deep sea, but we have a limited understanding of how habitat influences the molecular evolution of mitochondrial energy metabolism. Subterranean amphipods likely experience different evolutionary pressures on energy management compared to surface-dwelling taxa that generally encounter higher levels of predation and energy resources and live in more variable environments. In this study, we compared the mitogenomes, including the 13 protein-coding genes involved in the oxidative phosphorylation (OXPHOS) pathway, of surface and subterranean amphipods to uncover potentially different molecular signals of energy metabolism between surface and subterranean environments in this diverse crustacean group. We compared base composition, codon usage, gene order rearrangement, conducted comparative mitogenomic and phylogenomic analyses, and examined evolutionary signals of 35 amphipod mitogenomes representing 13 families, with an emphasis on Crangonyctidae. Mitogenome size, AT content, GC-skew, gene order, uncommon start codons, location of putative control region (CR), length of rrnL and intergenic spacers differed between surface and subterranean amphipods. Among crangonyctid amphipods, the spring-dwelling Crangonyx forbesi exhibited a unique gene order, a long nad5 locus, longer rrnL and rrnS loci, and unconventional start codons. Evidence of directional selection was detected in several protein-encoding genes of the OXPHOS pathway in the mitogenomes of surface amphipods, while a signal of purifying selection was more prominent in subterranean species, which is consistent with the hypothesis that the mitogenome of surface-adapted species has evolved in response to a more energy demanding environment compared to subterranean amphipods. Overall, gene order, locations of non-coding regions, and base-substitution rates points to habitat as an important factor influencing the evolution of amphipod mitogenomes.

Unusual rearrangements of mitogenomes in Diptera revealed by comparative analysis of 135 tachinid species (Insecta, Diptera, Tachinidae)

The Tachinidae is one of the most speciose families in Diptera, and the exclusively parasitoid species play an important role in regulating populations of many herbivorous insects in ecosystems, including many agricultural pests. To better comprehend the characteristics and evolution of the mitochondrial genome for the Tachinidae, we are adding a massive amount of new molecular data by assembling the mitogenomes for 71 genera and 135 species from all four tachinid subfamilies through next-generation sequencing, and we are presenting the most comprehensive mitogenomic phylogenetic analysis of this family so far. Extensive rearrangements observed in the mitogenome of Admontia podomyia (Exoristinae) are unique for the entire suborder Cyclorrhapha. The rearrangement pattern suggests that the process involved a tandem duplication of the complete mitogenome, followed by both random and nonrandom loss of one copy of each gene. Additionally, five minor mitogenome rearrangements are discovered and described in three subfamilies. We present the largest species-level phylogenetic hypothesis for Tachinidae to date, based on mitogenomes of 152 species of Tachinidae, representing all four subfamilies and with five non-tachinid outgroups. Our analyses support the monophyly of the Tachinidae and most tribes and genera were recovered with good support, but the higher-level phylogenetic relationships within Tachinidae were poorly resolved, indicating that mitogenome data alone are not enough to unambiguously resolve the deeper phylogenetic relationships within Tachinidae.

Comparative mitogenome analyses of twelve non-biting flies and provide insights into the phylogeny of Chironomidae (Diptera: Culicomorpha)

The family Chironomidae is represented by seven subfamilies in China, among which Chironominae and Orthocladiinae are the most diverse. To gain a better understanding of the architecture and evolution of the mitogenomes of Chironomidae, we sequenced mitogenomes of twelve species (including two published species) of the two subfamilies Chironominae and Orthocladiinae, and comparative mitogenomic analyses were performed. Thus, we identified highly conserved genome organization of twelve species with regard to genome content, nucleotide and amino acid composition, codon usage, and gene characteristics. The K_a/K_s values of most protein-coding genes were far smaller than 1, indicating that these genes were evolving under purifying selection. Phylogenetic relationships between the family Chironomidae were reconstructed using 23 species representing six subfamilies, based on protein-coding genes and rRNAs using Bayesian Inference and Maximum Likelihood. Our results suggested the following relationship within the Chironomidae: (Podonominae + Tanypodinae) + (Diamesinae + (Prodiamesinae + (Orthocladiinae + Chironominae))). This study contributes to the mitogenomic database of Chironomidae, which will be significant for studing the mitogenome evolution of Chironomidae.

DNA Barcoding of Culicoides Latreille (Diptera: Ceratopogonidae) From Thailand Reveals Taxonomic Inconsistencies and Novel Diversity Among Reported Sequences

Recent focus on Culicoides species diversity in Thailand was prompted by a need to identify vectors responsible for the transmission of African Horse Sickness in that country. To assist rapid genetic identification of species, we sampled mitochondrial cytochrome c oxidase subunit I (COI) DNA barcodes (N = 78) from 40 species of Culicoides biting midge from Thailand, including 17 species for which DNA barcodes were previously unavailable. The DNA barcodes were assigned to 39 Barcode Identification Numbers (BINs) representing terminal genetic clusters at the Barcode of Life Data systems (BOLD). BINs assisted with comparisons to published conspecific DNA barcodes and allowed partial barcodes obtained from seven specimens to be associated with BINs by their similarity. Some taxonomic issues were revealed and attributed to the possible misidentification of earlier reported specimens as well as a potential synonymy of C. elbeli Wirth & Hubert and C. menglaensis Chu & Liu. Comparison with published BINs also revealed genetic evidence of divergent population processes and or potentially cryptic species in 16 described taxa, flagged by their high levels of COI sequence difference among conspecifics. We recommend the BOLD BIN system to researchers preparing DNA barcodes of vouchered species for public release. This will alert them to taxonomic incongruencies between their records and publicly released DNA barcodes, and also flag genetically deep and potentially novel diversity in described species.

Characterization of the Complete Mitochondrial Genome of Eight Diurnal Hawkmoths (Lepidoptera: Sphingidae): New Insights into the Origin and Evolution of Diurnalism in Sphingids

In this study, the mitochondrial genomes of 22 species from three subfamilies in the Sphingidae were sequenced, assembled, and annotated. Eight diurnal hawkmoths were included, of which six were newly sequenced (Hemaris radians, Macroglossum bombylans, M. fritzei, M. pyrrhosticta, Neogurelca himachala, and Sataspes xylocoparis) and two were previously published (Cephonodes hylas and Macroglossum stellatarum). The mitochondrial genomes of these eight diurnal hawkmoths were comparatively analyzed in terms of sequence length, nucleotide composition, relative synonymous codon usage, non-synonymous/synonymous substitution ratio, gene spacing, and repeat sequences. The mitogenomes of the eight species, ranging in length from 15,201 to 15,461 bp, encode the complete set of 37 genes usually found in animal mitogenomes. The base composition of the mitochondrial genomes showed A+T bias. The most commonly used codons were UUA (Leu), AUU (Ile), UUU (Phe), AUA (Met), and AAU (Asn), whereas GCG (Ala) and CCG (Pro) were rarely used. A phylogenetic tree of Sphingidae was constructed based on both maximum likelihood and Bayesian methods. We verified the monophyly of the four current subfamilies of Sphingidae, all of which had high support. In addition, we performed divergence time estimation and ancestral character reconstruction analyses. Diurnal behavior in hawkmoths originated 29.19 million years ago (Mya). It may have been influenced by the combination of herbaceous flourishing, which occurred 26-28 Mya, the uplift of the Tibetan Plateau, and the large-scale evolution of bats in the Oligocene to Pre-Miocene. Moreover, diurnalism in hawkmoths had multiple independent origins in Sphingidae.

Mitochondrial genome sequencing, mapping, and assembly benchmarking for Culicoides species (Diptera: Ceratopogonidae)

BACKGROUND

Mitochondrial genomes are the most sequenced genomes after bacterial and fungal genomic DNA. However, little information on mitogenomes is available for multiple metazoan taxa, such as Culicoides, a globally distributed, megadiverse genus containing 1,347 species.

AIM

 Generating novel mitogenomic information from single Culicoides sonorensis and C. biguttatus specimens, comparing available mitogenome mapping and de novo assembly tools, and identifying the best performing strategy and tools for Culicoides species.

RESULTS

We present two novel and fully annotated mitochondrial haplotypes for two Culicoides species, C. sonorensis and C. biguttatus. We also annotated or re-annotated the only available reference mitogenome for C. sonorensis and C. arakawae. All species present a high similarity in mitogenome organization. The general gene arrangement for all Culicoides species was identical to the ancestral insect mitochondrial genome. Only short spacers were found in C. sonorensis (up to 30 bp), contrary to C. biguttatus (up to 114 bp). The mitochondrial genes ATP8, NAD2, NAD6, and LSU rRNA exhibited the highest nucleotide diversity and pairwise interspecific p genetic distance, suggesting that these genes might be suitable and complementary molecular barcodes for Culicoides identification in addition to the commonly utilized COI gene. We observed performance differences between the compared mitogenome generation strategies. The mapping strategy outperformed the de novo assembly strategy, but mapping results were partially biased in the absence of species-specific reference mitogenome. Among the utilized tools, BWA performed best for C. sonorensis while SPAdes, MEGAHIT, and MitoZ were among the best for C. biguttatus. The best-performing mitogenome annotator was MITOS2. Additionally, we were able to recover exogenous mitochondrial DNA from Bos taurus (biting midges host) from a C. biguttatus blood meal sample.

CONCLUSIONS

Two novel annotated mitogenome haplotypes for C. sonorensis and C. biguttatus using High-Throughput Sequencing are presented. Current results are useful as the baseline for mitogenome reconstruction of the remaining Culicoides species from single specimens to HTS and genome annotation. Mapping to a species-specific reference mitogenome generated better results for Culicoides mitochondrial genome reconstruction than de novo assembly, while de novo assembly resulted better in the absence of a closely related reference mitogenome. These results have direct implications for molecular-based identification of these vectors of human and zoonotic diseases, setting the basis for using the whole mitochondrial genome as a marker in Culicoides identification.

First complete mitogenomes of Diamesinae, Orthocladiinae, Prodiamesinae, Tanypodinae (Diptera: Chironomidae) and their implication in phylogenetics

BACKGROUND

The mitochondrial genome (mitogenome) has been extensively used for phylogenetic and evolutionary analysis in Diptera, but the study of mitogenome is still scarce in the family Chironomidae.

METHODS

Here, the first complete mitochondrial genomes of four Chironomid species representing Diamesinae, Orthocladiinae, Prodiamesinae and Tanypodinae are presented. Coupled with published mitogenomes of two, a comparative mitochondrial genomic analysis between six subfamilies of Chironomidae was carried out.

RESULTS

Mitogenomes of Chironomidae are conserved in structure, each contains 37 typical genes and a control region, and all genes arrange the same gene order as the ancestral insect mitogenome. Nucleotide composition is highly biased, the control region displayed the highest A + T content. All protein coding genes are under purifying selection, and the ATP8 evolves at the fastest rate. In addition, the phylogenetic analysis covering six subfamilies within Chironomidae was conducted. The monophyly of Chironomidae is strongly supported. However, the topology of six subfamilies based on mitogenomes in this study is inconsistent with previous morphological and molecular studies. This may be due to the high mutation rate of the mitochondrial genetic markers within Chironomidae. Our results indicate that mitogenomes showed poor signals in phylogenetic reconstructions at the subfamily level of Chironomidae.

Comparative Mitogenomics in Hyalella (Amphipoda: Crustacea)

We present the sequencing and comparative analysis of 17 mitochondrial genomes of Nearctic and Neotropical amphipods of the genus Hyalella, most from the Andean Altiplano. The mitogenomes obtained comprised the usual 37 gene-set of the metazoan mitochondrial genome showing a gene rearrangement (a reverse transposition and a reversal) between the North and South American Hyalella mitogenomes. Hyalella mitochondrial genomes show the typical AT-richness and strong nucleotide bias among codon sites and strands of pancrustaceans. Protein-coding sequences are biased towards AT-rich codons, with a preference for leucine and serine amino acids. Numerous base changes (539) were found in tRNA stems, with 103 classified as fully compensatory, 253 hemi-compensatory and the remaining base mismatches and indels. Most compensatory Watson–Crick switches were AU -> GC linked in the same haplotype, whereas most hemi-compensatory changes resulted in wobble GU and a few AC pairs. These results suggest a pairing fitness increase in tRNAs after crossing low fitness valleys. Branch-site level models detected positive selection for several amino acid positions in up to eight mitochondrial genes, with atp6 and nad5 as the genes displaying more sites under selection.

Rapid identification of the invasive fall armyworm Spodoptera frugiperda (Lepidoptera, Noctuidae) using species-specific primers in multiplex PCR

The fall armyworm (FAW), Spodoptera frugiperda (Smith), is a major pest native to the Americas. A recent invasion of FAWs from Africa eastward to South Asia, the Indochina Peninsula, and mainland China has received much attention due to the considerable economic losses in agriculture. FAWs can rapidly colonise a new area, likely due to the wide range of host plants, good flying capability, and high egg production. Therefore, a convenient, quick, and accurate tool for FAW identification is urgently required to establish a FAW invasion management strategy. In this study, FAW-specific primers were designed to recognise FAWs on the basis of internal transcribed spacer 1 (ITS1). The results revealed the accurate FAW recognition of the three congeneric species and eight common corn lepidopteran pests, especially at their larval stage. Furthermore, species-specific primers have confirmed their efficacy by using 69 FAW specimens from Taiwan, Thailand, and the United States, with a 96% success rate, excluding 3 decayed specimens. By using the simple, reliable, and convenient FAW-specific primers, a pest management programme can be developed not only to reduce sequencing costs and experimental time from 2 days to 4 h, but eradicate the FAW as soon as it enters a new area.

A phylogenetic analysis of the biting midges belonging to Culicoides Latreille (Diptera: Ceratopogonidae) subgenus Avaritia using molecular data

Background

Within the genus Culicoides (Diptera: Ceratopogonidae), the subgenus Avaritia is of particular interest as it contains a significant number of economically important vector species. Disagreements about the systematic classification of species within this subgenus have resulted in a taxonomic imbroglio.

Methods

A molecular phylogeny of the subgenus Avaritia was conducted to test the existing systematic classification, which is based on phenetic assessment of morphological characters. Three nuclear ribosomal markers, internal transcribed spacer 1 and 2 (ITS1, ITS2), 5.8S, and three mitochondrial markers, cytochrome c oxidase subunit 1 and 2, and cytochrome b (cox1, cox2 and cytb), were obtained for 37 species of the subgenus Avaritia from all six biogeographical regions. Phylogenetic reconstructions using these genes independently and in combination were implemented using Bayesian inference analysis and maximum likelihood methods.

Results

Phylogenetic reconstructions gave strong support to several monophyletic groups within the subgenus Avaritia. Both C. actoni and C. pusillus formed a single clade with C. grahamii so their respective groups, the Actoni and Pusillus groups, have been merged with the Grahamii group. Some support was provided for the Boophagus and Jacobsoni groups. A group of species currently placed into the Orientalis group clustered in a clade with poor support. The Obsoletus group was defined as a sister clade to all other Avaritia groups. The clade including the Imicola group was well supported based on phylogenetic criteria.

Conclusions

This phylogenetic study combining five distinct molecular markers has provided meaningful insights into the systematic relationships of Culicoides (Avaritia) and highlighted future directions to continue the study of this subgenus. While the cox2 marker appeared to be useful to investigate closely related species, the 5.8S marker was highly conserved and uninformative. Further investigations including species absent from this work are needed to confirm the proposed systematic scheme. However, this systematic scheme can now serve as a foundation to investigate cryptic species affiliation within the subgenus. We advocate that future studies employ a combination of morphological and molecular analyses.

Surveillance of Culicoides biting midges in northern Honshu, Japan, during the period of Akabane virus spread

A surveillance of Culicoides biting midges with light suction traps was conducted in the northern region of Honshu, main island of Japan, during the summers and autumns of 2009 and 2010. A total of 106 trap collections across 37 cattle farms were investigated for the structure and distribution of Culicoides species. Forty-thousand and one hundred forty-nine specimens of Culicoides biting midges were identified at the species level, and ≥19 species were included in the specimens. Culicoides oxystoma, which is a known major vector of Akabane virus (AKAV), appeared not to have expanded in northern Honshu during the surveillance. Of the potential AKAV vectors suggested by a previous laboratory experiment, C. tainanus and C. punctatus widely infested cowsheds across northern Honshu. The AKAV circulation was confirmed by serological surveillance of sentinel cattle in northern Honshu during the summer and autumn of 2010 and, consequently, >200 calves affected by the virus were identified as of spring 2011. Our surveillance demonstrated that C. tainanus and C. punctatus were widely spread and often dominated at cattle farms in/around the seroconverted regions, and our results thus suggest that these species played a critical role in the AKAV transmission in 2010. Because the distribution ranges of C. tainanus and C. punctatus cover almost all of mainland Japan, a potential risk of AKAV transmission might be expected even in areas outside the range of C. oxystoma.

Complete mitochondrial genomes of eight seahorses and pipefishes (Syngnathiformes: Syngnathidae): insight into the adaptive radiation of syngnathid fishes

Background

The evolution of male pregnancy is the most distinctive characteristic of syngnathids, and their specialized life history traits make syngnathid species excellent model species for many issues in biological evolution. However, the origin of syngnathids and the evolutionary divergence time of different syngnathid species remain poorly resolved. Comprehensive phylogenetic studies of the Syngnathidae will provide critical evidence to elucidate their origin, evolution, and dispersal patterns.

Results

We sequenced the mitochondrial genomes of eight syngnathid species in this study, and the estimated divergence times suggested that syngnathids diverged from other teleosts approximately 48.8 Mya during the Eocene period. Selection analysis showed that many mitochondrial genes of syngnathids exhibited significantly lower Ka/Ks values than those of other teleosts. The two most frequently used codons in syngnathid fishes were different from those in other teleosts, and a greater proportion of the mitochondrial simple sequence repeats (SSRs) were distributed in non-coding sequences in syngnathids compared with other teleosts.

Conclusions

Our study indicated that syngnathid fishes experienced an adaptive radiation process during the early explosion of species. Syngnathid mitochondrial OXPHOS genes appear to exhibit depressed Ka/Ks ratios compared with those of other teleosts, and this may suggest that their mitogenomes have experienced strong selective constraints to eliminate deleterious mutations.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1430-3) contains supplementary material, which is available to authorized users.

Oral Susceptibility of Japanese Culicoides (Diptera: Ceratopogonidae) Species to Akabane Virus

To test their virus susceptibility and capacity as transmission vectors, Japanese Culicoides species were artificially fed a mixture of Akabane virus (AKAV) and bovine blood, and virus recovery was attempted from infected midges 9-11 d post-exposure. Culicoides tainanus, C. punctatus, C. humeralis, C. jacobsoni, C. oxystoma, and C. asiana were found to be orally susceptible to AKAV. Virus titers in single infected midges of C. tainanus, C. oxystoma, C. punctatus, and C. jacobsoni ranged from 100.75 to 104.0 TCID50 (tissue culture infectious dose). The titers in the infected C. oxystoma were significantly higher than those in the other infected species. Viral RNA was detected from both midges testing positive and those testing negative for infectious virus particles, but the viral RNA copies in the infectious virus-negative midges were significantly lower than those in the infectious virus-positive midges. Lower viral amplification, limited dissemination or both caused by tissue barriers might occur in infected midges from which infectious viruses were undetectable. A fully disseminated infection was developed in orally infected C. oxystoma and C. tainanus. This finding indicates their capacity to transmit AKAV, assuming that salivary gland barriers have limited effects on viral entry to and replication in salivary gland tissue. This result also suggests that the other orally susceptible species are potentially competent for AKAV transmission and would be considered active vectors of its spread.

Molecular differentiation and species composition of genus Culicoides biting midges (Diptera: Ceratopogonidae) in different habitats in southern China

Culicoides biting midges (Diptera:Ceratopogonidae) cause a significant biting nuisance to humans, livestock, which are the biological vectors of a range of risky pathogens. Accurate illustration of vector play a key role in arthropod borne diseases surveillance. However, few studies have focused on the Culicoides, which caused bluetongue disease in 29 provinces of China since 1979. In this study, we assessed cytochrome oxidase subunit I (COI) mtDNA molecular marker for identification of ten major vector species and analyzed the Culicoides species community and diversity in different habitats. A total of 20,795 Culicoides samples collected from 11 sample sites were identified as 23 species belonging to 7 subgenera. Sequences of COI gene worked well as barcodes for identifying all the determined specimen in this study and were comparable with the existing sequence data from GenBank. We first reported COI barcode sequences of C. morisitai, C. insignipennis and C. homotomus. Morphological identification of Culicoides spp. samples within southern China appears relatively robust and some unidentified species were required further study. Our study shows that the COI sequence data can be used as a tool to identify species of Culicoides in Jiangxi Province. In our sampled area, the most abundant species was C. arakawae (61.89%), followed by C. oxystoma (13.77%), C. punctatus (10.10%), C. nipponensis (8.82%), C. homotomus (3.19%) and C.morisitai (1.17%) in this study. C. punctatus was the dominant species of Park habitat (62.22%), C. arakawae was the predominant species of Chicken habitat (96.66%), Vegetable plot habitat (92.0%), and Peasant household habitat (83.21%), respectively. C. oxystoma was the abundant species of Residential area habitat (40.11%), Hospital habitat (56.65%), and Pig & Cow habitat (48.77%), respectively. Results also show that the potential Bluetongue virus vectors belong to the Obsoletus and Pulicaris groups, are also included and relatively abundant, notably: Culicoides punctatus. These findings expand the current knowledge of Culicoides population composition in the southern part of China.

Molecular phylogeny of 42 species of Culicoides (Diptera, Ceratopogonidae) from three continents

The genus Culicoides includes vectors of important animal diseases such as bluetongue and Schmallenberg virus (BTV and SBV). This genus includes 1300 species classified in 32 subgenera and 38 unclassified species. However, the phylogenetic relationships between different subgenera of Culicoides have never been studied. Phylogenetic analyses of 42 species belonging to 12 subgenera and 8 ungrouped species of genus Culicoides from Ecuador, France, Gabon, Madagascar and Tunisia were carried out using two molecular markers (28S rDNA D1 and D2 domains and COI mtDNA). Sequences were subjected to non-probabilistic (maximum parsimony) and probabilistic (Bayesian inference (BI)) approaches. The subgenera Monoculicoides, Culicoides, Haematomyidium, Hoffmania, Remmia and Avaritia (including the main vectors of bluetongue disease) were monophyletic, whereas the subgenus Oecacta was paraphyletic. Our study validates the subgenus Remmia (= Schultzei group) as a valid subgenus, outside of the subgenus Oecacta. In Europe, Culicoides obsoletus, Culicoides scoticus and Culicoides chiopterus should be part of the Obsoletus complex whereas Culicoides dewulfi should be excluded from this complex. Our study suggests that the current Culicoides classification needs to be revisited with modern tools.

Evolution of mitochondrial genomes in Baikalian amphipods

BACKGROUND

Amphipods (Crustacea) of Lake Baikal are a very numerous and diverse group of invertebrates generally believed to have originated by adaptive radiation. The evolutionary history and phylogenetic relationships in Baikalian amphipods still remain poorly understood. Sequencing of mitochondrial genomes is a relatively feasible way for obtaining a set of gene sequences suitable for robust phylogenetic inferences. The architecture of mitochondrial genomes also may provide additional information on the mechanisms of evolution of amphipods in Lake Baikal.

RESULTS

Three complete and four nearly complete mitochondrial genomes of Baikalian amphipods were obtained by high-throughput sequencing using the Illumina platform. A phylogenetic inference based on the nucleotide sequences of all mitochondrial protein coding genes revealed the Baikalian species to be a monophyletic group relative to the nearest non-Baikalian species with a completely sequenced mitochondrial genome - Gammarus duebeni. The phylogeny of Baikalian amphipods also suggests that the shallow-water species Eulimnogammarus has likely evolved from a deep-water ancestor, however many other species have to be added to the analysis to test this hypothesis. The gene order in all mitochondrial genomes of studied Baikalian amphipods differs from the pancrustacean ground pattern. Mitochondrial genomes of four species possess 23 tRNA genes, and in three genomes the extra tRNA gene copies have likely undergone remolding. Widely varying lengths of putative control regions and other intergenic spacers are typical for the mitochondrial genomes of Baikalian amphipods.

CONCLUSIONS

The mitochondrial genomes of Baikalian amphipods display varying organization suggesting an intense rearrangement process during their evolution. Comparison of complete mitochondrial genomes is a potent approach for studying the amphipod evolution in Lake Baikal.

DNA barcoding and surveillance sampling strategies for Culicoides biting midges (Diptera: Ceratopogonidae) in southern India

Background

Culicoides spp. biting midges transmit bluetongue virus (BTV), the aetiological agent of bluetongue (BT), an economically important disease of ruminants. In southern India, hyperendemic outbreaks of BT exert high cost to subsistence farmers in the region, impacting on sheep production. Effective Culicoides spp. monitoring methods coupled with accurate species identification can accelerate responses for minimising BT outbreaks. Here, we assessed the utility of sampling methods and DNA barcoding for detection and identification of Culicoides spp. in southern India, in order to provide an informed basis for future monitoring of their populations in the region.

Methods

Culicoides spp. collected from Tamil Nadu and Karnataka were used to construct a framework for future morphological identification in surveillance, based on sequence comparison of the DNA barcode region of the mitochondrial cytochrome c oxidase I (COI) gene and achieving quality standards defined by the Barcode of Life initiative. Pairwise catches of Culicoides spp. were compared in diversity and abundance between green (570 nm) and ultraviolet (UV) (390 nm) light emitting diode (LED) suction traps at a single site in Chennai, Tamil Nadu over 20 nights of sampling in November 2013.

Results

DNA barcode sequences of Culicoides spp. were mostly congruent both with existing DNA barcode data from other countries and with morphological identification of major vector species. However, sequence differences symptomatic of cryptic species diversity were present in some groups which require further investigation. While the diversity of species collected by the UV LED Center for Disease Control (CDC) trap did not significantly vary from that collected by the green LED CDC trap, the UV CDC significantly outperformed the green LED CDC trap with regard to the number of Culicoides individuals collected.

Conclusions

Morphological identification of the majority of potential vector species of Culicoides spp. samples within southern India appears relatively robust; however, potential cryptic species diversity was present in some groups requiring further investigation. The UV LED CDC trap is recommended for surveillance of Culicoides in southern India.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1722-z) contains supplementary material, which is available to authorized users.

Comparative Mt Genomics of the Tipuloidea (Diptera: Nematocera: Tipulomorpha) and Its Implications for the Phylogeny of the Tipulomorpha

A traditionally controversial taxon, the Tipulomorpha has been frequently discussed with respect to both its familial composition and relationships with other Nematocera. The interpretation of internal relationships within the Tipuloidea, which include the Tipulidae sensu stricto, Cylindrotomidae, Pediciidae and Limoniidae, is also problematic. We sequenced the first complete mitochondrial (mt) genome of Symplecta hybrida (Meigen, 1804), which belongs to the subfamily Chioneinae of family Limoniidae, and another five nearly complete mt genomes from the Tipuloidea. We did a comparative analysis of these mt genomics and used them, along with some other representatives of the Nematocera to construct phylogenetic trees. Trees inferred by Bayesian methods strongly support a sister-group relationship between Trichoceridae and Tipuloidea. Tipulomorpha are not supported as the earliest branch of the Diptera. Furthermore, phylogenetic trees indicate that the family Limoniidae is a paraphyletic group.

The complete mitochondrial genome of the Atylotus miser (Diptera: Tabanomorpha: Tabanidae), with mitochondrial genome phylogeny of lower Brachycera (Orthorrhapha)

Brachycera is a clade with over 80,000 described species and originated from the Mesozoic, and its larvae employ comprehensive feeding strategies. The phylogeny of the lower Brachycera has been studied intensively over the past decades. In order to supplement the lack of genetic data in this important group, we sequenced the complete mitochondrial (mt) genome of Atylotus miser as well as the nearly complete mt genomes of another 11 orthorrhaphous flies. The mt genome of A. miser is 15,858bp, which is typical of Diptera, with 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and a 993bp control region. The rest of the orthorrhaphous mt genomes in our study have the similar structure with A. miser. Additionally, we conducted a phylogenetic analysis of 20 mt genomes using Maximum-likelihood and Bayesian methods in order to reconstruct the evolutionary relationship of Orthorrhapha. The results show that all infraorders of Brachycera are monophyletic, and a relationship of Tabanomorpha+((Xylophagomorpha+Stratiomyomorpha)+Muscomorpha) has been proposed. Within Xylophagomorpha, Nemestrinoidae forms the sister group of Xylophagidae.

Population Genetic Structure and Potential Incursion Pathways of the Bluetongue Virus Vector Culicoides brevitarsis (Diptera: Ceratopogonidae) in Australia

Culicoides brevitarsis is a vector of the bluetongue virus (BTV), which infects sheep and cattle. It is an invasive species in Australia with an assumed Asian/South East Asian origin. Using one mitochondrial marker (i.e., part of the cytochrome oxidase subunit I gene) and six nuclear markers, we inferred population genetic structure and possible incursion pathways for Australian C. brevitarsis. Nine mitochondrial haplotypes, with low nucleotide sequence diversity (0.0–0.7%) among these, were identified in a sample of 70 individuals from seven sites. Both sets of markers revealed a homogeneous population structure, albeit with evidence of isolation by distance and two genetically distinct clusters distributed along a north-to-south cline. No evidence of a cryptic species complex was found. The geographical distribution of the mitochondrial haplotypes is consistent with at least two incursion pathways into Australia since the arrival of suitable livestock hosts. By contrast, 15 mitochondrial haplotypes, with up to four times greater nucleotide sequence diversity (0.0–2.9%) among these, were identified in a sample of 16 individuals of the endemic C. marksi (sampled from a site in South Australia and another in New South Wales). A phylogenetic tree inferred using the mitochondrial marker revealed that the Australian and Japanese samples of C. brevitarsis are as evolutionarily different from one another as some of the other Australian species (e.g., C. marksi, C. henryi, C. pallidothorax) are. The phylogenetic tree placed four of the species endemic to Australia (C. pallidothorax, C. bundyensis, C. marksi, C. henryi) in a clade, with a fifth such species (C. bunrooensis) sharing a common ancestor with that clade and a clade comprising two Japanese species (C. verbosus, C. kibunensis).

Mitochondrial genomes of two phlebotomine sand flies, Phlebotomus chinensis and Phlebotomus papatasi (Diptera: Nematocera), the first representatives from the family Psychodidae

BACKGROUND

Leishmaniasis is a worldwide but neglected disease of humans and animal transmitted by sand flies, vectors that also transmit other important diseases. Mitochondrial genomes contain abundant information for population genetic and phylogenetic studies, important in disease management. However, the available mitochondrial sequences of these crucial vectors are limited, emphasizing the need for developing more mitochondrial genetic markers.

METHODS

The complete mitochondrial genome of Phlebotomus chinensis was amplified in eight fragments and sequenced using primer walking. The mitochondrial genome of Phlebotomus papatasi was reconstructed from whole-genome sequencing data available on Genbank. The phylogenetic relationship of 24 selected representatives of Diptera was deduced from codon positions 1 and 2 for 13 protein coding genes, using Bayesian inference (BI) and maximum likelihood (ML) methods.

RESULTS

We provide the first Phlebotomus (P. chinensis and P. papatasi) mitochondrial genomes. Both genomes contain 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and an A + T-rich region. The gene order of Phlebotomus mitochondrial genomes is identical with the ancestral gene order of insect. Phylogenetic analyses demonstrated that Psychodidae and Tanyderidae are sister taxa. Potential markers for population genetic study of Phlebotomus species were also revealed.

CONCLUSION

The generated mitochondrial genomes of P. chinensis and P. papatasi represent a useful resource for comparative genomic studies and provide valuable future markers for the population genetic study of these important Leishmania vectors. Our results also preliminary demonstrate the phylogenetic placement of Psychodidae based on their mitochondrial genomes.

Mitochondrial genome rearrangements at low taxonomic levels: three distinct mitogenome gene orders in the genus Pseudoniphargus (Crustacea: Amphipoda)

A comparison of mitochondrial genomes of three species of the amphipod Pseudoniphargus revealed the occurrence of a surprisingly high level of gene rearrangement involving protein-coding genes that is a rare phenomenon at low taxonomic levels. The three Pseudoniphargus mitogenomes also display a unique gene arrangement with respect to either the presumed Pancrustacean order or those known for other amphipods. Relative long non-coding sequences appear adjacent to the putative breakage points involved in gene rearrangements of protein coding genes. Other details of the newly obtained mitochondrial genomes - e.g., gene content, nucleotide composition and codon usage - are similar to those found in the mitogenomes of other amphipod species studied. They all contain the typical mitochondrial genome set consisting of 13 protein-coding genes, 22 tRNAs, and two rRNAS, as well as a large control region. The secondary structures and characteristics of tRNA and ribosomal mitochondrial genes of these three species are also discussed.

Culicoides Latreille (Diptera: Ceratopogonidae) taxonomy: current challenges and future directions

Culicoides Latreille biting midges (Diptera: Ceratopogonidae) cause a significant biting nuisance to humans, livestock and equines, and are the biological vectors of a range of internationally important pathogens of both veterinary and medical importance. Despite their economic significance, the delimitation and identification of species and evolutionary relationships between species within this genus remains at best problematic. To date no phylogenetic study has attempted to validate the subgeneric classification of the genus and the monophyly of many of the subgenera remains doubtful. Many informal species groupings are also known to exist but few are adequately described, further complicating accurate identification. Recent contributions to Culicoides taxonomy at the species level have revealed a high correlation between morphological and molecular analyses although molecular analyses are revealing the existence of cryptic species. This review considers the methods for studying the systematics of Culicoides using both morphological and genetic techniques, with a view to understanding the factors limiting our current understanding of Culicoides biology and hence arbovirus epidemiology. In addition, we examine the global status of Culicoides identification, highlighting areas that are poorly addressed, including the potential implementation of emerging technologies.

Biting midges monitoring (Diptera: Ceratopogonidae: Culicoides Latreille) in the governate of Monastir (Tunisia): species composition and molecular investigations

The results of entomological studies carried out in the governate of Monastir (Tunisia) in 2009-2010 (captures and emergences from muds) focusing on Culicoides species are presented in the present study. Identification of Culicoides at the species level is based on morphological characters, and a molecular study has been carried out based on mitochondrial DNA cytochrome C oxidase I gene (COI) and D1 and D2 domains of the 28S rDNA. The DNA sequences reported here are related to 10 species (on 25 known) of Culicoides described in Tunisia: Culicoides cataneii-gejgelensis, Culicoides circumscriptus, Culicoides imicola, Culicoides jumineri, Culicoides kingi, Culicoides langeroni, Culicoides newsteadi, Culicoides paolae, Culicoides puncticollis and Culicoides sahariensis. DNA sequencing of the COI gene and D1D2 domains discriminated all morphologically determined species. The choice of D1D2 domains considered as a conserved region is informative for Culicoides species identification. The molecular analyses of COI has grouped both C. circumscriptus, C. puncticollis within two clusters and C. newsteadi within five subclusters. However, C. newsteadi shows relatively deep intraspecific divergence using COI sequences.

Barcoding of biting midges in the genus Culicoides: a tool for species determination

Biting midges of the genus Culicoides (Diptera: Ceratopogonidae) are insect vectors of economically important veterinary diseases such as African horse sickness virus and bluetongue virus. However, the identification of Culicoides based on morphological features is difficult. The sequencing of mitochondrial cytochrome oxidase subunit I (COI), referred to as DNA barcoding, has been proposed as a tool for rapid identification to species. Hence, a study was undertaken to establish DNA barcodes for all morphologically determined Culicoides species in Swedish collections. In total, 237 specimens of Culicoides representing 37 morphologically distinct species were used. The barcoding generated 37 supported clusters, 31 of which were in agreement with the morphological determination. However, two pairs of closely related species could not be separated using the DNA barcode approach. Moreover, Culicoides obsoletus Meigen and Culicoides newsteadi Austen showed relatively deep intraspecific divergence (more than 10 times the average), which led to the creation of two cryptic species within each of C. obsoletus and C. newsteadi. The use of COI barcodes as a tool for the species identification of biting midges can differentiate 95% of species studied. Identification of some closely related species should employ a less conserved region, such as a ribosomal internal transcribed spacer.

Molecular identification of field-collected Culicoides larvae in the southern part of Japan

Although Culicoides biting midges act as a vector of important human and domestic animal diseases, their ecology is poorly understood. The lack of proper identification systems of Culicoides larvae is one of the main obstacles to progress in research. Based on mitochondrial sequences of 19 Japanese Culicoides species, we designed a universal primer set to amplify the partial sequence of the mitochondrial cytochrome c oxidase I (cox 1). The polymerase chain reaction product amplified from extracted DNA of Culicoides larvae using the primer set was directly sequenced, and species identification based on the variation at cox1 was conducted. Using the molecular identification system, we sorted 243 specimens of field-collected larvae from the southern part of Japan into 10 species including Culicoides arakawae (Arakawa), Culicoides oxystoma Kieffer, and Culicoides brevitarsis Kieffer, which are regarded as vectors of important livestock animal diseases. Eight species of Culicoides larvae, including C. arakawae and C. oxystoma, were recovered from active paddy fields and an abandoned paddy field. The result suggests that paddy fields contribute to breeding a variety of Culicoides species and maintenance and spread of Culicoides-borne pathogens. In contrast, larvae of C. brevitarsis were collected from cattle dung in pastures. The molecular identification system described herein using nucleotide sequences successfully achieved larval identification and will be useful for a better understanding of larval habitats of Culicoides biting midges.

The mitochondrial genome of Elodia flavipalpis Aldrich (Diptera: Tachinidae) and the evolutionary timescale of Tachinid flies

Tachinid flies are natural enemies of many lepidopteran and coleopteran pests of forests, crops, and fruit trees. In order to address the lack of genetic data in this economically important group, we sequenced the complete mitochondrial genome of the Palaearctic tachinid fly Elodia flavipalpis Aldrich, 1933. Usually found in Northern China and Japan, this species is one of the primary natural enemies of the leaf-roller moths (Tortricidae), which are major pests of various fruit trees. The 14,932-bp mitochondrial genome was typical of Diptera, with 13 protein-coding genes, 22 tRNA genes, and 2 rRNA genes. However, its control region is only 105 bp in length, which is the shortest found so far in flies. In order to estimate dipteran evolutionary relationships, we conducted a phylogenetic analysis of 58 mitochondrial genomes from 23 families. Maximum-likelihood and Bayesian methods supported the monophyly of both Tachinidae and superfamily Oestroidea. Within the subsection Calyptratae, Muscidae was inferred as the sister group to Oestroidea. Within Oestroidea, Calliphoridae and Sarcophagidae formed a sister clade to Oestridae and Tachinidae. Using a Bayesian relaxed clock calibrated with fossil data, we estimated that Tachinidae originated in the middle Eocene.

Molecular identification, phylogenetic status, and geographic distribution of Culicoides oxystoma (Diptera: Ceratopogonidae) in Israel

Culicoides oxystoma (Diptera: Ceratopogonidae) is an important vector species, reported mainly from Asia, with high potential to transmit viral diseases affecting livestock. In Japan, many arboviruses have been isolated from C. oxystoma, suggesting it as a key player in the epidemiology of several Culicoides-borne diseases. Over the years, C. oxystoma has also been reported in the Middle East region, including Israel. In this region, however, C. oxystoma cannot be easily distinguished morphologically from its sibling species included in the Culicoides schultzei complex. We therefore used genomic data for species identification and phylogeny resolution. Phylogenetic analyses based on internal transcribed spacer 1 (ITS-1) of ribosomal DNA and the mitochondrial gene encoding cytochrome oxidase subunit I (COI) showed that C. oxystoma from Israel is closely related to C. oxystoma from Japan. Using differential probing PCR, we showed that C. oxystoma is distributed all over the country, especially in Mediterranean climate regions. Culicoides oxystoma is less common or even absent in arid regions, while the other genetic cluster of C. schultzei complex was found only in the east of the country (mostly arid and semiarid regions). The molecular finding of C. oxystoma in wide geographical regions, together with its high proportion in the general Culicoides population and its vectoring potential, imply that it may be an important vector species in the Middle East.

The first mitochondrial genome for the wasp superfamily Platygastroidea: the egg parasitoid Trissolcus basalis

The nearly complete mitochondrial (mt) genome of an egg parasitoid, Trissolcus basalis (Wollaston), was sequenced using both 454 and Illumina next-generation sequencing technologies. A portion of the noncoding region remained unsequenced, possibly owing to the presence of repeats. The sequenced portion of the genome is 15,768 bp and has a high A+T content (84.2%), as is typical for hymenopteran mt genomes. A total of 36 of the 37 genes normally present in animal mt genomes were located. The one exception was trnR; a truncated version of this gene is present between trnS(1) and nd5, but it is unclear whether this gene fragment could code for the entire trnR gene. The mt gene arrangement of T. basalis is different from other Proctotrupomorpha mt genomes, with a number of trn genes in different positions. However, no shared derived gene rearrangements were identified in the present study. Bayesian analyses of mt genomes from 29 hymenopteran taxa and seven other orders of holometabolous insects support some uncontroversial evolutionary relationships, but indicate that much higher levels of taxonomic sampling are necessary for the resolution of family and superfamily relationships.

The mitochondrial genome of the 'twisted-wing parasite' Mengenilla australiensis (Insecta, Strepsiptera): a comparative study

BACKGROUND

Strepsiptera are an unusual group of sexually dimorphic, entomophagous parasitoids whose evolutionary origins remain elusive. The lineage leading to Mengenilla australiensis (Family Mengenillidae) is the sister group to all remaining extant strepsipterans. It is unique in that members of this family have retained a less derived condition, where females are free-living from pupation onwards, and are structurally much less simplified. We sequenced almost the entire mitochondrial genome of M. australiensis as an important comparative data point to the already available genome of its distant relative Xenos vesparum (Family Xenidae). This study represents the first in-depth comparative mitochondrial genomic analysis of Strepsiptera.

RESULTS

The partial genome of M. australiensis is presented as a 13421 bp fragment, across which all 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes and 18 transfer RNA (tRNA) sequences are identified. Two tRNA translocations disrupt an otherwise ancestral insect mitochondrial genome order. A+T content is measured at 84.3%, C-content is also very skewed. Compared with M. australiensis, codon bias in X. vesparum is more balanced. Interestingly, the size of the protein coding genome is truncated in both strepsipterans, especially in X. vesparum which, uniquely, has 4.3% fewer amino acids than the average holometabolan complement. A revised assessment of mitochondrial rRNA secondary structure based on comparative structural considerations is presented for M. australiensis and X. vesparum.

CONCLUSIONS

The mitochondrial genome of X. vesparum has undergone a series of alterations which are probably related to an extremely derived lifestyle. Although M. australiensis shares some of these attributes; it has retained greater signal from the hypothetical most recent common ancestor (MRCA) of Strepsiptera, inviting the possibility that a shift in the mitochondrial selective environment might be related to the specialization accompanying the evolution of a small, morphologically simplified completely host-dependent lifestyle. These results provide useful insights into the nature of the evolutionary transitions that accompanied the emergence of Strepsiptera, but we emphasize the need for adequate sampling across the order in future investigations concerning the extraordinary developmental and evolutionary origins of this group.

Prevalence of Cardinium bacteria in planthoppers and spider mites and taxonomic revision of "Candidatus Cardinium hertigii" based on detection of a new Cardinium group from biting midges

Cardinium bacteria, members of the phylum Cytophaga-Flavobacterium-Bacteroides (CFB), are intracellular bacteria in arthropods that are capable of inducing reproductive abnormalities in their hosts, which include parasitic wasps, mites, and spiders. A high frequency of Cardinium infection was detected in planthoppers (27 out of 57 species were infected). A high frequency of Cardinium infection was also found in spider mites (9 out of 22 species were infected). Frequencies of double infection by Cardinium and Wolbachia bacteria (Alphaproteobacteria capable of manipulating reproduction of their hosts) were disproportionately high in planthoppers but not in spider mites. A new group of bacteria, phylogenetically closely related to but distinct from previously described Cardinium bacteria (based on 16S rRNA and gyrB genes) was found in 4 out of 25 species of Culicoides biting midges. These bacteria possessed a microfilament-like structure that is a morphological feature previously found in Cardinium and Paenicardinium. The bacteria close to the genus Cardinium consist of at least three groups, A, B, and C. Group A is present in various species of arthropods and was previously referred to as "Candidatus Cardinium hertigii," group B is present in plant parasitic nematodes and was previously referred to as "Candidatus Paenicardinium endonii," and group C is present in Culicoides biting midges. On the basis of morphological and molecular data, we propose that the nomenclature of these three groups be integrated into a single species, "Candidatus Cardinium hertigii."

Characterization of internal transcribed spacer (ITS1)-ITS2 region of ribosomal RNA gene from 25 species of Culicoides biting midges (Diptera: Ceratopogonidae) in Japan

We determined nucleotide sequences of the nuclear rDNA internal transcribed spacer (ITS)1-5.8S-ITS2a-2S-ITS2 region in 103 individuals of 25 Culicoides species (Diptera: Ceratopogonidae) from 11 locations in Japan. Ribosomal RNA genes, 5.8S and 2S rDNA, were highly conserved among the species with few variations. The ITS2a region showed length variation among species. Both ITS1 and ITS2 showed highly varied sequences among species. The noticeable indel regions among ITS1 sequences are present in some Culicoides species, separating species into two types having long or short ITS1 region. However, Culicoides cylindratus Kitaoka possesses both types of ITS1 in each individual; these results seem to indicate that the ITS1-long type was the prototype and the short type was produced through deletion in many Culicoides species. One species, belonging to subgenus Avaritia, possessed an Avaritia-specific sequence in ITS1 and phylogenetically formed a monophyletic group. Geographical genotypes in a species were not clear. Species-specific sequence features were observed, enabling molecular identification of Culicoides species.