A Near-Chromosome Level Genome Assembly of Anopheles stephensi

Interspecies genome divergence is predominantly due to frequent small scale rearrangements in Eucalyptus

Synteny, the ordering of sequences within homologous chromosomes, must be maintained within the genomes of sexually reproducing species for the sharing of alleles and production of viable, reproducing offspring. However, when the genomes of closely related species are compared, a loss of synteny is often observed. Unequal homologous recombination is the primary mechanism behind synteny loss, occurring more often in transposon rich regions, and resulting in the formation of chromosomal rearrangements. To examine patterns of synteny among three closely related, interbreeding, and wild Eucalyptus species, we assembled their genomes using long-read DNA sequencing and de novo assembly. We identify syntenic and rearranged regions between these genomes and estimate that ~48% of our genomes remain syntenic while ~36% is rearranged. We observed that rearrangements highly fragment microsynteny. Our results suggest that synteny between these species is primarily lost through small-scale rearrangements, not through sequence loss, gain, or sequence divergence. Further examination of identified rearrangements suggests that rearrangements may be altering the phenotypes of Eucalyptus species. Our study also underscores that the use of single reference genomes in genomic variation studies could lead to reference bias, especially given the scale at which we show potentially adaptive loci have highly diverged, deleted, duplicated and/or rearranged. This study provides an unbiased framework to look at potential speciation and adaptive loci among a rapidly radiating foundation species of woodland trees that are free from selective breeding seen in most crop species.

Targetome Analysis of Malaria Sporozoite Transcription Factor AP2-Sp Reveals Its Role as a Master Regulator

Malaria transmission to humans begins with sporozoite infection of the liver. The elucidation of gene regulation during the sporozoite stage will promote the investigation of mechanisms of liver infection by this parasite and contribute to the development of strategies for preventing malaria transmission. AP2-Sp is a transcription factor (TF) essential for the formation of sporozoites or sporogony, which takes place in oocysts in the midguts of infected mosquitoes. To understand the role of this TF in the transcriptional regulatory system of this stage, we performed chromatin immunoprecipitation sequencing (ChIP-seq) analyses using whole mosquito midguts containing late oocysts as starting material and explored its genome-wide target genes. We identified 697 target genes, comprising those involved in distinct processes parasites experience during this stage, from sporogony to development into the liver stage and representing the majority of genes highly expressed in the sporozoite stage. These results suggest that AP2-Sp determines basal patterns of gene expression by targeting a broad range of genes directly. The ChIP-seq analyses also showed that AP2-Sp maintains its own expression by a transcriptional autoactivation mechanism (positive-feedback loop) and induces all TFs reported to be transcribed at this stage, including AP2-Sp2, AP2-Sp3, and SLARP. The results showed that AP2-Sp exists at the top of the transcriptional cascade of this stage and triggers the formation of this stage as a master regulator. IMPORTANCE The sporozoite stage plays a central role in malaria transmission from a mosquito to vertebrate host and is an important target for antimalarial strategies. AP2-Sp is a candidate master transcription factor for the sporozoite stage. However, study of its role in gene regulation has been hampered because of difficulties in performing genome-wide studies of gene regulation in this stage. Here, we conquered this problem and revealed that AP2-Sp has the following prominent features as a master transcription factor. First, it determines the repertory of gene expression during this stage. Second, it maintains its own expression through a transcriptional positive-feedback loop and induces all other transcription factors specifically expressed in this stage. This study represents a major breakthrough in fully understanding gene regulation in this important malarial stage.

Special Collection: Highlights of Medical, Urban and Veterinary Entomology. Highlights in Medical Entomology, 2021

Life remained far from normal as we completed the first year of the Covid-19 pandemic and entered a second year. Despite the challenges faced worldwide, together we continue to move the field of Medical Entomology forward. Here, I reflect on parallels between control of Covid-19 and vector-borne disease control, discuss the advantages and caveats of using new genotyping technologies for the study of invasive species, and proceed to highlight papers that were published between 2020 and 2021 with a focus on those related to mosquito surveillance and population genetics of mosquito vectors.

Identification of a TNF-TNFR-like system in malaria vectors (Anopheles stephensi) likely to influence Plasmodium resistance

