Not all isolates are equal: linkage disequilibrium analysis on Xq13.3 reveals different patterns in Sardinian sub-populations

Recent studies indicate that, whereas the Sardinian population as a whole is comparable to outbred populations for linkage disequilibrium (LD) mapping of common variants, LD in Sardinian sub-isolates is more extended, making these populations particularly suitable for this approach. To evaluate the extent of LD between microsatellite markers, we compared different sub-populations within Sardinia selected on the basis of their geographical position and isolation: two small isolated villages (Talana, Urzulei), two larger but remote areas (Ogliastra, Nuoro province) and a cohort of samples representing the wider Sardinian population. LD analysis was carried out by using six microsatellite markers that are located on Xq13.3 and that have been extensively studied in different populations. We found different extents and patterns of LD in the sub-population samples depending on their degree of isolation and demographic history. All LD measurements and haplotype analyses indicate that there is a decreasing trend from Talana (the most inbred population, LD up to 9.5-11.5 Mb) to the more outbred Sardinian population (LD only for intervals <2 Mb). In one village (Talana), five haplotype classes accounting for 80% of the entire sample perfectly matched five Ogliastra clusters, supporting the origin of the village from the Ogliastra genetic pool. In contrast, the other village (Urzulei) showed a different pattern of haplotypes with a closer relationship to the Nuoro region sub-population. LD analyses therefore show that even neighbouring isolate villages may differ in their genetic background. Here, we highlight the importance of selecting appropriate populations and/or sub-populations for the analysis of complex traits. Isolated sub-populations showing different extents of LD can provide a powerful method for mapping complex traits by LD scanning at relatively low marker density.     

Striking differentiation of sub-populations within a genetically homogeneous isolate (Ogliastra) in Sardinia as revealed by mtDNA analysis

Since the reduced genetic diversity found in isolates should simplify the study of complex traits, analyses of patterns of homogeneity within populations are of particular interest. We analysed the mtDNA haplogroups and hypervariable segment I (HVS-I) sequences of 475 individuals from a geographically restricted and isolated area (Ogliastra) within Sardinia, comprehending 175 random samples from 20 out of 23 villages. The remaining 300 subjects were chosen from the other three villages, Talana, Urzulei and Perdasdefogu, by sampling all maternal lineages. A comparison with other European populations reveals that Ogliastra ranks among the most genetically homogenous population and that it has been small and isolated throughout its history. The lack of variation and the high genetic homogeneity indicate that an important founder event and a demographic expansion took place during the Neolithic (~ 7,700 years before present) in Ogliastra's mtDNA gene pool. We present highly resolved phylogenetic networks for Ogliastra and for the three sub-isolates. MtDNA differentiation in the sub-populations versus Ogliastra is revealed by a strong demarcation in their genetic pools due to distinctive founder effects and genetic drift. We found that genetic homogeneity strictly depends on a scale factor in population size and on sampling methodology. The outstanding homogeneity and the reduced female gene pool observed in Ogliastra, in the European context, hide an extremely marked differentiation in sub-isolates originated from the same archaic population. Although Ogliastra can be considered a genetically homogeneous isolate, small villages' divergent genetic histories underline the importance of more systematic analysis of DNA variation between and within populations.C. Fraumene and E. Petretto contributed equally to this work     

Archival,demographic and genetic studies define a Sardinian sub-isolate as a suitable model for mapping complex traits

Genetic isolates represent exceptional resources for the mapping of complex traits but not all isolates are similar. We have selected a genetic and cultural isolate, the village of Talana from an isolated area of Sardinia, and propose that this population is suitable for the mapping of complex traits. A wealth of historical and archive data allowed the reconstruction of the demographic and genealogical history of the village. Key features of the population, which has grown slowly with no significant immigration, were defined by using a combination of historical, demographic and genetic studies. The genealogy of each Talana inhabitant was reconstructed and the main maternal and paternal lineages of the village were defined. Haplotype and phylogenetic analyses of the Y chromosome and characterisation of mitochondrial DNA haplogroups were used to determine the number of ancestral village founders. The extent of linkage disequilibrium (LD) was evaluated by the analysis of several microsatellites in chromosomal region Xq13.3, which was previously used to asses the extension of LD. Genealogical reconstructions were confirmed and reinforced by the genetic analyses, since some lineages were found to have merged prior to the beginning of the archival records, suggesting an even smaller number of founders than initially predicted. About 80% of the present-day population appears to derive from eight paternal and eleven maternal ancestral lineages. LD was found to span, on average, a 5-Mb region in Xq13.3. This suggests the possibility of identifying identical-by-descent regions associated with complex traits in a genome-wide search by using a low-density marker map. The present study emphasises the importance of combining genetic studies with genealogical and historical information.     

