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

Low coverage of species constrains the use of DNA barcoding to assess mosquito biodiversity

Mosquitoes (Culicidae) represent the main vector insects globally, and they also inhabit many of the terrestrial and aquatic habitats of the world. DNA barcoding and metabarcoding are now widely used in both research and routine practices involving mosquitoes. However, these methodologies rely on information available in databases consisting of barcode sequences representing taxonomically identified voucher specimens. In this study, we assess the availability of public data for mosquitoes in the main online databases, focusing specifically on the two most widely used DNA barcoding markers in Culicidae: COI and ITS2. In addition, we test hypotheses on possible factors affecting species coverage (i.e., the percentage of species covered in the online databases) for COI in different countries and the occurrence of the DNA barcode gap for COI. Our findings showed differences in the data publicly available in the repositories, with a taxonomic or species coverage of 28.4-30.11% for COI in BOLD + GenBank, and 12.32% for ITS2 in GenBank. Afrotropical, Australian and Oriental biogeographic regions had the lowest coverages, while Nearctic, Palearctic and Oceanian had the highest. The Neotropical region had an intermediate coverage. In general, countries with a higher diversity of mosquitoes and higher numbers of medically important species had lower coverage. Moreover, countries with a higher number of endemic species tended to have a higher coverage. Although our DNA barcode gap analyses suggested that the species boundaries need to be revised in half of the mosquito species available in the databases, additional data must be gathered to confirm these results and to allow explaining the occurrence of the DNA barcode gap. We hope this study can help guide regional species inventories of mosquitoes and the completion of a publicly available reference library of DNA barcodes for all mosquito species.

Systematic Review on Diversity and Distribution of Anopheles Species in Gabon: A Fresh Look at the Potential Malaria Vectors and Perspectives

Gabon is located in the malaria hyper-endemic zone, where data concerning malaria vector distribution remains fragmentary, making it difficult to implement an effective vector control strategy. Thus, it becomes crucial and urgent to undertake entomological surveys that will allow a better mapping of the Anopheles species present in Gabon. In this review, we examined different articles dealing with Anopheles in Gabon from ProQuest, Web of Science, PubMed, and Google scholar databases. After applying the eligibility criteria to 7543 articles collected from four databases, 42 studies were included that covered a 91-year period of study. The review revealed a wide diversity of Anopheles species in Gabon with a heterogeneous distribution. Indeed, our review revealed the presence of 41 Anopheles species, of which the most abundant were members of the Gambiae and Nili complexes and those of the Funestus and Moucheti groups. However, our review also revealed that the major and minor vectors of malaria in Gabon are present in both sylvatic, rural, and urban environments. The observation of human malaria vectors in sylvatic environments raises the question of the role that the sylvatic environment may play in maintaining malaria transmission in rural and urban areas. Ultimately, it appears that knowledge of biodiversity and spatial distribution of Anopheles mosquitoes is fragmentary in Gabon, suggesting that additional studies are necessary to complete and update these entomological data, which are useful for the implementation of vector control strategies.

Sympatric occurrence of deeply diverged mitochondrial DNA lineages in Siberian geometrid moths (Lepidoptera: Geometridae): cryptic speciation, mitochondrial introgression, secondary admixture or effect of Wolbachia?

The divergent sympatric mitochondrial lineages within traditionally recognized species present a challenge regularly faced by taxonomists and evolutionary biologists. We encountered this problem when studying the Siberian geometrid moths, Alcis deversata and Thalera chlorosaria. Within each of these species we found two deeply diverged mitochondrial lineages that demonstrated a level of genetic differentiation exceeding the standard interspecific DNA barcode threshold. Using analyses of nuclear genes, morphology, ecological preferences and Wolbachia endosymbionts, we tested five hypotheses that might explain the mitochondrial pattern observed: cryptic speciation, ancestral polymorphism, interspecific mitochondrial introgression, secondary admixture of allopatrically evolved populations and an effect of intracellular Wolbachia endosymbionts. We demonstrate that in A. deversata and Th. chlorosaria the mitochondrial differences are not correlated with differences in nuclear genes, morphology, ecology and Wolbachia infection status, thus not supporting the hypothesis of cryptic species and an effect of Wolbachia. Mitochondrial introgression can lead to a situation in which one species has both its own mitochondrial lineage and the lineage obtained from another species. We found this situation in the species pair Alcis repandata and Alcis extinctaria. We conclude that the mitochondrial heterogeneity in A. deversata and Th. chlorosaria is most likely to be attributable to the secondary admixture of allopatrically evolved populations.

Laurent B, Zubaidah S, Rozi IE, Kosasih S, Shinta, Sukowati S, Hakim L, Haryanto E, Mangunwardoyo W, Bangs MJ, Lobo NF. Anopheles sundaicus complex and the presence of Anopheles epiroticus in Indonesia

Anopheles sundaicus s.l. is an important malaria vector primarily found in coastal landscapes of western and central Indonesia. The species complex has a wide geographical distribution in South and Southeast Asia and exhibits ecological and behavioural variability over its range. Studies on understanding the distribution of different members in the complex and their bionomics related to malaria transmission might be important guiding more effective vector intervention strategies. Female An. sundaicus s.l. were collected from seven provinces, 12 locations in Indonesia representing Sumatra: North Sumatra, Bangka-Belitung, South Lampung, and Bengkulu; in Java: West Java; and the Lesser Sunda Islands: West Nusa Tenggara and East Nusa Tenggara provinces. Sequencing of ribosomal DNA ITS2 gene fragments and two mitochondrial DNA gene markers, COI and cytb, enabled molecular identification of morphologically indistinguishable members of the complex. Findings allowed inference on the distribution of the An. sundaicus s.l. present in Indonesia and further illustrate the phylogenetic relationships of An. epiroticus within the complex. A total of 370 An. sundaicus s.l specimens were analysed for the ITS2 fragment. The ITS2 sequence alignment revealed two consistent species-specific point mutations, a T>C transition at base 479 and a G>T transversion at base 538 that differentiated five haplotypes: TG, CG, TT, CT, and TY. The TG haplotype matched published An. epiroticus-indicative sequences from Thailand, Vietnam and peninsular Malaysia. The previously described insertion event (base 603) was observed in all identified specimens. Analysis of the COI and cytb genes revealed no consistent nucleotide variations that could definitively distinguish An. epiroticus from other members in the Sundaicus Complex. The findings indicate and support the existence of An. epiroticus in North Sumatra and Bangka-Belitung archipelago. Further studies are recommended to determine the full distributional extent of the Sundaicus complex in Indonesia and investigate the role of these species in malaria transmission.

Sequence-based identification of Anopheles species in eastern Ethiopia

BACKGROUND

The recent finding of a typically non-African Anopheles species in eastern Ethiopia emphasizes the need for detailed species identification and characterization for effective malaria vector surveillance. Molecular approaches increase the accuracy and interoperability of vector surveillance data. To develop effective molecular assays for Anopheles identification, it is important to evaluate different genetic loci for the ability to characterize species and population level variation. Here the utility of the internal transcribed spacer 2 (ITS2) and cytochrome oxidase I (COI) loci for detection of Anopheles species from understudied regions of eastern Ethiopia was investigated.

METHODS

Adult mosquitoes were collected from the Harewe locality (east) and Meki (east central) Ethiopia. PCR and Sanger sequencing were performed for portions of the ITS2 and COI loci. Both NCBI's Basic Local Alignment Search tool (BLAST) and phylogenetic analysis using a maximum-likelihood approach were performed to identify species of Anopheles specimens.

RESULTS

Two species from the east Ethiopian collection, Anopheles arabiensis and Anopheles pretoriensis were identified. Analyses of ITS2 locus resulted in delineation of both species. In contrast, analysis of COI locus could not be used to delineate An. arabiensis from other taxa in Anopheles gambiae complex, but could distinguish An. pretoriensis sequences from sister taxa.

CONCLUSION

The lack of clarity from COI sequence analysis highlights potential challenges of species identification within species complexes. These results provide supporting data for the development of molecular assays for delineation of Anopheles in east Ethiopia.

