Phylogeographic Inference Research Articles

Testing the Effect of Sampling Effort on Inferring Phylogeographic History in Psolodesmus mandarinus (Calopterygidae, Odonata)

Phylogeographic studies have revealed spatial genetic structure and inferred geographical processes that may have generated genetic diversity and divergence. These study results have implications not only on the processes that generate intraspecific and interspecific diversity but also on the essential integrals for defining evolutionary entities (e.g., species). However, the resulting phylogeographic inferences might be impacted by the sampling design, i.e., the number of individuals per population and the number of geographic populations studied. The effect of sampling bias on phylogeographic inferences remains poorly explored. With a comprehensive sampling design (including 186 samples from 56 localities), we studied the phylogeographic history of a Taiwanese endemic damselfly, Psolodesmus mandarinus, with a specific focus on testing the impact of the sampling design on phylogeographic inference. We found a significant difference in the genetic structure of eastern and western populations separated by the Central Mountain Range (CMR) of Taiwan. However, isolation by the CMR did not lead to reciprocally monophyletic geographic populations. We further showed that, when only a subset of individuals was randomly included in the study, monophyletic geographic populations were obtained. Furthermore, historical demographic expansion could become undetectable when only a subset of samples was used in the analyses. Our results demonstrate the impact of sampling design on phylogeographic inferences. Future studies need to be cautious when inferring the effect of isolation by a physical barrier.

SPREAD 4: online visualisation of pathogen phylogeographic reconstructions.

Phylogeographic analyses aim to extract information about pathogen spread from genomic data, and visualising spatio-temporal reconstructions is a key aspect of this process. Here we present SPREAD 4, a feature-rich web-based application that visualises estimates of pathogen dispersal resulting from Bayesian phylogeographic inference using BEAST on a geographic map, offering zoom-and-filter functionality and smooth animation over time. SPREAD 4 takes as input phylogenies with both discrete and continuous location annotation and offers customised visualisation as well as generation of publication-ready figures. SPREAD 4 now features account-based storage and easy sharing of visualisations by means of unique web addresses. SPREAD 4 is intuitive to use and is available online at https://spreadviz.org, with an accompanying web page containing answers to frequently asked questions at https://beast.community/spread4.

Importation and circulation of rubella virus lineages 1E-L2 and 2B-L2c between 2018 and 2021 in China: Virus evolution and spatial-temporal transmission characteristics.

To better understand the importation and circulation patterns of rubella virus lineages 1E-L2 and 2B-L2c circulating in China since 2018, 3,312 viral strains collected from 27 out of 31 provinces in China between 2018 and 2021 were sequenced and analyzed with the representative international strains of lineages 1E-L2 and 2B-L2c based on genotyping region. Time-scale phylogenetic analysis revealed that the global lineages 1E-L2 and 2B-L2c presented distinct evolutionary patterns. Lineage 1E-L2 circulated in relatively limited geographical areas (mainly Asia) and showed geographical and temporal clustering, while lineage 2B-L2c strains circulated widely throughout the world and exhibited a complicated topology with several independently evolved branches. Furthermore, both lineages showed extensive international transmission activities, and phylogeographic inference provided evidence that lineage 1E-L2 strains circulating in China possibly originated from Japan, while the source of lineage 2B-L2c isolated since 2018 is still unclear. After importation into China in 2018, the spread of lineage 1E-L2 presented a three-stage transmission pattern from southern to northern China, whereas lineage 2B-L2c spread from a single point in western China to all the other four regions. These two transmission patterns allowed both imported lineages to spread rapidly across China during the 2018-9 rubella epidemic and eventually established endemic circulations. This study provides critical scientific data for rubella control and elimination in China and worldwide.

Recombination-aware phylogeographic inference using the structured coalescent with ancestral recombination.

