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

Genomics of urban adaptation and exaptation in mosquitoes and consequences for vectorial capacity

As urbanization accelerates around the world, mosquitoes that are capable of surviving and thriving in urban habitats increasingly spread mosquito-borne diseases. Across the >3,500 known species of mosquitoes, only a few rapidly adapted to the novel (on an evolutionary timescale) urban environments. In this review we highlight several emerging themes and testable hypotheses from recent literature. First, apparent urban adaptations can be roughly divided into newer adaptations arising in an urban context and exaptations - traits that evolved in a different context, prior to modern urbanization. Second, variants involved in urban adaptation are often partitioned among species complexes and cryptic lineages, and the history of gene flow-selection balance may be related to the evolution of compact genomic architectures that could facilitate rapid urban adaptation. Third, urban adaptation often has consequences for vectorial capacity - the ability of mosquitoes to serve as effective vectors of a particular pathogen - though the selective drivers and genetic mechanisms underlying these differences are incompletely understood. To fully understand urban adaptation in mosquitoes, we advocate for a coordinated effort to increase linkages between evolutionary ecology, population genomics, and medical entomology research. We discuss the two traits for which all three perspectives are the most developed - host preference and insecticide resistance - before reviewing several other less studied traits.

Unraveling the extensive phylogenetic discordance and evolutionary history of spurless taxa within the Aquilegia ecalcarata complex

Parallel evolution of the same, or at least very similar, phenotype(s) in different lineages is often interpreted as evidence for the action of natural selection. However, caution is required when inferring parallel evolution based on uncertain or potentially incorrect phylogenetic frameworks. Here, by conducting extensive phylogenomic and population genetic analyses, we aim to clarify the evolutionary history of spurless taxa within the Aquilegia ecalcarata complex. We observed substantial discordance in the phylogenetic patterns across the entire genome, primarily attributed to ancient introgression and incomplete lineage sorting. Additionally, we identified several spurless lineages whose phylogenetic positions were distorted by admixture events. Using a backbone tree and demographic modeling, we determined that these spurless taxa independently originated twice within this group. Intriguingly, our investigation revealed that the spurless taxa experienced population expansion during global cooling, while their spurred sister groups underwent population contraction. The parallel losses of petal spurs, therefore, may be linked to adaptations for low-temperature conditions. These findings emphasize the importance of comprehensive population-level analyses in phylogenetic inference and provide valuable insights into the dynamics of trait loss and its implications for the adaptive strategies.

Genomic Analysis Reveals a New Cryptic Taxon Within the Anopheles gambiae Complex With a Distinct Insecticide Resistance Profile in the Coast of East Africa

Anopheles mosquitoes are major malaria vectors, encompassing several species complexes with diverse life histories, transmission risks and insecticide resistance profiles that challenge malaria control efforts. This study investigated the genetic structure and insecticide resistance profiles of Anopheles gambiae complex mosquitoes in Tanzania. We analysed whole-genome sequence data of 300 mosquitoes collected between 2012 and 2015 across four regions in northern Tanzania and identified An. gambiae s.s., An. arabiensis and a distinct taxonomic group that was previously unknown. This distinct taxon has a unique profile of genetic diversity and appears restricted to the coastal region, and we refer to it as the Pwani molecular form. Analysis of insecticide resistance based on target-site mutations and copy number variations (CNV) showed that these markers were strikingly absent from the Pwani molecular form in contrast to other taxa. Our analysis also revealed a pattern of geographical isolation in the An. gambiae s.s. populations, with samples from the north-western site (Muleba) clustering separately from those collected in the north-eastern site (Muheza). These geographically isolated subpopulations also had differing resistance and selection profiles, with An. gambiae s.s. from the north-western site showing genomic evidence of higher resistance to pyrethroids compared with the north-eastern population. Conversely, An. arabiensis showed no geographical population structuring, with a similar insecticide resistance profile across all sampling locations, suggesting unrestricted gene flow. Our findings underscore the need to incorporate genetic data into malaria vector surveillance and control decisions and could inform the development and deployment of new interventions.

Reconciling Gene Tree Discordance and Biogeography in European Crows

Reconstructing the evolutionary history of young lineages diverging with gene flow is challenging due to factors like incomplete lineage sorting, introgression, and selection causing gene tree discordance. The European crow hybrid zone between all-black carrion crows and grey-coated hooded crows exemplifies this challenge. Most of the genome in Western and Central European carrion crow populations is near-identical to hooded crows, but differs substantially from their Iberian congeners. A notable exception is a single major-effect colour-locus under sexual selection aligning with the 'species' tree. To understand the underlying evolutionary processes, we reconstructed the biogeographic history of the species complex. During the Pleistocene carrion and hooded crows took refuge in the Iberian Peninsula and the Middle East, respectively. Allele-sharing of all-black Western European populations with likewise black Iberian crows at the colour-locus represents the last trace of carrion crow ancestry, resisting gene flow from expanding hooded crow populations that have homogenised most of the genome. A model of colour-locus introgression from an Iberian ancestor into hooded crow populations near the Pyrenées was significantly less supported. We found no positive relationship between introgression and recombination rate consistent with the absence of genome-wide, polygenic barriers in this young species complex. Overall, this study portrays a scenario where few large-effect loci, subject to divergent sexual selection, resist rampant and asymmetric gene exchange. This study underscores the importance of integrating population demography and biogeography to accurately interpret patterns of gene tree discordance following population divergence.

Sweeps in space: leveraging geographic data to identify beneficial alleles in Anopheles gambiae

As organisms adapt to environmental changes, natural selection modifies the frequency of non-neutral alleles. For beneficial mutations, the outcome of this process may be a selective sweep, in which an allele rapidly increases in frequency and perhaps reaches fixation within a population. Selective sweeps have well-studied effects on patterns of local genetic variation in panmictic populations, but much less is known about the dynamics of sweeps in continuous space. In particular, because limited movement across a landscape leads to unique patterns of population structure, spatial dynamics may influence the trajectory of selected mutations. Here, we use forward-in-time, individual-based simulations in continuous space to study the impact of space on beneficial mutations as they sweep through a population. In particular, we show that selection changes the joint distribution of allele frequency and geographic range occupied by a focal allele and demonstrate that this signal can be used to identify selective sweeps. We then leverage this signal to identify in-progress selective sweeps within the malaria vector Anopheles gambiae , a species under strong selection pressure from vector control measures. By considering space, we identify multiple previously undescribed variants with potential phenotypic consequences, including mutations impacting known IR-associated genes and altering protein structure and properties. Our results demonstrate a novel signal for detecting selection in spatial population genetic data that may have implications for genomic surveillance and understanding geographic patterns of genetic variation.

Phylogenetic taxonomy of the Zambian Anopheles coustani group using a mitogenomics approach

Background

Mosquito species belonging to the Anopheles coustani group have been implicated in driving residual malaria transmission in sub-Saharan Africa and are regarded as an established primary vector in Madagascar. The morphological identification of mosquitoes in this group is challenging due to cryptic features and their molecular confirmation is difficult due to a paucity of reference sequence data representing all members of the group. Conventional molecular barcoding with the cytochrome oxidase I (COI) gene and the internal transcribed spacer 2 (ITS2) region targets is limited in their discrimination and conclusive identification of members of species complexes. In contrast, complete mitochondrial genomes (mitogenomes) have demonstrated much improved power over barcodes to be useful in rectifying taxonomic discrepancies in Culicidae.

Methods

We utilized a genome skimming approach via shallow shotgun sequencing on individual mosquito specimens to generate sequence reads for mitogenome assembly. Bayesian inferred phylogenies and molecular dating estimations were perfomed on the concatenated protein coding genes using the Bayesian Evolutionary Analysis by Sampling Trees 2 (BEAST 2) platform. Divergence estimates were calibrated on published calucations for Anopheles-Aedes.

Results

This study generated 17 new complete mitogenomes which were comprable to reference An. coustani mitogenomes in the GenBank repository by having 13 protein coding, 22 transfer RNA and 2 ribosomal RNA genes, with an average length of 15,400 bp and AT content of 78.3%. Bayesian inference using the concatenated protein coding genes from the generated and publicly available mitogenomes yielded six clades: one for each of the four taxa targeted in this study, the GenBank references, and a currently unknown species. Divergence times estimated that the An. coustani group separated from the An. gambiae complex approximately 110 million years ago (MYA), and members within the complex diverged at times points ranging from~34 MYA to as recent as ~7 MYA.

Conclusions

These findings demonstrate the value of mitochondrial genomes in differentiating cryptic taxa and help to confirm morphological identities of An. coustani s.s., An. paludis, An. zeimanni and An. tenebrosus. Divergence estimates with the An. coustani group are similar to those for well-studied anopheline vector groups. These analyses also highlight the likely prescence of other cryptic An. coustani group members circulating in Zambia.

Species Delimitation Using Genomic Data: Options and Limitations

