Genetic Discontinuity Research Articles

Genomic data reveal a north-south split and introgression history of blood fluke (Schistosoma haematobium) populations from across Africa.

The human parasitic fluke, Schistosoma haematobium hybridizes with the livestock parasite S. bovis in the laboratory, but the extent of hybridization in nature is unclear. We analyzed 34.6 million single nucleotide variants in 162 samples from 18 African countries, revealing a sharp genetic discontinuity between northern and southern S. haematobium. We found no evidence for recent hybridization. Instead the data reveal admixture events that occurred 257-879 generations ago in northern S. haematobium populations. Fifteen introgressed S. bovis genes are approaching fixation in northern S. haematobium with four genes potentially driving adaptation. We identified 19 regions that were resistant to introgression; these were enriched on the sex chromosomes. These results (i) demonstrate strong barriers to gene flow between these species, (ii) indicate that hybridization may be less common than currently envisaged, but (iii) reveal profound genomic consequences of interspecific hybridization between schistosomes of medical and veterinary importance.

Multi-proxy bioarchaeological analysis ofskeletal remains shows genetic discontinuity in a Medieval Sicilian community.

The medieval period in Sicily was turbulent, involving successiveregime changes, from Byzantine (Greek Christian), Aghlabid (Sunni Muslim), Fatimid (Shī'a Muslim), to Normans and Swabians (Latin Christian). To shed new light on the localimplications of regime changes, we conducted amultidisciplinary analysis of 27 individuals buried in adjacentMuslim and Christian cemeteries at the site of Segesta,western Sicily. By combining radiocarbon dating, genome-wide sequencing, stable and radiogenic isotopic data, and archaeological records, we uncover genetic differences between the two communities but find evidence of continuity in other aspects of life. Historical and archaeological evidence shows a Muslim community was present by the 12th century during Norman governance, with the Christian settlement appearing in the 13th century under Swabian governance. ABayesian analysis of radiocarbon dates from the burials finds the abandonment of the Muslim cemetery likely occurred after the establishment of the Christian cemetery, indicating that individuals of both faiths were present in the area in the first half of the 13th century. The biomolecular results suggest the Christians remained genetically distinct from the Muslim community at Segesta while following a substantially similar diet. This study demonstrates that medieval regime changes had major impacts beyond the political core, leading to demographic changes while economic systems persisted and new social relationships emerged.

Network of large pedigrees reveals social practices of Avar communities.

From AD 567-568, at the onset of the Avar period, populations from the Eurasian Steppe settled in the Carpathian Basin for approximately 250 years1. Extensive sampling for archaeogenomics (424 individuals) and isotopes, combined with archaeological, anthropological and historical contextualization of four Avar-period cemeteries, allowed for a detailed description of the genomic structure of these communities and their kinship and social practices. We present a set of large pedigrees, reconstructed using ancient DNA, spanning nine generations and comprising around 300 individuals. We uncover a strict patrilineal kinship system, in which patrilocality and female exogamy were the norm and multiple reproductive partnering and levirate unions were common. The absence of consanguinity indicates that this society maintained a detailed memory of ancestry over generations. These kinship practices correspond with previous evidence from historical sources and anthropological research on Eurasian Steppe societies2. Network analyses of identity-by-descent DNA connections suggest that social cohesion between communities was maintained via female exogamy. Finally, despite the absence of major ancestry shifts, the level of resolution of our analyses allowed us to detect genetic discontinuity caused by the replacement of a community at one of the sites. This was paralleled with changes in the archaeological record and was probably a result of local political realignment.

Population structure and microscale morphological differentiation in a freshwater snail from the Chilean Altiplano

BackgroundThe diversity and population genetic structure of many species have been shaped by historical and contemporary climatic changes. For the species of the South American Altiplano, the historical climatic changes are mainly related to the wet events of great magnitude and regional influence that occurred during the Pleistocene climatic oscillations (PCOs). In contrast, contemporary climate changes are associated with events of lesser magnitude and local influence related to intensifications of the South American Summer Monsoon (SASM). Although multiple studies have analyzed the effect of PCOs on the genetic patterns of highland aquatic species, little is known about the impact of contemporary climate changes in recent evolutionary history. Therefore, in this study, we investigated the change in population structure and connectivity using nuclear and mitochondrial markers throughout the distribution range of Heleobia ascotanensis, a freshwater Cochliopidae endemic to the Ascotán Saltpan. In addition, using geometric morphometric analyses, we evaluated the concomitance of genetic divergence and morphological differentiation.ResultsThe mitochondrial sequence analysis results revealed the presence of highly divergent co-distributed and geographically nested haplotypes. This pattern reflects an extension in the distribution of groups that previously would have differentiated allopatrically. These changes in distribution would have covered the entire saltpan and would be associated with the large-scale wet events of the PCOs. On the other hand, the microsatellite results defined five spatially isolated populations, separated primarily by geographic barriers. Contemporary gene flow analyses suggest that post-PCO, climatic events that would have connected all populations did not occur. The morphometric analyses results indicate that there is significant morphological differentiation in the populations that are more isolated and that present the greatest genetic divergence.ConclusionsThe contemporary population structure and morphological variation of H. ascotanensis mainly reflect the post-PCO climatic influence. Although both markers exhibit high genetic structuring, the microsatellite and morphology results show the preponderant influence of fragmentation in recent evolutionary history. The contemporary genetic pattern shows that in species that have limited dispersal capabilities, genetic discontinuities can appear rapidly, erasing signs of historical connectivity.

Microgeographic speciation in a complex of Anatolian bush crickets facilitated by fast evolution of reproductive isolation.

Identifying the drivers of microgeographic speciation (i.e., speciation over small, local geographic scales) is key to understand the origin of speciose groups. Here, we use genomic data to infer the demographic processes underlying diversification in Poecilimon luschani (Orthoptera: Tettigoniidae), a species complex belonging to the most diverse genus of bush crickets from the Mediterranean region (>170 taxa) that comprises three recognized subspecies with small allopatric distributions in the topographically complex Teke Peninsula, southwestern Anatolia. Phylogenomic reconstructions that include all other taxa within the species group confirmed that subspecies of P. luschani originated from a common ancestor during the Pleistocene, supporting recent (<1 Ma) diversification within a small geographical area (ca. 120 × 80 km). Genetic clustering analyses corroborated the distinctiveness of each subspecies and the cohesiveness of their respective populations, with abrupt genetic discontinuities coinciding with contemporary range boundaries. Indeed, our analyses uncovered the presence of two sympatric cryptic sister lineages that diverged <300 ka ago and do not admix despite being co-distributed. Collectively, these results support that all lineages within the complex represent independently evolving entities corresponding to full-fledged species. Statistical evaluation of alternative models of speciation strongly supports a scenario of divergence in isolation followed by a period of limited gene flow during the last glacial period, when all lineages experienced marked expansions according to demographic reconstructions. Our study exemplifies how localized allopatric divergence and fast evolution of reproductive isolation can promote microgeographic speciation and explain the high rates of endemism characterizing biodiversity hotspots.

The fragrant orchid &lt;i&gt;Vanilla chamissonis&lt;/i&gt;: an appraisal of the genetic structure of wild populations from Osununú Natural Reserve (Misiones, Argentina)

Due to diverse factors the commercial crops of vanilla are highly vulnerable to genetic erosion. The species Vanilla chamissonis Klotzsch has been shown to produce interesting bioeconomic compounds that could diversify the market. In Argentina, the distribution of this fragrant orchid is restricted to specific north-eastern regions, like the Osununú Natural Reserve (at San Ignacio, Misiones), where it grows wildly in sectors (nuclei). This germplasm represents a locally adapted endemism, and thus is a valuable resource worth preserving. Herein, the genetic characterization of unexplored V. chamissonis´ germplasm was accomplished; as well, the specimens´ molecular identification was assessed in a phylogenetic context using nuclear and chloroplastidial markers. The genotyping exposed that, at least, four gene pools contribute to the genetic diversity of these plants, and that ca. 81% of total genetic variation is allocated within populations. Present results suggest a predominance of sexual reproduction, that population N1 emerges as a candidate unit worth conserving, and that a genetic discontinuity was evidenced between northeastern and southwestern nuclei due to their different genetic constitutions. In concurrence, a tenuous isolation by distance was also unveiled. Our phylogenetic results suggest that it is appropriate to consider the vanilla germplasm surveyed as V. chamissonis. The basic knowledge generated in this work will aid in guiding the in-situ conservation and management initiative, ongoing at the Osununú´s reserve.

Contrasting patterns of genetic diversity and population discontinuity in the common red seaweed Gracilaria salicornia (Gracilariaceae) along the coasts of Thailand

ABSTRACT For long-term management and conservation of marine resources an understanding of their genetic diversity and geographic distribution is necessary. The red alga Gracilaria salicornia is widely distributed along the coasts of Thailand and has been used as food and as fertilizer. However, there have been no studies examining the genetic diversity of this alga. The present research aimed to employ two different genetic markers, mitochondrial (COI-5P) and chloroplast (RuBisCO spacer) for exploring the genetic diversity and phylogeography of G. salicornia in Thailand. Our analyses, inferred from both markers, revealed different levels of genetic diversity between the two coasts of Thailand, and the recognition of seven new COI haplotypes. Gracilaria salicornia in the Andaman Sea had higher genetic diversity than in the Gulf of Thailand. Significant genetic divergence was found between populations in the Andaman Sea and the Gulf of Thailand, likely due to the geographic barrier of the Thai-Malay Peninsula and the complex ocean currents. No genetic structuring was detected in the Gulf of Thailand, whereas low to moderate genetic differentiation was found between populations in the Andaman Sea, suggesting a possible break between the northern and southern Andaman Sea. This break could be explained by ocean circulation patterns around the Malacca Strait. These findings provide further evidence supporting the Thai-Malay Peninsula and oceanographic currents as strong barriers to dispersal, leading to genetic differentiation and population structuring of G. salicornia and other marine species. These results can be used to guide conservation management actions and exploitation of this marine resource.

Regional (south‐eastern Pacific Ocean) population genetics and global phylogeography of two endangered highly migratory pelagic sharks, the blue shark Prionace glauca and shortfin mako Isurus oxyrinchus

Abstract The blue shark Prionace glauca and the shortfin mako Isurus oxyrinchus are two large and highly migratory sharks inhabiting temperate and tropical waters worldwide that are heavily targeted by artisanal and industrial fisheries. The International Union for Nature Conservation classifies the blue shark and shortfin mako as ‘Near Threatened’ and ‘Vulnerable’, respectively, in the Nature Red List of Threatened Species v. 2019‐2. This study examined the population genetics of the shortfin mako and blue sharks at a regional (south‐eastern Pacific Ocean) and global scale. The null hypothesis of no genetic discontinuities among ocean basins and/or between hemispheres was tested using two mitochondrial markers suitable for population genetics inference in these species: the non‐coding control region and the protein‐coding gene cytochrome c oxidase I in I. oxyrinchus, and the control region and cytochrome b in P. glauca. Spatial genetic analyses suggested a single and two genetic clusters co‐occurring along the south‐eastern Pacific Ocean in the shortfin mako and blue shark, respectively. Phylogeographic analyses, migration estimates, haplotype networks and AMOVAs demonstrated that the two species exhibit an overall pattern of high genetic connectivity among hemispheres and across ocean basins with a signature of shallow genetic structuring worldwide. This study has generated valuable information for the management and conservation of heavily exploited sharks and highlights the need for additional inclusive research programmes assessing inter‐regional genomic discontinuities using more statistically powerful genetic markers to determine with precision population genetic discontinuities (if any) in these and other highly migratory sharks.

The role of the Sunda shelf biogeographic barrier in the cryptic differentiation of Conus litteratus (Gastropoda: Conidae) across the Indo-Pacific region.

Geographical and oceanographic processes have influenced the speciation of marine organisms. Cone snails are marine mollusks that show high levels of endemism and a wide distributional range across the Indian and Pacific Oceans. Discontinuities in distributions caused by biogeographic barriers can affect genetic connectivity. Here we analysed the connectivity within Conus litteratus using samples from the Lakshadweep archipelago (Arabian Sea, Indian Ocean) and from the Pacific Ocean. Maximum likelihood analyses based on the mitochondrial cytochrome C oxidase subunit I (COI) and on the non-coding 16S ribosomal RNA (16S rRNA) genes revealed cryptic diversity within C. literatus occupying distinct oceanographic regions. The intraspecific genetic distances between the two distinct clades of C. literatus from the Arabian Sea and the Pacific Ocean ranged from 7.4% to 7.6% for COI and from 2.4% to 2.8% for 16S rRNA genes, which is larger than the threshold limit for interspecific differentiation. The haplotype network analysis also corroborated the existence of two different lineages within C. litteratus. The detected genetic discontinuities reflect the effect of the Sunda shelf biogeographic barrier on the allopatric divergence of C. litteratus.

Phylogeography of the Asian green mussel, Perna viridis: pronounced regional genetic differentiation is explained by major freshwater input into the Bay of Bengal

The Asian green mussel (AGM), Perna viridis, is a mussel of ecological and economic importance throughout much of the tropical and subtropical regions of Asia, from China and Korea in the east to Oman on the Arabian Peninsula in the west. We collated published DNA sequence data and then analysed 467 bp of cytochrome c oxidase I gene (COI) sequence variation from 649 mussels (576 from SE Asia and 73 from India) to test for regional differentiation. Analysis of molecular variation, a haplotype network and a neighbour-joining tree all revealed significant differentiation between mussels from India and those from SE Asia. We observed a greater number of segregating sites (haplotypes) and private haplotypes than expected for Indian mussels compared to SE Asian mussels, based on the respective sample sizes. Tests of neutrality revealed population expansion or a recent selective sweep for only the mussels from India. We suggest that the differentiation of the Indian from the SE Asian mussels is explained by the pronounced seasonal input of freshwater into the northern Bay of Bengal area, which acts as a (semi-permeable) barrier to gene flow between the two regions (i.e., east versus west of the Bay). This suggestion is based on the oceanography of the region and is consistent with the biogeography and reports of genetic breaks in a range of taxa in this region. Further targetted sampling of AGMs from the east coast of India, Bangladesh and Myanmar is required to test this hypothesis, and additional sampling of AGMs from west of India (e.g., Pakistan, Iran, the Arabian Peninsula) will also be informative. Finally, this study is based on the analysis of a single marker (COI) because this is all that is presently available from published sources but we note that the application of new molecular markers such as single nucleotide polymorphisms to newly collected mussel samples will greatly advance our understanding of the AGM genetic discontinuity and its age, and help test its geographical location and its nature (e.g., a steep/shallow cline or a stepped cline).

Biogeographic provinces and genomically delineated stocks are congruent in snapper (Chrysophrys auratus) from southeastern Australia

AbstractMarine species often exhibit genetic discontinuities concordant with biogeographic boundaries, frequently occurring due to changes in ocean circulation, bathymetry, coastline topography, and temperature. Here, we used 10 916 single nucleotide polymorphisms (SNPs) to assess the concordance between population genomic differentiation and coastal biogeography in the fishery-important snapper (Chrysophrys auratus) across southeastern Australia. Additionally, we investigated whether spatial scales of assessment and management of snapper align with evidence from population genomics. We detected genomic structure concordant with the region's three biogeographic provinces across snapper from 11 localities (n = 488) between the west coast of South Australia and the south coast of New South Wales. We also detected fine-scale genetic structuring relating to spatial variation in spawning and recruitment dynamics, as well as temporal stability in the genomic signal associated with two important spawning areas. The current management boundaries in the region coincided with either the genetic breaks at bioregional boundaries or with local-scale variation. Our study highlights the value of population genomic surveys in species with high dispersal potential for uncovering stock boundaries and demographic variation related to spawning and recruitment. It also illustrates the importance of marine biogeography in shaping population structure in commercial species with high dispersal potential.

An HLA map of the world: A comparison of HLA frequencies in 200 worldwide populations reveals diverse patterns for class I and class II.

HLA frequencies show widespread variation across human populations. Demographic factors as well as selection are thought to have shaped HLA variation across continents. In this study, a worldwide comparison of HLA class I and class II diversity was carried out. Multidimensional scaling techniques were applied to 50 HLA-A and HLA-B (class I) as well as 13 HLA-DRB1 (class II) first-field frequencies in 200 populations from all continents. Our results confirm a strong effect of geography on the distribution of HLA class I allele groups, with principal coordinates analysis closely resembling geographical location of populations, especially those of Africa-Eurasia. Conversely, class II frequencies stratify populations along a continuum of differentiation less clearly correlated to actual geographic location. Double clustering analysis revealed finer intra-continental sub-clusters (e.g., Northern and Western Europe vs. South East Europe, North Africa and Southwest Asia; South and East Africa vs. West Africa), and HLA allele group patterns characteristic of these clusters. Ancient (Austronesian expansion) and more recent (Romani people in Europe) migrations, as well as extreme differentiation (Taiwan indigenous peoples, Native Americans), and interregional gene flow (Sámi, Egyptians) are also reflected by the results. Barrier analysis comparing DST and geographic location identified genetic discontinuities caused by natural barriers or human behavior explaining inter and intra-continental HLA borders for class I and class II. Overall, a progressive reduction in HLA diversity from African to Oceanian and Native American populations is noted. This analysis of HLA frequencies in a unique set of worldwide populations confirms previous findings on the remarkable similarity of class I frequencies to geography, but also shows a more complex development for class II, with implications for both human evolutionary studies and biomedical research.

It's a small world for parasites: evidence supporting the North American invasion of European Echinococcus multilocularis.

Echinococcus multilocularis (Em), the causative agent of human alveolar echinococcosis (AE), is present in the Holarctic region, and several genetic variants deem to have differential infectivity and pathogenicity. An unprecedented outbreak of human AE cases in Western Canada infected with a European-like strain circulating in wild hosts warranted assessment of whether this strain was derived from a recent invasion or was endemic but undetected. Using nuclear and mitochondrial markers, we investigated the genetic diversity of Em in wild coyotes and red foxes from Western Canada, compared the genetic variants identified to global isolates and assessed their spatial distribution to infer possible invasion dynamics. Genetic variants from Western Canada were closely related to the original European clade, with lesser genetic diversity than that expected for a long-established strain and spatial genetic discontinuities within the study area, supporting the hypothesis of a relatively recent invasion with various founder events.

Habitat connectivity and host relatedness influence virus spread across an urbanising landscape in a fragmentation-sensitive carnivore.

Spatially heterogeneous landscape factors such as urbanisation can have substantial effects on the severity and spread of wildlife diseases. However, research linking patterns of pathogen transmission to landscape features remains rare. Using a combination of phylogeographic and machine learning approaches, we tested the influence of landscape and host factors on feline immunodeficiency virus (FIVLru) genetic variation and spread among bobcats (Lynx rufus) sampled from coastal southern California. We found evidence for increased rates of FIVLru lineage spread through areas of higher vegetation density. Furthermore, single-nucleotide polymorphism (SNP) variation among FIVLru sequences was associated with host genetic distances and geographic location, with FIVLru genetic discontinuities precisely correlating with known urban barriers to host dispersal. An effect of forest land cover on FIVLru SNP variation was likely attributable to host population structure and differences in forest land cover between different populations. Taken together, these results suggest that the spread of FIVLru is constrained by large-scale urban barriers to host movement. Although urbanisation at fine spatial scales did not appear to directly influence virus transmission or spread, we found evidence that viruses transmit and spread more quickly through areas containing higher proportions of natural habitat. These multiple lines of evidence demonstrate how urbanisation can change patterns of contact-dependent pathogen transmission and provide insights into how continued urban development may influence the incidence and management of wildlife disease.

Drivers of population divergence and genetic variation in Elymus breviaristatus (Keng) Keng f. (Poaceae: Triticeae), an endemic perennial herb of the Qinghai-Tibet plateau

Elymus breviaristatus, a rare grass species with excellent resistance and ecological importance, is narrowly distributed on the Qinghai-Tibet plateau. Populations of E. breviaristatus are declining due to habitat fragmentation, and thus far, characteristics of genetic differentiation and adaptive responses to climate change remain poorly understood in this species. Here, we explored the genetic structure of 18 natural populations (269 individuals) in the transition zone between Tibet and the Hengduan Mountains using 15 expressed sequence tag (EST)-SSR primer pairs and identified possible barriers to gene flow that might have caused genetic discontinuities. Additional analyses were performed to identify the environmental factors affecting genetic diversity and to test whether the patterns of genetic variation among populations were more consistent with the isolation by distance (IBD) or isolation by environment (IBE) model. Multiple measures of genetic diversity revealed that intra-population genetic variation was low, while inter-population genetic variation was high. Clustering, structure, and principal coordinate analyses identified three genetic groups: (a) Eastern Qamdo, (b) Nagqu and Western Qamdo, and (c) Lhasa and Nyingchi. A clear physical barrier to gene flow was formed by the Yarlung Zangbo Grand Canyon and the Tanggula Mountains. We found that both IBD and IBE contributed to the observed patterns of genetic variation, and the IBE model played a leading role. In addition, precipitation-related variables, soil phosphorus content and soil K:P ratio significantly affected population genetic variation. Overall, our results emphasized the genetic fragility of E. breviaristatus populations and showed that this species requires attention, as future climate changes and human activities may further threaten its survival. In addition, the genetic differences among E. breviaristatus populations should be considered when formulating conservation measures for E. breviaristatus populations in the study area.

Fisheries genomics of snapper (Chrysophrys auratus) along the west Australian coast.

The efficacy of fisheries management strategies depends on stock assessment and management actions being carried out at appropriate spatial scales. This requires understanding of spatial and temporal population structure and connectivity, which is challenging in weakly structured and highly connected marine populations. We carried out a population genomics study of the heavily exploited snapper (Chrysophrys auratus) along ~2600 km of the Australian coastline, with a focus on Western Australia (WA). We used 10,903 filtered SNPs in 341 individuals from eight sampling locations to characterize population structure and connectivity in snapper across WA and to assess if current spatial scales of stock assessment and management agree with evidence from population genomics. Our dataset also enabled us to investigate temporal stability in population structure as well as connectivity between WA and its nearest, eastern jurisdictional neighbour. As expected for a species influenced by the extensive ocean boundary current in the region, low genetic differentiation and high connectivity were uncovered across WA. However, we did detect strong isolation by distance and genetic discontinuities in the mid‐west and south‐east. The discontinuities correlate with boundaries between biogeographic regions, influenced by on‐shelf oceanography, and the sites of important spawning aggregations. We also detected temporal instability in genetic structure at one of our sites, possibly due to interannual variability in recruitment in adjacent regions. Our results partly contrast with the current spatial management of snapper in WA, indicating the likely benefits of a review. This study supports the value of population genomic surveys in informing the management of weakly structured and wide‐ranging marine fishery resources.

Stepped Geomorphology Shaped the Phylogeographic Structure of a Widespread Tree Species (Toxicodendron vernicifluum, Anacardiaceae) in East Asia

Species’ phylogeographic patterns reflect the interplay between landscape features, climatic forces, and evolutionary processes. Here, we used two chloroplast DNA (cpDNA) markers (trnL and trnL-F) to explore the role of stepped geomorphology in shaping the phylogeographic structure of Toxicodendron vernicifluum, an economically important tree species widely distributed in East Asia. The range-wide pattern of sequence variation was analyzed based on a dataset including 357 individuals from China, together with published sequences of 92 individuals mainly from Japan and South Korea. We identified five chloroplast haplotypes based on seven substitutions across the 717-bp alignment. A clear east-west phylogeographic break was recovered according to the stepped landforms of mainland China. The wild trees of the western clade were found to be geographically restricted to the “middle step”, which is characterized by high mountains and plateaus, while those of the eastern clade were confined to the “low step”, which is mainly made up of hills and plains. The two major clades were estimated to have diverged during the Early Pleistocene, suggesting that the cool glacial climate may have caused the ancestral population to retreat to at least two glacial refugia, leading to allopatric divergence in response to long-term geographic isolation. Migration vector analyses based on the outputs of ecological niche models (ENMs) supported a gradual range expansion since the Last Interglacial. Mountain ranges in western China and the East China Sea land bridge were inferred to be dispersal corridors in the western and eastern distributions of T. vernicifluum, respectively. Overall, our study provides solid evidence for the role of stepped geomorphology in shaping the phylogeographic patterns of T. vernicifluum. The resulting east-west genetic discontinuities could persist for a long time, and could occur at a much larger scale than previously reported, extending from subtropical (e.g., the Xuefeng Mountain) to warm-temperate China (e.g., the Taihang Mountain).

Genomic and dietary discontinuities during the Mesolithic and Neolithic in Sicily

SummarySicily is a key region for understanding the agricultural transition in the Mediterranean because of its central position. Here, we present genomic and stable isotopic data for 19 prehistoric Sicilians covering the Mesolithic to Bronze Age periods (10,700–4,100 yBP). We find that Early Mesolithic hunter-gatherers (HGs) from Sicily are a highly drifted lineage of the Early Holocene western European HGs, whereas Late Mesolithic HGs carry ∼20% ancestry related to northern and (south) eastern European HGs, indicating substantial gene flow. Early Neolithic farmers are genetically most similar to farmers from the Balkans and Greece, with only ∼7% of ancestry from local Mesolithic HGs. The genetic discontinuities during the Mesolithic and Early Neolithic match the changes in material culture and diet. Three outlying individuals dated to ∼8,000 yBP; however, suggest that hunter-gatherers interacted with incoming farmers at Grotta dell’Uzzo, resulting in a mixed economy and diet for a brief interlude at the Mesolithic-Neolithic transition.

Main roads and land cover shaped the genetic structure of a Mediterranean island wild boar population.

Patterns of genetic differentiation within and among animal populations might vary due to the simple effect of distance or landscape features hindering gene flow. An assessment of how landscape connectivity affects gene flow can help guide management, especially in fragmented landscapes. Our objective was to analyze population genetic structure and landscape genetics of the native wild boar (Sus scrofa meridionalis) population inhabiting the island of Sardinia (Italy), and test for the existence of Isolation‐by‐Distance (IBD), Isolation‐by‐Barrier (IBB), and Isolation‐by‐Resistance (IBR). A total of 393 Sardinian wild boar samples were analyzed using a set of 16 microsatellite loci. Signals of genetic introgression from introduced non‐native wild boars or from domestic pigs were revealed by a Bayesian cluster analysis including 250 reference individuals belonging to European wild populations and domestic breeds. After removal of introgressed individuals, genetic structure in the population was investigated by different statistical approaches, supporting a partition into five discrete subpopulations, corresponding to five geographic areas on the island: north‐west (NW), central west (CW), south‐west (SW), north‐central east (NCE), and south‐east (SE). To test the IBD, IBB, and IBR hypotheses, we optimized resistance surfaces using genetic algorithms and linear mixed‐effects models with a maximum likelihood population effects parameterization. Landscape genetics analyses revealed that genetic discontinuities between subpopulations can be explained by landscape elements, suggesting that main roads, urban settings, and intensively cultivated areas are hampering gene flow (and thus individual movements) within the Sardinian wild boar population. Our results reveal how human‐transformed landscapes can affect genetic connectivity even in a large‐sized and highly mobile mammal such as the wild boar, and provide crucial information to manage the spread of pathogens, including the African Swine Fever virus, endemic in Sardinia.

Bovine landscape genetics in Brazil

The aim of this study was to explore spatial genetic structure patterns in cattle breeds adapted to local conditions in Brazil. We georeferenced 876 animals of ten breeds raised in Brazil kept in the Genebank of Brazilian Agricultural Research Corporation (Embrapa) by sample locations using the QGIS 2.4.0 software. The Mantel tests, spatial autocorrelation, and Monmonier tests were performed. The distances for spatial correlation tests ranged from 5 to 15 classes. The results indicated genetic discontinuities in cattle breeds from the Midwest, South, and Southeast of the country. Correlation between genetic and geographic distance was low, but significant. The Monmonier Maximum Distance Algorithm indicated an initial subdivision of Curraleiro and then Pantaneiro from the other breeds. In another subdivision, Criollo, Mocho Nacional, and Caracu were grouped. Genetic discontinuity was observed beyond 431 km, the minimum sampling distance between populations for conservation purposes.