Population genetic structure and demography analysis in palinurus elephas populations from atlantic and mediterranean waters

Report for the scientific sojourn at the University of Reading, United Kingdom, from January until May 2008. The main objectives have been firstly to infer population structure and parameters in demographic models using a total of 13 microsatellite loci for genotyping approximately 30 individuals per population in 10 Palinurus elephas populations both from Mediterranean and Atlantic waters. Secondly, developing statistical methods to identify discrepant loci, possibly under selection and implement those methods using the R software environment. It is important to consider that the calculation of the probability distribution of the demographic and mutational parameters for a full genetic data set is numerically difficult for complex demographic history (Stephens 2003). The Approximate Bayesian Computation (ABC), based on summary statistics to infer posterior distributions of variable parameters without explicit likelihood calculations, can surmount this difficulty. This would allow to gather information on different demographic prior values (i.e. effective population sizes, migration rate, microsatellite mutation rate, mutational processes) and assay the sensitivity of inferences to demographic priors by assuming different priors.

Y-chromosome diversity in native mexicans reveals continental transition of genetic structure in the Americas

The genetic characterization of Native Mexicans is important to understand multiethnic based features influencing the medical genetics of present Mexican populations, as well as to the reconstruct the peopling of the Americas. We describe the Y-chromosome genetic diversity of 197 Native Mexicans from 11 populations and 1,044 individuals from 44 Native American populations after combining with publicly available data. We found extensive heterogeneity among Native Mexican populations and ample segregation of Q-M242* (46%) and Q-M3 (54%) haplogroups within Mexico. The northernmost sampled populations falling outside Mesoamerica (Pima and Tarahumara) showed a clear differentiation with respect to the other populations, which is in agreement with previous results from mtDNA lineages. However, our results point toward a complex genetic makeup of Native Mexicans whose maternal and paternal lineages reveal different narratives of their population history, with sex-biased continental contributions and different admixture proportions. At a continental scale, we found that Arctic populations and the northernmost groups from North America cluster together, but we did not find a clear differentiation within Mesoamerica and the rest of the continent, which coupled with the fact that the majority of individuals from Central and South American samples are restricted to the Q-M3 branch, supports the notion that most Native Americans from Mesoamerica southwards are descendants from a single wave of migration. This observation is compatible with the idea that present day Mexico might have constituted an area of transition in the diversification of paternal lineages during the colonization of the Americas.

A genomic analysis identifies a novel component in the genetic structure of sub-Saharan African populations

Studies of large sets of SNP data have proven to be a powerful tool in the analysis of the genetic structure of human populations. In this work, we analyze genotyping data for 2,841 SNPs in 12 Sub-Saharan African populations, including a previously unsampled region of south-eastern Africa (Mozambique). We show that robust results in a world-wide perspective can be obtained when analyzing only 1,000 SNPs. Our main results both confirm the results of previous studies, and show new and interesting features in Sub-Saharan African genetic complexity. There is a strong differentiation of Nilo-Saharans, much beyond what would be expected by geography. Hunter-gatherer populations (Khoisan and Pygmies) show a clear distinctiveness with very intrinsic Pygmy (and not only Khoisan) genetic features. Populations of the West Africa present an unexpected similarity among them, possibly the result of a population expansion. Finally, we find a strong differentiation of the south-eastern Bantu population from Mozambique, which suggests an assimilation of a pre-Bantu substrate by Bantu speakers in the region.

Time and space: genetic structure of the common sea urchin paracentrotus lividus in the south-eastern iberian coast

Spatio-temporal variability in settlement and recruitment, high mortality during the first life-history stages, and selection may determine the genetic structure of cohorts of long-lived marine invertebrates at small scales. We conducted a spatial and temporal analysis of the common Mediterranean Sea urchin Paracentrotus lividus to determine the genetic structure of cohorts at different scales. In Tossa de Mar (NW Mediterranean), recruitment was followed over 5 consecutive springs (2006-2010). In spring 2008, recruits and two-year-old individuals were collected at 6 locations along East and South Iberian coasts separated from 200 to over 1,100 km. All cohorts presented a high genetic diversity based on a fragment of mtCOI. Our results showed a marked genetic homogeneity in the temporal monitoring and a low degree of spatial structure in 2006. In 2008, coupled with an abnormality in the usual circulation patterns in the area, the genetic structure of the southern populations studied changed markedly, with arrival of many private haplotypes. This fact highlights the importance of point events in renewing the genetic makeup of populations, which can only be detected through analysis of the cohort structure coupling temporal and spatial perspectives.

