The genome sequence of the malaria mosquito Anopheles gambiae

De novo transcriptome analysis and identification of reproduction control genes from the red palm weevil Rhynchophorus ferrugineus

Recent attacks by the red palm weevil, Rhynchophorus ferrugineus (Olivier), have become a severe problem for palm species. In present work, fat body transcriptome of adult female red palm weevil was analyzed, focusing on the identification of reproduction control genes. Transcriptome study was completed by means of next-generation sequencing (NGS) using Illumina Hiseq 2000 sequencing system. A total of 105,938,182 raw reads, 102,645,544 clean reads, and 9,238,098,960 clean nucleotides with a guanine-cytosine content of 40.31%, were produced. The processed transcriptome data resulted in 43,789 unique transcripts (with mean lengths of 1,172 bp). It was found that 20% of total unique transcripts shared up to 80%-100% sequence identity with homologous species, mainly the mountain pine beetle Dendroctonus ponderosae (59.9%) and red flour beetle Tribolium castaneum (26.9%). Nearly 25 annotated genes were predicted to be involved in red palm weevil reproduction, including five vitellogenin (Vg) transcripts. Among the five Vg gene transcripts, one was highly expressed compared with the other four (FPKM values of 1.963, 1.471, 1.028, and 1.017, respectively), and the five Vg gene transcripts were designated as RfVg, RfVg-equivalent1, RfVg-equivalent2, RfVg-equivalent3, and RfVg-equivalent4, respectively. The high expression level of RfVg verified by RT-polymerase chain reaction analysis suggested that RfVg is the primary functional Vg gene in red palm weevil. A high similarity of RfVg with other Coleopterans was also reflected in a phylogenetic tree, where RfVg was placed within the clade of the order Coleoptera. Awareness of the major genes that play critical roles in reproduction and proliferation of red palm weevil is valuable to understand their reproduction mechanism at a molecular level. In addition, for future molecular studies, the NGS dataset obtained will be useful and will promote the exploration of biotech-based control strategies against red palm weevil, a primary pest of palm trees.

G-Protein Coupled Receptors (GPCRs) in Insects-A Potential Target for New Insecticide Development

G-protein coupled receptors (GPCRs) play important roles in cell biology and insects' physiological processes, toxicological response and the development of insecticide resistance. New information on genome sequences, proteomic and transcriptome analysis and expression patterns of GPCRs in organs such as the central nervous system in different organisms has shown the importance of these signaling regulatory GPCRs and their impact on vital cell functions. Our growing understanding of the role played by GPCRs at the cellular, genome, transcriptome and tissue levels is now being utilized to develop new targets that will sidestep many of the problems currently hindering human disease control and insect pest management. This article reviews recent work on the expression and function of GPCRs in insects, focusing on the molecular complexes governing the insect physiology and development of insecticide resistance and examining the genome information for GPCRs in two medically important insects, mosquitoes and house flies, and their orthologs in the model insect species Drosophila melanogaster. The tissue specific distribution and expression of the insect GPCRs is discussed, along with fresh insights into practical aspects of insect physiology and toxicology that could be fundamental for efforts to develop new, more effective, strategies for pest control and resistance management.

The 100 most-cited articles on malaria: a bibliometric analysis

Purpose

The use of citation analysis to identify the first 100 articles on malaria offers unique insights into understanding the disease and subsequent follow-up treatment innovations over time. In this study, the 100 most-cited articles on malaria were analyzed, and key studies were highlighted.

Design/methodology/approach

The data of the most-cited articles for the period of four decades were extracted from the Web of Science database. The search terms malaria, plasmodium infection and remittent fever were used to identify the related articles for the study.

Findings

The preliminary data of the 100 most-cited articles were recorded and analyzed. The total number of retrieved articles was 55,517. Dondrop and colleagues wrote the most-cited articles focusing on a new treatment for falciparum malaria resistant to existing medications. The author, with the most publications, was N. J. White. The most-cited articles on malaria were published in 35 journals. The USA had published most of the influential articles, while the Mahidol Oxford Tropical Medicine Research Unit was the leading institute active in malaria research. The experimental method was the most frequent method used by the most-cited articles. There was a significant relationship between the number of authors, the presence of the corresponding author on Research Gate, the availability of the full text of the article on Research Gate, the impact factor of the journal in which the article was published and the international collaboration of authors and the number of citations on malaria. The most influential authors, countries, institutes, journals and articles were specified. Most of the notable articles on malaria were published in journals with high impact factors. A group of journals was introduced as the core journals.

Originality/value

It was found that having co-authors, the presence of the corresponding author on Research Gate, the availability of the full text of the article on Research Gate, the impact factor of the journal in which the article was published and international collaboration contributed to the publication of high-quality scientific products. Updated information on malaria is needed to present and expand the screening strategies to improve health and reduce burden of malaria.

Regulators and signalling in insect antimicrobial innate immunity: Functional molecules and cellular pathways

Insects possess an immune system that protects them from attacks by various pathogenic microorganisms that would otherwise threaten their survival. Immune mechanisms may deal directly with the pathogens by eliminating them from the host organism or disarm them by suppressing the synthesis of toxins and virulence factors that promote the invasion and destructive action of the intruder within the host. Insects have been established as outstanding models for studying immune system regulation because innate immunity can be explored as an integrated system at the level of the whole organism. Innate immunity in insects consists of basal immunity that controls the constitutive synthesis of effector molecules such as antimicrobial peptides, and inducible immunity that is activated after detection of a microbe or its product(s). Activation and coordination of innate immune defenses in insects involve evolutionary conserved immune factors. Previous research in insects has led to the identification and characterization of distinct immune signalling pathways that modulate the response to microbial infections. This work has not only advanced the field of insect immunology, but it has also rekindled interest in the innate immune system of mammals. Here we review the current knowledge on key molecular components of insect immunity and discuss the opportunities they present for confronting infectious diseases in humans.

Systematic identification of plausible pathways to potential harm via problem formulation for investigational releases of a population suppression gene drive to control the human malaria vector Anopheles gambiae in West Africa

BACKGROUND

Population suppression gene drive has been proposed as a strategy for malaria vector control. A CRISPR-Cas9-based transgene homing at the doublesex locus (dsxFCRISPRh) has recently been shown to increase rapidly in frequency in, and suppress, caged laboratory populations of the malaria mosquito vector Anopheles gambiae. Here, problem formulation, an initial step in environmental risk assessment (ERA), was performed for simulated field releases of the dsxFCRISPRh transgene in West Africa.

METHODS

Building on consultative workshops in Africa that previously identified relevant environmental and health protection goals for ERA of gene drive in malaria vector control, 8 potentially harmful effects from these simulated releases were identified. These were stratified into 46 plausible pathways describing the causal chain of events that would be required for potential harms to occur. Risk hypotheses to interrogate critical steps in each pathway, and an analysis plan involving experiments, modelling and literature review to test each of those risk hypotheses, were developed.

RESULTS

Most potential harms involved increased human (n = 13) or animal (n = 13) disease transmission, emphasizing the importance to subsequent stages of ERA of data on vectorial capacity comparing transgenics to non-transgenics. Although some of the pathways (n = 14) were based on known anatomical alterations in dsxFCRISPRh homozygotes, many could also be applicable to field releases of a range of other transgenic strains of mosquito (n = 18). In addition to population suppression of target organisms being an accepted outcome for existing vector control programmes, these investigations also revealed that the efficacy of population suppression caused by the dsxFCRISPRh transgene should itself directly affect most pathways (n = 35).

CONCLUSIONS

Modelling will play an essential role in subsequent stages of ERA by clarifying the dynamics of this relationship between population suppression and reduction in exposure to specific potential harms. This analysis represents a comprehensive identification of plausible pathways to potential harm using problem formulation for a specific gene drive transgene and organism, and a transparent communication tool that could inform future regulatory studies, guide subsequent stages of ERA, and stimulate further, broader engagement on the use of population suppression gene drive to control malaria vectors in West Africa.

