Comparative Population Genetics Research Articles

Genetic differentiation of wild boar populations in a region affected by African swine fever

In the European Union, African swine fever (ASF) affects wild boar (Sus scrofa) populations in several Member States. Knowledge of population connectivity is important for the implementation of control measures, in particular the establishment of effective barriers. Population genetic comparisons of neighbouring populations can be very helpful in this respect. The present study investigated the genetic differentiation of wild boar in eastern Germany. This region has been affected by ASF since September 2020. A total of 1,262 wild boars from 31 hunting grounds (populations) in ASF-affected and ASF-free districts were sampled over a total area of almost 100,000 km². The study area encompassed a network of geographical factors that promote (roads, rivers, cities) or inhibit (natural areas, habitat corridors) genetic differentiation between wild boar populations. The genetic differentiation of the areas was based on 12 microsatellite markers. Three different Bayesian algorithms were used to analyse the data. The results were combined into a common approach with 9 clusters. Based on the cluster distribution in each population, the connectivity between the areas was quantified. The strongest differentiation was found along an imaginary line along the lower Elbe valley through Berlin and the A11 freeway to the Szczecin Lagoon. In contrast, the Mecklenburg Lake District and the south-east of the study area showed strong connectivity between areas. The special features of the landscapes along the lower Elbe valley, which was assessed as highly connective, and the high barrier effect of the A11 freeway in contrast to the other freeways in the study area show that barrier effects cannot be generalised in principle, but are actually determined by the circumstances of individual structures. The results of the connectivity analysis were compared with the distribution of viral lineages and variants. The genotypes of the wild boar populations and the ASFV lineages and variants showed a good explanatory approach for the observed disease dynamics in the study area. The newly gained knowledge on barriers and regionally different connectivity between wild boar populations can support considerations and measures for the containment of ASF in the affected areas by improving the understanding of wild boar dispersal dynamics.

Unravelling human migrations: DISPERSALS Project sheds light on early Homo sapiens' journey

Unravelling human migrations: DISPERSALS Project sheds light on early Homo sapiens' journey DISPERSALS will compare the human occupation and ecology between central Mozambique and eastern and southern Africa using a multi-scale approach based on the study of regional diachronic cultural traits. It will reconstruct regional population patterns, followed by comparative quantitative population genetics combined with GIS computational network analyses. The results will then be integrated through agent-based modelling, based on the incremental creation, elimination, or reorientation of network links to simulate a quantitative framework to study the evolution of population dispersal across southern-eastern Africa. The project will be crucial in providing groundbreaking high-resolution archaeological, chronological and paleoenvironmental data. DISPERSALS will deliver a fundamental perspective on the key processes that triggered migrations and dispersals within Africa and out of Africa, ultimately resulting in the human diaspora over the entire planet.

PCSK9 single nucleotide variants might affect susceptibility to malaria through regulation of CD36

Genetic polymorphism of proprotein convertase subtilisin/kexin type 9 (PCSK9), identified as single nucleotide variant (SNV), when locating at the Cys-His rich domain (CHRD), for example, E670G and I474V, is herein hypothesized to affect the domain’s binding affinity to cluster of differentiation 36 (CD36), which plays a major role in the phagocytosis of Plasmodium parasitized erythrocytes, and thus bringing on divergent susceptibility to malaria, with supporting evidence coming from epidemiologic studies of Malian children infected by malaria and population genetic comparison of European-Americans and African-Americans whose ancestors originated in the land ravaged by malaria. Starting from this proposal, the hypothesis extends further to the possible influence of PCSK9 CHRD SNVs on the development, therapeutic effects, and prognosis of cancer based on recently discovered roles of CD36 in tumour biology. It is hoped that more laboratory, epidemiology, and computation analyses in the near future can prove the whole speculation and broaden our knowledge about the physiology of PCSK9.

The maternal phylogenetic insights of Yunnan Miao group revealed by complete mitogenomes

