Analysis of the complete mitochondrial DNA from Anopheles funestus: an improved dipteran mitochondrial genome annotation and a temporal dimension of mosquito evolution

Phylogenetic taxonomy of the Zambian Anopheles coustani group using a mitogenomics approach

Background

Mosquito species belonging to the Anopheles coustani group have been implicated in driving residual malaria transmission in sub-Saharan Africa and are regarded as an established primary vector in Madagascar. The morphological identification of mosquitoes in this group is challenging due to cryptic features and their molecular confirmation is difficult due to a paucity of reference sequence data representing all members of the group. Conventional molecular barcoding with the cytochrome oxidase I (COI) gene and the internal transcribed spacer 2 (ITS2) region targets is limited in their discrimination and conclusive identification of members of species complexes. In contrast, complete mitochondrial genomes (mitogenomes) have demonstrated much improved power over barcodes to be useful in rectifying taxonomic discrepancies in Culicidae.

Methods

We utilized a genome skimming approach via shallow shotgun sequencing on individual mosquito specimens to generate sequence reads for mitogenome assembly. Bayesian inferred phylogenies and molecular dating estimations were perfomed on the concatenated protein coding genes using the Bayesian Evolutionary Analysis by Sampling Trees 2 (BEAST 2) platform. Divergence estimates were calibrated on published calucations for Anopheles-Aedes.

Results

This study generated 17 new complete mitogenomes which were comprable to reference An. coustani mitogenomes in the GenBank repository by having 13 protein coding, 22 transfer RNA and 2 ribosomal RNA genes, with an average length of 15,400 bp and AT content of 78.3%. Bayesian inference using the concatenated protein coding genes from the generated and publicly available mitogenomes yielded six clades: one for each of the four taxa targeted in this study, the GenBank references, and a currently unknown species. Divergence times estimated that the An. coustani group separated from the An. gambiae complex approximately 110 million years ago (MYA), and members within the complex diverged at times points ranging from~34 MYA to as recent as ~7 MYA.

Conclusions

These findings demonstrate the value of mitochondrial genomes in differentiating cryptic taxa and help to confirm morphological identities of An. coustani s.s., An. paludis, An. zeimanni and An. tenebrosus. Divergence estimates with the An. coustani group are similar to those for well-studied anopheline vector groups. These analyses also highlight the likely prescence of other cryptic An. coustani group members circulating in Zambia.

Unravelling the complete mitochondrial genomes of Thrips tabaci Lindeman and Thrips parvispinus Karny (Thysanoptera: Thripidae) and their phylogenetic implications

Onion (Allium cepa Linnaeus) is an important vegetable crop valued for its nutritional properties and economics worldwide. Onion cultivation faces serious threats from pests and diseases, particularly onion thrips (Thrips tabaci), which cause substantial yield losses. Recently, Black thrips (Thrips parvispinus), an invasive key pest of chili, have been reported to cause severe damage in onion crop and is likely to devastate the onion cultivation in near future. Therefore, this study was conducted to address the knowledge gap concerning the genetic basis and evolutionary history of T. tabaci and T. parvispinus through sequencing of their mitochondrial genomes. T. tabaci and T. parvispinus were collected from different locations in Maharashtra, India, and reared in the laboratory. The mitochondrial genomes of T. tabaci and T. parvispinus were sequenced to a length of 15,277 and 15,285 bp, respectively. Both genomes exhibited similar gene organization with regard to thirteen protein-coding genes and two rRNA genes. T. tabaci contained 19 tRNA genes whereas T. parvispinus contained 18 tRNA genes. The evolutionary positions of T. tabaci and T. parvispinus within the Thysanoptera order were elucidated through phylogenetic analysis of the mitogenomes of 15 thrips species. These findings provide crucial insights into the genetic makeup and evolutionary dynamics of both the thrips species, thereby aiding the development of novel and sustainable pest management strategies to mitigate their impacts on crops in the changing climate scenario.

Seven mitochondrial genomes of tribe Hylurgini (Coleoptera: Curculionidae: Scolytinae) in Eurasia and their phylogenetic analysis

The Hylurgini tribe (Coleoptera: Curculionidae: Scolytinae) comprises commercially significant bark beetles, including invasive species within the genera Dendroctonus and Hylurgus. These invasive species coexist with native Tomicus species of Hylurgini and cooperatively infest host trees in China. However, we lack sufficient mitochondrial genome data of Hylurgini to conduct phylogenetic studies, clarify the phylogenetic relationships of the above species, and improve the understanding of niche divergence and common hazards. Here, we sequenced and analyzed the mitochondrial genomes of seven Hylurgini species, including Dendroctonus valens, Hylurgus ligniperda, Hylurgus micklitzi, Tomicus piniperda, Tomicus brevipilosus, Tomicus minor and Tomicus yunnanensis. All sequenced mitochondrial genomes ranged from 15,339 bp to 17,545 bp in length, and their AT contents ranged from 73.24% to 78.81%. The structure of the seven mitochondrial genomes was consistent with that of ancestral insects. Based on 13 protein-coding genes from the reported mitochondrial genomes of 29 species of bark beetles, we constructed phylogenetic trees using maximum likelihood and Bayesian inference methods. The topology of the two phylogenetic trees was almost consistent. The findings elucidated the taxonomy classification of Hylurgini and the evolutionary connections of its sister taxa within the Scolytinae. This study offers insights for examining the evolutionary connections between invasive and native bark beetles, as well as the molecular identification and detection of newly invading species.

The first two complete mitochondrial genomes for the genus Anagyrus (Hymenoptera, Encyrtidae) and their phylogenetic implications

Anagyrus, a genus of Encyrtidae (Hymenoptera, Chalcidoidea), represents a successful group of parasitoid insects that attack various mealybug pests of agricultural and forestry plants. Until now, only 20 complete mitochondrial genomes have been sequenced, including those in this study. To enrich the diversity of mitochondrial genomes in Encyrtidae and to gain insights into their phylogenetic relationships, the mitochondrial genomes of two species of Anagyrus were sequenced, and the mitochondrial genomes of these species were compared and analyzed. Encyrtid mitochondrial genomes exhibit similarities in nucleotide composition, gene organization, and control region patterns. Comparative analysis of protein-coding genes revealed varying molecular evolutionary rates among different genes, with six genes (ATP8, ND2, ND4L, ND6, ND4 and ND5) showing higher rates than others. A phylogenetic analysis based on mitochondrial genome sequences supports the monophyly of Encyrtidae; however, the two subfamilies, Encyrtinae and Tetracneminae, are non-monophyletic. This study provides valuable insights into the phylogenetic relationships within the Encyrtidae and underscores the utility of mitochondrial genomes in the systematics of this family.

Characterization of the complete mitochondrial genome of Paecilomyces variotii and comparative evolutionary mitochondriomics of 36 fungi

BACKGROUD

Paecilomyces variotii has important economic value in stimulating crop growth, biodegradation, and other aspects. Up to now, there are no research reports on its mitochondrial genome.

METHODS AND RESULTS

The mitochondrial genome of Paecilomyces variotii was determined with the next-generation sequencing method (Illumina, NovaSeq), and its characteristics were analyzed using various bioinformatics approaches. The length of complete mitochondrial genome sequence of P. variotii is 40,965 bp and consists of 14 protein-coding genes, 2 ribosomal RNA genes, 1 ribosomal protein S3 gene, 26 transport RNA genes. The results of phylogenetics analysis using Bayesian inference and Maximum likelihood methods showed that P. variotii belongs to the Eurotiales order in the Thermoascaceae family, and 9 genera within the Eurotiomycetes class were effectively distinguished with high support rates (bootstrap value > 92% and posterior probabilities > 99%). The analysis of synonymous substitution rates and nonsynonymous substitution rates indicated that the Ka/Ks values of the 14 PCGs in the mitochondrial genomes of the two orders in the Eurotiomycetes class ranged from 0 to 0.4333.

CONCLUSIONS

This study revealed the structural and sequence information characteristics of the mitochondrial genome of P. variotii, and the phylogenetic results strongly support its classification within the family Thermoascaceae, consistent with traditional morphological taxonomy studies. The 14 PCGs in the mitochondrial genomes of the two orders in the Eurotiomycetes class are subject to strong purifying (negative) selection. The results of this research provides an important molecular basis for the development of genomics, evolutionary genetics and molecular markers of P. variotii in the future.

Serine protease homolog pairs CLIPA4-A6, A4-A7Δ, and A4-A12 act as cofactors for proteolytic activation of prophenoloxidase-2 and -7 in Anopheles gambiae