Movement of individuals between populations or demes is often restricted, especially between geographically isolated populations. The structured coalescent provides an elegant theoretical framework for describing how movement between populations shapes the genealogical history of sampled individuals and thereby structures genetic variation within and between populations. However, in the presence of recombination an individual may inherit different regions of their genome from different parents, resulting in a mosaic of genealogical histories across the genome, which can be represented by an Ancestral Recombination Graph (ARG). In this case, different genomic regions may have different ancestral histories and so different histories of movement between populations. Recombination therefore poses an additional challenge to phylogeographic methods that aim to reconstruct the movement of individuals from genealogies, although also a potential benefit in that different loci may contain additional information about movement. Here, we introduce the Structured Coalescent with Ancestral Recombination (SCAR) model, which builds on recent approximations to the structured coalescent by incorporating recombination into the ancestry of sampled individuals. The SCAR model allows us to infer how the migration history of sampled individuals varies across the genome from ARGs, and improves estimation of key population genetic parameters such as population sizes, recombination rates and migration rates. Using the SCAR model, we explore the potential and limitations of phylogeographic inference using full ARGs. We then apply the SCAR to lineages of the recombining fungus Aspergillus flavus sampled across the United States to explore patterns of recombination and migration across the genome.

Comparative mitogenomics elucidates the population genetic structure of Amblyomma testudinarium in Japan and a closely related Amblyomma species in Myanmar.

Ticks are the second most important vector capable of transmitting diseases affecting the health of both humans and animals. Amblyomma testudinarium Koch 1844 (Acari: Ixodidae), is a hard tick species having a wide geographic distribution in Asia. In this study, we analyzed the composition of A. testudinarium whole mitogenomes from various geographical regions in Japan and investigated the population structure, demographic patterns, and phylogeographic relationship with other ixodid species. In addition, we characterized a potentially novel tick species closely related to A. testudinarium from Myanmar. Phylogeographic inference and evolutionary dynamics based on the 15 mitochondrial coding genes supported that A. testudinarium population in Japan is resolved into a star‐like haplogroup and suggested a distinct population structure of A. testudinarium from Amami island in Kyushu region. Correlation analysis using Mantel test statistics showed that no significant correlation was observed between the genetic and geographic distances calculated between the A. testudinarium population from different localities in Japan. Finally, demographic analyses, including mismatch analysis and Tajima’s D test, suggested a possibility of recent population expansion occurred within Japanese haplogroup after a bottleneck event. Although A. testudinarium has been considered widespread and common in East and Southeast Asia, the current study suggested that potentially several cryptic Amblyomma spp. closely related to A. testudinarium are present in Asia.

Accommodating sampling location uncertainty in continuous phylogeography.

Phylogeographic inference of the dispersal history of viral lineages offers key opportunities to tackle epidemiological questions about the spread of fast-evolving pathogens across human, animal and plant populations. In continuous space, i.e. when locations are specified by longitude and latitude, these reconstructions are however often limited by the availability or accessibility of precise sampling locations required for such spatially explicit analyses. We here review the different approaches that can be considered when genomic sequences are associated with a geographic area of sampling instead of precise coordinates. In particular, we describe and compare the approaches to define homogeneous and heterogeneous prior ranges of sampling coordinates.

Human‐mediated eco‐evolutionary processes of the herbivorous insect Hyalopterus arundiniformis during the Holocene

AbstractAimPleistocene climatic fluctuations have been identified as a dominant factor in evolutionary and ecological processes underlying contemporary species distribution. In contrast, anthropogenic activity during the Holocene has been largely neglected in phylogeographic inference. Here, by integrating archaeological evidence, phylogeographic analyses and a climatic niche‐based approach, the role of Holocene human activity in determining the evolutionary trajectory of insect pests is investigated using Hyalopterus arundiniformis (HAR) as a model.LocationEurasia.MethodsWe compiled archaeological records and historical data for the host plant (i.e. peach) of HAR. Phylogeographic structure, genetic diversity and gene flow were estimated. Historical demographic dynamics were explored using coalescent‐based analyses. The origin and dispersal history were assessed using an approximate Bayesian computation (ABC) approach. The changes in climatic niche breadth during geographic expansion were evaluated via niche comparison analyses.ResultsOur results revealed that HAR originated in Southwest China and dispersed to Southeast and Northwest China, and it was then introduced to western Eurasia from Northwest China. The climatic niche shifted significantly with range expansions, and niche expansion explained the largest proportion of niche dynamics. Population size expansion occurred in the Holocene and coincided with the domestication and trans‐Eurasian trade of peach.Main conclusionsThe geographic range, effective population size and climatic niche breadth of HAR expanded following the human‐mediated dispersal of its host plant during the Holocene. Our results highlight that the expansion of Neolithic agriculture and trans‐Eurasian crop exchange along ancient communication routes may have played an important role in shaping the distribution and demographic dynamics of herbivorous insects. This study also provides essential information for pest management and biological invasions.