The Miao group is one of the representative Hmong-Mien-speaking populations and primarily scattered in southern China and Southeast Asia, which has experienced massive migrations in history and thus forms distinctive evolutionary genetics. Yet, the genetic explorations of Miao group are relatively limited based on complete mitochondrial genome (mitogenome), especially for the Miao group from Yunnan Province (YNM). Here, we sequenced complete mitogenomes of 132 Miao individuals from Yunnan Province using massively parallel sequencing method. Total 132 Miao individuals could be allocated to 119 various haplotypes, which were mainly dominated by haplogroups prevalent in southern East Asia (B, F, M7 and R9), and rarely occupied by northern lineages (A, D, G and M8). In order to dissect the genetic background of YNM more comprehensively, we introduced 99 published population data with 7135 complete mitochondrial sequences for population genetic comparisons. YNM exhibited closer genetic relationships with Hmong-Mien, Tai-Kadai, Sino-Tibetan and Austroasiatic populations, especially for Hmong-Mien populations; we further speculated that Miao group might have certain direct or indirect gene exchanges with ancient Baiyue groups. Several maternal lineages, such as B5a1c1a, F1g1, B4a5 and D4e1a3, were found to be specifically shared by YNM and other Hmong-Mien populations, and these matrilineal expansions occurred roughly during the Neolithic period. Eventually, according to the population dynamic analyses of YNM, the population size began to emerge recovery ∼1–0.5 kya after a long-term population reduction ∼1–5 kya, during which the B5a1c1a haplogroup manifested relatively apparent lineage expansion.

Genetic diversity and microevolution in clinical Cryptococcus isolates from Cameroon.

Cryptococcal meningitis is the second most common cause of death in people living with HIV/AIDS, yet we have a limited understanding of how cryptococcal isolates change over the course of infection. Cryptococcal infections are environmentally acquired, and the genetic diversity of these infecting isolates can also be geographically linked. Here, we employ whole genome sequences for 372 clinical Cryptococcus isolates from 341 patients with HIV-associated cryptococcal meningitis obtained via a large clinical trial, across both Malawi and Cameroon, to enable population genetic comparisons of isolates between countries. We see that isolates from Cameroon are highly clonal, when compared to those from Malawi, with differential rates of disruptive variants in genes with roles in DNA binding and energy use. For a subset of patients (22) from Cameroon, we leverage longitudinal sampling, with samples taken at days 7 and 14 post-enrollment, to interrogate the genetic changes that arise over the course of infection, and the genetic diversity of isolates within patients. We see disruptive variants arising over the course of infection in several genes, including the phagocytosis-regulating transcription factor GAT204. In addition, in 13% of patients sampled longitudinally, we see evidence for mixed infections. This approach identifies geographically linked genetic variation, signatures of microevolution, and evidence for mixed infections across a clinical cohort of patients affected by cryptococcal meningitis in Central Africa.

Population genetic analyses reveal host association and genetically distinct populations of social parasite Solenopsis daguerrei (Hymenoptera: Formicidae)

Inquiline ant social parasites exploit other ant species for their reproductive benefit because they do not possess a worker caste. Due to their relative rarity in nature, the biology and natural history of inquilines are largely unknown. Likewise, not much research exists that details the close relationship between inquilines and their host(s), and how each organism influences the genetic structure of the other. Here, we conducted a comparative population genetics study to assess patterns of genetic structure within and among populations of inquiline Solenopsis daguerrei and its known fire ant hosts, which includes invasive Solenopsis invicta. Using nuclear and mitochondrial markers, we show that four genetically distinct groups of S. daguerrei likely exist, each with different degrees of host association. Consistent with previous inferences of the inquiline lifestyle, we find that inbreeding is common in S. daguerrei, presumably a result of intranidal mating and restricted dispersal. Results from this study, specifically host association patterns, may inform future biological control strategies to mitigate invasive S. invicta populations.

Comparative population genetics of habitat-forming octocorals in two marine protected areas: eco-evolutionary and management implications