The most effective approaches for species discovery and species validation with genomic data remain underexplored. This study evaluates the merits and limitations of phylogenetic approaches based on the multispecies coalescent model and population genetic approaches for species discovery, i.e., species delimitation in the absence of prior knowledge, using genomic datasets from four well-known radiations. Furthermore, it demonstrates how geographic data can be integrated with the genomic data for species validation, i.e., for testing primary species hypotheses. The multispecies coalescent model-based approaches tr2 and soda resulted in high over-splitting of species, low percentages of species delimited according to the current classification, and low percentages of individuals assigned to the same species as in the current classification across all four species complexes studied. Conversely, the species numbers were slightly underestimated based on the structure results. Although the proportion of species delimited according to the current classification and the proportion of individuals assigned to the same species as in the current classification in the classifications based on the structure results is approximately twice that of the classifications proposed by the multispecies coalescent model-based approaches, it remains unsatisfactory. A slight over-splitting of species into population groups can be corrected by species validation with isolation-by-distance tests if a sufficient number of populations have been sampled for each species. Sampling design is an essential step in any taxonomic study, as it has a significant impact on the delimitation of the species and the possibility of their validation.

Molecular evolution of gustatory receptors in the Anopheles gambiae complex

BACKGROUND

Mosquitoes in the Anopheles (An.) gambiae species complex are major vectors of Plasmodium falciparum malaria. One reason for this is the high anthropophily of the constituent species An. coluzzii, An. gambiae sensu stricto, and An. arabiensis. In contrast, their sister species An. quadriannulatus is highly zoophilic. Anopheles mosquitoes largely rely on chemical cues for host-seeking, which are primarily detected by four chemosensory gene families: olfactory receptors (Ors), ionotropic receptors (Irs), gustatory receptors (Grs), and odorant binding proteins (Obps). Genes from these families that have been implicated in host adaptation show evidence of positive selection in other insect species, including other mosquitoes. As such, we analyzed the molecular evolutionary patterns of the gustatory receptors within the Anopheles gambiae complex, with a particular interest in identifying Grs that show evidence of positive selection in highly anthropophilic species.

RESULTS

We identified sixteen Grs that show evidence of potential positive selection using the McDonald-Kreitman test, including four putative sugar receptors and two Grs with unknown ligands that are relatively highly expressed in chemosensory organs of either An. coluzzii or An. quadriannulatus. In addition, we identified twelve Grs that show evidence of potential purifying selection using the McDonald-Kreitman test, and twelve Grs that may have experienced a selective sweep using the DH test, including three putative sugar receptors and the carbon dioxide receptor Gr24. We also identified both positive and purifying selection in the coastal species An. melas (West Africa) and An. merus (East Africa).

CONCLUSIONS

Our results, together with transcriptomic data, identify four Grs as possible candidates for involvement in the evolution of vertebrate host preference in the An. gambiae complex, as may have occurred in the An. farauti complex. They also point to sugar receptors as playing a role in recent adaptation of some of these species. As the vast majority of Grs have unknown functions and much is still unknown about the role of Grs in these species, a more complete interpretation of our data necessitates further characterization of these genes.

Population genomics of Anopheles darlingi, the principal South American malaria vector mosquito

Malaria in South America remains a serious public health problem. Anopheles (Nyssorhynchus) darlingi is the most important malaria vector across tropical Latin America. Vector-targeted disease control efforts require a thorough understanding of mosquito demographic and evolutionary patterns. We present and analyze whole genomes of 1094 A. darlingi (median depth 18x) from six South American countries. We observe deep geographic population structure, high genetic diversity including thirteen putative segregating inversions, and no evidence for cryptic sympatric taxa despite high interpopulation divergence. Strong signals of selection are plausibly driven by insecticides, especially on cytochrome P450 genes, one of which we validated experimentally. Our results will facilitate effective mosquito surveillance and control, while highlighting ongoing challenges that a diverse vector poses for malaria elimination in the western hemisphere.

Comparative Chloroplast Genomics Reveals Intrageneric Divergence in Salix

As the most diverse genus of Salicaceae, Salix is primarily distributed in the temperate zone of the Northern Hemisphere, encompassing 350-500 species worldwide. The genus's evolutionary history is complex due to significant genetic differentiation. Chloroplast genes, being highly conserved, serve as effective tools for studying uniparental inheritance and evolution. In this study, we sequenced and assembled the chloroplast genomes of five representative Salix species. Phylogenetic relationships were constructed using chloroplast genome data, and structural differences among lineages were compared. These Salix chloroplast genomes exhibited a typical quadripartite structure, with lengths ranging from 154,444 to 155,725 bp. We successfully annotated 131 genes, including 88 protein-coding genes, 35 tRNA genes, and 8 rRNA genes. Clade I showed higher variability in the SSC region, identifying five highly variable regions: petA-psbJ, rps16-rps3, ndhD, ccsA-ndhD, and ndhG-ndhI. Two rapidly evolving genes, ndhI and ycf4, were also identified. The total length of insertions and deletions (InDels) in Clade I was 1046 bp. Clade II exhibited greater variability in the LSC region, with four highly variable regions being identified: trnK-trnQ, ndhC-trnV, trnV, and psdE-petL. Four rapidly evolving genes-infA, rpoC1, rps18, and ycf1-were identified. The total length of InDels in Clade II was 1282 bp. Therefore, this study elucidated the chloroplast genome evolution across different Salix lineages, thereby providing deeper insights into intrageneric phylogenetic relationships.

Adaptive evolution of traits for parasitism and pathogen transmission potential in bat flies

Deciphering the mechanisms underlying the transmission and spillover of zoonoses from reservoir hosts is essential in preventing future global pandemics. Bat flies-obligate blood-feeding ectoparasites of bats-are known carriers of diverse viruses. Here, we conducted a de novo assembly of a chromosome-level genome for the bat fly species Phthiridium sp. Comparative genomic analysis unveiled genes associated with specialized traits, such as the loss of eyes and wings, as well as elongated legs, which have adapted to parasitism on the dense fur of bats. Utilizing small RNA sequencing, we identified a spectrum of known and previously unclassified viruses in bat flies. Notably, experimental evidence indicated that bat flies can also feed on mammalian hosts other than bats, suggesting the potential for the spillover of bat-borne viruses. Furthermore, we demonstrated the role of the bat fly's RNA interference pathway in influencing the diversity and evolution of viruses. In summary, this study not only presents a new genome catalog to unveil the evolutionary mechanisms underpinning bat fly parasitism, but also provides a novel research system that can be used to investigate the mechanisms of cross-species transmission of bat-borne viruses and the co-evolution of bats and viruses.

Phylogenomic Analysis Reveals Evolutionary Relationships of Tropical Drosophilidae: From Drosophila to Scaptodrosophila

Species radiation in the family Drosophilidae has led to a diversity of species occupying a wide range of ecological niches. Despite the high diversity within this family, with over thousands of species and more than a hundred species recorded in Thailand and the Malay Peninsula, taxonomic classifications remain complicated due to morphological plasticity and inconsistent phylogenetic reconstructions based on limited genetic data. In this study, we assemble new mitochondrial genomes from Drosophila and Scaptodrosophila species collected in Thailand, expanding the genomic resources for these underexplored tropical regions. Phylogenetic analyses of 13 mitochondrial protein-coding genes revealed well-supported evolutionary relationships, with Scaptodrosophila forming a distinct lineage and several Drosophila subgroups, such as ananassae, montium, and melanogaster, exhibiting monophyly. Notable discrepancies were observed in the placement of the suzukii subgroup, which was not recovered as monophyletic, and the position of the punjabiensis subgroup, reflecting the complexities of lineage sorting and hybridization events. This comprehensive genomic analysis provides a more accurate understanding of evolutionary relationships within the family Drosophilidae's diversification in the tropics.

The hybrid history of zebrafish

Since the description of zebrafish (Danio rerio) in 1822, the identity of its closest living relative has been unclear. To address this problem, we sequenced the exomes of 10 species in genus Danio, using the closely related Devario aequipinnatus as outgroup, to infer relationships across the 25 chromosomes of the zebrafish genome. The majority of relationships within Danio were remarkably consistent across all chromosomes. Relationships of chromosome segments, however, depended systematically upon their genomic location within zebrafish chromosomes. Regions near chromosome centers identified Danio kyathit and/or Danio aesculapii as the closest relative of zebrafish, while segments near chromosome ends supported only D. aesculapii as the zebrafish sister species. Genome-wide comparisons of derived character states revealed that danio relationships are inconsistent with a simple bifurcating species history but support an ancient hybrid origin of the D. rerio lineage by homoploid hybrid speciation. We also found evidence of more recent gene flow limited to the high recombination ends of chromosomes and several megabases of chromosome 20 with a history distinct from the rest of the genome. Additional insights gained from incorporating genome structure into a phylogenomic study demonstrate the utility of such an approach for future studies in other taxa. The multiple genomic histories of species in the genus Danio have important implications for comparative studies in these morphologically varied and beautiful species and for our understanding of the hybrid evolutionary history of zebrafish.

