Search

Author(s): S. T. Cole [1]; K. Eiglmeier [1]; J. Parkhill [2]; K. D. James [2]; N. R. Thomson [2]; P. R. Wheeler [3]; N. Honoré [1]; T. Garnier [1]; C. [...]

The nucleotide sequence of cosmid B1790, carrying the Rif-Str regions of the Mycobacterium leprae chromosome, has been determined. Twelve open reading frames were identified in the 36716bp sequence, representing 40% of the coding capacity. Five ribosomal proteins, two elongation factors and the β and β' subunits of RNA polymerase have been characterized and two novel genes were found. One of these encodes a member of the so-called ABC family of ATP-binding proteins while the other appears to encode an enzyme involved in repairing genomic lesions caused by free radicals. This finding may well be significant as M. leprae, an intracellular pathogen, lives within macrophages. [ABSTRACT FROM AUTHOR]

In an attempt to unify the genetic and biological research on Mycobacterium leprae, the aetiological agent of leprosy, a cosmid library was constructed and then ordered by a combination of fingerprinting and hybridization techniques. The genome of M. leprae is represented by four contigs of overlapping clones which, together, account for nearly 2.8Mb of DNA. Several arguments suggest that the gaps between the contigs are small in size and that virtually complete coverage of the chromosome has been obtained. All of the cloned M. leprae genes have been positioned on the contig maps together with the 29 copies of the dispersed repetitive element, RLEP. These have been classified into four groups on the basis of differences in their organization. Several key housekeeping genes were identified and mapped by hybridization with heterologous probes, and the current genome map of this uncultivable pathogen comprises 72 loci. [ABSTRACT FROM AUTHOR]

In enteric bacteria, the expression of many genes encoding various anaerobic electron transfer functions is controlled by FNR, the product of the autoregulated fnr gene. FNR is structurally and functionally homologous to CAP, the catabolite gene activator protein, and increased FNR production strongly stimulates transcription of its target genes. By analysis of RNA produced in vivo the promoters of four FNR-dependent genes were localized and shown to display a common arrangement. A 22bp dyad symmetry was found about 30 nucleotides upstream of the transcriptional startpoints and a similar sequence was shown to overlap the site of transcription initiation in the negatively controlled fnr gene. The consensus sequence for the half site recognized by FNR (AAATTGAT) is only slightly different from that of CAP (AA-TGTGA). Studies with two mutant frd promoters from Escherichia coli, displaying altered regulation and FNR response, provided additional evidence for recognition of this sequence by FNR. [ABSTRACT FROM AUTHOR]

Background: Anopheles cell lines are used in a variety of ways to better understand the major vectors of malaria in sub-Saharan Africa. Despite this, commonly used cell lines are not well characterized, and no tools are available for cell line identification and authentication., Methods: Utilizing whole genome sequencing, genomes of 4a-3A and 4a-3B 'hemocyte-like' cell lines were characterized for insertions and deletions (indels) and SNP variation. Genomic locations of distinguishing sequence variation and species origin of the cell lines were also examined. Unique indels were targeted to develop a PCR-based cell line authentication assay. Mitotic chromosomes were examined to survey the cytogenetic landscape for chromosome structure and copy number in the cell lines., Results: The 4a-3A and 4a-3B cell lines are female in origin and primarily of Anopheles coluzzii ancestry. Cytogenetic analysis indicates that the two cell lines are essentially diploid, with some relatively minor chromosome structural rearrangements. Whole-genome sequence was generated, and analysis indicated that SNPs and indels which differentiate the cell lines are clustered on the 2R chromosome in the regions of the 2Rb, 2Rc and 2Ru chromosomal inversions. A PCR-based authentication assay was developed to fingerprint three indels unique to each cell line. The assay distinguishes between 4a-3A and 4a-3B cells and also uniquely identifies two additional An. coluzzii cell lines tested, Ag55 and Sua4.0. The assay has the specificity to distinguish four cell lines and also has the sensitivity to detect cellular contamination within a sample of cultured cells., Conclusions: Genomic characterization of the 4a-3A and 4a-3B Anopheles cell lines was used to develop a simple diagnostic assay that can distinguish these cell lines within and across research laboratories. A cytogenetic survey indicated that the 4a-3A and Sua4.0 cell lines carry essentially normal diploid chromosomes, which makes them amenable to CRISPR/Cas9 genome editing. The presented simple authentication assay, coupled with screening for mycoplasma, will allow validation of the integrity of experimental resources and will promote greater experimental reproducibility of results., (© 2022. The Author(s).)

Entomopathogenic fungi have been explored as a potential biopesticide to counteract the insecticide resistance issue in mosquitoes. However, little is known about the possibility that genetic resistance to fungal biopesticides could evolve in mosquito populations. Here, we detected an important genetic component underlying Anopheles coluzzii survival after exposure to the entomopathogenic fungus Metarhizium anisopliae. A familiality study detected variation for survival among wild mosquito isofemale pedigrees, and genetic mapping identified two loci that significantly influence mosquito survival after fungus exposure. One locus overlaps with a previously reported locus for Anopheles susceptibility to the human malaria parasite Plasmodium falciparum. Candidate gene studies revealed that two LRR proteins encoded by APL1C and LRIM1 genes in this newly mapped locus are required for protection of female A. coluzzii from M. anisopliae, as is the complement-like factor Tep1. These results indicate that natural Anopheles populations already segregate frequent genetic variation for differential mosquito survival after fungal challenge and suggest a similarity in Anopheles protective responses against fungus and Plasmodium. However, this immune similarity raises the possibility that fungus-resistant mosquitoes could also display enhanced resistance to Plasmodium, suggesting an advantage of selecting for fungus resistance in vector populations to promote naturally diminished malaria vector competence., (© 2022. The Author(s).)

