The chromosome-scale genome assembly for the West Nile vector Culex quinquefasciatus uncovers patterns of genome evolution in mosquitoes.

BACKGROUND: Understanding genome organization and evolution is important for species involved in transmission of human diseases, such as mosquitoes. Anophelinae and Culicinae subfamilies of mosquitoes show striking differences in genome sizes, sex chromosome arrangements, behavior, and ability to transmit pathogens. However, the genomic basis of these differences is not fully understood. METHODS: In this study, we used a combination of advanced genome technologies such as Oxford Nanopore Technology sequencing, Hi-C scaffolding, Bionano, and cytogenetic mapping to develop an improved chromosome-scale genome assembly for the West Nile vector Culex quinquefasciatus. RESULTS: We then used this assembly to annotate odorant receptors, odorant binding proteins, and transposable elements. A genomic region containing male-specific sequences on chromosome 1 and a polymorphic inversion on chromosome 3 were identified in the Cx. quinquefasciatus genome. In addition, the genome of Cx. quinquefasciatus was compared with the genomes of other mosquitoes such as malaria vectors An. coluzzi and An. albimanus, and the vector of arboviruses Ae. aegypti. Our work confirms significant expansion of the two chemosensory gene families in Cx. quinquefasciatus, as well as a significant increase and relocation of the transposable elements in both Cx. quinquefasciatus and Ae. aegypti relative to the Anophelines. Phylogenetic analysis clarifies the divergence time between the mosquito species. Our study provides new insights into chromosomal evolution in mosquitoes and finds that the X chromosome of Anophelinae and the sex-determining chromosome 1 of Culicinae have a significantly higher rate of evolution than autosomes. CONCLUSION: The improved Cx. quinquefasciatus genome assembly uncovered new details of mosquito genome evolution and has the potential to speed up the development of novel vector control strategies.

On the Origin and Evolution of the Mosquito Male-determining Factor Nix.

The mosquito family Culicidae is divided into 2 subfamilies named the Culicinae and Anophelinae. Nix, the dominant male-determining factor, has only been found in the culicines Aedes aegypti and Aedes albopictus, 2 important arboviral vectors that belong to the subgenus Stegomyia. Here we performed sex-specific whole-genome sequencing and RNAseq of divergent mosquito species and explored additional male-inclusive datasets to investigate the distribution of Nix. Except for the Culex genus, Nix homologs were found in all species surveyed from the Culicinae subfamily, including 12 additional species from 3 highly divergent tribes comprising 4 genera, suggesting Nix originated at least 133 to 165 million years ago (MYA). Heterologous expression of 1 of 3 divergent Nix open reading frames (ORFs) in Ae. aegypti resulted in partial masculinization of genetic females as evidenced by morphology and doublesex splicing. Phylogenetic analysis suggests Nix is related to femaleless (fle), a recently described intermediate sex-determining factor found exclusively in anopheline mosquitoes. Nix from all species has a conserved structure, including 3 RNA-recognition motifs (RRMs), as does fle. However, Nix has evolved at a much faster rate than fle. The RRM3 of both Nix and fle are distantly related to the single RRM of a widely distributed and conserved splicing factor transformer-2 (tra2). The RRM3-based phylogenetic analysis suggests this domain in Nix and fle may have evolved from tra2 or a tra2-related gene in a common ancestor of mosquitoes. Our results provide insights into the evolution of sex determination in mosquitoes and will inform broad applications of mosquito-control strategies based on manipulating sex ratios toward nonbiting males.

Marker-assisted mapping enables forward genetic analysis in Aedes aegypti, an arboviral vector with vast recombination deserts.

Aedes aegypti is a major vector of arboviruses that cause dengue, chikungunya, yellow fever, and Zika. Although recent success in reverse genetics has facilitated rapid progress in basic and applied research, integration of forward genetics with modern technologies remains challenging in this important species, as up to 47% of its chromosome is refractory to genetic mapping due to extremely low rate of recombination. Here, we report the development of a marker-assisted mapping strategy to readily screen for and genotype only the rare but informative recombinants, drastically increasing both the resolution and signal-to-noise ratio. Using marker-assisted mapping, we mapped a transgene that was inserted in a >100-Mb recombination desert and a sex-linked spontaneous red-eye (re) mutation just outside the region. We subsequently determined, by CRISPR/Cas9-mediated knockout, that cardinal is the causal gene of re, which is the first forward genetic identification of a causal gene in Ae. aegypti. The identification of the causal gene of the sex-linked re mutation provides the molecular foundation for using gene editing to develop versatile and stable genetic sexing methods. To facilitate genome-wide forward genetics in Ae. aegypti, we generated and compiled a number of lines with markers throughout the genome. Thus, by overcoming the challenges presented by the vast recombination deserts and the scarcity of markers, we have shown that effective forward genetic analysis is increasingly feasible in this important arboviral vector species.

The Beginning of the End: A Chromosomal Assembly of the New World Malaria Mosquito Ends with a Novel Telomere.

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

Activation of SsoPK4, an Archaeal eIF2α Kinase Homolog, by Oxidized CoA.

The eukaryotic protein kinase (ePK) paradigm provides integral components for signal transduction cascades throughout nature. However, while so-called typical ePKs permeate the Eucarya and Bacteria, atypical ePKs dominate the kinomes of the Archaea. Intriguingly, the catalytic domains of the handful of deduced typical ePKs from the archaeon Sulfolobus solfataricus P2 exhibit significant resemblance to the protein kinases that phosphorylate translation initiation factor 2α (eIF2α) in response to cellular stresses. We cloned and expressed one of these archaeal eIF2α protein kinases, SsoPK4. SsoPK4 exhibited protein-serine/threonine kinase activity toward several proteins, including the S. solfataricus homolog of eIF2α, aIF2α. The activity of SsoPK4 was inhibited in vitro by 3',5'-cyclic AMP (Ki of ~23 µM) and was activated by oxidized Coenzyme A, an indicator of oxidative stress in the Archaea. Activation enhanced the apparent affinity for protein substrates, Km, but had little effect on Vmax. Autophosphorylation activated SsoPK4 and rendered it insensitive to oxidized Coenzyme A.