Search

A major insecticide resistance mechanism in insect pests is knock-down resistance (kdr) caused by mutations in the voltage-gated sodium channel (Vgsc) gene. Despite being common in most malaria Anopheles vector species, kdr mutations have never been observed in Anopheles funestus, the principal malaria vector in Eastern and Southern Africa, with resistance mainly being conferred by detoxification enzymes. In a parallel study, we monitored 10 populations of An. funestus in Tanzania for insecticide resistance unexpectedly identified resistance to a banned insecticide, DDT, in the Morogoro region. Through whole-genome sequencing of 333 An. funestus samples from these populations, we found eight novel amino acid substitutions in the Vgsc gene, including the kdr variant, L976F (equivalent to L995F in An. gambiae), in tight linkage disequilibrium with another (P1842S). The mutants were found only at high frequency in one region and were accompanied by weak signatures of a selective sweep, with a significant decline between 2017 and 2023. Notably, kdr L976F was strongly associated with survivorship to exposure to DDT insecticide, while no clear association was noted with a pyrethroid insecticide (deltamethrin). The WHO prequalifies no DDT products for vector control, and the chemical is banned in Tanzania. Widespread DDT contamination and a legacy of extensive countrywide stockpiles may have selected for this mutation. Continued monitoring is necessary to understand the origin of kdr in An. funestus, and the threat posed to insecticide-based vector control in Africa., (Molecular Ecology© 2024 The Author(s). Molecular Ecology published by John Wiley & Sons Ltd.)

Most insects encountered in the field are initially entomological dark matter in that they cannot be identified to species while alive. This explains the enduring quest for efficient ways to identify collected specimens. Morphological tools came first but are now routinely replaced or complemented with DNA barcodes. Initially too expensive for widespread use, these barcodes have since evolved into powerful tools for specimen identification and sorting, given that the evolution of sequencing approaches has dramatically reduced the cost of barcodes, thus enabling decentralized deployment across the planet. In this article, we review how DNA barcodes have become a key tool for accelerating biodiversity discovery and analyzing insect communities through both megabarcoding and metabarcoding in an era of insect decline. We predict that DNA barcodes will be particularly important for assembling image training sets for deep learning algorithms, global biodiversity genomics, and functional analysis of insect communities.

We present genome assembly from individual female An. coustani (African malaria mosquito; Arthropoda; Insecta; Diptera; Culicidae) from Lopé, Gabon. The genome sequence is 270 megabases in span. Most of the assembly is scaffolded into three chromosomal pseudomolecules with the X sex chromosome assembled for both species. The complete mitochondrial genome was also assembled and is 15.4 kilobases in length., Competing Interests: No competing interests were disclosed., (Copyright: © 2024 Bouafou LBA et al.)

Biodiversity genomics research requires reliable organismal identification, which can be difficult based on morphology alone. DNA-based identification using DNA barcoding can provide confirmation of species identity and resolve taxonomic issues but is rarely used in studies generating reference genomes. Here, we describe the development and implementation of DNA barcoding for the Darwin Tree of Life Project (DToL), which aims to sequence and assemble high quality reference genomes for all eukaryotic species in Britain and Ireland. We present a standardised framework for DNA barcode sequencing and data interpretation that is then adapted for diverse organismal groups. DNA barcoding data from over 12,000 DToL specimens has identified up to 20% of samples requiring additional verification, with 2% of seed plants and 3.5% of animal specimens subsequently having their names changed. We also make recommendations for future developments using new sequencing approaches and streamlined bioinformatic approaches., Competing Interests: No competing interests were disclosed., (Copyright: © 2024 Twyford AD et al.)

Large-scale reference genome sequencing projects for all of biodiversity are underway and common standards have been in place for some years to enable the understanding and sharing of sequence data. However, the metadata that describes the collection, processing and management of samples, and link to the associated sequencing and genome data, are not yet adequately developed and standardised for these projects. At the time of writing, the Darwin Tree of Life (DToL) Project is over two years into its ten-year ambition to sequence all described eukaryotic species in Britain and Ireland. We have sought consensus from a wide range of scientists across taxonomic domains to determine the minimal set of metadata that we collectively deem as critically important to accompany each sequenced specimen. These metadata are made available throughout the subsequent laboratory processes, and once collected, need to be adequately managed to fulfil the requirements of good data management practice. Due to the size and scale of management required, software tools are needed. These tools need to implement rigorous development pathways and change management procedures to ensure that effective research data management of key project and sample metadata is maintained. Tracking of sample properties through the sequencing process is handled by Lab Information Management Systems (LIMS), so publication of the sequenced data is achieved via technical integration of LIMS and data management tools. Discussions with community members on how metadata standards need to be managed within large-scale programmes is a priority in the planning process. Here we report on the standards we developed with respect to a robust and reusable mechanism of metadata collection, in the hopes that other projects forthcoming or underway will adopt these practices for metadata., Competing Interests: No competing interests were disclosed., (Copyright: © 2024 Shaw F et al.)