Phenoloxidase (PO) catalyzed melanization and other insect immune responses are mediated by serine proteases (SPs) and their noncatalytic homologs (SPHs). Many of these SP-like proteins have a regulatory clip domain and are called CLIPs. In most insects studied so far, PO precursors are activated by a PAP (i.e., PPO activating protease) and its cofactor of clip-domain SPHs. Although melanotic encapsulation is a well-known refractory mechanism of mosquitoes against malaria parasites, it is unclear if a cofactor is required for PPO activation. In Anopheles gambiae, CLIPA4 is 1:1 orthologous to Manduca sexta SPH2; CLIPs A5-7, A12-14, A26, A31, A32, E6, and E7 are 11:4 orthologous to M. sexta SPH1a, 1b, 4, and 101, SPH2 partners in the cofactors. Here we produced proCLIPs A4, A6, A7Δ, A12, and activated them with CLIPB9 or M. sexta PAP3. A. gambiae PPO2 and PPO7 were expressed in Escherichia coli for use as PAP substrates. CLIPB9 was mutated to CLIPB9Xa by including a Factor Xa cleavage site. CLIPA7Δ was a deletion mutant with a low complexity region removed. After PAP3 or CLIPB9Xa processing, CLIPA4 formed a high Mr complex with CLIPA6, A7Δ or A12, which assisted PPO2 and PPO7 activation. High levels of specific PO activity (55-85 U/μg for PO2 and 1131-1630 U/μg for PO7) were detected in vitro, indicating that cofactor-assisted PPO activation also occurs in this species. The cleavage sites and mechanisms for complex formation and cofactor function are like those reported in M. sexta and Drosophila melanogaster. In conclusion, these data suggest that the three (and perhaps more) SPHI-II pairs may form cofactors for CLIPB9-mediated activation of PPOs for melanotic encapsulation in A. gambiae.

Recognition of polymorphic Csd proteins determines sex in the honeybee

Sex in honeybees, Apis mellifera, is genetically determined by heterozygous versus homo/hemizygous genotypes involving numerous alleles at the single complementary sex determination locus. The molecular mechanism of sex determination is however unknown because there are more than 4950 known possible allele combinations, but only two sexes in the species. We show how protein variants expressed from complementary sex determiner (csd) gene determine sex. In females, the amino acid differences between Csd variants at the potential-specifying domain (PSD) direct the selection of a conserved coiled-coil domain for binding and protein complexation. This recognition mechanism activates Csd proteins and, thus, the female pathway. In males, the absence of polymorphisms establishes other binding elements at PSD for binding and complexation of identical Csd proteins. This second recognition mechanism inactivates Csd proteins and commits male development via default pathway. Our results demonstrate that the recognition of different versus identical variants of a single protein is a mechanism to determine sex.

Advances in the genetic characterization of the malaria vector, Anopheles funestus, and implications for improved surveillance and control

Anopheles mosquitoes present a major public health challenge in sub-Saharan Africa; notably, as vectors of malaria that kill over half a million people annually. In parts of the east and southern Africa region, one species in the Funestus group, Anopheles funestus, has established itself as an exceptionally dominant vector in some areas, it is responsible for more than 90% of all malaria transmission events. However, compared to other malaria vectors, the species is far less studied, partly due to difficulties in laboratory colonization and the unresolved aspects of its taxonomy and systematics. Control of An. funestus is also increasingly difficult because it has developed widespread resistance to public health insecticides. Fortunately, recent advances in molecular techniques are enabling greater insights into species identity, gene flow patterns, population structure, and the spread of resistance in mosquitoes. These advances and their potential applications are reviewed with a focus on four research themes relevant to the biology and control of An. funestus in Africa, namely: (i) the taxonomic characterization of different vector species within the Funestus group and their role in malaria transmission; (ii) insecticide resistance profile; (iii) population genetic diversity and gene flow, and (iv) applications of genetic technologies for surveillance and control. The research gaps and opportunities identified in this review will provide a basis for improving the surveillance and control of An. funestus and malaria transmission in Africa.

The transcriptional response in mosquitoes distinguishes between fungi and bacteria but not Gram types

Mosquitoes are prolific vectors of human pathogens; a clear and accurate understanding of the organization of their antimicrobial defenses is crucial for informing the development of transmission control strategies. The canonical infection response in insects, as described in the insect model Drosophila melanogaster , is pathogen type-dependent, with distinct stereotypical responses to Gram-negative bacteria and Gram-positive bacteria/fungi mediated by the activation of the Imd and Toll pathways, respectively. To determine whether this pathogen-specific discrimination is shared by mosquitoes, we used RNAseq to capture the genome-wide transcriptional response of Aedes aegypti and Anopheles gambiae ( s.l. ) to systemic infection with Gram-negative bacteria, Gram-positive bacteria, yeasts, and filamentous fungi, as well as challenge with heat-killed Gram-negative, Gram-positive, and fungal pathogens. From the resulting data, we found that Ae. aegypti and An. gambiae both mount a core response to all categories of infection, and this response is highly conserved between the two species with respect to both function and orthology. When we compared the transcriptomes of mosquitoes infected with different types of bacteria, we observed that the intensity of the transcriptional response was correlated with both the virulence and growth rate of the infecting pathogen. Exhaustive comparisons of the transcriptomes of Gram-negative-challenged versus Gram-positive-challenged mosquitoes yielded no difference in either species. In Ae. aegypti , however, we identified transcriptional signatures specific to bacterial infection and to fungal infection. The bacterial infection response was dominated by the expression of defensins and cecropins, while the fungal infection response included the disproportionate upregulation of an uncharacterized family of glycine-rich proteins. These signatures were also observed in Ae. aegypti challenged with heat-killed bacteria and fungi, indicating that this species can discriminate between molecular patterns that are specific to bacteria and to fungi.

Characterization of the complete mitochondrial genome of Culex vishnui (Diptera: Culicidae), one of the major vectors of Japanese encephalitis virus

Culex mosquitoes (Diptera: Culicidae) can transmit a variety of arthropod-borne viruses (arboviruses), causing human and animal diseases. Cx. vishnui, Cx. pseudovishnui, and Cx. tritaeniorhynchus are three representative species in Culex vishnui subgroup, which are widely distributed in southeast Asia, and they have been proved as the main vectors transmitting Japanese encephalitis virus (JEV) that could cause human infectious mosquito-borne disease across Asia. However, the epidemiology, biology, and even molecular information of those mosquitos remain poorly understood, and only the mitochondrial genome (mitogenome) of Cx. tritaeniorhynchus has been reported in these species. In the present study, we sequenced and annotated the complete mitogenome sequence of Cx. vishnui which was 15,587 bp in length, comprising 37 genes. Comparisons of nucleotide and amino acid sequences between Cx. vishnui and Cx. tritaeniorhynchus revealed that most genes within Culex vishnui subgroup were conserved, except atp8, nad1, atp6, and nad6, with differences of 0.4 (rrnS) - 15.1% (tRNAs) and 0 (nad4L) - 9.4% (atp8), respectively, interestingly suggesting the genes nad4L and rrnS were the most conserved but atp8 gene was the least. The results based on nucleotide diversity also supported a relatively uniform distribution of the intraspecific differences in Cx. vishnui and Cx. tritaeniorhynchus with only one highly pronounced peak of divergence centered at the control region. Phylogenetic analyses using concatenated amino acid sequences of 13 protein-coding genes supported the previous taxonomic classification of the family Culicidae and the monophyly of tribes Aedini, Culicini, Mansoniini, and Sabethini. The present study revealed detailed information on the subgroup Culex vishnui, reanalyzed the relationships within the family Culicidae, provided better markers to identify and distinguish Culex species, and offered more markers for studying the molecular epidemiology, population genetics, and molecular phylogenetics of Cx. vishnui.

Identification and functional characterization of olfactory indolergic receptors in Musca domestica

In mosquitoes, indolic compounds are detected by a group of olfactory indolergic Odorant Receptors (indolORs). The ancient origin of indole and 3-methylindole as chemical signals suggest that they may be detected by insects outside the Culicidae clade. To test this hypothesis, we have identified potential indolOR genes in brachyceran flies based on sequence homology. Because of the crucial roles of indolic compounds in oviposition and foraging, we have focused our attention on the housefly Musca domestica. Using a heterologous expression system, we have identified indolOR transcript expression in the female antennae, and have characterized MdomOR30a and MdomOR49b as 3-methylindole and indole receptors, respectively. We have identified a set of 92 putative indolOR genes encoded in the genomes of Culicoidea, Psychodidae and brachycera, described their phylogenetic relationships, and exon/intron structures. Further characterization of indolORs will impact our understanding of insect chemical ecology and will provide targets for the development of novel odor-based tools that can be integrated into existing vector surveillance and control programs.

Interspecific Hybridization and Complete Mitochondrial Genome Analysis of Two Ghost Moth Species

Simple Summary

The Chinese cordyceps is a valuable parasitic Ophiocordyceps sinensis fungus–Thitarodes/Hepialus larva complex. In view of culturing this complex, a method for the artificial rearing of the Thitarodes/Hepialus ghost moth hosts was established. Deterioration of the host insect population and low mummification rates in infected larvae constrain effective cultivation. Hybridization of Thitarodes/Hepialus populations may overcome this problem. Thitarodes shambalaensis and Thitarodes sp. were inbred or hybridized, and the biological parameters, larval sensitivity to the fungal infection and mitochondrial genomes of the resulting populations were investigated. Hybridization of T. shambalaensis and Thitarodes sp. allowed producing a new generation. One hybrid population (T. shambalaensis females mated with Thitarodes sp. males) showed increased population growth as compared with the parental Thitarodes sp. population. The sensitivity of the inbred larval populations to four fungal isolates of O. sinensis differed. The complete mitochondrial genomes of T. shambalaensis, Thitarodes sp. and the hybrid population were 15,612 bp, 15,389 bp and 15,496 bp in length, respectively. A + T-rich regions were variable in sizes and repetitive sequences. The hybrid population was located in the same clade with T. shambalaensis, implying the maternal inheritance of mitochondrial DNA.

