Search

Background Phylogenetic analyses of closely related species of mosquitoes are important for better understanding the evolution of traits contributing to transmission of vector-borne diseases. Six out of 41 dominant malaria vectors of the genus Anopheles in the world belong to the Maculipennis Group, which is subdivided into two Nearctic subgroups (Freeborni and Quadrimaculatus) and one Palearctic (Maculipennis) subgroup. Although previous studies considered the Nearctic subgroups as ancestral, details about their relationship with the Palearctic subgroup, and their migration times and routes from North America to Eurasia remain controversial. The Palearctic species An. beklemishevi is currently included in the Nearctic Quadrimaculatus subgroup adding to the uncertainties in mosquito systematics. Results To reconstruct historic relationships in the Maculipennis Group, we conducted a phylogenomic analysis of 11 Palearctic and 2 Nearctic species based on sequences of 1271 orthologous genes. The analysis indicated that the Palearctic species An. beklemishevi clusters together with other Eurasian species and represents a basal lineage among them. Also, An. beklemishevi is related more closely to An. freeborni, which inhabits the Western United States, rather than to An. quadrimaculatus, a species from the Eastern United States. The time-calibrated tree suggests a migration of mosquitoes in the Maculipennis Group from North America to Eurasia about 20–25 million years ago through the Bering Land Bridge. A Hybridcheck analysis demonstrated highly significant signatures of introgression events between allopatric species An. labranchiae and An. beklemishevi. The analysis also identified ancestral introgression events between An. sacharovi and its Nearctic relative An. freeborni despite their current geographic isolation. The reconstructed phylogeny suggests that vector competence and the ability to enter complete diapause during winter evolved independently in different lineages of the Maculipennis Group. Conclusions Our phylogenomic analyses reveal migration routes and adaptive radiation timing of Holarctic malaria vectors and strongly support the inclusion of An. beklemishevi into the Maculipennis Subgroup. Detailed knowledge of the evolutionary history of the Maculipennis Subgroup provides a framework for examining the genomic changes related to ecological adaptation and susceptibility to human pathogens. These genomic variations may inform researchers about similar changes in the future providing insights into the patterns of disease transmission in Eurasia.

BackgroundUnderstanding the evolutionary relationships between closely related taxa is important for mosquitoes that transmit human diseases. Six out of 41 dominant malaria vectors in the world belong to the Maculipennis group, which is subdivided into two North American subgroups (Freeborni and Quadrimaculatus), and one Eurasian (Maculipennis) subgroup. Although previous studies considered the Nearctic subgroups as ancestral, details about their relationship with the Palearctic subgroup, and their migration times and routes from North America to Eurasia remain controversial. The Eurasian species An. beklemishevi is currently included in the North American Quadrimaculatus subgroup adding to the uncertainties in mosquito systematics.ResultsTo reconstruct historic relationships between the North American and Eurasian mosquitoes, we conducted a phylogenomic analysis of 11 Palearctic and 2 Nearctic species based on 1271 orthologous genes using their transcriptomic or genomic sequences. The analysis indicated that the Palearctic species An. beklemishevi clusters together with other Eurasian species and represents a basal lineage among them. Also, An. beklemishevi is related more closely to An. freeborni, which inhabits the Western United States, rather than to An. quadrimaculatus, a species from the Eastern United States. The time-calibrated tree suggests a migration of mosquitoes in the Maculipennis group from North America to Eurasia about 20-25 million years ago through the Bering Land Bridge. A Hybridcheck analysis demonstrated highly significant signatures of introgression events between allopatric species An. labranchiae and An. beklemishevi. The analysis also identified ancestral introgression events between An. sacharovi and its Nearctic relative An. freeborni despite their current geographic isolation.ConclusionsOur phylogenomic analyses reveal migration routes and adaptive radiation timing of Holarctic malaria vectors and strongly support inclusion of An. beklemishevi into the Maculipennis subgroup. The vectorial capacity and the ability to diapause during winter evolved multiple times in Maculipennis evolution. Detailed knowledge of the evolutionary history in the Maculipennis subgroup will help us better understand the current and future patterns of malaria transmission in Eurasia.

