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

HM, Antoshechkin I, Marshall JM, Akbari OS. Eliminating malaria vectors with precision-guided sterile males

Controlling the principal African malaria vector, the mosquito Anopheles gambiae, is considered essential to curtail malaria transmission. However, existing vector control technologies rely on insecticides, which are becoming increasingly ineffective. Sterile insect technique (SIT) is a powerful suppression approach that has successfully eradicated a number of insect pests, yet the A. gambiae toolkit lacks the requisite technologies for its implementation. SIT relies on iterative mass releases of nonbiting, nondriving, sterile males which seek out and mate with monandrous wild females. Once mated, females are permanently sterilized due to mating-induced refractoriness, which results in population suppression of the subsequent generation. However, sterilization by traditional methods renders males unfit, making the creation of precise genetic sterilization methods imperative. Here, we introduce a vector control technology termed precision-guided sterile insect technique (pgSIT), in A. gambiae for inducible, programmed male sterilization and female elimination for wide-scale use in SIT campaigns. Using a binary CRISPR strategy, we cross separate engineered Cas9 and gRNA strains to disrupt male-fertility and female-essential genes, yielding >99.5% male sterility and >99.9% female lethality in hybrid progeny. We demonstrate that these genetically sterilized males have good longevity, are able to induce sustained population suppression in cage trials, and are predicted to eliminate wild A. gambiae populations using mathematical models, making them ideal candidates for release. This work provides a valuable addition to the malaria genetic biocontrol toolkit, enabling scalable SIT-like confinable, species-specific, and safe suppression in the species.

HM, Antoshechkin I, Marshall JM, Akbari OS. Eliminating Malaria Vectors with Precision Guided Sterile Males

Controlling the principal African malaria vector, the mosquito Anopheles gambiae, is considered essential to curtail malaria transmission. However existing vector control technologies rely on insecticides, which are becoming increasingly ineffective. Sterile insect technique (SIT) is a powerful suppression approach that has successfully eradicated a number of insect pests, yet the A. gambiae toolkit lacks the requisite technologies for its implementation. SIT relies on iterative mass-releases of non-biting, non-driving, sterile males which seek out and mate with monandrous wild females. Once mated, females are permanently sterilized due to mating-induced refractoriness, which results in population suppression of the subsequent generation. However, sterilization by traditional methods renders males unfit, making the creation of precise genetic sterilization methods imperative. Here we develop precision guided Sterile Insect Technique (pgSIT) in the mosquito A. gambiae for inducible, programmed male-sterilization and female-elimination for wide scale use in SIT campaigns. Using a binary CRISPR strategy, we cross separate engineered Cas9 and gRNA strains to disrupt male-fertility and female-essential genes, yielding >99.5% male-sterility and >99.9% female-lethality in hybrid progeny. We demonstrate that these genetically sterilized males have good longevity, are able to induce population suppression in cage trials, and are predicted to eliminate wild A. gambiae populations using mathematical models, making them ideal candidates for release. This work provides a valuable addition to the malaria genetic biocontrol toolkit, for the first time enabling scalable SIT-like confinable suppression in the species.

Balancing selection on an MYB transcription factor maintains the twig trichome color variation in Melastoma normale

BACKGROUND

The factors that maintain phenotypic and genetic variation within a population have received long-term attention in evolutionary biology. Here the genetic basis and evolution of the geographically widespread variation in twig trichome color (from red to white) in a shrub Melastoma normale was investigated using Pool-seq and evolutionary analyses.

RESULTS

The results show that the twig trichome coloration is under selection in different light environments and that a 6-kb region containing an R2R3 MYB transcription factor gene is the major region of divergence between the extreme red and white morphs. This gene has two highly divergent groups of alleles, one of which likely originated from introgression from another species in this genus and has risen to high frequency (> 0.6) within each of the three populations under investigation. In contrast, polymorphisms in other regions of the genome show no sign of differentiation between the two morphs, suggesting that genomic patterns of diversity have been shaped by homogenizing gene flow. Population genetics analysis reveals signals of balancing selection acting on this gene, and it is suggested that spatially varying selection is the most likely mechanism of balancing selection in this case.

CONCLUSIONS

This study demonstrate that polymorphisms on a single transcription factor gene largely confer the twig trichome color variation in M. normale, while also explaining how adaptive divergence can occur and be maintained in the face of gene flow.

Fine-scale spatial distribution of deltamethrin resistance and population structure of Anopheles funestus and Anopheles arabiensis populations in Southern Mozambique

BACKGROUND

Insecticide resistance in malaria vectors can be spatially highly heterogeneous, yet population structure analyses frequently find relatively high levels of gene flow among mosquito populations. Few studies have contemporaneously assessed phenotypic, genotypic and population structure analysis on mosquito populations and none at fine geographical scales. In this study, genetic diversity, population structure, and insecticide resistance profiles of Anopheles funestus and Anopheles arabiensis were examined across mosquito populations from and within neighbouring villages.

METHODS

Mosquitoes were collected from 11 towns in southern Mozambique, as well as from different neighbourhoods within the town of Palmeira, during the peak malaria transmission season in 2016. CDC bottle bioassay and PCR assays were performed with Anopheles mosquitoes at each site to determine phenotypic and molecular insecticide resistance profiles, respectively. Microsatellite analysis was conducted on a subsample of mosquitoes to estimate genetic diversity and population structure.

RESULTS

Phenotypic insecticide resistance to deltamethrin was observed in An. funestus sensu stricto (s.s.) throughout the area, though a high level of mortality variation was seen. However, 98% of An. funestus s.s. were CYP6P9a homozygous resistant. An. arabiensis was phenotypically susceptible to deltamethrin and 99% were kdr homozygous susceptible. Both Anopheles species exhibited high allelic richness and heterozygosity. Significant deviations from Hardy-Weinberg equilibrium were observed, and high linkage disequilibrium was seen for An. funestus s.s., supporting population subdivision. However, the FST values were low for both anophelines (- 0.00457 to 0.04213), Nm values were high (9.4-71.8 migrants per generation), AMOVA results showed almost 100% genetic variation among and within individuals, and Structure analysis showed no clustering of An. funestus s.s. and An. arabiensis populations. These results suggest high gene flow among mosquito populations.

