The suitability of restriction fragment length polymorphism markers for evaluating genetic diversity among and synteny between mosquito species

Molecular phylogenetic study of Culex quinquefasciatus mosquito from different geographical regions of India using 16S rRNA gene sequences

Culex quinquefasciatus is a major vector of filariasis and various encephalitis in India and worldwide. Vector control remains the most successful strategy for the suppression of mosquito borne diseases. The genetic structure of vector populations in terms of insecticide resistance and susceptibility or refractoriness to infection may possibly vary. To exploit the genetic variability in vector population could pave the path for the alternative strategies in vector management. The sequences of ribosomal RNA molecules have been widely used for such studies. Here, we examined the molecular phylogenetic relationship among the Cx. quinquefasciatus collected from different geographical regions of India, using 16S ribosomal RNA (16S rRNA) gene nucleotide sequences. The distances among the species were measured using Pearson correlation; the Neighbor-Joining (NJ) method was used for the clustering with appropriate bootstrap values using Data Analysis in Molecular Biology and Evolution (DAMBE) software. The results revealed that the populations are genetically diverse. Based on the distance values and the tree topology on the basis of 16S rRNA sequences reflected the clear biogeographical and geoclimatic pattern among the different geographical populations from India.

Genome-based polymorphic microsatellite development and validation in the mosquito Aedes aegypti and application to population genetics in Haiti

Background

Microsatellite markers have proven useful in genetic studies in many organisms, yet microsatellite-based studies of the dengue and yellow fever vector mosquito Aedes aegypti have been limited by the number of assayable and polymorphic loci available, despite multiple independent efforts to identify them. Here we present strategies for efficient identification and development of useful microsatellites with broad coverage across the Aedes aegypti genome, development of multiplex-ready PCR groups of microsatellite loci, and validation of their utility for population analysis with field collections from Haiti.

Results

From 79 putative microsatellite loci representing 31 motifs identified in 42 whole genome sequence supercontig assemblies in the Aedes aegypti genome, 33 microsatellites providing genome-wide coverage amplified as single copy sequences in four lab strains, with a range of 2-6 alleles per locus. The tri-nucleotide motifs represented the majority (51%) of the polymorphic single copy loci, and none of these was located within a putative open reading frame. Seven groups of 4-5 microsatellite loci each were developed for multiplex-ready PCR. Four multiplex-ready groups were used to investigate population genetics of Aedes aegypti populations sampled in Haiti. Of the 23 loci represented in these groups, 20 were polymorphic with a range of 3-24 alleles per locus (mean = 8.75). Allelic polymorphic information content varied from 0.171 to 0.867 (mean = 0.545). Most loci met Hardy-Weinberg expectations across populations and pairwise F_ST comparisons identified significant genetic differentiation between some populations. No evidence for genetic isolation by distance was observed.

Conclusion

Despite limited success in previous reports, we demonstrate that the Aedes aegypti genome is well-populated with single copy, polymorphic microsatellite loci that can be uncovered using the strategy developed here for rapid and efficient screening of genome supercontig assemblies. These loci are suitable for genetic and population studies using multiplex-PCR.

Analysis of the phylogenetic relationship of Anopheles species, subgenus Cellia (Diptera: Culicidae) and using it to define the relationship of morphologically similar species

Studies on the relationship of various vectors and non-vectors of malaria from the evolutionary point of view are important. Use of molecular methods to define phylogeny helps to understand the interrelationship among the members of the anophelines and elucidate the ambiguity that has arisen from improper classification. It could also help to design molecular markers for species differentiation, particularly in those which pose difficulty when classified, based on morphological features. In the present study, the phylogenetic relationships among the species of the anophelines of subgenus Cellia are inferred from the mitochondrial genes COI and COII, the ribosomal RNA gene, in particular the D3 region, and Internal Transcribed Spacer 2 (ITS2) region. The molecular phylogeny obtained in this work matches with that of the classical morphological taxonomy reasonably well, and was useful in properly defining species positions and resolving the ambiguity that normally arises due to morphological taxonomy. The correct arrangement of the various anopheline taxa as per the traditional morphological character-based classification of anophelines was there when we considered the D3 region of 28S rRNA gene and ITS2 region of rDNA. However, the arrangement of the taxa did not match with that of the morphological classification in some aspects, when we considered the COI and COII region of mitochondrial DNA. It may have been due to the variable degree of the rate of evolution of the different genes within the organism. Thus, a proper selection of those particular genes that evolve at the rate that is reflected at the species differentiation level, could help to construct the correct phylogenetic relationship among the anophelines and could be used to correlate with the grouping pattern done from the morphological perspective.