Geographic patterns of genetic variation in a broadly distributed marine vertebrate: new insights into loggerhead turtle stock structure from expanded mitochondrial DNA sequences

Previous genetic studies have demonstrated that natal homing shapes the stock structure of marine turtle nesting populations. However, widespread sharing of common haplotypes based on short segments of the mitochondrial control region often limits resolution of the demographic connectivity of populations. Recent studies employing longer control region sequences to resolve haplotype sharing have focused on regional assessments of genetic structure and phylogeography. Here we synthesize available control region sequences for loggerhead turtles from the Mediterranean Sea, Atlantic, and western Indian Ocean basins. These data represent six of the nine globally significant regional management units (RMUs) for the species and include novel sequence data from Brazil, Cape Verde, South Africa and Oman. Genetic tests of differentiation among 42 rookeries represented by short sequences (380 bp haplotypes from 3,486 samples) and 40 rookeries represented by long sequences (~800 bp haplotypes from 3,434 samples) supported the distinction of the six RMUs analyzed as well as recognition of at least 18 demographically independent management units (MUs) with respect to female natal homing. A total of 59 haplotypes were resolved. These haplotypes belonged to two highly divergent global lineages, with haplogroup I represented primarily by CC-A1, CC-A4, and CC-A11 variants and haplogroup II represented by CC-A2 and derived variants. Geographic distribution patterns of haplogroup II haplotypes and the nested position of CC-A11.6 from Oman among the Atlantic haplotypes invoke recent colonization of the Indian Ocean from the Atlantic for both global lineages. The haplotypes we confirmed for western Indian Ocean RMUs allow reinterpretation of previous mixed stock analysis and further suggest that contemporary migratory connectivity between the Indian and Atlantic Oceans occurs on a broader scale than previously hypothesized. This study represents a valuable model for conducting comprehensive international cooperative data management and research in marine ecology.

Microsatellite markers reveal shallow genetic differentiation between cohorts of the common sea urchin Paracentrotus lividus (Lamarck) in northwest Mediterranean

Temporal variability was studied in the common sea urchin Paracentrotus lividus through the analysis of the genetic composition of three yearly cohorts sampled over two consecutive springs in a locality in northwestern Mediterranean. Individuals were aged using growth ring patterns observed in tests and samples were genotyped for five microsatellite loci. No reduction of genetic diversity was observed relative to a sample of the adult population from the same location or within cohorts across years. FST and amova results indicated that the differentiation between cohorts is rather shallow and not significant, as most variability is found within cohorts and within individuals. This mild differentiation translated into estimates of effective population size of 90100 individuals. When the observed excess of homozygotes was taken into account, the estimate of the average number of breeders increased to c. 300 individuals. Given our restricted sampling area and the known small-scale heterogeneity in recruitment in this species, our results suggest that at stretches of a few kilometres of shoreline, large numbers of progenitors are likely to contribute to the larval pool at each reproduction event. Intercohort variation in our samples is six times smaller than spatial variation between adults of four localities in the western Mediterranean. Our results indicate that, notwithstanding the stochastic events that take place during the long planktonic phase and during the settlement and recruitment processes, reproductive success in this species is high enough to produce cohorts genetically diverse and with little differentiation between them. Further research is needed before the link between genetic structure and underlying physical and biological processes can be well established.

Population structure in a highly pelagic seabird, the Cory's shearwater (Calonectris diomedea): an examination of genetics, morphology and ecology