The adaptation of Anopheles malaria vectors to domestic settings is directly linked to their ability to feed on humans. The strength of this species-habitat association is unequal across the species within the genus, with the major vectors being particularly dependent on humans. However, our understanding of how blood-feeding behavior interacts with and adapts to environmental settings, including the presence of humans, remains limited. Using a field-based approach, we first investigated Anopheles community structure and feeding behavior patterns in domestic and sylvatic settings in La Lopé National Park in Gabon, Central Africa. We characterized the preference indices using a dual-host choice sampling approach across mosquito species, habitats, and seasons. We then quantified the plastic biting behavior of mosquito species in each habitat. We collected individuals from 16 Anopheles species that exhibited significant differences in species composition and abundance between sylvatic and domestic settings. The host-seeking behavior also varied among the seven most abundant species. The general attractiveness to each host, human or animal, remained relatively constant for each species, but with significant variations between habitats across species. These variations, to more generalist and to more anthropophilic behavior, were related to seasonal changes and distance from the village, respectively. Finally, we pointed out that the host choice of major malaria vectors changed in the absence of humans, revealing a plastic feeding behavior of these species. This study highlights the effect of humans on Anopheles distribution and feeding evolution. The characterization of feeding behavior in wild and domestic settings provides opportunities to better understand the interplay between genetic determinants of host preference and ecological factors. Our findings suggest that protected areas may offer alternative thriving conditions to major malaria vectors., Competing Interests: None declared., (© 2024 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

The developmental decision made by malaria parasites to become sexual underlies all malaria transmission. Here, we describe a rich atlas of short- and long-read single-cell transcriptomes of over 37,000 Plasmodium falciparum cells across intraerythrocytic asexual and sexual development. We used the atlas to explore transcriptional modules and exon usage along sexual development and expanded it to include malaria parasites collected from four Malian individuals naturally infected with multiple P. falciparum strains. We investigated genotypic and transcriptional heterogeneity within and among these wild strains at the single-cell level, finding differential expression between different strains even within the same host. These data are a key addition to the Malaria Cell Atlas interactive data resource, enabling a deeper understanding of the biology and diversity of transmission stages.

We present a genome assembly from an individual female Anopheles gambiae (the malaria mosquito; Arthropoda; Insecta; Diptera; Culicidae), Ifakara strain. The genome sequence is 264 megabases in span. Most of the assembly is scaffolded into three chromosomal pseudomolecules with the X sex chromosome assembled. The complete mitochondrial genome was also assembled and is 15.4 kilobases in length., Competing Interests: No competing interests were disclosed., (Copyright: © 2024 Habtewold T et al.)

Microsporidia are obligate intracellular eukaryotic parasites with an extremely broad host range. They have both economic and public health importance. Ploidy in microsporidia is variable, with a few species formally identified as diploid and one as polyploid. Given the increase in the number of studies sequencing microsporidian genomes, it is now possible to assess ploidy levels across all currently explored microsporidian diversity. We estimate ploidy for all microsporidian data sets available on the Sequence Read Archive using k-mer-based analyses, indicating that polyploidy is widespread in Microsporidia and that ploidy change is dynamic in the group. Using genome-wide heterozygosity estimates, we also show that polyploid microsporidian genomes are relatively homozygous, and we discuss the implications of these findings on the timing of polyploidization events and their origin.IMPORTANCEMicrosporidia are single-celled intracellular parasites, distantly related to fungi, that can infect a broad range of hosts, from humans all the way to protozoans. Exploiting the wealth of microsporidian genomic data available, we use k-mer-based analyses to assess ploidy status across the group. Understanding a genome's ploidy is crucial in order to assemble it effectively and may also be relevant for better understanding a parasite's behavior and life cycle. We show that tetraploidy is present in at least six species in Microsporidia and that these polyploidization events are likely to have occurred independently. We discuss why these findings may be paradoxical, given that Microsporidia, like other intracellular parasites, have extremely small, reduced genomes., Competing Interests: The authors declare no conflict of interest.

We present a genome assembly from an individual male Bombylius major (the dark-edged bee fly; Arthropoda; Insecta; Diptera; Bombyliidae). The genome sequence is 304.3 megabases in span. The whole assembly is scaffolded into 7 chromosomal pseudomolecules, including the X and Y sex chromosomes. The mitochondrial genome has also been assembled and is 17.8 kilobases in length. Gene annotation of this assembly on Ensembl identified 10,852 protein coding genes., Competing Interests: No competing interests were disclosed., (Copyright: © 2023 Lawniczak MKN et al.)

We present a genome assembly from an individual female Shargacucullia verbasci (the Mullein moth; Arthropoda; Insecta; Lepidoptera; Noctuidae). The genome sequence is 422.7 megabases in span. Most of the assembly is scaffolded into 32 chromosomal pseudomolecules, including the W and Z sex chromosomes. The mitochondrial genome has also been assembled and is 15.32 kilobases in length., Competing Interests: No competing interests were disclosed., (Copyright: © 2023 Lawniczak MKN et al.)

Insects are crucial for ecosystem health but climate change and pesticide use are driving massive insect decline. To mitigate this loss, we need new and effective monitoring techniques. Over the past decade there has been a shift to DNA-based techniques. We describe key emerging techniques for sample collection. We suggest that the selection of tools should be broadened, and that DNA-based insect monitoring data need to be integrated more rapidly into policymaking. We argue that there are four key areas for advancement, including the generation of more complete DNA barcode databases to interpret molecular data, standardisation of molecular methods, scaling up of monitoring efforts, and integrating molecular tools with other technologies that allow continuous, passive monitoring based on images and/or laser imaging, detection, and ranging (LIDAR)., Competing Interests: Declaration of interests The authors declare no conflicts of interest., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Malaria transmission to mosquitoes requires a developmental switch in asexually dividing blood-stage parasites to sexual reproduction. In Plasmodium berghei, the transcription factor AP2-G is required and sufficient for this switch, but how a particular sex is determined in a haploid parasite remains unknown. Using a global screen of barcoded mutants, we here identify genes essential for the formation of either male or female sexual forms and validate their importance for transmission. High-resolution single-cell transcriptomics of ten mutant parasites portrays the developmental bifurcation and reveals a regulatory cascade of putative gene functions in the determination and subsequent differentiation of each sex. A male-determining gene with a LOTUS/OST-HTH domain as well as the protein interactors of a female-determining zinc-finger protein indicate that germ-granule-like ribonucleoprotein complexes complement transcriptional processes in the regulation of both male and female development of a malaria parasite., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

The ANOSPP amplicon panel is a genus-wide targeted sequencing panel to facilitate large-scale monitoring of Anopheles species diversity. Combining information from the 62 nuclear amplicons present in the ANOSPP panel allows for a more senstive and specific species assignment than single gene (e.g. COI) barcoding, which is desirable in the light of permeable species boundaries. Here, we present NNoVAE, a method using Nearest Neighbours (NN) and Variational Autoencoders (VAE), which we apply to k- mers resulting from the ANOSPP amplicon sequences in order to hierarchically assign species identity. The NN step assigns a sample to a species-group by comparing the k -mers arising from each haplotype's amplicon sequence to a reference database. The VAE step is required to distinguish between closely related species, and also has sufficient resolution to reveal population structure within species. In tests on independent samples with over 80% amplicon coverage, NNoVAE correctly classifies to species level 98% of samples within the An. gambiae complex and 89% of samples outside the complex. We apply NNoVAE to over two thousand new samples from Burkina Faso and Gabon, identifying unexpected species in Gabon. NNoVAE presents an approach that may be of value to other targeted sequencing panels, and is a method that will be used to survey Anopheles species diversity and Plasmodium transmission patterns through space and time on a large scale, with plans to analyse half a million mosquitoes in the next five years., Competing Interests: MB, AM, DA, LB, AD, UE, MK, GL, BM, MN, OO, NR, FT, RD, ML No competing interests declared, (© 2022, Boddé et al.)

We present a genome assembly from an individual Patella pellucida (the blue-rayed limpet; Mollusca; Gastropoda; Patellidae). The genome sequence is 712 megabases in span. The majority of the assembly (99.85%) is scaffolded into 9 chromosomal pseudomolecules. The mitochondrial genome was assembled and is 14.9 kilobases in length., Competing Interests: No competing interests were disclosed., (Copyright: © 2022 Lawniczak MKN et al.)

Early diverging lineages such as trypanosomes can provide clues to the evolution of sexual reproduction in eukaryotes. In Trypanosoma brucei, the pathogen that causes Human African Trypanosomiasis, sexual reproduction occurs in the salivary glands of the insect host, but analysis of the molecular signatures that define these sexual forms is complicated because they mingle with more numerous, mitotically-dividing developmental stages. We used single-cell RNA-sequencing (scRNAseq) to profile 388 individual trypanosomes from midgut, proventriculus, and salivary glands of infected tsetse flies allowing us to identify tissue-specific cell types. Further investigation of salivary gland parasite transcriptomes revealed fine-scale changes in gene expression over a developmental progression from putative sexual forms through metacyclics expressing variant surface glycoprotein genes. The cluster of cells potentially containing sexual forms was characterized by high level transcription of the gamete fusion protein HAP2, together with an array of surface proteins and several genes of unknown function. We linked these expression patterns to distinct morphological forms using immunofluorescence assays and reporter gene expression to demonstrate that the kinetoplastid-conserved gene Tb927.10.12080 is exclusively expressed at high levels by meiotic intermediates and gametes. Further experiments are required to establish whether this protein, currently of unknown function, plays a role in gamete formation and/or fusion., Competing Interests: The authors have declared that no competing interests exist.

Life on Earth has evolved from initial simplicity to the astounding complexity we experience today. Bacteria and archaea have largely excelled in metabolic diversification, but eukaryotes additionally display abundant morphological innovation. How have these innovations come about and what constraints are there on the origins of novelty and the continuing maintenance of biodiversity on Earth? The history of life and the code for the working parts of cells and systems are written in the genome. The Earth BioGenome Project has proposed that the genomes of all extant, named eukaryotes-about 2 million species-should be sequenced to high quality to produce a digital library of life on Earth, beginning with strategic phylogenetic, ecological, and high-impact priorities. Here we discuss why we should sequence all eukaryotic species, not just a representative few scattered across the many branches of the tree of life. We suggest that many questions of evolutionary and ecological significance will only be addressable when whole-genome data representing divergences at all of the branchings in the tree of life or all species in natural ecosystems are available. We envisage that a genomic tree of life will foster understanding of the ongoing processes of speciation, adaptation, and organismal dependencies within entire ecosystems. These explorations will resolve long-standing problems in phylogenetics, evolution, ecology, conservation, agriculture, bioindustry, and medicine., Competing Interests: The authors declare no competing interest., (Copyright © 2022 the Author(s). Published by PNAS.)

A global international initiative, such as the Earth BioGenome Project (EBP), requires both agreement and coordination on standards to ensure that the collective effort generates rapid progress toward its goals. To this end, the EBP initiated five technical standards committees comprising volunteer members from the global genomics scientific community: Sample Collection and Processing, Sequencing and Assembly, Annotation, Analysis, and IT and Informatics. The current versions of the resulting standards documents are available on the EBP website, with the recognition that opportunities, technologies, and challenges may improve or change in the future, requiring flexibility for the EBP to meet its goals. Here, we describe some highlights from the proposed standards, and areas where additional challenges will need to be met., Competing Interests: Competing interest statement: P.F. is a member of the scientific advisory boards of Fabric Genomics, Inc. and Eagle Genomics, Ltd., (Copyright © 2022 the Author(s). Published by PNAS.)

Competing Interests: The authors declare no competing interest.

Anopheles is a diverse genus of mosquitoes comprising over 500 described species, including all known human malaria vectors. While a limited number of key vector species have been studied in detail, the goal of malaria elimination calls for surveillance of all potential vector species. Here, we develop a multilocus amplicon sequencing approach that targets 62 highly variable loci in the Anopheles genome and two conserved loci in the Plasmodium mitochondrion, simultaneously revealing both the mosquito species and whether that mosquito carries malaria parasites. We also develop a cheap, nondestructive, and high-throughput DNA extraction workflow that provides template DNA from single mosquitoes for the multiplex PCR, which means specimens producing unexpected results can be returned to for morphological examination. Over 1000 individual mosquitoes can be sequenced in a single MiSeq run, and we demonstrate the panel's power to assign species identity using sequencing data for 40 species from Africa, Southeast Asia, and South America. We also show that the approach can be used to resolve geographic population structure within An. gambiae and An. coluzzii populations, as the population structure determined based on these 62 loci from over 1000 mosquitoes closely mirrors that revealed through whole genome sequencing. The end-to-end approach is quick, inexpensive, robust, and accurate, which makes it a promising technique for very large-scale mosquito genetic surveillance and vector control., (© 2021 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)

We present a genome assembly from an individual Marthasterias glacialis (the spiny starfish; Echinodermata; Asteroidea; Forcipulatida; Asteriidae). The genome sequence is 521 megabases in span. The majority of the assembly, 99.44%, is scaffolded into 22 chromosomal pseudomolecules. The mitochondrial genome has also been assembled, and is 16 kb in span., Competing Interests: No competing interests were disclosed., (Copyright: © 2021 Lawniczak MKN et al.)

Museum collections contain enormous quantities of insect specimens collected over the past century, covering a period of increased and varied insecticide usage. These historic collections are therefore incredibly valuable as genomic snapshots of organisms before, during, and after exposure to novel selective pressures. However, these samples come with their own challenges compared with present-day collections, as they are fragile and retrievable DNA is low yield and fragmented. In this article, we tested several DNA extraction procedures across pinned historic Diptera specimens from four disease vector genera: Anopheles, Aedes, Culex, and Glossina. We identify an approach that minimizes morphological damage while maximizing DNA retrieval for Illumina library preparation and sequencing that can accommodate the fragmented and low yield nature of historic DNA. We identify several key points in retrieving sufficient DNA while keeping morphological damage to a minimum: an initial rehydration step, a short incubation without agitation in a modified low salt Proteinase K buffer (referred to as "lysis buffer C" throughout), and critical point drying of samples post-extraction to prevent tissue collapse caused by air drying. The suggested method presented here provides a solid foundation for exploring the genomes and morphology of historic Diptera collections., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Malaria parasites have a complex life cycle featuring diverse developmental strategies, each uniquely adapted to navigate specific host environments. Here we use single-cell transcriptomics to illuminate gene usage across the transmission cycle of the most virulent agent of human malaria - Plasmodium falciparum. We reveal developmental trajectories associated with the colonization of the mosquito midgut and salivary glands and elucidate the transcriptional signatures of each transmissible stage. Additionally, we identify both conserved and non-conserved gene usage between human and rodent parasites, which point to both essential mechanisms in malaria transmission and species-specific adaptations potentially linked to host tropism. Together, the data presented here, which are made freely available via an interactive website, provide a fine-grained atlas that enables intensive investigation of the P. falciparum transcriptional journey. As well as providing insights into gene function across the transmission cycle, the atlas opens the door for identification of drug and vaccine targets to stop malaria transmission and thereby prevent disease.

Background: The emergence of insecticide resistance is a major threat to malaria control programmes in Africa, with many different factors contributing to insecticide resistance in its vectors, Anopheles mosquitoes. CYP6M2 has previously been recognized as an important candidate in cytochrome P450-mediated detoxification in Anopheles. As it has been implicated in resistance against pyrethroids, organochlorines and carbamates, its broad metabolic activity makes it a potential agent in insecticide cross-resistance. Currently, allelic variation within the Cyp6m2 gene remains unknown., Methods: Here, Illumina whole-genome sequence data from Phase 2 of the Anopheles gambiae 1000 Genomes Project (Ag1000G) was used to examine genetic variation in the Cyp6m2 gene across 16 populations in 13 countries comprising Anopheles gambiae and Anopheles coluzzii mosquitoes. To identify whether these alleles show evidence of selection either through potentially modified enzymatic function or by being linked to variants that change the transcriptional profile of the gene, hierarchical clustering of haplotypes, linkage disequilibrium, median joining networks and extended haplotype homozygosity analyses were performed., Results: Fifteen missense biallelic substitutions at high frequency (defined as > 5% frequency in one or more populations) are found, which fall into five distinct haplotype groups that carry the main high frequency variants: A13T, D65A, E328Q, Y347F, I359V and A468S. Despite consistent reports of Cyp6m2 upregulation and metabolic activity in insecticide resistant Anophelines, no evidence of directional selection is found occurring on these variants or on the haplotype clusters in which they are found., Conclusion: These results imply that emerging resistance associated with Cyp6m2 is potentially driven by distant regulatory loci such as transcriptional factors rather than by its missense variants, or that other genes are playing a more significant role in conferring metabolic resistance.

The crossing of the mosquito midgut epithelium by the malaria parasite motile ookinete form represents the most extreme population bottleneck in the parasite life cycle and is a prime target for transmission blocking strategies. However, we have little understanding of the clonal variation that exists in a population of ookinetes in the vector, partially because the parasites are difficult to access and are found in low numbers. Within a vector, variation may result as a response to specific environmental cues or may exist independent of those cues as a potential bet-hedging strategy. Here we use single-cell RNA-seq to profile transcriptional variation in Plasmodium berghei ookinetes across different vector species, and between and within individual midguts. We then compare our results to low-input transcriptomes from individual Anopheles coluzzii midguts infected with the human malaria parasite Plasmodium falciparum . Although the vast majority of transcriptional changes in ookinetes are driven by development, we have identified candidate genes that may be responding to environmental cues or are clonally variant within a population. Our results illustrate the value of single-cell and low-input technologies in understanding clonal variation of parasite populations., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Witmer, Dahalan, Metcalf, Talman, Howick and Lawniczak.)