Identification of Plasmodium-resistance genes in malaria vectors remains an elusive goal despite the recent availability of high-quality genomes of several mosquito vectors. Anopheles stephensi, with its three distinctly-identifiable forms at the egg stage, correlating with varying vector competence, offers an ideal species to discover functional mosquito genes implicated in Plasmodium resistance. Recently, the genomes of several strains of An. stephensi of the type-form, known to display high vectorial capacity, were reported. Here, we report a chromosomal-level assembly of an intermediate-form of An. stephensi strain (IndInt), shown to have reduced vectorial capacity relative to a strain of type-form (IndCh). The contig level assembly with a L50 of 4 was scaffolded into chromosomes by using the genome of IndCh as the reference. The final assembly shows a heterozygous paracentric inversion, 3Li, involving 8 Mbp, which is syntenic to the extensively-studied 2La inversion implicated in Plasmodium resistance in An. gambiae involving 21 Mbp. Deep annotation of genes within the 3Li region in the IndInt assembly using the state-of-the-art protein-fold prediction and other annotation tools reveals the presence of a tumor necrosis factor-alpha (TNF-alpha) like gene, which is the homolog of the Eiger gene in Drosophila. Subsequent chromosome-wide searches revealed homologs of Wengen (Wgn) and Grindelwald (Grnd) genes, which are known to be the receptors for Eiger in Drosophila. We have identified all the genes in IndInt required for Eiger-mediated signaling by analogy to the TNF-alpha system, suggesting the presence of a functionally-active Eiger signaling pathway in IndInt. Comparative genomics of the three type-forms with that of IndInt, reveals structurally disruptive mutations in Eiger gene in all three strains of the type-form, suggesting compromised innate immunity in the type-form as the likely cause of high vectorial capacity in these strains. This is the first report of the presence of a homolog of Eiger in malaria vectors, known to be involved in cell death in Drosophila, within an inversion region in IndInt syntenic to an inversion associated with Plasmodium resistance in An. gambiae.

Anopheles mosquitoes reveal new principles of 3D genome organization in insects

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.

The genome trilogy of Anopheles stephensi, an urban malaria vector, reveals structure of a locus associated with adaptation to environmental heterogeneity

Anopheles stephensi is the most menacing malaria vector to watch for in newly urbanising parts of the world. Its fitness is reported to be a direct consequence of the vector adapting to laying eggs in over-head water tanks with street-side water puddles polluted by oil and sewage. Large frequent inversions in the genome of malaria vectors are implicated in adaptation. We report the genome assembly of a strain of An. stephensi of the type-form, collected from a construction site from Chennai (IndCh) in 2016. The genome reported here with a L50 of 4, completes the trilogy of high-resolution genomes of strains with respect to a 16.5 Mbp 2Rb genotype in An. stephensi known to be associated with adaptation to environmental heterogeneity. Unlike the reported genomes of two other strains, STE2 (2R+^b/2Rb) and UCI (2Rb/2Rb), IndCh is found to be homozygous for the standard form (2R+^b/2R+^b). Comparative genome analysis revealed base-level details of the breakpoints and allowed extraction of 22,650 segregating SNPs for typing this inversion in populations. Whole genome sequencing of 82 individual mosquitoes from diverse geographical locations reveal that one third of both wild and laboratory populations maintain the heterozygous genotype of 2Rb. The large number of SNPs can be tailored to 1740 exonic SNPs enabling genotyping directly from transcriptome sequencing. The genome trilogy approach accelerated the study of fine structure and typing of an important inversion in An. stephensi, putting the genome resources for this understudied species on par with the extensively studied malaria vector, Anopheles gambiae. We argue that the IndCh genome is relevant for field translation work compared to those reported earlier by showing that individuals from diverse geographical locations cluster with IndCh, pointing to significant convergence resulting from travel and commerce between cities, perhaps, contributing to the survival of the fittest strain.

A chromosome-scale genome assembly and karyotype of the ctenophore Hormiphora californensis

Here, we present a karyotype, a chromosome-scale genome assembly, and a genome annotation from the ctenophore Hormiphora californensis (Ctenophora: Cydippida: Pleurobrachiidae). The assembly spans 110 Mb in 44 scaffolds and 99.47% of the bases are contained in 13 scaffolds. Chromosome micrographs and Hi-C heatmaps support a karyotype of 13 diploid chromosomes. Hi-C data reveal three large heterozygous inversions on chromosome 1, and one heterozygous inversion shares the same gene order found in the genome of the ctenophore Pleurobrachia bachei. We find evidence that H. californensis and P. bachei share thirteen homologous chromosomes, and the same karyotype of 1n = 13. The manually curated PacBio Iso-Seq-based genome annotation reveals complex gene structures, including nested genes and trans-spliced leader sequences. This chromosome-scale assembly is a useful resource for ctenophore biology and will aid future studies of metazoan evolution and phylogenetics.