Significance The presence of a dominant male-determining locus (M-locus) in one of a pair of autosomes establishes the male sex in the dengue fever mosquito Aedes aegypti. The Ae. aegypti M-locus contains 30 genes, including Nix, a previously reported male-determining factor. Here we show that the Nix transgene alone was sufficient to convert females into fertile males, which continued to produce sex-converted progeny. We also show that a second M-locus gene named myo-sex was needed for male flight. Nix-mediated sex conversion was 100% penetrant, heritable, and stable, indicating great potential for developing mosquito-control strategies to reduce vector populations by female-to-male conversion. This work also sheds lights into the molecular basis of the function of the M-locus., A dominant male-determining locus (M-locus) establishes the male sex (M/m) in the yellow fever mosquito, Aedes aegypti. Nix, a gene in the M-locus, was shown to be a male-determining factor (M factor) as somatic knockout of Nix led to feminized males (M/m) while transient expression of Nix resulted in partially masculinized females (m/m), with male reproductive organs but retained female antennae. It was not clear whether any of the other 29 genes in the 1.3-Mb M-locus are also needed for complete sex-conversion. Here, we report the generation of multiple transgenic lines that express Nix under the control of its own promoter. Genetic and molecular analyses of these lines provided insights unattainable from previous transient experiments. We show that the Nix transgene alone, in the absence of the M-locus, was sufficient to convert females into males with all male-specific sexually dimorphic features and male-like gene expression. The converted m/m males are flightless, unable to perform the nuptial flight required for mating. However, they were able to father sex-converted progeny when presented with cold-anesthetized wild-type females. We show that myo-sex, a myosin heavy-chain gene also in the M-locus, was required for male flight as knockout of myo-sex rendered wild-type males flightless. We also show that Nix-mediated female-to-male conversion was 100% penetrant and stable over many generations. Therefore, Nix has great potential for developing mosquito control strategies to reduce vector populations by female-to-male sex conversion, or to aid in a sterile insect technique that requires releasing only non-biting males.

Chromosomal inversions are important drivers of genome evolution. The Eurasian malaria vector Anopheles messeae has five polymorphic inversions. A cryptic species, An. daciae, has been discriminated from An. messeae based on five fixed nucleotide substitutions in the internal transcribed spacer 2 (ITS2) of ribosomal DNA. However, the inversion polymorphism in An. daciae and the genome divergence between these species remain unexplored. In this study, we sequenced the ITS2 region and analyzed the inversion frequencies of 289 Anopheles larvae specimens collected from three locations in the Moscow region. Five individual genomes for each of the two species were sequenced. We determined that An. messeae and An. daciae differ from each other by the frequency of polymorphic inversions. Inversion X1 was fixed in An. messeae but polymorphic in An. daciae populations. The genome sequence comparison demonstrated genome-wide divergence between the species, especially pronounced on the inversion-rich X chromosome (mean Fst = 0.331). The frequency of polymorphic autosomal inversions was higher in An. messeae than in An. daciae. We conclude that the X chromosome inversions play an important role in the genomic differentiation between the species. Our study determined that An. messeae and An. daciae are closely related species with incomplete reproductive isolation. Published version

Chromosomal inversions are important drivers of genome evolution. The Eurasian malaria vector Anopheles messeae has five polymorphic inversions. A cryptic species, An. daciae, has been discriminated from An. messeae based on five fixed nucleotide substitutions in the internal transcribed spacer 2 (ITS2) of ribosomal DNA. However, the inversion polymorphism in An. daciae and the genome divergence between these species remain unexplored. In this study, we sequenced the ITS2 region and analyzed the inversion frequencies of 289 Anopheles larvae specimens collected from three locations in the Moscow region. Five individual genomes for each of the two species were sequenced. We determined that An. messeae and An. daciae differ from each other by the frequency of polymorphic inversions. Inversion X1 was fixed in An. messeae but polymorphic in An. daciae populations. The genome sequence comparison demonstrated genome-wide divergence between the species, especially pronounced on the inversion-rich X chromosome (mean Fst = 0.331). The frequency of polymorphic autosomal inversions was higher in An. messeae than in An. daciae. We conclude that the X chromosome inversions play an important role in the genomic differentiation between the species. Our study determined that An. messeae and An. daciae are closely related species with incomplete reproductive isolation.

Background A dominant malaria vector, Anopheles messeae, is a highly polymorphic species with wide distribution throughout Eurasia. Five highly polymorphic inversions associated with the geographical distribution of the species have been reported. A sister species, An. daciae, was described and discriminated from An. messeae based on five fixed nucleotide substitutions in the internal transcribed spacer 2 (ITS2) of ribosomal DNA. However, the levels of genomic divergence, chromosomal variation, and ecological differentiation between these two cryptic species remain unexplored. Results In this study, we sequenced ITS2 and analyzed the inversion frequencies of 289 Anopheles larvae specimens collected from three locations in the Moscow region. We identified a high abundance of both An. messeae and An. daciae in all three locations. Five individual genomes for each species of An. messeae and An. daciae from one location were sequenced. Our study confirmed five previously described nucleotide substitutions in the ITS2 of An. messeae. However, we found that the ITS2 sequence in An. daciae was heterogenic in three of the five positions. Fixed nucleotide differences between An. messeae and An. daciae were found only in the last two positions. One mosquito was identified as a hybrid between An. messeae and An. daciae based on heterogeneous substitutions in all five positions. Although, the genomic sequence comparison demonstrated genome-wide divergence between the two species, which is especially pronounced on the X chromosome, an ADMIXTURE cluster analysis demonstrated the presence of two admixed individuals suggesting ongoing hybridization. Cytogenetic analysis demonstrated that An. messeae and An. daciae significantly differ from each other by their frequency of polymorphic inversions. Inversion X1 was fixed in An. messeae but was polymorphic in all An. daciae populations. The frequency of polymorphic autosomal inversions was higher in An. messeae than in An. daciae. The species composition was different among the studied locations suggesting species-specific ecological preferences. Conclusions Our study demonstrated that An. messeae and An. daciae represent closely related cryptic species with incomplete reproductive isolation that are able to maintain genomic differentiation in sympatry despite ongoing genetic introgression. The X chromosome plays an important role in the reproductive isolation between the species.