Endogamy and variation in blood pressure levels in Croatian island isolates

Blood pressure variation was investigated among populations inhabiting islands and peninsula of Middle Dalmatia, Croatia. The number of previous anthropological studies pointed to isolation and different genetic population structure in this environmentally fairly homogeneous area. Variation in blood pressure (systolic and diastolic) among the populations of the islands of Brac, Hvar, Korcula, and the Peljesac peninsula was assessed at three levels involving village populations, regional (western and eastern) populations and the entire island populations. The blood pressure data were collected from 3834 adult individuals inhabiting 37 rural communities and were adjusted for age and body mass index. Variation in blood pressure levels existed among regions and villages. Due to the history of differential settlement, small village sizes and high levels of reproductive isolation, the observed blood pressure variation could be attributed to founder effect, genetic drift and inbreeding. The involvement of genetic factors was tested by relating blood pressure variation among villages to degree of isolation among them. Blood pressure means and proportions of hypertensives increased with endogamy levels in males. In females, this effect could not be observed. However, in both sexes the highest proportions of hypertensives (more than 40%) were found in villages that are most reproductively closed (endogamy greater than 80%). These populations are considered particularly promising for further genetic epidemiological research.     

Molecular and quantitative genetic differentiation across Europe in yellow dung flies

Relating geographic variation in quantitative traits to underlying population structure is crucial for understanding processes driving population differentiation, isolation and ultimately speciation. Our study represents a comprehensive population genetic survey of the yellow dung fly Scathophaga stercoraria, an important model organism for evolutionary and ecological studies, over a broad geographic scale across Europe (10 populations from the Swiss Alps to Iceland). We simultaneously assessed differentiation in five quantitative traits (body size, development time, growth rate, proportion of diapausing individuals and duration of diapause), to compare differentiation in neutral marker loci (F(ST)) to that of quantitative traits (Q(ST)). Despite long distances and uninhabitable areas between sampled populations, population structuring was very low but significant (F(ST) = 0.007, 13 microsatellite markers; F(ST) = 0.012, three allozyme markers; F(ST) = 0.007, markers combined). However, only two populations (Iceland and Sweden) showed significant allelic differentiation to all other populations. We estimated high levels of gene flow [effective number of migrants (Nm) = 6.2], there was no isolation by distance, and no indication of past genetic bottlenecks (i.e. founder events) and associated loss of genetic diversity in any northern or island population. In contrast to the low population structure, quantitative traits were strongly genetically differentiated among populations, following latitudinal clines, suggesting that selection is responsible for life history differentiation in yellow dung flies across Europe.     

Latent structure of dermatoglyphs in the population of Selska Valley.

The historical records of Selska Valley reveal that the eastern part of this area was first settled by Slovene agrarian colonists, the western part by German colonists and the central part by Friulians. These were later followed by Slovene and Slovenized settlers, who penetrated the valley from north to south. Because of its reproductive isolation, the population of Selska Valley is highly suitable for the study of population structures. The quantitative traits of the digital and palmar dermatoglyphs are polygenetically determined characteristics, which, due to their selective inertness to changes, may provide an insight into microevolutionary processes. The purpose of our study was to identify the possible differences between the populations of villages in the valley and the mountain villages attributable to various migration flows through history. Altogether 340 finger and palm prints of 163 males and 177 females were collected in two groups of villages: (1) the lowland villages (Praprotno, Bukovica, Sevlje, Dolenja vas, Selca, Zelezniki and Zali log), and (2) the mountain villages (Podlonk, Prtovc, Spodnje Danje, Zgornja Sorica and Spodnja Sorica). The 18 dermatoglyphic variables were analyzed. A statistical analysis using standard methods was performed and the latent structure evaluated using factor analysis. The discriminant analysis and latent structure of the quantitative properties of dermatoglyphs suggest the presence of certain differences in gene pools of two studied populations (the group of villages in the valley and the group of mountain villages). It is highly probable that these differences can be attributed to low migration in the Selska Valley and to the 'selective inertness' of quantitative dermatoglyphic traits. In a previous study, no significant biological differences between the studied populations were found in qualitative dermatoglyphic traits. This indicates that Selska Valley and its village populations represent a specific isolate, and therefore expanded studies of this type could significantly contribute to a better understanding of the population concerned and isolates in general.     