Comparative proteomics reveals mechanisms that underlie insecticide resistance in Culex pipiens pallens Coquillett

Mosquito control based on chemical insecticides is considered as an important element of the current global strategies for the control of mosquito-borne diseases. Unfortunately, the development of insecticide resistance of important vector mosquito species jeopardizes the effectiveness of insecticide-based mosquito control. In contrast to target site resistance, other mechanisms are far from being fully understood. Global protein profiles among cypermethrin-resistant, propoxur-resistant, dimethyl-dichloro-vinyl-phosphate-resistant and susceptible strain of Culex pipiens pallens were obtained and proteomic differences were evaluated by using isobaric tags for relative and absolute quantification labeling coupled with liquid chromatography/tandem mass spectrometric analysis. A susceptible strain of Culex pipiens pallens showed elevated resistance levels after 25 generations of insecticide selection, through iTRAQ data analysis detected 2,502 proteins, of which 1,513 were differentially expressed in insecticide-selected strains compared to the susceptible strain. Finally, midgut differential protein expression profiles were analyzed, and 62 proteins were selected for verification of differential expression using iTRAQ and parallel reaction monitoring strategy, respectively. iTRAQ profiles of adaptation selection to three insecticide strains combined with midgut profiles revealed that multiple insecticide resistance mechanisms operate simultaneously in resistant insects of Culex pipiens pallens. Significant molecular resources were developed for Culex pipiens pallens, potential candidates were involved in metabolic resistance and reducing penetration or sequestering insecticide. Future research that is targeted towards RNA interference of the identified metabolic targets, such as cuticular proteins, cytochrome P450s, glutathione S-transferases and ribosomal proteins proteins and biological pathways (drug metabolism—cytochrome P450, metabolism of xenobiotics by cytochrome P450, oxidative phosphorylation, ribosome) could lay the foundation for a better understanding of the genetic basis of insecticide resistance in Culex pipiens pallens.

Global protein profiles were compared among a susceptible strain of Cx. pipiens pallens and strains that were cypermethrin-resistant, propoxur-resistant, and dimethyl-dichloro-vinyl-phosphate-resistant after 25 generations of selection by distinct chemical insecticide families, multiple mechanisms were found to operate simultaneously in resistant mosquitoes of Cx. pipiens pallens, including mechanisms to lower penetration of or sequester the insecticide or to increase biodegradation of the insecticide via subtle alterations in either the cuticular protein levels or the activities of detoxification enzymes (P450s and glutathione S-transferases).

Ultra-conserved sequences in the genomes of highly diverse Anopheles mosquitoes, with implications for malaria vector control

DNA sequences that are exactly conserved over long evolutionary time scales have been observed in a variety of taxa. Such sequences are likely under strong functional constraint and they have been useful in the field of comparative genomics for identifying genome regions with regulatory function. A potential new application for these ultra-conserved elements (UCEs) has emerged in the development of gene drives to control mosquito populations. Many gene drives work by recognizing and inserting at a specific target sequence in the genome, often imposing a reproductive load as a consequence. They can therefore select for target sequence variants that provide resistance to the drive. Focusing on highly conserved, highly constrained sequences lowers the probability that variant, gene drive-resistant alleles can be tolerated. Here, we search for conserved sequences of 18 bp and over in an alignment of 21 Anopheles genomes, spanning an evolutionary timescale of 100 million years, and characterize the resulting sequences according to their location and function. Over 8000 UCEs were found across the alignment, with a maximum length of 164 bp. Length-corrected gene ontology analysis revealed that genes containing Anopheles UCEs were over-represented in categories with structural or nucleotide-binding functions. Known insect transcription factor binding sites were found in 48% of intergenic Anopheles UCEs. When we looked at the genome sequences of 1142 wild-caught mosquitoes, we found that 15% of the Anopheles UCEs contained no polymorphisms. Our list of Anopheles UCEs should provide a valuable starting point for the selection and testing of new targets for gene-drive modification in the mosquitoes that transmit malaria.

Comprehensive Quantitative Proteome Analysis of Aedes aegypti Identifies Proteins and Pathways Involved in Wolbachia pipientis and Zika Virus Interference Phenomenon

Zika virus (ZIKV) is a global public health emergency due to its association with microcephaly, Guillain-Barré syndrome, neuropathy, and myelitis in children and adults. A total of 87 countries have had evidence of autochthonous mosquito-borne transmission of ZIKV, distributed across four continents, and no antivirus therapy or vaccines are available. Therefore, several strategies have been developed to target the main mosquito vector, Aedes aegypti, to reduce the burden of different arboviruses. Among such strategies, the use of the maternally-inherited endosymbiont Wolbachia pipientis has been applied successfully to reduce virus susceptibility and decrease transmission. However, the mechanisms by which Wolbachia orchestrate resistance to ZIKV infection remain to be elucidated. In this study, we apply isobaric labeling quantitative mass spectrometry (MS)-based proteomics to quantify proteins and identify pathways altered during ZIKV infection; Wolbachia infection; co-infection with Wolbachia/ZIKV in the A. aegypti heads and salivary glands. We show that Wolbachia regulates proteins involved in reactive oxygen species production, regulates humoral immune response, and antioxidant production. The reduction of ZIKV polyprotein in the presence of Wolbachia in mosquitoes was determined by MS and corroborates the idea that Wolbachia helps to block ZIKV infections in A. aegypti. The present study offers a rich resource of data that may help to elucidate mechanisms by which Wolbachia orchestrate resistance to ZIKV infection in A. aegypti, and represents a step further on the development of new targeted methods to detect and quantify ZIKV and Wolbachia directly in complex tissues.

Chromosome-level genome assemblies of the malaria vectors Anopheles coluzzii and Anopheles arabiensis

BACKGROUND

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.

FINDINGS

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.

CONCLUSION

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.

Recent progress in the research of exosomes and Dscam regulated crab antiviral immunity

Crustaceans, including crab and shrimp, generally lack lymphocytes or adaptive immunity, and they rely solely on innate immunity for pathogen defense. The white spot syndrome virus (WSSV) causes the most prevalent viral disease in penaeid shrimps, which are widely cultured species in coastal waters worldwide. Numerous studies have elucidated the role of the immune system in protecting shrimps from WSSV infection for the development of safe and effective defensive strategies against WSSV. Although WSSV has a wide host range, it appears to exhibit high pathogenicity and virulence in only penaeid shrimps. Crabs are interesting models for studying immune responses after WSSV infection. Therefore, we reviewed recent information on the innate immune responses of crabs to WSSV and mainly focused on the antiviral functions of exosome-mediated apoptosis and alternatively spliced Down syndrome cell adhesion molecule. Our review may provide novel insights into antiviral management for crustaceans, especially penaeid shrimps.

A mosquito small RNA genomics resource reveals dynamic evolution and host responses to viruses and transposons

Although mosquitoes are major transmission vectors for pathogenic arboviruses, viral infection has little impact on mosquito health. This immunity is caused in part by mosquito RNA interference (RNAi) pathways that generate antiviral small interfering RNAs (siRNAs) and Piwi-interacting RNAs (piRNAs). RNAi also maintains genome integrity by potently repressing mosquito transposon activity in the germline and soma. However, viral and transposon small RNA regulatory pathways have not been systematically examined together in mosquitoes. Therefore, we developed an integrated mosquito small RNA genomics (MSRG) resource that analyzes the transposon and virus small RNA profiles in mosquito cell cultures and somatic and gonadal tissues across four medically important mosquito species. Our resource captures both somatic and gonadal small RNA expression profiles within mosquito cell cultures, and we report the evolutionary dynamics of a novel Mosquito-Conserved piRNA Cluster Locus (MCpiRCL) made up of satellite DNA repeats. In the larger culicine mosquito genomes we detected highly regular periodicity in piRNA biogenesis patterns coinciding with the expansion of Piwi pathway genes. Finally, our resource enables detection of cross talk between piRNA and siRNA populations in mosquito cells during a response to virus infection. The MSRG resource will aid efforts to dissect and combat the capacity of mosquitoes to tolerate and spread arboviruses.