Current efforts to halt the decline of biodiversity are based primarily on protecting species richness. This narrow focus overlooks key components of biological diversity, particularly the infra-species genetic diversity, which is critical to consider with respect to genetic adaptation in changing environments. While comparative population genetics is recognized as a relevant approach to improve biodiversity management, it is still barely considered in practice. Here, a comparative population genetics study was conducted on two key habitat-forming octocoral species, Corallium rubrum and Paramuricea clavata, to contribute to management of two Marine Protected Areas (MPAs) in the northwestern Mediterranean. Contrasting patterns of genetic diversity and structure were observed in the two species, although they share many common biological features and live in similar habitats. Differential genetic drift effects induced by species-specific reproductive strategies and demographic histories most likely explain these differences. The translation of our results into management strategies supports the definition of four management units. We identified a coldspot of genetic diversity, with genetically isolated populations, and a hotspot of genetic diversity that has a central role in the system’s connectivity. Interestingly, they corresponded to the most recent and the oldest protected areas, respectively. This case study shows how moving from a “species pattern” perspective to an “eco-evolutionary processes” perspective can help assess and contribute to the effectiveness of biodiversity management plans.

Comparative Phylogeography, Historical Demography, and Population Genetics of Three Common Coastal Fauna in Spartina Marshes of the Northwestern Gulf of Mexico

Coastal wetlands worldwide are experiencing high rates of loss and degradation that may lead to a reduction in diversity in faunal populations. Since salt marsh habitats are subject to a multitude of stressors, evaluations of the genetic diversity, connectivity, and potential resilience of faunal communities within salt marsh habitats are relevant. This study characterizes mitochondrial DNA (mtDNA) diversity for three common faunal residents of salt marshes along the northern Gulf of Mexico. Gulf Killifish (Fundulus grandis) samples were characterized for 1077 bp of the concatenated nucleotide sequence corresponding to the Control Region and Nitrogen Dehydrogenase, Subunits 2 and 5. Daggerblade grass shrimp (Palaemon pugio) samples were characterized using 466 bp of 16sRNA sequence, and phloem-feeding planthoppers (Prokelisia marginata) were characterized using 372 bp of Cytochrome c Oxidase Subunit I (COI) sequence. For F. grandis, our data revealed high levels of haplotypic diversity, evidence of isolation by distance (IBD), and regional population structuring associated with the distribution of two distinct phylogroups and distinct historical demography signatures. P. pugio and P. marginata displayed low levels of haplotypic diversity and evidence of population structure, but both appear to contain only snapshots of the total potential diversity for these species in the Gulf of Mexico. Greater resolution of the patterns of historical demography of Gulf Killifish may be obtained in future studies by including localities from Florida and Mexico. For both P. pugio and planthoppers, future studies would benefit from the characterization of genetic markers with a higher degree of polymorphism. We conclude that despite these three species inhabiting the same habitats along the same stretch of coast, each is subject to a different combination of evolutionary forces, and this study was able to reconstruct differences in how the genetic variation in each of these species emerged, and how it is maintained.

Evolution and biogeography of apple stem grooving virus

BackgroundApple stem grooving virus (ASGV) has a wide host range, notably including apples, pears, prunes and citrus. It is found worldwide.MethodIn this study, two near complete genomes, and seven coat protein (CP) sequences of Iranian isolates from apple were determined. Sequences added from GenBank provided alignments of 120 genomic sequences (54 of which were recombinant), and 276 coat protein genes (none of them recombinant).ResultThe non-recombinant genomes gave a well supported phylogeny with isolates from diverse hosts in China forming the base of the phylogeny, and a monophyletic clade of at least seven clusters of isolates from around the world with no host or provenace groupings among them, and all but one including isolates from China. The six regions of the ASGV genome (five in one frame, one − 2 overlapping) gave significantly correlated phylogenies, but individually had less statistical support. The largest cluster of isolates contained those from Iran and had isolates with worldwide provenances, and came from a wide range of mono- and dicotyledonous hosts. Population genetic comparisons of the six regions of the ASGV genome showed that four were under strong negative selection, but two of unknown function were under positive selection.ConclusionASGV most likely originated and spread in East Asia in one or more of various plant species, but not in Eurasia; the ASGV population of China had the greatest overall nucleotide diversity and largest number of segregating sites.