Genetic diversity and connectivity of the Ostreid herpesvirus 1 populations in France: A first attempt to phylogeographic inference for a marine mollusc disease.

The genetic diversity of viral populations is a key driver of the spatial and temporal diffusion of viruses; yet, studying the diversity of whole genomes from natural populations still remains a challenge. Phylodynamic approaches are commonly used for RNA viruses harboring small genomes but have only rarely been applied to DNA viruses with larger genomes. Here, we used the Pacific oyster mortality syndrome (a disease that affects oyster farms around the world) as a model to study the genetic diversity of its causative agent, the Ostreid herpesvirus 1 (OsHV-1) in the three main French oyster-farming areas. Using ultra-deep sequencing on individual moribund oysters and an innovative combination of bioinformatics tools, we de novo assembled twenty-one OsHV-1 new genomes. Combining quantification of major and minor genetic variations, phylogenetic analysis, and ancestral state reconstruction of discrete traits approaches, we assessed the connectivity of OsHV-1 viral populations between the three oyster-farming areas. Our results suggest that the Marennes-Oléron Bay represents the main source of OsHV-1 diversity, from where the virus has dispersed to other farming areas, a scenario consistent with current practices of oyster transfers in France. We demonstrate that phylodynamic approaches can be applied to aquatic DNA viruses to determine how epidemiological, immunological, and evolutionary processes act and potentially interact to shape their diversity patterns.

Brine shrimps adrift: historical species turnover in Western Mediterranean Artemia (Anostraca)

Brine shrimps (Artemia) have undergone geographic range and demographic expansions as a result of their interaction with humans since the beginning of salt harvesting. This interaction has favoured the expansion of some species but compromising the survival of others. Mediterranean native populations of Artemia salina from coastal salterns and lagoons are facing the presence and expansion of the introduced and invasive American species Artemia monica (= A. franciscana). However, this species could not be the only threat. Parthenogenetic populations of the Asian species A. urmiana and A. sinica are widespread along the Mediterranean and other areas of the world. In this work, with the use of large cox1 and mitogenomic datasets, phylogenetic and phylogeographic inferences, and a time calibrated tree, we confirmed the Asian origin and recent arrival of the current Western Mediterranean parthenogenetic populations of Artemia. In addition, the replacement of Iberian populations of A. salina by Asiatic parthenogenetic populations lead us to recognize parthenogens as invasive. Current salterns development and commercial importance of Artemia make human-mediated introduction probable. These results demonstrate again the impact that changing human interests have on population expansion or decline of species adapted to anthropogenic habitats. Artemia salina decline makes urgent the implementation of conservation measures such as its use in fish farming and salt production or its inoculation in inland salterns.

Molecular Phylogeny and Historical Biogeography of Petroleuciscus (Teleostei: Leuciscidae) Species in Turkey

In this study, phylogenetic relationships of Petroleuciscus species in Turkey were investigated by using mitochondrial cytochrome b (cyt b) (1140 bp). For this purpose, 84 samples from 13 populations for the cyt b gene were sequenced. The result of of genetic analysis, twenty-one haplotypes were identified and overall haplotype diversity was calculated as Hd: 0.921±0.011. The result of phylogenetic analysis of Petroleuciscus species shows the presence of two clades and four haplogroups. First clade includes Petroleuciscus smyrneaus and P. ninae from the Western Anatolia. Second clade indicated that P. borysthenicus consists of Black sea and Marmara Sea drainages and northern of Aegean Sea drainages. Haplotype network analysis and phylogeographic inference show that Petroleuciscus genus was monophyletic in Turkey.