Interspecific gamete compatibility and hybrid larval fitness in reef-building corals: Implications for coral reef restoration

Climate warming is a major cause of the global decline of coral reefs. Active reef restoration, although still in its infancy, is one of several possible ways to help restore coral cover and reef ecosystem function. The deployment of mature coral larvae onto depauperate reef substratum has been shown to significantly increase larval recruitment, providing a novel option for the delivery of ex situ bred coral stock to the reef for restoration purposes. The success of such reef restoration approaches may be improved by the use of coral larval stock augmented for climate resilience. Here we explore whether coral climate resilience can be enhanced via interspecific hybridization through hybrid vigour. Firstly, we assessed cross-fertility of four pairs of Acropora species from the Great Barrier Reef. Temporal isolation in gamete release between the Acropora species was limited, but gametic incompatibility was present with varying strength between species pairs and depending on the direction of the hybrid crosses. We subsequently examined the fitness of hybrid and purebred larvae under heat stress by comparing their survival and settlement success throughout 10 days of exposure to 28 °C, 29.5 °C and 31 °C. Fitness of the majority of Acropora hybrid larvae was similar to that of the purebred larvae of both parental species, and in some instances it was higher than that of the purebred larvae of one of the parental species. Lower hybrid fertilization success did not affect larval fitness. These findings indicate that high hybrid fitness can be achieved after overcoming partial prezygotic barriers, and that interspecific hybridization may be a tool to enhance coral recruitment and climate resilience.

Spatio-temporal genetic structure of Anopheles gambiae in the Northwestern Lake Victoria Basin, Uganda: implications for genetic control trials in malaria endemic regions

Background

Understanding population genetic structure in the malaria vector Anopheles gambiae (s.s.) is crucial to inform genetic control and manage insecticide resistance. Unfortunately, species characteristics such as high nucleotide diversity, large effective population size, recent range expansion, and high dispersal ability complicate the inference of genetic structure across its range in sub-Saharan Africa. The ocean, along with the Great Rift Valley, is one of the few recognized barriers to gene flow in this species, but the effect of inland lakes, which could be useful sites for initial testing of genetic control strategies, is relatively understudied. Here we examine Lake Victoria as a barrier between the Ugandan mainland and the Ssese Islands, which lie up to 60 km offshore. We use mitochondrial DNA (mtDNA) from populations sampled in 2002, 2012 and 2015, and perform Bayesian cluster analysis on mtDNA combined with microsatellite data previously generated from the same 2002 mosquito DNA samples.

Results

Hierarchical analysis of molecular variance and Bayesian clustering support significant differentiation between the mainland and lacustrine islands. In an mtDNA haplotype network constructed from this and previous data, haplotypes are shared even between localities separated by the Rift Valley, a result that more likely reflects retention of shared ancestral polymorphism than contemporary gene flow.

Conclusions

The relative genetic isolation of An. gambiae on the Ssese Islands, their small size, level terrain and ease of access from the mainland, the relative simplicity of the vectorial system, and the prevalence of malaria, are all attributes that recommend these islands as possible sites for the testing of genetic control strategies.

DNA barcoding mosquitoes: advice for potential prospectors

Mosquitoes’ importance as vectors of pathogens that drive disease underscores the importance of precise and comparable methods of taxa identification among their species. While several molecular targets have been used to study mosquitoes since the initiation of PCR in the 1980s, its application to mosquito identification took off in the early 1990s. This review follows the research's recent journey into the use of mitochondrial DNA (mtDNA) cytochrome oxidase 1 (COI or COX1) as a DNA barcode target for mosquito species identification – a target whose utility for discriminating mosquitoes is now escalating. The pros and cons of using a mitochondrial genome target are discussed with a broad sweep of the mosquito literature suggesting that nuclear introgressions of mtDNA sequences appear to be uncommon and that the COI works well for distantly related taxa and shows encouraging utility in discriminating more closely related species such as cryptic/sibling species groups. However, the utility of COI in discriminating some closely related groups can be problematic and investigators are advised to proceed with caution as problems with incomplete lineage sorting and introgression events can result in indistinguishable COI sequences appearing in reproductively independent populations. In these – if not all – cases, it is advisable to run a nuclear marker alongside the mtDNA and thus the utility of the ribosomal DNA – and in particular the internal transcribed spacer 2 – is also briefly discussed as a useful counterpoint to the COI.

Introgression among three rockfish species (Sebastes spp.) in the Salish Sea, northeast Pacific Ocean

Interspecific hybridization is often seen as a major conservation issue, potentially threatening endangered species and decreasing biodiversity. In natural populations, the conservation implications of hybridization depends on both on anthropogenic factors and the evolutionary processes maintaining the hybrid zone. However, the timeline and patterns of hybridization in the hybrid zone are often not known. Therefore, species conservation becomes a concern when recent anthropogenic changes influence hybridization and not if hybridization is part of a long-term process. Here, we use sequence data from one mitochondrial gene, three nuclear introns and one nuclear exon to estimate the direction, geographic extent, frequency and possible timeline of hybridization between three rockfish species (Sebastes auriculatus, S. caurinus, S. maliger) in the Salish Sea, Washington, USA. We show that (i) introgression occurred much more frequently in the Salish Sea than on the outer coast, (ii) introgression was highly asymmetrical from S. maliger into the other two species, (iii) almost 40% of individuals in the Salish Sea were hybrids, with frequency of hybrids increasing with isolation from the coast, and (iv) all hybrids were later generation backcrosses rather than F1 hybrids. Our results suggest long-standing low-level hybridization rather than recent onset of interbreeding because of human induced environmental change, possibly facilitated by specific environmental conditions in the sub-basins of the Salish Sea, and by differences in population sizes during recolonization of the area after the last glaciation. This rockfish hybrid system, with asymmetrical introgression and the maintenance of parental species, may prove useful to study both mechanisms that maintain species boundaries and that facilitate speciation in the presence of rapid environmental change.

Mitochondrial genomes of Anopheles arabiensis, An. gambiae and An. coluzzii show no clear species division

Here we report the complete mitochondrial sequences of 70 individual field collected mosquito specimens from throughout Sub-Saharan Africa. We generated this dataset to identify species specific markers for the following Anopheles species and chromosomal forms: An. arabiensis, An. coluzzii (The Forest and Mopti chromosomal forms) and An. gambiae (The Bamako and Savannah chromosomal forms).  The raw Illumina sequencing reads were mapped to the NC_002084 reference mitogenome sequence. A total of 783 single nucleotide polymorphisms (SNPs) were detected on the mitochondrial genome, of which 460 are singletons (58.7%). None of these SNPs are suitable as molecular markers to distinguish among An. arabiensis, An. coluzzii and An. gambiae or any of the chromosomal forms. The lack of species or chromosomal form specific markers is also reflected in the constructed phylogenetic tree, which shows no clear division among the operational taxonomic units considered here.

Phylogeny with introgression in Habronattus jumping spiders (Araneae: Salticidae)

Background

Habronattus is a diverse clade of jumping spiders with complex courtship displays and repeated evolution of Y chromosomes. A well-resolved species phylogeny would provide an important framework to study these traits, but has not yet been achieved, in part because the few genes available in past studies gave conflicting signals. Such discordant gene trees could be the result of incomplete lineage sorting (ILS) in recently diverged parts of the phylogeny, but there are indications that introgression could be a source of conflict.

Results

