Continent-wide genomic analysis of the African buffalo (Syncerus caffer).

The African buffalo (Syncerus caffer) is a wild bovid with a historical distribution across much of sub-Saharan Africa. Genomic analysis can provide insights into the evolutionary history of the species, and the key selective pressures shaping populations, including assessment of population level differentiation, population fragmentation, and population genetic structure. In this study we generated the highest quality de novo genome assembly (2.65 Gb, scaffold N50 69.17 Mb) of African buffalo to date, and sequenced a further 195 genomes from across the species distribution. Principal component and admixture analyses provided little support for the currently described four subspecies. Estimating Effective Migration Surfaces analysis suggested that geographical barriers have played a significant role in shaping gene flow and the population structure. Estimated effective population sizes indicated a substantial drop occurring in all populations 5-10,000 years ago, coinciding with the increase in human populations. Finally, signatures of selection were enriched for key genes associated with the immune response, suggesting infectious disease exert a substantial selective pressure upon the African buffalo. These findings have important implications for understanding bovid evolution, buffalo conservation and population management.

Recurrent gene flow between Neanderthals and modern humans over the past 200,000 years.

Although it is well known that the ancestors of modern humans and Neanderthals admixed, the effects of gene flow on the Neanderthal genome are not well understood. We develop methods to estimate the amount of human-introgressed sequences in Neanderthals and apply it to whole-genome sequence data from 2000 modern humans and three Neanderthals. We estimate that Neanderthals have 2.5 to 3.7% human ancestry, and we leverage human-introgressed sequences in Neanderthals to revise estimates of Neanderthal ancestry in modern humans, show that Neanderthal population sizes were significantly smaller than previously estimated, and identify two distinct waves of modern human gene flow into Neanderthals. Our data provide insights into the genetic legacy of recurrent gene flow between modern humans and Neanderthals.

The first complete genome of the extinct European wild ass (Equus hemionus hydruntinus).

We present palaeogenomes of three morphologically unidentified Anatolian equids dating to the first millennium BCE, sequenced to a coverage of 0.6-6.4×. Mitochondrial DNA haplotypes of the Anatolian individuals clustered with those of Equus hydruntinus (or Equus hemionus hydruntinus), the extinct European wild ass, secular name 'hydruntine'. Further, the Anatolian wild ass whole genome profiles fell outside the genomic diversity of other extant and past Asiatic wild ass (E. hemionus) lineages. These observations suggest that the three Anatolian wild asses represent hydruntines, making them the latest recorded survivors of this lineage, about a millennium later than the latest observations in the zooarchaeological record. Our mitogenomic and genomic analyses indicate that E. h. hydruntinus was a clade belonging to ancient and present-day E. hemionus lineages that radiated possibly between 0.6 and 0.8 Mya. We also find evidence consistent with recent gene flow between hydruntines and Middle Eastern wild asses. Analyses of genome-wide heterozygosity and runs of homozygosity suggest that the Anatolian wild ass population may have lost genetic diversity by the mid-first millennium BCE, a possible sign of its eventual demise.

Impact of patrilocality on contrasting patterns of paternal and maternal heritage in Central-West Africa.

Despite their ancient past and high diversity, African populations are the least represented in human population genetic studies. In this study, uniparental markers (mtDNA and Y chromosome) were used to investigate the impact of sociocultural factors on the genetic diversity and inter-ethnolinguistic gene flow in the three major Nigerian groups: Hausa (n = 89), Yoruba (n = 135) and Igbo (n = 134). The results show a distinct history from the maternal and paternal perspectives. The three Nigerian groups present a similar substrate for mtDNA, but not for the Y chromosome. The two Niger-Congo groups, Yoruba and Igbo, are paternally genetically correlated with populations from the same ethnolinguistic affiliation. Meanwhile, the Hausa is paternally closer to other Afro-Asiatic populations and presented a high diversity of lineages from across Africa. When expanding the analyses to other African populations, it is observed that language did not act as a major barrier to female-mediated gene flow and that the differentiation of paternal lineages is better correlated with linguistic than geographic distances. The results obtained demonstrate the impact of patrilocality, a common and well-established practice in populations from Central-West Africa, in the preservation of the patrilineage gene pool and in the affirmation of identity between groups.