Exposure of mosquitoes to numerous eukaryotic and prokaryotic microbes in their associated microbiomes has probably helped drive the evolution of the innate immune system. To our knowledge, a metagenomic catalog of the eukaryotic microbiome has not been reported from any insect. Here we employ a novel approach to preferentially deplete host 18S ribosomal RNA gene amplicons to reveal the composition of the eukaryotic microbial communities of Anopheles larvae sampled in Kenya, Burkina Faso and Republic of Guinea (Conakry). We identified 453 eukaryotic operational taxonomic units (OTUs) associated with Anopheles larvae in nature, but an average of 45% of the 18S rRNA sequences clustered into OTUs that lacked a taxonomic assignment in the Silva database. Thus, the Anopheles microbiome contains a striking proportion of novel eukaryotic taxa. Using sequence similarity matching and de novo phylogenetic placement, the fraction of unassigned sequences was reduced to an average of 4%, and many unclassified OTUs were assigned as relatives of known taxa. A novel taxon of the genus Ophryocystis in the phylum Apicomplexa (which also includes Plasmodium ) is widespread in Anopheles larvae from East and West Africa. Notably, Ophryocystis is present at fluctuating abundance among larval breeding sites, consistent with the expected pattern of an epidemic pathogen. Species richness of the eukaryotic microbiome was not significantly different across sites from East to West Africa, while species richness of the prokaryotic microbiome was significantly lower in West Africa. Laboratory colonies of Anopheles coluzzii harbor 26 eukaryotic OTUs, of which 38% ( n = 10) are shared with wild populations, while 16 OTUs are unique to the laboratory colonies. Genetically distinct An. coluzzii colonies co-housed in the same facility maintain different prokaryotic microbiome profiles, suggesting a persistent host genetic influence on microbiome composition. These results provide a foundation to understand the role of the Anopheles eukaryotic microbiome in vector immunity and pathogen transmission. We hypothesize that prevalent apicomplexans such as Ophryocystis associated with Anopheles could induce interference or competition against Plasmodium within the vector. This and other members of the eukaryotic microbiome may offer candidates for new vector control tools., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Cuesta, Coulibaly, Bukhari, Eiglmeier, Kone, Coulibaly, Zongo, Barry, Gneme, Guelbeogo, Beavogui, Traore, Sagnon, Vernick and Riehle.)

In the Greater Mekong Subregion, malaria cases have significantly decreased but little is known about the vectors or mechanisms responsible for residual malaria transmission. We analysed a total of 3920 Anopheles mosquitoes collected during the rainy and dry seasons from four ecological settings in Cambodia (villages, forested areas near villages, rubber tree plantations and forest sites). Using odor-baited traps, 81% of the total samples across all sites were collected in cow baited traps, although 67% of the samples attracted by human baited traps were collected in forest sites. Overall, 20% of collected Anopheles were active during the day, with increased day biting during the dry season. 3131 samples were identified morphologically as 14 different species, and a subset was also identified by DNA amplicon sequencing allowing determination of 29 Anopheles species. The investigation of well characterized insecticide mutations (ace-1, kdr, and rdl genes) indicated that individuals carried mutations associated with response to all the different classes of insecticides. There also appeared to be a non-random association between mosquito species and insecticide resistance with Anopheles peditaeniatus exhibiting nearly fixed mutations. Molecular screening for Plasmodium sp. presence indicated that 3.6% of collected Anopheles were positive, most for P. vivax followed by P. falciparum. These results highlight some of the key mechanisms driving residual human malaria transmission in Cambodia, and illustrate the importance of diverse collection methods, sites and seasons to avoid bias and better characterize Anopheles mosquito ecology in Southeast Asia.

The commensal gut microbiome is contained by the enteric epithelial barrier, but little is known about the degree of specificity of host immune barrier interactions for particular bacterial taxa. Here, we show that depletion of leucine-rich repeat immune factor APL1 in the Asian malaria mosquito Anopheles stephensi is associated with higher midgut abundance of just the family Enterobacteraceae , and not generalized dysbiosis of the microbiome. The effect is explained by the response of a narrow clade containing two main taxa related to Klebsiella and Cedecea . Analysis of field samples indicate that these two taxa are recurrent members of the wild Anopheles microbiome. Triangulation using sequence and functional data incriminated relatives of C. neteri and Cedecea NFIX57 as candidates for the Cedecea component, and K. michiganensis , K. oxytoca , and K.sp. LTGPAF-6F as candidates for the Klebsiella component. APL1 presence is associated with host ability to specifically constrain the abundance of a narrow microbiome clade of the Enterobacteraceae , and the immune factor may promote homeostasis of this clade in the enteric microbiome for host benefit., (Copyright © 2020 Mitri, Bischoff, Belda Cuesta, Volant, Ghozlane, Eiglmeier, Holm, Dieme, Brito-Fravallo, Guelbeogo, Sagnon, Riehle and Vernick.)