To infer Habronattus phylogeny and investigate the cause of gene tree discordance, we assembled transcriptomes for 34 Habronattus species and 2 outgroups. The concatenated 2.41 Mb of nuclear data (1877 loci) resolved phylogeny by Maximum Likelihood (ML) with high bootstrap support (95-100%) at most nodes, with some uncertainty surrounding the relationships of H. icenoglei, H. cambridgei, H. oregonensis, and Pellenes canadensis. Species tree analyses by ASTRAL and SVDQuartets gave almost completely congruent results. Several nodes in the ML phylogeny from 12.33 kb of mitochondrial data are incongruent with the nuclear phylogeny and indicate possible mitochondrial introgression: the internal relationships of the americanus and the coecatus groups, the relationship between the altanus, decorus, banksi, and americanus group, and between H. clypeatus and the coecatus group. To determine the relative contributions of ILS and introgression, we analyzed gene tree discordance for nuclear loci longer than 1 kb using Bayesian Concordance Analysis (BCA) for the americanus group (679 loci) and the VCCR clade (viridipes/clypeatus/coecatus/roberti groups) (517 loci) and found signals of introgression in both. Finally, we tested specifically for introgression in the concatenated nuclear matrix with Patterson’s D statistics and D_FOIL. We found nuclear introgression resulting in substantial admixture between americanus group species, between H. roberti and the clypeatus group, and between the clypeatus and coecatus groups.

Conclusions

Our results indicate that the phylogenetic history of Habronattus is predominantly a diverging tree, but that hybridization may have been common between phylogenetically distant species, especially in subgroups with complex courtship displays.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1137-x) contains supplementary material, which is available to authorized users.

Extensive shared polymorphism at non-MHC immune genes in recently diverged North American prairie grouse

Gene polymorphisms shared between recently diverged species are thought to be widespread and most commonly reflect introgression from hybridization or retention of ancestral polymorphism through incomplete lineage sorting. Shared genetic diversity resulting from incomplete lineage sorting is usually maintained for a relatively short period of time, but under strong balancing selection it may persist for millions of years beyond species divergence (balanced trans-species polymorphism), as in the case of the major histocompatibility complex (MHC) genes. However, balancing selection is much less likely to act on non-MHC immune genes. The aim of this study was to investigate the patterns of shared polymorphism and selection at non-MHC immune genes in five grouse species from Centrocercus and Tympanuchus genera. For this purpose, we genotyped five non-MHC immune genes that do not interact directly with pathogens, but are involved in signaling and regulate immune cell growth. In contrast to previous studies with MHC, we found no evidence for balancing selection or balanced trans-species polymorphism among the non-MHC immune genes. No haplotypes were shared between genera and in most cases more similar allelic variants sorted by genus. Between species within genera, however, we found extensive shared polymorphism, which was most likely attributable to introgression or incomplete lineage sorting following recent divergence and large ancestral effective population size (i.e., weak genetic drift). Our study suggests that North American prairie grouse may have attained relatively low degree of reciprocal monophyly at nuclear loci and reinforces the rarity of balancing selection in non-MHC immune genes.

An evaluation of the suitability of COI and COII gene variation for reconstructing the phylogeny of, and identifying cryptic species in, anopheline mosquitoes (Diptera Culicidae)

We assessed the practicality and effectiveness of using variation in the mitochondrial COI and COII genes to discriminate species and reconstruct the phylogeny of anophelene mosquitoes. Phylogenetic relationships among the subfamily Anophelinae were inferred from portions of the mitochondrial COI (92 species) and COII genes (108 species). Phylogenetic trees were reconstructed on the basis of parsimony, maximum likelihood and Bayesian methods. The suitability of COI and COII gene variation for identifying cryptic species was compared by comparing the sequence divergence within species groups and complexes. The results show that the COI gene was more useful for identifying sibling and cryptic species, but that phylogenetic relationships reconstructed using the COII gene were more similar to those based on morphological data. We conclude that: (1) there is a significant molecular divergence among An. sinensis; (2) the COI and COII are valid genetic markers for resolving taxonomic relationships among anopheline mosquitoes and the resultant phylogeny raises some questions about the taxonomic status of anopheline species groups and complexes; (3) the genus Anopheles is not demonstrably monophyletic with regard to the genus Bironella; (4) the subgenera Kerteszia and Nyssorhynchus are monophyletic; (5) below the group-level, COI data support the existence of monophyletic taxa within the Anopheles funestus, Anopheles maculipennis and Anopheles strode and Anopheles barbirostris subgroups, and within the Anopheles nuneztovari complex, whereas COII data support the monophyletic taxa within the Anopheles minimus and Anopheles oswaldoi subgroups, and Anopheles hyrcanus group. The monophyletic taxa within the Anopheles gambiae and Anopheles albitarsis complexes are supported by both COI and COII data.

DNA barcoding of crickets, katydids and grasshoppers (Orthoptera) from Central Europe with focus on Austria, Germany and Switzerland

We present a DNA barcoding study on the insect order Orthoptera that was generated in collaboration between four barcoding projects in three countries, viz. Barcoding Fauna Bavarica (Germany), German Barcode of Life, Austrian Barcode of Life and Swiss Barcode of Life. Our data set includes 748 COI sequences from 127 of the 162 taxa (78.4%) recorded in the three countries involved. Ninety-three of these 122 species (76.2%, including all Ensifera) can be reliably identified using DNA barcodes. The remaining 26 caeliferan species (families Acrididae and Tetrigidae) form ten clusters that share barcodes among up to five species, in three cases even across different genera, and in six cases even sharing individual barcodes. We discuss incomplete lineage sorting and hybridization as most likely causes of this phenomenon, as the species concerned are phylogenetically young and hybridization has been previously observed. We also highlight the problem of nuclear mitochondrial pseudogenes (numts), a known problem in the barcoding of orthopteran species, and the possibility of Wolbachia infections. Finally, we discuss the possible taxonomic implications of our barcoding results and point out future research directions.

Molecular taxonomy provides new insights into anopheles species of the neotropical arribalzagia series

Phylogenetic analysis of partial mitochondrial cytochrome oxidase c subunit I (COI) and nuclear internal transcribed spacer 2 (ITS2) sequences were used to evaluate initial identification and to investigate phylogenetic relationships of seven Anopheles morphospecies of the Arribalzagia Series from Colombia. Phylogenetic trees recovered highly supported clades for An. punctimaculas.s., An. calderoni, An. malefactor s.l., An. neomaculipalpus, An. apicimacula s.l., An. mattogrossensis and An. peryassui. This study provides the first molecular confirmation of An. malefactorfrom Colombia and discovered conflicting patterns of divergence for the molecular markers among specimens from northeast and northern Colombia suggesting the presence of two previously unrecognized Molecular Operational Taxonomic Units (MOTUs). Furthermore, two highly differentiated An. apicimacula MOTUs previously found in Panama were detected. Overall, the combined molecular dataset facilitated the detection of known and new Colombian evolutionary lineages, and constitutes the baseline for future research on their bionomics, ecology and potential role as malaria vectors.

High morphological and genetic variabilities of Ochlerotatus scapularis, a potential vector of filarias and arboviruses

Background

Ochlerotatus scapularis is a potential vector of filarias and arboviruses in the Neotropics. This species was once typically associated with sylvatic environments; however, cases of synanthropy and urbanization of this species have been increasingly reported in southeast Brazil. Despite the medical relevance of Oc. scapularis, its populational variability is not yet known. To our knowledge, this is the first report describing the morphological and genetic variabilities of this species.

Methods

Population samples were characterized using the cytochrome oxidase subunit I (COI) mitochondrial gene and wing geometrics. Adult mosquitoes were collected from five sampling sites from remnants of the Atlantic forest embedded in the urban or rural areas of southeast Brazil.

Results

In the 130 individuals analyzed, 46 COI haplotypes were detected. Haplotype diversity was high and ranged from 0.66 to 0.97. Six haplotypes were present in 61% of the individuals, whereas the remaining haplotypes were less frequent (39%). Wing shape was also highly polymorphic. Differentiation of populations across sampling sites according to genetic distances (F_st = −0.009 to 0.060) and morphological distances (Q_st = 0.47) indicated that populations were not identical. No correlations were noted for phenetic and genetic diversities (p = 0.19) or for genetic or phenetic distances with geographical distances (p = 0.2 and p = 0.18, respectively).

Conclusions

