Evolution of the mitochondrial DNA control region in the Anopheles gambiae complex

Mitochondrial Variation in Anopheles gambiae and Anopheles coluzzii: Phylogeographic Legacy and Mitonuclear Associations With Metabolic Resistance to Pathogens and Insecticides

Mitochondrial DNA has been a popular marker in phylogeography, phylogeny, and molecular ecology, but its complex evolution is increasingly recognized. Here, we investigated mitochondrial DNA variation in Anopheles gambiae and Anopheles coluzzii, in relation to other species in the Anopheles gambiae complex, by assembling the mitogenomes of 1,219 mosquitoes across Africa. The mitochondrial DNA phylogeny of the Anopheles gambiae complex was consistent with previously reported highly reticulated evolutionary history, revealing important discordances with the species tree. The three most widespread species (An. gambiae, An. coluzzii, and Anopheles arabiensis), known for extensive historical introgression, could not be discriminated based on mitogenomes. Furthermore, a monophyletic clustering of the three saltwater-tolerant species (Anopheles merus, Anopheles melas, and Anopheles bwambae) in the Anopheles gambiae complex also suggested that introgression and possibly selection shaped mitochondrial DNA evolution. Mitochondrial DNA variation in An. gambiae and An. coluzzii across Africa revealed significant partitioning among populations and species. A peculiar mitochondrial DNA lineage found predominantly in An. coluzzii and in the hybrid taxon of the African "far-west" exhibited divergence comparable to the interspecies divergence in the Anopheles gambiae complex, with a geographic distribution matching closely An. coluzzii's geographic range. This phylogeographic relict of the An. coluzzii and An. gambiae split was associated with population and species structure, but not with the rare Wolbachia occurrence. The lineage was significantly associated with single nucleotide polymorphisms in the nuclear genome, particularly in genes associated with pathogen and insecticide resistance. These findings underline potential mitonuclear coevolution history and the role played by mitochondria in shaping metabolic responses to pathogens and insecticides in Anopheles.

In silico identification of multiple conserved motifs within the control region of Culicidae mitogenomes

Mosquitoes are important vectors for human and animal diseases. Genetic markers, like the mitochondrial COI gene, can facilitate the taxonomic classification of disease vectors, vector-borne disease surveillance, and prevention. Within the control region (CR) of the mitochondrial genome, there exists a highly variable and poorly studied non-coding AT-rich area that contains the origin of replication. Although the CR hypervariable region has been used for species differentiation of some animals, few studies have investigated the mosquito CR. In this study, we analyze the mosquito mitogenome CR sequences from 125 species and 17 genera. We discovered four conserved motifs located 80 to 230 bp upstream of the 12S rRNA gene. Two of these motifs were found within all 392 Anopheles (An.) CR sequences while the other two motifs were identified in all 37 Culex (Cx.) CR sequences. However, only 3 of the 304 non-Culicidae Dipteran mitogenome CR sequences contained these motifs. Interestingly, the short motif found in all 37 Culex sequences had poly-A and poly-T stretch of similar length that is predicted to form a stable hairpin. We show that supervised learning using the frequency chaos game representation of the CR can be used to differentiate mosquito genera from their dipteran relatives.

Sequence analysis of mtDNA COI barcode region revealed three haplotypes within Culex pipiens assemblage

Members of the Culex (Culex) pipiens assemblage are known vectors of deadly encephalitides, periodic filariasis, and West Nile virus throughout the world. However, members of this assemblage are morphologically indistinguishable or hard to distinguish and play distinct roles in transmission of the diseases. The current study aimed to provide further evidence on utility of the two most popular nuclear (ITS2-rDNA) and mitochondrial (COI barcode region) genetic markers to identify members of the assemblage. Culex pipiens assemblage specimens from different climate zones of Iran were collected and identified to species level based on morphological characteristics. Nucleotide sequences of the loci for the specimens plus available data in the GenBank were analyzed to find species specific genetic structures useful for diagnosis purposes. ITS2 region was highly divergent within species or populations suggesting lack of consistency as a reliable molecular marker. In contrast, sequence analysis of 710 bp of COI gene revealed three fixed haplotypes named here "C, T, H" within the assemblage which can be distinguished by HaeIII and AluI enzymes. There were a correlation between the haplotypes and the world climate regions, where the haplotypes H/T and C are present mainly in temperate and tropical regions of the world, respectively. In the New world, Australia, and Japan only haplotype H is found. In conjunction between tropical and temperate regions such Iran, China, and Turkey, a mix of C/H or C/H/T are present. Although, the haplotypes are not strictly species-specific, however, Cx. quinquefasciatus was mainly of haplotype C. Due to the lack of mating barrier and questionable taxonomic situation of the complex members, the mentioned haplotypes in combination with other morphological and molecular characters might be used to address the genetic structure of the studied populations.

Population Genetic Structure of the Malaria Vector Anopheles sinensis (Diptera: Culicidae) Sensu Stricto and Evidence for Possible Introgression in the Republic of Korea