Anopheles atroparvus (Diptera: Culicidae) is one of the main malaria vectors of the Maculipennis group in Europe. Cytogenetic analysis based on salivary gland chromosomes has been used in taxonomic and population genetic studies of mosquitoes from this group. However, a high-resolution cytogenetic map that could be used in physical genome mapping in An. atroparvus is still lacking. In the present study, a high-quality photomap of the polytene chromosomes from ovarian nurse cells of An. atroparvus was developed. Using fluorescent in situ hybridization, 10 genes from the five largest genomic supercontigs on the polytene chromosome were localized and 28% of the genome was anchored to the cytogenetic map. The study established chromosome arm homology between An. atroparvus and the major African malaria vector Anopheles gambiae, suggesting a whole-arm translocation between autosomes of these two species. The standard photomap constructed for ovarian nurse cell chromosomes of An. atroparvus will be useful for routine physical mapping. This map will assist in the development of a fine-scale chromosome-based genome assembly for this species and will also facilitate comparative and evolutionary genomics studies in the genus Anopheles.

Background Malaria mosquitoes have had a remarkable stability in the number of chromosomes in their karyotype (2n = 6) during 100 million years of evolution. Moreover, autosomal arms were assumed to maintain their integrity even if their associations with each other changed via whole-arm translocations. Here we use high-coverage comparative physical genome mapping of three Anopheles species to test the extent of evolutionary conservation of chromosomal arms in malaria mosquitoes. Results In this study, we developed a physical genome map for Anopheles atroparvus, one of the dominant malaria vectors in Europe. Using fluorescence in situ hybridization (FISH) of DNA probes with the ovarian nurse cell polytene chromosomes and synteny comparison, we anchored 56 genomic scaffolds to the An. atroparvus chromosomes. The obtained physical map represents 89.6% of the An. atroparvus genome. This genome has the second highest mapping coverage among Anophelinae assemblies after An. albimanus, which has 98.2% of the genome assigned to its chromosomes. A comparison of the An. atroparvus, An. albimanus, and An. gambiae genomes identified partial-arm translocations between the autosomal arms that break down the integrity of chromosome elements in evolution affecting the structure of the genetic material in the pericentromeric regions. Unlike An. atroparvus and An. albimanus, all chromosome elements of An. gambiae are fully syntenic with chromosome elements of the putative ancestral Anopheles karyotype. We also detected nonrandom distribution of large conserved synteny blocks and confirmed a higher rate of inversion fixation in the X chromosome compared with autosomes. Conclusions Our study demonstrates the power of physical mapping for understanding the genome evolution in malaria mosquitoes. The results indicate that syntenic relationships among chromosome elements of Anopheles species have not been fully preserved because of multiple partial-arm translocations. Electronic supplementary material The online version of this article (10.1186/s12864-018-4663-4) contains supplementary material, which is available to authorized users.

There is an error in the legend for Fig 3, “A landmark-guided two-step physical mapping approach on Cx. quinquefasciatus chromosomes.” The complete, correct Fig 3 legend is given below. Fig 3 A landmark-guided two-step physical mapping approach on Cx. quinquefasciatus chromosomes. There is an error in the legend for Fig 4, “Chromosome idiograms with positions of supercontigs and genetic markers.” The complete, correct Fig 4 legend is given below. Fig 4 Chromosome idiograms with positions of supercontigs and genetic markers.

The genome assembly of southern house mosquito Cx. quinquefasciatus is represented by a high number of supercontigs with no order or orientation on the chromosomes. Although cytogenetic maps for the polytene chromosomes of this mosquito have been developed, their utilization for the genome mapping remains difficult because of the low number of high quality spreads in chromosome preparations. Therefore, a simple and robust mitotic-chromosome-based approach for the genome mapping of Cx. quinquefasciatus still needs to be developed. In this study, we performed physical mapping of 37 genomic supercontigs using fluorescent in situ hybridization on mitotic chromosomes from imaginal discs of 4th instar larvae. The genetic linkage map nomenclature was adopted for the chromosome numbering based on the direct positioning of 58 markers that were previously genetically mapped. The smallest, largest, and intermediate chromosomes were numbered as 1, 2, and 3, respectively. For idiogram development, we analyzed and described in detail the morphology and proportions of the mitotic chromosomes. Chromosomes were subdivided into 19 divisions and 72 bands of four different intensities. These idiograms were used for mapping the genomic supercontigs/genetic markers. We also determined the presence of length polymorphism in the q arm of sex-determining chromosome 1 in Cx. quinquefasciatus related to the size of ribosomal locus. Our physical mapping and previous genetic linkage mapping resulted in the chromosomal assignment of 13% of the total genome assembly to the chromosome bands.We provided the first detailed description, nomenclature, and idiograms for the mitotic chromosomes of Cx. quinquefasciatus. Further application of the approach developed in this study will help to improve the quality of the southern house mosquito genome. Published version