Two Nested Inversions in the X Chromosome Differentiate the Dominant Malaria Vectors in Europe, Anopheles atroparvus and Anopheles messeae.

The Maculipennis subgroup of malaria mosquitoes includes both dominant malaria vectors and non-vectors in Eurasia. Understanding the genetic factors, particularly chromosomal inversions, that differentiate Anopheles species can provide valuable insights for vector control strategies. Although autosomal inversions between the species in this subgroup have been characterized based on the chromosomal banding patterns, the number and positions of rearrangements in the X chromosome remain unclear due to the divergent banding patterns. Here, we identified two large X chromosomal inversions, approximately 13 Mb and 10 Mb in size, using fluorescence in situ hybridization. The inversion breakpoint regions were mapped by hybridizing 53 gene markers with polytene chromosomes of An. messeae. The DNA probes were designed based on gene sequences from the annotated An. atroparvus genome. The two nested inversions resulted in five syntenic blocks. Only two small syntenic blocks, which encompass 181 annotated genes in the An. atroparvus genome, changed their position and orientation in the X chromosome. The analysis of the An. atroparvus genome revealed an enrichment of gene ontology terms associated with immune system and mating behavior in the rearranged syntenic blocks. Additionally, the enrichment of DNA transposons was found in sequences homologous to three of the four breakpoint regions. This study demonstrates the successful application of the physical genome mapping approach to identify rearrangements that differentiate species in insects with polytene chromosomes.

Patterns of genetic differentiation imply distinct phylogeographic history of the mosquito species Anopheles messeae and Anopheles daciae in Eurasia.

Detailed knowledge of phylogeography is important for control of mosquito species involved in the transmission of human infectious diseases. Anopheles messeae is a geographically widespread and genetically diverse dominant vector of malaria in Eurasia. A closely related species, An. daciae, was originally distinguished from An. messeae based on five nucleotide substitutions in its ribosomal DNA (rDNA). However, the patterns of phylogeographic history of these species in Eurasia remain poorly understood. Here, using internal transcribed spacer 2 (ITS2) of rDNA and karyotyping for the species identification we determined the composition of five Anopheles species in 28 locations in Eurasia. Based on the frequencies of 11 polymorphic chromosomal inversions used as genetic markers, a large-scale population genetics analysis was performed of 1932 mosquitoes identified as An. messeae, An. daciae and their hybrids. The largest genetic differences between the species were detected in the X sex chromosome suggesting a potential involvement of this chromosome in speciation. The frequencies of autosomal inversions in the same locations differed by 13%-45% between the species demonstrating a restricted gene flow between the species. Overall, An. messeae was identified as a diverse species with a more complex population structure than An. daciae. The clinal gradients in frequencies of chromosomal inversions were determined in both species implicating their possible involvement in climate adaptations. The frequencies of hybrids were low ~1% in northern Europe but high up to 50% in south-eastern populations. Thus, our study revealed critical differences in patterns of phylogeographic history between An. messeae and An. daciae in Eurasia. This knowledge will help to predict the potential of the malaria transmission in the northern territories of the continent.

Phylogenomics revealed migration routes and adaptive radiation timing of Holarctic malaria mosquito species of the Maculipennis Group.

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.

Karyotypes and Physical Mapping of Ribosomal DNA with Oligo-Probes in Eranthis sect. Eranthis (Ranunculaceae).

A comparative karyotype analysis of four species of yellow-flowered Eranthis sect. Eranthis, i.e., E. bulgarica, E. cilicica, E. hyemalis, and E. longistipitata from different areas, has been carried out for the first time. All the studied specimens had somatic chromosome number 2n = 16 with basic chromosome number x = 8. Karyotypes of the investigated plants included five pairs of metacentric chromosomes and three pairs of submetacentric/subtelocentric chromosomes. The chromosome sets of the investigated species differ mainly in the ratio of submetacentric/subtelocentric chromosomes, their relative lengths, and arm ratios. A new oligonucleotide probe was developed and tested to detect 45S rDNA clusters. Using this probe and an oligonucleotide probe to 5S rDNA, 45S and 5S rDNA clusters were localized for the first time on chromosomes of E. cilicica, E. hyemalis, and E. longistipitata. Major 45S rDNA clusters were identified on satellite chromosomes in all the species; in E. cilicica, minor clusters were also identified in the terminal regions of one metacentric chromosome pair. The number and distribution of 5S rDNA clusters is more specific. In E. cilicica, two major clusters were identified in the pericentromeric region of a pair of metacentric chromosomes. Two major clusters in the pericentromeric region of a pair of submetacentric chromosomes and two major clusters in the interstitial region of a pair of metacentric chromosomes were observed in E. longistipitata. E. hyemalis has many clusters of different sizes, localized mainly in the pericentromeric regions. Summarizing new data on the karyotype structure of E. sect. Eranthis and previously obtained data on E. sect. Shibateranthis allowed conclusions to be formed about the clear interspecific karyological differences of the genus Eranthis.