Abstract

The Chinese cordyceps, a parasitic Ophiocordyceps sinensis fungus–Thitarodes/Hepialus larva complex, is a valuable biological resource endemic to the Tibetan Plateau. Protection of the Plateau environment and huge market demand make it necessary to culture this complex in an artificial system. A method for the large-scale artificial rearing of the Thitarodes/Hepialus insect host has been established. However, the deterioration of the insect rearing population and low mummification of the infected larvae by the fungus constrain effective commercial cultivation. Hybridization of Thitarodes/Hepialus populations may be needed to overcome this problem. The species T. shambalaensis (GG♂ × GG♀) and an undescribed Thitarodes species (SD♂ × SD♀) were inbred or hybridized to evaluate the biological parameters, larval sensitivity to the fungal infection and mitochondrial genomes of the resulting populations. The two parental Thitarodes species exhibited significant differences in adult fresh weights and body lengths but not in pupal emergence rates. Hybridization of T. shambalaensis and Thitarodes sp. allowed producing a new generation. The SD♂ × GG♀ population showed a higher population trend index than the SD♂ × SD♀ population, implying increased population growth compared with the male parent. The sensitivity of the inbred larval populations to four fungal isolates of O. sinensis also differed. This provides possibilities to create Thitarodes/Hepialus populations with increased growth potential for the improved artificial production of the insect hosts. The mitochondrial genomes of GG♂ × GG♀, SD♂ × SD♀ and SD♂ × GG♀ were 15,612 bp, 15,389 bp and 15,496 bp in length, with an A + T content of 80.92%, 82.35% and 80.87%, respectively. The A + T-rich region contains 787 bp with two 114 bp repetitive sequences, 554 bp without repetitive sequences and 673 bp without repetitive sequences in GG♂ × GG♀, SD♂ × SD♀ and SD♂ × GG♀, respectively. The hybrid population (SD♂ × GG♀) was located in the same clade with GG♂ × GG♀, based on the phylogenetic tree constructed by 13 PCGs, implying the maternal inheritance of mitochondrial DNA.

Updating the bionomy and geographical distribution of Anopheles (Nyssorhynchus) albitarsis F: A vector of malaria parasites in northern South America

Anopheles albitarsis F is a putative species belonging to the Albitarsis Complex, recognized by rDNA, mtDNA, partial white gene, and microsatellites sequences. It has been reported from the island of Trinidad, Venezuela and Colombia, and incriminated as a vector of malaria parasites in the latter. This study examined mitochondrially encoded cytochrome c oxidase I (MT-CO1) sequences of An. albitarsis F from malaria-endemic areas in Colombia and Venezuela to understand its relations with other members of the Complex, revised and update the geographical distribution and bionomics of An. albitarsis F and explore hypotheses to explain its phylogenetic relationships and geographical expansion. Forty-five MT-CO1 sequences obtained in this study were analyzed to estimate genetic diversity and possible evolutionary relationships. Sequences generated 37 haplotypes clustered in a group where the genetic divergence of Venezuelan populations did not exceed 1.6% with respect to Colombian samples. Anopheles albitarsis F (π = 0.013) represented the most recent cluster located closer to An. albitarsis I (π = 0.009). Barcode gap was detected according to Albitarsis Complex lineages previously reported (threshold 0.014–0.021). Anopheles albitarsis F has a wide distribution in northern South America and might play an important role in the transmission dynamics of malaria due to its high expansion capacity. Future studies are required to establish the southern distribution of An. albitarsis F in Venezuela, and its occurrence in Guyana and Ecuador.

The complete mitochondrial genome of Trifida elongate and comparative analysis of 43 leafhoppers

In this study, the mitochondrial genome of Trifida elongate was sequenced, and comparative analysis of T. elongate and other 43 leafhoppers was performed based on the mitochondrial genome. The mitochondrial genome sequence length of T. elongate was 14,924 bp. It comprised 13 protein-coding genes, 2 ribosomal RNA (rRNA) genes, 22 transport RNA (tRNA) genes, and 1 non-coding control region. The control region is located between the rrnS and trnI genes, is characterized by two tandem repeats and three simple sequence repeats. Phylogenetic analysis showed that T. elongate is closely related with Bolanusoides shaanxiensis and Limassolla lingchuanensis (bootstrap value = 92% and posterior probabilities = 1). Analysis of synonymous and non-synonymous nucleotide substitutions showed that Ka/Ks value of the 13 protein-coding genes of 8 subfamily leafhoppers were less than 1 ranging from 0.0315 to 0.9928. atp8 had the highest Ka/Ks value whereas cox1 had the lowest Ka/Ks value. This study provides information on the structure and sequence characteristics of the mitochondrial genome of T. elongata. Typhlocybinae is clustered with (Cicadellinae+(Idiocerinae+(Mileewinae+(Nirvaninae+(Evacanthinae+Ledrinae))))).

Suppressed expression of oxidoreductin-like protein, Oxidor, increases follicle degeneration and decreases survival during the overwintering diapause of the mosquito Culex pipiens

Throughout diapause in mosquitoes, stress resistance and subsequent prolonged lifespan are a few important features of diapause that are crucial for overwintering success. In the mosquito Culex pipiens, we suggest that oxidoreductin-like protein is involved with these diapause characteristics for overwintering survival. Expression of oxidor was more than two-fold higher in early stage diapausing females compared to their non-diapausing counterparts. Suppression of the gene that encodes oxidoreductin-like protein by RNAi significantly increased the proportion of degenerating follicles in early-stage adult diapausing females. Inhibition of oxidor also significantly reduced the survivability of diapausing females which indicates that this protein plays a key role in protecting multiple tissues during early diapause.

Utility of Complete Mitochondrial Genomes in Phylogenetic Classification of the Species of Anopheles (Culicidae: Anophelinae)

Background:

Among the blood-sucking insects, Anopheles mosquitoes have a very special position, because they transmit parasites of the genus Plasmodium, which cause malaria as one of the main vector-borne disease worldwide. The aim of this review study was to evaluate utility of complete mitochondrial genomes in phylogenetic classification of the species of Anopheles.

Methods:

The complete mitochondrial genome sequences belonging to 28 species of the genus Anopheles (n=32) were downloaded from NCBI. The phylogenetic trees were constructed using the ML, NJ, ME, and Bayesian inference methods.

Results:

In general, the results of the present survey revealed that the complete mitochondrial genomes act very accurately in recognition of the taxonomic and phylogenetic status of these species and provide a higher level of support than those based on individual or partial mitochondrial genes so that by using them, we can meticulously reconstruct and modify Anopheles classification.

Conclusion:

Understanding the taxonomic position of Anopheles, can be a very effective step in better planning for controlling these malaria vectors in the world and will improve our knowledge of their evolutionary biology.

femaleless Controls Sex Determination and Dosage Compensation Pathways in Females of Anopheles Mosquitoes

The insect sex determination and the intimately linked dosage compensation pathways represent a challenging evolutionary puzzle that has been solved only in Drosophila melanogaster. Analyses of orthologs of the Drosophila genes identified in non-drosophilid taxa^1,2 revealed that evolution of sex determination pathways is consistent with a bottom-up mode,³ where only the terminal genes within the pathway are well conserved. doublesex (dsx), occupying a bottom-most position and encoding sex-specific proteins orchestrating downstream sexual differentiation processes, is an ancient sex-determining gene present in all studied species.^2,4,5 With the exception of lepidopterans, its female-specific splicing is known to be regulated by transformer (tra) and its co-factor transformer-2 (tra2).6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20 Here we show that in the African malaria mosquito Anopheles gambiae, a gene, which likely arose in the Anopheles lineage and which we call femaleless (fle), controls sex determination in females by regulating splicing of dsx and fruitless (fru; another terminal gene within a branch of the sex determination pathway). Moreover, fle represents a novel molecular link between the sex determination and dosage compensation pathways. It is necessary to suppress activation of dosage compensation in females, as demonstrated by the significant upregulation of the female X chromosome genes and a correlated female-specific lethality, but no negative effect on males, in response to fle knockdown. This unexpected property, combined with a high level of conservation in sequence and function in anopheline mosquitoes, makes fle an excellent target for genetic control of all major vectors of human malaria.

•

fle is a new sex determination pathway element conserved in Anopheles mosquitoes

•

fle may have originated in the Anopheles lineage and is highly conserved in Anopheles

•

fle suppresses activation of dosage compensation in females

•

Depletion of fle transcripts is lethal or otherwise deleterious to females

Ultraconserved Non-coding DNA Within Diptera and Hymenoptera