CONCLUSION

Despite a relatively high level of phenotypic variation in the An. funestus population, molecular analysis shows the population is admixed. These data indicate that CYP6P9a resistance markers do not capture all phenotypic variation in the area, but also that resistance genes of high impact are likely to easily spread in the area. Conversely, other strategies, such as transgenic mosquito release programmes will likely not face challenges in this locality.

Gene drive in species complexes: defining target organisms

Engineered gene drives, which bias their own inheritance to increase in frequency in target populations, are being developed to control mosquito malaria vectors. Such mosquitoes can belong to complexes of both vector and nonvector species that can produce fertile interspecific hybrids, making vertical gene drive transfer (VGDT) to sibling species biologically plausible. While VGDT to other vectors could positively impact human health protection goals, VGDT to nonvectors might challenge biodiversity ones. Therefore, environmental risk assessment of gene drive use in species complexes invites more nuanced considerations of target organisms and nontarget organisms than for transgenes not intended to increase in frequency in target populations. Incorporating the concept of target species complexes offers more flexibility when assessing potential impacts from VGDT.

Ghost lineages can invalidate or even reverse findings regarding gene flow

Introgression, endosymbiosis, and gene transfer, i.e., horizontal gene flow (HGF), are primordial sources of innovation in all domains of life. Our knowledge on HGF relies on detection methods that exploit some of its signatures left on extant genomes. One of them is the effect of HGF on branch lengths of constructed phylogenies. This signature has been formalized in statistical tests for HGF detection and used for example to detect massive adaptive gene flows in malaria vectors or to order evolutionary events involved in eukaryogenesis. However, these studies rely on the assumption that ghost lineages (all unsampled extant and extinct taxa) have little influence. We demonstrate here with simulations and data reanalysis that when considering the more realistic condition that unsampled taxa are legion compared to sampled ones, the conclusion of these studies become unfounded or even reversed. This illustrates the necessity to recognize the existence of ghosts in evolutionary studies.

Recommendations for environmental risk assessment of gene drive applications for malaria vector control

Building on an exercise that identified potential harms from simulated investigational releases of a population suppression gene drive for malaria vector control, a series of online workshops identified nine recommendations to advance future environmental risk assessment of gene drive applications.

Genomic insights into variation in thermotolerance between hybridizing swordtail fishes

Understanding how organisms adapt to changing environments is a core focus of research in evolutionary biology. One common mechanism is adaptive introgression, which has received increasing attention as a potential route to rapid adaptation in populations struggling in the face of ecological change, particularly global climate change. However, hybridization can also result in deleterious genetic interactions that may limit the benefits of adaptive introgression. Here, we used a combination of genome-wide quantitative trait locus mapping and differential gene expression analyses between the swordtail fish species Xiphophorus malinche and X. birchmanni to study the consequences of hybridization on thermotolerance. While these two species are adapted to different thermal environments, we document a complicated architecture of thermotolerance in hybrids. We identify a region of the genome that contributes to reduced thermotolerance in individuals heterozygous for X. malinche and X. birchmanni ancestry, as well as widespread misexpression in hybrids of genes that respond to thermal stress in the parental species, particularly in the circadian clock pathway. We also show that a previously mapped hybrid incompatibility between X. malinche and X. birchmanni contributes to reduced thermotolerance in hybrids. Together, our results highlight the challenges of understanding the impact of hybridization on complex ecological traits and its potential impact on adaptive introgression.

Unravelling population structure heterogeneity within the genome of the malaria vector Anopheles gambiae

Background

Whole genome re-sequencing provides powerful data for population genomic studies, allowing robust inferences of population structure, gene flow and evolutionary history. For the major malaria vector in Africa, Anopheles gambiae, other genetic aspects such as selection and adaptation are also important. In the present study, we explore population genetic variation from genome-wide sequencing of 765 An. gambiae and An. coluzzii specimens collected from across Africa. We used t-SNE, a recently popularized dimensionality reduction method, to create a 2D-map of An. gambiae and An. coluzzii genes that reflect their population structure similarities.

Results

The map allows intuitive navigation among genes distributed throughout the so-called “mainland” and numerous surrounding “island-like” gene clusters. These gene clusters of various sizes correspond predominantly to low recombination genomic regions such as inversions and centromeres, and also to recent selective sweeps. Because this mosquito species complex has been studied extensively, we were able to support our interpretations with previously published findings. Several novel observations and hypotheses are also made, including selective sweeps and a multi-locus selection event in Guinea-Bissau, a known intense hybridization zone between An. gambiae and An. coluzzii.

Conclusions

Our results present a rich dataset that could be utilized in functional investigations aiming to shed light onto An. gambiae s.l genome evolution and eventual speciation. In addition, the methodology presented here can be used to further characterize other species not so well studied as An. gambiae, shortening the time required to progress from field sampling to the identification of genes and genomic regions under unique evolutionary processes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07722-y.

Chromosomal evolution in Raphicerus antelope suggests divergent X chromosomes may drive speciation through females, rather than males, contrary to Haldane's rule

Chromosome structural change has long been considered important in the evolution of post-zygotic reproductive isolation. The premise that karyotypic variation can serve as a possible barrier to gene flow is founded on the expectation that heterozygotes for structurally distinct chromosomal forms would be partially sterile (negatively heterotic) or show reduced recombination. We report the outcome of a detailed comparative molecular cytogenetic study of three antelope species, genus Raphicerus, that have undergone a rapid radiation. The species are largely conserved with respect to their euchromatic regions but the X chromosomes, in marked contrast, show distinct patterns of heterochromatic amplification and localization of repeats that have occurred independently in each lineage. We argue a novel hypothesis that postulates that the expansion of heterochromatic blocks in the homogametic sex can, with certain conditions, contribute to post-zygotic isolation. i.e., female hybrid incompatibility, the converse of Haldane’s rule. This is based on the expectation that hybrids incur a selective disadvantage due to impaired meiosis resulting from the meiotic checkpoint network’s surveillance of the asymmetric expansions of heterochromatic blocks in the homogametic sex. Asynapsis of these heterochromatic regions would result in meiotic silencing of unsynapsed chromatin and, if this persists, germline apoptosis and female infertility.