De Novo Whole-Genome Assembly of the Swede Midge (Contarinia nasturtii), a Specialist of Brassicaceae, Using Linked-Read Sequencing

The swede midge, Contarinia nasturtii, is a cecidomyiid fly that feeds specifically on plants within the Brassicaceae. Plants in this family employ a glucosinolate-myrosinase defense system, which can be highly toxic to nonspecialist feeders. Feeding by C. nasturtii larvae induces gall formation, which can cause substantial yield losses thus making it a significant agricultural pest. A lack of genomic resources, in particular a reference genome, has limited deciphering the mechanisms underlying glucosinolate tolerance in C. nasturtii, which is of particular importance for managing this species. Here, we present an annotated, scaffolded reference genome of C. nasturtii using linked-read sequencing from a single individual and explore systems involved in glucosinolate detoxification. The C. nasturtii genome is similar in size and annotation completeness to that of the Hessian fly, Mayetiola destructor, but has greater contiguity. Several genes encoding enzymes involved in glucosinolate detoxification in other insect pests, including myrosinases, sulfatases, and glutathione S-transferases, were found, suggesting that C. nasturtii has developed similar strategies for feeding on Brassicaceae. The C. nasturtii genome will, therefore, be integral to continued research on plant-insect interactions in this system and contribute to effective pest management strategies.

Physical Mapping of the Anopheles (Nyssorhynchus) darlingi Genomic Scaffolds

The genome assembly of Anopheles darlingi consists of 2221 scaffolds (N50 = 115,072 bp) and has a size spanning 136.94 Mbp. This assembly represents one of the smallest genomes among Anopheles species. Anopheles darlingi genomic DNA fragments of ~37 Kb were cloned, end-sequenced, and used as probes for fluorescence in situ hybridization (FISH) with salivary gland polytene chromosomes. In total, we mapped nine DNA probes to scaffolds and autosomal arms. Comparative analysis of the An. darlingi scaffolds with homologous sequences of the Anopheles albimanus and Anopheles gambiae genomes identified chromosomal rearrangements among these species. Our results confirmed that physical mapping is a useful tool for anchoring genome assemblies to mosquito chromosomes.

The evolution and history of gene editing technologies

Scientific enquiry must be the driving force of research. This sentiment is manifested as the profound impact gene editing technologies are having in our current world. There exist three main gene editing technologies today: Zinc Finger Nucleases, TALENs and the CRISPR-Cas system. When these systems were being uncovered, none of the scientists set out to design tools to engineer genomes. They were simply trying to understand the mechanisms existing in nature. If it was not for this simple sense of wonder, we probably would not have these breakthrough technologies. In this chapter, we will discuss the history, applications and ethical issues surrounding these technologies, focusing on the now predominant CRISPR-Cas technology. Gene editing technologies, as we know them now, are poised to have an overwhelming impact on our world. However, it is impossible to predict the route they will take in the future or to comprehend the full impact of its repercussions.

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.

A Gene-Based Method for Cytogenetic Mapping of Repeat-Rich Mosquito Genomes

Long-read sequencing technologies have opened up new avenues of research on the mosquito genome biology, enabling scientists to better understand the remarkable abilities of vectors for transmitting pathogens. Although new genome mapping technologies such as Hi-C scaffolding and optical mapping may significantly improve the quality of genomes, only cytogenetic mapping, with the help of fluorescence in situ hybridization (FISH), connects genomic scaffolds to a particular chromosome and chromosome band. This mapping approach is important for creating and validating chromosome-scale genome assemblies for mosquitoes with repeat-rich genomes, which can potentially be misassembled. In this study, we describe a new gene-based physical mapping approach that was optimized using the newly assembled Aedes albopictus genome, which is enriched with transposable elements. To avoid amplification of the repetitive DNA, 15 protein-coding gene transcripts were used for the probe design. Instead of using genomic DNA, complementary DNA was utilized as a template for development of the PCR-amplified probes for FISH. All probes were successfully amplified and mapped to specific chromosome bands. The genome-unique probes allowed to perform unambiguous mapping of genomic scaffolds to chromosome regions. The method described in detail here can be used for physical genome mapping in other insects.

Transposable Elements and the Evolution of Insects

Insects are major contributors to our understanding of the interaction between transposable elements (TEs) and their hosts, owing to seminal discoveries, as well as to the growing number of sequenced insect genomes and population genomics and functional studies. Insect TE landscapes are highly variable both within and across insect orders, although phylogenetic relatedness appears to correlate with similarity in insect TE content. This correlation is unlikely to be solely due to inheritance of TEs from shared ancestors and may partly reflect preferential horizontal transfer of TEs between closely related species. The influence of insect traits on TE landscapes, however, remains unclear. Recent findings indicate that, in addition to being involved in insect adaptations and aging, TEs are seemingly at the cornerstone of insect antiviral immunity. Thus, TEs are emerging as essential insect symbionts that may have deleterious or beneficial consequences on their hosts, depending on context.

A population genomic unveiling of a new cryptic mosquito taxon within the malaria-transmitting Anopheles gambiae complex

The Anopheles gambiae complex consists of multiple morphologically indistinguishable mosquito species including the most important vectors of the malaria parasite Plasmodium falciparum in sub-Saharan Africa. Nine cryptic species have been described so far within the complex. The ecological, immunological and reproductive differences among these species will critically impact population responses to disease control strategies and environmental changes. Here, we examine whole-genome sequencing data from a longitudinal study of putative A. coluzzii in western Burkina Faso. Surprisingly, many specimens are genetically divergent from A. coluzzii and all other Anopheles species and represent a new taxon, here designated Anopheles TENGRELA (AT). Population genetic analysis suggests that the cryptic GOUNDRY subgroup, previously collected as larvae in central Burkina Faso, represents an admixed population descended from both A. coluzzii and AT. AT harbours low nucleotide diversity except for the 2La inversion polymorphism which is maintained by overdominance. It shows numerous fixed differences with A. coluzzii concentrated in several regions reflecting selective sweeps, but the two taxa are identical at standard diagnostic loci used for taxon identification, and thus, AT may often go unnoticed. We present an amplicon-based genotyping assay for identifying AT which could be usefully applied to numerous existing samples. Misidentified cryptic taxa could seriously confound ongoing studies of Anopheles ecology and evolution in western Africa, including phenotypic and genotypic surveys of insecticide resistance. Reproductive barriers between cryptic species may also complicate novel vector control efforts, for example gene drives, and hinder predictions about evolutionary dynamics of Anopheles and Plasmodium.

A Near-Chromosome Level Genome Assembly of Anopheles stephensi

Malaria remains a major healthcare risk to growing economies like India, and a chromosome-level reference genome of Anopheles stephensi is critical for successful vector management and understanding of vector evolution using comparative genomics. We report chromosome-level assemblies of an Indian strain, STE2, and a Pakistani strain SDA-500 by combining draft genomes of the two strains using a homology-based iterative approach. The resulting assembly IndV3/PakV3 with L50 of 9/12 and N50 6.3/6.9 Mb had scaffolds long enough for building 90% of the euchromatic regions of the three chromosomes, IndV3s/PakV3s, using low-resolution physical markers and enabled the generation of the next version of genome assemblies, IndV4/PakV4, using HiC data. We have validated these assemblies using contact maps against publicly available HiC raw data from two strains including STE2 and another lab strain of An. stephensi from UCI and compare the quality of the assemblies with other assemblies made available as preprints since the submission of the manuscript. We show that the IndV3s and IndV4 assemblies are sensitive in identifying a homozygous 2Rb inversion in the UCI strain and a 2Rb polymorphism in the STE2 strain. Multiple tandem copies of CYP6a14, 4c1, and 4c21 genes, implicated in insecticide resistance, lie within this inversion locus. Comparison of assembled genomes suggests a variation of 1 in 81 positions between the UCI and STE2 lab strains, 1 in 82 between SDA-500 and UCI strain, and 1 in 113 between SDA-500 and STE2 strains of An. stephensi, which are closer than 1 in 68 variations among individuals from two other lab strains sequenced and reported here. Based on the developmental transcriptome and orthology of all the 54 olfactory receptors (ORs) to those of other Anopheles species, we identify an OR with the potential for host recognition in the genus Anopheles. A comparative analysis of An. stephensi genomes with the completed genomes of a few other Anopheles species suggests limited inter-chromosomal gene flow and loss of synteny within chromosomal arms even among the closely related species.