Anopheles sinensis Wiedemann sensu stricto (s.s.) is a dominant mosquito and considered a secondary malaria vector in the Republic of Korea (ROK). Despite the potential significance for malaria control, population genetics studies have been conducted using only mitochondrial DNA (mtDNA), and studies of the genetics of hybridization have never been attempted. In this study, 346 specimens from 23 localities were subject to experiments. Among them, 305 An. sinensis s.s. specimens from 20 localities were used for mtDNA analysis, and 346 specimens comprising 341 An. sinensis s.s. from 22 localities and five Anopheles kleini Rueda from one locality were examined in the microsatellite study. Neighbor-joining analysis of pairwise FST and RST based on microsatellite results showed that the populations are divided into two groups, as did the mtDNA results. However, the Bayesian analysis and factorial correspondence analysis plots showed three distinct clusters. Among the mtDNA and microsatellite results, only microsatellites represented small but positive and significant isolation-by-distance patterns. Both molecular markers show the Taebaek and Sobaek Mountain ranges as barriers between the northern and southern parts of the ROK. The newly recognized third group suggests possible introgressive hybridization of An. sinensis s.s. with closely related species. The slightly different composition of populations in each group based on different markers is probably because of different population dynamics in each group. These results imply that there is restricted gene flow of epidemiologically important malaria-related genes between the northern and southern parts of the ROK.

Mitochondrial genomes and comparative analyses of Culex camposi, Culex coronator, Culex usquatus and Culex usquatissimus (Diptera:Culicidae), members of the coronator group

BACKGROUND

The Coronator Group currently encompasses six morphologically similar species (Culex camposi Dyar, Culex coronator Dyar and Knab, Culex covagarciai Forattini, Culex usquatus Dyar, Culex usquatissimus Dyar, and Culex ousqua Dyar). Culex coronator has been incriminated as a potential vector of West Nile Virus (WNV), Saint Louis Encephalitis Virus (SLEV), and Venezuelan Equine Encephalitis Virus (VEEV). The complete mitochondrial genome of Cx. coronator, Cx. usquatus, Cx.usquatissimus, and Cx. camposi was sequenced, annotated, and analyzed to provide genetic information about these species.

RESULTS

The mitochondrial genomes of Cx. coronator, Cx. usquatus, Cx.usquatissimus, and Cx. camposi varied from 15,573 base pairs in Cx. usquatus to 15,576 in Cx. coronator. They contained 37 genes (13 protein-encoding genes, 2 rRNA genes, and 22 tRNA genes) and the AT-rich control region. Comparative analyses of the 37 genes demonstrated the mitochondrial genomes to be composed of variable and conserved genes. Despite the small size, the ATP8, ATP6 plus NADH5 protein-encoding genes were polymorphic, whereas tRNAs and rRNAs were conserved. The control region contained some poly-T stretch. The Bayesian phylogenetic tree corroborated that both the Coronator Group and the Culex pipens complex are monophyletic taxa.

CONCLUSIONS

The mitochondrial genomes of Cx. coronator, Cx. usquatus, Cx. usquatissimus and Cx. camposi share the same gene composition and arrangement features that match to those reported for most Culicidae species. They are composed of the same 37 genes and the AT-rich control region, which contains poly-T stretches that may be involved in the functional role of the mitochondrial genome. Taken together, results of the dN/dS ratios, the sliding window analyses and the Bayesian phylogenetic analyses suggest that ATP6, ATP8 and NADH5 are promising genes to be employed in phylogenetic studies involving species of the Coronator Group, and probably other species groups of the subgenus Culex. Bayesian topology corroborated the morphological hypothesis of the Coronator Group as monophyletic lineage within the subgenus Culex.

The complete mitochondrial genomes of three grasshoppers, Asiotmethis zacharjini, Filchnerella helanshanensis and Pseudotmethis rubimarginis (Orthoptera: Pamphagidae)