Natural skip oviposition of the mosquito Aedes aegypti indicated by codominant genetic markers

This study examines the use of codominant restriction fragment length polymorphism (RFLP) markers to estimate the number of sibling families found within and among oviposition sites used by the mosquito Aedes aegypti (L) (Diptera: Culicidae). Estimates were made using pairwise relatedness (rxy) calculations based on alleles shared between individuals. Genotypes for eight laboratory mosquito families were determined at six RFLP loci and the observed allele frequencies were used to generate simulated distributions of rxy from full-sibling and unrelated pairs of individuals. The midpoint (mp) between the means of the pairwise rxy distributions was used to discriminate full-sibling families from unrelated families. Clusters of individuals with rxy values higher than the mp value were grouped as putative sibling families. This method was tested by calculating actual rxy for all pairwise comparisons of the known laboratory full-sibling and paternal half-sibling families, followed by upgma cluster analysis to group sibling families. The technique was then used for sibling estimations on wild caught mosquitoes collected at three locations in Trinidad, West Indies. From field populations, 35 families were estimated among 122 individuals tested with an average of 6.2 families per container. Members of 19 predicted families clustered as groups across multiple containers, providing molecular evidence for skip-oviposition behaviour in Ae. aegypti females, whereby individual females oviposit in more than one container.

Linkage map organization of expressed sequence tags and sequence tagged sites in the mosquito, Aedes aegypti

A composite genetic linkage map for the yellow fever mosquito Aedes aegypti was constructed based on restriction fragment length polymorphism (RFLP), single nucleotide polymorphism (SNP) and single strand conformation polymorphism (SSCP) markers. The map consists of 146 marker loci distributed across 205 cM, and includes several morphological mutant marker loci. Most of the genetic markers are derived from random cDNAs or Ae. aegypti genes of known function. A number of markers are derived from random genomic DNAs, including several cloned RAPD-PCR fragments, and also several cDNAs from Drosophila melanogaster. Most of the random cDNAs (80.2%) have high BlastX sequence identities to known genes, with the majority of matches to genes from D. melanogaster. Access to sequence data for all markers will facilitate their continued development for use in high-throughput SNP marker analyses and also provides additional physical anchor points for an anticipated genome sequencing effort.

A comparative genomic analysis of two distant diptera, the fruit fly, Drosophila melanogaster, and the malaria mosquito, Anopheles gambiae

Genome evolution entails changes in the DNA sequence of genes and intergenic regions, changes in gene numbers, and also changes in gene order along the chromosomes. Genes are reshuffled by chromosomal rearrangements such as deletions/insertions, inversions, translocations, and transpositions. Here we report a comparative study of genome organization in the main African malaria vector, Anopheles gambiae, relative to the recently determined sequence of the Drosophila melanogaster genome. The ancestral lines of these two dipteran insects are thought to have separated approximately 250 Myr, a long period that makes this genome comparison especially interesting. Sequence comparisons have identified 113 pairs of putative orthologs of the two species. Chromosomal mapping of orthologous genes reveals that each polytene chromosome arm has a homolog in the other species. Between 41% and 73% of the known orthologous genes remain linked in the respective homologous chromosomal arms, with the remainder translocated to various nonhomologous arms. Within homologous arms, gene order is extensively reshuffled, but a limited degree of conserved local synteny (microsynteny) can be recognized.

Comparative linkage map development and identification of an autosomal locus for insensitive acetylcholinesterase-mediated insecticide resistance in Culex tritaeniorhynchus