Mosquito genomes are frequently invaded by transposable elements through horizontal transfer

Transposable elements (TEs) are mobile genetic elements that parasitize basically all eukaryotic species genomes. Due to their complexity, an in-depth TE characterization is only available for a handful of model organisms. In the present study, we performed a de novo and homology-based characterization of TEs in the genomes of 24 mosquito species and investigated their mode of inheritance. More than 40% of the genome of Aedes aegypti, Aedes albopictus, and Culex quinquefasciatus is composed of TEs, while it varied substantially among Anopheles species (0.13%-19.55%). Class I TEs are the most abundant among mosquitoes and at least 24 TE superfamilies were found. Interestingly, TEs have been extensively exchanged by horizontal transfer (172 TE families of 16 different superfamilies) among mosquitoes in the last 30 million years. Horizontally transferred TEs represents around 7% of the genome in Aedes species and a small fraction in Anopheles genomes. Most of these horizontally transferred TEs are from the three ubiquitous LTR superfamilies: Gypsy, Bel-Pao and Copia. Searching more than 32,000 genomes, we also uncovered transfers between mosquitoes and two different Phyla-Cnidaria and Nematoda-and two subphyla-Chelicerata and Crustacea, identifying a vector, the worm Wuchereria bancrofti, that enabled the horizontal spread of a Tc1-mariner element among various Anopheles species. These data also allowed us to reconstruct the horizontal transfer network of this TE involving more than 40 species. In summary, our results suggest that TEs are frequently exchanged by horizontal transfers among mosquitoes, influencing mosquito's genome size and variability.

Genome streamlining in a minute herbivore that manipulates its host plant

The tomato russet mite, Aculops lycopersici, is among the smallest animals on earth. It is a worldwide pest on tomato and can potently suppress the host's natural resistance. We sequenced its genome, the first of an eriophyoid, and explored whether there are genomic features associated with the mite's minute size and lifestyle. At only 32.5 Mb, the genome is the smallest yet reported for any arthropod and, reminiscent of microbial eukaryotes, exceptionally streamlined. It has few transposable elements, tiny intergenic regions, and is remarkably intron-poor, as more than 80% of coding genes are intronless. Furthermore, in accordance with ecological specialization theory, this defense-suppressing herbivore has extremely reduced environmental response gene families such as those involved in chemoreception and detoxification. Other losses associate with this species' highly derived body plan. Our findings accelerate the understanding of evolutionary forces underpinning metazoan life at the limits of small physical and genome size.

Expression, Purification, Refolding, and Characterization of a Neverland Protein From Caenorhabditis elegans

Steroid hormones that serve as vital compounds are necessary for the development and metabolism of a variety of organisms. The neverland (NVD) family genes encode the conserved Rieske-type oxygenases, which are accountable for the dehydrogenation during the synthesis and regulation of steroid hormones. However, the His-tagged NVD protein from Caenorhabditis elegans expresses as inclusion bodies in Escherichia coli BL21 (DE3). This bottleneck can be solved through refolding by urea or the introduction of a maltose-binding protein (MBP) tag at the N-terminus. Through further research on purification after the introduction of a MBP tag at the N-terminus, the CD measurement and fluorescence-based thermal shift assay indicated that MBP was favorable for the NVD proteins' solubility and stability, which may be beneficial for the large-scale manufacture of NVD protein for further research. The structural model contained the Rieske [2Fe-2S] domain and non-heme iron-binding motif, which were similar to 3-ketosteroid 9 α-hydroxylase.

Genome variation and population structure among 1142 mosquitoes of the African malaria vector species Anopheles gambiae and Anopheles coluzzii

Mosquito control remains a central pillar of efforts to reduce malaria burden in sub-Saharan Africa. However, insecticide resistance is entrenched in malaria vector populations, and countries with a high malaria burden face a daunting challenge to sustain malaria control with a limited set of surveillance and intervention tools. Here we report on the second phase of a project to build an open resource of high-quality data on genome variation among natural populations of the major African malaria vector species Anopheles gambiae and Anopheles coluzzii We analyzed whole genomes of 1142 individual mosquitoes sampled from the wild in 13 African countries, as well as a further 234 individuals comprising parents and progeny of 11 laboratory crosses. The data resource includes high-confidence single-nucleotide polymorphism (SNP) calls at 57 million variable sites, genome-wide copy number variation (CNV) calls, and haplotypes phased at biallelic SNPs. We use these data to analyze genetic population structure and characterize genetic diversity within and between populations. We illustrate the utility of these data by investigating species differences in isolation by distance, genetic variation within proposed gene drive target sequences, and patterns of resistance to pyrethroid insecticides. This data resource provides a foundation for developing new operational systems for molecular surveillance and for accelerating research and development of new vector control tools. It also provides a unique resource for the study of population genomics and evolutionary biology in eukaryotic species with high levels of genetic diversity under strong anthropogenic evolutionary pressures.

Proteases and pseudoproteases in parasitic arthropods of clinical importance

Parasitic arthropods feed on blood or skin tissue and share comparable repertoires of proteases involved in haematophagy, digestion, egg development and immunity. While proteolytically active proteases of multiple classes dominate, an increasing number of pseudoproteases have been discovered that have no proteolytic function but are pharmacologically active biomolecules, evolved to carry out alternative functions as regulatory, antihaemostatic, anti-inflammatory or immunomodulatory compounds. In this review, we provide an overview of proteases and pseudoproteases from clinically important arthropod parasites. Many of these act in central biological pathways of parasite survival and host-parasite interaction and may be potential targets for therapeutic interventions.

Large-Scale Comparative Analyses of Tick Genomes Elucidate Their Genetic Diversity and Vector Capacities

Among arthropod vectors, ticks transmit the most diverse human and animal pathogens, leading to an increasing number of new challenges worldwide. Here we sequenced and assembled high-quality genomes of six ixodid tick species and further resequenced 678 tick specimens to understand three key aspects of ticks: genetic diversity, population structure, and pathogen distribution. We explored the genetic basis common to ticks, including heme and hemoglobin digestion, iron metabolism, and reactive oxygen species, and unveiled for the first time that genetic structure and pathogen composition in different tick species are mainly shaped by ecological and geographic factors. We further identified species-specific determinants associated with different host ranges, life cycles, and distributions. The findings of this study are an invaluable resource for research and control of ticks and tick-borne diseases.

Ultraconserved Non-coding DNA Within Diptera and Hymenoptera

This study has taken advantage of the availability of the assembled genomic sequence of flies, mosquitos, ants and bees to explore the presence of ultraconserved sequence elements in these phylogenetic groups. We compared non-coding sequences found within and flanking Drosophila developmental genes to homologous sequences in Ceratitis capitata and Musca domestica. Many of the conserved sequence blocks (CSBs) that constitute Drosophila cis-regulatory DNA, recognized by EvoPrinter alignment protocols, are also conserved in Ceratitis and Musca. Also conserved is the position but not necessarily the orientation of many of these ultraconserved CSBs (uCSBs) with respect to flanking genes. Using the mosquito EvoPrint algorithm, we have also identified uCSBs shared among distantly related mosquito species. Side by side comparison of bee and ant EvoPrints of selected developmental genes identify uCSBs shared between these two Hymenoptera, as well as less conserved CSBs in either one or the other taxon but not in both. Analysis of uCSBs in these dipterans and Hymenoptera will lead to a greater understanding of their evolutionary origin and function of their conserved non-coding sequences and aid in discovery of core elements of enhancers.