The Genomic Landscape, Causes, and Consequences of Extensive Phylogenomic Discordance in Murine Rodents

A species tree is a central concept in evolutionary biology whereby a single branching phylogeny reflects relationships among species. However, the phylogenies of different genomic regions often differ from the species tree. Although tree discordance is widespread in phylogenomic studies, we still lack a clear understanding of how variation in phylogenetic patterns is shaped by genome biology or the extent to which discordance may compromise comparative studies. We characterized patterns of phylogenomic discordance across the murine rodents-a large and ecologically diverse group that gave rise to the laboratory mouse and rat model systems. Combining recently published linked-read genome assemblies for seven murine species with other available rodent genomes, we first used ultraconserved elements (UCEs) to infer a robust time-calibrated species tree. We then used whole genomes to examine finer-scale patterns of discordance across ∼12 million years of divergence. We found that proximate chromosomal regions tended to have more similar phylogenetic histories. There was no clear relationship between local tree similarity and recombination rates in house mice, but we did observe a correlation between recombination rates and average similarity to the species tree. We also detected a strong influence of linked selection whereby purifying selection at UCEs led to appreciably less discordance. Finally, we show that assuming a single species tree can result in substantial deviation from the results with gene trees when testing for positive selection under different models. Collectively, our results highlight the complex relationship between phylogenetic inference and genome biology and underscore how failure to account for this complexity can mislead comparative genomic studies.

Thrushes in Love: Extensive Gene Flow, With Differential Resistance and Selection, Obscures and Reveals the Evolutionary History of a Songbird Clade

The application of high-throughput sequencing to phylogenetic analyses is allowing authors to reconstruct the true evolutionary history of species. This work can illuminate specific mechanisms underlying divergence when combined with analyses of gene flow, recombination and selection. We conducted a phylogenomic analysis of Catharus, a songbird genus with considerable potential for gene flow, variation in migratory behaviour and genomic resources. We documented discordance among trees constructed for mitochondrial, autosomal and sex (Z) chromosome partitions. Two trees were recovered on the Z. Both trees differed from the autosomes, one matched the mitochondria, and the other was unique to the Z. Gene flow with one species likely generated much of this discordance; substantial admixture between ustulatus and the remaining species was documented and linked to at least two historic events. The tree unique to the Z likely reflects the true history of Catharus; local genomic analyses recovered the same tree in autosomal regions with reduced admixture and recombination. Genes previously connected to migration were enriched in these regions suggesting transitions between migratory and non-migratory states helped generate divergence. Migratory (vs. nonmigratory) Catharus formed a monophyletic clade in a subset of genomic regions. Gene flow was elevated in some of these regions suggesting adaptive introgression may have occurred, but the dominant pattern was of balancing selection maintaining ancestral polymorphisms important for olfaction and perhaps, by extension, adaptation to temperate climates. This work illuminates the evolutionary history of an important model in speciation and demonstrates how differential resistance to gene flow can affect local genomic patterns.

Headaches and adolescents: why so many failures in their management

The management of headaches in children and adolescents is still a challenge, with patients experiencing pain for years and polypharmacy. We reviewed the medical history of 31 patients referred to our pain clinic for chronic headaches between April 2023 and April 2024. There were more female than male patients (73%). Patients have been reporting headaches for 52 ± 44 months on average. Twenty-nine patients (94%) were experiencing different types of pain besides headaches. The most common medication prescribed in this group of patients was topiramate (70%), followed by tricyclic antidepressants (35%) and triptans (21%). Patients had been prescribed and tried, on average, 4.5 ± 2 (range 1-10) different medications to manage headaches and concomitant psychiatric disorders. Twenty-two patients (71%) had been diagnosed with a psychiatric disorder, including depression, anxiety, and PTSD, and 16% had a history of attempted suicide/self-harm. Fourteen of them (45%) had been prescribed antidepressants or benzodiazepines. They had been prescribed, on average, 2 ± 1 (range 1-4) psychiatric stabilizer medications. We agreed with the referral diagnosis in 39% of the patients. We attributed the headaches to more complex chronic pain conditions, including fibromyalgia (15%) and AMPS (15%), autism with sensory integration problems (9%), and major depression (9%). Patients had seen an average of 3 ± 1 (range 1-5) different specialists; none consulted a pain specialist. Patients underwent between 0 (12%) and four tests (6%), including MRI (52%) and CT of the brain (8%). These neuroimaging studies did not demonstrate any brain pathology. We prescribed new medications and treatments, including nerve blocks, in 19 (61%) patients. In 47% of the cases, patients reported improved headaches, while 22% did not feel our recommendations were effective. Twenty-one percent of patients never came back to the clinic for a follow-up. Significant catastrophizing was present in 57% of the patients; 52% of patients had mild to severe anxiety, and 57% had symptoms of depressive disorder.

CONCLUSIONS

Headaches are often the manifestation of more complex pain syndromes that require a more holistic approach, different from conventional pharmacological management.

WHAT IS KNOWN

• Headache is one of the most disabling diseases. • Prevalence of headaches in hildrens and adolescents can be as high as 58%.

WHAT IS NEW

• Conventional pharmacological management often fails to help young patients. • A relationship between chronic headaches and psychopathology should be investigated in these young patients.

Contrast and Genomic Characterisation of Ancient and Recent Interspecific Introgression Between Deeply Diverged Moustache Toads (Leptobrachium)

Recent genomic analyses have provided new insights into the process of interspecific introgression and its consequences on species evolution. Most recent studies, however, focused on hybridization between recently radiated species, with few examining the genomic outcomes of ancient hybridization across deeply diverged species. Using whole genome data of moustache toads (Leptobrachium), we identified signals of three hybridization events among nine species that diverged at the Eocene. An ancient introgression from L. leishanense to the ancestral branch (C1) of L. liui introduced adaptive variants. The highly introgressed regions include genes with important functions in odorant detection and immune responses. These genes are preserved in all three descendent populations of L. liui_C1, and these regions likely have been positively selected over a long filtering process. A recent introgression occurred from L. huashen to L. tengchongense, with the introgressed regions being mostly neutral. Furthermore, one F1 hybrid individual was detected between sympatric L. ailaonicum and L. promustache. The signals of introgression largely disappeared after removing the hybrid individual, indicating an occasional hybridization but minimal introgression. Further examination of highly divergent but low introgressed genomic regions revealed both pre-mating isolation and genetic incompatibility as potential mechanisms of resisting introgression and maintaining species boundaries. Additionally, no large X-effect was found in these introgression events. Hybridization between deeply diverged amphibian species may be common, but detectable introgressions are likely less so, with recent introgression being mostly neutral and the rare ancient one potentially adaptive. Our findings complement recent genomic work, and together they provide a better understanding of the genomic characteristics of interspecific introgression and its significance in species adaptation and evolution.

Hierarchical Heuristic Species Delimitation Under the Multispecies Coalescent Model with Migration

The multispecies coalescent (MSC) model accommodates genealogical fluctuations across the genome and provides a natural framework for comparative analysis of genomic sequence data from closely related species to infer the history of species divergence and gene flow. Given a set of populations, hypotheses of species delimitation (and species phylogeny) may be formulated as instances of MSC models (e.g., MSC for 1 species versus MSC for 2 species) and compared using Bayesian model selection. This approach, implemented in the program bpp, has been found to be prone to over-splitting. Alternatively, heuristic criteria based on population parameters (such as population split times, population sizes, and migration rates) estimated from genomic data may be used to delimit species. Here, we develop hierarchical merge and split algorithms for heuristic species delimitation based on the genealogical divergence index (gdi) and implement them in a Python pipeline called hhsd. We characterize the behavior of the gdi under a few simple scenarios of gene flow. We apply the new approaches to a dataset simulated under a model of isolation by distance as well as 3 empirical datasets. Our tests suggest that the new approaches produced sensible results and were less prone to oversplitting. We discuss possible strategies for accommodating paraphyletic species in the hierarchical algorithm, as well as the challenges of species delimitation based on heuristic criteria.

Unravelling the role of mitochondrial DNA in hybrid incompatibility within species of the Anopheles gambiae complex

Isolation mechanisms between mosquito species of the Anopheles gambiae complex, which includes major malaria vectors, remain poorly understood. In some cases, pre-zygotic barriers have been shown to limit gene flow between species of the complex, leading to a low level of hybridisation in nature. Post-zygotic mechanisms manifest with F1 hybrid males fully sterile and F1 hybrid females with reduced fertility. Genetic approaches combined with DNA sequencing techniques have highlighted the involvement of genomic regions in hybrid incompatibility with a predominant role of the X chromosome. In addition, differences in the phenotype of F1 hybrid males have been identified depending on the directionality of the parental cross used to generate them. All these studies have focused on the interaction of nuclear DNA elements in hybrid individuals. Given the role that mitochondrial DNA plays in genetic incompatibilities within other organisms and its unique inheritance pattern, commonly maternal, we conducted a genetic study that relied on the introgression of mitochondrial DNA between Anopheles gambiae and Anopheles arabiensis. The findings indicate that the mitochondrial switch does not appear to restore the fertility of F1 hybrid males, suggesting that mitochondrial DNA may not play a role in hybrid incompatibilities in these Anopheles species.