Comparative population genetics of congeneric limpets across a biogeographic transition zone reveals common patterns of genetic structure and demographic history.

The distribution of genetic diversity is often heterogeneous in space, and it usually correlates with environmental transitions or historical processes that affect demography. The coast of Chile encompasses two biogeographic provinces and spans a broad environmental gradient together with oceanographic processes linked to coastal topography that can affect species' genetic diversity. Here, we evaluated the genetic connectivity and historical demography of four Scurria limpets, S. scurra, S. variabilis, S. ceciliana and S. araucana, between ca. 19° S and 53° S in the Chilean coast using genome-wide SNPs markers. Genetic structure varied among species which was evidenced by species-specific breaks together with two shared breaks. One of the shared breaks was located at 22-25° S and was observed in S. araucana and S. variabilis, while the second break around 31-34° S was shared by three Scurria species. Interestingly, the identified genetic breaks are also shared with other low-disperser invertebrates. Demographic histories show bottlenecks in S. scurra and S. araucana populations and recent population expansion in all species. The shared genetic breaks can be linked to oceanographic features acting as soft barriers to dispersal and also to historical climate, evidencing the utility of comparing multiple and sympatric species to understand the influence of a particular seascape on genetic diversity.

Emergence and clonal expansion of Vibrio aestuarianus lineages pathogenic for oysters in Europe.

Crassostrea gigas oysters represent a significant global food source with 4.7 million tons harvested per year. In 2001, the bacterium Vibrio aestuarianus subsp. francensis emerged as a pathogen that causes adult oyster mortality in France and Ireland. Its impact on oyster aquaculture has increased in Europe since its re-emergence in 2012. To better understand the evolutionary mechanisms leading to the emergence and persistence over time of this pathogen, we conducted a survey of mollusc diseases through national reference laboratories across Europe. We analysed 54 new genomes of Vibrio aestuarianus (Va) isolated from multiple environmental compartments since 2001, in areas with and without bivalve mortalities. We used a combination of comparative genomics and population genetics approaches and show that Va has a classical epidemic population structure from which the pathogenic Va francensis subspecies emerged and clonally expanded. Furthermore, we identified a specific cus-cop-containing island conferring copper resistance to Va francensis whose acquisition may have favoured the emergence of pathogenic lineages adapted and specialized to oysters.

Numt Parser: Automated identification and removal of nuclear mitochondrial pseudogenes (numts) for accurate mitochondrial genome reconstruction in Panthera.

Nuclear mitochondrial pseudogenes (numts) may hinder the reconstruction of mtDNA genomes and affect the reliability of mtDNA datasets for phylogenetic and population genetic comparisons. Here, we present the program Numt Parser, which allows for the identification of DNA sequences that likely originate from numt pseudogene DNA. Sequencing reads are classified as originating from either numt or true cytoplasmic mitochondrial (cymt) DNA by direct comparison against cymt and numt reference sequences. Classified reads can then be parsed into cymt or numt datasets. We tested this program using whole genome shotgun-sequenced data from 2 ancient Cape lions (Panthera leo), because mtDNA is often the marker of choice for ancient DNA studies and the genus Panthera is known to have numt pseudogenes. Numt Parser decreased sequence disagreements that were likely due to numt pseudogene contamination and equalized read coverage across the mitogenome by removing reads that likely originated from numts. We compared the efficacy of Numt Parser to 2 other bioinformatic approaches that can be used to account for numt contamination. We found that Numt Parser outperformed approaches that rely only on read alignment or Basic Local Alignment Search Tool (BLAST) properties, and was effective at identifying sequences that likely originated from numts while having minimal impacts on the recovery of cymt reads. Numt Parser therefore improves the reconstruction of true mitogenomes, allowing for more accurate and robust biological inferences.