A comparative linkage map for Culex tritaeniorhynchus was constructed based on restriction fragment length polymorphism markers using cDNAs from Aedes aegypti. Linear orders of marker loci in Cx. tritaeniorhynchus were identical to Culex pipiens wherein chromosomes 2 and 3 reflect whole-arm rearrangements compared to A. aegypti. However, the sex determination locus in Cx. tritaeniorhynchus maps to chromosome 3, in contrast to Cx. pipiens and Ae. aegypti where it is located on chromosome 1. Our results indicate that insensitive acetylcholinesterase (AChE)-mediated organophosphate resistance is controlled by a single major gene (AChE) on chromosome 2, while the AChE structural gene (Ace) is located on chromosome 1. No evidence for a second Ace gene was observed, even under very low stringency hybridization conditions.

SSCP analysis of cDNA markers provides a dense linkage map of the Aedes aegypti genome

An intensive linkage map of the yellow fever mosquito, Aedes aegypti, was constructed using single-strand conformation polymorphism (SSCP) analysis of cDNA markers to identify single nucleotide polymorphisms (SNPs). A total of 94 A. aegypti cDNAs were downloaded from GenBank and primers were designed to amplify fragments <500 bp in size. These primer pairs amplified 94 loci, 57 (61%) of which segregated in a single F(1) intercross family among 83 F(2) progeny. This allowed us to produce a dense linkage map of one marker every 2 cM distributed over a total length of 134 cM. Many A. aegypti cDNAs were highly similar to genes in the Drosophila melanogaster genome project. Comparative linkage analysis revealed areas of synteny between the two species. SNP polymorphisms are abundant in A. aegypti genes and should prove useful in both population genetics and mapping studies.

Mannose phosphate isomerase isoenzymes in Plutella xylostella support common genetic bases of resistance to Bacillus thuringiensis toxins in Llpidopteran species

A strong correlation between two mannose phosphate isomerase (MPI) isoenzymes and resistance to Cry1A toxins from Bacillus thuringiensis has been found in a Plutella xylostella population. MPI linkage to Cry1A resistance had previously been reported for a Heliothis virescens population. The fact that the two populations share similar biochemical, genetic, and cross-resistance profiles of resistance suggests the occurrence of homologous resistance loci in both species.

Genetic and physical mapping in mosquitoes: molecular approaches

The genetic background of individual mosquito species and populations within those species influences the transmission of mosquito-borne pathogens to humans. Technical advances in contemporary genomics are contributing significantly to the detailed genetic analysis of this mosquito-pathogen interaction as well as all other aspects of mosquito biology, ecology, and evolution. A variety of DNA-based marker types are being used to develop genetic maps for a number of mosquito species. Complex phenotypic traits such as vector competence are being dissected into their discrete genetic components, with the intention of eventually using this information to develop new methods to prevent disease transmission. Both genetic- and physical-mapping techniques are being used to define and compare genome architecture among and within mosquito species. The integration of genetic- and physical-map information is providing a sound framework for map-based positional cloning of target genes of interest. This review focuses on advances in genome-based analysis and their specific applications to mosquitoes.

Sequence analysis of mitochondrial 16S ribosomal RNA gene fragment from seven mosquito species

Mosquitoes are vectors for the transmission of many human pathogens that include viruses, nematodes and protozoa. For the understanding of their vectorial capacity, identification of disease carrying and refractory strains is essential. Recently, molecular taxonomic techniques have been utilized for this purpose. Sequence analysis of the mitochondrial 16S rRNA gene has been used for molecular taxonomy in many insects. In this paper, we have analysed a 450 bp hypervariable region of the mitochondrial 16S rRNA gene in three major genera of mosquitoes, Aedes, Anopheles and Culex. The sequence was found to be unusually A+T rich and in substitutions the rate of transversions was higher than the transition rate. A phylogenetic tree was constructed with these sequences. An interesting feature of the sequences was a stretch of Ts that distinguished between Ae-des and Culex on the one hand, and Anopheles on the other. This is the first report of mitochondrial rRNA sequences from these medically important genera of mosquitoes.

Evidence for genetic hitchhiking effect associated with insecticide resistance in Aedes aegypti