Cas9/guide RNA-based gene-drive dynamics following introduction and introgression into diverse anopheline mosquito genetic backgrounds

BACKGROUND

Novel technologies are needed to combat anopheline vectors of malaria parasites as the reductions in worldwide disease incidence has stalled in recent years. Gene drive-based approaches utilizing Cas9/guide RNA (gRNA) systems are being developed to suppress anopheline populations or modify them by increasing their refractoriness to the parasites. These systems rely on the successful cleavage of a chromosomal DNA target site followed by homology-directed repair (HDR) in germline cells to bias inheritance of the drive system. An optimal drive system should be highly efficient for HDR-mediated gene conversion with minimal error rates. A gene-drive system, AgNosCd-1, with these attributes has been developed in the Anopheles gambiae G3 strain and serves as a framework for further development of population modification strains. To validate AgNosCd-1 as a versatile platform, it must perform well in a variety of genetic backgrounds.

RESULTS

We introduced or introgressed AgNosCd-1 into different genetic backgrounds, three in geographically-diverse Anopheles gambiae strains, and one each in an An. coluzzii and An. arabiensis strain. The overall drive inheritance, determined by presence of a dominant marker gene in the F2 hybrids, far exceeded Mendelian inheritance ratios in all genetic backgrounds that produced viable progeny. Haldane's rule was confirmed for AgNosCd-1 introgression into the An. arabiensis Dongola strain and sterility of the F1 hybrid males prevented production of F2 hybrid offspring. Back-crosses of F1 hybrid females were not performed to keep the experimental design consistent across all the genetic backgrounds and to avoid maternally-generated mutant alleles that might confound the drive dynamics. DNA sequencing of the target site in F1 and F2 mosquitoes with exceptional phenotypes revealed drive system-generated mutations resulting from non-homologous end joining events (NHEJ), which formed at rates similar to AgNosCd-1 in the G3 genetic background and were generated via the same maternal-effect mechanism.

CONCLUSIONS

These findings support the conclusion that the AgNosCd-1 drive system is robust and has high drive inheritance and gene conversion efficiency accompanied by low NHEJ mutation rates in diverse An. gambiae s.l. laboratory strains.

The Meaning and Measure of Concordance Factors in Phylogenomics

As phylogenomic datasets have grown in size, researchers have developed new ways to measure biological variation and to assess statistical support for specific branches. Larger datasets have more sites and loci and therefore less sampling variance. While we can more accurately measure the mean signal in these datasets, lower sampling variance is often reflected in uniformly high measures of branch support-such as the bootstrap and posterior probability-limiting their utility. Larger datasets have also revealed substantial biological variation in the topologies found across individual loci, such that the single species tree inferred by most phylogenetic methods represents a limited summary of the data for many purposes. In contrast to measures of statistical support, the degree of underlying topological variation among loci should be approximately constant regardless of the size of the dataset. "Concordance factors" (CFs) and similar statistics have therefore become increasingly important tools in phylogenetics. In this review, we explain why CFs should be thought of as descriptors of topological variation rather than as measures of statistical support, and argue that they provide important information about the predictive power of the species tree not contained in measures of support. We review a growing suite of statistics for measuring concordance, compare them in a common framework that reveals their interrelationships, and demonstrate how to calculate them using an example from birds. We also discuss how measures of topological variation might change in the future as we move beyond estimating a single "tree of life" toward estimating the myriad evolutionary histories underlying genomic variation.

Detection and quantification of introgression using Bayesian inference based on conjugate priors

SUMMARY

Introgression (the flow of genes between species) is a major force structuring the evolution of genomes, potentially providing raw material for adaptation. Here, we present a versatile Bayesian model selection approach for detecting and quantifying introgression, df-BF, that builds upon the recently published distance-based df statistic. Unlike df, df-BF accounts for the number of variant sites within a genomic region. The underlying model parameter of our df-BF method, here denoted as dfθ, accurately quantifies introgression, and the corresponding Bayes Factors (df-BF) enables weighing the strength of evidence for introgression. To ensure fast computation, we use conjugate priors with no need for computationally demanding MCMC iterations. We compare our method with other approaches including df, fd, Dp, and Patterson's D using a wide range of coalescent simulations. Furthermore, we showcase the applicability of df-BF and dfθ using whole-genome mosquito data. Finally, we integrate the new method into the powerful genomics R-package PopGenome.

AVAILABILITY AND IMPLEMENTATION