Evolutionary history and introduction of SARS-CoV-2 Alpha VOC/B.1.1.7 in Pakistan through international travelers.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants continue to emerge, and their identification is important for the public health response to coronavirus disease 2019 (COVID-19). Genomic sequencing provides robust information but may not always be accessible, and therefore, mutation-based polymerase chain reaction (PCR) approaches can be used for rapid identification of known variants. International travelers arriving in Karachi between December 2020 and February 2021 were tested for SARS-CoV-2 by PCR. A subset of positive samples was tested for S-gene target failure (SGTF) on TaqPathTM COVID-19 (Thermo Fisher Scientific) and for mutations using the GSD NovaType SARS-CoV-2 (Eurofins Technologies) assays. Sequencing was conducted on the MinION platform (Oxford Nanopore Technologies). Bayesian phylogeographic inference was performed integrating the patients’ travel history information. Of the thirty-five COVID-19 cases screened, thirteen had isolates with SGTF. The travelers transmitted infection to sixty-eight contact cases. The B.1.1.7 lineage was confirmed through sequencing and PCR. The phylogenetic analysis of sequence data available for six cases included four B.1.1.7 strains and one B.1.36 and B.1.1.212 lineage isolate. Phylogeographic modeling estimated at least three independent B.1.1.7 introductions into Karachi, Pakistan, originating from the UK. B.1.1.212 and B.1.36 were inferred to be introduced either from the UK or the travelers’ layover location. We report the introduction of SARS-CoV-2 B.1.1.7 and other lineages in Pakistan by international travelers arriving via different flight routes. This highlights SARS-CoV-2 transmission through travel, importance of testing, and quarantine post-travel to prevent transmission of new strains, as well as recording detailed patients’ metadata. Such results help inform policies on restricting travel from destinations where new highly transmissible variants have emerged.

Phylogeographic Assessment Reveals Geographic Sources of HIV-1 Dissemination Among Men Who Have Sex With Men in Kenya.

HIV-1 transmission dynamics involving men who have sex with men (MSM) in Africa are not well understood. We investigated the rates of HIV-1 transmission between MSM across three regions in Kenya: Coast, Nairobi, and Nyanza. We analyzed 372 HIV-1 partial pol sequences sampled during 2006–2019 from MSM in Coast (N = 178, 47.9%), Nairobi (N = 137, 36.8%), and Nyanza (N = 57, 15.3%) provinces in Kenya. Maximum-likelihood (ML) phylogenetics and Bayesian inference were used to determine HIV-1 clusters, evolutionary dynamics, and virus migration rates between geographic regions. HIV-1 sub-subtype A1 (72.0%) was most common followed by subtype D (11.0%), unique recombinant forms (8.9%), subtype C (5.9%), CRF 21A2D (0.8%), subtype G (0.8%), CRF 16A2D (0.3%), and subtype B (0.3%). Forty-six clusters (size range 2–20 sequences) were found—half (50.0%) of which had evidence of extensive HIV-1 mixing among different provinces. Data revealed an exponential increase in infections among MSM during the early-to-mid 2000s and stable or decreasing transmission dynamics in recent years (2017–2019). Phylogeographic inference showed significant (Bayes factor, BF > 3) HIV-1 dissemination from Coast to Nairobi and Nyanza provinces, and from Nairobi to Nyanza province. Strengthening HIV-1 prevention programs to MSM in geographic locations with higher HIV-1 prevalence among MSM (such as Coast and Nairobi) may reduce HIV-1 incidence among MSM in Kenya.

Pliocene-Pleistocene dispersal bring along low inter species diversity between Vimba species based on multilocus analysis

This study investigates phylogenetic and phylogeographic relationships of Vimba species using mitochondrial cytochrome b (cyt b) (1023 bp) and cytochrome c oxidase subunit I (COI) barcoding region (652 bp) genes. Ninety-one samples from 36 populations for the cyt b gene and 67 samples from 20 populations for the COI were analyzed. We identified 29 haplotypes and calculated overall haplotype diversity as Hd: 0.907 ± 0.015 for cyt b. We also identified 13 COI haplotypes and calculated overall haplotype diversity as 0.826 ± 0.026 for this marker. The phylogenetic analysis of Vimba species reveals the presence of four clades, based on concatenated cyt b and COI sequences. The first and second clade consist of Vimba vimba Western lineage, and Vimba vimba Caspian lineage, while the third and fourth clade consist of Vimba mirabilis and Vimba melanops. Based on haplotype network analyses and phylogeographic inferences, the Vimba genus is monophyletic, and its species dispersed in the Pleistocene era.

The phylogeographic history of Megistostegium (Malvaceae) in the dry, spiny thickets of southwestern Madagascar using RAD-seq data and ecological niche modeling.