Background: The recent reference genome assembly and annotation of the Asian malaria vector Anopheles stephensi detected only one gene encoding the leucine-rich repeat immune factor APL1, while in the Anopheles gambiae and sibling Anopheles coluzzii, APL1 factors are encoded by a family of three paralogs. The phylogeny and biological function of the unique APL1 gene in An. stephensi have not yet been specifically examined., Methods: The APL1 locus was manually annotated to confirm the computationally predicted single APL1 gene in An. stephensi. APL1 evolution within Anopheles was explored by phylogenomic analysis. The single or paralogous APL1 genes were silenced in An. stephensi and An. coluzzii, respectively, followed by mosquito survival analysis, experimental infection with Plasmodium and expression analysis., Results: APL1 is present as a single ancestral gene in most Anopheles including An. stephensi but has expanded to three paralogs in an African lineage that includes only the Anopheles gambiae species complex and Anopheles christyi. Silencing of the unique APL1 copy in An. stephensi results in significant mosquito mortality. Elevated mortality of APL1-depleted An. stephensi is rescued by antibiotic treatment, suggesting that pathology due to bacteria is the cause of mortality, and indicating that the unique APL1 gene is essential for host survival. Successful Plasmodium development in An. stephensi depends upon APL1 activity for protection from high host mortality due to bacteria. In contrast, silencing of all three APL1 paralogs in An. coluzzii does not result in elevated mortality, either with or without Plasmodium infection. Expression of the single An. stephensi APL1 gene is regulated by both the Imd and Toll immune pathways, while the two signaling pathways regulate different APL1 paralogs in the expanded APL1 locus., Conclusions: APL1 underwent loss and gain of functions concomitant with expansion from a single ancestral gene to three paralogs in one lineage of African Anopheles. We infer that activity of the unique APL1 gene promotes longevity in An. stephensi by conferring protection from or tolerance to an effect of bacterial pathology. The evolution of an expanded APL1 gene family could be a factor contributing to the exceptional levels of malaria transmission mediated by human-feeding members of the An. gambiae species complex in Africa.

Following the publication of this article [1], the authors reported that the original shading in columns 3 and 4 of Table 3, which indicated the presence or absence of viruses in each library, had been removed during typesetting.

Background: Mosquitoes are colonized by a large but mostly uncharacterized natural virome of RNA viruses, and the composition and distribution of the natural RNA virome may influence the biology and immunity of Anopheles malaria vector populations., Results: Anopheles mosquitoes were sampled in malaria endemic forest village sites in Senegal and Cambodia, including Anopheles funestus, Anopheles gambiae group sp., and Anopheles coustani in Senegal, and Anopheles hyrcanus group sp., Anopheles maculatus group sp., and Anopheles dirus in Cambodia. The most frequent mosquito species sampled at both study sites are human malaria vectors. Small and long RNA sequences were depleted of mosquito host sequences, de novo assembled and clustered to yield non-redundant contigs longer than 500 nucleotides. Analysis of the assemblies by sequence similarity to known virus families yielded 115 novel virus sequences, and evidence supports a functional status for at least 86 of the novel viral contigs. Important monophyletic virus clades in the Bunyavirales and Mononegavirales orders were found in these Anopheles from Africa and Asia. The remaining non-host RNA assemblies that were unclassified by sequence similarity to known viruses were clustered by small RNA profiles, and 39 high-quality independent contigs strongly matched a pattern of classic RNAi processing of viral replication intermediates, suggesting they are entirely undescribed viruses. One thousand five hundred sixty-six additional high-quality unclassified contigs matched a pattern consistent with Piwi-interacting RNAs (piRNAs), suggesting that strand-biased piRNAs are generated from the natural virome in Anopheles. To functionally query piRNA effect, we analyzed piRNA expression in Anopheles coluzzii after infection with O'nyong nyong virus (family Togaviridae), and identified two piRNAs that appear to display specifically altered abundance upon arbovirus infection., Conclusions: Anopheles vectors of human malaria in Africa and Asia are ubiquitously colonized by RNA viruses, some of which are monophyletic but clearly diverged from other arthropod viruses. The interplay between small RNA pathways, immunity, and the virome may represent part of the homeostatic mechanism maintaining virome members in a commensal or nonpathogenic state, and could potentially influence vector competence.

Background: In many African countries malaria has declined sharply due to a synergy of actions marked by the introduction of vector control strategies, but the disease remains a leading cause of morbidity and mortality in Central African Republic (CAR). An entomological study was initiated with the aim to characterize the malaria vectors in Bangui, the capital of CAR, and determine their vector competence., Methods: A cross-sectional entomological study was conducted in 15 sites of the district of Bangui, the capital of CAR, in September-October 2013 and a second collection was done in four of those sites between November and December 2013. Mosquitoes were collected by human landing catch (HLC) indoors and outdoors and by pyrethrum spray catch of indoor-resting mosquitoes. Mosquitoes were analysed for species and multiple other attributes, including the presence of Plasmodium falciparum circumsporozoite protein or DNA, blood meal source, 2La inversion karyotype, and the L1014F kdr insecticide resistance mutation., Results: Overall, 1292 anophelines were analysed, revealing a predominance of Anopheles gambiae and Anopheles funestus, with a small fraction of Anopheles coluzzii. Molecular typing of the An. gambiae complex species showed that An. gambiae was predominant (95.7 %) as compared to An. coluzzii (2.1 %), and Anopheles arabiensis was not present. In some areas the involvement of secondary vectors, such as Anopheles coustani, expands the risk of infection. By HLC sampling, An. funestus displayed a stronger endophilic preference than mosquitoes from the An. gambiae sister taxa, with a mean indoor-capture rate of 54.3 % and 67.58 % for An. gambiae sister taxa and An. funestus, respectively. Human biting rates were measured overall for each of the species with 28 or 29 bites/person/night, respectively. Both vectors displayed a strong human feeding preference as determined by blood meal source, which was not different between the different sampling sites. An. coustani appears to be highly exophilic, with 92 % of HLC samples captured outdoors. The mean CSP rate in head-thorax sections of all Anopheles was 5.09 %, and was higher in An. gambiae s.l. (7.4 %) than in An. funestus (3.3 %). CSP-positive An. coustani were also detected in outdoor HLC samples. In the mosquitoes of the An. gambiae sister taxa the kdr-w mutant allele was nearly fixed, with 92.3 % resistant homozygotes, and no susceptible homozygotes detected., Conclusions: This study collected data on anopheline populations in CAR, behaviour of vectors and transmission levels. Further studies should investigate the biting behaviour and susceptibility status of the anophelines to different insecticides to allow the establishment of appropriate vector control based on practical entomological knowledge.