This study has taken advantage of the availability of the assembled genomic sequence of flies, mosquitos, ants and bees to explore the presence of ultraconserved sequence elements in these phylogenetic groups. We compared non-coding sequences found within and flanking Drosophila developmental genes to homologous sequences in Ceratitis capitata and Musca domestica. Many of the conserved sequence blocks (CSBs) that constitute Drosophila cis-regulatory DNA, recognized by EvoPrinter alignment protocols, are also conserved in Ceratitis and Musca. Also conserved is the position but not necessarily the orientation of many of these ultraconserved CSBs (uCSBs) with respect to flanking genes. Using the mosquito EvoPrint algorithm, we have also identified uCSBs shared among distantly related mosquito species. Side by side comparison of bee and ant EvoPrints of selected developmental genes identify uCSBs shared between these two Hymenoptera, as well as less conserved CSBs in either one or the other taxon but not in both. Analysis of uCSBs in these dipterans and Hymenoptera will lead to a greater understanding of their evolutionary origin and function of their conserved non-coding sequences and aid in discovery of core elements of enhancers.

This study applies the phylogenetic footprinting program EvoPrinter to detection of ultraconserved non-coding sequence elements in Diptera, including flies and mosquitos, and Hymenoptera, including ants and bees. EvoPrinter outputs an interspecies comparison as a single sequence in terms of the input reference sequence. Ultraconserved sequences flanking known developmental genes were detected in Ceratitis and Musca when compared with Drosophila species, in Aedes and Culex when compared with Anopheles, and between ants and bees. Our methods are useful in detecting and understanding the core evolutionarily hardened sequences required for gene regulation.

Mitochondrial genes associated with pyrethroid resistance revealed by mitochondrial genome and transcriptome analyses in the malaria vector Anopheles sinensis (Diptera: Culicidae)

BACKGROUND

Insecticides are still the main method of mosquito control, but mosquito resistance presents a large obstacle. The function of mitochondrial genes in the evolution of insecticide resistance is still poorly understood. Pyrethroid is the most commonly used insecticide, and Anopheles sinensis is an important malaria vector in China and Southeast Asia. In this study, we investigated the mitochondrial genes associated with pyrethroid resistance through their genetic and expression variation based on analyses of transcriptomes and 36 individuals with resequencing in three geographical populations in China.

RESULTS

The nucleotide diversity (Pi) in 18 resistant individuals was much lower than that in 18 susceptible individuals, which suggests that some sites experienced purifying selection subject to pyrethroid stress. Ka/Ks and amino acid analyses showed that ND4 experienced positive selection and had 23 amino acid mutations due to pyrethroid stress. These mutations might change the ND4 structure and function and thus alter the efficiency of the respiratory chain. ND5 was significantly upregulated, and ATP8 was significantly downregulated in these three pyrethroid resistant populations, which suggests that these two genes function in the production and maintenance of pyrethroid resistance. There are differences in mitochondrial genes involved in pyrethroid resistance among these three populations.

CONCLUSION

This is the first study to reveal the association of mitochondrial genes in the evolution of insecticide resistance through amino acid mutation and expression patterns and can help us further understand insecticide resistance mechanisms. © 2019 Society of Chemical Industry.

Five mitochondrial genomes of black fungus gnats (Sciaridae) and their phylogenetic implications

Sciaridae is a family of great species diversity, distributed worldwide, that includes important agricultural pests of cultivated mushrooms and plants produced in greenhouses. Here we sequenced five nearly complete mitochondrial genomes representing three subfamilies of Sciaridae. The lengths of these mitogenomes range from 13,849 bp to 16,923 bp with 13 protein-coding genes (PCGs), 20-22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a control region (CR). Compared with other dipteran species, rearrangements in Sciaridae are more common. Inversion or transition is observed frequently of trnL2, and in the tRNA clusters trnI-trnQ-trnM, trnW-trnC-trnY, and trnA-trnR-trnN-trnS1-trnE-trnF. Phylogenetic relationships within the family were reconstructed based on these newly sequenced species, combined with the published mitogenomes of related families, and recovered the topology within Sciaroidea as Cecidomyiidae + (Sciaridae + Keroplatidae). Relationships recovered within Sciaridae were Sciarinae + ('Pseudolycoriella group' + Megalosphyinae).

Mitochondrial genome characterization of the family Trigonidiidae (Orthoptera) reveals novel structural features and nad1 transcript ends

The Trigonidiidae, a family of crickets, comprises 981 valid species with only one mitochondrial genome (mitogenome) sequenced to date. To explore mitogenome features of Trigonidiidae, six mitogenomes from its two subfamilies (Nemobiinae and Trigonidiinae) were determined. Two types of gene rearrangements involving a trnN-trnS1-trnE inversion and a trnV shuffling were shared by Trigonidiidae. A long intergenic spacer was observed between trnQ and trnM in Trigonidiinae (210-369 bp) and Nemobiinae (80-216 bp), which was capable of forming extensive stem-loop secondary structures in Trigonidiinae but not in Nemobiinae. The anticodon of trnS1 was TCT in Trigonidiinae, rather than GCT in Nemobiinae and other related subfamilies. There was no overlap between nad4 and nad4l in Dianemobius, as opposed to a conserved 7-bp overlap commonly found in insects. Furthermore, combined comparative analysis and transcript verification revealed that nad1 transcripts ended with a U, corresponding to the T immediately preceding a conserved motif GAGAC in the superfamily Grylloidea, plus poly-A tails. The resultant UAA served as a stop codon for species lacking full stop codons upstream of the motif. Our findings gain novel understanding of mitogenome structural diversity and provide insight into accurate mitogenome annotation.

The complete mitochondrial genome of a field cricket Turanogryllus eous (Insecta: Orthoptera)

The complete mitochondrial genome (mitogenome) of a field cricket Turanogryllus eous Bey-Bienko, 1956 was determined using next-generation sequencing. The mitogenome was 16,045 bp in length comprising 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a control region. Relative to the ancestral insect gene order, T. eous possessed an inversion of trnN-trnS1-trnE. The control region contained 3.4 tandem copies of a 194-bp sequence. Phylogenetic analysis supported that T. eous was sister to the clade comprising Teleogryllus and Velarifictorus. This study provides essential genetic information for genetic diversity analysis of T. eous.

Biological Adaptations Associated with Dehydration in Mosquitoes

Diseases that are transmitted by mosquitoes are a tremendous health and socioeconomic burden with hundreds of millions of people being impacted by mosquito-borne illnesses annually. Many factors have been implicated and extensively studied in disease transmission dynamics, but knowledge regarding how dehydration impacts mosquito physiology, behavior, and resulting mosquito-borne disease transmission remain underdeveloped. The lapse in understanding on how mosquitoes respond to dehydration stress likely obscures our ability to effectively study mosquito physiology, behavior, and vectorial capabilities. The goal of this review is to develop a profile of factors underlying mosquito biology that are altered by dehydration and the implications that are related to disease transmission.

Characterization of the complete mitogenome of Anopheles aquasalis, and phylogenetic divergences among Anopheles from diverse geographic zones

Whole mitogenome sequences (mtDNA) have been exploited for insect ecology studies, using them as molecular markers to reconstruct phylogenies, or to infer phylogeographic relationships and gene flow. Recent Anopheles phylogenomic studies have provided information regarding the time of deep lineage divergences within the genus. Here we report the complete 15,393 bp mtDNA sequences of Anopheles aquasalis, a Neotropical human malaria vector. When comparing its structure and base composition with other relevant and available anopheline mitogenomes, high similarity and conserved genomic features were observed. Furthermore, 22 mtDNA sequences comprising anopheline and Dipteran sibling species were analyzed to reconstruct phylogenies and estimate dates of divergence between taxa. Phylogenetic analysis using complete mtDNA sequences suggests that A. aquasalis diverged from the Anopheles albitarsis complex ~28 million years ago (MYA), and ~38 MYA from Anopheles darlingi. Bayesian analysis suggests that the most recent ancestor of Nyssorhynchus and Anopheles + Cellia was extant ~83 MYA, corroborating current estimates of ~79–100 MYA. Additional sampling and publication of African, Asian, and North American anopheline mitogenomes would improve the resolution of the Anopheles phylogeny and clarify early continental dispersal routes.

Is Genetic Continuity Between Anopheles sinensis (Diptera: Culicidae) and its Sibling Species Due to Gene Introgression or Incomplete Speciation?

The Anopheles mosquito Hyrcanus Group is widely distributed geographically across both Palearctic and Oriental regions and comprises 26 valid species. Although the species Anopheles sinensis Wiedemann (1828) is the most common in China and has a low potential vector rank, it has nevertheless long been thought to be an important natural malaria vector within the middle and lower reaches of the Yangtze River. A number of previous research studies have found evidence to support the occurrence of natural hybridization between An. sinensis and Anopheles kleini Rueda, 2005 (a competent malaria vector). We, therefore, collected a sample series of An. sinensis and morphologically similar species across China and undertook ribosomal and mitochondrial DNA analyses in order to assess genetic differentiation (Fst) and gene flow (Nm) amongst different groups. This enabled us to evaluate divergence times between morphologically similar species using the cytochrome oxidase I (COI) gene. The results of this study reveal significant genetic similarities between An. sinensis, An. kleini, and Anopheles belenrae Rueda, 2005 and therefore imply that correct molecular identifications will require additional molecular markers. As results also reveal the presence of gene flow between these three species, their taxonomic status will require further work. Data suggest that An. kleini is the most basal of the three species, while An. sinensis and An. belenrae share the closest genetic relationship.