An updated list of mosquito species in Armenia and Transcaucasian region responsible for Dirofilaria transmission: A review.

Canine dirofilariasis is a world-wide distributed mosquito-borne helminthiasis that has a potential zoonotic impact. This disease is a serious problem in Southern Caucasus region, including Armenia. Average extensity of Dirofilaria infection in canids in Armenia is 26.9%. At present, 77 mosquito species (Diptera: Culicidae) are presumed to have a role in the transmission of dirofilariasis. 80 species of Culicidae family mosquitoes are registered in Southern Caucasus Region, and 31 of them are described in Armenia. The following species of Culicidae can serve as vectors of Dirofilaria spp. in Armenia: Aedes vexans, Ae. caspius, Ae. albopictus, Ae. geniculatus, Culex pipiens, Cx. theileri, Cx. modestus, Anopheles maculipennis s.s., An. claviger, An. hyrcanus, An. superpictus, Coquillettidia richardi, Culiseta annulata and Uranotaenia unguiculta. We see a risk of transmission in active foci of dirofilariasis, as well as for human infection in Armenia. Comprehensive research should be conducted on detection of microfilariae in the target mosquito species and in target localities.

New Cytogenetic Photomap and Molecular Diagnostics for the Cryptic Species of the Malaria Mosquitoes Anopheles messeae and Anopheles daciae from Eurasia.

The Eurasian malaria vector Anopheles messeae is a widely spread and genetically diverse species. Five widespread polymorphic chromosomal inversions were found in natural populations of this mosquito. A cryptic species, Anopheles daciae, was differentiated from An. messeae by the presence of several nucleotide substitutions in the Internal Transcribed Spacer 2 (ITS2) region of ribosomal DNA. However, because of the absence of a high-quality reference cytogenetic map, the inversion polymorphisms in An. daciae and An. messeae remain poorly understood. Moreover, a recently determined heterogeneity in ITS2 in An. daciae questioned the accuracy of the previously used Restriction Fragment Length Polymorphism (RFLP) assay for species diagnostics. In this study, a standard-universal cytogenetic map was constructed based on orcein stained images of chromosomes from salivary glands for population studies of the chromosomal inversions that can be used for both An. messeae and An. daciae. In addition, a new ITS2-RFLP approach for species diagnostics was developed. Both methods were applied to characterize inversion polymorphism in populations of An. messeae and An. daciae from a single location in Western Siberia in Russia. The analysis demonstrates that cryptic species are remarkably different in their frequencies of chromosomal inversion variants. Our study supports previous observations that An. messeae has higher inversion polymorphism in all autosomes than the cryptic species An. daciae.

Physical Genome Mapping Using Fluorescence In Situ Hybridization with Mosquito Chromosomes.

The development of genomic resources and tools is an important step in designing novel approaches to genetic control of mosquitoes. Physical genome maps enhance the quality of the genome assemblies, improve gene annotation, and provide a better framework for comparative and population genomics studies in mosquitoes. In this chapter, we describe protocols for an important procedure in physical genome mapping-fluorescence in situ hybridization (FISH). We provide details on (1) dissection of salivary glands, ovaries, and imaginal discs for obtaining high-quality polytene or mitotic chromosome preparations; (2) DNA-labeling procedures and extraction of repetitive DNA fractions; and (3) approaches to FISH on polytene and mitotic chromosomes.

The Development of Cytogenetic Maps for Malaria Mosquitoes.

