Membracis mexicana (Hemiptera: Membracidae) is distributed in four biogeographic provinces of Mexico, field observations indicate that there are different forms of this species but so far it has not been recorded how the different forms are and how they are distributed. The aim of this study was to quantify the phenotypic and genetic variation of M. mexicana and determine if the config-uration of biogeographic provinces impacts the distribution of this variation. To do this we ana-lyzed 307 photographs using 19 landmarks and geometric morphometrics to quantify the pheno-typic variation in helmets. We sequenced five molecular markers for 205 individuals to describe the phylogeographic pattern. As a result, we identified three morphological configurations of the helmet of M. mexicana and two genetic lineages. The morphotypes are 1) large and wide helmet with small dorsal spots, 2) small and narrow helmet with large dorsal spots, and 3) small and narrow helmet with small spots. Genetic lineages are distributed at southeast and western of Mexico. The western lineage corresponds to two helmet morphotypes (1 and 2) and the southeast lineage to morphotype 3. We found that the larger helmets correspond the western lineage and are distributed in Transmexican Volcanic Belt and Pacific lowlands provinces whereas the smallest helmets correspond to the southeast lineage and are present in Veracruzan and Yucatan Peninsula provinces.

Tularemia is a severe infectious disease caused by the Gram-negative bacteria Francisella tularensis. F. tularensis is currently divided into three subspecies, holarctica, tularensis and mediasiatica which differ in their virulence and geographic distribution. Subspecies mediasiatica is the least studied because of its very low documented virulence for humans and limited geographic distribution. It was discovered in sparsely populated regions of Central Asia. Since 2011, a new subsp. mediasiatica lineage was identified in Altai (Russia). In 2021, we isolated one subsp. mediasiatica strain in Krasnoyarsk Territory. In spite of its geographic origin, 500 km east from Altai, this strain belongs to the Altai lineage and contributes surprisingly little genetic diversity to previous knowledge. This improved knowledge of the phylogeography of subsp. mediasiatica led us to propose a scenario in which the two zoonotic lineages, holarctica and tularensis, independently emerged from the mediasiatica lineage in Siberia or Central Asia, and make predictions which will allow to challenge this hypothesis

The geographic variability of the dorsal pattern (DP) of the Italian wall lizard, Podarcis siculus, across its native range was studied with the aim to understand whether the distributions of this phenotypic trait were more shaped by allopatric differentiation rather than adaptive processes. A total of 1298 georeferenced observations scattered across the Italian peninsula and the main islands (Sicily, Corsica and Sardinia) were obtained from citizen science databases and five DPs were characterized by different shapes of the dark pattern (“reticulated”, “campestris”, “reticulated/campestris” and “striped”) or by absence of it (“concolor”). Frequencies of different DP phenotypes differ between the two main mtDNA lineages settled in central-northern and in southern Italy respectively. This pattern may be indicative of a role of long-term allopatric historical processes in determining the observed pattern. The analysis also identified a putative wide area of secondary contact, in central southern Italy, characterized by high diversity of the DP. Generalized Linear Models (GLMs), used to estimate a possible association between bioclimatic variables and the observed phenotypic variation, showed that each of the five DPs is correlated to different environmental factors and show different distribution of areas with high probability of occurrence. However, for all but one of the DPs, the area with the greatest probability does not correspond exactly to the real distribution of the DP. Conversely, the “concolor” phenotype does not seem related to any particular mtDNA lineage and it shows a preference for areas with high temperature and low rainfall. This is in agreement with the expectation of low amount of melanin of the dorsal pattern that, in the study areas, is characterized by a light uniform coloration which could confer a better thermoregulation ability in high temperatures environments avoiding overheating.

Chikungunya virus (CHIKV), a mosquito-borne alphavirus of the family Togaviridae, has recently emerged in the Americas from lineages from two continents, Asia and Africa. Historically, CHIKV circulated as at least four lineages worldwide with both enzootic and epidemic transmission cycles. To understand the recent patterns of emergence and the current status of the CHIKV spread, updated analyses of the viral genetic data and metadata are needed. Here, we performed phylogenetic and comparative genomics screens of CHIKV genomes, taking advantage of the public availability of many recently sequenced isolates. Based on these new data and analyses, we derive a revised phylogeny from nucleotide sequences in coding regions. Using this phylogeny, we uncover the presence of several distinct lineages in Africa that were previously considered a single one. In parallel, we performed thermodynamic modeling of CHIKV untranslated regions (UTRs), which revealed evolutionarily conserved structured and unstructured RNA elements in the 3'UTR. We provide evidence for duplication events in recently emerged American isolates of the Asian CHIKV lineage and propose the existence of a flexible 3'UTR architecture among different CHIKV lineages.