Genetic analyses and dispersal patterns unveil the Amazonian origin of guava domestication.

Guava (Psidium guajava L.) is a semi-domesticated fruit tree of moderate importance in the Neotropics, utilized for millennia due to its nutritional and medicinal benefits, but its origin of domestication remains unknown. In this study, we examine genetic diversity and population structure in 215 plants from 11 countries in Mesoamerica, the Andes, and Amazonia using 25 nuclear microsatellite loci to propose an origin of domestication. Genetic analyses reveal one gene pool in Mesoamerica (Mexico) and four in South America (Brazilian Amazonia, Peruvian Amazonia and Andes, and Colombia), indicating greater differentiation among localities, possibly due to isolation between guava populations, particularly in the Amazonian and Andean regions. Moreover, Mesoamerican populations show high genetic diversity, with moderate genetic structure due to gene flow from northern South American populations. Dispersal scenarios suggest that Brazilian Amazonia is the probable origin of guava domestication, spreading from there to the Peruvian Andes, northern South America, Central America, and Mexico. These findings present the first evidence of guava domestication in the Americas, contributing to a deeper understanding of its evolutionary history.

Coastal dunes as drivers of genetic differentiation in the honeypot ant Myrmecocystus baja (Formicidae: Formicinae).

Coastal dunes are unique habitats, threatened by human activities. In biogeographical terms, coastal dunes are habitat islands, being discrete and distinct patches of similar habitat among themselves, separated from each other by a different type of habitat. Furthermore, coastal dunes harbor endemic species, adapted to living solely in the habitats found on specific dune systems. For example, the honeypot ant Myrmecocystus baja is endemic and restricted to coastal dunes of Mexico's Baja California Pacific coast. This ecological and biogeographical scenario led to the questions whether their geographical isolation is reflected in their genetic diversity and structuring, and how their demographic history is related with the formation of the dune system habitats. To answer these questions, population genetic, isolation-with-migration, and phylogeographical analyses were carried out, based on mitochondrial and five nuclear intronic markers. Minimal gene flow was detected only between two of the dune systems sampled; otherwise, the M. baja populations were found to be isolated and genetically structured, and their divergence generally pre-dated the modern-day dune systems. It is therefore highly likely that these ants were already present in paleodunes and that each of the populations was established from founder populations as the dunes formed. These findings highlight the importance of coastal dunes for species such as the honeypot ant from Baja California, in promoting genetic differentiation.

Genetic diversity, determinants, and dissemination of Burkholderia pseudomallei lineages implicated in melioidosis in Northeast Thailand.

Melioidosis is an often-fatal neglected tropical disease caused by an environmental bacterium Burkholderia pseudomallei. However, our understanding of the disease-causing bacterial lineages, their dissemination, and adaptive mechanisms remains limited. To address this, we conduct a comprehensive genomic analysis of 1,391 B. pseudomallei isolates collected from nine hospitals in northeast Thailand between 2015 and 2018, and contemporaneous isolates from neighbouring countries, representing the most densely sampled collection to date. Our study identifies three dominant lineages, each with unique gene sets potentially enhancing bacterial fitness in the environment. We find that recombination drives lineage-specific gene flow. Transcriptome analyses of representative clinical isolates from each dominant lineage reveal increased expression of lineage-specific genes under environmental conditions in two out of three lineages. This underscores the potential importance of environmental persistence for these dominant lineages. The study also highlights the influence of environmental factors such as terrain slope, altitude, and river direction on the geographical dispersal of B. pseudomallei. Collectively, our findings suggest that environmental persistence may play a role in facilitating the spread of B. pseudomallei, and as a prerequisite for exposure and infection, thereby providing useful insights for informing melioidosis prevention and control strategies.

Southern marsh deer (Blastocerus dichotomus) populations assessed using Amplicon Sequencing on fecal samples.