Anopheline mosquitoes are important vectors of human malaria. Next-generation sequencing opens new opportunities for studies of mosquito genomes to uncover the genetic basis of a Plasmodium transmission. Physical mapping of genome sequences to polytene chromosomes significantly improves reference assemblies. High-resolution cytogenetic maps are essential for anchoring genome sequences to chromosomes as well as for studying breakpoints of chromosome rearrangements and chromatin protein localization. Here we describe a detailed pipeline for the development of high-resolution cytogenetic maps using polytene chromosomes of malaria mosquitoes. We apply this workflow to the refinement of the cytogenetic map developed for Anopheles beklemishevi.

A standard photomap of the ovarian nurse cell chromosomes for the dominant malaria vector in Europe and Middle East Anopheles sacharovi.

BACKGROUND: Anopheles sacharovi is a dominant malaria vector species in South Europe and the Middle East which has a highly plastic behaviour at both adult and larval stages. Such plasticity has prevented this species from eradication by several anti-vector campaigns. The development of new genome-based strategies for vector control will benefit from genome sequencing and physical chromosome mapping of this mosquito. Although a cytogenetic photomap for chromosomes from salivary glands of An. sacharovi has been developed, no cytogenetic map suitable for physical genome mapping is available. METHODS: Mosquitoes for this study were collected at adult stage in animal shelters in Armenia. Polytene chromosome preparations were prepared from ovarian nurse cells. Fluorescent in situ hybridization (FISH) was performed using PCR amplified probes. RESULTS: This study constructed a high-quality standard photomap for polytene chromosomes from ovarian nurse cells of An. sacharovi. Following the previous nomenclature, chromosomes were sub-divided into 39 numbered and 119 lettered sub-divisions. Chromosomal landmarks for the chromosome recognition were described. Using FISH, 4 PCR-amplified genic probes were mapped to the chromosomes. The positions of the probes demonstrated gene order reshuffling between An. sacharovi and Anopheles atroparvus which has not been seen cytologically. In addition, this study described specific chromosomal landmarks that can be used for the cytotaxonomic diagnostics of An. sacharovi based on the banding pattern of its polytene chromosomes. CONCLUSIONS: This study constructed a high-quality standard photomap for ovarian nurse cell chromosomes of An. sacharovi and validated its utility for physical genome mapping. Based on the map, cytotaxonomic features for identification of An. sacharovi have been described. The cytogenetic map constructed in this study will assist in creating a chromosome-based genome assembly for this mosquito and in developing cytotaxonomic tools for identification of other species from the Maculipennis group.

Partial-arm translocations in evolution of malaria mosquitoes revealed by high-coverage physical mapping of the Anopheles atroparvus genome.

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.

A standard photomap of ovarian nurse cell chromosomes and inversion polymorphism in Anopheles beklemishevi.

BACKGROUND: Anopheles beklemishevi is a member of the Maculipennis group of malaria mosquitoes that has the most northern distribution among other members of the group. Although a cytogenetic map for the larval salivary gland chromosomes of this species has been developed, a high-quality standard cytogenetic photomap that enables genomics and population genetics studies of this mosquito at the adult stage is still lacking. METHODS: In this study, a cytogenetic map for the polytene chromosomes of An. beklemishevi from ovarian nurse cells was developed using high-resolution digital imaging from field collected mosquitoes. PCR-amplified DNA probes for fluorescence in situ hybridization (FISH) were designed based on the genome of An. atroparvus. The DNA probe obtained by microdissection procedures from the breakpoint region was labelled in a DOP-PCR reaction. Population analysis was performed on 371 specimens collected in 18 locations. RESULTS: We report the development of a high-quality standard photomap for the polytene chromosomes from ovarian nurse cells of An. beklemishevi. To confirm the suitability of the map for physical mapping, several PCR-amplified probes were mapped to the chromosomes of An. beklemishevi using FISH. In addition, we identified and mapped DNA probes to flanking regions of the breakpoints of two inversions on chromosome X of this species. Inversion polymorphism was determined in 13 geographically distant populations of An. beklemishevi. Four polymorphic inversions were detected. The positions of common chromosomal inversions were indicated on the map. CONCLUSIONS: The study constructed a standard photomap for ovarian nurse cell chromosomes of An. beklemishevi and tested its suitability for physical genome mapping and population studies. Cytogenetic analysis determined inversion polymorphism in natural populations of An. beklemishevi related to this species' adaptation.