Our study results suggest that Oc. scapularis has a rich genetic patrimony, even though its habitat is fragmented. Implications of such genetic richness with respect to vectorial competence, plasticity, and ability to exploit urbanized areas need to be further investigated.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-0740-6) contains supplementary material, which is available to authorized users.

Adaptive introgression in an African malaria mosquito coincident with the increased usage of insecticide-treated bed nets

Animal species adapt to changes in their environment, including man-made changes such as the introduction of insecticides, through selection for advantageous genes already present in populations or newly arisen through mutation. A possible alternative mechanism is the acquisition of adaptive genes from related species via a process known as adaptive introgression. Differing levels of insecticide resistance between two African malaria vectors, Anopheles coluzzii and Anopheles gambiae, have been attributed to assortative mating between the two species. In a previous study, we reported two bouts of hybridization observed in the town of Selinkenyi, Mali in 2002 and 2006. These hybridization events did not appear to be directly associated with insecticide-resistance genes. We demonstrate that during a brief breakdown in assortative mating in 2006, A. coluzzii inherited the entire A. gambiae-associated 2L divergence island, which includes a suite of insecticide-resistance alleles. In this case, introgression was coincident with the start of a major insecticide-treated bed net distribution campaign in Mali. This suggests that insecticide exposure altered the fitness landscape, favoring the survival of A. coluzzii/A. gambiae hybrids, and provided selection pressure that swept the 2L divergence island through A. coluzzii populations in Mali. We propose that the work described herein presents a unique description of the temporal dynamics of adaptive introgression in an animal species and represents a mechanism for the rapid evolution of insecticide resistance in this important vector of human malaria in Africa.

Mosquito genomics. Extensive introgression in a malaria vector species complex revealed by phylogenomics

Introgressive hybridization is now recognized as a widespread phenomenon, but its role in evolution remains contested. Here, we use newly available reference genome assemblies to investigate phylogenetic relationships and introgression in a medically important group of Afrotropical mosquito sibling species. We have identified the correct species branching order to resolve a contentious phylogeny and show that lineages leading to the principal vectors of human malaria were among the first to split. Pervasive autosomal introgression between these malaria vectors means that only a small fraction of the genome, mainly on the X chromosome, has not crossed species boundaries. Our results suggest that traits enhancing vectorial capacity may be gained through interspecific gene flow, including between nonsister species.

A genotypically distinct, melanic variant of Anopheles arabiensis in Sudan is associated with arid environments

BACKGROUND

Anopheles arabiensis, an important malaria vector in Sudan and other countries in sub-Saharan Africa, exhibits considerable ecological and behavioural plasticity allowing it to survive in the harsh conditions of arid regions. It has been shown that adult populations of An. arabiensis in the semi-desert habitat of western Khartoum State survive through the long dry season in a state of partial aestivation, characterized by limited feeding activity and a degree of arrested ovarian development. Anopheles arabiensis in these sites occurs in two phenotypic forms. One is large and heavily melanized, the other has the typical characteristics of An. arabiensis as found elsewhere in Africa. The extent of genetic variation in these forms was examined in widely separated locations in Sudan, including Kassala, Gedaref and the Northern States between 1998 and 1999 and 2004 and 2006.

METHODS

Each mosquito specimen was identified using standard morphological keys and a species-specific PCR test. Sequence variation in a 660 bp fragment of the mtDNA ND5 coding region was examined and the extent of genetic divergence between the forms was estimated from FST values using DNASP version 4.9. TCS 1.13 software was used to determine the genealogical relationships and to reflect clustering among mtDNA haplotypes.

RESULTS

The melanic and normal forms were found in sympatry in Kassala, Gedaref and Khartoum states, with the melanic form commonest in the hottest and most arid areas. Both forms were encountered in the periods of study: 1998-1999, and 2004-2006. Only ten specimens of An. arabiensis were collected from the Northern State in February 2006, all of which were of the normal form.Based on the ND5 analysis, there was a marked subdivision between the normal and melanic forms (FST = 0.59). Furthermore, the melanic form showed more genetic variability, as measured by haplotype diversity (0.95) compared with the normal form (0.57), suggesting larger effective population.

CONCLUSIONS

This is the first demonstration of correspondent phenotypic and genetic structuring in An. arabiensis. The high level of genetic differentiation shown by the mtDNA ND5 locus suggests that the two forms may represent separate species. It is hypothesized that the melanic form is better adapted to hot and arid environments.

Islands and stepping-stones: comparative population structure of Anopheles gambiae sensu stricto and Anopheles arabiensis in Tanzania and implications for the spread of insecticide resistance

Population genetic structures of the two major malaria vectors Anopheles gambiae s.s. and An. arabiensis, differ markedly across Sub-Saharan Africa, which could reflect differences in historical demographies or in contemporary gene flow. Elucidation of the degree and cause of population structure is important for predicting the spread of genetic traits such as insecticide resistance genes or artificially engineered genes. Here the population genetics of An. gambiae s.s. and An. arabiensis in the central, eastern and island regions of Tanzania were compared. Microsatellite markers were screened in 33 collections of female An. gambiae s.l., originating from 22 geographical locations, four of which were sampled in two or three years between 2008 and 2010. An. gambiae were sampled from six sites, An. arabiensis from 14 sites, and both species from two sites, with an additional colonised insectary sample of each species. Frequencies of the knock-down resistance (kdr) alleles 1014S and 1014F were also determined. An. gambiae exhibited relatively high genetic differentiation (average pairwise FST = 0.131), significant even between nearby samples, but without clear geographical patterning. In contrast, An. arabiensis exhibited limited differentiation (average FST = 0.015), but strong isolation-by-distance (Mantel test r = 0.46, p = 0.0008). Most time-series samples of An. arabiensis were homogeneous, suggesting general temporal stability of the genetic structure. An. gambiae populations from Dar es Salaam and Bagamoyo were found to have high frequencies of kdr 1014S (around 70%), with almost 50% homozygote but was at much lower frequency on Unguja Island, with no. An. gambiae population genetic differentiation was consistent with an island model of genetic structuring with highly restricted gene flow, contrary to An. arabiensis which was consistent with a stepping-stone model of extensive, but geographically-restricted gene flow.

Contemporary gene flow between wild An. gambiae s.s. and An. arabiensis

Background

In areas where the morphologically indistinguishable malaria mosquitoes Anopheles gambiae Giles and An. arabiensis Patton are sympatric, hybrids are detected occasionally via species-diagnostic molecular assays. An. gambiae and An. arabiensis exhibit both pre- and post-reproductive mating barriers, with swarms largely species-specific and male F1 (first-generation) hybrids sterile. Consequently advanced-stage hybrids (back-crosses to parental species), which would represent a route for potentially-adaptive introgression, are expected to be very rare in natural populations. Yet the use of one or two physically linked single-locus diagnostic assays renders them indistinguishable from F1 hybrids and levels of interspecific gene flow are unknown.

Methods

We used data from over 350 polymorphic autosomal SNPs to investigate post F1 gene flow via patterns of genomic admixture between An. gambiae and An. arabiensis from eastern Uganda. Simulations were used to investigate the statistical power to detect hybrids with different levels of crossing and to identify the hybrid category significantly admixed genotypes could represent.

Results

A range of admixture proportions were detected for 11 field-collected hybrids identified via single-locus species-diagnostic PCRs. Comparison of admixture data with simulations indicated that at least seven of these hybrids were advanced generation crosses, with backcrosses to each species identified. In addition, of 36 individuals typing as An. gambiae or An. arabiensis that exhibited outlying admixture proportions, ten were identified as significantly mixed backcrosses, and at least four of these were second or third generation crosses.

Conclusions

Our results show that hybrids detected using standard diagnostics will often be hybrid generations beyond F1, and that in our study area around 5% (95% confidence intervals 3%-9%) of apparently ‘pure’ species samples may also be backcrosses. This is likely an underestimate because of rapidly-declining detection power beyond the first two backcross generations. Post-F1 gene flow occurs at a far from inconsequential rate between An. gambiae and An. arabiensis, and, especially for traits under strong selection, could readily lead to adaptive introgression of genetic variants relevant for vector control.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-3305-7-345) contains supplementary material, which is available to authorized users.