This study applies the phylogenetic footprinting program EvoPrinter to detection of ultraconserved non-coding sequence elements in Diptera, including flies and mosquitos, and Hymenoptera, including ants and bees. EvoPrinter outputs an interspecies comparison as a single sequence in terms of the input reference sequence. Ultraconserved sequences flanking known developmental genes were detected in Ceratitis and Musca when compared with Drosophila species, in Aedes and Culex when compared with Anopheles, and between ants and bees. Our methods are useful in detecting and understanding the core evolutionarily hardened sequences required for gene regulation.

HiCanu: accurate assembly of segmental duplications, satellites, and allelic variants from high-fidelity long reads

Complete and accurate genome assemblies form the basis of most downstream genomic analyses and are of critical importance. Recent genome assembly projects have relied on a combination of noisy long-read sequencing and accurate short-read sequencing, with the former offering greater assembly continuity and the latter providing higher consensus accuracy. The recently introduced Pacific Biosciences (PacBio) HiFi sequencing technology bridges this divide by delivering long reads (>10 kbp) with high per-base accuracy (>99.9%). Here we present HiCanu, a modification of the Canu assembler designed to leverage the full potential of HiFi reads via homopolymer compression, overlap-based error correction, and aggressive false overlap filtering. We benchmark HiCanu with a focus on the recovery of haplotype diversity, major histocompatibility complex (MHC) variants, satellite DNAs, and segmental duplications. For diploid human genomes sequenced to 30× HiFi coverage, HiCanu achieved superior accuracy and allele recovery compared to the current state of the art. On the effectively haploid CHM13 human cell line, HiCanu achieved an NG50 contig size of 77 Mbp with a per-base consensus accuracy of 99.999% (QV50), surpassing recent assemblies of high-coverage, ultralong Oxford Nanopore Technologies (ONT) reads in terms of both accuracy and continuity. This HiCanu assembly correctly resolves 337 out of 341 validation BACs sampled from known segmental duplications and provides the first preliminary assemblies of nine complete human centromeric regions. Although gaps and errors still remain within the most challenging regions of the genome, these results represent a significant advance toward the complete assembly of human genomes.

Reference genes for proximal femoral epiphysiolysis expression studies in broilers cartilage

The use of reference genes is required for relative quantification in gene expression analysis and the stability of these genes can be variable depending on the experimental design. Therefore, it is indispensable to test the reliability of endogenous genes previously to their use. This study evaluated nine candidate reference genes to select the most stable genes to be used as reference in gene expression studies with the femoral cartilage of normal and epiphysiolysis-affected broilers. The femur articular cartilage of 29 male broilers with 35 days of age was collected, frozen and further submitted to RNA extraction and quantitative PCR (qPCR) analysis. The candidate reference genes evaluated were GAPDH, HMBS, HPRT1, MRPS27, MRPS30, RPL30, RPL4, RPL5, and RPLP1. For the gene stability evaluation, three software were used: GeNorm, BestKeeper and NormFinder, and a global ranking was generated using the function RankAggreg. In this study, the RPLP1 and RPL5 were the most reliable endogenous genes being recommended for expression studies with femur cartilage in broilers with epiphysiolysis and possible other femur anomalies.

The Genetic Architecture of Post-Zygotic Reproductive Isolation Between Anopheles coluzzii and An. quadriannulatus

The Anopheles gambiae complex is comprised of eight morphologically indistinguishable species and has emerged as a model system for the study of speciation genetics due to the rapid radiation of its member species over the past two million years. Male hybrids between most An. gambiae complex species pairs are sterile, and some genotype combinations in hybrid males cause inviability. We investigated the genetic basis of hybrid male inviability and sterility between An. coluzzii and An. quadriannulatus by measuring segregation distortion and performing a QTL analysis of sterility in a backcross population. Hybrid males were inviable if they inherited the An. coluzzii X chromosome and were homozygous at one or more loci in 18.9 Mb region of chromosome 3. The An. coluzzii X chromosome has a disproportionately large effect on hybrid sterility when introgressed into an An. quadriannulatus genetic background. Additionally, an epistatic interaction between the An. coluzzii X and a 1.12 Mb, pericentric region of the An. quadriannulatus 3L chromosome arm has a statistically significant contribution to the hybrid sterility phenotype. This same epistatic interaction occurs when the An. coluzzii X is introgressed into the genetic background of An. arabiensis, the sister species of An. quadriannulatus, suggesting that this may represent one of the first Dobzhansky-Muller incompatibilities to evolve early in the radiation of the Anopheles gambiae species complex. We describe the additive effects of each sterility QTL, epistatic interactions between them, and genes within QTL with protein functions related to mating behavior, reproduction, spermatogenesis, and microtubule morphogenesis, whose divergence may contribute to post-zygotic reproductive isolation between An. coluzzii and An. quadriannulatus.

Recent advances and future perspectives in vector-omics

We have reviewed recent progress and the remaining challenges in vector-omics. We have highlighted several technologies and applications that facilitate novel biological insights beyond achieving a reference-quality genome assembly. Among other topics, we have discussed the applications of chromatin conformation capture, chromatin accessibility assays, optical mapping, full-length RNA sequencing, single cell RNA analysis, proteomics, and population genomics. We anticipate that we will witness a great expansion in vector-omics research not only in its application in a broad range of species, but also its ability to uncover novel genetic elements and tackle previously inaccessible regions of the genome. It is our hope that the continued innovation in device portability, cost reduction, and informatics support will in the foreseeable future bring vector-omics to every vector laboratory and field station in the world, which will have an unparalleled impact on basic research and the control of vector-borne infectious diseases.

TickSialoFam (TSFam): A Database That Helps to Classify Tick Salivary Proteins, a Review on Tick Salivary Protein Function and Evolution, With Considerations on the Tick Sialome Switching Phenomenon

Tick saliva contains a complex mixture of peptides and non-peptides that counteract their hosts' hemostasis, immunity, and tissue-repair reactions. Recent transcriptomic studies have revealed over one thousand different transcripts coding for secreted polypeptides in a single tick species. Not only do these gene products belong to many expanded families, such as the lipocalins, metalloproteases, Antigen-5, cystatins, and apyrases, but also families that are found exclusively in ticks, such as the evasins, Isac, DAP36, and many others. Phylogenetic analysis of the deduced protein sequences indicate that the salivary genes exhibit an increased rate of evolution due to a lower evolutionary constraint and/or positive selection, allowing for a large diversity of tick salivary proteins. Thus, for each new tick species that has its salivary transcriptome sequenced and assembled, a formidable task of annotation of these transcripts awaits. Currently, as of November 2019, there are over 287 thousand coding sequences deposited at the National Center for Biotechnology Information (NCBI) that are derived from tick salivary gland mRNA. Here, from these 287 thousand sequences we identified 45,264 potential secretory proteins which possess a signal peptide and no transmembrane domains on the mature peptide. By using the psiblast tools, position-specific matrices were constructed and assembled into the TickSialoFam (TSF) database. The TSF is a rpsblastable database that can help with the annotation of tick sialotranscriptomes. The TSA database identified 136 tick salivary secreted protein families, as well as 80 families of endosomal-related products, mostly having a protein modification function. As the number of sequences increases, and new annotation details become available, new releases of the TSF database may become available.

Genetic engineering and bacterial pathogenesis against the vectorial capacity of mosquitoes