Tissue-specific features of the X chromosome and nucleolus spatial dynamics in a malaria mosquito, Anopheles atroparvus.

Spatial organization of chromosome territories is important for maintenance of genomic stability and regulation of gene expression. Recent studies have shown tissue-specific features of chromosome attachments to the nuclear envelope in various organisms including malaria mosquitoes. However, other spatial characteristics of nucleus organization, like volume and shape of chromosome territories, have not been studied in Anopheles. We conducted a thorough analysis of tissue-specific features of the X chromosome and nucleolus volume and shape in follicular epithelium and nurse cells of the Anopheles atroparvus ovaries using a modern open-source software. DNA of the polytene X chromosome from ovarian nurse cells was obtained by microdissection and was used as a template for amplification with degenerate oligo primers. A fluorescently labeled X chromosome painting probe was hybridized with formaldehyde-fixed ovaries of mosquitoes using a 3D-FISH method. The nucleolus was stained by immunostaining with an anti-fibrillarin antibody. The analysis was conducted with TANGO-a software for a chromosome spatial organization analysis. We show that the volume and position of the X chromosome have tissue-specific characteristics. Unlike nurse cell nuclei, the growth of follicular epithelium nuclei is not accompanied with the proportional growth of the X chromosome. However, the shape of the X chromosome does not differ between the tissues. The dynamics of the X chromosome attachment regions location is tissue-specific and it is correlated with the process of nucleus growth in follicular epithelium and nurse cells.

The Physical Genome Mapping of Anopheles albimanus Corrected Scaffold Misassemblies and Identified Interarm Rearrangements in Genus Anopheles.

The genome of the Neotropical malaria vector Anopheles albimanus was sequenced as part of the 16 Anopheles Genomes Project published in 2015. The draft assembly of this species consisted of 204 scaffolds with an N50 scaffold size of 18.1 Mb and a total assembly size of 170.5 Mb. It was among the smallest genomes with the longest scaffolds in the 16 Anopheles species cluster, making An. albimanus the logical choice for anchoring the genome assembly to chromosomes. In this study, we developed a high-resolution cytogenetic photomap with completely straightened polytene chromosomes from the salivary glands of the mosquito larvae. Based on this photomap, we constructed a chromosome-based genome assembly using fluorescent in situ hybridization of PCR-amplified DNA probes. Our physical mapping, assisted by an ortholog-based bioinformatics approach, identified and corrected nine misassemblies in five large genomic scaffolds. Misassemblies mostly occurred in junctions between contigs. Our comparative analysis of scaffolds with the An. gambiae genome detected multiple genetic exchanges between pericentromeric regions of chromosomal arms caused by partial-arm translocations. The final map consists of 40 ordered genomic scaffolds and corrected fragments of misassembled scaffolds. The An. albimanus physical map comprises 98.2% of the total genome assembly and represents the most complete genome map among mosquito species. This study demonstrates that physical mapping is a powerful tool for correcting errors in draft genome assemblies and for creating chromosome-anchored reference genomes.

Chromosome evolution in malaria mosquitoes inferred from physically mapped genome assemblies.

Polymorphic inversions in mosquitoes are distributed nonrandomly among chromosomes and are associated with ecological, behavioral, and physiological adaptations related to pathogen transmission. Despite their significance, the patterns and mechanism of genome rearrangements are not well understood. Recent sequencing and physical mapping of the genomes for 16 Anopheles mosquito species provided an opportunity to study chromosome evolution at the highest resolution. New studies revealed that fixed rearrangement accumulated [Formula: see text]3 times faster on the X chromosome than on autosomes. The highest densities of transposable elements (TEs) and satellites of different sizes have also been found on the X chromosome, suggesting a mechanism for the inversion generation. The high rate of X chromosome rearrangements is in sharp contrast with the paucity of polymorphic inversions on the X in the majority of anopheline species. This paper highlights the advances in understanding chromosome evolution in malaria vectors and discusses possible future directions in studying mechanisms and biological roles of genome rearrangements.