Resistance to artemisinin-based combination therapy (ACT) in the Plasmodium falciparum parasite is threatening to reverse recent gains in reducing global deaths from malaria. While resistance manifests as delayed parasite clearance in patients, the phenotype can only spread geographically via the sexual stages and mosquito transmission. In addition to their asexual killing properties, artemisinin and its derivatives sterilize sexual male gametocytes. Whether resistant parasites overcome this sterilizing effect has not, however, been fully tested. Here, we analyzed P. falciparum clinical isolates from the Greater Mekong Subregion, each demonstrating delayed clinical clearance and known resistance-associated polymorphisms in the Kelch13 (PfK13 var ) gene. As well as demonstrating reduced asexual sensitivity to drug, certain PfK13 var isolates demonstrated a marked reduction in sensitivity to artemisinin in an in vitro male gamete formation assay. Importantly, this same reduction in sensitivity was observed when the most resistant isolate was tested directly in mosquito feeds. These results indicate that, under artemisinin drug pressure, while sensitive parasites are blocked, resistant parasites continue transmission. This selective advantage for resistance transmission could favor acquisition of additional host-specificity or polymorphisms affecting partner drug sensitivity in mixed infections. Favored resistance transmission under ACT coverage could have profound implications for the spread of multidrug-resistant malaria beyond Southeast Asia., (Copyright © 2020 Witmer et al.)