Increasing evidence suggests oceanic traits may play a key role in the genetic structuring of marine organisms. Whereas genetic breaks in the open ocean are well known in fishes and marine invertebrates, the importance of marine habitat characteristics in seabirds remains less certain. We investigated the role of oceanic transitions versus population genetic processes in driving population differentiation in a highly vagile seabird, the Cory"s shearwater, combining molecular, morphological and ecological data from 27 breeding colonies distributed across the Mediterranean (Calonectris diomedea diomedea) and the Atlantic (C. d. borealis). Genetic and biometric analyses showed a clear differentiation between Atlantic and Mediterranean Cory"s shearwaters. Ringing-recovery data indicated high site fidelity of the species, but we found some cases of dispersal among neighbouring breeding sites (<300 km) and a few long distance movements (>1000 km) within and between each basin. In agreement with this, comparison of phenotypic and genetic data revealed both current and historical dispersal events. Within each region, we did not detect any genetic substructure among archipelagos in the Atlantic, but we found a slight genetic differentiation between western and eastern breeding colonies in the Mediterranean. Accordingly, gene flow estimates suggested substantial dispersal among colonies within basins. Overall, genetic structure of the Cory"s shearwater matches main oceanographic breaks (Almería-Oran Oceanic Front and Siculo-Tunisian Strait), but spatial analyses suggest that patterns of genetic differentiation are better explained by geographic rather than oceanographic distances. In line with previous studies, genetic, phenotypic and ecological evidence supported the separation of Atlantic and Mediterranean forms, suggesting the 2 taxa should be regarded as different species.

Combining genetic and demographic data for the conservation of a Mediterranean habitat-forming species

The integration of ecological and evolutionary data is highly valuable for conservation planning. However, it has been rarely used in the marine realm, where the adequate design of marine protected areas (MPAs) is urgently needed. Here, we examined the interacting processes underlying the patterns of genetic structure and demographic strucuture of a highly vulnerable Mediterranean habitat-forming species (i.e. Paramuricea clavata (Risso, 1826)), with particular emphasis on the processes of contemporary dispersal, genetic drift, and colonization of a new population. Isolation by distance and genetic discontinuities were found, and three genetic clusters were detected; each submitted to variations in the relative impact of drift and gene flow. No founder effect was found in the new population. The interplay of ecology and evolution revealed that drift is strongly impacting the smallest, most isolated populations, where partial mortality of individuals was highest. Moreover, the eco-evolutionary analyses entailed important conservation implications for P. clavata. Our study supports the inclusion of habitat-forming organisms in the design of MPAs and highlights the need to account for genetic drift in the development of MPAs. Moreover, it reinforces the importance of integrating genetic and demographic data in marine conservation.

Origen i diversificació de les comunitats de rèptils dels deserts del Nord d’Àfrica i Aràbia: els paràsits com a eines filogenètiques

Els deserts constitueixen aproximadament un terç de la superfície terrestre i estan caracteritzats per la seva aridesa extrema. Els Sàhara, la zona àrida més gran del Nord d’Àfrica és, amb diferencia, el desert més gran del mon. Ocupa una extensió de més de 9 milions de kilòmetres quadrats, expandint-se 5500 km a través del nord d’Àfrica, des de l’oceà Atlàntic fins el mar Roig. El desert d’Aràbia es troba a l’est del Sàhara i és aproximadament una octava part més petit. Tot i el seu clima tan extrem, les flores i faunes dels deserts acostumen a ser relativament riques, el que fa que ens preguntem com aquestes biotes han estat adquirides i com es mantenen. En aquest projecte utilitzem els rèptils com a model d’estudi, un dels habitants més comuns dels deserts. El propòsit d’aquest projecte era utilitzar filogènies moleculars de diversos taxons de rèptils dels deserts del nord d’Àfrica i Aràbia per respondre Quan i Com els deserts han adquirit la seva fauna de rèptils endèmica (origen i diversificació), i de quina manera aquestes faunes s’han mantingut fins ara (adaptació). Però més enllà de les dades morfològiques i les eines de filogènia molecular, els paràsits representen una alternativa excitant i innovadora dins del camp de la biologia evolutiva. En aquest sentit, aquest projecte de beca proposava anar més lluny i utilitzar els paràsits com a eina (biological tags) per millor entendre l’historia evolutiva dels seus hostes en base a una aproximació biogeogràfica i co-evolutiva. Els objectius durant aquest primer any de projecte han sigut: 1) Estudiar l’origen i diversificació de la fauna de rèptils endèmica dels deserts del Sàhara i Aràbia, 2) caracteritzar la comunitat d’ecto- haemoparàsits de les espècies d’estudi, 3) posar a punt un estudi comparatiu de la filogeografia i estructuració genètica d’hostes i paràsits.

O-chromosome lethal frequencies in Serbian and Montenegrin Drosophila subobscura populations