Population genetics of the invasive Asian bush mosquito Aedes japonicus (Diptera, Culicidae) in Germany-a re-evaluation in a time period of separate populations merging

The Asian bush mosquito Aedes japonicus, endemic to East Asia, is one of the most expansive mosquito species in the world and has as yet established in 15 countries of Europe. Within Germany, the species has been spreading tremendously during the last years, and its four once geographically isolated populations were on the verge of merging in 2017. To reveal relationships and carry-over ways between the various populations, and thus, migration and displacement routes, the genetic make-up of Ae. japonicus from ten different locations throughout its German distribution area was investigated. For this purpose, a part of the mitochondrial DNA (nad4 gene) of collected specimens was sequenced and seven loci of short tandem repeats (microsatellites) were genotyped. When related to similar genetic studies carried out between 2012 and 2015, the results suggest that admixtures had since occurred, but no complete genetic mixture of populations had taken place. At the time of sampling for the present study, the western collection sites were still uniform in their genetic make-up; however, a carry-over of individuals from the southeastern to the northern and southwestern German populations was determined. Further introductions from abroad are possible. In summary, the genetic diversity of Ae. japonicus in Germany had grown considerably, thus increasing ecological variability and adaptability of the species. At this point (10 years after the first detection), it is not possible anymore to draw conclusions on the origins of the populations.

Evolution of sexually-transferred steroids and mating-induced phenotypes in Anopheles mosquitoes

Human malaria, which remains a major public health problem, is transmitted by a subset of Anopheles mosquitoes belonging to only three out of eight subgenera: Anopheles, Cellia and Nyssorhynchus. Unlike almost every other insect species, males of some Anopheles species produce steroid hormones which are transferred to females during copulation to influence their reproduction. Steroids are consequently a potential target for malaria vector control. Here, we analysed the evolution of sexually-transferred steroids and their effects on female reproductive traits across Anopheles by using a set of 16 mosquito species (five Anopheles, eight Cellia, and three Nyssorhynchus), including malaria vector and non-vector species. We show that male steroid production and transfer are specific to the Cellia and therefore represent a synapomorphy of this subgenus. Furthermore, we show that mating-induced effects in females are variable across species and differences are not correlated with sexually-transferred steroids or with Anopheles ability to transmit human malaria. Overall, our findings highlight that Anopheles mosquitoes have evolved different reproductive strategies, independently of being a malaria vector or not.

Mode and Rate of Evolution of Haemosporidian Mitochondrial Genomes: Timing the Radiation of Avian Parasites

Haemosporidians are a diverse group of vector-borne parasitic protozoa that includes the agents of human malaria; however, most of the described species are found in birds and reptiles. Although our understanding of these parasites’ diversity has expanded by analyses of their mitochondrial genes, there is limited information on these genes’ evolutionary rates. Here, 114 mitochondrial genomes (mtDNA) were studied from species belonging to four genera: Leucocytozoon, Haemoproteus, Hepatocystis, and Plasmodium. Contrary to previous assertions, the mtDNA is phylogenetically informative. The inferred phylogeny showed that, like the genus Plasmodium, the Leucocytozoon and Haemoproteus genera are not monophyletic groups. Although sensitive to the assumptions of the molecular dating method used, the estimated times indicate that the diversification of the avian haemosporidian subgenera/genera took place after the Cretaceous–Paleogene boundary following the radiation of modern birds. Furthermore, parasite clade differences in mtDNA substitution rates and strength of negative selection were detected. These differences may affect the biological interpretation of mtDNA gene lineages used as a proxy to species in ecological and parasitological investigations. Given that the mitochondria are critically important in the parasite life cycle stages that take place in the vector and that the transmission of parasites belonging to particular clades has been linked to specific insect families/subfamilies, this study suggests that differences in vectors have affected the mode of evolution of haemosporidian mtDNA genes. The observed patterns also suggest that the radiation of haemosporidian parasites may be the result of community-level evolutionary processes between their vertebrate and invertebrate hosts.

Supersensitive Odorant Receptor Underscores Pleiotropic Roles of Indoles in Mosquito Ecology

Mosquitoes exhibit highly diverse and fast evolving odorant receptors (ORs). The indole-sensitive OR gene clade, comprised of Or2 and Or10 is a notable exception on account of its conservation in both mosquito subfamilies. This group of paralogous genes exhibits a complex developmental expression pattern in Aedes aegypti: AaegOr2 is expressed in both adults and larvae, AaegOr10 is adult-specific and a third member named AaegOr9 is larva-specific. OR2 and OR10 have been deorphanized and are selectively activated by indole and skatole, respectively. Using the two-electrode voltage clamp of Xenopus oocytes expressing Ae. aegypti ORs, we show that AaegOR9 is supersensitive and narrowly tuned to skatole. Our findings suggest that Ae. aegypti has evolved two distinct molecular strategies to detect skatole in aquatic and terrestrial environments, highlighting the central ecological roles of indolic compounds in the evolutionary and life histories of these insects.

The PAGODAS protocol: pediatric assessment group of dengue and Aedes saliva protocol to investigate vector-borne determinants of Aedes-transmitted arboviral infections in Cambodia

BACKGROUND

Mosquito-borne arboviruses, like dengue virus, continue to cause significant global morbidity and mortality, particularly in Southeast Asia. When the infectious mosquitoes probe into human skin for a blood meal, they deposit saliva containing a myriad of pharmacologically active compounds, some of which alter the immune response and influence host receptivity to infection, and consequently, the establishment of the virus. Previous reports have highlighted the complexity of mosquito vector-derived factors and immunity in the success of infection. Cumulative evidence from animal models and limited data from humans have identified various vector-derived components, including salivary components, that are co-delivered with the pathogen and play an important role in the dissemination of infection. Much about the roles and effects of these vector-derived factors remain to be discovered.

METHODS/DESIGN

We describe a longitudinal, pagoda (community)-based pediatric cohort study to evaluate the burden of dengue virus infection and document the immune responses to salivary proteins of Aedes aegypti, the mosquito vector of dengue, Zika, and chikungunya viruses. The study includes community-based seroprevalence assessments in the peri-urban town of Chbar Mon in Kampong Speu Province, Cambodia. The study aims to recruit 771 children between the ages of 2 and 9 years for a three year period of longitudinal follow-up, including twice per year (rainy and dry season) serosurveillance for dengue seroconversion and Ae. aegypti salivary gland homogenate antibody intensity determinations by ELISA assays. Diagnostic tests for acute dengue, Zika and chikungunya viral infections will be performed by RT-PCR.

DISCUSSION

This study will serve as a foundation for further understanding of mosquito saliva immunity and its impact on Aedes-transmitted arboviral diseases endemic to Cambodia.

TRIAL REGISTRATION

NCT03534245 registered on 23 May 2018.

Complete mitogenome of Anopheles sinensis and mitochondrial insertion segments in the nuclear genomes of 19 mosquito species

Anopheles sinensis is a major malarial vector in China and Southeast Asia. The mitochondria is involved in many important biological functions. Nuclear mitochondrial DNA segments (NUMTs) are common in eukaryotic organisms, but their characteristics are poorly understood. We sequenced and analyzed the complete mitochondrial (mt) genome of An. sinensis. The mt genome is 15,418 bp long and contains 13 protein-coding genes (PCGs), two rRNAs, 22 tRNAs and a large non-coding region. Its gene arrangement is similar to previously published mosquito mt genomes. We identified and analyzed the NUMTs of 19 mosquito species with both nuclear genomes and mt genome sequences. The number, total length and density of NUMTs are significantly correlated with genome size. About half of NUMTs are short (< 200 bp), but larger genomes can house longer NUMTs. NUMTs may help explain genome size expansion in mosquitoes. The expansion due to mitochondrial insertion segments is variable in different insect groups. PCGs are transferred to nuclear genomes at a higher frequency in mosquitoes, but NUMT origination is more different than in mammals. Larger-sized nuclear genomes have longer mt genome sequences in both mosquitoes and mammals. The study provides a foundation for the functional research of mitochondrial genes in An. sinensis and helps us understand the characteristics and origin of NUMTs and the potential contribution to genome expansion.

Identification of Uranotaenia sapphirina as a specialist of annelids broadens known mosquito host use patterns

Feeding upon vertebrate blood by mosquitoes permits transmission of diverse pathogens, including viruses, protozoa, and nematodes. Despite over a century of intensive study, no mosquito species is known to specialize on non-vertebrate hosts. Using molecular analyses and field observations, we provide the first evidence, to our knowledge, that a mosquito, Uranotaenia sapphirina, specializes on annelid hosts (earthworms and leeches) while its sympatric congener, Uranotaenia lowii, feeds only on anurans (frogs and toads). Our results demonstrate that Ur. sapphirina feeds on annelid hosts (100% of identified blood meals; n = 72; collected throughout Florida), findings that are supported by field observations of these mosquitoes feeding on Sparganophilus worms and freshwater leeches. These findings indicate that adult mosquitoes utilize a much broader range of host taxa than previously recognized, with implications for epidemiology and the evolution of host use patterns in mosquitoes.