Mosquitoes of the Anopheles gambiae complex display strong preference for human bloodmeals and are major malaria vectors in Africa. However, their interaction with viruses or role in arbovirus transmission during epidemics has been little examined, with the exception of O'nyong-nyong virus, closely related to Chikungunya virus. Deep-sequencing has revealed different RNA viruses in natural insect viromes, but none have been previously described in the Anopheles gambiae species complex. Here, we describe two novel insect RNA viruses, a Dicistrovirus and a Cypovirus, found in laboratory colonies of An. gambiae taxa using small-RNA deep sequencing. Sequence analysis was done with Metavisitor, an open-source bioinformatic pipeline for virus discovery and de novo genome assembly. Wild-collected Anopheles from Senegal and Cambodia were positive for the Dicistrovirus and Cypovirus, displaying high sequence identity to the laboratory-derived virus. Thus, the Dicistrovirus (Anopheles C virus, AnCV) and Cypovirus (Anopheles Cypovirus, AnCPV) are components of the natural virome of at least some anopheline species. Their possible influence on mosquito immunity or transmission of other pathogens is unknown. These natural viruses could be developed as models for the study of Anopheles-RNA virus interactions in low security laboratory settings, in an analogous manner to the use of rodent malaria parasites for studies of mosquito anti-parasite immunity.

Members of the Anopheles gambiae species complex are primary vectors of human malaria in Africa. Population heterogeneities for ecological and behavioral attributes expand and stabilize malaria transmission over space and time, and populations may change in response to vector control, urbanization and other factors. There is a need for approaches to comprehensively describe the structure and characteristics of a sympatric local mosquito population, because incomplete knowledge of vector population composition may hinder control efforts. To this end, we used a genome-wide custom SNP typing array to analyze a population collection from a single geographic region in West Africa. The combination of sample depth (n = 456) and marker density (n = 1536) unambiguously resolved population subgroups, which were also compared for their relative susceptibility to natural genotypes of Plasmodium falciparum malaria. The population subgroups display fluctuating patterns of differentiation or sharing across the genome. Analysis of linkage disequilibrium identified 19 new candidate genes for association with underlying population divergence between sister taxa, A. coluzzii (M-form) and A. gambiae (S-form).

Nucleotide variation patterns across species are shaped by the processes of natural selection, including exposure to environmental pathogens. We examined patterns of genetic variation in two sister species, Anopheles gambiae and Anopheles coluzzii, both efficient natural vectors of human malaria in West Africa. We used the differentiation signature displayed by a known coordinate selective sweep of immune genes APL1 and TEP1 in A. coluzzii to design a population genetic screen trained on the sweep, classified a panel of 26 potential immune genes for concordance with the signature, and functionally tested their immune phenotypes. The screen results were strongly predictive for genes with protective immune phenotypes: genes meeting the screen criteria were significantly more likely to display a functional phenotype against malaria infection than genes not meeting the criteria (p = 0.0005). Thus, an evolution-based screen can efficiently prioritize candidate genes for labor-intensive downstream functional testing, and safely allow the elimination of genes not meeting the screen criteria. The suite of immune genes with characteristics similar to the APL1-TEP1 selective sweep appears to be more widespread in the A. coluzzii genome than previously recognized. The immune gene differentiation may be a consequence of adaptation of A. coluzzii to new pathogens encountered in its niche expansion during the separation from A. gambiae, although the role, if any of natural selection by Plasmodium is unknown. Application of the screen allowed identification of new functional immune factors, and assignment of new functions to known factors. We describe biochemical binding interactions between immune proteins that underlie functional activity for malaria infection, which highlights the interplay between pathogen specificity and the structure of immune complexes. We also find that most malaria-protective immune factors display phenotypes for either human or rodent malaria, with broad specificity a rarity.

Background: The genome-wide association study (GWAS) techniques that have been used for genetic mapping in other organisms have not been successfully applied to mosquitoes, which have genetic characteristics of high nucleotide diversity, low linkage disequilibrium, and complex population stratification that render population-based GWAS essentially unfeasible at realistic sample size and marker density., Methods: We designed a novel mapping strategy for the mosquito system that combines the power of linkage mapping with the resolution afforded by genetic association. We established founder colonies from West Africa, controlled for diversity, linkage disequilibrium and population stratification. Colonies were challenged by feeding on the infectious stage of the human malaria parasite, Plasmodium falciparum, mosquitoes were phenotyped for parasite load, and DNA pools for phenotypically similar mosquitoes were Illumina sequenced. Phenotype-genotype mapping was carried out in two stages, coarse and fine., Results: In the first mapping stage, pooled sequences were analysed genome-wide for intervals displaying relativereduction in diversity between phenotype pools, and candidate genomic loci were identified for influence upon parasite infection levels. In the second mapping stage, focused genotyping of SNPs from the first mapping stage was carried out in unpooled individual mosquitoes and replicates. The second stage confirmed significant SNPs in a locus encoding two Toll-family proteins. RNAi-mediated gene silencing and infection challenge revealed that TOLL 11 protects mosquitoes against P. falciparum infection., Conclusions: We present an efficient and cost-effective method for genetic mapping using natural variation segregating in defined recent Anopheles founder colonies, and demonstrate its applicability for mapping in a complex non-model genome. This approach is a practical and preferred alternative to population-based GWAS for first-pass mapping of phenotypes in Anopheles. This design should facilitate mapping of other traits involved in physiology, epidemiology, and behaviour.