Population structure inferred by local spatial autocorrelation: an example from an Amerindian tribal population

Spatial autocorrelation (SA) methods were recently extended to detect local spatial autocorrelation (LSA) at individual localities. LSA statistics serve as useful indicators of local genetic population structure. We applied this method to 15 allele frequencies from 43 villages of a South American tribe, the Yanomama. Based on a network of links 相似文献   

Historical Environment Is Reflected in Modern Population Genetics and Biogeography of an Island Endemic Lizard (Xantusia riversiana reticulata)

The restricted distribution and isolation of island endemics often produces unique genetic and phenotypic diversity of conservation interest to management agencies. However, these isolated species, especially those with sensitive life history traits, are at high risk for the adverse effects of genetic drift and habitat degradation by non-native wildlife. Here, we study the population genetic diversity, structure, and stability of a classic “island giant” (Xantusia riversiana, the Island Night Lizard) on San Clemente Island, California following the removal of feral goats. Using DNA microsatellites, we found that this population is reasonably genetically robust despite historical grazing, with similar effective population sizes and genetic diversity metrics across all sampling locations irrespective of habitat type and degree of degradation. However, we also found strong site-specific patterns of genetic variation and low genetic diversity compared to mainland congeners, warranting continued special management as an island endemic. We identify both high and low elevation areas that remain valuable repositories of genetic diversity and provide a case study for other low-dispersal coastal organisms in the face of future climate change.     

Dispersal ability and habitat requirements determine landscape‐level genetic patterns in desert aquatic insects

Species occupying the same geographic range can exhibit remarkably different population structures across the landscape, ranging from highly diversified to panmictic. Given limitations on collecting population‐level data for large numbers of species, ecologists seek to identify proximate organismal traits—such as dispersal ability, habitat preference and life history—that are strong predictors of realized population structure. We examined how dispersal ability and habitat structure affect the regional balance of gene flow and genetic drift within three aquatic insects that represent the range of dispersal abilities and habitat requirements observed in desert stream insect communities. For each species, we tested for linear relationships between genetic distances and geographic distances using Euclidean and landscape‐based metrics of resistance. We found that the moderate‐disperser Mesocapnia arizonensis (Plecoptera: Capniidae) has a strong isolation‐by‐distance pattern, suggesting migration–drift equilibrium. By contrast, population structure in the flightless Abedus herberti (Hemiptera: Belostomatidae) is influenced by genetic drift, while gene flow is the dominant force in the strong‐flying Boreonectes aequinoctialis (Coleoptera: Dytiscidae). The best‐fitting landscape model for M. arizonensis was based on Euclidean distance. Analyses also identified a strong spatial scale‐dependence, where landscape genetic methods only performed well for species that were intermediate in dispersal ability. Our results highlight the fact that when either gene flow or genetic drift dominates in shaping population structure, no detectable relationship between genetic and geographic distances is expected at certain spatial scales. This study provides insight into how gene flow and drift interact at the regional scale for these insects as well as the organisms that share similar habitats and dispersal abilities.     

Campora: a young genetic isolate in South Italy

Genetic isolates have been successfully used in the study of complex traits, mainly because due to their features, they allow a reduction in the complexity of the genetic models underlying the trait. The aim of the present study is to describe the population of Campora, a village in the South of Italy, highlighting its properties of a genetic isolate. Both historical evidence and multi-locus genetic data (genomic and mitochondrial DNA polymorphisms) have been taken into account in the analyses. The extension of linkage disequilibrium (LD) regions has been evaluated on autosomes and on a region of the X chromosome. We defined a study sample population on the basis of the genealogy and exogamy data. We found in this population a few different mitochondrial and Y chromosome haplotypes and we ascertained that, similarly to other isolated populations, in Campora LD extends over wider region compared to large and genetically heterogeneous populations. These findings indicate a conspicuous genetic homogeneity in the genome. Finally, we found evidence for a recent population bottleneck that we propose to interpret as a demographic crisis determined by the plague of the 17th century. Overall our findings demonstrate that Campora displays the genetic characteristics of a young isolate.     