Montiaceae (Portulacineae) comprise a clade of at least 268 species plus ca. 27 subspecific taxa primarily of western North America, the Chilean Floristic Region, and temperate Australasia. This work uses existing phylogenetic metadata to elaborate a new cladistic taxonomic synthesis. A total of 24 taxa are validated, nine new and 15 necessary nomenclatural recombinations. Hypotheses of Montiaceae historical biogeographical, ecological, and phenotypic evolution are evaluated in light of recent metadata and in terms of classical, contemporary, and novel systematic and evolutionary epistemology.

Despite constituting the western-most edge of the population distributions for several native European flora, Ireland has largely been left out of key Europe-wide phylogeographic studies. This is true for birch (Betula pubescens Ehrh. and Betula pendula Roth.), for which the genetic diversity has yet to be mapped for Ireland. Here, we used eight cpDNA (two RFLP and six SSR) markers to map the genetic diversity of B. pubescens, B. pendula and putative hybrid individuals sampled from 19 populations spread cross most of the island of Ireland. Within Ireland, 11 distinct haplotypes were detected, the most common of which (H1) was also detected in England, Scotland, France and Norway. A moderate level of population structuring (GST = 0.282) was found across Ireland and the genetic diversity of its northern populations was twice that of its southern populations. This indicates that, unlike other native Irish trees, such as oak and alder, post-glacial recolonisation by birch did not begin in the south (i.e., from Iberia). Rather, and in agreement with palynological data, birch most likely migrated in from eastern populations in Britain. Finally, we highlight Irish populations with comparatively unique genetic structure which may be included as part of European genetic conservation networks.

Asian giant hornet, Vespa mandarinia, is an invasive species that could potentially destroy local honeybee industry in North America. It has been observed to nest in the coastal regions of British Columbia in Canada and Washington State in USA. What is the source population of the immigrant hornets? The identification of the source population can not only shed light on the route of immigration, but also on the similarity between the native habitat and the potential new habitat in the Pacific Northwest. We analyzed mitochondrial COX1 sequences of specimens sampled from multiple populations in China, South Korea, Japan and the Russian Far East. V. mandarinia exhibits phylogeographic patterns, forming monophyletic clades for 16 specimens from China, six specimens from South Korea and two specimens from Japan. The two mitochondrial COX1 sequences from Nanaimo, British Columbia are identical to the two sequences from Japan. The COX1 sequence from Blaine, Washington State, clustered with those from South Korea, and is identical to one sequence from South Korea (GenBank accession MN716828). Our geophylogeny, which allows visualization of genetic variation over time and space, provides evolutionary insights on the evolution and speciation of three closely related vespine species (V. tropica, V. soror and V. mandarinia). The geophylogeny also highlights a strong insufficiency in sample collection in China. The existing sequenced specimens from mainland China with geographic coordinates were represent only four provinces, missing populations along the coastal regions in East China. One therefore cannot exclude the possibility that the Asian giant hornet found in Canada and USA may also be identical to V. mandarinia in Eastern China. We strongly recommend more extensive DNA barcoding data with geographic information and longer DNA barcode (e.g., the entire mitochondrial genome) so that critical questions on invasive species and environmental conservation can be addressed more accurately than what is possible with existing data.

Dengue virus (DENV) was detected in Senegal in 1979 for the first time. Since 2017, unprecedented frequent occurrence of DENV outbreaks was noticed yearly. In this context, epidemiological and molecular evolution data are paramount to decipher virus diffusion route. In the current study, a dengue outbreak occurred in Senegal in 2018 in the context of the largest religious gathering with 263 confirmed DENV cases out of 832 collected samples, including 25 life-threatening cases and 2 deaths. It was characterized by a co-circulation of dengue serotypes 1 and 3. Phylogenetic analysis based on the E gene revealed that the main detected serotype in Touba city was DENV-3 and belonged to Genotype III. Bayesian phylogeographic analysis was performed and suggested one viral introduction around 2017.07 (95 % HPD = 2016.61 – 2017.57) followed by cryptic circulation before identification of the first case on 01 October 2018. DENV-3 strains are phylogenetically related, with strong phylogenetic links between strains retrieved from Burkina Faso and West African countries. These phylogenetic data substantiate epidemiological data of the DENV-3 origin and spread between African countries and subsequent diffusion after religious mass events. The study also highlighted the usefulness of the Mobile Laboratory during outbreak response allowing rapid diagnosis resulting in improved patient management.