Lethal chromosomal frequencies were obtained from three Drosophila subobscura samples from the Mt. Avala (Serbia) population in September 2003 (0.218), June 2004 (0.204) and September 2004 (0.250). These values and those from other Balkan populations studied previously (Petnica, Kamariste, Zanjic and Djerdap) were used to analyze the possible effect of population, year, month and altitude above sea level on lethal chromosomal frequencies. According to ANOVAS no effect were observed. Furthermore, the lethal frequencies of the Balkan populations did not vary according to latitude. This is probably due to the relative proximity and high gene flow between these populations. From a joint study of all the Palearctic D. subobscura populations so far analyzed, it can be deduced that the Balkan populations are located in the central area of the species distribution. Finally, it seems that lethal chromosomal frequencies are a consequence of the genetic structure of the populations.

O-chromosome lethal frequencies in Serbian and Montenegrin Drosophila subobscura populations

Lethal chromosomal frequencies were obtained from three Drosophila subobscura samples from the Mt. Avala (Serbia) population in September 2003 (0.218), June 2004 (0.204) and September 2004 (0.250). These values and those from other Balkan populations studied previously (Petnica, Kamariste, Zanjic and Djerdap) were used to analyze the possible effect of population, year, month and altitude above sea level on lethal chromosomal frequencies. According to ANOVAS no effect were observed. Furthermore, the lethal frequencies of the Balkan populations did not vary according to latitude. This is probably due to the relative proximity and high gene flow between these populations. From a joint study of all the Palearctic D. subobscura populations so far analyzed, it can be deduced that the Balkan populations are located in the central area of the species distribution. Finally, it seems that lethal chromosomal frequencies are a consequence of the genetic structure of the populations.

Natural or naturalized? phylogeography suggests that the abundant sea urchin arbacia lixula Is a recent Colonizer of the Mediterranean

We present the global phylogeography of the black sea urchin Arbacia lixula, an amphi-Atlantic echinoid with potential to strongly impact shallow rocky ecosystems. Sequences of the mitochondrial cytochrome c oxidase gene of 604 specimens from 24 localities were obtained, covering most of the distribution area of the species, including the Mediterranean and both shores of the Atlantic. Genetic diversity measures, phylogeographic patterns, demographic parameters and population differentiation were analysed. We found high haplotype diversity but relatively low nucleotide diversity, with 176 haplotypes grouped within three haplogroups: one is shared between Eastern Atlantic (including Mediterranean) and Brazilian populations, the second is found in Eastern Atlantic and the Mediterranean and the third is exclusively from Brazil. Significant genetic differentiation was found between Brazilian, Eastern Atlantic and Mediterranean regions, but no differentiation was found among Mediterranean sub-basins or among Eastern Atlantic sub-regions. The star-shaped topology of the haplotype network and the unimodal mismatch distributions of Mediterranean and Eastern Atlantic samples suggest that these populations have suffered very recent demographic expansions. These expansions could be dated 94-205 kya in the Mediterranean, and 31-67 kya in the Eastern Atlantic. In contrast, Brazilian populations did not show any signature of population expansion. Our results indicate that all populations of A. lixula constitute a single species. The Brazilian populations probably diverged from an Eastern Atlantic stock. The present-day genetic structure of the species in Eastern Atlantic and the Mediterranean is shaped by very recent demographic processes. Our results support the view (backed by the lack of fossil record) that A. lixula is a recent thermophilous colonizer which spread throughout the Mediterranean during a warm period of the Pleistocene, probably during the last interglacial. Implications for the possible future impact of A. lixula on shallow Mediterranean ecosystems in the context of global warming trends must be considered.

The whereabouts of an ancient wonderer: Global phylogeography of the solitary ascidian Styela plicata