High resolution analysis and phylogenetic network construction using complete mtDNA sequences in sardinian genetic isolates

For mitochondrial phylogenetic analysis, the best result comes from complete sequences. We therefore decided to sequence the entire mitochondrial DNA (mtDNA) (coding and D-loop regions) of 63 individuals selected in 3 small Ogliastra villages, an isolated area of eastern Sardinia: Talana, Urzulei, and Perdasdefogu. We studied at least one individual for each of the most frequent maternal genealogical lineages belonging to haplogroups H, V, J, K, T, U, and X. We found in our 63 samples, 172 and 69 sequence changes in the coding and in the D-loop region, respectively. Thirteen out of 172 sequence changes in the coding region are novel. It is our hypothesis that some of them are characteristic of the Ogliastra region and/or Sardinia. We reconstructed the phylogenetic network of the 63 complete mtDNA sequences for the 3 villages. We also drew a network including a large number of European sequences and calculated various indices of genetic diversity in Ogliastra. It appears that these small populations remained extremely isolated and genetically differentiated compared with other European populations. We also identified in our samples a never previously described subhaplogroup, U5b3, which seems peculiar to the Ogliastra region.     

Demographic history of the trace metal hyperaccumulator Noccaea caerulescens (J. Presl and C. Presl) F. K. Mey. in Western Europe

Noccaea caerulescens (Brassicaceae) is a major pseudometallophyte model for the investigation of the genetics and evolution of metal hyperaccumulation in plants. We studied the population genetics and demographic history of this species to advance the understanding of among‐population differences in metal hyperaccumulation and tolerance abilities. Sampling of seven to 30 plants was carried out in 62 sites in Western Europe. Genotyping was carried out using a combination of new chloroplast and nuclear neutral markers. A strong genetic structure was detected, allowing the definition of three genetic subunits. Subunits showed a good geographic coherence. Accordingly, distant metallicolous populations generally belonged to distinct subunits. Approximate Bayesian computation analysis of demographic scenarios among subunits further supported a primary isolation of populations from the southern Massif Central prior to last glacial maximum, whereas northern populations may have derived during postglacial recolonization events. Estimated divergence times among subunits were rather recent in comparison with the species history, but certainly before the establishment of anthropogenic metalliferous sites. Our results suggest that the large‐scale genetic structure of N. caerulescens populations pre‐existed to the local adaptation to metalliferous sites. The population structure of quantitative variation for metal‐related adaptive traits must have established independently in isolated gene pools. However, features of the most divergent genetic unit (e.g. extreme levels of Cd accumulation observed in previous studies) question the putative relationships between adaptive evolution of metal‐related traits and subunits isolation. Finally, admixture signals among distant metallicolous populations suggest a putative role of human activities in facilitating long‐distance genetic exchanges.     

Multivariate classification of human populations. I. Allocation of Yanomama indians to villages.

A set of 12 anthropometric measures and six genetic traits, available for 520 Yanomama Indians from 19 villages in nine clusters, were used to allocate individuals to villages. On the basis of anthropometrics alone, 36% of the individuals were allocated to the right village and 60% to the right cluster. On the basis of genetic traits alone, 16% were allocated to the right village and 26% to the right cluster. A combination of all 18 characters yielded 41% allocation to the right village and 63% to the right cluster. Of the 924 possible combinations of six anthropometric measures, only one provided poorer resolution than did the six genetic traits. We explain the better resolution of the anthropometric traits by noting that the anthropometric traits are not totally heritable and that genetic traits are not continuously distributed. Randomization studies indicated that all of the observed correct-allocation fractions are far in excess of random expectation. We infer that the village phenotype distributions overlap only partially, and that they represent real and substantial population differentiation.     

Sky island bird populations isolated by ancient genetic barriers are characterized by different song traits than those isolated by recent deforestation