The agamid lizard Phrynocephalus melanurus, is restricted to the Northwest China (Dzungar Basin) and adjacent Eastern Kazakhstan (Zaisan and Alakol basins). To elucidate the phylogeography of P. melanurus, we obtained the mitochondrial DNA COI segments of 175 sampled lizards from 44 localities across the whole distribution range. Phylogenetic analyses revealed two main clades comprising six geographically structured lineages (I, IIa, IIb, IIc1, IIc2, and IIc3) that fits an isola-tion by distance (IBD) model. The divergence from the most recent ancestor was dated to ~1.86 million years ago (Ma). Demographic analyses demonstrated lineage-specific response to past climate change: stable population for Clade I, IIb and IIc2; while past population expansion for IIc1 and IIc3 since 0.02 Ma and 0.18 Ma, respectively. Bayesian phylogeographic diffusion anal-yses detected initial spreading at the Saur Mount vicinity approximately 1.8 Ma. Historical species distribution modeling (SDM) projected expansion of the suitable habitat in the last interglacial period and stepwise contraction in the Last Glacial Maximum and Holocene epochs. The SDM predicts a shrinkage and shift of suitable area eastward as a response to future climate change. Our findings suggest that the evolution of P. melanurus followed a parapatric divergence with subsequent dispersal and adaptation to cold and dry environments during the Quaternary. Moreover, future conservation strategies need to be taken to maintain the population of lineage IIc1.

Montiaceae comprise a clade of at least 270 species plus about 20 accepted subspecific taxa, primarily of western America and Australia. The present paper is the first of a two-part work that seeks to evaluate evolutionary theory via metadata analysis of Montiaceae. In particular, it uses metadata analysis to evaluate the theory in theory-laden methods that have been applied in evolutionary analyses of Montiaceae. This part focuses on phylogeny and phylogeography. The second part focuses on phenotypic and ecological diversification. An emergent theme in this paper is the degree to which historical idiosyncrasy during Montiaceae evolution misleads quantitative methods of evolutionary reconstruction and phylogeographic interpretation. This suggests that idiosyncraticity itself is a fundamental property of evolution. The second part of this work elaborates this notion as the Principle of Evolutionary Idiosyncraticity. The present part describes idiosyncraticity in molecular phylogenetic and phylogeographic data and uses this notion to refine ideas on Montiaceae evolution. Phylogenetic metadata conflicts and conflicting phylogeographic interpretations are discussed. I conclude that, owing to PEI, quantitative methods of evolutionary analysis cannot be globally accurate, though they are useful heuristically. In contrast, classical narrative analysis is robust in the face of PEI.

Temperature is one of the most important range-limiting factors for many seaweeds. Driven by the recent climatic changes, rapid northward shifts of species’ distribution ranges can potentially modify the phylogeographic signature of Last Glacial Maximum such as increased genetic diversity at lower latitudes. To explore this question in detail in the kelp Saccharina latissima, we used microsatellites and double digest RAD-sequencing derived SNPs on S. latissima sporophytes sampled within 11 sites spanning the entire European Atlantic latitudinal range of this species. In addition, we checked for statistical correlation between genetic marker allele frequencies and three environmental proxies (sea surface temperature, salinity, and water turbidity). Our findings revealed that genetic diversity was significantly higher for the northern localities compared to the southern ones in contrast to the expected phylogeographic pattern. This suggests that the southernmost S. latissima populations are negatively affected by the recent climatic changes but also that the recolonization of S. latissima range following the LGM may have occurred from northerly refugia. Seven SNPs and 12 microsatellite alleles were found to be significantly associated with at least one of the three environmental variables. We discuss the potential adaptive functions of the genes associated with the outlier markers and the importance of these markers for successful conservation and aquaculture strategies for S. latissima in this age of rapid global change.

