Chromosomes are hierarchically folded within cell nuclei into territories, domains and subdomains, but the functional importance and evolutionary dynamics of these hierarchies are poorly defined. Here, we comprehensively profile genome organizations of five Anopheles mosquito species and show how different levels of chromatin architecture influence each other. Patterns observed on Hi-C maps are associated with known cytological structures, epigenetic profiles, and gene expression levels. Evolutionary analysis reveals conservation of chromatin architecture within synteny blocks for tens of millions of years and enrichment of synteny breakpoints in regions with increased genomic insulation. However, in-depth analysis shows a confounding effect of gene density on both insulation and distribution of synteny breakpoints, suggesting limited causal relationship between breakpoints and regions with increased genomic insulation. At the level of individual loci, we identify specific, extremely long-ranged looping interactions, conserved for ~100 million years. We demonstrate that the mechanisms underlying these looping contacts differ from previously described Polycomb-dependent interactions and clustering of active chromatin.

Chromosome level assemblies are accumulating in various taxonomic groups including mosquitoes. However, even in the few reference-quality mosquito assemblies, a significant portion of the heterochromatic regions including telomeres remain unresolved. Here we produce a de novo assembly of the New World malaria mosquito, Anopheles albimanus by integrating Oxford Nanopore sequencing, Illumina, Hi-C and optical mapping. This 172.6 Mbps female assembly, which we call AalbS3, is obtained by scaffolding polished large contigs (contig N50=13.7 Mbps) into three chromosomes. All chromosome arms end with telomeric repeats, which is the first in mosquito assemblies and represents a significant step towards the completion of a genome assembly. These telomeres consist of tandem repeats of a novel 30-32 bp telomeric repeat unit (TRU) and are confirmed by analysing the termini of long reads and through both chromosomal in situ hybridization and a Bal31 sensitivity assay. The AalbS3 assembly included previously uncharacterized centromeric and rDNA clusters and more than doubled the content of transposable elements and other repetitive sequences. This telomere-to-telomere assembly, although still containing gaps, represents a significant step towards resolving biologically important but previously hidden genomic components. The comparison of different scaffolding methods will also inform future efforts to obtain reference-quality genomes for other mosquito species.

Anopheles coluzzii and Anopheles arabiensis belong to the Anopheles gambiae complex and are among the major malaria vectors in sub-Saharan Africa. However, chromosome-level reference genome assemblies are still lacking for these medically important mosquito species. In this study, we produced de novo chromosome-level genome assemblies for A. coluzzii and A. arabiensis using the long-read Oxford Nanopore sequencing technology and the Hi-C scaffolding approach. We obtained 273.4 and 256.8 Mb of the total assemblies for A. coluzzii and A. arabiensis, respectively. Each assembly consists of 3 chromosome-scale scaffolds (X, 2, 3), complete mitochondrion, and unordered contigs identified as autosomal pericentromeric DNA, X pericentromeric DNA, and Y sequences. Comparison of these assemblies with the existing assemblies for these species demonstrated that we obtained improved reference-quality genomes. The new assemblies allowed us to identify genomic coordinates for the breakpoint regions of fixed and polymorphic chromosomal inversions in A. coluzzii and A. arabiensis. The new chromosome-level assemblies will facilitate functional and population genomic studies in A. coluzzii and A. arabiensis. The presented assembly pipeline will accelerate progress toward creating high-quality genome references for other disease vectors.

Hybrid male sterility (HMS) contributes to speciation by restricting gene flow between related taxa. Detailed cytological characterization of reproductive organs in hybrid males is important for identifying phenotypes that can help guide searches of speciation genes. To investigate possible cellular causes of HMS, we performed crosses between closely related species of theAnopheles gambiaecomplex:An. meruswithAn. gambiaeorAn. coluzzii. We demonstrate that HMS in African malaria mosquitoes involves two defects in the reciprocal crosses: a premeiotic arrest of germline stem cells in degenerate testes and a failure of the reductional meiotic division of primary spermatocytes in normal-like testes. The premeiotic arrest in degenerate testes of hybrids is accompanied by a strong suppression of meiotic and postmeiotic genes. Unlike pure species, sex chromosomes in normal-like testes of F1 hybrids are largely unpaired during meiotic prophase I and all chromosomes show various degrees of insufficient condensation. Instead of entering reductional division in meiosis I, primary spermatocytes prematurely undergo an equational mitotic division producing non-motile diploid sperm. Thus, our study identified cytogenetic errors in interspecies hybrids that arise during the early stages of postzygotic isolation.

Genes that originate during evolution are an important source of novel biological functions. Retrogenes are functional copies of genes produced by retroduplication and as such are located in different genomic positions. To investigate retroposition patterns and retrogene expression, we computationally identified interchromosomal retroduplication events in nine portions of the phylogenetic history of malaria mosquitoes, making use of species that do or do not have classical sex chromosomes to test the roles of sex-linkage. We found 40 interchromosomal events and a significant excess of retroduplications from the X chromosome to autosomes among a set of young retrogenes. These young retroposition events occurred within the last 100 million years in lineages where all species possessed differentiated sex chromosomes. An analysis of available microarray and RNA-seq expression data for Anopheles gambiae showed that many of the young retrogenes evolved male-biased expression in the reproductive organs. Young autosomal retrogenes with increased meiotic or postmeiotic expression in the testes tend to be male biased. In contrast, older retrogenes, i.e., in lineages with undifferentiated sex chromosomes, do not show this particular chromosomal bias and are enriched for female-biased expression in reproductive organs. Our reverse-transcription PCR data indicates that most of the youngest retrogenes, which originated within the last 47.6 million years in the subgenus Cellia, evolved non-uniform expression patterns across body parts in the males and females of An. coluzzii. Finally, gene annotation revealed that mitochondrial function is a prominent feature of the young autosomal retrogenes. We conclude that mRNA-mediated gene duplication has produced a set of genes that contribute to mosquito reproductive functions and that different biases are revealed after the sex chromosomes evolve. Overall, these results suggest potential roles for the evolution of meiotic sex chromosome inactivation in males and of sexually antagonistic conflict related to mitochondrial energy function as the main selective pressures for X-to-autosome gene reduplication and testis-biased expression in these mosquito lineages.

Spatial organization of chromosome territories and interactions between interphase chromosomes themselves, as well as with the nuclear periphery, play important roles in epigenetic regulation of the genome function. However, the interplay between inter-chromosomal contacts and chromosome-nuclear envelope attachments in an organism’s development is not well-understood. To address this question, we conducted microscopic analyses of the three-dimensional chromosome organization in malaria mosquitoes. We employed multi-colored oligonucleotide painting probes, spaced 1 Mb apart along the euchromatin, to quantitatively study chromosome territories in larval salivary gland cells and adult ovarian nurse cells of Anopheles gambiae, An. coluzzii, and An. merus. We found that the X chromosome territory has a significantly smaller volume and is more compact than the autosomal arm territories. The number of inter-chromosomal, and the percentage of the chromosome–nuclear envelope, contacts were conserved among the species within the same cell type. However, the percentage of chromosome regions located at the nuclear periphery was typically higher, while the number of inter-chromosomal contacts was lower, in salivary gland cells than in ovarian nurse cells. The inverse correlation was considerably stronger for the autosomes. Consistent with previous theoretical arguments, our data indicate that, at the genome-wide level, there is an inverse relationship between chromosome-nuclear envelope attachments and chromosome–chromosome interactions, which is a key feature of the cell type-specific nuclear architecture.