Background: Members of the Anopheles gambiae species complex are primary vectors of human malaria in Africa. It is known that a large haplotype shared between An. gambiae and Anopheles coluzzii by introgression carries point mutations of the voltage-gated sodium channel gene para, including the L1014F kdr mutation associated with insensitivity to pyrethroid insecticides. Carriage of L1014F kdr is also correlated with higher susceptibility to infection with Plasmodium falciparum. However, the genetic mechanism and causative gene(s) underlying the parasite susceptibility phenotype are not known., Methods: Mosquitoes from the wild Burkina Faso population were challenged by feeding on natural P. falciparum gametocytes. Oocyst infection phenotypes were determined and were tested for association with SNP genotypes. Candidate genes in the detected locus were prioritized and RNAi-mediated gene silencing was used to functionally test for gene effects on P. falciparum susceptibility., Results: A genetic locus, Pfin6, was identified that influences infection levels of P. falciparum in mosquitoes. The locus segregates as a ~3 Mb haplotype carrying 65 predicted genes including the para gene. The haplotype carrying the kdr allele of para is linked to increased parasite infection prevalence, but many single nucleotide polymorphisms on the haplotype are also equally linked to the infection phenotype. Candidate genes in the haplotype were prioritized and functionally tested. Silencing of para did not influence P. falciparum infection, while silencing of a predicted immune gene, serine protease ClipC9, allowed development of significantly increased parasite numbers., Conclusions: Genetic variation influencing Plasmodium infection in wild Anopheles is linked to a natural ~3 megabase haplotype on chromosome 2L that carries the kdr allele of the para gene. Evidence suggests that para gene function does not directly influence parasite susceptibility, and the association of kdr with infection may be due to tight linkage of kdr with other gene(s) on the haplotype. Further work will be required to determine if ClipC9 influences the outcome of P. falciparum infection in nature, as well as to confirm the absence of a direct influence by para.

Background: Anopheles 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., Results: Here, 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., Conclusions: The 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.

Background: The Anopheles gambiae sensu lato (s.l.) species complex in Burkina Faso consists of Anopheles arabiensis, and molecular forms M and S of Anopheles gambiae sensu stricto (s.s.). Previous studies comparing the M and S forms for level of infection with Plasmodium falciparum have yielded conflicting results., Methods: Mosquito larvae were sampled from natural pools, reared to adulthood under controlled conditions, and challenged with natural P. falciparum by experimental feeding with blood from gametocyte carriers. Oocyst infection prevalence and intensity was determined one week after infection. DNA from carcasses was genotyped to identify species and molecular form., Results: In total, 7,400 adult mosquitoes grown from wild-caught larvae were challenged with gametocytes in 29 experimental infections spanning four transmission seasons. The overall infection prevalence averaged 40.7% for A. gambiae M form, 41.4% for A. gambiae S form, and 40.1% for A. arabiensis. There was no significant difference in infection prevalence or intensity between the three population groups. Notably, infection experiments in which the population groups were challenged in parallel on the same infective blood displayed less infection difference between population groups, while infections with less balanced composition of population groups had lower statistical power and displayed apparent differences that fluctuated more often from the null average., Conclusion: The study clearly establishes that, at the study site in Burkina Faso, there is no difference in genetic susceptibility to P. falciparum infection between three sympatric population groups of the A. gambiae s.l. complex. Feeding the mosquito groups on the same infective blood meal greatly increases statistical power. Conversely, comparison of the different mosquito groups between, rather than within, infections yields larger apparent difference between mosquito groups, resulting from lower statistical power and greater noise, and could lead to false-positive results. In making infection comparisons between population groups, it is more accurate to compare the different groups after feeding simultaneously upon the same infective blood.

Host-pathogen interactions can be powerful drivers of adaptive evolution, shaping the patterns of molecular variation at the genes involved. In this study, we sequenced alleles from 28 immune-related loci in wild samples of multiple genetic subpopulations of the African malaria mosquito Anopheles gambiae, obtaining unprecedented sample sizes and providing the first opportunity to contrast patterns of molecular evolution at immune-related loci in the recently discovered GOUNDRY population to those of the indoor-resting M and S molecular forms. In contrast to previous studies that focused on immune genes identified in laboratory studies, we centered our analysis on genes that fall within a quantitative trait locus associated with resistance to Plasmodium falciparum in natural populations of A. gambiae. Analyses of haplotypic and genetic diversity at these 28 loci revealed striking differences among populations in levels of genetic diversity and allele frequencies in coding sequence. Putative signals of positive selection were identified at 11 loci, but only one was shared by two subgroups of A. gambiae. Striking patterns of linkage disequilibrium were observed at several loci. We discuss these results with respect to ecological differences among these strata as well as potential implications for disease transmission.

Functional studies have demonstrated a role for the Anopheles gambiae APL1A gene in resistance against the human malaria parasite, Plasmodium falciparum. Here, we exhaustively characterize the structure of the APL1 locus and show that three structurally different APL1A alleles segregate in the Ngousso colony. Genetic association combined with RNAi-mediated gene silencing revealed that APL1A alleles display distinct protective profiles against P. falciparum. One APL1A allele is sufficient to explain the protective phenotype of APL1A observed in silencing experiments. Epitope-tagged APL1A isoforms expressed in an in vitro hemocyte-like cell system showed that under assay conditions, the most protective APL1A isoform (APL1A(2)) localizes within large cytoplasmic vesicles, is not constitutively secreted, and forms only one protein complex, while a less protective isoform (APL1A(1)) is constitutively secreted in at least two protein complexes. The tested alleles are identical to natural variants in the wild A. gambiae population, suggesting that APL1A genetic variation could be a factor underlying natural heterogeneity of vector susceptibility to P. falciparum.