Mosquitoes are the main vector of multiple diseases worldwide and transmit viral (malaria, chikungunya, encephalitis, yellow fever, as well as dengue fever), as well as bacterial diseases (tularemia). To manage the outbreak of mosquito populations, various management programs include the application of chemicals, followed by biological and genetic control. Here we aimed to focus on the role of bacterial pathogenesis and molecular tactics for the management of mosquitoes and their vectorial capacity. Bacterial pathogenesis and molecular manipulations have a substantial impact on the biology of mosquitoes, and both strategies change the gene expression and regulation of disease vectors. The strategy for genetic modification is also proved to be excellent for the management of mosquitoes, which halt the development of population via incompatibility of different sex. Therefore, the purpose of the present discussion is to illustrate the impact of both approaches against the vectorial capacity of mosquitoes. Moreover, it could be helpful to understand the relationship of insect-pathogen and to manage various insect vectors as well as diseases.

Molecular aspects of lipid metabolism in the midgut gland of the brown shrimp Crangon crangon

The brown shrimp, Crangon crangon, is well adapted to the variable environmental conditions in the southern North Sea. It is very abundant, has high reproduction rates, and holds a key position in coastal ecosystems. This species has very low lipid deposits in the midgut gland, suggesting that the main function of the midgut gland is metabolic turnover rather than energy storage. Based on seasonal gene expression studies and established transcriptome data, we investigated key components of lipid metabolic pathways. Gene expression of triacylglycerol lipase, phospholipase, and fatty acid desaturase were analyzed and compared with that of other digestive enzymes involved in lipid, carbohydrate, and protein catabolism. Our results suggest that gene expression of digestive enzymes involved in lipid metabolism is modulated by the lipid content in the midgut gland and is related to food availability. Brown shrimp seem to be capable of using cellular phospholipids during periods of food paucity but high energetic (lipid) requirements. Two of three isoforms of fatty acid binding proteins (FABPs) from the midgut gland involved in fatty acid transport showed specific mutations of the binding site. We hypothesize that the mutations in FABPs and deficiencies in anabolic pathways limit lipid storage capacities in the midgut gland of C. crangon. In turn, food utilization, including lipid catabolism, has to be efficient to fulfill the energetic requirements of brown shrimp.

Delivery Methods for RNAi in Mosquito Larvae

Mosquito-transmitted diseases pose a threat for a great portion of the world population. Chemical insecticides are the main tool for mosquito control. Heavy dependence on chemicals created several problems such as resistance development in many mosquito species, environmental effects, and human health issues. Other tools for mosquito control were developed and used in some parts of the world. Ribonucleic acid interference (RNAi) is a reverse genetic mechanism that was recently introduced as a new tool for pest control. Regarding mosquito, RNAi was used to study gene function and to discover genes that can be used as targets for control purposes. Several delivery methods are used to induce RNAi in mosquito larvae. Some methods such as injection and soaking are used routinely in RNAi research but have no application in the field. Other methods such as nanoparticles and microbes have some characteristics that make them good candidates for field application. In this report, we will focus on delivery methods for RNAi in mosquito larvae and will give examples for each method.

Genome-wide identification and analysis of nuclear receptors genes for lethal screening against Bemisia tabaci Q

BACKGROUND

Nuclear receptors (NRs) play an essential role in diverse biological processes, such as insect metamorphosis. Here, transcriptome analysis and functional studies were used to determine whether NRs are involved in metamorphosis of whitefly Bemisia tabaci Q, a serious pest to crops, and to find some potential insecticide targets.

RESULTS

Twenty NRs were identified in the Bemisia tabaci Q genome and categorized into the NR0-NR6 subfamilies. The phylogenetic tree of NRs from Bemisia tabaci Q and other representative species was constructed, which provided evolutionary insight into their genetic distances. The results of spatiotemporal gene expression indicated that the majority of NR gene expression was higher in the head than the abdomen and higher in eggs than adults. Further functional analysis using RNA interference (RNAi) showed that NR genes play an important role in Bemisia tabaci Q pupation and eclosion. With respect to high mortality and effects on growth, this was reflected in the unable to become pupa when the third-stage nymph treated with double-stranded RNA (dsRNA) and the developmental time delay (4-7 days) when pupae were treated with dsRNA for the 12 NR genes during molting compared with the development time in the control.

CONCLUSION

This study provides insight into NR functions during the metamorphosis stages of Bemisia tabaci Q. Several candidate genes could be potential insecticide targets for whitefly pest control due to their important roles in insect development. © 2020 Society of Chemical Industry.

Characterization of an adulticidal and larvicidal interfering RNA pesticide that targets a conserved sequence in mosquito G protein-coupled dopamine 1 receptor genes

G protein-coupled receptors (GPCRs), key regulators of a variety of critical biological processes, are attractive targets for insecticide development. Given the importance of these receptors in many organisms, including humans, it is critical that novel pesticides directed against GPCRs are designed to be species-specific. Here, we present characterization of an interfering RNA pesticide (IRP) targeting the mosquito GPCR-encoding dopamine 1 receptor (dop1) genes. A small interfering RNA corresponding to dop1 was identified in a screen for IRPs that kill Aedes aegypti during both the adult and larval stages. The 25 bp sequence targeted by this IRP is conserved in the dop1 genes of multiple mosquito species, but not in non-target organisms, indicating that it could function as a biorational mosquito insecticide. Aedes aegypti adults treated through microinjection or attractive toxic sugar bait delivery of small interfering RNA corresponding to the target site exhibited severe neural and behavioral defects and high levels of adult mortality. Likewise, A. aegypti larval consumption of dried inactivated yeast tablets prepared from a Saccharomyces cerevisiae strain engineered to express short hairpin RNA corresponding to the dop1 target site resulted in severe neural defects and larval mortality. Aedes albopictus and Anopheles gambiae adult and larval mortality was also observed following treatment with dop1 IRPs, which were not toxic to non-target arthropods. The results of this investigation indicate that dop1 IRPs can be used for species-specific targeting of dop1 GPCRs and may represent a new biorational strategy for control of both adult and larval mosquitoes.

Chromosome-level analysis of the Crassostrea hongkongensis genome reveals extensive duplication of immune-related genes in bivalves

Crassostrea hongkongensis is a popular and important native oyster species that is cultured mainly along the coast of the South China Sea. However, the absence of a reference genome has restricted genetic studies and the development of molecular breeding schemes for this species. Here, we combined PacBio and 10 × Genomics technologies to create a C. hongkongensis genome assembly, which has a size of 610 Mb, and is close to that estimated by flow cytometry (~650 Mb). Contig and scaffold N50 are 2.57 and 4.99 Mb, respectively, and BUSCO analysis indicates that 95.8% of metazoan conserved genes are completely represented. Using a high-density linkage map of its closest related species, C. gigas, a total of 521 Mb (85.4%) was anchored to 10 haploid chromosomes. Comparative genomic analyses with other molluscs reveal that several immune- or stress response-related genes extensively expanded in bivalves by tandem duplication, including C1q, Toll-like receptors and Hsp70, which are associated with their adaptation to filter-feeding and sessile lifestyles in shallow sea and/or deep-sea ecosystems. Through transcriptome sequencing, potential genes and pathways related to sex determination and gonad development were identified. The genome and transcriptome of C. hongkongensis provide valuable resources for future molecular studies, genetic improvement and genome-assisted breeding of oysters.

De novo assembly of the olive fruit fly (Bactrocera oleae) genome with linked-reads and long-read technologies minimizes gaps and provides exceptional Y chromosome assembly

BACKGROUND

The olive fruit fly, Bactrocera oleae, is the most important pest in the olive fruit agribusiness industry. This is because female flies lay their eggs in the unripe fruits and upon hatching the larvae feed on the fruits thus destroying them. The lack of a high-quality genome and other genomic and transcriptomic data has hindered progress in understanding the fly's biology and proposing alternative control methods to pesticide use.

RESULTS