Characterization of Anopheles stephensi Odorant Receptor 8, an Abundant Component of the Mouthpart Chemosensory Transcriptome

Several mosquito species within the genus Anopheles are vectors for human malaria, and the spread of this disease is driven by the propensity of certain species to feed preferentially on humans. The study of olfaction in mosquitoes is important to understand dynamics of host-seeking and host-selection; however, the majority of these studies focus on Anopheles gambiae or An. coluzzii, both vectors of malaria in Sub-Saharan Africa. Other malaria vectors may recognize different chemical cues from potential hosts; therefore, in this study, we investigated An. stephensi, the south Asian malaria mosquito. We specifically focused on the mouthparts (primarily the maxillary palp and labella) that have been much less investigated compared to the antennae but are also important for host-seeking. To provide a broad view of chemoreceptor expression, RNAseq was used to examine the transcriptomes from the mouthparts of host-seeking females, blood-fed females, and males. Notably, AsOr8 had a high transcript abundance in all transcriptomes and was, therefore, cloned and expressed in the Drosophila empty neuron system. This permitted characterization with a panel of odorants that were selected, in part, for their presence in the human odor profile. The responsiveness of AsOr8 to odorants was highly similar to An. gambiae Or8 (AgOr8), except for sulcatone, which was detected by AsOr8 but not AgOr8. Subtle differences in the receptor sensitivity to specific odorants may provide clues to species- or strain-specific approaches to host-seeking and host selection. Further exploration of the profile of An. stephensi chemosensory proteins may yield a better understanding of how different malaria vectors navigate host-finding and host-choice.

Improving mosquito control strategies with population genomics

Mosquito control strategies increasingly apply knowledge from population genomics research. This review highlights recent applications to three research domains: mosquito invasions, insecticide resistance evolution, and rear and release programs. Current research trends follow developments in reference assemblies, either as improvements to existing assemblies (particularly Aedes) or assemblies for new taxa (particularly Anopheles). With improved assemblies, studies of invasive and rear and release target populations are better able to incorporate adaptive as well as demographic hypotheses. New reference assemblies are aiding comparisons of insecticide resistance across sister taxa while helping resolve taxon boundaries amidst frequent introgression. Anopheles gene drive deployments and improved Aedes genome assemblies should lead to a convergence in research aims for Anopheles and Aedes in the coming years.

Hidden genomic features of an invasive malaria vector, Anopheles stephensi, revealed by a chromosome-level genome assembly

BACKGROUND

The mosquito Anopheles stephensi is a vector of urban malaria in Asia that recently invaded Africa. Studying the genetic basis of vectorial capacity and engineering genetic interventions are both impeded by limitations of a vector's genome assembly. The existing assemblies of An. stephensi are draft-quality and contain thousands of sequence gaps, potentially missing genetic elements important for its biology and evolution.

RESULTS

To access previously intractable genomic regions, we generated a reference-grade genome assembly and full transcript annotations that achieve a new standard for reference genomes of disease vectors. Here, we report novel species-specific transposable element (TE) families and insertions in functional genetic elements, demonstrating the widespread role of TEs in genome evolution and phenotypic variation. We discovered 29 previously hidden members of insecticide resistance genes, uncovering new candidate genetic elements for the widespread insecticide resistance observed in An. stephensi. We identified 2.4 Mb of the Y chromosome and seven new male-linked gene candidates, representing the most extensive coverage of the Y chromosome in any mosquito. By tracking full-length mRNA for > 15 days following blood feeding, we discover distinct roles of previously uncharacterized genes in blood metabolism and female reproduction. The Y-linked heterochromatin landscape reveals extensive accumulation of long-terminal repeat retrotransposons throughout the evolution and degeneration of this chromosome. Finally, we identify a novel Y-linked putative transcription factor that is expressed constitutively throughout male development and adulthood, suggesting an important role.

CONCLUSION

Collectively, these results and resources underscore the significance of previously hidden genomic elements in the biology of malaria mosquitoes and will accelerate the development of genetic control strategies of malaria transmission.