The presented methods are implemented within the R-package PopGenome (https://github.com/pievos101/PopGenome) and the simulation as the application results can be reproduced from the source code available from a dedicated GitHub repository (https://github.com/pievos101/Introgression-Simulation).

Towards Reliable Detection of Introgression in the Presence of Among-Species Rate Variation

The role of interspecific hybridization has recently seen increasing attention, especially in the context of diversification dynamics. Genomic research has now made it abundantly clear that both hybridization and introgression-the exchange of genetic material through hybridization and backcrossing-are far more common than previously thought. Besides cases of ongoing or recent genetic exchange between taxa, an increasing number of studies report "ancient introgression"- referring to results of hybridization that took place in the distant past. However, it is not clear whether commonly used methods for the detection of introgression are applicable to such old systems, given that most of these methods were originally developed for analyses at the level of populations and recently diverged species, affected by recent or ongoing genetic exchange. In particular, the assumption of constant evolutionary rates, which is implicit in many commonly used approaches, is more likely to be violated as evolutionary divergence increases. To test the limitations of introgression detection methods when being applied to old systems, we simulated thousands of genomic datasets under a wide range of settings, with varying degrees of among-species rate variation and introgression. Using these simulated datasets, we showed that some commonly applied statistical methods, including the D-statistic and certain tests based on sets of local phylogenetic trees, can produce false-positive signals of introgression between divergent taxa that have different rates of evolution. These misleading signals are caused by the presence of homoplasies occurring at different rates in different lineages. To distinguish between the patterns caused by rate variation and genuine introgression, we developed a new test that is based on the expected clustering of introgressed sites along the genome and implemented this test in the program Dsuite.

Common misconceptions of speciation

Speciation is a complex process that can unfold in many different ways. Speciation researchers sometimes simplify core principles in their writing in a way that implies misconceptions about the speciation process. While we think that these misconceptions are usually inadvertently implied (and not actively believed) by the researchers, they nonetheless risk warping how external readers understand speciation. Here we highlight six misconceptions of speciation that are especially widespread. First, species are implied to be clearly and consistently defined entities in nature, whereas in reality species boundaries are often fuzzy and semipermeable. Second, speciation is often implied to be 'good', which is two-fold problematic because it implies both that evolution has a goal and that speciation universally increases the chances of lineage persistence. Third, species-poor clades with species-rich sister clades are considered 'primitive' or 'basal', falsely implying a ladder of progress. Fourth, the evolution of species is assumed to be strictly tree-like, but genomic findings show widespread hybridization more consistent with network-like evolution. Fifth, a lack of association between a trait and elevated speciation rates in macroevolutionary studies is often interpreted as evidence against its relevance in speciation-even if microevolutionary case studies show that it is relevant. Sixth, obvious trait differences between species are sometimes too readily assumed to be (i) barriers to reproduction, (ii) a stepping-stone to inevitable speciation, or (iii) reflective of the species' whole divergence history. In conclusion, we call for caution, particularly when communicating science, because miscommunication of these ideas provides fertile ground for misconceptions to spread.

Whole-genome sequencing of major malaria vectors reveals the evolution of new insecticide resistance variants in a longitudinal study in Burkina Faso

BACKGROUND

Intensive deployment of insecticide based malaria vector control tools resulted in the rapid evolution of phenotypes resistant to these chemicals. Understanding this process at the genomic level is important for the deployment of successful vector control interventions. Therefore, longitudinal sampling followed by whole genome sequencing (WGS) is necessary to understand how these evolutionary processes evolve over time. This study investigated the change in genetic structure and the evolution of the insecticide resistance variants in natural populations of Anopheles gambiae over time and space from 2012 to 2017 in Burkina Faso.

METHODS

New genomic data have been generated from An. gambiae mosquitoes collected from three villages in the western part of Burkina Faso between 2012 and 2017. The samples were whole-genome sequenced and the data used in the An. gambiae 1000 genomes (Ag1000G) project as part of the Vector Observatory. Genomic data were analysed using the analysis pipeline previously designed by the Ag1000G project.

RESULTS

The results showed similar and consistent nucleotide diversity and negative Tajima's D between An. gambiae sensu stricto (s.s.) and Anopheles coluzzii. Principal component analysis (PCA) and the fixation index (FST) showed a clear genetic structure in the An. gambiae sensu lato (s.l.) species. Genome-wide FST and H12 scans identified genomic regions under divergent selection that may have implications in the adaptation to ecological changes. Novel voltage-gated sodium channel pyrethroid resistance target-site alleles (V402L, I1527T) were identified at increasing frequencies alongside the established alleles (Vgsc-L995F, Vgsc-L995S and N1570Y) within the An. gambiae s.l.

POPULATIONS

Organophosphate metabolic resistance markers were also identified, at increasing frequencies, within the An. gambiae s.s. populations from 2012 to 2017, including the SNP Ace1-G280S and its associated duplication. Variants simultaneously identified in the same vector populations raise concerns about the long-term efficacy of new generation bed nets and the recently organophosphate pirimiphos-methyl indoor residual spraying in Burkina Faso.

CONCLUSION

These findings highlighted the benefit of genomic surveillance of malaria vectors for the detection of new insecticide resistance variants, the monitoring of the existing resistance variants, and also to get insights into the evolutionary processes driving insecticide resistance.

Population genomic evidence of a putative 'far-west' African cryptic taxon in the Anopheles gambiae complex

The two main Afrotropical malaria vectors - Anopheles coluzzii and An. gambiae - are genetically distinct and reproductively isolated across West Africa. However, populations at the western extreme of their range are assigned as "intermediate" between the two species by whole genome sequence (WGS) data, and as hybrid forms by conventional molecular diagnostics. By exploiting WGS data from 1190 specimens collected across west Africa via the Anopheles gambiae 1000 Genomes network, we identified a putative taxon in the far-west (provisionally named Bissau molecular form), which did not arise by admixture but rather may have originated at the same time as the split between An. coluzzii and An. gambiae. Intriguingly, this taxon lacks insecticide resistance mechanisms commonly observed in the two main species. These findings lead to a change of perspective on malaria vector species in the far-west region with potential for epidemiological implications, and a new challenge for genetic-based mosquito control approaches.

Museum Skins Enable Identification of Introgression Associated with Cytonuclear Discordance

Increased sampling of genomes and populations across closely related species has revealed that levels of genetic exchange during and after speciation are higher than previously thought. One obvious manifestation of such exchange is strong cytonuclear discordance, where the divergence in mitochondrial DNA (mtDNA) differs from that for nuclear genes more (or less) than expected from differences between mtDNA and nuclear DNA (nDNA) in population size and mutation rate. Given genome-scale data sets and coalescent modeling, we can now confidently identify cases of strong discordance and test specifically for historical or recent introgression as the cause. Using population sampling, combining exon capture data from historical museum specimens and recently collected tissues we showcase how genomic tools can resolve complex evolutionary histories in the brachyotis group of rock-wallabies (Petrogale). In particular, applying population and phylogenomic approaches we can assess the role of demographic processes in driving complex evolutionary patterns and assess a role of ancient introgression and hybridization. We find that described species are well supported as monophyletic taxa for nDNA genes, but not for mtDNA, with cytonuclear discordance involving at least 4 operational taxonomic units across 4 species which diverged 183-278 kya. ABC modeling of nDNA gene trees supports introgression during or after speciation for some taxon pairs with cytonuclear discordance. Given substantial differences in body size between the species involved, this evidence for gene flow is surprising. Heterogenous patterns of introgression were identified but do not appear to be associated with chromosome differences between species. These and previous results suggest that dynamic past climates across the monsoonal tropics could have promoted reticulation among related species.

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

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

Complex patterns of gene flow and convergence in the evolutionary history of the spiral-horned antelopes (Tragelaphini)

The Tragelaphini, also known as spiral-horned antelope, is a phenotypically diverse mammalian tribe comprising a single genus, Tragelaphus. The evolutionary history of this tribe has attracted the attention of taxonomists and molecular geneticists for decades because its diversity is characterised by conflicts between morphological and molecular data as well as between mitochondrial, nuclear and chromosomal DNA. These inconsistencies point to a complex history of ecological diversification, coupled by either phenotypic convergence or introgression. Therefore, to unravel the phylogenetic relationships among spiral-horned antelopes, and to further investigate the role of divergence and gene flow in trait evolution, we sequenced genomes for all nine accepted species of the genus Tragelaphus, including a genome each for the highly divergent bushbuck lineages (T. s. scriptus and T. s. sylvaticus). We successfully reconstructed the Tragelaphus species tree, providing genome-level support for the early Pliocene divergence and monophyly of the nyala (T. angasii) and lesser kudu (T. imberbis), the monophyly of the two eland species (T. oryx and T. derbianus) and, importantly, the monophyly of kéwel (T. s. scriptus) and imbabala (T. s. sylvaticus) bushbuck. We found strong evidence for gene flow in at least four of eight nodes on the species tree. Among the six phenotypic traits assessed here, only habitat type mapped onto the species tree without homoplasy, showing that trait evolution was the result of complex patterns of divergence, introgression and convergent evolution.

Introgression and incomplete lineage sorting blurred phylogenetic relationships across the genomes of sclerophyllous oaks from southwest China

Resolving evolutionary relationships among closely related species with interspecific gene flow is challenging. Genome-scale data provide opportunities to clarify complex evolutionary relationships in closely related species and to observe variations in species relationships across the genomes of such species. The Himalayan-Hengduan subalpine oaks have a nearly completely sympatric distribution in southwest China and probably constitute a syngameon. In this study, we mapped resequencing data from different species in this group to the Quercus aquifolioides reference genome to obtain a high-quality filtered single nucleotide polymorphism (SNP) dataset. We also assembled their plastomes. We reconstructed their phylogenetic relationships, explored the level and pattern of introgression among these species and investigated gene tree variation in the genomes of these species using sliding windows. The same or closely related plastomes were found to be shared extensively among different species within a specific geographical area. Phylogenomic analyses of genome-wide SNP data found that most oaks in the Himalayan-Hengduan subalpine clade showed genetic coherence, but several species were found to be connected by introgression. The gene trees obtained using sliding windows showed that the phylogenetic relationships in the genomes of oaks are highly heterogeneous and therefore highly obscured. Our study found that all the oaks of the Himalayan-Hengduan subalpine clade from southwest China form a syngameon. The obscured phylogenetic relationships observed empirically across the genome are best explained by interspecific gene flow in conjunction with incomplete lineage sorting.

No sexual pheromones in Anopheles mosquitoes?

Swarming behavior is the cornerstone of the anopheline mating system. At dusk, males congregate in monospecific swarms in which females come to find a mate once in their lives. Although many Anopheles species coexist in sympatry, hybrids are infrequent, suggesting the existence of strong premating reproductive barriers. Chemical cues, particularly pheromones, often play a crucial role in bringing sexes together in a species-specific manner among insects. While the existence of sexual pheromones in Anopheles species has been postulated, only a few studies developed experimental designs to investigate their presence. Here, we discuss the contrasting and debatable findings regarding both long-range and contact sex pheromones in the context of swarm ecology in Anopheles species.

Phylogenetic Conflict Between Species Tree and Maternally Inherited Gene Trees in a Clade of Emberiza Buntings (Aves: Emberizidae)

Different genomic regions may reflect conflicting phylogenetic topologies primarily due to incomplete lineage sorting and/or gene flow. Genomic data are necessary to reconstruct the true species tree and explore potential causes of phylogenetic conflict. Here, we investigate the phylogenetic relationships of 4 Emberiza species (Aves: Emberizidae) and discuss the potential causes of the observed mitochondrial non-monophyly of Emberiza godlewskii (Godlewski's bunting) using phylogenomic analyses based on whole genome resequencing data from 41 birds. Analyses based on both the whole mitochondrial genome and ~39 kilobases from the non-recombining W chromosome reveal sister relationships between each the northern and southern populations of E. godlewskii with E. cioides and E. cia, respectively. In contrast, the monophyly of E. godlewskii is reflected by the phylogenetic signal of autosomal and Z chromosomal sequence data as well as demographic inference analyses, which-in combination-support the following tree topology: ([{E. godlewskii, E. cia}, E. cioides], E. jankowskii). Using D-statistics, we detected multiple gene flow events among different lineages, indicating pervasive introgressive hybridization within this clade. Introgression from an unsampled lineage that is sister to E. cioides or introgression from an unsampled mitochondrial + W chromosomal lineage of E. cioides into northern E. godlewskii may explain the phylogenetic conflict between the species tree estimated from genome-wide data versus mtDNA/W tree topologies. These results underscore the importance of using genomic data for phylogenetic reconstruction and species delimitation.

Complex Polyploids: Origins, Genomic Composition, and Role of Introgressed Alleles

Introgression allows polyploid species to acquire new genomic content from diploid progenitors or from other unrelated diploid or polyploid lineages, contributing to genetic diversity and facilitating adaptive allele discovery. In some cases, high levels of introgression elicit the replacement of large numbers of alleles inherited from the polyploid's ancestral species, profoundly reshaping the polyploid's genomic composition. In such complex polyploids, it is often difficult to determine which taxa were the progenitor species and which taxa provided additional introgressive blocks through subsequent hybridization. Here, we use population-level genomic data to reconstruct the phylogenetic history of Betula pubescens (downy birch), a tetraploid species often assumed to be of allopolyploid origin and which is known to hybridize with at least four other birch species. This was achieved by modeling polyploidization and introgression events under the multispecies coalescent and then using an approximate Bayesian computation rejection algorithm to evaluate and compare competing polyploidization models. We provide evidence that B. pubescens is the outcome of an autoploid genome doubling event in the common ancestor of B. pendula and its extant sister species, B. platyphylla, that took place approximately 178,000-188,000 generations ago. Extensive hybridization with B. pendula, B. nana, and B. humilis followed in the aftermath of autopolyploidization, with the relative contribution of each of these species to the B. pubescens genome varying markedly across the species' range. Functional analysis of B. pubescens loci containing alleles introgressed from B. nana identified multiple genes involved in climate adaptation, while loci containing alleles derived from B. humilis revealed several genes involved in the regulation of meiotic stability and pollen viability in plant species.

Inference of Locus-Specific Population Mixtures from Linked Genome-Wide Allele Frequencies

Admixture between populations and species is common in nature. Since the influx of new genetic material might be either facilitated or hindered by selection, variation in mixture proportions along the genome is expected in organisms undergoing recombination. Various graph-based models have been developed to better understand these evolutionary dynamics of population splits and mixtures. However, current models assume a single mixture rate for the entire genome and do not explicitly account for linkage. Here, we introduce TreeSwirl, a novel method for inferring branch lengths and locus-specific mixture proportions by using genome-wide allele frequency data, assuming that the admixture graph is known or has been inferred. TreeSwirl builds upon TreeMix that uses Gaussian processes to estimate the presence of gene flow between diverged populations. However, in contrast to TreeMix, our model infers locus-specific mixture proportions employing a hidden Markov model that accounts for linkage. Through simulated data, we demonstrate that TreeSwirl can accurately estimate locus-specific mixture proportions and handle complex demographic scenarios. It also outperforms related D- and f-statistics in terms of accuracy and sensitivity to detect introgressed loci.

Genetic surveillance of insecticide resistance in African Anopheles populations to inform malaria vector control

Insecticide resistance in malaria vector populations poses a major threat to malaria control, which relies largely on insecticidal interventions. Contemporary vector-control strategies focus on combatting resistance using multiple insecticides with differing modes of action within the mosquito. However, diverse genetic resistance mechanisms are present in vector populations, and continue to evolve. Knowledge of the spatial distribution of these genetic mechanisms, and how they impact the efficacy of different insecticidal products, is critical to inform intervention deployment decisions. We developed a catalogue of genetic-resistance mechanisms in African malaria vectors that could guide molecular surveillance. We highlight situations where intervention deployment has led to resistance evolution and spread, and identify challenges in understanding and mitigating the epidemiological impacts of resistance.

Recombination rate variation shapes genomic variability of phylogeographic structure in a widespread North American songbird (Aves: Certhia americana)

The nonrandom distribution of chromosomal characteristics and functional elements-genomic architecture-impacts the relative strengths and impacts of population genetic processes across the genome. Due to this relationship, genomic architecture has the potential to shape variation in population genetic structure across the genome. Population genetic structure has been shown to vary across the genome in a variety of taxa, but this body of work has largely focused on pairwise population genomic comparisons between closely related taxa. Here, we used whole genome sequencing of seven phylogeographically structured populations of a North American songbird, the Brown Creeper (Certhia americana), to determine the impacts of genomic architecture on phylogeographic structure variation across the genome. Using multiple methods to infer phylogeographic structure-ordination, clustering, and phylogenetic methods-we found that recombination rate variation explained a large proportion of phylogeographic structure variation. Genomic regions with low recombination showed phylogeographic structure consistent with the genome-wide pattern. In regions with high recombination, we found strong phylogeographic structure, but with discordant patterns relative to the genome-wide pattern. In regions with high recombination rate, we found that populations with small effective population sizes evolve relatively more rapidly than larger populations, leading to discordant signatures of phylogeographic structure. These results suggest that the interplay between recombination rate variation and effective population sizes shape the relative impacts of selection and genetic drift in different parts of the genome. Overall, the combined interactions of population genetic processes, genomic architecture, and effective population sizes shape patterns of variability in phylogeographic structure across the genome of the Brown Creeper.

A new species of Populus and the extensive hybrid speciation arising from it on the Qinghai-Tibet Plateau

While the diversity of species formation is broadly acknowledged, significant debate exists regarding the universal nature of hybrid species formation. Through an 18-year comprehensive study of all Populus species on the Qinghai-Tibet Plateau, 23 previously recorded species and 8 new species were identified. Based on morphological characteristics, these can be classified into three groups: species in section Leucoides, species with large leaves, and species with small leaves in section Tacamahaca. By conducting whole-genome re-sequencing of 150 genotypes from these 31 species, 2.28 million single nucleotide polymorphisms (SNPs) were identified. Phylogenetic analysis utilizing these SNPs not only revealed a highly intricate evolutionary network within the large-leaf species of section Tacamahaca but also confirmed that a new species, P. curviserrata, naturally hybridized with P. cathayana, P. szechuanica, and P. ciliata, resulting in 11 hybrid species. These findings indicate the widespread occurrence of hybrid species formation within this genus, with hybridization serving as a key evolutionary mechanism for Populus on the plateau. A novel hypothesis, "Hybrid Species Exterminating Their Ancestral Species (HSEAS)," is introduced to explain the mechanisms of hybrid species formation at three different scales: the entire plateau, the southeastern mountain region, and individual river valleys.

Schistosoma haematobium and Schistosoma bovis first generation hybrids undergo gene expressions changes consistent with species compatibility and heterosis

When two species hybridize, the two parental genomes are brought together and some alleles might interact for the first time. To date, the extent of the transcriptomic changes in first hybrid generations, along with their functional outcome constitute an important knowledge gap, especially in parasite species. Here we explored the molecular and functional outcomes of hybridization in first-generation hybrids between the blood fluke parasites Schistosoma haematobium and S. bovis. Through a transcriptomic approach, we measured gene expression in both parental species and hybrids. We described and quantified expression profiles encountered in hybrids along with the main biological processes impacted. Up to 7,100 genes fell into a particular hybrid expression profile (intermediate between the parental expression levels, over-expressed, under-expressed, or expressed like one of the parental lines). Most of these genes were different depending on the direction of the parental cross (S. bovis mother and S. haematobium father or the reverse) and depending on the sex. For a given sex and cross direction, the vast majority of genes were hence unassigned to a hybrid expression profile: either they were differentially expressed genes but not typical of any hybrid expression profiles or they were not differentially expressed neither between hybrids and parental lines nor between parental lines. The most prevalent profile of gene expression in hybrids was the intermediate one (24% of investigated genes). These results suggest that transcriptomic compatibility between S. haematobium and S. bovis remains quite high. We also found support for an over-dominance model (over- and under-expressed genes in hybrids compared to parental lines) potentially associated with heterosis. In females in particular, processes such as reproductive processes, metabolism and cell interactions as well as signaling pathways were indeed affected. Our study hence provides new insight on the biology of Schistosoma hybrids with evidences supporting compatibility and heterosis.

Evolution of Chromosomal Inversions across an Avian Radiation

Chromosomal inversions are structural mutations that can play a prominent role in adaptation and speciation. Inversions segregating across species boundaries (trans-species inversions) are often taken as evidence for ancient balancing selection or adaptive introgression, but can also be due to incomplete lineage sorting. Using whole-genome resequencing data from 18 populations of 11 recognized munia species in the genus Lonchura (N = 176 individuals), we identify four large para- and pericentric inversions ranging in size from 4 to 20 Mb. All four inversions cosegregate across multiple species and predate the numerous speciation events associated with the rapid radiation of this clade across the prehistoric Sahul (Australia, New Guinea) and Bismarck Archipelago. Using coalescent theory, we infer that trans-specificity is improbable for neutrally segregating variation despite substantial incomplete lineage sorting characterizing this young radiation. Instead, the maintenance of all three autosomal inversions (chr1, chr5, and chr6) is best explained by selection acting along ecogeographic clines not observed for the collinear parts of the genome. In addition, the sex chromosome inversion largely aligns with species boundaries and shows signatures of repeated positive selection for both alleles. This study provides evidence for trans-species inversion polymorphisms involved in both adaptation and speciation. It further highlights the importance of informing selection inference using a null model of neutral evolution derived from the collinear part of the genome.

Host preference patterns in domestic and wild settings: Insights into Anopheles feeding behavior

The adaptation of Anopheles malaria vectors to domestic settings is directly linked to their ability to feed on humans. The strength of this species-habitat association is unequal across the species within the genus, with the major vectors being particularly dependent on humans. However, our understanding of how blood-feeding behavior interacts with and adapts to environmental settings, including the presence of humans, remains limited. Using a field-based approach, we first investigated Anopheles community structure and feeding behavior patterns in domestic and sylvatic settings in La Lopé National Park in Gabon, Central Africa. We characterized the preference indices using a dual-host choice sampling approach across mosquito species, habitats, and seasons. We then quantified the plastic biting behavior of mosquito species in each habitat. We collected individuals from 16 Anopheles species that exhibited significant differences in species composition and abundance between sylvatic and domestic settings. The host-seeking behavior also varied among the seven most abundant species. The general attractiveness to each host, human or animal, remained relatively constant for each species, but with significant variations between habitats across species. These variations, to more generalist and to more anthropophilic behavior, were related to seasonal changes and distance from the village, respectively. Finally, we pointed out that the host choice of major malaria vectors changed in the absence of humans, revealing a plastic feeding behavior of these species. This study highlights the effect of humans on Anopheles distribution and feeding evolution. The characterization of feeding behavior in wild and domestic settings provides opportunities to better understand the interplay between genetic determinants of host preference and ecological factors. Our findings suggest that protected areas may offer alternative thriving conditions to major malaria vectors.

Introgression Underlies Phylogenetic Uncertainty But Not Parallel Plumage Evolution in a Recent Songbird Radiation

Instances of parallel phenotypic evolution offer great opportunities to understand the evolutionary processes underlying phenotypic changes. However, confirming parallel phenotypic evolution and studying its causes requires a robust phylogenetic framework. One such example is the "black-and-white wagtails," a group of 5 species in the songbird genus Motacilla: 1 species, Motacilla alba, shows wide intra-specific plumage variation, while the 4r others form 2 pairs of very similar-looking species (M. aguimp + M. samveasnae and M. grandis + M. maderaspatensis, respectively). However, the 2 species in each of these pairs were not recovered as sisters in previous phylogenetic inferences. Their relationships varied depending on the markers used, suggesting that gene tree heterogeneity might have hampered accurate phylogenetic inference. Here, we use whole genome resequencing data to explore the phylogenetic relationships within this group, with a special emphasis on characterizing the extent of gene tree heterogeneity and its underlying causes. We first used multispecies coalescent methods to generate a "complete evidence" phylogenetic hypothesis based on genome-wide variants, while accounting for incomplete lineage sorting (ILS) and introgression. We then investigated the variation in phylogenetic signal across the genome to quantify the extent of discordance across genomic regions and test its underlying causes. We found that wagtail genomes are mosaics of regions supporting variable genealogies, because of ILS and inter-specific introgression. The most common topology across the genome, supporting M. alba and M. aguimp as sister species, appears to be influenced by ancient introgression. Additionally, we inferred another ancient introgression event, between M. alba and M. grandis. By combining results from multiple analyses, we propose a phylogenetic network for the black-and-white wagtails that confirms that similar phenotypes evolved in non-sister lineages, supporting parallel plumage evolution. Furthermore, the inferred reticulations do not connect species with similar plumage coloration, suggesting that introgression does not underlie parallel plumage evolution in this group. Our results demonstrate the importance of investing genome-wide patterns of gene tree heterogeneity to help understand the mechanisms underlying phenotypic evolution. [Gene tree heterogeneity; incomplete lineage sorting; introgression; parallel evolution; phylogenomics; plumage evolution; wagtails.].

Detection of Ghost Introgression Requires Exploiting Topological and Branch Length Information

In recent years, the study of hybridization and introgression has made significant progress, with ghost introgression-the transfer of genetic material from extinct or unsampled lineages to extant species-emerging as a key area for research. Accurately identifying ghost introgression, however, presents a challenge. To address this issue, we focused on simple cases involving 3 species with a known phylogenetic tree. Using mathematical analyses and simulations, we evaluated the performance of popular phylogenetic methods, including HyDe and PhyloNet/MPL, and the full-likelihood method, Bayesian Phylogenetics and Phylogeography (BPP), in detecting ghost introgression. Our findings suggest that heuristic approaches relying on site-pattern counts or gene-tree topologies struggle to differentiate ghost introgression from introgression between sampled non-sister species, frequently leading to incorrect identification of donor and recipient species. The full-likelihood method BPP uses multilocus sequence alignments directly-hence taking into account both gene-tree topologies and branch lengths, by contrast, is capable of detecting ghost introgression in phylogenomic datasets. We analyzed a real-world phylogenomic dataset of 14 species of Jaltomata (Solanaceae) to showcase the potential of full-likelihood methods for accurate inference of introgression.

Geogenomic Predictors of Genetree Heterogeneity Explain Phylogeographic and Introgression History: A Case Study in an Amazonian Bird (Thamnophilus aethiops)

Can knowledge about genome architecture inform biogeographic and phylogenetic inference? Selection, drift, recombination, and gene flow interact to produce a genomic landscape of divergence wherein patterns of differentiation and genealogy vary nonrandomly across the genomes of diverging populations. For instance, genealogical patterns that arise due to gene flow should be more likely to occur on smaller chromosomes, which experience high recombination, whereas those tracking histories of geographic isolation (reduced gene flow caused by a barrier) and divergence should be more likely to occur on larger and sex chromosomes. In Amazonia, populations of many bird species diverge and introgress across rivers, resulting in reticulated genomic signals. Herein, we used reduced representation genomic data to disentangle the evolutionary history of 4 populations of an Amazonian antbird, Thamnophilus aethiops, whose biogeographic history was associated with the dynamic evolution of the Madeira River Basin. Specifically, we evaluate whether a large river capture event ca. 200 Ka, gave rise to reticulated genealogies in the genome by making spatially explicit predictions about isolation and gene flow based on knowledge about genomic processes. We first estimated chromosome-level phylogenies and recovered 2 primary topologies across the genome. The first topology (T1) was most consistent with predictions about population divergence and was recovered for the Z-chromosome. The second (T2), was consistent with predictions about gene flow upon secondary contact. To evaluate support for these topologies, we trained a convolutional neural network to classify our data into alternative diversification models and estimate demographic parameters. The best-fit model was concordant with T1 and included gene flow between non-sister taxa. Finally, we modeled levels of divergence and introgression as functions of chromosome length and found that smaller chromosomes experienced higher gene flow. Given that (1) genetrees supporting T2 were more likely to occur on smaller chromosomes and (2) we found lower levels of introgression on larger chromosomes (and especially the Z-chromosome), we argue that T1 represents the history of population divergence across rivers and T2 the history of secondary contact due to barrier loss. Our results suggest that a significant portion of genomic heterogeneity arises due to extrinsic biogeographic processes such as river capture interacting with intrinsic processes associated with genome architecture. Future phylogeographic studies would benefit from accounting for genomic processes, as different parts of the genome reveal contrasting, albeit complementary histories, all of which are relevant for disentangling the intricate geogenomic mechanisms of biotic diversification. [Amazonia; biogeography; demographic modeling; gene flow; gene tree; genome architecture; geogenomics; introgression; linked selection; neural network; phylogenomic; phylogeography; reproductive isolation; speciation; species tree.].

Genomic Signatures of Microgeographic Adaptation in Anopheles coluzzii Along an Anthropogenic Gradient in Gabon

Species distributed across heterogeneous environments often evolve locally adapted populations, but understanding how these persist in the presence of homogenizing gene flow remains puzzling. In Gabon, Anopheles coluzzii, a major African malaria mosquito is found along an ecological gradient, including a sylvatic population, away of any human presence. This study identifies into the genomic signatures of local adaptation in populations from distinct environments including the urban area of Libreville, and two proximate sites 10km apart in the La Lopé National Park (LLP), a village and its sylvatic neighborhood. Whole genome re-sequencing of 96 mosquitoes unveiled ∼ 5.7millions high-quality single nucleotide polymorphisms. Coalescent-based demographic analyses suggest an ∼ 8,000-year-old divergence between Libreville and La Lopé populations, followed by a secondary contact ( ∼ 4,000 ybp) resulting in asymmetric effective gene flow. The urban population displayed reduced effective size, evidence of inbreeding, and strong selection pressures for adaptation to urban settings, as suggested by the hard selective sweeps associated with genes involved in detoxification and insecticide resistance. In contrast, the two geographically proximate LLP populations showed larger effective sizes, and distinctive genomic differences in selective signals, notably soft-selective sweeps on the standing genetic variation. Although neutral loci and chromosomal inversions failed to discriminate between LLP populations, our findings support that microgeographic adaptation can swiftly emerge through selection on standing genetic variation despite high gene flow. This study contributes to the growing understanding of evolution of populations in heterogeneous environments amid ongoing gene flow and how major malaria mosquitoes adapt to human.

Significance

Anopheles coluzzii , a major African malaria vector, thrives from humid rainforests to dry savannahs and coastal areas. This ecological success is linked to its close association with domestic settings, with human playing significant roles in driving the recent urban evolution of this mosquito. Our research explores the assumption that these mosquitoes are strictly dependent on human habitats, by conducting whole-genome sequencing on An. coluzzii specimens from urban, rural, and sylvatic sites in Gabon. We found that urban mosquitoes show de novo genetic signatures of human-driven vector control, while rural and sylvatic mosquitoes exhibit distinctive genetic evidence of local adaptations derived from standing genetic variation. Understanding adaptation mechanisms of this mosquito is therefore crucial to predict evolution of vector control strategies.

The contours of evolution: In defence of Darwin's tree of life paradigm

Both the concept of a Darwinian tree of life (TOL) and the possibility of its accurate reconstruction have been much criticized. Criticisms mostly revolve around the extensive occurrence of lateral gene transfer (LGT), instances of uptake of complete organisms to become organelles (with the associated subsequent gene transfer to the nucleus), as well as the implications of more subtle aspects of the biological species concept. Here we argue that none of these criticisms are sufficient to abandon the valuable TOL concept and the biological realities it captures. Especially important is the need to conceptually distinguish between organismal trees and gene trees, which necessitates incorporating insights of widely occurring LGT into modern evolutionary theory. We demonstrate that all criticisms, while based on important new findings, do not invalidate the TOL. After considering the implications of these new insights, we find that the contours of evolution are best represented by a TOL.

The Anopheles coluzzii range extends into Kenya: detection, insecticide resistance profiles and population genetic structure in relation to conspecific populations in West and Central Africa

BACKGROUND

Anopheles coluzzii is a primary vector of malaria found in West and Central Africa, but its presence has hitherto never been documented in Kenya. A thorough understanding of vector bionomics is important as it enables the implementation of targeted and effective vector control interventions. Malaria vector surveillance efforts in the country have tended to focus on historically known primary vectors. The current study sought to determine the taxonomic status of samples collected from five different malaria epidemiological zones in Kenya as well as describe the population genetic structure and insecticide resistance profiles in relation to other An. coluzzii populations.

METHODS

Mosquitoes were sampled as larvae from Busia, Kwale, Turkana, Kirinyaga and Kiambu counties, representing the range of malaria endemicities in Kenya, in 2019 and 2021 and emergent adults analysed using Whole Genome Sequencing (WGS) data processed in accordance with the Anopheles gambiae 1000 Genomes Project phase 3. Where available, historical samples from the same sites were included for WGS. Comparisons were made with An. coluzzii cohorts from West and Central Africa.

RESULTS

This study reports the detection of An. coluzzii for the first time in Kenya. The species was detected in Turkana County across all three time points from which samples were analyzed and its presence confirmed through taxonomic analysis. Additionally, there was a lack of strong population genetic differentiation between An. coluzzii from Kenya and those from the more northerly regions of West and Central Africa, suggesting they represent a connected extension to the known species range. Mutations associated with target-site resistance to DDT and pyrethroids and metabolic resistance to DDT were found at high frequencies up to 64%. The profile and frequencies of the variants observed were similar to An. coluzzii from West and Central Africa but the ace-1 mutation linked to organophosphate and carbamate resistance present in An. coluzzii from coastal West Africa was absent in Kenya.

CONCLUSIONS

These findings emphasize the need for the incorporation of genomics in comprehensive and routine vector surveillance to inform on the range of malaria vector species, and their insecticide resistance status to inform the choice of effective vector control approaches.

Two Nested Inversions in the X Chromosome Differentiate the Dominant Malaria Vectors in Europe, Anopheles atroparvus and Anopheles messeae

The Maculipennis subgroup of malaria mosquitoes includes both dominant malaria vectors and non-vectors in Eurasia. Understanding the genetic factors, particularly chromosomal inversions, that differentiate Anopheles species can provide valuable insights for vector control strategies. Although autosomal inversions between the species in this subgroup have been characterized based on the chromosomal banding patterns, the number and positions of rearrangements in the X chromosome remain unclear due to the divergent banding patterns. Here, we identified two large X chromosomal inversions, approximately 13 Mb and 10 Mb in size, using fluorescence in situ hybridization. The inversion breakpoint regions were mapped by hybridizing 53 gene markers with polytene chromosomes of An. messeae. The DNA probes were designed based on gene sequences from the annotated An. atroparvus genome. The two nested inversions resulted in five syntenic blocks. Only two small syntenic blocks, which encompass 181 annotated genes in the An. atroparvus genome, changed their position and orientation in the X chromosome. The analysis of the An. atroparvus genome revealed an enrichment of gene ontology terms associated with immune system and mating behavior in the rearranged syntenic blocks. Additionally, the enrichment of DNA transposons was found in sequences homologous to three of the four breakpoint regions. This study demonstrates the successful application of the physical genome mapping approach to identify rearrangements that differentiate species in insects with polytene chromosomes.

Metazoa-level USCOs as markers in species delimitation and classification

Metazoa-level universal single-copy orthologs (mzl-USCOs) are universally applicable markers for DNA taxonomy in animals that can replace or supplement single-gene barcodes. Previously, mzl-USCOs from target enrichment data were shown to reliably distinguish species. Here, we tested whether USCOs are an evenly distributed, representative sample of a given metazoan genome and therefore able to cope with past hybridization events and incomplete lineage sorting. This is relevant for coalescent-based species delimitation approaches, which critically depend on the assumption that the investigated loci do not exhibit autocorrelation due to physical linkage. Based on 239 chromosome-level assembled genomes, we confirmed that mzl-USCOs are genetically unlinked for practical purposes and a representative sample of a genome in terms of reciprocal distances between USCOs on a chromosome and of distribution across chromosomes. We tested the suitability of mzl-USCOs extracted from genomes for species delimitation and phylogeny in four case studies: Anopheles mosquitos, Drosophila fruit flies, Heliconius butterflies and Darwin's finches. In almost all instances, USCOs allowed delineating species and yielded phylogenies that corresponded to those generated from whole genome data. Our phylogenetic analyses demonstrate that USCOs may complement single-gene DNA barcodes and provide more accurate taxonomic inferences. Combining USCOs from sources that used different versions of ortholog reference libraries to infer marker orthology may be challenging and, at times, impact taxonomic conclusions. However, we expect this problem to become less severe as the rapidly growing number of reference genomes provides a better representation of the number and diversity of organismal lineages.

Lessons from assembling UCEs: A comparison of common methods and the case of Clavinomia (Halictidae)

Sequence data assembly is a foundational step in high-throughput sequencing, with untold consequences for downstream analyses. Despite this, few studies have interrogated the many methods for assembling phylogenomic UCE data for their comparative efficacy, or for how outputs may be impacted. We study this by comparing the most commonly used assembly methods for UCEs in the under-studied bee lineage Nomiinae and a representative sampling of relatives. Data for 63 UCE-only and 75 mixed taxa were assembled with five methods, including ABySS, HybPiper, SPAdes, Trinity and Velvet, and then benchmarked for their relative performance in terms of locus capture parameters and phylogenetic reconstruction. Unexpectedly, Trinity and Velvet trailed the other methods in terms of locus capture and DNA matrix density, whereas SPAdes performed favourably in most assessed metrics. In comparison with SPAdes, the guided-assembly approach HybPiper generally recovered the highest quality loci but in lower numbers. Based on our results, we formally move Clavinomia to Dieunomiini and render Epinomia once more a subgenus of Dieunomia. We strongly advise that future studies more closely examine the influence of assembly approach on their results, or, minimally, use better-performing assembly methods such as SPAdes or HybPiper. In this way, we can move forward with phylogenomic studies in a more standardized, comparable manner.

Speciation within the Anopheles gambiae complex: high-throughput whole genome sequencing reveals evidence of a putative new cryptic taxon in 'far-west' Africa

The two main Afrotropical malaria vectors - Anopheles coluzzii and An. gambiae - are genetically distinct and reproductively isolated across West Africa. However, populations at the western extreme of their range are assigned as "intermediate" between the two species by whole genome sequence (WGS) data, and as hybrid forms by conventional molecular diagnostics. By exploiting WGS data from 1,190 specimens collected across west Africa via the Anopheles gambiae 1000 Genomes network, we identify a novel putative taxon in the far-west (provisionally named Bissau molecular form), which did not arise by admixture but rather originated at the same time as the split between An. coluzzii and An. gambiae. Intriguingly, these populations lack insecticide resistance mechanisms commonly observed in the two main species. These findings lead to a change of perspective on malaria vector species in the far-west region with potential for epidemiological implications, and a new challenge for genetic-based mosquito control approaches.