Zhijian Tu

A male-determining factor in the mosquito Aedes aegypti

Manipulating M factor alters mosquito sex Female mosquitoes feed on blood and in so doing transmit pathogens to millions annually. Although the molecular mechanism for determining sex in many animals is known, the specific factors in mosquitoes have been elusive. This is because sex determination in insects involves a section of the genome that is highly repetitive. Hall et al. now identify a male-determining factor (M factor) in Aedes aegypti. Manipulation of the M factor produced sex-change phenotypes. Knocking out the gene Nix resulted in feminized males, and ectopic expression gave masculinized females. These findings should help to advance strategies for converting female mosquitoes into nonbiting males. Science, this issue p. 1268 An M-locus gene is necessary and sufficient for male development in the mosquito that transmits dengue and yellow fever. Sex determination in the mosquito Aedes aegypti is governed by a dominant male-determining factor (M factor) located within a Y chromosome–like region called the M locus. Here, we show that an M-locus gene, Nix, functions as an M factor in A. aegypti. Nix exhibits persistent M linkage and early embryonic expression, two characteristics required of an M factor. Nix knockout with clustered regularly interspaced short palindromic repeats (CRISPR)–Cas9 resulted in largely feminized genetic males and the production of female isoforms of two key regulators of sexual differentiation: doublesex and fruitless. Ectopic expression of Nix resulted in genetic females with nearly complete male genitalia. Thus, Nix is both required and sufficient to initiate male development. This study provides a foundation for mosquito control strategies that convert female mosquitoes into harmless males.

Genome analysis of a major urban malaria vector mosquito, Anopheles stephensi

BackgroundAnopheles stephensi is the key vector of malaria throughout the Indian subcontinent and Middle East and an emerging model for molecular and genetic studies of mosquito-parasite interactions. The type form of the species is responsible for the majority of urban malaria transmission across its range.ResultsHere, we report the genome sequence and annotation of the Indian strain of the type form of An. stephensi. The 221 Mb genome assembly represents more than 92% of the entire genome and was produced using a combination of 454, Illumina, and PacBio sequencing. Physical mapping assigned 62% of the genome onto chromosomes, enabling chromosome-based analysis. Comparisons between An. stephensi and An. gambiae reveal that the rate of gene order reshuffling on the X chromosome was three times higher than that on the autosomes. An. stephensi has more heterochromatin in pericentric regions but less repetitive DNA in chromosome arms than An. gambiae. We also identify a number of Y-chromosome contigs and BACs. Interspersed repeats constitute 7.1% of the assembled genome while LTR retrotransposons alone comprise more than 49% of the Y contigs. RNA-seq analyses provide new insights into mosquito innate immunity, development, and sexual dimorphism.ConclusionsThe genome analysis described in this manuscript provides a resource and platform for fundamental and translational research into a major urban malaria vector. Chromosome-based investigations provide unique perspectives on Anopheles chromosome evolution. RNA-seq analysis and studies of immunity genes offer new insights into mosquito biology and mosquito-parasite interactions.

Six novel Y chromosome genes in Anopheles mosquitoes discovered by independently sequencing males and females

BackgroundY chromosomes are responsible for the initiation of male development, male fertility, and other male-related functions in diverse species. However, Y genes are rarely characterized outside a few model species due to the arduous nature of studying the repeat-rich Y.ResultsThe chromosome quotient (CQ) is a novel approach to systematically discover Y chromosome genes. In the CQ method, genomic DNA from males and females is sequenced independently and aligned to candidate reference sequences. The female to male ratio of the number of alignments to a reference sequence, a parameter called the chromosome quotient (CQ), is used to determine whether the sequence is Y-linked. Using the CQ method, we successfully identified known Y sequences from Homo sapiens and Drosophila melanogaster. The CQ method facilitated the discovery of Y chromosome sequences from the malaria mosquitoes Anopheles stephensi and An. gambiae. Comparisons to transcriptome sequence data with blastn led to the discovery of six Anopheles Y genes, three from each species. All six genes are expressed in the early embryo. Two of the three An. stephensi Y genes were recently acquired from the autosomes or the X. Although An. stephensi and An. gambiae belong to the same subgenus, we found no evidence of Y genes shared between the species.ConclusionsThe CQ method can reliably identify Y chromosome sequences using the ratio of alignments from male and female sequence data. The CQ method is widely applicable to species with fragmented genome assemblies produced from next-generation sequencing data. Analysis of the six Y genes characterized in this study indicates rapid Y chromosome evolution between An. stephensi and An. gambiae. The Anopheles Y genes discovered by the CQ method provide unique markers for population and phylogenetic analysis, and opportunities for novel mosquito control measures through the manipulation of sexual dimorphism and fertility.

Control of Mosquito-Borne Infectious Diseases: Sex and Gene Drive