Complete Anopheles funestus mitogenomes reveal an ancient history of mitochondrial lineages and their distribution in southern and central Africa

Anopheles funestus s.s. is a primary vector of malaria in sub-Saharan Africa. Despite its important role in human Plasmodium transmission, evolutionary history, genetic diversity, and population structure of An. funestus in southern and central Africa remains understudied. We deep sequenced, assembled, and annotated the complete mitochondrial genome of An. funestus s.s. for the first time, providing a foundation for further genetic research of this important malaria vector species. We further analyzed the complete mitochondrial genomes of 43 An. funestus s.s. from three sites in Zambia, Democratic Republic of the Congo, and Tanzania. From these 43 mitogenomes we identified 41 unique haplotypes that comprised 567 polymorphic sites. Bayesian phylogenetic reconstruction confirmed the co-existence of two highly divergent An. funestus maternal lineages, herein defined as lineages I and II, in Zambia and Tanzania. The estimated coalescence time of these two mitochondrial lineages is ~500,000 years ago (95% HPD 426,000–594,000 years ago) with subsequent independent diversification. Haplotype network and phylogenetic analysis revealed two major clusters within lineage I, and genetic relatedness of samples with deep branching in lineage II. At this time, data suggest that the lineages are partially sympatric. This study illustrates that accurate retrieval of full mitogenomes of Anopheles vectors enables fine-resolution studies of intraspecies genetic relationships, population differentiation, and demographic history. Further investigations on whether An. funestus mitochondrial lineages represent biologically meaningful populations and their potential implications for malaria vector control are warranted.

Characterization of the adult Aedes aegypti early midgut peritrophic matrix proteome using LC-MS

The Aedes aegypti mosquito is the principal vector of arboviruses such as dengue, chikungunya, yellow fever, and Zika virus. These arboviruses are transmitted during adult female mosquito bloodfeeding. While these viruses must transverse the midgut to replicate, the blood meal must also reach the midgut to be digested, absorbed, or excreted, as aggregation of blood meal metabolites can be toxic to the female mosquito midgut. The midgut peritrophic matrix (PM), a semipermeable extracellular layer comprised of chitin fibrils, glycoproteins, and proteoglycans, is one such mechanism of protection for the mosquito midgut. However, this structure has not been characterized for adult female Ae. aegypti. We conducted a mass spectrometry based proteomic analysis to identify proteins that comprise or are associated with the adult female Ae. aegypti early midgut PM. Altogether, 474 unique proteins were identified, with 115 predicted as secreted. GO-term enrichment analysis revealed an abundance of serine-type proteases and several known and novel intestinal mucins. In addition, approximately 10% of the peptides identified corresponded to known salivary proteins, indicating Ae. aegypti mosquitoes extensively swallow their own salivary secretions. However, the physiological relevance of this remains unclear, and further studies are needed to determine PM proteins integral for midgut protection from blood meal derived toxicity and pathogen protection. Finally, we describe substantial discordance between previously described transcriptionally changes observed in the midgut in response to a bloodmeal and the presence of the corresponding protein in the PM. Data are available via ProteomeXchange with identifier PXD007627.

Phylogeny of Anophelinae using mitochondrial protein coding genes

Malaria is a vector-borne disease that is a great burden on the poorest and most marginalized communities of the tropical and subtropical world. Approximately 41 species of Anopheline mosquitoes can effectively spread species of Plasmodium parasites that cause human malaria. Proposing a natural classification for the subfamily Anophelinae has been a continuous effort, addressed using both morphology and DNA sequence data. The monophyly of the genus Anopheles, and phylogenetic placement of the genus Bironella, subgenera Kerteszia, Lophopodomyia and Stethomyia within the subfamily Anophelinae, remain in question. To understand the classification of Anophelinae, we inferred the phylogeny of all three genera (Anopheles, Bironella, Chagasia) and major subgenera by analysing the amino acid sequences of the 13 protein coding genes of 150 newly sequenced mitochondrial genomes of Anophelinae and 18 newly sequenced Culex species as outgroup taxa, supplemented with 23 mitogenomes from GenBank. Our analyses generally place genus Bironella within the genus Anopheles, which implies that the latter as it is currently defined is not monophyletic. With some inconsistencies, Bironella was placed within the major clade that includes Anopheles, Cellia, Kerteszia, Lophopodomyia, Nyssorhynchus and Stethomyia, which were found to be monophyletic groups within Anophelinae. Our findings provided robust evidence for elevating the monophyletic groupings Kerteszia, Lophopodomyia, Nyssorhynchus and Stethomyia to genus level; genus Anopheles to include subgenera Anopheles, Baimaia, Cellia and Christya; Anopheles parvus to be placed into a new genus; Nyssorhynchus to be elevated to genus level; the genus Nyssorhynchus to include subgenera Myzorhynchella and Nyssorhynchus; Anopheles atacamensis and Anopheles pictipennis to be transferred from subgenus Nyssorhynchus to subgenus Myzorhynchella; and subgenus Nyssorhynchus to encompass the remaining species of Argyritarsis and Albimanus Sections.

Phylogeny of Anophelinae using mitochondrial protein coding genes

Malaria is a vector-borne disease that is a great burden on the poorest and most marginalized communities of the tropical and subtropical world. Approximately 41 species of Anopheline mosquitoes can effectively spread species of Plasmodium parasites that cause human malaria. Proposing a natural classification for the subfamily Anophelinae has been a continuous effort, addressed using both morphology and DNA sequence data. The monophyly of the genus Anopheles, and phylogenetic placement of the genus Bironella, subgenera Kerteszia, Lophopodomyia and Stethomyia within the subfamily Anophelinae, remain in question. To understand the classification of Anophelinae, we inferred the phylogeny of all three genera (Anopheles, Bironella, Chagasia) and major subgenera by analysing the amino acid sequences of the 13 protein coding genes of 150 newly sequenced mitochondrial genomes of Anophelinae and 18 newly sequenced Culex species as outgroup taxa, supplemented with 23 mitogenomes from GenBank. Our analyses generally place genus Bironella within the genus Anopheles, which implies that the latter as it is currently defined is not monophyletic. With some inconsistencies, Bironella was placed within the major clade that includes Anopheles, Cellia, Kerteszia, Lophopodomyia, Nyssorhynchus and Stethomyia, which were found to be monophyletic groups within Anophelinae. Our findings provided robust evidence for elevating the monophyletic groupings Kerteszia, Lophopodomyia, Nyssorhynchus and Stethomyia to genus level; genus Anopheles to include subgenera Anopheles, Baimaia, Cellia and Christya; Anopheles parvus to be placed into a new genus; Nyssorhynchus to be elevated to genus level; the genus Nyssorhynchus to include subgenera Myzorhynchella and Nyssorhynchus; Anopheles atacamensis and Anopheles pictipennis to be transferred from subgenus Nyssorhynchus to subgenus Myzorhynchella; and subgenus Nyssorhynchus to encompass the remaining species of Argyritarsis and Albimanus Sections.

Sequencing and analysis of the complete mitochondrial genome in Anopheles sinensis (Diptera: Culicidae)

Background

Anopheles sinensis (Diptera: Culicidae) is a primary vector of Plasmodium vivax and Brugia malayi in most regions of China. In addition, its phylogenetic relationship with the cryptic species of the Hyrcanus Group is complex and remains unresolved. Mitochondrial genome sequences are widely used as molecular markers for phylogenetic studies of mosquito species complexes, of which mitochondrial genome data of An. sinensis is not available.

Methods

An. sinensis samples was collected from Shandong, China, and identified by molecular marker. Genomic DNA was extracted, followed by the Illumina sequencing. Two complete mitochondrial genomes were assembled and annotated using the mitochondrial genome of An. gambiae as reference. The mitochondrial genomes sequences of the 28 known Anopheles species were aligned and reconstructed phylogenetic tree by Maximum Likelihood (ML) method.

Findings

The length of complete mitochondrial genomes of An. sinensis was 15,076 bp and 15,138 bp, consisting of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and an AT-rich control region. As in other insects, most mitochondrial genes are encoded on the J strand, except for ND5, ND4, ND4L, ND1, two rRNA and eight tRNA genes, which are encoded on the N strand. The bootstrap value was set as 1000 in ML analyses. The topologies restored phylogenetic affinity within subfamily Anophelinae. The ML tree showed four major clades, corresponding to the subgenera Cellia, Anopheles, Nyssorhynchus and Kerteszia of the genus Anopheles.

Conclusions

The complete mitochondrial genomes of An. sinensis were obtained. The number, order and transcription direction of An. sinensis mitochondrial genes were the same as in other species of family Culicidae.

Electronic supplementary material

The online version of this article (10.1186/s40249-017-0362-7) contains supplementary material, which is available to authorized users.

Complete mitochondrial genomes of Anopheles stephensi and An. dirus and comparative evolutionary mitochondriomics of 50 mosquitoes