Ancient lakes are among the most interesting models for evolution studies, because their biodiversity is the result of a complex combination of migration and speciation. Here, we investigate the origin of single celled planktonic eukaryotes from the oldest lake in the world – Lake Baikal. By using 18S rDNA metabarcoding we recovered 1,427 Operational Taxonomic Units (OTUs) belonging to protists populating surface waters (1-50 meters) and belonging to pico- and nano-plankton size fractions. The recovered communities resembled other lacustrine freshwater assemblages found elsewhere, especially the unclassified protists. However, our results suggest that a fraction of Baikal protists originated evolutionary recently from marine/brackish ancestors. Moreover, our results suggest that rapid radiation may have occurred among some protist taxa, partially mirroring what was already shown for multicellular organisms in Lake Baikal. We found potential species flocks in Stramenopiles, Alveolata, Opisthokonta, Archaeplastida, Rhizaria and Hacrobia. Putative flocks predominated in Chrysophytes, which are highly diverse in Lake Baikal. Some of the species, especially those from these flocks, may be endemic from Lake Baikal, because their 18S rDNA differed > 10% from known DNA. Overall, our study points to novel diversity of planktonic protists in Lake Baikal, some of which may have emerged in situ after evolutionary diversification.

Classic mitochondrial DNA phylogenetic and phylogeographic studies on the aboriginal, historical, and extant human populations of the Canary Islands have contributed to reconstructing the origin and settlement on the islands of the pre-European colonizers. The recent use of new ancient DNA targeted enrichment and next-generation sequencing techniques on new Canary Islands samples has greatly improved these molecular results, but it has also revealed significant contamination in the islands' previous aboriginal genetic pool. Following a thorough review of these cases, new phylogeographic analysis revealed the existence of a heterogeneous aboriginal population, asymmetrically distributed across the various islands, and most likely descended from a unique main settlement. These new results and new proposed coalescent ages are compatible with a Roman-mediated arrival driven by the exploitation of the purple on the Canary Islands.

The eastern subterranean termite, Reticulitermes flavipes, currently inhabits previously glaciated regions of the northeastern U.S., as well as the unglaciated southern Appalachian Mountains and surrounding areas. We hypothesized that Pleistocene climatic fluctuations have influenced the distribution of R. flavipes, and thus the evolutionary history of the species. We estimated contemporary and historical geographic distributions of R. flavipes by constructing Species Distribution Models (SDM). We also inferred the evolutionary and demographic history of the species using mitochondrial (cytochrome oxidase I and II) and nuclear (endo-beta-1,4-glucanase) DNA sequence data. To do this, genetic populations were delineated using Bayesian spatial genetic clustering, competing hypotheses about population divergence were assessed using approximate Bayesian computation (ABC), and changes in population size were estimated using Bayesian skyline plots. SDMs identified areas in the north with suitable habitat during the transition from the Last Interglacial to the Last Glacial Maximum, as well as an expanding distribution from the mid-Holocene to the present. Genetic analyses identified three geographically cohesive populations, corresponding with northern, central, and southern portions of the study region. Based on ABC analyses, divergence between the Northern and Southern populations was the oldest, estimated to have occurred 64.80 thousand years ago (kya), which corresponds with the timing of available habitat in the north. The Central and Northern populations diverged in the mid-Holocene, 8.63 kya, after which the Central population continued to expand. Accordingly, phylogeographic patterns of R. flavipes in the southern Appalachians appear to have been strongly influenced by glacial-interglacial climate change.

The field of mitochondrial genomics has advanced rapidly and revolutionized disciplines from molecular anthropology, population genetics, to medical/oncogenetics. However, mtDNA next-generation sequencing (NGS) analysis for matrilineal haplotyping and phylogeographic inference remains hindered by the lack of a consolidated, mitogenome database and efficient bioinformatics pipeline. To address this, we developed a customized human mitogenome database (hMITO DB) embedded in a CLC Genomics workflow for read mapping, variant analysis, haplotyping, and geo-mapping. The database was constructed from 4,286 mitogenomes. The macro-haplogroup (A to Z) distribution and representative phylogenetic tree were found consistent with published literature. The hMITO DB automated workflow was tested using mtDNA-NGS sequences derived from Pap smears and cervical cancer cell lines. The auto-generated read mapping, variants track, and table of haplotypes and geo-origins were completed in 15 min for 47 samples. The mtDNA workflow proved to be a rapid, efficient and accurate means of sequence analysis for translational mitogenomics.