Sterile male releases have successfully reduced local populations of the dengue vector, Aedes aegypti, but challenges remain in scale and in separating sexes before release. The recent discovery of the first mosquito male determining factor (M factor) will facilitate our understanding of the genetic programs that initiate sexual development in mosquitoes. Manipulation of the M factor and possible intermediary factors may result in female-to-male conversion or female killing, enabling efficient sex separation and effective reduction of target mosquito populations. Given recent breakthroughs in the development of CRISPR-Cas9 reagents as a source of gene drive, more advanced technologies at driving maleness, the ultimate disease refractory phenotype, become possible and may represent efficient and self-limiting methods to control mosquito populations.

Identification of Early Zygotic Genes in the Yellow Fever Mosquito Aedes aegypti and Discovery of a Motif Involved in Early Zygotic Genome Activation

During early embryogenesis the zygotic genome is transcriptionally silent and all mRNAs present are of maternal origin. The maternal-zygotic transition marks the time over which embryogenesis changes its dependence from maternal RNAs to zygotically transcribed RNAs. Here we present the first systematic investigation of early zygotic genes (EZGs) in a mosquito species and focus on genes involved in the onset of transcription during 2–4 hr. We used transcriptome sequencing to identify the “pure” (without maternal expression) EZGs by analyzing transcripts from four embryonic time ranges of 0–2, 2–4, 4–8, and 8–12 hr, which includes the time of cellular blastoderm formation and up to the start of gastrulation. Blast of 16,789 annotated transcripts vs. the transcriptome reads revealed evidence for 63 (P<0.001) and 143 (P<0.05) nonmaternally derived transcripts having a significant increase in expression at 2–4 hr. One third of the 63 EZG transcripts do not have predicted introns compared to 10% of all Ae. aegypti genes. We have confirmed by RT-PCR that zygotic transcription starts as early as 2–3 hours. A degenerate motif VBRGGTA was found to be overrepresented in the upstream sequences of the identified EZGs using a motif identification software called SCOPE. We find evidence for homology between this motif and the TAGteam motif found in Drosophila that has been implicated in EZG activation. A 38 bp sequence in the proximal upstream sequence of a kinesin light chain EZG (KLC2.1) contains two copies of the mosquito motif. This sequence was shown to support EZG transcription by luciferase reporter assays performed on injected early embryos, and confers early zygotic activity to a heterologous promoter from a divergent mosquito species. The results of these studies are consistent with the model of early zygotic genome activation via transcriptional activators, similar to what has been found recently in Drosophila.

The expression profile of Aedes albopictus miRNAs is altered by dengue virus serotype-2 infection.

BackgroundAedes albopictus is an important vector of Dengue virus (DENV) and it has quickly invaded the tropical and temperate environments worldwide. A few studies have shown that, microRNAs (miRNAs) regulate mosquito defense against pathogens. However, there is no systematic analysis of the impact of DENV infection on miRNA expression in Ae. albopictus. We conducted this study to investigate the miRNA expression of Ae. albopictus upon DENV-2 infection using Illumina RNA sequencing.ResultsA total of 103 known and 5 novel candidate miRNAs were identified in DENV-2 infected and non-infected adult female Ae. albopictus. Comparative analysis indicated that 52 miRNAs were significantly down-regulated and 18 were up-regulated significantly after infection. Furthermore, RT-qPCR validated the expression patterns of eleven of these differentially expressed miRNAs. Targets prediction and functional analysis of these regulated miRNAs suggested that miR-34-5p and miR-87 might be involved in the anti-pathogen and immune responses.ConclusionThis is the first systematic study on the impact of DENV infection on miRNA expression in Ae. albopictus. Complex changes in miRNA expression suggest a potential role of miRNAs in antiviral responses by regulating immune-related genes. This investigation provides information concerning DENV-induced miRNAs and offers clues for identifying potential candidates for vector based antiviral strategies.

miRNA Genes of an Invasive Vector Mosquito, Aedes albopictus

Aedes albopictus, a vector of Dengue and Chikungunya viruses, is a robust invasive species in both tropical and temperate environments. MicroRNAs (miRNAs) regulate gene expression and biological processes including embryonic development, innate immunity and infection. While a number of miRNAs have been discovered in some mosquitoes, no comprehensive effort has been made to characterize them from different developmental stages from a single species. Systematic analysis of miRNAs in Ae. albopictus will improve our understanding of its basic biology and inform novel strategies to prevent virus transmission. Between 10–14 million Illumina sequencing reads per sample were obtained from embryos, larvae, pupae, adult males, sugar-fed and blood-fed adult females. A total of 119 miRNA genes represented by 215 miRNA or miRNA star (miRNA*) sequences were identified, 15 of which are novel. Eleven, two, and two of the newly-discovered miRNA genes appear specific to Aedes, Culicinae, and Culicidae, respectively. A number of miRNAs accumulate predominantly in one or two developmental stages and the large number that showed differences in abundance following a blood meal likely are important in blood-induced mosquito biology. Gene Ontology (GO) analysis of the targets of all Ae. albopictus miRNAs provides a useful starting point for the study of their functions in mosquitoes. This study is the first systematic analysis of miRNAs based on deep-sequencing of small RNA samples of all developmental stages of a mosquito species. A number of miRNAs are related to specific physiological states, most notably, pre- and post-blood feeding. The distribution of lineage-specific miRNAs is consistent with mosquito phylogeny and the presence of a number of Aedes-specific miRNAs likely reflects the divergence between the Aedes and Culex genera.