To better understand the phylogeny and evolution of mosquitoes, the complete mitochondrial genome (mitogenome) of Anopheles stephensi and An. dirus were sequenced and annotated, and a total of 50 mosquito mitogenomes were comparatively analyzed. The complete mitogenome of An. stephensi and An. dirus is 1,5371 bp and 1,5406 bp long, respectively. The main features of the 50 mosquito mitogenomes are conservative: 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes, positive AT-skew and negative GC-skew. The gene order trnA-trnR in ancestral insects is rearranged. All tRNA genes have the typical clover leaf secondary structure but tRNA Ser . The control regions are highly variable in size. PCGs show signals of purifying selection, but evidence for positive selection in ND2, ND4 and ND6 is found. Bayesian and Maximum Likelihood phylogenetic analyses based on all PCG nucleotides produce an identical tree topology and strongly support the monophyly of subgenera Cellia, Anopheles, Keterszia and Nyssorhynchus, the sister relationship of the subgenera Nyssorhynchus and Keterszia, and Cellia and Anopheles. The most recent ancestor of the genus Anopheles and Culicini + Aedini exited ~145 Mya ago. This is the first comprehensive study of mosquito mitogenomes, which are effective for mosquito phylogeny at various taxonomic levels.

Phylogenetic Relationships Among Malaria Vectors and Closely Related Species in Thailand Using Multilocus DNA Sequences

The evolutionary and taxonomic status is important for understanding speciation events and phylogenetic relationships between closely related vector and nonvector species. This information is useful for targeting important disease vector species groups for the development of novel genetic-based vector and pathogen control methods. In this study, different phylogenetic analyses were performed to reconstruct phylogenetic trees for the primary malaria vectors in Thailand based on sequence information of 4 DNA fragments from the nuclear and mitochondrial regions. The primary Anopheles species in the subgenus Cellia involved in malaria transmission in Thailand separate clearly into 3 distinct clades: the Leucosphyrus group, Minimus subgroup, and Maculatus group. The phylogenetic trees based on different reconstructed algorithms and different gene regions provided congruent phylogenetic status of the mosquito species studied. The phylogenetic relationships of malaria vector species examined followed similar patterns based on morphological characters. An estimate of the divergence time among the Anopheles species infers that they were present during the Eocene and Miocene periods (>41 million years ago). Congruent phylogenetic analysis of malaria vectors is presented with different algorithms and gene regions. The nuclear TOLL6 fragment appears useful for molecular phylogenetic, species DNA barcode, and Anopheles population genetic analyses.

Retention of duplicated long-wavelength opsins in mosquito lineages by positive selection and differential expression

Background

Opsins are light sensitive receptors associated with visual processes. Insects typically possess opsins that are stimulated by ultraviolet, short and long wavelength (LW) radiation. Six putative LW-sensitive opsins predicted in the yellow fever mosquito, Aedes aegypti and malaria mosquito, Anopheles gambiae, and eight in the southern house mosquito, Culex quinquefasciatus, suggest gene expansion in the Family Culicidae (mosquitoes) relative to other insects. Here we report the first detailed molecular and evolutionary analyses of LW opsins in three mosquito vectors, with a goal to understanding the molecular basis of opsin-mediated visual processes that could be exploited for mosquito control.

Results

Time of divergence estimates suggest that the mosquito LW opsins originated from 18 or 19 duplication events between 166.9/197.5 to 1.07/0.94 million years ago (MY) and that these likely occurred following the predicted divergence of the lineages Anophelinae and Culicinae 145–226 MY. Fitmodel analyses identified nine amino acid residues in the LW opsins that may be under positive selection. Of these, eight amino acids occur in the N and C termini and are shared among all three species, and one residue in TMIII was unique to culicine species. Alignment of 5′ non-coding regions revealed potential Conserved Non-coding Sequences (CNS) and transcription factor binding sites (TFBS) in seven pairs of LW opsin paralogs.

Conclusions

Our analyses suggest opsin gene duplication and residues possibly associated with spectral tuning of LW-sensitive photoreceptors. We explore two mechanisms - positive selection and differential expression mediated by regulatory units in CNS – that may have contributed to the retention of LW opsin genes in Culicinae and Anophelinae. We discuss the evolution of mosquito LW opsins in the context of major Earth events and possible adaptation of mosquitoes to LW-dominated photo environments, and implications for mosquito control strategies based on disrupting vision-mediated behaviors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0910-6) contains supplementary material, which is available to authorized users.

Genes and Odors Underlying the Recent Evolution of Mosquito Preference for Humans

Mosquito species that specialize in biting humans are few but dangerous. They include the African malaria vectors Anopheles gambiae and Anopheles coluzzii, as well as Aedes aegypti, the cosmopolitan vector of dengue, chikungunya, and yellow fever. These mosquitoes have evolved a remarkable innate preference for human odor that helps them find and bite us. Here I review what is known about this important evolutionary adaptation, from its historical documentation to its chemical and molecular basis.

The mitochondrial genome of the lepidopteran host cadaver (Thitarodes sp.) of Ophiocordyceps sinensis and related phylogenetic analysis

To understand the phylogeny of the host insect (Thitarodes sp.) of the fungus Ophiocordyceps sinensis, we sequenced, annotated and characterized the complete mitochondrial (mt) genome of the host cadaver of a natural O. sinensis. Further, we compared the Thitarodes sp. mt genome with those of the other 7 sequenced Hepialidae and examined the phylogenetic relationships using a constructed Maximum Likelihood (ML) phylogenetic tree and mt genomic features (genetic distances and intergenic spacers). The mt genome is a circular molecule of 16,280bp in length with a high A+T content (81.20%) and contains 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes and an AT-rich region. The gene arrangement is identical to the ancestral arrangement but differs from those of other lepidopteran mt genomes because of the arrangement of tRNA genes. The tRNA region, which is located between the AT-rich region and nad2, is trnI/trnQ/trnM (IQM) in Thitarodes sp., rather than the trnM/trnI/trnQ (MIQ) of the Lepidoptera-specific rearrangement. All PCGs begin with the canonical start codons ATN or NTG, except for cox1, which starts with CGA. Most PCGs terminate with the typical stop codon TAA, although some have an incomplete stop codon (T). The 1473bp AT-rich region is located between the rrnS (12S rRNA) and trnI, which is the longest sequenced in a Thitarodes mt genome to date, containing nine 112bp copies and one partial copy of a 55bp sequence. The results derived from the phylogenetic tree, the genetic distances and the intergenic spacers of the mt genome show that the host insect of O. sinensis belongs to the Thitarodes, while Endoclita signifer and Napialus hunanensis form a relatively distinct lineage from Thitarodes. The sequence and full annotation of this moth mt genome will provide more molecular information about the Exoporia within the Lepidoptera, and the clarification of its phylogeny will improve the management of this insect resource and the conservation and sustainable use of this endangered medicinal species in China.

The Mitochondrial Cytochrome Oxidase Subunit I Gene Occurs on a Minichromosome with Extensive Heteroplasmy in Two Species of Chewing Lice, Geomydoecus aurei and Thomomydoecus minor

In animals, mitochondrial DNA (mtDNA) typically occurs as a single circular chromosome with 13 protein-coding genes and 22 tRNA genes. The various species of lice examined previously, however, have shown mitochondrial genome rearrangements with a range of chromosome sizes and numbers. Our research demonstrates that the mitochondrial genomes of two species of chewing lice found on pocket gophers, Geomydoecus aurei and Thomomydoecus minor, are fragmented with the 1,536 base-pair (bp) cytochrome-oxidase subunit I (cox1) gene occurring as the only protein-coding gene on a 1,916-1,964 bp minicircular chromosome in the two species, respectively. The cox1 gene of T. minor begins with an atypical start codon, while that of G. aurei does not. Components of the non-protein coding sequence of G. aurei and T. minor include a tRNA (isoleucine) gene, inverted repeat sequences consistent with origins of replication, and an additional non-coding region that is smaller than the non-coding sequence of other lice with such fragmented mitochondrial genomes. Sequences of cox1 minichromosome clones for each species reveal extensive length and sequence heteroplasmy in both coding and noncoding regions. The highly variable non-gene regions of G. aurei and T. minor have little sequence similarity with one another except for a 19-bp region of phylogenetically conserved sequence with unknown function.

The mitochondrial genomes of Culex tritaeniorhynchus and Culex pipiens pallens (Diptera: Culicidae) and comparison analysis with two other Culex species

Background

Culex tritaeniorhynchus and Culex pipiens pallens are the major vectors of the Japanese encephalitis virus and Wuchereria bancrofti, the causative agent of filariasis. The knowledge of mitochondrial genomes has been widely useful for the studies on molecular evolution, phylogenetics and population genetics.

Methods

In this study, we sequenced and annotated the mitochondrial (mt) genomes of Cx. tritaeniorhynchus and Cx. p. pallens, and performed a comparative analysis including four known mt genomes of species of the subgenus Culex (Culex). The phylogenetic relationships of Cx. tritaeniorhynchus, Cx. p. pallens and four known Culex mt genome sequences were reconstructed by maximum likelihood based on concatenated protein-coding gene sequences.

Results