Methods to deconvolve single-cell RNA-sequencing (scRNA-seq) data are necessary for samples containing a mixture of genotypes, whether they are natural or experimentally combined. Multiplexing across donors is a popular experimental design that can avoid batch effects, reduce costs and improve doublet detection. By using variants detected in scRNA-seq reads, it is possible to assign cells to their donor of origin and identify cross-genotype doublets that may have highly similar transcriptional profiles, precluding detection by transcriptional profile. More subtle cross-genotype variant contamination can be used to estimate the amount of ambient RNA. Ambient RNA is caused by cell lysis before droplet partitioning and is an important confounder of scRNA-seq analysis. Here we develop souporcell, a method to cluster cells using the genetic variants detected within the scRNA-seq reads. We show that it achieves high accuracy on genotype clustering, doublet detection and ambient RNA estimation, as demonstrated across a range of challenging scenarios.

A possible malaria control approach involves the dissemination in mosquitoes of inherited symbiotic microbes to block Plasmodium transmission. However, in the Anopheles gambiae complex, the primary African vectors of malaria, there are limited reports of inherited symbionts that impair transmission. We show that a vertically transmitted microsporidian symbiont (Microsporidia MB) in the An. gambiae complex can impair Plasmodium transmission. Microsporidia MB is present at moderate prevalence in geographically dispersed populations of An. arabiensis in Kenya, localized to the mosquito midgut and ovaries, and is not associated with significant reductions in adult host fecundity or survival. Field-collected Microsporidia MB infected An. arabiensis tested negative for P. falciparum gametocytes and, on experimental infection with P. falciparum, sporozoites aren't detected in Microsporidia MB infected mosquitoes. As a microbe that impairs Plasmodium transmission that is non-virulent and vertically transmitted, Microsporidia MB could be investigated as a strategy to limit malaria transmission.