Populations in isolated and small fragments lose genetic variability very fast and are usually of conservation concern because they are at greater risk of local extinction. The largest native deer in South America, Blastocerus dichotomus (Illiger, 1815), is a Vulnerable species according to the IUCN categorization, which inhabits tropical and subtropical swampy areas. In Argentina, its presence has been restricted to four isolated fragments. Here we examine the genetic diversity and differentiation among three of them, including the three different patches that form the southernmost population, using 18 microsatellite markers genotyped by Amplicon Sequencing of DNA extracted from fecal samples. Genetic diversity was low (HE < 0.45) in all three populations studied. We found three genetic clusters compatible with the geographic location of the samples. We also found a metapopulation dynamics that involves the patches that make up the southernmost population, with evidence of a barrier to gene flow between two of them. Our results point to the creation of a corridor as a necessary and urgent management action. This is the first study, at the population level, employing microsatellite genotyping by Amplicon Sequencing with non-invasive samples in an endangered species.

Plasmodium vivax populations in the western Greater Mekong Subregion evaluated using a genetic barcode.

An improved understanding of the Plasmodium vivax populations in the Great Mekong Subregion (GMS) is needed to monitor the progress of malaria elimination. This study aimed to use a P. vivax single nucleotide polymorphism (SNP) barcode to evaluate the population dynamics and explore the gene flow among P. vivax parasite populations in the western GMS (China, Myanmar and Thailand). A total of 315 P. vivax patient samples collected in 2011 and 2018 from four regions of the western GMS were genotyped for 42 SNPs using the high-throughput MassARRAY SNP genotyping technology. Population genetic analysis was conducted to estimate the genetic diversity, effective population size, and population structure among the P. vivax populations. Overall, 291 samples were successfully genotyped at 39 SNPs. A significant difference was observed in the proportion of polyclonal infections among the five P. vivax populations (P = 0.0012, Pearson Chi-square test, χ2 = 18.1), with western Myanmar having the highest proportion (96.2%, 50/52) in 2018. Likewise, the average complexity of infection was also highest in western Myanmar (1.31) and lowest in northeast Myanmar (1.01) in 2018. The older samples from western China in 2011 had the highest pairwise nucleotide diversity (π, 0.388 ± 0.046), expected heterozygosity (He, 0.363 ± 0.02), and the largest effective population size. In comparison, in the neighboring northeast Myanmar, the more recent samples in 2018 showed the lowest values (π, 0.224 ± 0.036; He, 0.220 ± 0.026). Furthermore, the 2018 northeast Myanmar parasites showed high and moderate genetic differentiation from other populations with FST values of 0.162-0.252, whereas genetic differentiation among other populations was relatively low (FST ≤ 0.059). Principal component analysis, phylogeny, and STRUCTURE analysis showed that the P. vivax population in northeast Myanmar in 2018 substantially diverged from other populations. Although the 42 SNP barcode is a valuable tool for tracking parasite origins of worldwide parasite populations, a more extended barcode with additional SNPs is needed to distinguish the more related parasite populations in the western GMS.

Genomics of hybrid parallel origin in Aquilegia ecalcarata.

BACKGROUND: The parallel evolution of similar traits or species provides strong evidence for the role of natural selection in evolution. Traits or species that evolved repeatedly can be driven by separate de novo mutations or interspecific gene flow. Although parallel evolution has been reported in many studies, documented cases of parallel evolution caused by gene flow are scarce by comparison. Aquilegia ecalcarata and A. kansuensis belong to the genus of Aquilegia, and are the closest related sister species. Mutiple origins of A. ecalcarata have been reported in previous studies, but whether they have been driven by separate de novo mutations or gene flow remains unclear. RESULTS: In this study, We conducted genomic analysis from 158 individuals of two repeatedly evolving pairs of A. ecalcarata and A. kansuensis. All samples were divided into two distinct clades with obvious geographical distribution based on phylogeny and population structure. Demographic modeling revealed that the origin of the A. ecalcarata in the Eastern of China was caused by gene flow, and the Eastern A. ecalcarata occurred following introgression from Western A. ecalcarata population. Analysis of Treemix and D-statistic also revealed that a strong signal of gene flow was detected from Western A. ecalcarata to Eastern A. ecalcarata. Genetic divergence and selective sweep analyses inferred parallel regions of genomic divergence and identified many candidate genes associated with ecologically adaptive divergence between species pair. Comparative analysis of parallel diverged regions and gene introgression confirms that gene flow contributed to the parallel evolution of A. ecalcarata. CONCLUSIONS: Our results further confirmed the multiple origins of A. ecalcarata and highlighted the roles of gene flow. These findings provide new evidence for parallel origin after hybridization as well as insights into the ecological adaptation mechanisms underlying the parallel origins of species.