Genomic analyses of three malaria vectors reveals extensive shared polymorphism but contrasting population histories

Anopheles gambiae s.l. are important malaria vectors, but little is known about their genomic variation in the wild. Here, we present inter- and intraspecies analysis of genome-wide RADseq data, in three Anopheles gambiae s.l. species collected from East Africa. The mosquitoes fall into three genotypic clusters representing described species (A. gambiae, A. arabiensis, and A. merus) with no evidence of cryptic breeding units. Anopheles merus is the most divergent of the three species, supporting a recent new phylogeny based on chromosomal inversions. Even though the species clusters are well separated, there is extensive shared polymorphism, particularly between A. gambiae and A. arabiensis. Divergence between A. gambiae and A. arabiensis does not vary across the autosomes but is higher in X-linked inversions than elsewhere on X or on the autosomes, consistent with the suggestion that this inversion (or a gene within it) is important in reproductive isolation between the species. The 2La/2L+(a) inversion shows no more evidence of introgression between A. gambiae and A. arabiensis than the rest of the autosomes. Population differentiation within A. gambiae and A. arabiensis is weak over approximately 190-270 km, implying no strong barriers to dispersal. Analysis of Tajima's D and the allele frequency spectrum is consistent with modest population increases in A. arabiensis and A. merus, but a more complex demographic history of expansion followed by contraction in A. gambiae. Although they are less than 200 km apart, the two A. gambiae populations show evidence of different demographic histories.

Multilocus analysis of divergence and introgression in sympatric and allopatric sibling species of the Lutzomyia longipalpis complex in Brazil

BACKGROUND

Lutzomyia longipalpis, the main vector of visceral leishmaniasis in Latin America, is a complex of sibling species. In Brazil, a number of very closely related sibling species have been revealed by the analyses of copulation songs, sex pheromones and molecular markers. However, the level of divergence and gene flow between the sibling species remains unclear. Brazilian populations of this vector can be divided in two main groups: one producing Burst-type songs and the Cembrene-1 pheromone and a second more diverse group producing various Pulse song subtypes and different pheromones.

METHODOLOGY/PRINCIPAL FINDINGS

We analyzed 21 nuclear loci in two pairs of Brazilian populations: two sympatric populations from the Sobral locality (1S and 2S) in northeastern Brazil and two allopatric populations from the Lapinha and Pancas localities in southeastern Brazil. Pancas and Sobral 2S are populations of the Burst/Cembrene-1 species while Lapinha and Sobral 1S are two putative incipient species producing the same pheromone and similar Pulse song subtypes. The multilocus analysis strongly suggests the occurrence of gene flow during the divergence between the sibling species, with different levels of introgression between loci. Moreover, this differential introgression is asymmetrical, with estimated gene flow being higher in the direction of the Burst/Cembrene-1 species.

CONCLUSIONS/SIGNIFICANCE

The results indicate that introgressive hybridization has been a crucial phenomenon in shaping the genome of the L. longipalpis complex. This has possible epidemiological implications and is particularly interesting considering the potential for increased introgression caused by man-made environmental changes and the current trend of leishmaniasis urbanization in Brazil.

Evidence for the occurrence of two sympatric sibling species within the Anopheles (Kerteszia) cruzii complex in southeast Brazil and the detection of asymmetric introgression between them using a multilocus analysis

BACKGROUND

Anopheles (Kerteszia) cruzii (Diptera: Culicidae) is a primary vector of human and simian malaria parasites in southern and southeastern Brazil. Earlier studies using chromosome inversions, isoenzymes and a number of molecular markers have suggested that An. cruzii is a species complex.

RESULTS

In this study, a multilocus approach using six loci, three circadian clock genes and three encoding ribosomal proteins, was carried out to investigate in more detail the genetic differentiation between the An. cruzii populations from Florianópolis-Santa Catarina (southern Brazil) and Itatiaia-Rio de Janeiro States (southeastern Brazil). The analyses were performed first comparing Florianópolis and Itatiaia, and then comparing the two putative sympatric incipient species from Itatiaia (Itatiaia A and Itatiaia B). The analysis revealed high FST values between Florianópolis and Itatiaia (considering Itatiaia A and B together) and also between the sympatric Itatiaia A and Itatiaia B, irrespective of their function. Also, using the IM program, no strong indication of migration was found between Florianópolis and Itatiaia (considering Itatiaia A and B together) using all loci together, but between Itatiaia A and Itatiaia B, the results show evidence of migration only in the direction of Itatiaia B.

CONCLUSIONS

The results of the multilocus analysis indicate that Florianópolis and Itatiaia represent different species of the An. cruzii complex that diverged around 0.6 Mya, and also that the Itatiaia sample is composed of two sympatric incipient species A and B, which diverged around 0.2 Mya. Asymmetric introgression was found between the latter two species despite strong divergence in some loci.

Insecticide resistance monitoring of field-collected Anopheles gambiae s.l. populations from Jinja, eastern Uganda, identifies high levels of pyrethroid resistance

Insecticide resistance in the malaria vector Anopheles gambiae s.l. (Diptera: Culicidae) threatens insecticide-based control efforts, necessitating regular monitoring. We assessed resistance in field-collected An. gambiae s.l. from Jinja, Uganda using World Health Organization (WHO) bioassays. Only An. gambiae s.s. and An. arabiensis (≈70%) were present. Female An. gambiae exhibited extremely high pyrethroid resistance (permethrin LT50 > 2 h; deltamethrin LT50 > 5 h). Female An. arabiensis were resistant to permethrin and exhibited reduced susceptibility to deltamethrin. However, while An. gambiae were DDT resistant, An. arabiensis were fully susceptible. Both species were fully susceptible to bendiocarb and fenitrothion. Kdr 1014S has increased rapidly in the Jinja population of An. gambiae s.s. and now approaches fixation (≈95%), consistent with insecticide-mediated selection, but is currently at a low frequency in An. arabiensis (0.07%). Kdr 1014F was also at a low frequency in An. gambiae. These frequencies preclude adequately-powered tests for an association with phenotypic resistance. PBO synergist bioassays resulted in near complete recovery of pyrethroid susceptibility suggesting involvement of CYP450s in resistance. A small number (0.22%) of An. gambiae s.s. ×An. arabiensis hybrids were found, suggesting the possibility of introgression of resistance alleles between species. The high levels of pyrethroid resistance encountered in Jinja threaten to reduce the efficacy of vector control programmes which rely on pyrethroid-impregnated bednets or indoor spraying of pyrethroids.

An analysis of two island groups as potential sites for trials of transgenic mosquitoes for malaria control

Considerable technological advances have been made towards the generation of genetically modified mosquitoes for vector control. In contrast, less progress has been made towards field evaluations of transformed mosquitoes which are critical for evaluating the success of, and hazards associated with, genetic modification. Oceanic islands have been highlighted as potentially the best locations for such trials. However, population genetic studies are necessary to verify isolation. Here, we used a panel of genetic markers to assess for evidence of genetic isolation of two oceanic island populations of the African malaria vector, Anopheles gambiae s.s. We found no evidence of isolation between the Bijagós archipelago and mainland Guinea-Bissau, despite separation by distances beyond the known dispersal capabilities of this taxon. Conversely, the Comoros Islands appear to be genetically isolated from the East African mainland, and thus represent a location worthy of further investigation for field trials. Based on assessments of gene flow within and between the Comoros islands, the island of Grande Comore was found to be genetically isolated from adjacent islands and also exhibited local population structure, indicating that it may be the most suitable site for trials with existing genetic modification technologies.

Wing geometric morphometrics and molecular assessment of members in the Albitarsis Complex from Colombia