The complete mitogenomes of Asiotmethis zacharjini, Filchnerella helanshanensis and Pseudotmethis rubimarginis are 15,660 bp, 15,657 bp and 15,661 bp in size, respectively. All three mitogenomes contain a standard set of 13 protein - coding genes, 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and an A+T-rich region in the same order as those of the other analysed caeliferan species, including the rearrangement of trnAsp and trnLys. The putative initiation codon for the cox1 gene in the three species is CCG. The long polythymine stretch (T-stretch) in the A+T-rich region of the three species is not adjacent to the trnIle but inside the stem-loop sequence in the majority strand. The mitogenomes of F. helanshanensis and P. rubimarginis have higher overall similarities. The characterization of the three mitogenomes will enrich our knowledge on the Pamphagidae mitogenome. The phylogenetic analyses indicated that within the Caelifera, Pyrgomorphoidea is a sister group to Acridoidea. The species from the Pamphagidae form a monophyletic group, as is the case for Acrididae. Furthermore, the two families cluster as sister groups, supporting the monophyly of Acridoidea. The relationships among eight acridid subfamilies were (Cyrtacanthacridinae+(Calliptaminae+(Catantopinae+(Oxyinae+(Melanopline+(Acridinae+(Oedipodinae+Gomphocerinae).

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.

Malaria causes more than one million deaths every year, mostly among children in Sub-Saharan Africa. Anopheles mosquitoes are exclusive vectors of human malaria. Many malaria vectors belong to species complexes, and members within these complexes can vary significantly in their ecological adaptations and ability to transmit the parasite. To better understand evolution of epidemiologically important traits, we studied relationships among nonvector and vector species of the African Anopheles gambiae complex. We analyzed gene orders at genomic regions where evolutionary breaks of chromosomal inversions occurred in members of the complex and compared them with gene orders in species outside the complex. This approach allowed us to identify ancient and recent gene orders for three chromosomal inversions. Surprisingly, the more ancestral chromosomal arrangements were found in mosquito species that are vectors of human malaria, while the more derived arrangements were found in both nonvectors and vectors. Our finding strongly suggests that the increased ability to transmit human malaria originated repeatedly during the recent evolution of these African mosquitoes. This knowledge can be used to identify specific genetic changes associated with the human blood choice and ecological adaptations.

Evolution of Anopheles gambiae in Relation to Humans and Malaria

The closely related and morphologically indistinguishable mosquito species in the Afrotropical Anopheles gambiae complex differ dramatically in their contribution to malaria transmission, ranging from major vectors through minor or locally important vectors and nonvectors. Radiation of the A. gambiae complex and ongoing diversification within its nominal species appears to be a product of recent and rapid adaptation to environmental heterogeneities, notably those of anthropogenic origin. Polytene chromosome and genomic analyses suggest that paracentric chromosomal inversions and possibly other low-recombination regions have played instrumental roles in this process by facilitating ecotypic differentiation both within and across semipermeable species boundaries. Forthcoming complete genome sequences from several members of the A. gambiae complex will provide powerful tools to accelerate ongoing investigation of how genetic diversification of populations and species has shaped behavioral and physiological traits, such as vector competence, that bear on vectorial importance.

Geographically widespread swordfish barcode stock identification: a case study of its application

Background

The swordfish (Xiphias gladius) is a cosmopolitan large pelagic fish inhabiting tempered and tropical waters and it is a target species for fisheries all around the world. The present study investigated the ability of COI barcoding to reliably identify swordfish and particularly specific stocks of this commercially important species.

Methodology

We applied the classical DNA barcoding technology, upon a 682 bp segment of COI, and compared swordfish sequences from different geographical sources (Atlantic, Indian Oceans and Mediterranean Sea). The sequences of the 5′ hyper-variable fragment of the control region (5′dloop), were also used to validate the efficacy of COI as a stock-specific marker.

Case Report

This information was successfully applied to the discrimination of unknown samples from the market, detecting in some cases mislabeled seafood products.

Conclusions

The NJ distance-based phenogram (K2P model) obtained with COI sequences allowed us to correlate the swordfish haplotypes to the different geographical stocks. Similar results were obtained with 5′dloop. Our preliminary data in swordfish Xiphias gladius confirm that Cytochrome Oxidase I can be proposed as an efficient species-specific marker that has also the potential to assign geographical provenance. This information might speed the samples analysis in commercial application of barcoding.

Complete Mitochondrial Genome Sequence of Acrida cinerea (Acrididae: Orthoptera) and Comparative Analysis of Mitochondrial Genomes in Orthoptera

The complete 15,599-bp mitogenome of Acrida cinerea was determined and compared with that of the other 20 orthopterans. It displays characteristic gene content, genome organization, nucleotide composition, and codon usage found in other Caelifera mitogenomes. Comparison of 21 orthopteran sequences revealed that the tRNAs encoded by the H-strand appear more conserved than those by the L-stand. All tRNAs form the typical clover-leaf structure except trnS (agn), and most of the size variation among tRNAs stemmed from the length variation in the arm and loop of TΨC and the loop of DHU. The derived secondary structure models of the rrnS and rrnL from 21 orthoptera species closely resemble those from other insects on CRW except a considerably enlarged loop of helix 1399 of rrnS in Caelifera, which is a potentially autapomorphy of Caelifera. In the A+T-rich region, tandem repeats are not only conserved in the closely related mitogenome but also share some conserved motifs in the same subfamily. A stem-loop structure, 16 bp or longer, is likely to be involved in replication initiation in Caelifera and Grylloidea. A long T-stretch (>17 bp) with conserved stem-loop structure next to rrnS on the H-strand, bounded by a purine at either end, exists in the three species from Tettigoniidae.

Complete mtDNA genomes of Anopheles darlingi and an approach to anopheline divergence time

BACKGROUND

The complete sequences of the mitochondrial genomes (mtDNA) of members of the northern and southern genotypes of Anopheles (Nyssorhynchus) darlingi were used for comparative studies to estimate the time to the most recent common ancestor for modern anophelines, to evaluate differentiation within this taxon, and to seek evidence of incipient speciation.

METHODS

The mtDNAs were sequenced from mosquitoes from Belize and Brazil and comparative analyses of structure and base composition, among others, were performed. A maximum likelihood approach linked with phylogenetic information was employed to detect evidence of selection and a Bayesian approach was used to date the split between the subgenus Nyssorhynchus and other Anopheles subgenera.

RESULTS

The comparison of mtDNA sequences within the Anopheles darlingi taxon does not provide sufficient resolution to establish different units of speciation within the species. In addition, no evidence of positive selection in any protein-coding gene of the mtDNA was detected, and purifying selection likely is the basis for this lack of diversity. Bayesian analysis supports the conclusion that the most recent ancestor of Nyssorhynchus and Anopheles+Cellia was extant ~94 million years ago.

CONCLUSION

Analyses of mtDNA genomes of Anopheles darlingi do not provide support for speciation in the taxon. The dates estimated for divergence among the anopheline groups tested is in agreement with the geological split of western Gondwana (95 mya), and provides additional support for explaining the absence of Cellia in the New World, and Nyssorhynchus in the Afro-Eurasian continents.

Two different routes of colonization of Aedes aegypti in Argentina from neighboring countries

Aedes aegypti L. (Diptera, Culicidae) is the main vector of dengue and yellow fever. In Argentina, the species was apparently eradicated approximately in 1964; by 1986, it was reintroduced. To identify different gene pools in geographical populations of the species and to ascertain the possible routes of colonization, we analyzed the diversity of mitochondrial DNA haplotypes in 572 specimens from Argentina and neighboring countries. We found that the restriction fragment length polymorphism-polymerase chain reaction screening of a large DNA fragment including the A+T-rich region was the best strategy to reconstruct the colonization pattern ofAe. aegypti in Argentina. Twenty haplotypes were recognized; levels of genetic similarity varied among populations from different geographical locations. The haplotype network constructed on the basis of genetic distances showed three well differentiated groups. Two of them exhibited a well defined spatial distribution and populations in these groups presented an isolation-by-distance pattern. The persistence of relictual populations after the last eradication campaigns would explain the high levels of haplotype diversity and the presence of exclusive haplotypes in urban centers from northwestern Argentina. Eastern Argentine populations showed one prevalent haplotype, also predominant in Brazil and Paraguay. Our results highlight the need for efficient surveys and control campaigns, given the strong effect of land trade on genetic exchange among mosquito populations from Argentina and neighboring countries where dengue is endemic.

Structure and evolution of the mitochondrial control region of the pollen beetleMeligethes thalassophilus(Coleoptera: Nitidulidae)

The organization of the mitochondrial DNA (mtDNA) control region (CR) of the pollen beetle Meligethes thalassophilus is described. This mtDNA CR represents the longest sequenced for beetles so far, since the entire nucleotide sequence ranges from ~5000 to ~5500 bp. The CR of M. thalassophilus is organized in three distinct domains: a conserved domain near the tRNAIlegene, a variable domain flanking the 12S rRNA gene, and a relatively large central tandem array made up of a variable number of ~170 bp repeats that is responsible for the intraspecific length variation observed. Like other CRs found in insects, the M. thalassophilus CR contains two long homopolymeric runs that may be involved in mtDNA replication. Furthermore, conserved stem-and-loop structures in the repetitive domain were identified and their possible role in generating length variation is examined. Intraspecific comparison of the tandem repeat elements of M. thalassophilus suggests mechanisms of concerted evolution leading to homogenization of the repetitive region. The utility of such an array of tandem repeats as a genetic marker for assessing population-level variability and evolutionary relationships among populations is discussed. Finally, the technical difficulties found in isolating the mtDNA CR in beetles are remarked upon.

Population structure and gene flow of Anopheles farauti s.s. (Diptera: Culicidae) among ten sites on five islands of Vanuatu: implications for malaria control

The Anopheles punctulatus (Diptera: Culicidae) group is the main vector for malaria and Bancroftian filariasis in Vanuatu. Anopheles larvae were collected from 10 localities on five islands of Vanuatu during the 2004 dry season for species identification as well as for estimating population structure and gene flow within and among islands. Species identification was determined using polymerase chain reaction-restriction fragment length polymorphism analysis of the internal transcribed spacer 2 region. Population structure and gene flow were examined by sequencing a portion of the ND4/ND5 region of the mitochondrial genome. Only one species of the An. punctulatus group, An. farauti s.s., was identified, consistent with previous studies in Vanuatu. A nonrandom distribution of An. farauti s.s. lineages was observed with one cosmopolitan lineage shared by eight sites on all five islands and a preponderance of island-specific lineages (36/40), indicating the introduction of a single main lineage into Vanuatu followed by dispersal, diversification, and limited lineage exchange between islands. Network analysis suggests a possible second introduction of An. farauti s.s. into the northern islands of Gaua and Malekula. Gene flow was high on three of the five islands, whereas Tanna and Santo have significant population structure. Among islands, gene flow was limited, indicating active mosquito dispersal only over short distances and a paucity of passive human-mediated dispersal over long distances. Minimal risk of active dispersal among these islands indicates that vector control can be effectively initiated at the island level within the archipelago of Vanuatu.

The mitochondrial DNA control region of Muscidae flies: evolution and structural conservation in a dipteran context

The structure and evolution of the mtDNA control region (CR) and its flanking genes in economically important dipterans from the family Muscidae (Brachycera: Calyptratae), Haematobia irritans, Musca domestica, Atherigona orientalis, and Stomoxys calcitrans are presented in this paper, along with the description of short noncoding intergenic regions possibly related to CR flanking sequences in Stomoxys calcitrans and Ophyra aenescens mtDNAs (ScIR and OaIR, respectively). S. calcitrans showed a large CR with an approximately 550-bp element tandemly repeated and a duplicated tRNA(Ile) gene. The characterization of H. irritans, M. domestica, A. orientalis, and S. calcitrans CR sequences led to the identification of seven conserved sequence blocks homologous to the elements previously described for Calliphoridae and Oestridae species (Brachycera: Calyptratae). Comparative analysis with Drosophila species (Brachycera: Acalyptratae) revealed four conserved regions. The putative functional roles of the conserved elements in the regulation of replication and transcription processes are addressed. The characterization of the structural organization of the mitochondrial genome CR demonstrates the plasticity of the mtDNA molecule in family Muscidae.

The complete mitochondrial genome sequence of the Mormon cricket (Anabrus simplex: Tettigoniidae: Orthoptera) and an analysis of control region variability

The Anabrus simplex is a swarming plague orthopteran found in western North America. The genome is 15 766 bp in length and genome organization follows the ancestral insect gene arrangement. atp6 lacked any readily identifiable stop codon. Examination of mRNA secondary structure for this gene suggested a stem/loop-mediated mRNA post-transcriptional processing to liberate a mature atp6 mRNA with a complete stop codon produced by polyadenylation. Comparison of similar protein with protein gene boundaries in other insect species reveal a general mechanism for mRNA excision and provide further supporting evidence for post-transcriptional mRNA processing in mitochondrial genomes. The A + T-rich region, or control region, was sequenced for 55 A. simplex individuals from 12 different populations. Variance studies between these individuals show that the A + T-rich region contains significant phylogenetic signal to be used in population studies.

Structural organization of the mitochondrial DNA control region in Aedes aegypti

The complete A+T - rich region of Aedes aegypti mitochondrial DNA has been cloned and sequenced. In Argentinean populations of the species, a polymorphism in the length of the amplified fragment was observed. Nucleotide sequence comparison of the shortest and longest A+T - rich amplified fragments detected revealed the presence of 2 types of tandemly repeated blocks. The size variation observed in natural populations is mainly due to the presence of a variable number of a 181 bp tandem repeat unit, located toward the 12S rRNA gene end. The size of the longest A+T - rich region was of 2070 bp, representing the largest control sequence reported for any mosquito species. Few relevant short blocks of primary-sequence similarity conserved in the control region of mosquitoes and other insects were detected scattered throughout the whole region. Five putative stem-loop secondary structures were found, one of them flanked by conserved sequences described in other insects. Our results suggest that there are no universal models of structure-function relations in the control region of insect mtDNA. In addition, we identified a short A+T - rich variable segment in the Ae. aegyti control region that would be suitable for population genetic studies.

The hypervariable domain of the mitochondrial control region in Atlantic spiny lobsters and its potential as a marker for investigating phylogeographic structuring

Atlantic spiny lobsters support major fisheries in northeastern Brazilian waters and in the Caribbean Sea. To avoid reduction in diversity and elimination of distinct stocks, understanding their population dynamics, including structuring of populations and genetic diversity, is critical. We here explore the potential of using the hypervariable domain in the control region of the mitochondrial DNA as a genetic marker to characterize population subdivision in spiny lobsters, using Panulirus argus as the species model. The primers designed on the neighboring conserved genes have amplified the entire control region (approx. 780 bases) of P. argus and other closely related species. Average nucleotide and haplotype diversity within P. argus were found to be high, and population structuring was hypothesized. The data suggest a division of P. argus into genetically different phylogeographic groups. The hypervariable domain seems to be useful for determining genetic differentiation of geographically distinct stocks of P. argus and other Atlantic spiny lobsters.

The classification of genus Anopheles (Diptera: Culicidae): a working hypothesis of phylogenetic relationships

The internal classification of genus Anopheles is updated to reflect taxonomic actions published since the classification was last reviewed in 1994. Both formal and informal taxa are included. The classification is intended to aid researchers and students who are interested in analysing species relationships, making group comparisons and testing phylogenetic hypotheses. The genus includes 444 formally named and 40 provisionally designated extant species divided between six subgenera: Anopheles, Cellia, Kerteszia, Lophopodomyia, Nyssorhynchus and Stethomyia. Subgenera Anopheles, Cellia and Nyssorhynchus are subdivided hierarchically into nested informal groups of morphologically similar species that are believed to represent monophyletic lineages based on morphological similarity. Changes to the classification include additional species, eliminated species and changes to the hierarchical organization and composition of supraspecific groups, some as a result of molecular studies.

The utility of the neglected mitochondrial control region for evolutionary studies in lepidoptera (insecta)

The insect mitochondrial control region (=AT-rich-region) is a rarely used genetic marker in phylogeographic studies and population genetic surveys. Reasons for this are that the high AT content and the presence of tandem repeats and indels pose technical and analytical problems. We provide a new pair of primers and the first taxonomically wide-scale description of control region (CR) structure in an insect order after sequencing it in 31 lepidopteran species. We assessed levels of variation occurring in the CR and cytochrome oxidase I (COI) by sequencing and comparison. Intrapopulation analyses in five species of butterflies showed that CR was more variable than COI. Interpopulation variation from three populations of Erebia triaria and E. palarica was slightly lower in the CR than COI with regard to single-nucleotide polymorphisms, but the results were concordant between both markers and highly congruent with regard to population differentiation. Using 15 species of Satyrinae we found that the CR has the same, or stronger, phylogenetic resolution as COI. Our results indicate that the CR can be of importance in addition to COI in population genetic studies. Alignments for the whole CR are direct and unambiguous at the intraspecific level. Indels show phylogenetic signal, but make this marker more complex to use than COI for higher phylogenetic analyses. Nevertheless, alignments at the generic level are straightforward over one part of the CR. The combination CR+COI appears to be a very promising phylogenetic tool to resolve fast-evolving species-level phylogenies.

Revisiting the role of introgression vs shared ancestral polymorphisms as key processes shaping genetic diversity in the recently separated sibling species of the Anopheles gambiae complex

The role of interspecific hybridisation in the evolution of pest species is poorly understood. In mosquito disease vectors this is of particular importance due to the evolution of insecticide resistance and the proposed release of transgenic strains that are refractory to the malaria parasite. In this study, we apply population genetic methods in a novel manner to determine whether mitochondrial DNA sequences have introgressed between the closely related African malaria vectors Anopheles gambiae and A. arabiensis. Our results suggest that speciation was geologically recent and ancestral haplotypes at the ND5 locus are retained in both species. In addition, comparing haplotype frequencies in allopatric and sympatric populations, suggest locale specific unidirectional introgression of mitochondria from A. arabiensis into A. gambiae.

Semipermeable species boundaries between Anopheles gambiae and Anopheles arabiensis: evidence from multilocus DNA sequence variation

Attempts to reconstruct the phylogenetic history of the Anopheles gambiae cryptic species complex have yielded strongly conflicting results. In particular, An. gambiae, the primary African malaria vector, is variously placed as a sister taxon to either Anopheles arabiensis or Anopheles merus. The recent divergence times for members of this complex complicate phylogenetic analysis, making it difficult to unambiguously implicate interspecific gene flow, versus retained ancestral polymorphism, as the source of conflict. Using sequences at four unlinked loci, which were determined from multiple specimens within each of five species in the complex, we found contrasting patterns of sequence divergence between the X chromosome and the autosomes. The isolation model of speciation assumes a lack of gene flow between species since their separation. This model could not be rejected for An. gambiae and An. arabiensis, although the data fit the model poorly. On the other hand, evidence from gene trees supports genetic introgression of chromosome 2 inversions between An. gambiae and An. arabiensis, and also points to more broad scale genetic exchange of autosomal sequences between this species pair. That such exchange has been relatively recent is suggested not only by the lack of fixed differences at three autosomal loci but also by the sharing of full haplotypes at two of the three loci, which is in contrast to several fixed differences and considerably deeper divergence on the X. The proposed acquisition by An. gambiae of sequences from the more arid-adapted An. arabiensis may have contributed to the spread and ecological dominance of this malaria vector.

Genetic markers for study of the anopheline vectors of human malaria

Human malaria is truly a disease of global proportions and is one of the most broadly distributed vector-borne infections. Anopheline mosquitoes are the exclusive vectors of human malaria. A handful of species predominate as the most notorious malaria vectors, but the species and forms involved in the transmission of human malaria world-wide are incredibly diverse. Many of the anophelines that vector malaria exist as members of species complexes that often contain vector and non-vector species. Additionally, single anopheline species often exhibit significant heterogeneity across the species' range. This phenotypic and genotypic plasticity exacerbates the difficulties in identification of vector populations and implementation of effective surveillance and control strategies. Polytene chromosome investigations were among the first to provide researchers with tangible genetic markers that could be used to differentiate between what are now recognised as species and chromosomal forms of anopheline mosquitoes. The advent of the polymerase chain reaction gave access to the molecular genetics of genomes and the techniques that followed have facilitated investigation of the genetics of individual specimens or population size samples. The variety and number of genetic markers available for the study of malaria vectors has literally exploded in the last 10 years. Markers have expanded from the 'traditional tools' to include a vast array of molecular markers. Contemporary markers range from what are now referred to as 'classical genetic markers' to methods used to detect and identify single nucleotide polymorphisms and finally to highly polymorphic markers. One of the greatest advantages of this wide variety of genetic markers is that researchers may choose to utilise any combination of markers or techniques to address multifaceted questions relating to malaria transmission. These molecular markers have proven useful in a wide variety of applications including molecular taxonomy, evolutionary systematics, population genetics, genetic mapping, and investigation of defined phenotypes.

Arrangement and structural conservation of the mitochondrial control region of two species of Plecoptera: utility of tandem repeat-containing regions in studies of population genetics and evolutionary history

Low levels of primary sequence similarity across insect taxa have led to the suggestion of conserved structural elements in the insect mitochondrial control region. Our aim was to determine whether previously described motifs and secondary structures exist in stoneflies (Plecoptera). Several motifs and structural elements previously described in Orthoptera and Diptera were found, including a conserved 'hairpin' structure that may play a role in the initiation of mtDNA replication. The repeat region had the highest percentage similarity, lowest A-T content and highest transition to transversion ratio, suggesting a unique evolutionary pattern for the repeats. Finally, we discuss the usefulness of the control region in population genetic and evolutionary studies.

Genetic differentiation in the African malaria vector, Anopheles gambiae s.s., and the problem of taxonomic status

Of the seven recognized species of the Anopheles gambiae complex, A. gambiae s.s. is the most widespread and most important vector of malaria. It is becoming clear that, in parts of West Africa, this nominal species is not a single panmictic unit. We found that the internal transcribed spacer (ITS) of the X-linked rDNA has two distinct sequences with three fixed nucleotide differences; we detected no heterozygotes at these three sites, even in areas of sympatry of the two ITS types. The intergenic spacer (IGS) of this region also displays two distinct sequences that are in almost complete linkage disequilibrium with the distinct ITS alleles. We have designated these two types as S/type I and M/type II. These rDNA types correspond at least partly to the previously recognized chromosomal forms. Here we expand the geographic range of sampling to 251 individuals from 38 populations. Outside of West Africa, a single rDNA type, S/type I, corresponds to the Savanna chromosomal form. In West Africa, both types are often found in a single local sample. To understand if these findings might be due to unusual behavior of the rDNA region, we sequenced the same region for 46 A. arabiensis, a sympatric sibling species. No such distinct discontinuity was observed for this species. Autosomal inversions in one chromosome arm (2R), an insecticide resistance gene on 2L, and this single X-linked region indicate at least two genetically differentiated subpopulations of A. gambiae. Yet, rather extensive studies of other regions of the genome have failed to reveal genetic discontinuity. Evidently, incomplete genetic isolation exists within this single nominal species.

Molecular systematics of Anopheles: from subgenera to subpopulations

The century-old discovery of the role of Anopheles in human malaria transmission precipitated intense study of this genus at the alpha taxonomy level, but until recently little attention was focused on the systematics of this group. The application of molecular approaches to systematic problems ranging from subgeneric relationships to relationships at and below the species level is helping to address questions such as anopheline phylogenetics and biogeography, the nature of species boundaries, and the forces that have structured genetic variation within species. Current knowledge in these areas is reviewed, with an emphasis on the Anopheles gambiae model. The recent publication of the genome of this anopheline mosquito will have a profound impact on inquiries at all taxonomic levels, supplying better tools for estimating phylogeny and population structure in the short term, and ultimately allowing the identification of genes and/or regulatory networks underlying ecological differentiation, speciation, and vectorial capacity.

Gene flow, dispersal, and nested clade analysis among populations of the stonefly Peltoperla tarteri in the southern Appalachians

We examined the genetic structure and phylogeography of populations of the stonefly Peltoperla tarteri in the Southern Appalachians to determine the extent and likely mechanism for dispersal of this stream insect. A 454-base-pair (bp) portion of the mitochondrial control region was sequenced from a minimum of 20 individuals from eight populations. Pairwise FST and exact tests showed high levels of differentiation among almost all populations except those on the same stream. amova analysis detected significant genetic differentiation between streams within drainages (phi(SD) = 0.14, P < 0.001), and there was a slight positive correlation between aquatic distance and genetic distance (r = 0.295, P = 0.03). According to nested clade analysis, the present day pattern of genetic variation in P. tarteri is the result of a historical range expansion coupled with restricted gene flow with isolation by distance. Together, these analyses suggest that adult dispersal is limited and that movement by larvae is the primary dispersal mechanism for P. tarteri.

Attempts to molecularly distinguish cryptic taxa in Anopheles gambiae s.s

Analyses of inversions in polytene chromosomes indicate that, in West Africa, Anopheles gambiae (sensu stricto) may be a complex of more than a single taxonomic unit, and these units have been called chromosomal forms. In order to determine whether this genetic discontinuity extends to the rest of the genome, as would be expected if reproductive isolation exists, we have sequenced several regions of both the nuclear and mitochondrial genomes. With one exception, we were unable to identify any nucleotide sites that differentiate the chromosomal forms. The exception was the internal transcribed spacer (ITS) of the ribosomal DNA (rDNA). Three sites in this region distinguish Mopti chromosomal form from Savanna and Bamako in Mali and Burkina Faso. However, outside these two countries, the association between chromosomal form and rDNA type does not always hold. Together with the variants in the rDNA intergenic spacer (IGS) described in the accompanying papers (della Torre et al., 2001; Favia et al., 2001), we can recognize two major types of rDNA, Type I and Type II (corresponding to molecular forms S and M in della Torre et al., 2001). Type I is widespread in West Africa and is the only type found outside of West Africa (i.e. Tanzania and Madagascar). Type II is confined to West Africa. We were unable to detect any heterozygosity for the ITS types even in five collections containing both types. A sample from the island of São Tomé could not be classified into either Type I or Type II as the rDNA had characteristics of both. In general, our results confirm that An. gambiae is not a single pan-mictic unit, but exactly how to define any new taxa remains problematic. Finally, we have found minor variants of the major rDNA types fixed in local populations; contrary to most previous studies, this suggests restricted gene flow among populations of this species.

Molecular systematics of European Hyalodaphnia: the role of contemporary hybridization in ancient species

We examined phylogenetic relationships among Daphnia using mitochondrial DNA (mtDNA) sequences from the small subunit ribosomal RNA (12S), cytochrome c oxidase subunit I and nuclear DNA sequences from the first and second internal transcribed spacer representing 1612 base positions. Phylogenetic analyses using several species of the three main Daphnia subgenera, Ctenodaphnia, Hyalodaphnia and Daphnia, revealed that the Hyalodaphnia are a monophyletic sister group of the Daphnia. Most Hyalodaphnia species occur on one continent, whereas only three are found in North America and Europe. Endemicity of species is associated with variation in thermal tolerance and habitat differentiation. Although many species of the Hyalodaphnia are known to hybridize in nature, mtDNA divergence is relatively high ca. 9%) compared to other hybridizing arthropods (ca. 3%). Reproductive isolation in Daphnia seems to evolve significantly slower than genetic isolation. We related these findings to what is known about the ecology and genetics of Daphnia in order to better understand the evolutionary diversification of lineages. The relationship of these data to phylogenetic patterns is discussed in the context of speciation processes in Daphnia.

Evidence for mitochondrial introgression between Anopheles bwambae and Anopheles gambiae

There is evidence for introgression between the malaria vectors Anopheles gambiae and Anopheles arabiensis, which belong to the Anopheles gambiae complex. Here we consider their relationship to another member of the complex, Anopheles bwambae. Species identifications were made using rDNA-PCR and revealed one An. bwambae/gambiae hybrid among 459 wild-caught mosquitoes. A 237-nucleotide region of the mitochondrial ND5 gene was sequenced in sixty-nine individuals. The results suggest that there is an An. bwambae specific group of haplotypes. However, two An. bwambae individuals carried haplotypes typical of An. gambiae and An. arabiensis. This is the first evidence for introgression between An. bwambae and other species of the An. gambiae complex and suggests that introgression may be more widespread, also occurring between other members of the complex.

Evolution and structural organisation of mitochondrial DNA control region of myiasis-causing flies

This study reports the molecular characterization of the mtDNA control region (called the A+T-rich region in insects) of five dipteran species which cause myiasis: Cochliomyia hominivorax Coquerel, Cochliomyia macellaria Fabricius, Chrysomya megacephala Fabricius, Lucilia eximia Wiedemann (Diptera: Calliphoridae) and Dermatobia hominis Linnaeus Jr (Diptera: Oestridae). The control region in these species varies in length from 1000 to 1600 bp. Two structural domains with specific evolutionary patterns were identified. These were (1) conserved sequence blocks containing primary sequence motifs, including dinucleotide pyrimidine-purine series and long T-stretches, located at the 5' end adjacent to the tRNA(Ile) gene and (2) a hypervariable domain at the 3' end characterized by increased nucleotide divergence and size variation. A high frequency of A<-->T transversions at nucleotide substitution level indicated directional mutation pressure. The phylogenetic usefulness of the insect control region is discussed.

The current state of insect molecular systematics: a thriving Tower of Babel

Insect molecular systematics has undergone remarkable recent growth. Advances in methods of data generation and analysis have led to the accumulation of large amounts of DNA sequence data from most major insect groups. In addition to reviewing theoretical and methodological advances, we have compiled information on the taxa and regions sequenced from all available phylogenetic studies of insects. It is evident that investigators have not usually coordinated their efforts. The genes and regions that have been sequenced differ substantially among studies and the whole of our efforts is thus little greater than the sum of its parts. The cytochrome oxidase I, 16S, 18S, and elongation factor-1 alpha genes have been widely used and are informative across a broad range of divergences in insects. We advocate their use as standards for insect phylogenetics. Insect molecular systematics has complemented and enhanced the value of morphological and ecological data, making substantial contributions to evolutionary biology in the process. A more coordinated approach focused on gathering homologous sequence data will greatly facilitate such efforts.

Multiple origins of cytologically identical chromosome inversions in the Anopheles gambiae complex

For more than 60 years, evolutionary cytogeneticists have been using naturally occurring chromosomal inversions to infer phylogenetic histories, especially in insects with polytene chromosomes. The validity of this method is predicated on the assumption that inversions arise only once in the history of a lineage, so that sharing a particular inversion implies shared common ancestry. This assumption of monophyly has been generally validated by independent data. We present the first clear evidence that naturally occurring inversions, identical at the level of light microscopic examination of polytene chromosomes, may not always be monophyletic. The evidence comes from DNA sequence analyses of regions within or very near the breakpoints of an inversion called the 2La that is found in the Anopheles gambiae complex. Two species, A. merus and A. arabiensis, which are fixed for the "same" inversion, do not cluster with each other in a phylogenetic analysis of the DNA sequences within the 2La. Rather, A. merus 2La is most closely related to strains of A. gambiae homozygous for the 2L+. A. gambiae and A. merus are sister taxa, the immediate ancestor was evidently homozygous 2L+, and A. merus became fixed for an inversion cytologically identical to that in A. arabiensis. A. gambiae is polymorphic for 2La/2L+, and the 2La in this species is nearly identical at the DNA level to that in A. arabiensis, consistent with the growing evidence that introgression has or is occurring between these two most important vectors of malaria in the world. The parallel evolution of the "same" inversion may be promoted by the presence of selectively important genes within the breakpoints.

Contrasting levels of variability between cytoplasmic genomes and incompatibility types in the mosquito Culex pipiens

Reproductive incompatibilities called cytoplasmic incompatibilities are known to affect a large number of arthropod species and are mediated by Wolbachia, a maternally transmitted microorganism. The crossing relationships between strains of potential hosts define their incompatibility types and it is generally assumed that differences between strains of Wolbachia induce different crossing types. Among all the described host species, the mosquito, Culex pipiens, displays the greatest variability of cytoplasmic incompatibility crossing types. We analysed mitochondrial and bacterial DNA variability in Culex pipiens in order to investigate some possible causes of incompatibility crossing type variability. We sequenced fragments of the ftsZ gene, and the A + T-rich control region of the mtDNA. We also sequenced the second subunit of the mitochondrial cytochrome oxidase (COII) gene, in Culex pipiens and a closely related species, C. torrentium, in order to verify the usefulness of the A + T-rich region for the present purposes. No variability was found in the Wolbachia ftsZ gene fragment, and very limited variation of the mitochondrial marker whatever the compatibility type or the origin of the host. A low variability was found in the A + T-rich region and comparison of divergence of the A + T-rich region and COII gene between C. pipiens and C. torrentium did not reveal any special constraints affecting this region. In contrast to observations in other host species, variability of incompatibility crossing types is not due to multiple infections by distantly related Wolbachia strains.