The spiny thicket of southwestern Madagascar represents an extreme and ancient landscape with extraordinary levels of biodiversity and endemism. Few hypotheses exist for explaining speciation in the region and few plant studies have explored hypotheses for species diversification. Here, we investigate three species in the endemic genus Megistostegium (Malvaceae) to evaluate phylogeographic structure and explore the roles of climate, soil, and paleoclimate oscillations on population divergence and speciation throughout the region. We combine phylogenetic and phylogeographic inference of RADseq data with ecological niche modeling across space and time. Population structure is concurrent with major rivers in the region and we identify a new, potentially important biogeographic break coincident with several landscape features. Our data further suggests that niches occupied by species and populations differ substantially across their distribution. Paleodistribution modeling provide evidence that past climatic change could be responsible for the current distribution, population structure, and maintenance of species in Megistostegium.

Bluetongue virus in Europe: the current epidemiological situation

The article reviews the history of BT occurrence in Europe and its present status. It describes the distribution of BT in Europe before 1998, the emergence of BTV in southern and eastern Europe in 1998-2006 and the epidemiology of BT in north-western Europe after 2006. Up to 1998, sporadic cases of BT were noted in Cyprus, on the Iberian Peninsula and on several Greek islands. However, since 1998, probably due to climatic changes, BTV has spread northwards into the Mediterranean Basin, where five BTV serotypes (1, 2, 4, 9 and 16) have been identified. In August 2006, BTV passed for the first time latitude 50°N, and BT outbreaks caused by BTV serotype 8 occurred in the Netherlands, Belgium, Germany, France and Luxembourg. Mass vaccination campaigns implemented in Europe in the spring of 2008 quickly limited the spread of disease caused by BTV-8, and it was eradicated by 2011. However, after a 3-year break, in September 2015, BTV-8 re-emerged in Europe, in central France, and subsequently spread throughout the entire country. In the following years, BTV-8 outbreaks were found in Switzerland, Germany, Belgium and Spain. In addition to BTV-8 outbreaks, BTV serotypes 1, 2, 4, 9 and 16 have recently circulated in Europe. As revealed by phylogeographic inference, the recent spread of BTV in Europe is a consequence of climatic, landscape and vertebrate host factors

Genomic epidemiology reveals multiple introductions of SARS-CoV-2 followed by community and nosocomial spread, Germany, February to May 2020.

BackgroundIn the SARS-CoV-2 pandemic, viral genomes are available at unprecedented speed, but spatio-temporal bias in genome sequence sampling precludes phylogeographical inference without additional contextual data.AimWe applied genomic epidemiology to trace SARS-CoV-2 spread on an international, national and local level, to illustrate how transmission chains can be resolved to the level of a single event and single person using integrated sequence data and spatio-temporal metadata.MethodsWe investigated 289 COVID-19 cases at a university hospital in Munich, Germany, between 29 February and 27 May 2020. Using the ARTIC protocol, we obtained near full-length viral genomes from 174 SARS-CoV-2-positive respiratory samples. Phylogenetic analyses using the Auspice software were employed in combination with anamnestic reporting of travel history, interpersonal interactions and perceived high-risk exposures among patients and healthcare workers to characterise cluster outbreaks and establish likely scenarios and timelines of transmission.ResultsWe identified multiple independent introductions in the Munich Metropolitan Region during the first weeks of the first pandemic wave, mainly by travellers returning from popular skiing areas in the Alps. In these early weeks, the rate of presumable hospital-acquired infections among patients and in particular healthcare workers was high (9.6% and 54%, respectively) and we illustrated how transmission chains can be dissected at high resolution combining virus sequences and spatio-temporal networks of human interactions.ConclusionsEarly spread of SARS-CoV-2 in Europe was catalysed by superspreading events and regional hotspots during the winter holiday season. Genomic epidemiology can be employed to trace viral spread and inform effective containment strategies.

Phylogenomics of Northeast AsianPungitiussticklebacks