Genomic DNA was sequenced from male and female Demokritos strain flies, maintained in the laboratory for over 45 years. We used short-, mate-pair-, and long-read sequencing technologies to generate a combined male-female genome assembly (GenBank accession GCA_001188975.2). Genomic DNA sequencing from male insects using 10x Genomics linked-reads technology followed by mate-pair and long-read scaffolding and gap-closing generated a highly contiguous 489 Mb genome with a scaffold N50 of 4.69 Mb and L50 of 30 scaffolds (GenBank accession GCA_001188975.4). RNA-seq data generated from 12 tissues and/or developmental stages allowed for genome annotation. Short reads from both males and females and the chromosome quotient method enabled identification of Y-chromosome scaffolds which were extensively validated by PCR.

CONCLUSIONS

The high-quality genome generated represents a critical tool in olive fruit fly research. We provide an extensive RNA-seq data set, and genome annotation, critical towards gaining an insight into the biology of the olive fruit fly. In addition, elucidation of Y-chromosome sequences will advance our understanding of the Y-chromosome's organization, function and evolution and is poised to provide avenues for sterile insect technique approaches.

Structural Variation of the X Chromosome Heterochromatin in the Anopheles gambiae Complex

Heterochromatin is identified as a potential factor driving diversification of species. To understand the magnitude of heterochromatin variation within the Anopheles gambiae complex of malaria mosquitoes, we analyzed metaphase chromosomes in An. arabiensis, An. coluzzii, An. gambiae, An. merus, and An. quadriannulatus. Using fluorescence in situ hybridization (FISH) with ribosomal DNA (rDNA), a highly repetitive fraction of DNA, and heterochromatic Bacterial Artificial Chromosome (BAC) clones, we established the correspondence of pericentric heterochromatin between the metaphase and polytene X chromosomes of An. gambiae. We then developed chromosome idiograms and demonstrated that the X chromosomes exhibit qualitative differences in their pattern of heterochromatic bands and position of satellite DNA (satDNA) repeats among the sibling species with postzygotic isolation, An. arabiensis, An. merus, An. quadriannulatus, and An. coluzzii or An. gambiae. The identified differences in the size and structure of the X chromosome heterochromatin point to a possible role of repetitive DNA in speciation of mosquitoes. We found that An. coluzzii and An. gambiae, incipient species with prezygotic isolation, share variations in the relative positions of the satDNA repeats and the proximal heterochromatin band on the X chromosomes. This previously unknown genetic polymorphism in malaria mosquitoes may be caused by a differential amplification of DNA repeats or an inversion in the sex chromosome heterochromatin.

Identification and expression profiling of neuropeptides and neuropeptide receptor genes in Atrijuglans hetaohei

Atrijuglans hetaohei Yang (Lepidoptera: Gelechioidea), is one of the major pests that can seriously damage the walnut fruits. Neuropeptides and their receptors regulate most physiological functions in insects and represent new targets for the development of control agents. To identify the neuropeptides and their receptors from A. hetaohei, we sequenced and analyzed its head transcriptomic data, identified 32 neuropeptides and 39 neuropeptide receptor genes. Sequence comparisons and phylogenetic analyses suggest that A. hetaohei neuropeptides and receptor genes have high homology with those in Bombyx mori, Chilo suppressalis, Plutella xylostella and Helicoverpa armigera. Moreover, gene expression patterns revealed that neuropeptide genes such as AKH1, CP, MS and PTTH were expressed specifically in male head, while CAP3, DH, NPLP1, PBAN and SIF showed higher expression in the female head. Bur showed abdomen biased expression in both male and female. Neuropeptide receptor genes such as A8, A11, A15 and LGR were highly expressed in male head, whereas A24 and LGR2 were preferentially expressed in female head. This is the first sequencing, identification and expression analyses of neuropeptides and neuropeptide receptor genes from A. hetaohei. Our results could provide a powerful background that will facilitate the further investigations using transcriptomics to determine neuropeptides and their receptors presence, functions, and indicates potential targets in A. hetaohei for a novel pest management strategy.

Linked-read sequencing identifies abundant microinversions and introgression in the arboviral vector Aedes aegypti

Background

Aedes aegypti is the principal mosquito vector of Zika, dengue, and yellow fever viruses. Two subspecies of Ae. aegypti exhibit phenotypic divergence with regard to habitat, host preference, and vectorial capacity. Chromosomal inversions have been shown to play a major role in adaptation and speciation in dipteran insects and would be of great utility for studies of Ae. aegypti. However, the large and highly repetitive genome of Ae. aegypti makes it difficult to detect inversions with paired-end short-read sequencing data, and polytene chromosome analysis does not provide sufficient resolution to detect chromosome banding patterns indicative of inversions.

Results

To characterize chromosomal diversity in this species, we have carried out deep Illumina sequencing of linked-read (10X Genomics) libraries in order to discover inversion loci as well as SNPs. We analyzed individuals from colonies representing the geographic limits of each subspecies, one contact zone between subspecies, and a closely related sister species. Despite genome-wide SNP divergence and abundant microinversions, we do not find any inversions occurring as fixed differences between subspecies. Many microinversions are found in regions that have introgressed and have captured genes that could impact behavior, such as a cluster of odorant-binding proteins that may play a role in host feeding preference.

Conclusions

Our study shows that inversions are abundant and widely shared among subspecies of Aedes aegypti and that introgression has occurred in regions of secondary contact. This library of 32 novel chromosomal inversions demonstrates the capacity for linked-read sequencing to identify previously intractable genomic rearrangements and provides a foundation for future population genetics studies in this species.

Diverse patterns of constitutive and inducible overexpression of detoxifying enzyme genes among resistant Aphis glycines populations

Understanding the mechanisms of pyrethroid resistance is essential to the effective management of pesticide resistance in Aphis glycines Matsumura (Hemiptera: Aphididae). We mined putative detoxifying enzyme genes in the draft genome sequence of A. glycines for cytochrome oxidase P450 (CYP), glutathione-S-transferase (GST) and esterases (E4 and carboxylesterases-CES). Aphids from clonal populations resistant to pyrethroids from three sites in Minnesota, USA, were screened against a diagnostic LC99 concentration of either λ-cyhalothrin or bifenthrin and detoxifying enzyme genes expression in survivors was analyzed by qPCR. Their expression profiles were compared relative to a susceptible clonal population. We found 61 CYP (40 full-length), seven GST (all full-length), seven E4 (five full-length) and three CES (two full-length) genes, including 24 possible pseudogenes. The detoxifying enzymes had different expression profiles across resistant aphid populations, possibly reflecting differences in the genetic background and pyrethroid selection pressures as the number of constitutively overexpressed detoxifying enzyme genes was correlated with the level of resistance. Our findings will strengthen the understanding of the pyrethroid resistance mechanisms in A. glycines.

Dynamics of prevalence and distribution pattern of avian Plasmodium species and its vectors in diverse zoogeographical areas - A review

Avian Plasmodium is of special interest to health care scientists and veterinarians due to the potency of causing avian malaria in non-adapted birds and their evolutionary phylogenetic relationship with human malaria species. This article aimed to provide a comprehensive list of the common avian Plasmodium parasites in the birds and mosquitoes, to specify the common Plasmodium species and lineages in the selected regions of West of Asia, East of Europe, and North of Africa/Middle East, and to determine the contribution of generalist and host-specific Plasmodium species and lineages. The final list of published infected birds includes 146 species, among which Passer domesticus was the most prevalent in the studied areas. The species of Acrocephalus arundinaceus and Sylvia atricapilla were reported as common infected hosts in the examined regions of three continents. The highest numbers of common species of infected birds between continent pairs were from Asia and Europe, and no common record was found from Europe and Africa. The species of Milvus migrans and Upupa epops were recorded as common species from Asia and Africa. The lineage of GRW11 and species of P. relictum were the most prevalent parasites among all the infection records in birds. The most prevalent genus of vectors of avian malaria belonged to Culex and species of Cx. pipiens. The lineage SGS1 with the highest number of occurrence has been found in various vectors comprising Cx. pipiens, Cx. modestus, Cx. theileri, Cx. sasai, Cx. perexiguus, Lutzia vorax, and Culicoides alazanicus. A total of 31 Plasmodium species and 59 Plasmodium lineages were recorded from these regions. SGS1, GRW04, and GRW11, and P. relictum and P. vaughani are specified as common generalist avian malaria parasites from these three geographic areas. The presence of avian Plasmodium parasites in distant geographic areas and various hosts may be explained by the movement of the infected birds through the migration routes. Although most recorded lineages were from Asia, investigating the distribution of lineages in some of the countries has not been done. Thus, the most important outcome of this review is the determination of the distribution pattern of parasite and vector species that shed light on gaps requiring further studies on the monitoring of avian Plasmodium and common vectors extension. This task could be achieved through scientific field and laboratory networking, performing active surveillance and designing regional/continental control programs of birds' malaria and other zoonotic diseases.

