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.

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

Analysis of the Metaphase Chromosome Karyotypes in Imaginal Discs of Aedes communis, Ae. punctor, Ae. intrudens, and Ae. rossicus (Diptera: Culicidae) Mosquitoes.

In this study, cytogenetic analysis of the metaphase chromosomes from imaginal discs of Aedes (Diptera: Culicidae) mosquitoes-Aedes communis, Ae. punctor, Ae. intrudens, and Ae. rossicus-was performed. The patterns of C-banding and DAPI staining of the heteroсhromatin and the length of the chromosomes demonstrate species specificity. In particular, the Ae. punctor chromosomes are the shortest compared with Ae. communis, Ae. intrudens, and Ae. rossicus, and they also carry additional C and DAPI bands in intercalary regions. The Ae. intrudens chromosomes are the longest, they have pericentromeric C bands, and they almost lack any DAPI bands near the centromere of chromosome 3 versus Ae. communis, which has the largest pericentromeric DAPI blocks in all three chromosome pairs. Ae. rossicus also possesses DAPI bands in the centromeric regions of all chromosomes, but their staining is weaker compared with those of Ae. communis. Therefore, the analysis of karyotypes is a tool for species-level identification of these mosquitoes.

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 comparative analysis of the metaphase karyotypes of Aedes excrucians, Ae. behningi, and Ae. euedes (Diptera: Culicidae) imaginal disсs.

Karyotypes of Aedes (Culicidae) mosquitoes (Ae. excrucians, Ae. behningi, and Ae. euedes) have been analyzed using the metaphase chromosomes of imaginal discs. Lacto-aceto-orcein, C-banding, and DAPI staining have detected species-specific features in the morphology and lengths of these chromosomes in the examined species. Species-specific features of chromosome 1 in the location of heterochromatin blocks have been shown. Thus, the metaphase chromosomes in the imaginal discs of Ae. excrucians, Ae. behningi, and Ae. euedes are a characteristic for species identification of mosquito species.

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.

Mitotic-chromosome-based physical mapping of the Culex quinquefasciatus genome.

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.

Cyst geometry in the egg chambers of Calliphora erythrocephala Mg. (Diptera: Calliphoridae) ovaries.

In the germarium of polytrophic ovarioles of Calliphora erythrocephala (Mg.) fly, four mitotic divisions of cystoblasts give rise to 16-cell germ-line cysts. One cell differentiates into an oocyte, while the remaining 15 cells become nurse cells. Concomitantly actin-rich ring canals are formed at the intercellular junctions. The present study considers a mutual arrangement of the ring canals formed after the second to fourth mitoses relative to the ring canal formed after the first mitotic division in different regions of the germarium and egg chambers. During the cyst formation and its movement to the posterior end of the germarium, the ring canals are displaced relative to one another, thereby giving different branching variants of the cyst. The pattern of cell interconnections becomes stable in germarium region 2b and does not change during the cyst movement along the ovariole despite the cyst polarizes and increases in size.

The changes in chromosome 6 spatial organization during chromatin polytenization in the Calliphora erythrocephala Mg. (Diptera: Calliphoridae) nurse cells.

Localization of Calliphora erythrocephala chromosome 6 in a 3D nuclear space at different stages of nurse cell chromatin polytenization was analyzed by fluorescence in situ hybridization and 3D microscopy. The obtained results suggest a large-scale chromatin relocation in the C. erythrocephala nurse cell nuclei, which is accompanied by a change in the chromosome territory of chromosome 6 associated with the change in expression activity of the nucleus and formation of reticular chromatin structure. It was revealed that the relocation of chromosome 6 (nucleolus organizer chromosome) is accompanied by fragmentation of the single large nucleolus into micronucleoli, which are spread over the entire nuclear space being associated with their nucleolar organizer regions. Presumably, the chromosome 6 material during transition to a highly polytenized structure is redistributed in the nucleus so that the inactive pericentromeric regions are displaced to the nuclear periphery, while the chromosome regions carrying rDNA sequences loop out beyond the chromosome territory. Being dispersed over the entire nuclear space, rDNA sequences are likely to be amplified, thereby providing numerous small signals from the chromosome 6-specific DNA probe. Micronucleoli are formed around the actively transcribed nucleolar organizer regions.