The genome sequence of the Mullein moth, Shargacucullia verbasci (Linnaeus, 1758)

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

Souporcell: robust clustering of single-cell RNA-seq data by genotype without reference genotypes.

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.We acknowledge the Wellcome Sanger Institute’s DNA Pipelines for construction of the 10x sequencing libraries. We thank Allan Muhwezi and Andrew Russell for assistance with parasite culture and 10x Single-cell 3’ RNA-seq respectively. In addition, we would like to thank Matthew Young for useful conversations about ambient RNA, Mirjana Efremova for providing information about the maternal/fetal data, and Katie Gray for assistance in interpreting the previously unannotated cluster. The Wellcome Sanger Institute is funded by the Wellcome Trust (grant 206194/Z/17/Z), which supports MKNL and MH. This work was supported by an MRC Career Development Award (G1100339) to MKNL. We would like to acknowledge the Wellcome Trust Sanger Institute as the source of the human induced pluripotent cell lines that were generated under the Human Induced Pluripotent Stem Cell Initiative funded by a grant from the Wellcome Trust and Medical Research Council, supported by the Wellcome Trust (WT098051) and the NIHR/Wellcome Trust Clinical Research Facility, and acknowledges Life Science Technologies Corporation as the provider of Cytotune (HipSci.org). The Cardiovascular Epidemiology Unit is supported by core funding from the UK Medical Research Council (MR/L003120/1), the British Heart Foundation (RG/13/13/30194; RG/18/13/33946) and the National Institute for Health Research [Cambridge Biomedical Research Centre at the Cambridge University Hospital’s NHS Foundation Trust]. The views expressed are those of the authors and not necessarily those of the NHS, the NIHR or the Department of Health and Social Care

Adaptive Gene Expression Divergence Inferred from Population Genomics

Detailed studies of individual genes have shown that gene expression divergence often results from adaptive evolution of regulatory sequence. Genome-wide analyses, however, have yet to unite patterns of gene expression with polymorphism and divergence to infer population genetic mechanisms underlying expression evolution. Here, we combined genomic expression data—analyzed in a phylogenetic context—with whole genome light-shotgun sequence data from six Drosophila simulans lines and reference sequences from D. melanogaster and D. yakuba. These data allowed us to use molecular population genetics to test for neutral versus adaptive gene expression divergence on a genomic scale. We identified recent and recurrent adaptive evolution along the D. simulans lineage by contrasting sequence polymorphism within D. simulans to divergence from D. melanogaster and D. yakuba. Genes that evolved higher levels of expression in D. simulans have experienced adaptive evolution of the associated 3′ flanking and amino acid sequence. Concomitantly, these genes are also decelerating in their rates of protein evolution, which is in agreement with the finding that highly expressed genes evolve slowly. Interestingly, adaptive evolution in 5′ cis-regulatory regions did not correspond strongly with expression evolution. Our results provide a genomic view of the intimate link between selection acting on a phenotype and associated genic evolution

A Minimally Morphologically Destructive Approach for DNA Retrieval and Whole-Genome Shotgun Sequencing of Pinned Historic Dipteran Vector Species

Abstract: 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.CopyrightThe Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited. The attached file is the published version of the article

Using scRNA-seq to identify transcriptional variation in the malaria parasite ookinete stage

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

Uptake of plasmodium falciparum gametocytes during mosquito bloodmeal by direct and membrane feeding

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

A High-Quality De novo Genome Assembly from a Single Mosquito Using PacBio Sequencing.

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

The Malaria Cell Atlas: single parasite transcriptomes across the complete Plasmodium life cycle

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

Sequence locally, think globally:The Darwin tree of life project

The goals of the Earth Biogenome Project—to sequence the genomes of all eukaryotic life on earth—are as daunting as they are ambitious. The Darwin Tree of Life Project was founded to demonstrate the credibility of these goals and to deliver at-scale genome sequences of unprecedented quality for a biogeographic region: the archipelago of islands that constitute Britain and Ireland. The Darwin Tree of Life Project is a collaboration between biodiversity organizations (museums, botanical gardens, and biodiversity institutes) and genomics institutes. Together, we have built a workflow that collects specimens from the field, robustly identifies them, performs sequencing, generates high-quality, curated assemblies, and releases these openly for the global community to use to build future science and conservation efforts.</jats:p

Whole-genome sequencing reveals high complexity of copy number variation at insecticide resistance loci in malaria mosquitoes

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