Culex tritaeniorhynchus and Cx. p. pallens mt genomes are 14,844 bp and 15,617 bp long, both consists of 13 PCGs, 22 tRNAs, 2 rRNAs and 1 CR (not sequenced for Cx. tritaeniorhynchus). The initiation and termination codons of PCGs are ATN and TAA, respectively, except for COI starting with TCG, and COI and COII terminated with T. tRNAs have the typical clover-leaf secondary structures except for trnS^(AGN) that is lacking the DHU stem. 16S rRNA and 12S rRNA secondary structures were drawn for the first time for mosquito mt genomes. The control region of Cx. p. pallens mt genome is 747 bp long and with four tandem repeat structures. Phylogenetic analyses demonstrated that the mt genome of Cx. tritaeniorhynchus was significantly separated from the remaining five mt genomes of Culex spp. Culex p. pipiens, Cx. p. pallens and Cx. p. quinquefasciatus formed a monophyletic clade with Cx. p. quinquefasciatus linked in the middle of the clade, and Cx. p. pallens should have the same taxonomic level as Culex p. pipiens and Cx. p. quinquefasciatus.

Conclusions

The mt genomes of Cx. tritaeniorhynchus and Cx. p. pallens share the same gene composition and order with those of two other Culex species. Culex p. pallens of the Pipiens complex should have the same taxonomic level as Culex p. pipiens and Cx. p. quinquefasciatus investigated. We enriched the Culex mt genome data and provided a reference basis for further Culex mt genome sequencing and analyses.

Electronic supplementary material

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

Mitochondrial Genomes of Anopheles (Kerteszia) (Diptera: Culicidae) From the Atlantic Forest, Brazil

Mitochondrial genome sequences are widely used as molecular markers for phylogenetic studies of mosquito species complexes, such as the Anopheles albitarsis complex. Except for a few studies that employed a limited number of nuclear or mitochondrial loci to address the genetic structure and species status of Anopheles cruzii, Anopheles bellator, and Anopheles homunculus, little is known about genetic markers that can be employed in studies focusing on Kerteszia species. The complete mitochondrial genomes of seven specimens of An. bellator, An. cruzii, An. homunculus, and Anopheles laneanus were sequenced using long-range polymerase chain reaction and Illumina sequencing. The mitochondrial genomes varied from 15,446 to 15,738 bp in length and contained 37 genes (13 protein-encoding genes, 2 rRNA genes [12S rRNA and 16S rRNA] and 22 tRNA genes), and the AT-rich control region, as all do other Anopheles mitochondrial genomes sequenced to date. Specimens from four populations of An. cruzii showed differences in codon composition.

Dosage Compensation in the African Malaria Mosquito Anopheles gambiae

Dosage compensation is the fundamental process by which gene expression from the male monosomic X chromosome and from the diploid set of autosomes is equalized. Various molecular mechanisms have evolved in different organisms to achieve this task. In Drosophila, genes on the male X chromosome are upregulated to the levels of expression from the two X chromosomes in females. To test whether a similar mechanism is operating in immature stages of Anopheles mosquitoes, we analyzed global gene expression in the Anopheles gambiae fourth instar larvae and pupae using high-coverage RNA-seq data. In pupae of both sexes, the median expression ratios of X-linked to autosomal genes (X:A) were close to 1.0, and within the ranges of expression ratios between the autosomal pairs, consistent with complete compensation. Gene-by-gene comparisons of expression in males and females revealed mild female bias, likely attributable to a deficit of male-biased X-linked genes. In larvae, male to female ratios of the X chromosome expression levels were more female biased than in pupae, suggesting that compensation may not be complete. No compensation mechanism appears to operate in male germline of early pupae. Confirmation of the existence of dosage compensation in A. gambiae lays the foundation for research into the components of dosage compensation machinery in this important vector species.

Molecular tools for studying the major malaria vector Anopheles funestus: improving the utility of the genome using a comparative poly(A) and Ribo-Zero RNAseq analysis

Background

Next-generation sequencing (NGS) offers great opportunities for studying the biology of insect vectors of disease. Prerequisites for successful analyses include high quality annotated genome assemblies and that techniques designed for use with model organisms be tested and optimised for use with these insects. We aimed to test and improve genomic tools for studying the major malaria vector Anopheles funestus.

Results

To guide future RNAseq transcriptomic studies of An. funestus, we compared two methods for enrichment of non-ribosomal RNA for analysis: enrichment of polyadenylated RNA and ribosomal RNA depletion using a kit designed to deplete human/rat/mouse rRNA. We found large differences between the two methods in the resulting transcriptomes, some of which is due to differential representation of polyadenylated and non-polyadenylated transcripts. We used the RNAseq data for validation and targeted manual editing of the draft An. funestus genome annotation, validating 62 % of annotated introns, manually improving the annotation of seven gene families involved in the detoxification of xenobiotics and integrated two published transcriptomic datasets with the recently published genome assembly.

Conclusions

The mRNA enrichment method makes a substantial, replicable difference to the transcriptome composition, at least partly due to the representation of non-polyadenylated transcripts in the final transcriptome. Therefore, great care should be taken in comparing gene expression data among studies. Ribosomal RNA depletion of total RNA using a kit designed to deplete human/rat/mouse rRNA works in mosquitoes and, we argue, results in a truer representation of the transcriptome than poly(A) selection. The An. funestus genome annotation can be considerably improved with the help of these new RNAseq data and further guided manual gene editing efforts will be of great benefit to the Anopheles research community for studies of this insect’s genome and transcriptome.

Electronic supplementary material

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

Mitochondrial genomes and comparative analyses of Culex camposi, Culex coronator, Culex usquatus and Culex usquatissimus (Diptera:Culicidae), members of the coronator group

BACKGROUND

The Coronator Group currently encompasses six morphologically similar species (Culex camposi Dyar, Culex coronator Dyar and Knab, Culex covagarciai Forattini, Culex usquatus Dyar, Culex usquatissimus Dyar, and Culex ousqua Dyar). Culex coronator has been incriminated as a potential vector of West Nile Virus (WNV), Saint Louis Encephalitis Virus (SLEV), and Venezuelan Equine Encephalitis Virus (VEEV). The complete mitochondrial genome of Cx. coronator, Cx. usquatus, Cx.usquatissimus, and Cx. camposi was sequenced, annotated, and analyzed to provide genetic information about these species.

RESULTS

The mitochondrial genomes of Cx. coronator, Cx. usquatus, Cx.usquatissimus, and Cx. camposi varied from 15,573 base pairs in Cx. usquatus to 15,576 in Cx. coronator. They contained 37 genes (13 protein-encoding genes, 2 rRNA genes, and 22 tRNA genes) and the AT-rich control region. Comparative analyses of the 37 genes demonstrated the mitochondrial genomes to be composed of variable and conserved genes. Despite the small size, the ATP8, ATP6 plus NADH5 protein-encoding genes were polymorphic, whereas tRNAs and rRNAs were conserved. The control region contained some poly-T stretch. The Bayesian phylogenetic tree corroborated that both the Coronator Group and the Culex pipens complex are monophyletic taxa.

CONCLUSIONS

The mitochondrial genomes of Cx. coronator, Cx. usquatus, Cx. usquatissimus and Cx. camposi share the same gene composition and arrangement features that match to those reported for most Culicidae species. They are composed of the same 37 genes and the AT-rich control region, which contains poly-T stretches that may be involved in the functional role of the mitochondrial genome. Taken together, results of the dN/dS ratios, the sliding window analyses and the Bayesian phylogenetic analyses suggest that ATP6, ATP8 and NADH5 are promising genes to be employed in phylogenetic studies involving species of the Coronator Group, and probably other species groups of the subgenus Culex. Bayesian topology corroborated the morphological hypothesis of the Coronator Group as monophyletic lineage within the subgenus Culex.

The Phylogeny and Evolutionary Timescale of Muscoidea (Diptera: Brachycera: Calyptratae) Inferred from Mitochondrial Genomes

Muscoidea is a significant dipteran clade that includes house flies (Family Muscidae), latrine flies (F. Fannidae), dung flies (F. Scathophagidae) and root maggot flies (F. Anthomyiidae). It is comprised of approximately 7000 described species. The monophyly of the Muscoidea and the precise relationships of muscoids to the closest superfamily the Oestroidea (blow flies, flesh flies etc) are both unresolved. Until now mitochondrial (mt) genomes were available for only two of the four muscoid families precluding a thorough test of phylogenetic relationships using this data source. Here we present the first two mt genomes for the families Fanniidae (Euryomma sp.) (family Fanniidae) and Anthomyiidae (Delia platura (Meigen, 1826)). We also conducted phylogenetic analyses containing of these newly sequenced mt genomes plus 15 other species representative of dipteran diversity to address the internal relationship of Muscoidea and its systematic position. Both maximum-likelihood and Bayesian analyses suggested that Muscoidea was not a monophyletic group with the relationship: (Fanniidae + Muscidae) + ((Anthomyiidae + Scathophagidae) + (Calliphoridae + Sarcophagidae)), supported by the majority of analysed datasets. This also infers that Oestroidea was paraphyletic in the majority of analyses. Divergence time estimation suggested that the earliest split within the Calyptratae, separating (Tachinidae + Oestridae) from the remaining families, occurred in the Early Eocene. The main divergence within the paraphyletic muscoidea grade was between Fanniidae + Muscidae and the lineage ((Anthomyiidae + Scathophagidae) + (Calliphoridae + Sarcophagidae)) which occurred in the Late Eocene.