Population subgroups of the African malaria vector Anopheles gambiae have not been comprehensively characterized owing to the lack of unbiased sampling methods. In the arid savanna zone of West Africa, where potential oviposition sites are scarce, widespread collection from larval pools in the peridomestic human habitat yielded a comprehensive genetic survey of local A. gambiae population subgroups, independent of adult resting behavior and ecological preference. A previously unknown subgroup of exophilic A. gambiae is sympatric with the known endophilic A. gambiae in this region. The exophilic subgroup is abundant, lacks differentiation into M and S molecular forms, and is highly susceptible to infection with wild Plasmodium falciparum. These findings might have implications for the epidemiology of malaria transmission and control.

We present a draft sequence of the genome of Aedes aegypti, the primary vector for yellow fever and dengue fever, which at approximately 1376 million base pairs is about 5 times the size of the genome of the malaria vector Anopheles gambiae. Nearly 50% of the Ae. aegypti genome consists of transposable elements. These contribute to a factor of approximately 4 to 6 increase in average gene length and in sizes of intergenic regions relative to An. gambiae and Drosophila melanogaster. Nonetheless, chromosomal synteny is generally maintained among all three insects, although conservation of orthologous gene order is higher (by a factor of approximately 2) between the mosquito species than between either of them and the fruit fly. An increase in genes encoding odorant binding, cytochrome P450, and cuticle domains relative to An. gambiae suggests that members of these protein families underpin some of the biological differences between the two mosquito species.

To understand the evolution, attenuation, and variable protective efficacy of bacillus Calmette-Guérin (BCG) vaccines, Mycobacterium bovis BCG Pasteur 1173P2 has been subjected to comparative genome and transcriptome analysis. The 4,374,522-bp genome contains 3,954 protein-coding genes, 58 of which are present in two copies as a result of two independent tandem duplications, DU1 and DU2. DU1 is restricted to BCG Pasteur, although four forms of DU2 exist; DU2-I is confined to early BCG vaccines, like BCG Japan, whereas DU2-III and DU2-IV occur in the late vaccines. The glycerol-3-phosphate dehydrogenase gene, glpD2, is one of only three genes common to all four DU2 variants, implying that BCG requires higher levels of this enzyme to grow on glycerol. Further amplification of the DU2 region is ongoing, even within vaccine preparations used to immunize humans. An evolutionary scheme for BCG vaccines was established by analyzing DU2 and other markers. Lesions in genes encoding sigma-factors and pleiotropic transcriptional regulators, like PhoR and Crp, were also uncovered in various BCG strains; together with gene amplification, these affect gene expression levels, immunogenicity, and, possibly, protection against tuberculosis. Furthermore, the combined findings suggest that early BCG vaccines may even be superior to the later ones that are more widely used.

Background: Complete genome annotation is a necessary tool as Anopheles gambiae researchers probe the biology of this potent malaria vector., Results: We reannotate the A. gambiae genome by synthesizing comparative and ab initio sets of predicted coding sequences (CDSs) into a single set using an exon-gene-union algorithm followed by an open-reading-frame-selection algorithm. The reannotation predicts 20,970 CDSs supported by at least two lines of evidence, and it lowers the proportion of CDSs lacking start and/or stop codons to only approximately 4%. The reannotated CDS set includes a set of 4,681 novel CDSs not represented in the Ensembl annotation but with EST support, and another set of 4,031 Ensembl-supported genes that undergo major structural and, therefore, probably functional changes in the reannotated set. The quality and accuracy of the reannotation was assessed by comparison with end sequences from 20,249 full-length cDNA clones, and evaluation of mass spectrometry peptide hit rates from an A. gambiae shotgun proteomic dataset confirms that the reannotated CDSs offer a high quality protein database for proteomics. We provide a functional proteomics annotation, ReAnoXcel, obtained by analysis of the new CDSs through the AnoXcel pipeline, which allows functional comparisons of the CDS sets within the same bioinformatic platform. CDS data are available for download., Conclusion: Comprehensive A. gambiae genome reannotation is achieved through a combination of comparative and ab initio gene prediction algorithms.

The only natural mechanism of malaria transmission in sub-Saharan Africa is the mosquito, generally Anopheles gambiae. Blocking malaria parasite transmission by stopping the development of Plasmodium in the insect vector would provide a useful alternative to the current methods of malaria control. Toward this end, it is important to understand the molecular basis of the malaria parasite refractory phenotype in An. gambiae mosquito strains. We have selected and sequenced six bacterial artificial chromosome (BAC) clones from the Pen-1 region that is the major quantitative trait locus involved in Plasmodium encapsulation. The sequence and the annotation of five overlapping BAC clones plus one adjacent, but not contiguous clone, totaling 585kb of genomic sequence from the centromeric end of the Pen-1 region of the PEST strain were compared to that of the genome sequence of the same strain produced by the whole genome shotgun technique. This project identified 23 putative mosquito genes plus putative copies of the retrotransposable elements BEL12 and TRANSIBN1_AG in the six BAC clones. Nineteen of the predicted genes are most similar to their Drosophila melanogaster homologs while one is more closely related to vertebrate genes. Comparison of these new BAC sequences plus previously published BAC sequences to the cognate region of the assembled genome sequence identified three retrotransposons present in one sequence version but not the other. One of these elements, Indy, has not been previously described. These observations provide evidence for the recent active transposition of these elements and demonstrate the plasticity of the Anopheles genome. The BAC sequences strongly support the public whole genome shotgun assembly and automatic annotation while also demonstrating the benefit of complementary genome sequences and of human curation. Importantly, the data demonstrate the differences in the genome sequence of an individual mosquito compared to that of a hypothetical, average genome sequence generated by whole genome shotgun assembly.