A New Chromosomal Phylogeny Supports the Repeated Origin of Vectorial Capacity in Malaria Mosquitoes of the Anopheles gambiae Complex

Understanding phylogenetic relationships within species complexes of disease vectors is crucial for identifying genomic changes associated with the evolution of epidemiologically important traits. However, the high degree of genetic similarity among sibling species confounds the ability to determine phylogenetic relationships using molecular markers. The goal of this study was to infer the ancestral–descendant relationships among malaria vectors and nonvectors of the Anopheles gambiae species complex by analyzing breakpoints of fixed chromosomal inversions in ingroup and several outgroup species. We identified genes at breakpoints of fixed overlapping chromosomal inversions 2Ro and 2Rp of An. merus using fluorescence in situ hybridization, a whole-genome mate-paired sequencing, and clone sequencing. We also mapped breakpoints of a chromosomal inversion 2La (common to An. merus, An. gambiae, and An. arabiensis) in outgroup species using a bioinformatics approach. We demonstrated that the “standard” 2R+p arrangement and “inverted” 2Ro and 2La arrangements are present in outgroup species Anopheles stephensi, Aedes aegypti, and Culex quinquefasciatus. The data indicate that the ancestral species of the An. gambiae complex had the 2Ro, 2R+p, and 2La chromosomal arrangements. The “inverted” 2Ro arrangement uniquely characterizes a malaria vector An. merus as the basal species in the complex. The rooted chromosomal phylogeny implies that An. merus acquired the 2Rp inversion and that its sister species An. gambiae acquired the 2R+o inversion from the ancestral species. The karyotype of nonvectors An. quadriannulatus A and B was derived from the karyotype of the major malaria vector An. gambiae. We conclude that the ability to effectively transmit human malaria had originated repeatedly in the complex. Our findings also suggest that saltwater tolerance originated first in An. merus and then independently in An. melas. The new chromosomal phylogeny will facilitate identifying the association of evolutionary genomic changes with epidemiologically important phenotypes.

Integrated proteomic and transcriptomic analysis of the Aedes aegypti eggshell

BackgroundMosquito eggshells show remarkable diversity in physical properties and structure consistent with adaptations to the wide variety of environments exploited by these insects. We applied proteomic, transcriptomic, and hybridization in situ techniques to identify gene products and pathways that participate in the assembly of the Aedes aegypti eggshell. Aedes aegypti population density is low during cold and dry seasons and increases immediately after rainfall. The survival of embryos through unfavorable periods is a key factor in the persistence of their populations. The work described here supports integrated vector control approaches that target eggshell formation and result in Ae. aegypti drought-intolerant phenotypes for public health initiatives directed to reduce mosquito-borne diseases.ResultsA total of 130 proteins were identified from the combined mass spectrometric analyses of eggshell preparations.ConclusionsClassification of proteins according to their known and putative functions revealed the complexity of the eggshell structure. Three novel Ae. aegypti vitelline membrane proteins were discovered. Odorant-binding and cysteine-rich proteins that may be structural components of the eggshell were identified. Enzymes with peroxidase, laccase and phenoloxidase activities also were identified, and their likely involvements in cross-linking reactions that stabilize the eggshell structure are discussed.

Identification and characterization of odorant-binding protein 1 gene from the Asian malaria mosquito, Anopheles stephensi

Insect odorant‐binding proteins (OBPs) are small, water‐soluble molecules that are thought to transport the hydrophobic odorants to their receptors in the chemosensory neurones. Here we report the identification and molecular characterization of the Anopheles stephensi odorant‐binding protein 1 gene (AsteObp1), an Obp1 gene in An. stephensi, a major malaria vector in Asia. We show that AsteObp1 and Anopheles gambiae Obp1 (AgamObp1) are orthologues. These two genes share similar coding sequences and conserved noncoding sequences (CNSs) that may be involved in their regulation. Transcript of AsteObp1 was observed in larvae and reached a relatively high level in late pupae. Quantitative real‐time PCR on female adult chemosensory tissues showed ∼900‐fold higher expression of AsteObp1 in antennae than in maxillary palp and proboscis. The amount of AsteObp1 in female legs was approximately 15‐fold lower than that of maxillary palp and proboscis. The level of AsteObp1 transcript was seven and 85‐fold higher in females than in males in the antennae, and maxillary palp and proboscis, respectively. Moreover, the AsteObp1 level was reduced by approximately 20‐fold in maxillary palp and proboscis 24 h after a bloodmeal. Our results indicate that AsteObp1 is likely to function in the female olfactory response and may also be involved in blood‐feeding behaviour.