Population structure and diversification of Gymnospermium kiangnanense, a plant species with extremely small populations endemic to eastern China.

Background: Gymnospermium kiangnanense is the only species distributed in the subtropical region within the spring ephemeral genus Gymnospermium. Extensive human exploitation and habitat destruction have resulted in a rapid shrink of G. kiangnanense populations. This study utilizes microsatellite markers to analyze the genetic diversity and structure and to deduce historical population events of extant populations of G. kiangnanense. Methods: A total of 143 individuals from eight extant populations of G. kiangnanense, including two populations from Anhui Province and six populations from Zhejiang Province, were analyzed with using 21 pairs of microsatellite markers. Genetic diversity indices were calculated using Cervus, GENEPOP, GenALEX. Population structure was assessed using genetic distance (UPGMA), principal coordinate analysis (PCoA), Bayesian clustering method (STRUCTURE), and molecular variation analysis of variance (AMOVA). Population history events were inferred using DIYABC. Results: The studied populations of G. kiangnanense exhibited a low level of genetic diversity (He = 0.179, I = 0.286), but a high degree of genetic differentiation (FST = 0.521). The mean value of gene flow (Nm ) among populations was 1.082, indicating prevalent gene exchange via pollen dispersal. Phylogeographic analyses suggested that the populations of G. kiangnanense were divided into two lineages, Zhejiang (ZJ) and Anhui (AH). These two lineages were separated by the Huangshan-Tianmu Mountain Range. AMOVA analysis revealed that 36.59% of total genetic variation occurred between the two groups. The ZJ lineage was further divided into the Hangzhou (ZJH) and Zhuji (ZJZ) lineages, separated by the Longmen Mountain and Fuchun River. DIYABC analyses suggested that the ZJ and AH lineages were separated at 5.592 ka, likely due to the impact of Holocene climate change and human activities. Subsequently, the ZJZ lineage diverged from the ZJH lineage around 2.112 ka. Given the limited distribution of G. kiangnanense and the significant genetic differentiation among its lineages, both in-situ and ex-situ conservation strategies should be implemented to protect the germplasm resources of G. kiangnanense.

Y Chromosome Story-Ancient Genetic Data as a Supplementary Tool for the Analysis of Modern Croatian Genetic Pool.

Due to its turbulent demographic history, marked by extensive settlement and gene flow from diverse regions of Eurasia, Southeastern Europe (SEE) has consistently served as a genetic crossroads between East and West and a junction for the migrations that reshaped Europe's population. SEE, including modern Croatian territory, was a crucial passage from the Near East and even more distant regions and human populations in this region, as almost any other European population represents a remarkable genetic mixture. Modern humans have continuously occupied this region since the Upper Paleolithic era, and different (pre)historical events have left a distinctive genetic signature on the historical narrative of this region. Our views of its history have been mostly renewed in the last few decades by extraordinary data obtained from Y-chromosome studies. In recent times, the international research community, bringing together geneticists and archaeologists, has steadily released a growing number of ancient genomes from this region, shedding more light on its complex past population dynamics and shaping the genetic pool in Croatia and this part of Europe.

Evolution: Genomes reveal the reticulated history of giraffes.

The diversification and taxonomy of modern giraffe lineages have been a riddle for more than 200 years. A new genomic study shows that divergence with gene flow has played a significant role in the history of this zoological icon.

Broken Ring Speciation in California Mygalomorph Spiders (Nemesiidae, Calisoga).

AbstractIdealized ring species, with approximately continuous gene flow around a geographic barrier but singular reproductive isolation at a ring terminus, are rare in nature. A broken ring species model preserves the geographic setting and fundamental features of an idealized model but accommodates varying degrees of gene flow restriction over complex landscapes through evolutionary time. Here we examine broken ring species dynamics in Calisoga spiders, which, like the classic ring species Ensatina salamanders, are distributed around the Central Valley of California. Using nuclear and mitogenomic data, we test key predictions of common ancestry, ringlike biogeography, biogeographic timing, population connectivity, and terminal overlap. We show that a ring complex of populations shares a single common ancestor, and from an ancestral area in the Sierra Nevada mountains, two distributional and phylogenomic arms encircle the Central Valley. Isolation by distance occurs along these distributional arms, although gene flow restriction is also evident. Where divergent lineages meet in the South Coast Ranges, we find rare lineage sympatry, without evidence for nuclear gene flow and with clear evidence for morphological and ecological divergence. We discuss general insights provided by broken ring species and how such a model could be explored and extended in other systems and future studies.