We describe the preliminary analysis of over 35,000 clones from a full-length enriched cDNA library from the malaria mosquito vector Anopheles gambiae. The clones define nearly 3,700 genes, of which around 2,600 significantly improve current gene definitions. An additional 17% of the genes were not previously annotated, suggesting that an equal percentage may be missing from the current Anopheles genome annotation.

We performed genome-wide sequence comparisons at the protein coding level between the genome sequences of Drosophila melanogaster and Anopheles gambiae. Such comparisons detect evolutionarily conserved regions (ecores) that can be used for a qualitative and quantitative evaluation of the available annotations of both genomes. They also provide novel candidate features for annotation. The percentage of ecores mapping outside annotations in the A. gambiae genome is about fourfold higher than in D. melanogaster. The A. gambiae genome assembly also contains a high proportion of duplicated ecores, possibly resulting from artefactual sequence duplications in the genome assembly. The occurrence of 4063 ecores in the D. melanogaster genome outside annotations suggests that some genes are not yet or only partially annotated. The present work illustrates the power of comparative genomics approaches towards an exhaustive and accurate establishment of gene models and gene catalogues in insect genomes.

Mycobacterium bovis is the causative agent of tuberculosis in a range of animal species and man, with worldwide annual losses to agriculture of $3 billion. The human burden of tuberculosis caused by the bovine tubercle bacillus is still largely unknown. M. bovis was also the progenitor for the M. bovis bacillus Calmette-Guérin vaccine strain, the most widely used human vaccine. Here we describe the 4,345,492-bp genome sequence of M. bovis AF2122/97 and its comparison with the genomes of Mycobacterium tuberculosis and Mycobacterium leprae. Strikingly, the genome sequence of M. bovis is >99.95% identical to that of M. tuberculosis, but deletion of genetic information has led to a reduced genome size. Comparison with M. leprae reveals a number of common gene losses, suggesting the removal of functional redundancy. Cell wall components and secreted proteins show the greatest variation, indicating their potential role in host-bacillus interactions or immune evasion. Furthermore, there are no genes unique to M. bovis, implying that differential gene expression may be the key to the host tropisms of human and bovine bacilli. The genome sequence therefore offers major insight on the evolution, host preference, and pathobiology of M. bovis.

Mycobacterium microti is a member of the Mycobacterium tuberculosis complex that causes tuberculosis in voles. Most strains of M. microti are harmless for humans, and some have been successfully used as live tuberculosis vaccines. In an attempt to identify putative virulence factors of the tubercle bacilli, genes that are absent from the avirulent M. microti but present in human pathogen M. tuberculosis or Mycobacterium bovis were searched for. A minimal set of 50 bacterial artificial chromosome (BAC) clones that covers almost all of the genome of M. microti OV254 was constructed, and individual BACs were compared to the corresponding BACs from M. bovis AF2122/97 and M. tuberculosis H37Rv. Comparison of pulsed-field gel-separated DNA digests of BAC clones led to the identification of 10 regions of difference (RD) between M. microti OV254 and M. tuberculosis. A 14-kb chromosomal region (RD1(mic)) that partly overlaps the RD1 deletion in the BCG vaccine strain was missing from the genomes of all nine tested M. microti strains. This region covers 13 genes, Rv3864 to Rv3876, in M. tuberculosis, including those encoding the potent ESAT-6 and CFP-10 antigens. In contrast, RD5(mic), a region that contains three phospholipase C genes (plcA to -C), was missing from only the vole isolates and was present in M. microti strains isolated from humans. Apart from RD1(mic) and RD5(mic) other M. microti-specific deleted regions have been identified (MiD1 to MiD3). Deletion of MiD1 has removed parts of the direct repeat region in M. microti and thus contributes to the characteristic spoligotype of M. microti strains.

The small size of their genomes made bacterial ideal model organisms for the emerging field of genomics. Elucidating the genome sequences of mycobacteria was particularly attractive owing to the difficulties inherent in their manipulation. The slow growth rate, clumping, and requirement for category III containment make manipulation of Mycobacterium tuberculosis-complex strains laborious. M. leprae presents even greater problems as it has resisted all attempts at axenic culture. Availability of genome sequence data promised to accelerate our knowledge of the fundamental biology of these organisms, and to offer clues to the basis for their virulence, tropism and persistence in the host. This article will focus on what the genome sequences of M. tuberculosis and M. leprae have taught us about these pathogens, and how comparative genomics has exposed some of the fundamental differences between the species.

Everything that we need to know about Mycobacterium leprae, a close relative of the tubercle bacillus, is encrypted in its genome. Inspection of the 3.27 Mb genome sequence of an armadillo-derived Indian isolate of the leprosy bacillus identified 1,605 genes encoding proteins and 50 genes for stable RNA species. Comparison with the genome sequence of Mycobacterium tuberculosis revealed an extreme case of reductive evolution, since less than half of the genome contains functional genes while inactivated or pseudogenes are highly abundant. The level of gene duplication was approximately 34% and, on classification of the proteins into families, the largest functional groups were found to be involved in the metabolism and modification of fatty acids and polyketides, transport of metabolites, cell envelope synthesis and gene regulation. Reductive evolution, gene decay and genome downsizing have eliminated entire metabolic pathways, together with their regulatory circuits and accessory functions, particularly those involved in catabolism. This may explain the unusually long generation time and account for our inability to culture the leprosy bacillus.