Malaria parasites are transmitted to humans by female mosquitoes of the genus Anopheles. The Albitarsis Complex harbours at least eight species not readily differentiable by morphology. This complicates the determination of those species involved in malaria transmission and the implementation of targeted and effective vector control strategies. In Colombia, there is little information about the identity and distribution of the Albitarsis Complex members. In this work, COI DNA barcoding was used to assign specimens Anopheles albitarsis s.l. to any of the previously designated species of the Albitarsis Complex. Two molecular operational taxonomic units (MOTUs), differentially distributed in Colombia, were detected, A. albitarsis I in the NW and NE, and A. albitarsis F, E and NE Colombia. In contrast, nuclear white gene and ITS2 sequence analyses did not allow differentiating between the MOTUs. Wing landmark-based geometric morphometrics applied to explore intertaxa phenotypic heterogeneity showed a subtle but significant difference in size, while shape did not allow the separation of the MOTUs. In general, the multiple marker analysis was not supportive of the existence in Colombia of more than one species of the Albitarsis Complex.

Population genetic structure of the major malaria vector Anopheles funestus s.s. and allied species in southern Africa

Background

Anopheles funestus s.s., one of the major malaria vectors in sub-Saharan Africa, belongs to a group of eleven African species that are morphologically similar at the adult stage, most of which do not transmit malaria. The population structure of An. funestus based on mitochondrial DNA data led to the description of two cryptic subdivisions, clade I widespread throughout Africa and clade II known only from Mozambique and Madagascar. In this study, we investigated five common members of the Anopheles funestus group in southern Africa in order to determine relationships within and between species.

Methods

A total of 155 specimens of An. funestus, An. parensis, An. vaneedeni, An. funestus-like and An. rivulorum from South Africa, Mozambique and Malawi were used for the study. The population genetic structure was assessed within and between populations using mitochondrial DNA.

Results

The phylogenetic trees revealed three main lineages: 1) An. rivulorum; 2) An. funestus-like clade I and An. parensis clade II; and 3) An. funestus clades I and II, An. funestus-like clade II, An. parensis clade I and An. vaneedeni clades I and II. Within An. funestus, 32 specimens from Mozambique consisted of 40.6% clade I and 59.4% clade II while all 21 individuals from Malawi were clade I. In the analysis of mitochondrial DNA sequences, there were 37 polymorphic sites and 9 fixed different nucleotides for ND5 and 21 polymorphic sites and 6 fixed different nucleotides for COI between the two An. funestus clades. The results for COI supported the ND5 analysis.

Conclusion

This is the first report comparing An. funestus group species including An. funestus clades I and II and the new species An. funestus-like. Anopheles funestus clade I is separated from the rest of the members of the An. funestus subgroup and An. funestus-like is distinctly distributed from the other species in this study. However, there were two clades for An. funestus-like, An. parensis and An. vaneedeni. Further investigations are needed to determine what these results mean in terms of the specific status of the clades within each taxon and whether this has any epidemiological implications for malaria transmission.

Next-generation hybridization and introgression

Hybridization has a major role in evolution-from the introgression of important phenotypic traits between species, to the creation of new species through hybrid speciation. Molecular studies of hybridization aim to understand the class of hybrids and the frequency of introgression, detect the signature of ancient hybridization, and understand the behaviour of introgressed loci in their new genomic background. This often involves a large investment in the design and application of molecular markers, leading to a compromise between the depth and breadth of genomic data. New techniques designed to assay a large sub-section of the genome, in association with next-generation sequencing (NGS) technologies, will allow genome-wide hybridization and introgression studies in organisms with no prior sequence data. These detailed genotypic data will unite the breadth of sampling of loci characteristic of population genetics with the depth of sequence information associated with molecular phylogenetics. In this review, we assess the theoretical and methodological constraints that limit our understanding of natural hybridization, and promote the use of NGS for detecting hybridization and introgression between non-model organisms. We also make recommendations for the ways in which emerging techniques, such as pooled barcoded amplicon sequencing and restriction site-associated DNA tags, should be used to overcome current limitations, and enhance our understanding of this evolutionary significant process.

Weak prezygotic isolating mechanisms in threatened Caribbean Acropora corals

The Caribbean corals, Acropora palmata and A. cervicornis, recently have undergone drastic declines primarily as a result of disease. Previous molecular studies have demonstrated that these species form a hybrid (A. prolifera) that varies in abundance throughout the range of the parental distribution. There is variable unidirectional introgression across loci and sites of A. palmata genes flowing into A. cervicornis. Here we examine the efficacy of prezygotic reproductive isolating mechanisms within these corals including spawning times and choice and no-choice fertilization crosses. We show that these species have subtly different mean but overlapping spawning times, suggesting that temporal isolation is likely not an effective barrier to hybridization. We found species-specific differences in gametic incompatibilities. Acropora palmata eggs were relatively resistant to hybridization, especially when conspecific sperm are available to outcompete heterospecific sperm. Acropora cervicornis eggs demonstrated no evidence for gametic incompatibility and no evidence of reduced viability after aging four hours. This asymmetry in compatibility matches previous genetic data on unidirectional introgression.

Asymmetric introgression between the M and S forms of the malaria vector, Anopheles gambiae, maintains divergence despite extensive hybridization

The suggestion that genetic divergence can arise and/or be maintained in the face of gene flow has been contentious since first proposed. This controversy and a rarity of good examples have limited our understanding of this process. Partially reproductively isolated taxa have been highlighted as offering unique opportunities for identifying the mechanisms underlying divergence with gene flow. The African malaria vector, Anopheles gambiae s.s., is widely regarded as consisting of two sympatric forms, thought by many to represent incipient species, the M and S molecular forms. However, there has been much debate about the extent of reproductive isolation between M and S, with one view positing that divergence may have arisen and is being maintained in the presence of gene flow, and the other proposing a more advanced speciation process with little realized gene flow because of low hybrid fitness. These hypotheses have been difficult to address because hybrids are typically rare (<1%). Here, we assess samples from an area of high hybridization and demonstrate that hybrids are fit and responsible for extensive introgression. Nonetheless, we show that strong divergent selection at a subset of loci combined with highly asymmetric introgression has enabled M and S to remain genetically differentiated despite extensive gene flow. We propose that the extent of reproductive isolation between M and S varies across West Africa resulting in a 'geographic mosaic of reproductive isolation'; a finding which adds further complexity to our understanding of divergence in this taxon and which has considerable implications for transgenic control strategies.

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.

Molecular characterization and evolution of a gene family encoding male-specific reproductive proteins in the African malaria vector Anopheles gambiae

Background

During copulation, the major Afro-tropical malaria vector Anopheles gambiae s.s. transfers male accessory gland (MAG) proteins to females as a solid mass (i.e. the "mating plug"). These proteins are postulated to function as important modulators of female post-mating responses. To understand the role of selective forces underlying the evolution of these proteins in the A. gambiae complex, we carried out an evolutionary analysis of gene sequence and expression divergence on a pair of paralog genes called AgAcp34A-1 and AgAcp34A-2. These encode MAG-specific proteins which, based on homology with Drosophila, have been hypothesized to play a role in sperm viability and function.

Results

Genetic analysis of 6 species of the A. gambiae complex revealed the existence of a third paralog (68-78% of identity), that we named AgAcp34A-3. FISH assays showed that this gene maps in the same division (34A) of chromosome-3R as the other two paralogs. In particular, immuno-fluorescence assays targeting the C-terminals of AgAcp34A-2 and AgAcp34A-3 revealed that these two proteins are localized in the posterior part of the MAG and concentrated at the apical portion of the mating plug. When transferred to females, this part of the plug lies in proximity to the duct connecting the spermatheca to the uterus, suggesting a potential role for these proteins in regulating sperm motility. AgAcp34A-3 is more polymorphic than the other two paralogs, possibly because of relaxation of purifying selection. Since both unequal crossing-over and gene conversion likely homogenized the members of this gene family, the interpretation of the evolutionary patterns is not straightforward. Although several haplotypes of the three paralogs are shared by most A. gambiae s.l. species, some fixed species-specific replacements (mainly placed in the N- and C-terminal portions of the secreted peptides) were also observed, suggesting some lineage-specific adaptation.

Conclusions