A phylogenomic study of Steganinae fruit flies (Diptera: Drosophilidae): strong gene tree heterogeneity and evidence for monophyly

BACKGROUND

The Drosophilidae family is traditionally divided into two subfamilies: Drosophilinae and Steganinae. This division is based on morphological characters, and the two subfamilies have been treated as monophyletic in most of the literature, but some molecular phylogenies have suggested Steganinae to be paraphyletic. To test the paraphyletic-Steganinae hypothesis, here, we used genomic sequences of eight Drosophilidae (three Steganinae and five Drosophilinae) and two Ephydridae (outgroup) species and inferred the phylogeny for the group based on a dataset of 1,028 orthologous genes present in all species (> 1,000,000 bp). This dataset includes three genera that broke the monophyly of the subfamilies in previous works. To investigate possible biases introduced by small sample sizes and automatic gene annotation, we used the same methods to infer species trees from a set of 10 manually annotated genes that are commonly used in phylogenetics.

RESULTS

Most of the 1,028 gene trees depicted Steganinae as paraphyletic with distinct topologies, but the most common topology depicted it as monophyletic (43.7% of the gene trees). Despite the high levels of gene tree heterogeneity observed, species tree inference in ASTRAL, in PhyloNet, and with the concatenation approach strongly supported the monophyly of both subfamilies for the 1,028-gene dataset. However, when using the concatenation approach to infer a species tree from the smaller set of 10 genes, we recovered Steganinae as a paraphyletic group. The pattern of gene tree heterogeneity was asymmetrical and thus could not be explained solely by incomplete lineage sorting (ILS).

CONCLUSIONS

Steganinae was clearly a monophyletic group in the dataset that we analyzed. In addition to ILS, gene tree discordance was possibly the result of introgression, suggesting complex branching processes during the early evolution of Drosophilidae with short speciation intervals and gene flow. Our study highlights the importance of genomic data in elucidating contentious phylogenetic relationships and suggests that phylogenetic inference for drosophilids based on small molecular datasets should be performed cautiously. Finally, we suggest an approach for the correction and cleaning of BUSCO-derived genomic datasets that will be useful to other researchers planning to use this tool for phylogenomic studies.

Evaluating the probability of CRISPR-based gene drive contaminating another species

The probability D that a given clustered regularly interspaced short palindromic repeats (CRISPR)-based gene drive element contaminates another, nontarget species can be estimated by the following Drive Risk Assessment Quantitative Estimate (DRAQUE) Equation: D = h y b + t r a n s f × e x p r e s s × c u t × f l a n k × i m m u n e × n o n e x t i n c t with hyb = probability of hybridization between the target species and a nontarget species; transf = probability of horizontal transfer of a piece of DNA containing the gene drive cassette from the target species to a nontarget species (with no hybridization); express = probability that the Cas9 and guide RNA genes are expressed; cut = probability that the CRISPR-guide RNA recognizes and cuts at a DNA site in the new host; flank = probability that the gene drive cassette inserts at the cut site; immune = probability that the immune system does not reject Cas9-expressing cells; nonextinct = probability of invasion of the drive within the population. We discuss and estimate each of the seven parameters of the equation, with particular emphasis on possible transfers within insects, and between rodents and humans. We conclude from current data that the probability of a gene drive cassette to contaminate another species is not insignificant. We propose strategies to reduce this risk and call for more work on estimating all the parameters of the formula.

Genetic analysis and population structure of the Anopheles gambiae complex from different ecological zones of Burkina Faso