Genetic diversity and population structure of Physella acuta (Gastropoda: Physidae) in Thailand using mitochondrial gene markers: COI and 16S rDNA.

Physella acuta is a freshwater snail native to North America. Understanding the phylogeography and genetic structure of P. acuta will help elucidate its evolution. In this study, we used mitochondrial (COI and 16S rDNA) and nuclear (ITS1) markers to identify the species and examine its genetic diversity, population structure, and demographic history of P. acuta in Thailand. Phylogenetic and network analyses of P. acuta in Thailand pertained to clade A, which exhibits a global distribution. Analysis of the genetic structure of the population revealed that the majority of pairwise comparisons showed no genetic dissimilarity. An isolation-by-distance test indicates no significant correlation between genetic and geographical distances among P. acuta populations, suggesting that gene flow is not restricted by distance. Demographic history and haplotype network analyses suggest a population expansion of P. acuta, as evidenced by the star-like structure detected in the median-joining network. Based on these results, we concluded that P. acuta in Thailand showed gene flow and recent population expansion. Our findings provide fundamental insights into the genetic variation of P. acuta in Thailand.

Population differentiation and dynamics of five pioneer species of Gaultheria from the secondary forests in subtropical China.

BACKGROUND: The influence of native secondary succession associated with anthropogenic disturbance on the biodiversity of the forests in subtropical China remains uncertain. In particular, the evolutionary response of small understory shrubs, particularly pioneer species inhabiting continuously disturbed habitats, to topographic heterogeneity and climate change is poorly understood. This study aimed to address this knowledge gap by focusing on the Gaultheria crenulata group, a clade of small pioneer shrubs in subtropical China. RESULTS: We examined the genetic structure and demographic history of all five species of the G. crenulata group with two maternally inherited chloroplast DNA (cpDNA) fragments and two biparentally inherited low-copy nuclear genes (LCG) over 89 natural populations. We found that the genetic differentiation of this group was influenced by the geomorphological boundary between different regions of China in association with Quaternary climatic events. Despite low overall genetic diversity, we observed an isolation-by-distance (IBD) pattern at a regional scale, rather than isolation-by-environment (IBE), which was attributed to ongoing human disturbance in the region. CONCLUSION: Our findings suggest that the genetic structure of the G. crenulata group reflects the interplay of geological topography, historical climates, and anthropogenic disturbance during the Pliocene-Pleistocene-Holocene periods in subtropical China. The observed IBD pattern, particularly prominent in western China, highlights the role of limited dispersal and gene flow, possibly influenced by physical barriers or decreased connectivity over geographic distance. Furthermore, the east-to-west trend of gene flow, potentially facilitated by the East Asian monsoon system, underscores the complex interplay of biotic and abiotic factors shaping the genetic dynamics of pioneer species in subtropical China's secondary forests. These findings can be used to assess the impact of environmental changes on the adaptation and persistence of biodiversity in subtropical forest ecosystems.

Genetic diversity and population structure of Bellamya purificata in Guangxi.

Bellamya purificata is an important medicinal value and economically farmed species in China. However, because little is known about the genetic characteristics of this species, the utilization of high-quality germplasm resources is hindered. The study examined the genetic differentiation between, and the structure of 12 B. purificata populations in Guangxi using 7 microsatellite DNA markers. High genetic diversity occurred in each population, with mean observed heterozygosity 0.655 and a mean expected heterozygosity 0.832. Analysis of molecular variance reveals genetic diversity to be greater within (95.2%) than among populations (4.8%). Genetic differentiation between populations is weak (Fst = 0.048, P < 0.001), with mixing of genetic clusters prevalent at the level of the individual. Genetic flow exists between populations (Nm = 3.084-11.778), with Longshui and Guilin populations exchanging frequently. A Mantel test reveals a low correlation between geographic and genetic distances (r = 0.2482, P < 0.071), suggesting that dispersal between neighboring populations facilitates population exchange. No significant heterozygosity excess was observed for any population (P > 0.05), indicating a lack of recent genetic bottlenecks. The results provide important genetic information for B. purificata, and data for potential germplasm discovery and aquaculture development.