Progress in understanding the signaling cascade in the A. gambiae reproductive pathway will elucidate the interaction of this MAG-specific protein family with their female counterparts. This knowledge will allow a better evaluation of the relative importance of genes involved in the reproductive isolation and fertility of A. gambiae species and could help the interpretation of the observed evolutionary patterns.

Gene flow and hybridization between numerically imbalanced populations of two duck species in the Falkland Islands

Interspecific hybridization is common in plants and animals, particularly in waterfowl (Anatidae). One factor shown to contribute to hybridization is restricted mate choice, which can occur when two species occur in sympatry but one is rare. The Hubbs principle, or “desperation hypothesis,” states that under such circumstances the rarer species is more likely to mate with heterospecifics. Here we report interspecific hybridization between two waterfowl species that coexist in broad sympatry and mixed flocks throughout southern South America. Speckled teal (Anas flavirostris) and yellow-billed pintails (Anas georgica) are abundant in continental South America, but in the Falkland Islands speckled teal outnumber yellow-billed pintails approximately ten to one. Using eight genetic loci (mtDNA and 7 nuclear introns) coupled with Bayesian assignment tests and relatedness analysis, we identified a speckled teal x yellow-billed pintail F₁ hybrid female and her duckling sired by a male speckled teal. Although our sample in the Falkland Islands was small, we failed to identify unequivocal evidence of hybridization or introgression in a much larger sample from Argentina using a three-population “isolation with migration” coalescent analysis. While additional data are needed to determine if this event in the Falkland Islands was a rare singular occurrence, our results provide further support for the “desperation hypothesis,” which states that scarcity in one population and abundance of another will often lead to hybridization.

Molecular evolution of a gene cluster of serine proteases expressed in the Anopheles gambiae female reproductive tract

Background

Genes involved in post-mating processes of multiple mating organisms are known to evolve rapidly due to coevolution driven by sexual conflict among male-female interacting proteins. In the malaria mosquito Anopheles gambiae - a monandrous species in which sexual conflict is expected to be absent or minimal - recent data strongly suggest that proteolytic enzymes specifically expressed in the female lower reproductive tissues are involved in the processing of male products transferred to females during mating. In order to better understand the role of selective forces underlying the evolution of proteins involved in post-mating responses, we analysed a cluster of genes encoding for three serine proteases that are down-regulated after mating, two of which specifically expressed in the atrium and one in the spermatheca of A. gambiae females.

Results

The analysis of polymorphisms and divergence of these female-expressed proteases in closely related species of the A. gambiae complex revealed a high level of replacement polymorphisms consistent with relaxed evolutionary constraints of duplicated genes, allowing to rapidly fix novel replacements to perform new or more specific functions. Adaptive evolution was detected in several codons of the 3 genes and hints of episodic selection were also found. In addition, the structural modelling of these proteases highlighted some important differences in their substrate specificity, and provided evidence that a number of sites evolving under selective pressures lie relatively close to the catalytic triad and/or on the edge of the specificity pocket, known to be involved in substrate recognition or binding. The observed patterns suggest that these proteases may interact with factors transferred by males during mating (e.g. substrates, inhibitors or pathogens) and that they may have differently evolved in independent A. gambiae lineages.

Conclusions

Our results - also examined in light of constraints in the application of selection-inference methods to the closely related species of the A. gambiae complex - reveal an unexpectedly intricate evolutionary scenario. Further experimental analyses are needed to investigate the biological functions of these genes in order to better interpret their molecular evolution and to assess whether they represent possible targets for limiting the fertility of Anopheles mosquitoes in malaria vector control strategies.

Evolutionary patterns of two major reproduction candidate genes (Zp2 and Zp3) reveal no contribution to reproductive isolation between bovine species

BACKGROUND

It has been established that mammalian egg zona pellucida (ZP) glycoproteins are responsible for species-restricted binding of sperm to unfertilized eggs, inducing the sperm acrosome reaction, and preventing polyspermy. In mammals, ZP apparently represents a barrier to heterospecific fertilization and thus probably contributes to reproductive isolation between species. The evolutionary relationships between some members of the tribe Bovini are complex and highly debatable, particularly, those involving Bos and Bison species for which interspecific hybridization is extensively documented. Because reproductive isolation is known to be a major precursor of species divergence, testing evolutionary patterns of ZP glycoproteins may shed some light into the speciation process of these species. To this end, we have examined intraspecific and interspecific genetic variation of two ZP genes (Zp2 and Zp3) for seven representative species (111 individuals) from the Bovini tribe, including five species from Bos and Bison, and two species each from genera Bubalus and Syncerus.

RESULTS

A pattern of low levels of intraspecific polymorphism and interspecific divergence was detected for the two sequenced fragments each for Zp2 and Zp3. At intraspecific level, none of neutrality tests detected deviations from neutral equilibrium expectations for the two genes. Several haplotypes in both genes were shared by multiple species from Bos and Bison.

CONCLUSIONS

Here we argue that neither ancestral polymorphism nor introgressive hybridization alone can fully account for haplotype sharing among species from Bos and Bison, and that both scenarios have contributed to such a pattern of haplotype sharing observed here. Additionally, codon-based tests revealed strong evidence for purifying selection in the Zp3 coding haplotype sequences and weak evidence for purifying selection in the Zp2 coding haplotype sequences. Contrary to a general genetic pattern that genes or genomic regions contributing to reproductive isolation between species often evolve rapidly and show little or no gene flow between species, these results demonstrate that, particularly, those sequenced exons of the Zp2 and the Zp3 did not show any contribution to reproductive isolation between the bovine species studied here.

Evolutionary dynamics of the Ty3/gypsy LTR retrotransposons in the genome of Anopheles gambiae

Ty3/gypsy elements represent one of the most abundant and diverse LTR-retrotransposon (LTRr) groups in the Anopheles gambiae genome, but their evolutionary dynamics have not been explored in detail. Here, we conduct an in silico analysis of the distribution and abundance of the full complement of 1045 copies in the updated AgamP3 assembly. Chromosomal distribution of Ty3/gypsy elements is inversely related to arm length, with densities being greatest on the X, and greater on the short versus long arms of both autosomes. Taking into account the different heterochromatic and euchromatic compartments of the genome, our data suggest that the relative abundance of Ty3/gypsy LTRrs along each chromosome arm is determined mainly by the different proportions of heterochromatin, particularly pericentric heterochromatin, relative to total arm length. Additionally, the breakpoint regions of chromosomal inversion 2La appears to be a haven for LTRrs. These elements are underrepresented more than 7-fold in euchromatin, where 33% of the Ty3/gypsy copies are associated with genes. The euchromatin on chromosome 3R shows a faster turnover rate of Ty3/gypsy elements, characterized by a deficit of proviral sequences and the lowest average sequence divergence of any autosomal region analyzed in this study. This probably reflects a principal role of purifying selection against insertion for the preservation of longer conserved syntenyc blocks with adaptive importance located in 3R. Although some Ty3/gypsy LTRrs show evidence of recent activity, an important fraction are inactive remnants of relatively ancient insertions apparently subject to genetic drift. Consistent with these computational predictions, an analysis of the occupancy rate of putatively older insertions in natural populations suggested that the degenerate copies have been fixed across the species range in this mosquito, and also are shared with the sibling species Anopheles arabiensis.

Reconstructing the evolutionary history of the radiation of the land snail genus Xerocrassa on Crete based on mitochondrial sequences and AFLP markers

BACKGROUND

A non-adaptive radiation triggered by sexual selection resulted in ten endemic land snail species of the genus Xerocrassa on Crete. Only five of these species and a more widespread species are monophyletic in a mitochondrial gene tree. The reconstruction of the evolutionary history of such closely related species can be complicated by incomplete lineage sorting, introgression or inadequate taxonomy. To distinguish between the reasons for the nonmonophyly of several species in the mitochondrial gene tree we analysed nuclear AFLP markers.

RESULTS