The Anopheles gambiae complex (Diptera: Culicidae) is the most important vector for malaria in Sub-Saharan Africa, besides other vectors such as Anopheles funestus. Malaria vector control should encompass specific identification, genetic diversity and population structure of An. gambiae to design vector control strategies. The aim of this study was to determine the distribution of sibling species of the An. gambiae complex according to climatic regions related to cotton-growing or cotton-free areas by using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). Then, variation in mitochondrial cytochrome c oxidase 1 (COI) was used to assess the genetic structure within and between populations from our selected ecological zones. At the sibling species level, the following proportions were found across all samples (n = 180): An. coluzzii 65.56%, An. gambiae stricto sensu (s.s). 21.11%, and An. arabiensis 3.33%. Hybrids between An. gambiae s.s. and An. coluzzii (7.78%) and hybrids between An. coluzzii and An. arabiensis (2.22%) were found. The phylogenetic tree and Integer Neighbour-Joining (IntNJ) haplotype network did not reveal any distinct genetic structure pattern related to climatic or agricultural conditions in Burkina Faso. The Fst (Wright's F-statistic) values close to zero showed a free gene flow and no differentiation in An. gambiae complex populations. Furthermore, neutrality indices calculated by Tajima's D, Fu and Li's D^⁎, Fu and Li's F^⁎, Fu's Fs tests suggested an excess of rare mutations in the populations. Overall, variation in the proportions of An. gambiae s.s., An. coluzzii and An. arabiensis was found according to climatic regions, but COI analysis did not evidence any population structuring of the An. gambiae complex. These scientific contributions can be used as a basis for further in-depth study of the genetic diversity of the An. gambiae complex for epidemiological risk assessment of malaria in Burkina Faso.

Agricultural insect hybridization and implications for pest management

Biological invasions, the expansion of agricultural frontiers, and climate change favor encounters of divergent lineages of animals and plants, increasing the likelihood of hybridization. However, hybridization of insect species and its consequences for agroecosystems have not received sufficient attention. Gene exchange between distinct and distant genetic pools can improve the survival and reproduction of insect pests, and threaten beneficial insects in disturbed agricultural environments. Hybridization may be the underlying explanation for the recurrent pest outbreaks and control failures in putative hybrid zones, as suspected for bollworm, corn borer, whiteflies, and stink bugs. Reliable predictions of the types of changes that can be expected in pest insect genomes and fitness, and of their impacts on the fate of species and populations remain elusive. Typical steps in pest management, such as insect identification, pest monitoring, and control are likely affected by gene flow and adaptive introgression mediated by hybridization, and we do not have ways to respond to or mitigate the problem. To address the adverse effects of farming intensification and global trade, we must ensure that current integrated pest management programs incorporate up-to-date monitoring and diagnostic tools. The rapid identification of hybrids, quantification of levels of introgression, and in-depth knowledge of what genes have been transferred may help to explain and predict insect population outbreaks and control failures in the future. © 2019 Society of Chemical Industry.

A new species in the major malaria vector complex sheds light on reticulated species evolution

Complexes of closely related species provide key insights into the rapid and independent evolution of adaptive traits. Here, we described and studied Anopheles fontenillei sp.n., a new species in the Anopheles gambiae complex that we recently discovered in the forested areas of Gabon, Central Africa. Our analysis placed the new taxon in the phylogenetic tree of the An. gambiae complex, revealing important introgression events with other members of the complex. Particularly, we detected recent introgression, with Anopheles gambiae and Anopheles coluzzii, of genes directly involved in vectorial capacity. Moreover, genome analysis of the new species allowed us to clarify the evolutionary history of the 3La inversion. Overall, An. fontenillei sp.n. analysis improved our understanding of the relationship between species within the An. gambiae complex, and provided insight into the evolution of vectorial capacity traits that are relevant for the successful control of malaria in Africa.

Host plant-related genomic differentiation in the European cherry fruit fly, Rhagoletis cerasi

Elucidating the mechanisms and conditions facilitating the formation of biodiversity are central topics in evolutionary biology. A growing number of studies imply that divergent ecological selection may often play a critical role in speciation by counteracting the homogenising effects of gene flow. Several examples involve phytophagous insects, where divergent selection pressures associated with host plant shifts may generate reproductive isolation, promoting speciation. Here, we use ddRADseq to assess the population structure and to test for host‐related genomic differentiation in the European cherry fruit fly, Rhagoletis cerasi (L., 1758) (Diptera: Tephritidae). This tephritid is distributed throughout Europe and western Asia, and has adapted to two different genera of host plants, Prunus spp. (cherries) and Lonicera spp. (honeysuckle). Our data imply that geographic distance and geomorphic barriers serve as the primary factors shaping genetic population structure across the species range. Locally, however, flies genetically cluster according to host plant, with consistent allele frequency differences displayed by a subset of loci between Prunus and Lonicera flies across four sites surveyed in Germany and Norway. These 17 loci display significantly higher F_ST values between host plants than others. They also showed high levels of linkage disequilibrium within and between Prunus and Lonicera flies, supporting host‐related selection and reduced gene flow. Our findings support the existence of sympatric host races in R. cerasi embedded within broader patterns of geographic variation in the fly, similar to the related apple maggot, Rhagoletis pomonella, in North America.

Is Genetic Continuity Between Anopheles sinensis (Diptera: Culicidae) and its Sibling Species Due to Gene Introgression or Incomplete Speciation?

The Anopheles mosquito Hyrcanus Group is widely distributed geographically across both Palearctic and Oriental regions and comprises 26 valid species. Although the species Anopheles sinensis Wiedemann (1828) is the most common in China and has a low potential vector rank, it has nevertheless long been thought to be an important natural malaria vector within the middle and lower reaches of the Yangtze River. A number of previous research studies have found evidence to support the occurrence of natural hybridization between An. sinensis and Anopheles kleini Rueda, 2005 (a competent malaria vector). We, therefore, collected a sample series of An. sinensis and morphologically similar species across China and undertook ribosomal and mitochondrial DNA analyses in order to assess genetic differentiation (Fst) and gene flow (Nm) amongst different groups. This enabled us to evaluate divergence times between morphologically similar species using the cytochrome oxidase I (COI) gene. The results of this study reveal significant genetic similarities between An. sinensis, An. kleini, and Anopheles belenrae Rueda, 2005 and therefore imply that correct molecular identifications will require additional molecular markers. As results also reveal the presence of gene flow between these three species, their taxonomic status will require further work. Data suggest that An. kleini is the most basal of the three species, while An. sinensis and An. belenrae share the closest genetic relationship.

Coalescent Analysis of Phylogenomic Data Confidently Resolves the Species Relationships in the Anopheles gambiae Species Complex

Deep coalescence and introgression make it challenging to infer phylogenetic relationships among closely related species that arose through radiative speciation events. Despite numerous phylogenetic analyses and the availability of whole genomes, the phylogeny in the Anopheles gambiae species complex has not been confidently resolved. Here we extract over 80, 000 coding and noncoding short segments (called loci) from the genomes of six members of the species complex and use a Bayesian method under the multispecies coalescent model to infer the species tree, which takes into account genealogical heterogeneity across the genome and uncertainty in the gene trees. We obtained a robust estimate of the species tree from the distal region of the X chromosome: (A. merus, ((A. melas, (A. arabiensis, A. quadriannulatus)), (A. gambiae, A. coluzzii))), with A. merus to be the earliest branching species. This species tree agrees with the chromosome inversion phylogeny and provides a parsimonious interpretation of inversion and introgression events. Simulation informed by the real data suggest that the coalescent approach is reliable while the sliding-window analysis used in a previous phylogenomic study generates artifactual species trees. Likelihood ratio test of gene flow revealed strong evidence of autosomal introgression from A. arabiensis into A. gambiae (at the average rate of ∼0.2 migrants per generation), but not in the opposite direction, and introgression of the 3 L chromosomal region from A. merus into A. quadriannulatus. Our results highlight the importance of accommodating incomplete lineage sorting and introgression in phylogenomic analyses of species that arose through recent radiative speciation events.

Forecasting Zoonotic Infectious Disease Response to Climate Change: Mosquito Vectors and a Changing Environment

Infectious diseases are changing due to the environment and altered interactions among hosts, reservoirs, vectors, and pathogens. This is particularly true for zoonotic diseases that infect humans, agricultural animals, and wildlife. Within the subset of zoonoses, vector-borne pathogens are changing more rapidly with climate change, and have a complex epidemiology, which may allow them to take advantage of a changing environment. Most mosquito-borne infectious diseases are transmitted by mosquitoes in three genera: Aedes, Anopheles, and Culex, and the expansion of these genera is well documented. There is an urgent need to study vector-borne diseases in response to climate change and to produce a generalizable approach capable of generating risk maps and forecasting outbreaks. Here, we provide a strategy for coupling climate and epidemiological models for zoonotic infectious diseases. We discuss the complexity and challenges of data and model fusion, baseline requirements for data, and animal and human population movement. Disease forecasting needs significant investment to build the infrastructure necessary to collect data about the environment, vectors, and hosts at all spatial and temporal resolutions. These investments can contribute to building a modeling community around the globe to support public health officials so as to reduce disease burden through forecasts with quantified uncertainty.

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.

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.

Demographic history influences spatial patterns of genetic diversityin recently expanded coyote (Canis latrans) populations

Human-mediated range expansions have increased in recent decades and represent unique opportunities to evaluate genetic outcomes of establishing peripheral populations across broad expansion fronts. Over the past century, coyotes (Canis latrans) have undergone a pervasive range expansion and now inhabit every state in the continental United States. Coyote expansion into eastern North America was facilitated by anthropogenic landscape changes and followed two broad expansion fronts. The northern expansion extended through the Great Lakes region and southern Canada, where hybridization with remnant wolf populations was common. The southern and more recent expansion front occurred approximately 40 years later and across territory where gray wolves have been historically absent and remnant red wolves were extirpated in the 1970s. We conducted a genetic survey at 10 microsatellite loci of 482 coyotes originating from 11 eastern U.S. states to address how divergent demographic histories influence geographic patterns of genetic diversity. We found that population structure corresponded to a north-south divide, which is consistent with the two known expansion routes. Additionally, we observed extremely high genetic diversity, which is atypical of recently expanded populations and is likely the result of multiple complex demographic processes, in addition to hybridization with other Canis species. Finally, we considered the transition of allele frequencies across geographic space and suggest the mid-Atlantic states of North Carolina and Virginia as an emerging contact zone between these two distinct coyote expansion fronts.

Hybridization as a facilitator of species range expansion

Explaining the evolution of species geographical ranges is fundamental to understanding how biodiversity is distributed and maintained. The solution to this classic problem in ecology and evolution remains elusive: we still do not fully know how species geographical ranges evolve and what factors fuel range expansions. Resolving this problem is now more crucial than ever with increasing biodiversity loss, global change and movement of species by humans. Here, we describe and evaluate the hypothesis that hybridization between species can contribute to species range expansion. We discuss how such a process can occur and the empirical data that are needed to test this hypothesis. We also examine how species can expand into new environments via hybridization with a resident species, and yet remain distinct species. Generally, hybridization may play an underappreciated role in influencing the evolution of species ranges. Whether-and to what extent-hybridization has such an effect requires further study across more diverse taxa.

Genetic variation during range expansion: effects of habitat novelty and hybridization

How species' ranges evolve remains an enduring problem in ecology and evolutionary biology. Species' range limits are potentially set by the inability of peripheral populations to adapt to range-edge habitat. Indeed, peripheral populations are often assumed to have reduced genetic diversity and population sizes, which limit evolvability. However, support for this assumption is mixed, possibly because the genetic effects of range expansion depend on two factors: the extent that habitat into which expansion occurs is novel and sources of gene flow. Here, we used spadefoot toads, Spea bombifrons, to contrast the population genetic effects of expansion into novel versus non-novel habitat. We further evaluated gene flow from conspecifics and from heterospecifics via hybridization with a resident species. We found that range expansion into novel habitat, relative to non-novel habitat, resulted in higher genetic differentiation, lower conspecific gene flow and bottlenecks. Moreover, we found that hybridizing with a resident species introduced genetic diversity in the novel habitat. Our results suggest the evolution of species' ranges can depend on the extent of differences in habitat between ancestral and newly occupied ranges. Furthermore, our results highlight the potential for hybridization with a resident species to enhance genetic diversity during expansions into novel habitat.

Anopheles darlingi polytene chromosomes: revised maps including newly described inversions and evidence for population structure in Manaus

Salivary gland polytene chromosomes of 4th instar Anopheles darlingi Root were examined from multiple locations in the Brazilian Amazon. Minor modifications were made to existing polytene photomaps. These included changes to the breakpoint positions of several previously described paracentric inversions and descriptions of four new paracentric inversions, two on the right arm of chromosome 3 and two on the left arm of chromosome 3 that were found in multiple locations. A total of 18 inversions on the X (n = 1) chromosome, chromosome 2 (n = 7) and 3 (n = 11) were scored for 83 individuals from Manaus, Macapá and Porto Velho municipalities. The frequency of 2Ra inversion karyotypes in Manaus shows significant deficiency of heterozygotes (p < 0.0009). No significant linkage disequilibrium was found between inversions on chromosome 2 and 3. We hypothesize that at least two sympatric subpopulations exist within the An. darlingi population at Manaus based on inversion frequencies.

Microsatellite markers reveal low frequency of natural hybridization between the white-footed mouse (Peromyscus leucopus) and deer mouse (Peromyscus maniculatus) in southern Quebec, Canada

In some parts of southern Quebec, two closely related rodent species - the white-footed mouse (Peromyscus leucopus) and the deer mouse (Peromyscus maniculatus) - have recently come in contact because of climate-driven changes in the distribution of the former. Both species share similar morphology, ecology, and life history traits, which suggests that natural hybridization may be possible. Hybridization among these two species can have important implications on the ecological roles these rodents play in disease transmission, yet few researchers have attempted to examine this phenomenon and results from previous hybridization experiments have remained inconclusive and conflicting. In this study, we attempt to investigate the occurrence of hybridization among white-footed mice and deer mice in southern Quebec by genotyping wild caught specimens with selectively neutral, polymorphic microsatellite markers. Our analyses suggest that hybridization may be occurring at extremely low frequency between both species in our study area. The presence of such hybridization events, even at low frequencies, may have implications on disease transmission risk in the region and further detailed studies are necessary.

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.

A CRISPR-Cas9 sex-ratio distortion system for genetic control

Genetic control aims to reduce the ability of insect pest populations to cause harm via the release of modified insects. One strategy is to bias the reproductive sex ratio towards males so that a population decreases in size or is eliminated altogether due to a lack of females. We have shown previously that sex ratio distortion can be generated synthetically in the main human malaria vector Anopheles gambiae, by selectively destroying the X-chromosome during spermatogenesis, through the activity of a naturally-occurring endonuclease that targets a repetitive rDNA sequence highly-conserved in a wide range of organisms. Here we describe a CRISPR-Cas9 sex distortion system that targets ribosomal sequences restricted to the member species of the Anopheles gambiae complex. Expression of Cas9 during spermatogenesis resulted in RNA-guided shredding of the X-chromosome during male meiosis and produced extreme male bias among progeny in the absence of any significant reduction in fertility. The flexibility of CRISPR-Cas9 combined with the availability of genomic data for a range of insects renders this strategy broadly applicable for the species-specific control of any pest or vector species with an XY sex-determination system by targeting sequences exclusive to the female sex chromosome.

The Evolution of Genome Structure by Natural and Sexual Selection

Progress on understanding how genome structure evolves is accelerating with the arrival of new genomic, comparative, and theoretical approaches. This article reviews progress in understanding how chromosome inversions and sex chromosomes evolve, and how their evolution affects species' ecology. Analyses of clines in inversion frequencies in flies and mosquitoes imply strong local adaptation, and roles for both over- and under dominant selection. Those results are consistent with the hypothesis that inversions become established when they capture locally adapted alleles. Inversions can carry alleles that are beneficial to closely related species, causing them to introgress following hybridization. Models show that this "adaptive cassette" scenario can trigger large range expansions, as recently happened in malaria mosquitoes. Sex chromosomes are the most rapidly evolving genome regions of some taxa. Sexually antagonistic selection may be the key force driving transitions of sex determination between different pairs of chromosomes and between XY and ZW systems. Fusions between sex-chromosomes and autosomes most often involve the Y chromosome, a pattern that can be explained if fusions are mildly deleterious and fix by drift. Sexually antagonistic selection is one of several hypotheses to explain the recent discovery that the sex determination system has strong effects on the adult sex ratios of tetrapods. The emerging view of how genome structure evolves invokes a much richer constellation of forces than was envisioned during the Golden Age of research on Drosophila karyotypes.

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.

A standard photomap of ovarian nurse cell chromosomes in the European malaria vector Anopheles atroparvus

Anopheles atroparvus (Diptera: Culicidae) is one of the main malaria vectors of the Maculipennis group in Europe. Cytogenetic analysis based on salivary gland chromosomes has been used in taxonomic and population genetic studies of mosquitoes from this group. However, a high-resolution cytogenetic map that could be used in physical genome mapping in An. atroparvus is still lacking. In the present study, a high-quality photomap of the polytene chromosomes from ovarian nurse cells of An. atroparvus was developed. Using fluorescent in situ hybridization, 10 genes from the five largest genomic supercontigs on the polytene chromosome were localized and 28% of the genome was anchored to the cytogenetic map. The study established chromosome arm homology between An. atroparvus and the major African malaria vector Anopheles gambiae, suggesting a whole-arm translocation between autosomes of these two species. The standard photomap constructed for ovarian nurse cell chromosomes of An. atroparvus will be useful for routine physical mapping. This map will assist in the development of a fine-scale chromosome-based genome assembly for this species and will also facilitate comparative and evolutionary genomics studies in the genus Anopheles.

Species delimitation in plants using the Qinghai-Tibet Plateau endemic Orinus (Poaceae: Tridentinae) as an example

BACKGROUND AND AIMS

Accurate identification of species is essential for the majority of biological studies. However, defining species objectively and consistently remains a challenge, especially for plants distributed in remote regions where there is often a lack of sufficient previous specimens. In this study, multiple approaches and lines of evidence were used to determine species boundaries for plants occurring in the Qinghai-Tibet Plateau, using the genus Orinus (Poaceae) as a model system for an integrative approach to delimiting species.

METHODS

A total of 786 individuals from 102 populations of six previously recognized species were collected for niche, morphological and genetic analyses. Three plastid DNA regions (matK, rbcL and trnH-psbA) and one nuclear DNA region [internal transcribed space (ITS)] were sequenced.

KEY RESULTS

Whereas six species had been previously recognized, statistical analyses based on character variation, molecular data and niche differentiation identified only two well-delimited clusters, together with a third possibly originating from relatively recent hybridization between, or historical introgression from, the other two.

CONCLUSIONS

Based on a principle of integrative species delimitation to reconcile different sources of data, the results provide compelling evidence that the six previously recognized species of the genus Orinus that were examined should be reduced to two, with new circumscriptions, and a third, identified in this study, should be described as a new species. This empirical study highlights the value of applying genetic differentiation, morphometric statistics and ecological niche modelling in an integrative approach to re-circumscribing species boundaries. The results produce relatively objective, operational and unbiased taxonomic classifications of plants occurring in remote regions.

Semi-permeable species boundaries in Iberian barbels (Barbus and Luciobarbus, Cyprinidae)

BACKGROUND

The evolution of species boundaries and the relative impact of selection and gene flow on genomic divergence are best studied in populations and species pairs exhibiting various levels of divergence along the speciation continuum. We studied species boundaries in Iberian barbels, Barbus and Luciobarbus, a system of populations and species spanning a wide degree of genetic relatedness, as well as geographic distribution and range overlap. We jointly analyze multiple types of molecular markers and morphological traits to gain a comprehensive perspective on the nature of species boundaries in these cyprinid fishes.

RESULTS

Intraspecific molecular and morphological differentiation is visible among many populations. Genomes of all sympatric species studied are porous to gene flow, even if they are not sister species. Compared to their allopatric counterparts, sympatric representatives of different species share alleles and show an increase in all measures of nucleotide polymorphism (S, Hd, K, π and θ). High molecular diversity is particularly striking in L. steindachneri from the Tejo and Guadiana rivers, which co-varies with other sympatric species. Interestingly, different nuclear markers introgress across species boundaries at various levels, with distinct impacts on population trees. As such, some loci exhibit limited introgression and population trees resemble the presumed species tree, while alleles at other loci introgress more freely and population trees reflect geographic affinities and interspecific gene flow. Additionally, extent of introgression decreases with increasing genetic divergence in hybridizing species pairs.

CONCLUSIONS

We show that reproductive isolation in Iberian Barbus and Luciobarbus is not complete and species boundaries are semi-permeable to (some) gene flow, as different species (including non-sister) are exchanging genes in areas of sympatry. Our results support a speciation-with-gene-flow scenario with heterogeneous barriers to gene flow across the genome, strengthening with genetic divergence. This is consistent with observations coming from other systems and supports the notion that speciation is not instantaneous but a gradual process, during which different species are still able to exchange some genes, while selection prevents gene flow at other loci. We also provide evidence for a hybrid origin of a barbel ecotype, L. steindachneri, suggesting that ecology plays a key role in species coexistence and hybridization in Iberian barbels. This ecotype with intermediate, yet variable, molecular, morphological, trophic and ecological characteristics is the local product of introgressive hybridization of L. comizo with up to three different species (with L. bocagei in the Tejo, with L. microcephalus and L. sclateri in the Guadiana). In spite of the homogenizing effects of ongoing gene flow, species can still be discriminated using a combination of morphological and molecular markers. Iberian barbels are thus an ideal system for the study of species boundaries, since they span a wide range of genetic divergences, with diverse ecologies and degrees of sympatry.

Phylogeny of anopheline (Diptera: Culicidae) species in southern Africa, based on nuclear and mitochondrial genes

A phylogeny of anthropophilic and zoophilic anopheline mosquito species was constructed, using the nuclear internal transcribed spacer 2 (ITS2) and mitochondrial cytochrome oxidase subunit I (COI) genes. The ITS2 alignment, typically difficult due to its noncoding nature and large size variations, was aided by using predicted secondary structure, making this phylogenetically useful gene more amenable to investigation. This phylogeny is unique in explicitly including zoophilic, non-vector anopheline species in order to illustrate their relationships to malaria vectors. Two new, cryptic species, Anopheles funestus-like and Anopheles rivulorum-like, were found to be present in Zambia for the first time. Sequences from the D3 region of the 28S rDNA suggest that the Zambian An. funestus-like may be a hybrid or geographical variant of An. funestus-like, previously reported in Malawi. This is the first report of An. rivulorum-like sympatric with An. rivulorum (Leeson), suggesting that these are separate species rather than geographic variants.

Genome-wide tests for introgression between cactophilic Drosophila implicate a role of inversions during speciation

Models of speciation‐with‐gene‐flow have shown that the reduction in recombination between alternative chromosome arrangements can facilitate the fixation of locally adaptive genes in the face of gene flow and contribute to speciation. However, it has proven frustratingly difficult to show empirically that inversions have reduced gene flow and arose during or shortly after the onset of species divergence rather than represent ancestral polymorphisms. Here, we present an analysis of whole genome data from a pair of cactophilic fruit flies, Drosophila mojavensis and D. arizonae, which are reproductively isolated in the wild and differ by several large inversions on three chromosomes. We found an increase in divergence at rearranged compared to colinear chromosomes. Using the density of divergent sites in short sequence blocks we fit a series of explicit models of species divergence in which gene flow is restricted to an initial period after divergence and may differ between colinear and rearranged parts of the genome. These analyses show that D. mojavensis and D. arizonae have experienced postdivergence gene flow that ceased around 270 KY ago and was significantly reduced in chromosomes with fixed inversions. Moreover, we show that these inversions most likely originated around the time of species divergence which is compatible with theoretical models that posit a role of inversions in speciation with gene flow.

Chromosome inversions, adaptive cassettes and the evolution of species' ranges

A chromosome inversion can spread when it captures locally adapted alleles or when it is introduced into a species by hybridization with adapted alleles that were previously absent. We present a model that shows how both processes can cause a species range to expand. Introgression of an inversion that carries novel, locally adapted alleles is a particularly powerful mechanism for range expansion. The model supports the earlier proposal that introgression of an inversion triggered a large range expansion of a malaria mosquito. These results suggest a role for inversions as cassettes of genes that can accelerate adaptation by crossing species boundaries, rather than protecting genomes from introgression.

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.

Biological invasions, climate change and genomics

The rate of biological invasions is expected to increase as the effects of climate change on biological communities become widespread. Climate change enhances habitat disturbance which facilitates the establishment of invasive species, which in turn provides opportunities for hybridization and introgression. These effects influence local biodiversity that can be tracked through genetic and genomic approaches. Metabarcoding and metagenomic approaches provide a way of monitoring some types of communities under climate change for the appearance of invasives. Introgression and hybridization can be followed by the analysis of entire genomes so that rapidly changing areas of the genome are identified and instances of genetic pollution monitored. Genomic markers enable accurate tracking of invasive species' geographic origin well beyond what was previously possible. New genomic tools are promoting fresh insights into classic questions about invading organisms under climate change, such as the role of genetic variation, local adaptation and climate pre-adaptation in successful invasions. These tools are providing managers with often more effective means to identify potential threats, improve surveillance and assess impacts on communities. We provide a framework for the application of genomic techniques within a management context and also indicate some important limitations in what can be achieved.

Hybrid sterility in crosses between two Brazilian sibling species of the Anopheles albitarsis complex

Background

Complexes of cryptic species are common in several taxa and this is also the case in the Anopheles genus, a group including all known human malaria vectors. The Anopheles albitarsis complex comprises at least nine cryptic species, some of which are implicated as vectors of human malaria. Several different types of data have been generated for this species complex such as cytogenetics, alloenzymes, morphological and feeding behavioral, hybridization experiments, RAPD-PCR and RFLP and mitochondrial and nuclear markers. Studies focused on its postzygotic isolation are still somewhat rare in the literature despite their importance to understand the speciation process and the level of gene flow potentially occurring among the different sibling species.

Methods

Hybridization experiments between Anopheles albitarsis s.s. and Anopheles marajoara, as well as backcrosses between hybrids and Anopheles albitarsis s.s., were performed using the induced mating technique. Results were compared to intraspecific crosses. Larva-to-adult viability and sex ratio were also assessed.

Results

Male hybrids show very low insemination rates and nearly complete sterility, apparently due to abnormalities in their reproductive organs. Evidence of partial sterility among the hybrid females was also observed.

Conclusions

Our data indicated that Anopheles albitarsis s.s. and Anopheles marajoara show a high level of postzygotic isolation with a strong hybrid male sterility. This result is consistent with the Haldane’s rule which states that in interspecific crosses the heterogametic sex is the first to be affected. However, the fact that the females are not completely sterile raises the possibility of introgression between these two siblings species.

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.

Adaptation through chromosomal inversions in Anopheles

Chromosomal inversions have been repeatedly involved in local adaptation in a large number of animals and plants. The ecological and behavioral plasticity of Anopheles species-human malaria vectors-is mirrored by high amounts of polymorphic inversions. The adaptive significance of chromosomal inversions has been consistently attested by strong and significant correlations between their frequencies and a number of phenotypic traits. Here, we provide an extensive literature review of the different adaptive traits associated with chromosomal inversions in the genus Anopheles. Traits having important consequences for the success of present and future vector control measures, such as insecticide resistance and behavioral changes, are discussed.

The speciation continuum: population structure, gene flow, and maternal ancestry in the Simulium arcticum complex (Diptera: Simuliidae)

Comparative analyses of populations at different stages of divergence can yield insights into the process of speciation. We assess population structure, gene flow, and maternal ancestry at five locations containing sympatric members of the Simulium arcticum complex at different stages of chromosome divergence. We analyze both nuclear and mitochondrial DNA markers, including 11 microsatellite loci, as well as COI, COII, cytb, and ND4 gene sequences. Simulium negativum, representing the later stages of divergence, shows both nuclear and mitochondrial differentiation when compared with allopatric and sympatric chromosomal forms, as well as both low contemporary and historical gene flow in sympatry. At intermediate stages of chromosome divergence, populations differ at nuclear, but not mitochondrial, loci in allopatry and sympatry. In one comparison of intermediate stage chromosomal forms (S. arcticum sensu stricto and S. apricarium), populations demonstrate low contemporary, but higher historical, gene flow in sympatry. In a second sympatric comparison (S. arcticum s. s. and S. brevicercum), both contemporary and historical gene flow are high. All analyses of sympatric populations at the earliest stages of chromosome divergence demonstrate panmixia; yet, some nuclear differentiation in allopatry is apparent. These findings suggest that molecular divergence is tracking chromosome divergence along a chromosomally-defined continuum of speciation in black flies.

Multigene phylogenetics reveals temporal diversification of major African malaria vectors

The major vectors of malaria in sub-Saharan Africa belong to subgenus Cellia. Yet, phylogenetic relationships and temporal diversification among African mosquito species have not been unambiguously determined. Knowledge about vector evolutionary history is crucial for correct interpretation of genetic changes identified through comparative genomics analyses. In this study, we estimated a molecular phylogeny using 49 gene sequences for the African malaria vectors An. gambiae, An. funestus, An. nili, the Asian malaria mosquito An. stephensi, and the outgroup species Culex quinquefasciatus and Aedes aegypti. To infer the phylogeny, we identified orthologous sequences uniformly distributed approximately every 5 Mb in the five chromosomal arms. The sequences were aligned and the phylogenetic trees were inferred using maximum likelihood and neighbor-joining methods. Bayesian molecular dating using a relaxed log normal model was used to infer divergence times. Trees from individual genes agreed with each other, placing An. nili as a basal clade that diversified from the studied malaria mosquito species 47.6 million years ago (mya). Other African malaria vectors originated more recently, and independently acquired traits related to vectorial capacity. The lineage leading to An. gambiae diverged 30.4 mya, while the African vector An. funestus and the Asian vector An. stephensi were the most closely related sister taxa that split 20.8 mya. These results were supported by consistently high bootstrap values in concatenated phylogenetic trees generated individually for each chromosomal arm. Genome-wide multigene phylogenetic analysis is a useful approach for discerning historic relationships among malaria vectors, providing a framework for the correct interpretation of genomic changes across species, and comprehending the evolutionary origins of this ubiquitous and deadly insect-borne disease.

Evidence to support natural hybridization between Anopheles sinensis and Anopheles kleini (Diptera: Culicidae): possibly a significant mechanism for gene introgression in sympatric populations

Background

Malaria caused by Plasmodium vivax is still a public health problem in the Republic of Korea (ROK), particularly regarding the recent re-emergence of this malarial species near the demilitarized zone in northwestern Paju City, Gyeonggi-do Province. Currently, at least 4 species (An. kleini, An. pullus, An. belenrae and An. lesteri) of the Hyrcanus Group are reported as possible natural vectors of vivax malaria in the ROK, and An. sinensis, which is the most dominant species, has long been incriminated as an important natural vector of this P. vivax. However, An. sinensis was ranked recently as a low potential vector. According to the discovery of natural hybrids between An. sinensis (a low potential vector for P. vivax) and An. kleini (a high potential vector for P. vivax) in Paju City, intensive investigation of this phenomenon is warranted under laboratory conditions.

Methods

Mosquitoes were collected during 2010-2012 from Paju City, ROK. Hybridization experiments used iso-female line colonies of these anophelines together with DNA analysis of ribosomal DNA [second internal transcribed spacer (ITS2)] and mitochondrial DNA [cytochrome c oxidase subunit I (COI)] of the parental colonies, F₁-hybrids and repeated backcross progenies were performed intensively by using a PCR-based assay and pyrosequencing technology.

Results

The results from hybridization experiments and molecular investigations revealed that the mitochondrial COI gene was introgressed from An. sinensis into An. kleini. The An. sinensis progenies obtained from consecutive repeated backcrosses in both directions, i.e., F₂-₁₁ progeny [(An. sinensis x An. kleini) x An. sinensis] and F₃-₅ progeny [(An. kleini x An. sinensis) x An. kleini] provided good supportive evidence.

Conclusions

This study revealed introgression of the mitochondrial COI gene between An. sinensis and An. kleini through consecutive repeated backcrosses under laboratory conditions. This new body of knowledge will be emphasized in reliable promising strategies in order to replace the population of An. kleini as a high potential vector for P. vivax, with that of a low potential vector, An. sinensis, through the mechanism of gene introgression in nature.