In cells, organs and whole organisms, nutrient sensing is key to maintaining homeostasis and adapting to a fluctuating environment 1 . In many animals, nutrient sensors are found within the enteroendocrine cells of the digestive system; however, less is known about nutrient sensing in their cellular siblings, the absorptive enterocytes 1 . Here we use a genetic screen in Drosophila melanogaster to identify Hodor, an ionotropic receptor in enterocytes that sustains larval development, particularly in nutrient-scarce conditions. Experiments in Xenopus oocytes and flies indicate that Hodor is a pH-sensitive, zinc-gated chloride channel that mediates a previously unrecognized dietary preference for zinc. Hodor controls systemic growth from a subset of enterocytes-interstitial cells-by promoting food intake and insulin/IGF signalling. Although Hodor sustains gut luminal acidity and restrains microbial loads, its effect on systemic growth results from the modulation of Tor signalling and lysosomal homeostasis within interstitial cells. Hodor-like genes are insect-specific, and may represent targets for the control of disease vectors. Indeed, CRISPR-Cas9 genome editing revealed that the single hodor orthologue in Anopheles gambiae is an essential gene. Our findings highlight the need to consider the instructive contributions of metals-and, more generally, micronutrients-to energy homeostasis.

Plasmodium falciparum remains one of the leading causes of child mortality, and nearly half of the world's population is at risk of contracting malaria. While pathogenesis results from replication of asexual forms in human red blood cells, it is the sexually differentiated forms, gametocytes, which are responsible for the spread of the disease. For transmission to succeed, both mature male and female gametocytes must be taken up by a female Anopheles mosquito during its blood meal for subsequent differentiation into gametes and mating inside the mosquito gut. Observed circulating numbers of gametocytes in the human host are often surprisingly low. A pre-fertilization behavior, such as skin sequestration, has been hypothesized to explain the efficiency of human-to-mosquito transmission but has not been sufficiently tested due to a lack of appropriate tools. In this study, we describe the optimization of a qPCR tool that enables the relative quantification of gametocytes within very small input samples. Such a tool allows for the quantification of gametocytes in different compartments of the host and the vector that could potentially unravel mechanisms that enable highly efficient malaria transmission. We demonstrate the use of our gametocyte quantification method in mosquito blood meals from both direct skin feeding on Plasmodium gametocyte carriers and standard membrane feeding assay. Relative gametocyte abundance was not different between mosquitoes fed through a membrane or directly on the skin suggesting that there is no systematic enrichment of gametocytes picked up in the skin., (Copyright © 2020 Talman, Ouologuem, Love, Howick, Mulamba, Haidara, Dara, Sylla, Sacko, Coulibaly, Dao, Sangare, Djimde and Lawniczak.)