Whereas six of the eleven morphologically delimited Xerocrassa species from Crete are monophyletic in the mitochondrial gene tree, nine of these species are monophyletic in the tree based on AFLP markers. Only two morphologically delimited species could not be distinguished with the multilocus data and might have diverged very recently or might represent extreme forms of a single species. The nonmonophyly of X. rhithymna with respect to X. kydonia is probably the result of incomplete lineage sorting, because there is no evidence for admixture in the AFLP data and the mitochondrial haplotype groups of these species coalesce deeply. The same is true for the main haplotype groups of X. mesostena. The nonmonophyly of X. franciscoi might be the result of mitochondrial introgression, because the coalescences of the haplotypes of this species with some X. mesostena haplotypes are shallow and there is admixture with neighbouring X. mesostena.

CONCLUSION

The most likely causes for the nonmonophyly of species in the mitochondrial gene tree of the Xerocrassa radiation on Crete could be inferred using AFLP data by a combination of several criteria, namely the depth of the coalescences in the gene tree, the geographical distribution of shared genetic markers, and concordance with results of admixture analyses of nuclear multilocus markers. The strongly subdivided population structure increases the effective population size of land snail species and, thus, the likelihood of a long persistence of ancestral polymorphisms. Our study suggests that ancestral polymorphisms are a frequent cause for nonmonophyly of species with a strongly subdivided population structure in gene trees.

Anopheles immune genes and amino acid sites evolving under the effect of positive selection

BACKGROUND

It has long been the goal of vector biology to generate genetic knowledge that can be used to "manipulate" natural populations of vectors to eliminate or lessen disease burden. While long in coming, progress towards reaching this goal has been made. Aiming to increase our understanding regarding the interactions between Plasmodium and the Anopheles immune genes, we investigated the patterns of genetic diversity of four anti-Plasmodium genes in the Anopheles gambiae complex of species.

METHODOLOGY/PRINCIPAL FINDINGS

Within a comparative phylogenetic and population genetics framework, the evolutionary history of four innate immunity genes within the An. gambiae complex (including the two most important human malaria vectors, An. gambiae and An. arabiensis) is reconstructed. The effect of natural selection in shaping the genes' diversity is examined. Introgression and retention of ancestral polymorphisms are relatively rare at all loci. Despite the potential confounding effects of these processes, we could identify sites that exhibited dN/dS ratios greater than 1.

CONCLUSIONS/SIGNIFICANCE

In two of the studied genes, CLIPB14 and FBN8, several sites indicated evolution under positive selection, with CLIPB14 exhibiting the most consistent evidence. Considering only the sites that were consistently identified by all methods, two sites in CLIPB14 are adaptively driven. However, the analysis inferring the lineage -specific evolution of each gene was not in favor of any of the Anopheles lineages evolving under the constraints imposed by positive selection. Nevertheless, the loci and the specific amino acids that were identified as evolving under strong evolutionary pressure merit further investigation for their involvement in the Anopheles defense against microbes in general.

Complete lack of mitochondrial divergence between two species of NE Atlantic marine intertidal gastropods

Some mitochondrial introgression is common between closely related species, but distinct species rarely show substantial introgression in their entire distribution range. In this study, however, we report a complete lack of mitochondrial divergence between two sympatric species of flat periwinkles (Littorina fabalis and Littorina obtusata) which, based on previous allozyme studies, diverged approximately 1 Ma. We re-examined their species status using both morphology (morphometric analysis) and neutral genetic markers (microsatellites) and our results confirmed that these species are well separated. Despite this, the two species shared all common cytochrome-b haplotypes throughout their NE Atlantic distribution and no deep split between typical L. fabalis and L. obtusata haplotypes could be found. We suggest that incomplete lineage sorting explains most of the lack of mitochondrial divergence between these species. However, coalescent-based analyses and the sympatric sharing of unique haplotypes suggest that introgressive hybridization also has occurred.

Paraphyly of organelle DNAs in Cycas Sect. Asiorientales due to ancient ancestral polymorphisms

Background

This study addresses the apportionment of genetic diversity between Cycas revoluta and C. taitungensis, species that constitute the section Asiorientales and represent a unique, basal lineage of the Laurasian genus Cycas. Fossil evidence indicates divergence of the section from the rest of Cycas at least 30 million years ago. Geographically, C. taitungensis is limited to Taiwan whereas C. revoluta is found in the Ryukyu Archipelago and on mainland China.

Results

The phylogenies of ribosomal ITS region of mtDNA and the intergenic spacer between atpB and rbcL genes of cpDNA were reconstructed. Phylogenetic analyses revealed paraphyly of both loci in the two species and also in the section Asiorientales. The lack of reciprocal monophyly between these long isolated sections is likely due to persistent shared ancestral polymorphisms. Molecular dating estimated that mt- and cp DNA lineages coalesced to the most recent common ancestors (TMRCA) about 327 (mt) and 204 MYA (cp), corresponding with the divergence of cycad sections in the Mesozoic.

Conclusion

Fates of newly derived mutations of cycads follow Klopfstein et al.'s surfing model where the majority of new mutations do not spread geographically and remain at low frequencies or are eventually lost by genetic drift. Only successful 'surfing mutations' reach very high frequencies and occupy a large portion of a species range. These mutations exist as dominant cytotypes across populations and species. Geographical subdivision is lacking in both species, even though recurrent gene flow by both pollen and seed is severely limited. In total, the contrasting levels between historical and ongoing gene flow, large population sizes, a long lifespan, and slow mutation rates in both organelle DNAs have all likely contributed to the unusually long duration of paraphyly in cycads.

Inferring selection in the Anopheles gambiae species complex: an example from immune-related serine protease inhibitors

Background

Mosquitoes of the Anopheles gambiae species complex are the primary vectors of human malaria in sub-Saharan Africa. Many host genes have been shown to affect Plasmodium development in the mosquito, and so are expected to engage in an evolutionary arms race with the pathogen. However, there is little conclusive evidence that any of these mosquito genes evolve rapidly, or show other signatures of adaptive evolution.

Methods

Three serine protease inhibitors have previously been identified as candidate immune system genes mediating mosquito-Plasmodium interaction, and serine protease inhibitors have been identified as hot-spots of adaptive evolution in other taxa. Population-genetic tests for selection, including a recent multi-gene extension of the McDonald-Kreitman test, were applied to 16 serine protease inhibitors and 16 other genes sampled from the An. gambiae species complex in both East and West Africa.

Results

Serine protease inhibitors were found to show a marginally significant trend towards higher levels of amino acid diversity than other genes, and display extensive genetic structuring associated with the 2La chromosomal inversion. However, although serpins are candidate targets for strong parasite-mediated selection, no evidence was found for rapid adaptive evolution in these genes.

Conclusion

It is well known that phylogenetic and population history in the An. gambiae complex can present special problems for the application of standard population-genetic tests for selection, and this may explain the failure of this study to detect selection acting on serine protease inhibitors. The pitfalls of uncritically applying these tests in this species complex are highlighted, and the future prospects for detecting selection acting on the An. gambiae genome are discussed.

Introgression as a likely cause of mtDNA paraphyly in two allopatric skippers (Lepidoptera: Hesperiidae)

Gene transfer between species during interspecific hybridization is a widely accepted reality in plants but is considered a relatively rare phenomenon among animals. Here we describe a unique case of mitochondrial DNA (mtDNA) paraphyly in the skipper genus, Erynnis, that involves well-diverged allopatric species. Using molecular evidence from both mitochondrial and nuclear genomes, we found high levels of intraspecific divergence in the mitochondrial genome within E. propertius (over 4% pair-wise sequence divergence) but no such differentiation in the nuclear genome. Sequence comparisons with related Erynnis suggest that past, but recent and infrequent introgression between E. propertius and E. horatius is the most reasonable explanation for the observed pattern of mtDNA paraphyly. This example of putative introgression highlights the complexity of mtDNA evolution and suggests that similar processes could be operating in other taxa that have not been extensively sampled. Our observations reinforce the importance of involving multiple genes with different modes of inheritance in the analysis of population history of congeneric taxa.