Population genetic analyses inferred a limited genetic diversity across the pvama-1 DI domain among Plasmodium vivax isolates from Khyber Pakhtunkhwa regions of Pakistan

BackgroundPlasmodium vivax apical membrane antigen-1 (pvama-1) is an important vaccine candidate against Malaria. The genetic composition assessment of pvama-1 from wide-range geography is vital to plan the antigen based vaccine designing against Malaria.MethodsThe blood samples were collected from 84 P. vivax positive malaria patients from different districts of Khyber Pakhtunkhwa (KP) province of Pakistan. The highly polymorphic and immunogenic domain-I (DI) region of pvama-1 was PCR amplified and DNA sequenced. The QC based sequences raw data filtration was done using DNASTAR package. The downstream population genetic analyses were performed using MEGA4, DnaSP, Arlequin v3.5 and Network.5 resources.ResultsThe analyses unveiled total 57 haplotypes of pvama-1 (DI) in KP samples with majorly prevalent H-14 and H-5 haplotypes. Pairwise comparative population genetics analyses identified limited to moderate genetic distinctions among the samples collected from different districts of KP, Pakistan. In context of worldwide available data, the KP samples depicted major genetic differentiation against the Korean samples with Fst = 0.40915 (P-value = 0.0001), while least distinction was observed against Indian and Iranian samples. The statistically significant negative values of Fu and Li’s D* and F* tests indicate the evidence of population expansion and directional positive selection signature. The slow LD decay across the nucleotide distance in KP isolates indicates low nucleotide diversity. In context of reference pvama-1 sequence, the KP samples were identified to have 09 novel non-synonymous single nucleotide polymorphisms (nsSNPs), including several trimorphic and tetramorphic substitutions. Few of these nsSNPs are mapped within the B-cell predicted epitopic motifs of the pvama-1, and possibly modulate the immune response mechanism.ConclusionLow genetic differentiation was observed across the pvama-1 DI among the P. vivax isolates acquired from widespread regions of KP province of Pakistan. The information may implicate in future vaccine designing strategies based on antigenic features of pvama-1.

Mitochondrial cytochrome b indicates the presence of two paraphyletic diverged lineages of the blue sheep Pseudois nayaur across the Indian Himalaya: conservation implications.

Populations exhibit signatures of local adaptive traits due to spatial and environmental heterogeneity resulting in microevolution. The blue sheep is widely distributed across the high Asian mountains and are the snow leopard's principal prey species. These mountains differ in their evolutionary history due to differential glaciation and deglaciation periods, orography, and rainfall patterns, and such factors causes diversification in species. Therefore, we assess the phylogeographic status of blue sheep using the mitochondrial cytochrome b gene (220bp) across the Indian Himalayan region (IHR) and its relationship with other populations. Of the observed five haplotypes, two and three were from the western Himalayas (WH) and eastern Himalayas (EH) respectively. One of the haplotypes from WH was shared with the population of Pamir plateau, suggesting historical maternal connectivity between these areas. The phylogenetic analyses split the blue sheep into two paraphyletic clades, and western and eastern populations of IHR were within the Pamir and Tibetan plateau clades, respectively. We observed a relatively higher mean sequence divergence in the EH population than in the WH. We propose five 'Evolutionary Significant Units' across the blue sheep distribution range based on observed variation in the species' ecological requirements, orography, climatic conditions, and maternal lineages, viz.; Western Himalaya-Pamir plateau (WHPP); Eastern Himalaya-Tibetan plateau (EHTP); Qilian mountains; Helan mountains and Hengduan mountains population. Despite the small sample size, population divergence was observed across the IHR, therefore, we suggest a transboundary, collaborative study on comparative morphology, anatomy, ecology, behaviour, and population genetics using harmonized different genetic markers for identifying the overall taxonomic status of the blue sheep across its range for planning effective conservation strategies.

Can sea snakes slither through seascape structure? Comparative phylogeography and population genetics of Hydrophis group sea snakes in Australia and Southeast Asia