Mating causes dramatic changes in female physiology, behaviour, and immunity in many insects, inducing oogenesis, oviposition, and refractoriness to further mating. Females from the Anopheles gambiae species complex typically mate only once in their lifetime during which they receive sperm and seminal fluid proteins as well as a mating plug that contains the steroid hormone 20-hydroxyecdysone. This hormone, which is also induced by blood-feeding, plays a major role in activating vitellogenesis for egg production. Here we show that female Anopheles coluzzii susceptibility to Plasmodium falciparum infection is significantly higher in mated females compared to virgins. We also find that mating status has a major impact on the midgut transcriptome, detectable only under sugar-fed conditions: once females have blood-fed, the transcriptional changes that are induced by mating are likely masked by the widespread effects of blood-feeding on gene expression. To determine whether increased susceptibility to parasites could be driven by the additional 20E that mated females receive from males, we mimicked mating by injecting virgin females with 20E, finding that these females are significantly more susceptible to human malaria parasites than virgin females injected with the control 20E carrier. Further RNAseq was carried out to examine whether the genes that change upon 20E injection in the midgut are similar to those that change upon mating. We find that 79 midgut-expressed genes are regulated in common by both mating and 20E, and 96% (n = 76) of these are regulated in the same direction (up vs down in 20E/mated). Together, these findings show that male Anopheles mosquitoes induce changes in the female midgut that can affect female susceptibility to P. falciparum. This implies that in nature, males might contribute to malaria transmission in previously unappreciated ways, and that vector control strategies that target males may have additional benefits towards reducing transmission., Competing Interests: The authors have declared that no competing interests exist.