Three-dimensional Organization of Polytene Chromosomes in Somatic and Germline Tissues of Malaria Mosquitoes

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.

Assessing connectivity despite high diversity in island populations of a malaria mosquito

Documenting isolation is notoriously difficult for species with vast polymorphic populations. High proportions of shared variation impede estimation of connectivity, even despite leveraging information from many genetic markers. We overcome these impediments by combining classical analysis of neutral variation with assays of the structure of selected variation, demonstrated using populations of the principal African malaria vector Anopheles gambiae. Accurate estimation of mosquito migration is crucial for efforts to combat malaria. Modeling and cage experiments suggest that mosquito gene drive systems will enable malaria eradication, but establishing safety and efficacy requires identification of isolated populations in which to conduct field testing. We assess Lake Victoria islands as candidate sites, finding one island 30 km offshore is as differentiated from mainland samples as populations from across the continent. Collectively, our results suggest sufficient contemporary isolation of these islands to warrant consideration as field-testing locations and illustrate shared adaptive variation as a useful proxy for connectivity in highly polymorphic species.

Transcriptome Sequencing and Analysis of Changes Associated with Insecticide Resistance in the Dengue Mosquito (Aedes aegypti) in Vietnam

The mosquito Aedes aegypti is a transmission vector for dangerous epidemic diseases in humans. Insecticides have been used as the most general vector control method in the world. However, Ae. aegypti have developed many resistant mechanisms such as reduced neuronal sensitivity to insecticides (target-site resistance), enhanced insecticide metabolism (metabolic resistance), altered transport, sequestration, and other mechanisms. It has become a major problem for vector control programs. Transcriptome sequencing and bioinformatic analysis were used to compare transcription levels between a susceptible strain (Bora7) and a resistant strain (KhanhHoa7) collected from the field. A total of 161 million Illumina reads, including 66,076,678 reads from the Bora7 strain and 69,606,654 reads from the KhanhHoa7 strain, were generated and assembled into 11,174 genes. A comparison of the KhanhHoa7 transcriptome to that of Bora7 showed 672 upregulated genes and 488 downregulated genes. We identified the highly upregulated genes: cytochrome P450 4C1, 4C3, 4C21, 4D1, 4D1 isoform X2, 4D2, 4D2 isoform X2, 4G15, 6A2, 6A8, 6D3, and 9E2; Glutathione S transferase (GST1), UGT1-3, 1-7, 2B15, and 2B37; binding cassette transporter (ABC) transporter F family member 4 and ABC transporter G family member 20. Interestingly, there was a significant increase in the expression of the genes such as CYP9E2 (8.3-fold), CYP6A8 (5.9-fold), CYP6D3 (5.4-fold), CYP4C21 (5.4-fold), CYP4G15 (5.2-fold), GST1 (3.5-fold), and ABC transporter 4 (2.1-fold). Our results suggested a potential relationship between the expression of the genes in metabolic processes and insecticide resistance in the studied strain. These results may contribute to the understanding of the mechanisms of insecticide resistance in Ae. aegypti.

A New Assessment of Thioester-Containing Proteins Diversity of the Freshwater Snail Biomphalaria glabrata

Thioester-containing proteins (TEPs) superfamily is known to play important innate immune functions in a wide range of animal phyla. TEPs are involved in recognition, and in the direct or mediated killing of several invading organisms or pathogens. While several TEPs have been identified in many invertebrates, only one TEP (named BgTEP) has been previously characterized in the freshwater snail, Biomphalaria glabrata. As the presence of a single member of that family is particularly intriguing, transcriptomic data and the recently published genome were used to explore the presence of other BgTEP related genes in B. glabrata. Ten other TEP members have been reported and classified into different subfamilies: Three complement-like factors (BgC3-1 to BgC3-3), one α-2-macroblobulin (BgA2M), two macroglobulin complement-related proteins (BgMCR1, BgMCR2), one CD109 (BgCD109), and three insect TEP (BgTEP2 to BgTEP4) in addition to the previously characterized BgTEP that we renamed BgTEP1. This is the first report on such a level of TEP diversity and of the presence of macroglobulin complement-related proteins (MCR) in mollusks. Gene structure analysis revealed alternative splicing in the highly variable region of three members (BgA2M, BgCD109, and BgTEP2) with a particularly unexpected diversity for BgTEP2. Finally, different gene expression profiles tend to indicate specific functions for such novel family members.

Evolutionary superscaffolding and chromosome anchoring to improve Anopheles genome assemblies

Background

New sequencing technologies have lowered financial barriers to whole genome sequencing, but resulting assemblies are often fragmented and far from ‘finished’. Updating multi-scaffold drafts to chromosome-level status can be achieved through experimental mapping or re-sequencing efforts. Avoiding the costs associated with such approaches, comparative genomic analysis of gene order conservation (synteny) to predict scaffold neighbours (adjacencies) offers a potentially useful complementary method for improving draft assemblies.

Results

We evaluated and employed 3 gene synteny-based methods applied to 21 Anopheles mosquito assemblies to produce consensus sets of scaffold adjacencies. For subsets of the assemblies, we integrated these with additional supporting data to confirm and complement the synteny-based adjacencies: 6 with physical mapping data that anchor scaffolds to chromosome locations, 13 with paired-end RNA sequencing (RNAseq) data, and 3 with new assemblies based on re-scaffolding or long-read data. Our combined analyses produced 20 new superscaffolded assemblies with improved contiguities: 7 for which assignments of non-anchored scaffolds to chromosome arms span more than 75% of the assemblies, and a further 7 with chromosome anchoring including an 88% anchored Anopheles arabiensis assembly and, respectively, 73% and 84% anchored assemblies with comprehensively updated cytogenetic photomaps for Anopheles funestus and Anopheles stephensi.

Conclusions

Experimental data from probe mapping, RNAseq, or long-read technologies, where available, all contribute to successful upgrading of draft assemblies. Our evaluations show that gene synteny-based computational methods represent a valuable alternative or complementary approach. Our improved Anopheles reference assemblies highlight the utility of applying comparative genomics approaches to improve community genomic resources.

Genomic landscape and genetic manipulation of the black soldier fly Hermetia illucens, a natural waste recycler

The black soldier fly (BSF), Hermetia illucens (Diptera: Stratiomyidae), is renowned for its bioconversion of organic waste into a sustainable source of animal feed. We report a high-quality genome of 1.1 Gb and a consensus set of 16,770 gene models for this beneficial species. Compared to those of other dipteran species, the BSF genome has undergone a substantial expansion in functional modules related to septic adaptation, including immune system factors, olfactory receptors, and cytochrome P450s. We further profiled midgut transcriptomes and associated microbiomes of BSF larvae fed with representative types of organic waste. We find that the pathways related to digestive system and fighting infection are commonly enriched and that Firmicutes bacteria dominate the microbial community in BSF across all diets. To extend its potential practical applications, we further developed an efficient CRISPR/Cas9-based gene editing approach and implemented this to yield flightless and enhanced feeding capacity phenotypes, both of which could expand BSF production capabilities. Our study provides valuable genomic and technical resources for optimizing BSF lines for industrialization.