The integrated map of the Mycobacterium leprae genome unveiled for the first time the genomic organization of this obligate intracellular parasite. Selected cosmid clones, isolated from a genomic library created in the cosmid vector Lorist6, were identified as representing nearly the complete genome and were subsequently used in the M. leprae genome sequencing project. Now a new version of the integrated map of M. leprae can be presented, combining the mapping results from the Lorist6 cosmids with data obtained from a second genomic library constructed in an Escherichia coli-mycobacterium shuttle cosmid, pYUB18. More than 98% of the M. leprae genome is now covered by overlapping large insert genomic clones representing a renewable source of well defined DNA segments and a powerful tool for functional genomics.

The imminent completion of the genome sequence of Mycobacterium bovis will reveal the genetic blueprint for this most successful pathogen. Comparative analysis with the genome sequences of M. tuberculosis and M. bovis BCG promises to expose the genetic basis for the phenotypic differences between the tubercle bacilli, offering unparalleled insight into the virulence factors of the M. tuberculosis complex. Initial analysis of the sequence data has already revealed a novel deletion from M. bovis, as well as identifying variation in members of the PPE family of proteins. As the study of bacterial pathogenicity enters the postgenomic phase, the genome sequence of M. bovis promises to serve as a cornerstone of mycobacterial genetics., (Copyright 2001 Harcourt Publishers Ltd.)

The genus mycobacteria includes two important human pathogens Mycobacterium tuberculosis and Mycobacterium lepra. The former is reputed to have the highest annual global mortality of all pathogens. Their slow growth, virulence for humans and particular physiology makes these organisms extremely difficult to work with. However the rapid development of mycobacterial genomics following the completion of the Mycobacterium tuberculosis genome sequence provides the basis for a powerful new approach for the understanding of these organisms. Five further genome sequencing projects of closely related mycobacterial species with differing host range, virulence for humans and physiology are underway. A comparative genomic analysis of these species has the potential to define the genetic basis of these phenotypes which will be invaluable for the development of urgently needed new vaccines and drugs. This minireview summarises the different techniques that have been employed to compare these genomes and gives an overview of the wealth of data that has already been generated by mycobacterial comparative genomics.

To achieve the quantum leap in understanding required to overcome two major human diseases, leprosy and tuberculosis, systematic and comparative genome analysis has been undertaken. New insight into the biology of their causative agents has been obtained and the principle findings are reported here.

Whole-genome comparisons of the tubercle bacilli were undertaken using ordered bacterial artificial chromosome (BAC) libraries of Mycobacterium tuberculosis and the vaccine strain, Mycobacterium bovis BCG-Pasteur, together with the complete genome sequence of M. tuberculosis H37Rv. Restriction-digested BAC arrays of M. tuberculosis H37Rv were used in hybridization experiments with radiolabelled M. bovis BCG genomic DNA to reveal the presence of 10 deletions (RD1-RD10) relative to M. tuberculosis. Seven of these regions, RD4-RD10, were also found to be deleted from M. bovis, with the three M. bovis BCG-specific deletions being identical to the RD1-RD3 loci described previously. The distribution of RD4-RD10 in Mycobacterium africanum resembles that of M. tuberculosis more closely than that of M. bovis, whereas an intermediate arrangement was found in Mycobacterium microti, suggesting that the corresponding genes may affect host range and virulence of the various tubercle bacilli. Among the known products encoded by these loci are a copy of the proposed mycobacterial invasin Mce, three phospholipases, several PE, PPE and ESAT-6 proteins, epoxide hydrolase and an insertion sequence. In a complementary approach, direct comparison of BACs uncovered a third class of deletions consisting of two M. tuberculosis H37Rv loci, RvD1 and RvD2, deleted from the genome relative to M. bovis BCG and M. bovis. These deletions affect a further seven genes, including a fourth phospholipase, plcD. In summary, the insertions and deletions described here have important implications for our understanding of the evolution of the tubercle complex.

The bacterial artificial chromosome (BAC) cloning system is capable of stably propagating large, complex DNA inserts in Escherichia coli. As part of the Mycobacterium tuberculosis H37Rv genome sequencing project, a BAC library was constructed in the pBeloBAC11 vector and used for genome mapping, confirmation of sequence assembly, and sequencing. The library contains about 5,000 BAC clones, with inserts ranging in size from 25 to 104 kb, representing theoretically a 70-fold coverage of the M. tuberculosis genome (4.4 Mb). A total of 840 sequences from the T7 and SP6 termini of 420 BACs were determined and compared to those of a partial genomic database. These sequences showed excellent correlation between the estimated sizes and positions of the BAC clones and the sizes and positions of previously sequenced cosmids and the resulting contigs. Many BAC clones represent linking clones between sequenced cosmids, allowing full coverage of the H37Rv chromosome, and they are now being shotgun sequenced in the framework of the H37Rv sequencing project. Also, no chimeric, deleted, or rearranged BAC clones were detected, which was of major importance for the correct mapping and assembly of the H37Rv sequence. The minimal overlapping set contains 68 unique BAC clones and spans the whole H37Rv chromosome with the exception of a single gap of approximately 150 kb. As a postgenomic application, the canonical BAC set was used in a comparative study to reveal chromosomal polymorphisms between M. tuberculosis, M. bovis, and M. bovis BCG Pasteur, and a novel 12.7-kb segment present in M. tuberculosis but absent from M. bovis and M. bovis BCG was characterized. This region contains a set of genes whose products show low similarity to proteins involved in polysaccharide biosynthesis. The H37Rv BAC library therefore provides us with a powerful tool both for the generation and confirmation of sequence data as well as for comparative genomics and other postgenomic applications. It represents a major resource for present and future M. tuberculosis research projects.