Various mechanisms of isolation can structure populations and result in cultural and genetic differentiation. Similar to genetic markers, for songbirds, culturally transmitted sexual signals such as breeding song can be used as a measure of differentiation as songs can also be impacted by geographic isolation resulting in population‐level differences in song structure. Several studies have found differences in song structure either across ancient geographic barriers or across contemporary habitat barriers owing to deforestation. However, very few studies have examined the effect of both ancient barriers and recent deforestation in the same system. In this study, we examined the geographic variation in song structure across six populations of the White‐bellied Shortwing, a threatened and endemic songbird species complex found on isolated mountaintops or “sky islands” of the Western Ghats. While some sky islands in the system are isolated by ancient valleys, others are separated by deforestation. We examined 14 frequency and temporal spectral traits and two syntax traits from 835 songs of 38 individuals across the six populations. We identified three major song clusters based on a discriminant model of spectral traits, degree of similarity of syntax features, as well as responses of birds to opportunistic playback. However, some traits like complex vocal mechanisms (CVM), relating to the use of syrinxes, clearly differentiated both ancient and recently fragmented populations. We suggest that CVMs may have a cultural basis and can be used to identify culturally isolated populations that cannot be differentiated using genetic markers or commonly used frequency‐based song traits. Our results demonstrate the use of bird songs to reconstruct phylogenetic groups and impacts of habitat fragmentation even in complex scenarios of historic and contemporary isolation.     

Distribution, endogamy, and isolation of Malas of Chittoor district, Andhra Pradesh, India

This study investigates the distribution of Malas, a scheduled caste population of Andhra Pradesh, their isolates in different eco-cultural zones, and their endogamy and isolation by marriage district. The Malas, formerly "untouchables," occupy the lowest status in the Hindu hierarchy. The sample consists of 10% of Malas from 10% of the villages in 2 taluks of Chitoor district of Andhra Pradesh. 6 Mala populations--Tangala, Maladasari, Pakanati, Rampala, Murikinati, and Bommanati--live in the area. These populations show a regionality in their distribution, with very little overlapping even when 2 populations inhabit the same village. Of 885 marriages in the 6 endogamous populations, all but 3 have been contracted between individuals belonging to the same Mala group. The 3 exogamous marriages took place between Mala men and women from another caste. Such small exceptions to the general rule do not mean that the Mala populations are not breeding isolates; these 6 populations satisfy Wright's island model. The high incidence of matings between closely related populations also contributes to their genetic and breeding isolation. Consanguineous marriages range from 26.76 to 38.75%. The distance between the birth place of spouses in miles, called marriage distance, shows a range from 7.72 to 15.71 miles. Lower values mean higher population densities. Groups within each population are isolated by distance and form small overlapping Mendelian populations, approaching a stepping stone model with continuous variation of genetic traits between adjacent groups of people.     

Physical and ecological isolation contribute to maintain genetic differentiation between fire salamander subspecies

Landscape features shape patterns of gene flow among populations, ultimately determining where taxa lay along the continuum between panmixia to complete reproductive isolation. Gene flow can be restricted, leading to population differentiation in two non-exclusive ways: “physical isolation”, in which geographic distance in combination with the landscape features restricts movement of individuals promoting genetic drift, and “ecological isolation”, in which adaptive mechanisms constrain gene flow between different environments via divergent natural selection. In central Iberia, two fire salamander subspecies occur in parapatry across elevation gradients along the Iberian Central System mountains, while in the adjacent Montes de Toledo Region only one of them occurs. By integrating population and landscape genetic analyses, we show a ubiquitous role of physical isolation between and within mountain ranges, with unsuitable landscapes increasing differentiation between populations. However, across the Iberian Central System, we found strong support for a significant contribution of ecological isolation, with low genetic differentiation in environmentally homogeneous areas, but high differentiation across sharp transitions in precipitation seasonality. These patterns are consistent with a significant contribution of ecological isolation in restricting gene flow among subspecies. Overall, our results suggest that ecological divergence contributes to reduce genetic admixture, creating an opportunity for lineages to follow distinct evolutionary trajectories.Subject terms: Evolutionary ecology, Speciation     

Evidence for Geographic Isolation and Signs of Endemism within a Protistan Morphospecies