Pleistocene sea level changes substantially shaped the biogeography of northern Australia and the Indo-Malayan Archipelago (IMA). For co-distributed species, their phylogeographic and population genetic patterns are expected to be concomitant with geological transformations of the Pleistocene. However, species-specific ecologies and life history traits may also be influential in generating patterns which depart from simple expectations arising from biogeographic features. Thus, comparative population genetic studies, which use taxa that reduces variation in taxonomy and geography, may refine our understanding of how biogeographic elements shape the populations of co-occurring species. Here, we sampled two sea snake species, Hydrophis curtus and H. elegans, throughout their known ranges in the IMA and northern Australia. These sea snakes have similar life history strategies and ecologies as well as overlapping distributions across the Torres Strait, a well-known biogeographic feature. We analysed two mitochondrial DNA (mtDNA) fragments and 10 microsatellite loci using traditional population genetic approaches and used Bayesian clustering methods to examine species- specific phylogenetic relationships, genetic diversities, and population genetic structures. For both species, we found a consistent lack of significant genetic variation among sampling sites across the Gulf of Carpentaria (GOC) and the Great Barrier Reef (GBR). Similarly, Bayesian clustering showed no to weak genetic partitioning across the historical Torres Strait land bridge. Both species sampled in Australia displayed population expansion signatures in tests using mtDNA and microsatellite markers. We conclude that the phylogeographic and population genetic patterns of these sea snake species do not align with the Torres Strait land bridge. This lack of population genetic structure departs from previous findings on Aipysurus sea snakes and may be linked to the association of Hydrophis species to soft sediment habitats typically found across northern Australia. These divergent patterns between the sea snake groups present the importance of considering taxon-specific attributes in formulating conservation strategies.

A genomic can of worms for schistosome host-specificity

Understanding the genetic underpinnings of schistosome host preferences is critical. Luo et al. recently identified genes associated with intermediate and definitive host-switching based on a new chromosome-level genome for Schistosoma japonicum, population genetic comparisons, and follow-up experiments. This represents a guide to fully map-selected schistosome genes using population genetics.

Covariance of pairwise differences on a multi-species coalescent tree and implications for FST.

The multi-species coalescent (MSC) provides a theoretical foundation for modern phylogenetics and comparative population genetics. Its theoretical properties have been heavily studied but there are still aspects of the MSC that are largely unknown, including the covariances in pairwise coalescence times, which are fundamental for understanding the properties of statistics that combine data from multiple species, such as the fixation index (FST). The major contribution of this study is the derivation and implementation of exact expressions for the covariances of pairwise coalescence times under phylogenetic models with piecewise constant changes in population size, assuming no gene flow after species divergence. We use these expressions to derive the variance in average pairwise differences within and between populations. We then derive approximations for the expectation and bias of a sequence-based estimator of FST, a commonly used genetic measurement of population differentiation, when it is applied to a non-recombining region of the genome. We show that the estimator of FST is generally biased downward. A freely available software package is provided, STCov, to calculate the mean, variances and covariances in coalescence times presented here under user-defined piecewise-constant species trees.This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.

High-resolution genomic surveillance elucidates a multilayered hierarchical transfer of resistance between WWTP- and human/animal-associated bacteria