Genetic tools have greatly aided in tracing the sources and colonization history of introduced species. However, recurrent introductions and repeated shuffling of populations may have blurred some of the genetic signals left by ancient introductions. Styela plicata is a solitary ascidian distributed worldwide. Although its origin remains unclear, this species is believed to have spread worldwide by travelling on ship's hulls. The goals of this study were to infer the genetic structure and global phylogeography of S. plicata and to look for present-day and historical genetic patterns. Two genetic markers were used: a fragment of the mitochondrial gene Cytochrome Oxidase subunit I (COI) and a fragment of the nuclear gene Adenine Nucleotide Transporter/ADP-ATP Translocase (ANT). A total of 368 individuals for COI and 315 for ANT were sequenced from 17 locations worldwide. The levels of gene diversity were moderate for COI to high for ANT. The Mediterranean populations showed the least diversity and allelic richness for both markers, while the Indian, Atlantic and Pacific Oceans had the highest gene and nucleotide diversities. Network and phylogenetic analyses with COI and ANT revealed two groups of alleles separated by 15 and 4 mutational steps, respectively. The existence of different lineages suggested an ancient population split. However, the geographic distributions of these groups did not show any consistent pattern, indicating different phylogeographic histories for each gene. Genetic divergence was significant for many population-pairs irrespective of the geographic distance among them. Stochastic introduction events are reflected in the uneven distribution of COI and ANT allele frequencies and groups among many populations. Our results confirmed that S. plicata has been present in all studied oceans for a long time, and that recurrent colonization events and occasional shuffling among populations have determined the actual genetic structure of this species.

Linguistic and maternal genetic diversity are not correlated in native mexicans

Mesoamerica, defined as the broad linguistic and cultural area from middle southern Mexico to Costa Rica, might have played a pivotal role during the colonization of theAmerican continent. It has been suggested that the Mesoamerican isthmus could have played an important role in severely restricting prehistorically gene flow between North and SouthAmerica. Although the Native American component has been already described in admixedMexican populations, few studies have been carried out in native Mexican populations. In thisstudy we present mitochondrial DNA (mtDNA) sequence data for the first hypervariable region (HVR-I) in 477 unrelated individuals belonging to eleven different native populations from Mexico. Almost all the Native Mexican mtDNAs could be classified into the four pan-Amerindian haplogroups (A2, B2, C1 and D1); only three of them could be allocated to the rare Native American lineage D4h3. Their haplogroup phylogenies are clearly star-like, as expected from relatively young populations that have experienced diverse episodes of genetic drift (e.g. extensive isolation, genetic drift and founder effects) and posterior population expansions. In agreement with this observation is the fact that Native Mexican populations show a high degree of heterogeneity in their patterns of haplogroup frequencies. HaplogroupX2a was absent in our samples, supporting previous observations where this clade was only detected in the American northernmost areas. The search for identical sequences in the American continent shows that, although Native Mexican populations seem to show a closer relationship to North American populations, they cannot be related to a single geographical region within the continent. Finally, we did not find significant population structure on the maternal lineages when considering the four main and distinct linguistic groups represented in our Mexican samples (Oto-Manguean, Uto-Aztecan, Tarascan, and Mayan), suggesting that genetic divergence predates linguistic diversification in Mexico.

Origin and genetic differentiation of Three Mexican Native Groups (Purépechas, Triquis and Mayas): contribution of CODIS-STRs to the history of human populations of Mesoamerica

BACKGROUND: CODIS-STRs in Native Mexican groups have rarely been analysed for human identification and anthropological purposes. AIM:To analyse the genetic relationships and population structure among three Native Mexican groups from Mesoamerica.SUBJECTS AND METHODS: 531 unrelated Native individuals from Mexico were PCR-typed for 15 and 9 autosomal STRs (Identifiler™ and Profiler™ kits, respectively), including five population samples: Purépechas (Mountain, Valley and Lake), Triquis and Yucatec Mayas. Previously published STR data were included in the analyses. RESULTS:Allele frequencies and statistical parameters of forensic importance were estimated by population. The majority of Native groups were not differentiated pairwise, excepting Triquis and Purépechas, which was attributable to their relative geographic and cultural isolation. Although Mayas, Triquis and Purépechas-Mountain presented the highest number of private alleles, suggesting recurrent gene flow, the elevated differentiation of Triquis indicates a different origin of this gene flow. Interestingly, Huastecos and Mayas were not differentiated, which is in agreement with the archaeological hypothesis that Huastecos represent an ancestral Maya group. Interpopulation variability was greater in Natives than in Mestizos, both significant.CONCLUSION: Although results suggest that European admixture has increased the similarity between Native Mexican groups, the differentiation and inconsistent clustering by language or geography stresses the importance of serial founder effect and/or genetic drift in showing their present genetic relationships.