Information on genetic variation within and between populations is critical for understanding the evolutionary history of mosquito populations and disease epidemiology. Previous studies with Drosophila suggest that genetic variation of selectively neutral loci in a large fraction of genome may be constrained by fixation of advantageous mutations associated with hitchhiking effect. This study examined restriction fragment length polymorphisms of four natural Aedes aegypti mosquito populations from Trinidad and Tobago, at 16 loci. These populations have been subjected to organophosphate (OP) insecticide treatments for more than two decades, while dichlor-diphenyltrichlor (DDT) was the insecticide of choice prior to this period. We predicted that genes closely linked to the OP target loci would exhibit reduced genetic variation as a result of the hitchhiking effect associated with intensive OP insecticide selection. We also predicted that genetic variability of the genes conferring resistance to DDT and loci near the target site would be similar to other unlinked loci. As predicted, reduced genetic variation was found for loci in the general chromosomal region of a putative OP target site, and these loci generally exhibited larger F(ST) values than other random loci. In contrast, the gene conferring resistance to DDT and its linked loci show polymorphisms and genetic differentiation similar to other random loci. The reduced genetic variability and apparent gene deletion in some regions of chromosome 1 likely reflect the hitchhiking effect associated with OP insecticide selection.

Development of a comparative genetic linkage map for Armigeres subalbatus using Aedes aegypti RFLP markers

One of the causative agents of lympahtic filariasis is the nematode parasite Brugia malayi that requires a competent mosquito vector for its development and transmission. Armigeres subalbatus mosquitoes rapidly destroy invading B. malayi microfilariae via a defense response known as melanotic encapsulation. We have constructed a genetic linkage map for this mosquito species using RFLP markers from Aedes aegypti. This heterologous approach was possible because of the conserved nature of the coding sequences used as markers and provided an experimental framework to evaluate the hypothesis that linkage and gene order are conserved between these mosquito species. Of the 56 Ae. aegypti markers tested, 77% hybridize to genomic DNA digests of Ar. subalbatus under stringent conditions, with 53% of these demonstrating strain-specific polymorphisms. Twenty-six Ae. aegypti markers have been mapped using an F2- segregating Ar. subalbatus population derived from a cross of strains originating in Japan and Malaysia. Linear order of these marker loci is highly conserved between the two species. Only 1 of these markers, LF92, was not linked in the manner predicted by the Ae. aegypti map. In addition, the autosomal sex-determination locus that occurs in linkage group 1 in Ae. aegypti resides in group 3 in Ar. subalbatus. The Ar. subalbatus map provides a basic genetic context that can be utilized in further genetic studies to clarify the genetic basis of parasite resistance in this mosquito and is a necessary precursor to the identification of genome regions that carry genes that determine the encapsulation phenotype. [The composite map and sequence database information for Ae. aegypti markers can be retrieved directly from the Ae. aegypti Genome Database through the World Wide Web: http://klab.agsci.colostate.edu.]

Comparative linkage maps for the mosquitoes, Aedes albopictus and Ae. aegypti, based on common RFLP loci

Aedes albopictus and Aedes aegypti are members of the mosquito family Culicidae and share a haploid chromosome complement of three. Although a genetic linkage map based on restriction fragment length polymorphism (RFLP), markers exists for Ae. aegypti, the extent of synteny and linkage order conservation between the two species was unknown. A comparative linkage map for Ae. albopictus was constructed based mainly on cDNA clones from Ae. aegypti. Nearly all Ae. aegypti probes hybridized to Ae. albopictus DNA at high stringency. For eighteen RFLP markers tested, the linkage group and linear order appears to be identical for the two species. 78% of the loci tested exhibited significant deviations from the expected segregation ratio in at least one of the test crosses. An excess of heterozygote genotypes was recovered with most loci. This probably reflects the effects of lethal loci on survival of F2 progeny homozygous for the parental genotypes. These results demonstrate that comparative linkage maps based on common DNA markers provide a basis for rapidly developing linkage maps for various mosquito species, and the opportunity to examine the significance and function of orthologous quantitative trait loci associated with mosquito vector competence for disease transmission.

Applications of molecular marker analysis to mosquito vector competence

A rapidly expanding cadre of molecular biology techniques is being developed for human and plant genetics, including development of the technology to identify large numbers of genetic markers and to evaluate these markers relative to phenotypic observations. In this review, David Severson discusses applications of these techniques for the analysis of mosquito vector competence.