BackgroundOur interconnected world and the ability of bacteria to quickly swap antibiotic resistance genes (ARGs) make it particularly important to establish the epidemiological links of multidrug resistance (MDR) transfer between wastewater treatment plant (WWTP)- and human/animal-associated bacteria, under the One Health framework. However, evidence of ARGs exchange and potential factors that contribute to this transfer remain limited.ResultsHere, by combining culture-based population genomics and genetic comparisons with publicly available datasets, we reconstructed the complete genomes of 82 multidrug-resistant isolates from WWTPs and found that most WWTP-associated isolates were genetically distinct from their closest human/animal-associated relatives currently available in the public database. Even in the minority of lineages that were closely related, WWTP-associated isolates were characterized by quite different plasmid compositions. We identified a high diversity of circular plasmids (264 in total, of which 141 were potentially novel), which served as the main source of resistance, and showed potential horizontal transfer of ARG-bearing plasmids between WWTP- and humans/animal-associated bacteria. Notably, the potentially transferred ARGs and virulence factors (VFs) with different genetic backgrounds were closely associated with flanking insertion sequences (ISs), suggesting the importance of synergy between plasmids and ISs in mediating a multilayered hierarchical transfer of MDR and potentiating the emergence of MDR-hypervirulent clones.ConclusionOur findings advance the current efforts to establish potential epidemiological links of MDR transmission between WWTP- and human/animal-associated bacteria. Plasmids play an important role in mediating the transfer of ARGs and the IS-associated ARGs that are carried by conjugative plasmids should be prioritized to tackle the spread of resistance.FH5Vi-SGsgm43KEHZNerU2Video

Circulation and Evolution of SARS-CoV-2 in India: Let the Data Speak.

The COVID-19 pandemic is a global challenge that impacted 200+ countries. India ranks in the second and third positions in terms of number of reported cases and deaths. Being a populous country with densely packed cities, SARS-CoV-2 spread exponentially. India sequenced ≈0.14% isolates from confirmed cases for pandemic surveillance and contributed ≈1.58% of complete genomes sequenced globally. This study was designed to map the circulating lineage diversity and to understand the evolution of SARS-CoV-2 in India using comparative genomics and population genetics approaches. Despite varied sequencing coverage across Indian States and Union Territories, isolates belonging to variants of concern (VoC) and variants of interest (VoI) circulated, persisted, and diversified during the first seventeen months of the pandemic. Delta and Kappa lineages emerged in India and spread globally. The phylogenetic tree shows lineage-wise monophyletic clusters of VoCs/VoIs and diversified tree topologies for non-VoC/VoI lineages designated as ‘Others’ in this study. Evolutionary dynamics analyses substantiate a lack of spatio-temporal clustering, which is indicative of multiple global and local introductions. Sites under positive selection and significant variations in spike protein corroborate with the constellation of mutations to be monitored for VoC/VoI as well as substitutions that are characteristic of functions with implications in virus–host interactions, differential glycosylation, immune evasion, and escape from neutralization.

Comparative population genetics of the invasive mosquito Aedes albopictus and the native mosquito Aedes flavopictus in the Korean peninsula

BackgroundAedes mosquitoes are important invasive species contributing to the spread of chikungunya, dengue fever, yellow fever, zika virus, and other dangerous vector-borne diseases. Aedes albopictus is native to southeast Asia, with rapid expansion due to human activity, showing a wide distribution in the Korean peninsula. Aedes flavopictus is considered to be native to East Asia, with a broad distribution in the region, including the Korean peninsula. A better understanding of the genetic diversity of these species is critical for establishing strategies for disease prevention and vector control.MethodsWe obtained DNA from 148 specimens of Ae. albopictus and 166 specimens of Ae. flavopictus in Korea, and amplified two mitochondrial genes (COI and ND5) to compare the genetic diversity and structure of the two species.ResultsWe obtained a 658-bp sequence of COI and a 423-bp sequence of ND5 from both mosquito species. We found low diversity and a nonsignificant population genetic structure in Ae. albopictus, and high diversity and a nonsignificant structure in Ae. flavopictus for these two mitochondrial genes. Aedes albopictus had fewer haplotypes with respect to the number of individuals, and a slight mismatch distribution was confirmed. By contrast, Ae. flavopictus had a large number of haplotypes compared with the number of individuals, and a large unimodal-type mismatch distribution was confirmed. Although the genetic structure of both species was nonsignificant, Ae. flavopictus exhibited higher genetic diversity than Ae. albopictus.ConclusionsAedes albopictus appears to be an introduced species, whereas Ae. flavopictus is endemic to the Korean peninsula, and the difference in genetic diversity between the two species is related to their adaptability and introduction history. Further studies on the genetic structure and diversity of these mosquitos will provide useful data for vector control.