Malaria parasites adopt a remarkable variety of morphological life stages as they transition through multiple mammalian host and mosquito vector environments. We profiled the single-cell transcriptomes of thousands of individual parasites, deriving the first high-resolution transcriptional atlas of the entire Plasmodium berghei life cycle. We then used our atlas to precisely define developmental stages of single cells from three different human malaria parasite species, including parasites isolated directly from infected individuals. The Malaria Cell Atlas provides both a comprehensive view of gene usage in a eukaryotic parasite and an open-access reference dataset for the study of malaria parasites., (Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Polymorphisms in genetic copy number can influence gene expression, coding sequence, and zygosity, making them powerful actors in the evolutionary process. Copy number variants (CNVs) are however understudied, being more difficult to detect than single-nucleotide polymorphisms. We take advantage of the intense selective pressures on the major malaria vector Anopheles gambiae , caused by the widespread use of insecticides for malaria control, to investigate the role of CNVs in the evolution of insecticide resistance. Using the whole-genome sequencing data from 1142 samples in the An. gambiae 1000 genomes project, we identified 250 gene-containing CNVs, encompassing a total of 267 genes of which 28 were in gene families linked to metabolic insecticide resistance, representing significant enrichment of these families. The five major gene clusters for metabolic resistance all contained CNVs, with 44 different CNVs being found across these clusters and multiple CNVs frequently covering the same genes. These 44 CNVs are widespread (45% of individuals carry at least one of them) and have been spreading through positive selection, indicated by their high local frequencies and extended haplotype homozygosity. Our results demonstrate the importance of CNVs in the response to selection, highlighting the urgent need to identify the contribution of each CNV to insecticide resistance and to track their spread as the use of insecticides in malaria endemic countries intensifies and as the operational deployment of next-generation bed nets targeting metabolic resistance gathers pace. Our detailed descriptions of CNVs found across the species range provide the tools to do so., (© 2019 Lucas et al.; Published by Cold Spring Harbor Laboratory Press.)

A high-quality reference genome is a fundamental resource for functional genetics, comparative genomics, and population genomics, and is increasingly important for conservation biology. PacBio Single Molecule, Real-Time (SMRT) sequencing generates long reads with uniform coverage and high consensus accuracy, making it a powerful technology for de novo genome assembly. Improvements in throughput and concomitant reductions in cost have made PacBio an attractive core technology for many large genome initiatives, however, relatively high DNA input requirements (~5 µg for standard library protocol) have placed PacBio out of reach for many projects on small organisms that have lower DNA content, or on projects with limited input DNA for other reasons. Here we present a high-quality de novo genome assembly from a single Anopheles coluzzii mosquito. A modified SMRTbell library construction protocol without DNA shearing and size selection was used to generate a SMRTbell library from just 100 ng of starting genomic DNA. The sample was run on the Sequel System with chemistry 3.0 and software v6.0, generating, on average, 25 Gb of sequence per SMRT Cell with 20 h movies, followed by diploid de novo genome assembly with FALCON-Unzip. The resulting curated assembly had high contiguity (contig N50 3.5 Mb) and completeness (more than 98% of conserved genes were present and full-length). In addition, this single-insect assembly now places 667 (>90%) of formerly unplaced genes into their appropriate chromosomal contexts in the AgamP4 PEST reference. We were also able to resolve maternal and paternal haplotypes for over 1/3 of the genome. By sequencing and assembling material from a single diploid individual, only two haplotypes were present, simplifying the assembly process compared to samples from multiple pooled individuals. The method presented here can be applied to samples with starting DNA amounts as low as 100 ng per 1 Gb genome size. This new low-input approach puts PacBio-based assemblies in reach for small highly heterozygous organisms that comprise much of the diversity of life.

Understanding how phenotypic differences between males and females arise from the sex-biased expression of nearly identical genomes can reveal important insights into the biology and evolution of a species. Among Anopheles mosquito species, these phenotypic differences include vectorial capacity, as it is only females that blood feed and thus transmit human malaria. Here, we use RNA-seq data from multiple tissues of four vector species spanning the Anopheles phylogeny to explore the genomic and evolutionary properties of sex-biased genes. We find that, in these mosquitoes, in contrast to what has been found in many other organisms, female-biased genes are more rapidly evolving in sequence, expression, and genic turnover than male-biased genes. Our results suggest that this atypical pattern may be due to the combination of sex-specific life history challenges encountered by females, such as blood feeding. Furthermore, female propensity to mate only once in nature in male swarms likely diminishes sexual selection of post-reproductive traits related to sperm competition among males. We also develop a comparative framework to systematically explore tissue- and sex-specific splicing to document its conservation throughout the genus and identify a set of candidate genes for future functional analyses of sex-specific isoform usage. Finally, our data reveal that the deficit of male-biased genes on the X Chromosomes in Anopheles is a conserved feature in this genus and can be directly attributed to chromosome-wide transcriptional regulation that de-masculinizes the X in male reproductive tissues., (© 2017 Papa et al.; Published by Cold Spring Harbor Laboratory Press.)

Mosquitoes in the Anopheles gambiae complex show rapid ecological and behavioral diversification, traits that promote malaria transmission and complicate vector control efforts. A high-density, genome-wide mosquito SNP-genotyping array allowed mapping of genomic differentiation between populations and species that exhibit varying levels of reproductive isolation. Regions near centromeres or within polymorphic inversions exhibited the greatest genetic divergence, but divergence was also observed elsewhere in the genomes. Signals of natural selection within populations were overrepresented among genomic regions that are differentiated between populations, implying that differentiation is often driven by population-specific selective events. Complex genomic differentiation among speciating vector mosquito populations implies that tools for genome-wide monitoring of population structure will prove useful for the advancement of malaria eradication.