The possible existence of endemism among microorganisms resulting from and preserved by geographic isolation is one of the most controversial topics in microbial ecology. We isolated 31 strains of “Spumella-like” flagellates from remote sampling sites from all continents, including Antarctica. These and another 23 isolates from a former study were characterized morphologically and by small-subunit rRNA gene sequence analysis and tested for the maximum temperature tolerance. Only a minority of the Spumella morpho- and phylotypes from the geographically isolated Antarctic continent follow the worldwide trend of a linear correlation between ambient (air) temperature during strain isolation and heat tolerance of the isolates. A high percentage of the Antarctic isolates, but none of the isolates from locations on all other continents, were obligate psychrophilic, although some of the latter were isolated at low ambient temperatures. The drastic deviation of Antarctic representatives of Spumella from the global trend of temperature adaptation of this morphospecies provides strong evidence for geographic transport restriction of a microorganism; i.e., Antarctic protistan communities are less influenced by transport of protists to and from the Antarctic continent than by local adaptation, a subtle form of endemism.     

Genetic variation in walnuts (Juglans regia and J. sigillata; Juglandaceae): Species distinctions, human impacts, and the conservation of agrobiodiversity in Yunnan, China

Walnuts are a major crop of many countries and mostly cultivated in large-scale plantations with few cultivars. Landraces provide important genetic reservoirs; thus, understanding factors influencing the geographic distribution of genetic variation in crop resources is a fundamental goal of agrobiodiversity conservation. Here, we investigated the role of human settlements and kinship on genetic variation and population structure of two walnut species: Juglans regia, an introduced species widely cultivated for its nuts, and J. sigillata, a native species cultivated locally in Yunnan. The objectives of this study were to characterize sympatric populations of J. regia and J. sigillata using 14 molecular markers and evaluate the role of Tibetan villages and kin groups (related households) on genotypic variation and population structure of J. regia and J. sigillata. Our results based on 220 walnut trees from six Tibetan villages show that although J. regia and J. sigillata are morphologically distinct, the two species are indistinguishable based on microsatellite data. Despite the lack of interspecific differences, AMOVAs partitioned among villages (5.41%, P = 0.0068) and kin groups within villages (3.34%, P = 0.0068) showed significant genetic variation. These findings suggest that village environments and familial relationships are factors contributing to the geographic structure of genetic variation in Tibetan walnuts.     

Clonal Evolution of Enterocytozoon bieneusi Populations in Swine and Genetic Differentiation in Subpopulations between Isolates from Swine and Humans

Enterocytozoon bieneusi is a widespread parasite with high genetic diversity among hosts. Its natural reservoir remains elusive and data on population structure are available only in isolates from primates. Here we describe a population genetic study of 101 E. bieneusi isolates from pigs using sequence analysis of the ribosomal internal transcribed spacer (ITS) and four mini- and microsatellite markers. The presence of strong linkage disequilibrium (LD) and limited genetic recombination indicated a clonal structure for the population. Bayesian inference of phylogeny, structural analysis, and principal coordinates analysis separated the overall population into three subpopulations (SP3 to SP5) with genetic segregation of the isolates at some geographic level. Comparative analysis showed the differentiation of SP3 to SP5 from the two known E. bieneusi subpopulations (SP1 and SP2) from primates. The placement of a human E. bieneusi isolate in pig subpopulation SP4 supported the zoonotic potential of some E. bieneusi isolates. Network analysis showed directed evolution of SP5 to SP3/SP4 and SP1 to SP2. The high LD and low number of inferred recombination events are consistent with the possibility of host adaptation in SP2, SP3, and SP4. In contrast, the reduced LD and high genetic diversity in SP1 and SP5 might be results of broad host range and adaptation to new host environment. The data provide evidence of the potential occurrence of host adaptation in some of E. bieneusi isolates that belong to the zoonotic ITS Group 1.     

Phylogeographical analysis of an estuarine fish, Salanx ariakensis (Osmeridae: Salanginae) in the north-western Pacific

This study extended the geographic coverage of a previous study to explore population genetic structure and demographic history in the Ariake icefish Salanx ariakensis from three populations of continental coastlines and one island population in the north-western Pacific based on a partial sequence of the mitochondrial cytochrome b gene. The S. ariakensis showed high genetic diversity and strong genetic structure. Phylogenetic analysis showed a shallow gene tree with no clear phylogeographical structure. Contiguous range expansion and restricted gene flow were inferred to be main population events by nested-clade analysis. Significant genetic differentiations between populations could be attributable to negligible gene flow by coalescent analysis. High nucleotide diversity of each population was due to geographic mixing of heterogenous haplotypes during lowering sea levels of the Pleistocene. These findings indicate that cycles of geographic isolation and secondary contact happened in the Pleistocene glacial–interglacial cycles shaping genetic structure and population demography of S. ariakensis .