Integrating Multiple Database Resources to Elucidate the Gene Flow in Southeast Asian Pig Populations.

The domestic pig (Sus scrofa) and its subfamilies have experienced long-term and extensive gene flow, particularly in Southeast Asia. Here, we analyzed 236 pigs, focusing on Yunnan indigenous, European commercial, East Asian, and Southeast Asian breeds, using the Pig Genomics Reference Panel (PGRP v1) of Pig Genotype-Tissue Expression (PigGTEx) to investigate gene flow and associated complex traits by integrating multiple database resources. In this study, we discovered evidence of admixtures from European pigs into the genome of Yunnan indigenous pigs. Additionally, we hypothesized that a potential conceptual gene flow route that may have contributed to the genetic composition of the Diannan small-ear pig is a gene exchange from the Vietnamese pig. Based on the most stringent gene introgression scan using the fd statistic, we identified three specific loci on chromosome 8, ranging from 51.65 to 52.45 Mb, which exhibited strong signatures of selection and harbored the NAF1, NPY1R, and NPY5R genes. These genes are associated with complex traits, such as fat mass, immunity, and litter weight, in pigs, as supported by multiple bio-functionalization databases. We utilized multiple databases to explore the potential dynamics of genetic exchange in Southeast Asian pig populations and elucidated specific gene functionalities.

City divided: Unveiling family ties and genetic structuring of coyotes in Seattle.

Linear barriers pose significant challenges for wildlife gene flow, impacting species persistence, adaptation, and evolution. While numerous studies have examined the effects of linear barriers (e.g., fences and roadways) on partitioning urban and non-urban areas, understanding their influence on gene flow within cities remains limited. Here, we investigated the impact of linear barriers on coyote (Canis latrans) population structure in Seattle, Washington, where major barriers (i.e., interstate highways and bodies of water) divide the city into distinct quadrants. Just under 1000 scats were collected to obtain genetic data between January 2021 and December 2022, allowing us to identify 73 individual coyotes. Notably, private allele analysis underscored limited interbreeding among quadrants. When comparing one quadrant to each other, there were up to 16 private alleles within a single quadrant, representing nearly 22% of the population allelic diversity. Our analysis revealed weak isolation by distance, and despite being a highly mobile species, genetic structuring was apparent between quadrants even with extremely short geographic distance between individual coyotes, implying that Interstate 5 and the Ship Canal act as major barriers. This study uses coyotes as a model species for understanding urban gene flow and its consequences in cities, a crucial component for bolstering conservation of rarer species and developing wildlife friendly cities.

Genomic insights into hybrid zone formation: The role of climate, landscape, and demography in the emergence of a novel hybrid lineage.

Population demographic changes, alongside landscape, geographic and climate heterogeneity, can influence the timing, stability and extent of introgression where species hybridise. Thus, quantifying interactions across diverged lineages, and the relative contributions of interspecific genetic exchange and selection to divergence at the genome-wide level is needed to better understand the drivers of hybrid zone formation and maintenance. We used seven latitudinally arrayed transects to quantify the contributions of climate, geography and landscape features to broad patterns of genetic structure across the hybrid zone of Populus trichocarpa and P. balsamifera and evaluated the demographic context of hybridisation over time. We found genetic structure differed among the seven transects. While ancestry was structured by climate, landscape features influenced gene flow dynamics. Demographic models indicated a secondary contact event may have influenced contemporary hybrid zone formation with the origin of a putative hybrid lineage that inhabits regions with higher aridity than either of the ancestral groups. Phylogenetic relationships based on chloroplast genomes support the origin of this hybrid lineage inferred from demographic models based on the nuclear data. Our results point towards the importance of climate and landscape patterns in structuring the contact zones between P. trichocarpa and P. balsamifera and emphasise the value whole genome sequencing can have to advancing our understanding of how neutral processes influence divergence across space and time.