AbstractAimUnderstanding the phylogeography of a species complex can provide important insights into its evolutionary history. However, phylogeographic inference often faces the dilemma of regionally inadequate sampling.Pungitiussticklebacks are a case in point: although the highest species diversity is found in Northeast Asia, their phylogeography in this region is still poorly understood.LocationNortheast Asia.MethodsWith the aid of whole‐genome resequencing data, we investigated the phylogeography of Northeast AsianPungitiussticklebacks, with newly sampled 83 worldwidePungitiusindividuals from 11 locations including eight Chinese locations reported to host onlyP. sinensis.ResultsWe discovered that three of the eight Chinese locations hosted populations ofP. kaibarae andP. bussei, species new to the fauna of China. Phylogeographic analyses further clarified the sequence and timing of colonization of Northeast Asia by differentPungitiusspecies, shedding new light on their origins and current distribution ranges. Colonization of inland Northeast Asia byPungitiussticklebacks occurred in multiple waves, and the widespreadP. sinensisexpanded its range relatively late in the Pleistocene.Main conclusionsThis study complements our understanding of the phylogeography ofPungitiussticklebacks by extending sampling to cover an area that comprises nearly half of the known distribution area of this genus in Northeast Asia. The discovery of threePungitusspecies from China is of particular interest, as translocations to support locally declining populations have occurred under the assumption that all sticklebacks in China—except the endangeredP. stenurus—areP. sinensis, raising conservation concerns associated with unintentional translocations and possible admixture.

Emergence and spread of SARS-CoV-2 lineage B.1.620 with variant of concern-like mutations and deletions

Distinct SARS-CoV-2 lineages, discovered through various genomic surveillance initiatives, have emerged during the pandemic following unprecedented reductions in worldwide human mobility. We here describe a SARS-CoV-2 lineage - designated B.1.620 - discovered in Lithuania and carrying many mutations and deletions in the spike protein shared with widespread variants of concern (VOCs), including E484K, S477N and deletions HV69Δ, Y144Δ, and LLA241/243Δ. As well as documenting the suite of mutations this lineage carries, we also describe its potential to be resistant to neutralising antibodies, accompanying travel histories for a subset of European cases, evidence of local B.1.620 transmission in Europe with a focus on Lithuania, and significance of its prevalence in Central Africa owing to recent genome sequencing efforts there. We make a case for its likely Central African origin using advanced phylogeographic inference methodologies incorporating recorded travel histories of infected travellers.

Genomic sequencing of SARS-CoV-2 in Rwanda reveals the importance of incoming travelers on lineage diversity

COVID-19 transmission rates are often linked to locally circulating strains of SARS-CoV-2. Here we describe 203 SARS-CoV-2 whole genome sequences analyzed from strains circulating in Rwanda from May 2020 to February 2021. In particular, we report a shift in variant distribution towards the emerging sub-lineage A.23.1 that is currently dominating. Furthermore, we report the detection of the first Rwandan cases of the B.1.1.7 and B.1.351 variants of concern among incoming travelers tested at Kigali International Airport. To assess the importance of viral introductions from neighboring countries and local transmission, we exploit available individual travel history metadata to inform spatio-temporal phylogeographic inference, enabling us to take into account infections from unsampled locations. We uncover an important role of neighboring countries in seeding introductions into Rwanda, including those from which no genomic sequences were available. Our results highlight the importance of systematic genomic surveillance and regional collaborations for a durable response towards combating COVID-19.

Rivers, not refugia, drove diversification in arboreal, sub-Saharan African snakes.

The relative roles of rivers versus refugia in shaping the high levels of species diversity in tropical rainforests have been widely debated for decades. Only recently has it become possible to take an integrative approach to test predictions derived from these hypotheses using genomic sequencing and paleo‐species distribution modeling. Herein, we tested the predictions of the classic river, refuge, and river‐refuge hypotheses on diversification in the arboreal sub‐Saharan African snake genus Toxicodryas. We used dated phylogeographic inferences, population clustering analyses, demographic model selection, and paleo‐distribution modeling to conduct a phylogenomic and historical demographic analysis of this genus. Our results revealed significant population genetic structure within both Toxicodryas species, corresponding geographically to river barriers and divergence times from the mid‐Miocene to Pliocene. Our demographic analyses supported the interpretation that rivers are indications of strong barriers to gene flow among populations since their divergence. Additionally, we found no support for a major contraction of suitable habitat during the last glacial maximum, allowing us to reject both the refuge and river‐refuge hypotheses in favor of the river‐barrier hypothesis. Based on conservative interpretations of our species delimitation analyses with the Sanger and ddRAD data sets, two new cryptic species are identified from east‐central Africa. This study highlights the complexity of diversification dynamics in the African tropics and the advantages of integrative approaches to studying speciation in tropical regions.