Search

The biological behavior of the Y chromosome, which is paternally inherited, implies that males sharing the same surname may also share a similar Y chromosome. However, socio-cultural factors, such as polyphyletism, non-paternity, adoption, or matrilineal surname transmission, may prevent the joint transmission of the surname and the Y chromosome. By genotyping 17 Y-STRs and 68 SNPs in ~2500 male samples that each carried one of the 50 selected Catalan surnames, we could determine sets of descendants of a common ancestor, the population of origin of the common ancestor, and the date when such a common ancestor lived. Haplotype diversity was positively correlated with surname frequency, that is, rarer surnames showed the strongest signals of coancestry. Introgression rates of Y chromosomes into a surname by non-paternity, adoption, and transmission of the maternal surname were estimated at 1.5-2.6% per generation, with some local variation. Average ages for the founders of the surnames were estimated at ~500 years, suggesting a delay between the origin of surnames (twelfth and thirteenth centuries) and the systematization of their paternal transmission. We have found that, in general, a foreign etymology for a surname does not often result in a non-indigenous origin of surname founders; however, bearers of some surnames with an Arabic etymology show an excess of North African haplotypes. Finally, we estimate that surname prediction from a Y-chromosome haplotype, which may have interesting forensic applications, has a ~60% sensitivity but a 17% false discovery rate.

The demographic history of human populations in North Africa has been characterized by complex processes of admixture and isolation that have modeled its current gene pool. Diverse genetic ancestral components with different origins (autochthonous, European, Middle Eastern, and sub-Saharan) and genetic heterogeneity in the region have been described. In this complex genetic landscape, Algeria, the largest country in Africa, has been poorly covered, with most of the studies using a single Algerian sample. In order to evaluate the genetic heterogeneity of Algeria, Y-chromosome, mtDNA and autosomal genome-wide makers have been analyzed in several Berber- and Arab-speaking groups. Our results show that the genetic heterogeneity found in Algeria is not correlated with geography or linguistics, challenging the idea of Berber groups being genetically isolated and Arab groups open to gene flow. In addition, we have found that external sources of gene flow into North Africa have been carried more often by females than males, while the North African autochthonous component is more frequent in paternally transmitted genome regions. Our results highlight the different demographic history revealed by different markers and urge to be cautious when deriving general conclusions from partial genomic information or from single samples as representatives of the total population of a region.

While most patients affected by the influenza A(H1N1) pandemic experienced mild symptoms, a small fraction required hospitalization, often without concomitant factors that could explain such a severe course. We hypothesize that host genetic factors could contribute to aggravate the disease. To test this hypothesis, we compared the allele frequencies of 547,296 genome-wide single nucleotide polymorphisms (SNPs) between 49 severe and 107 mild confirmed influenza A cases, as well as against a general population sample of 549 individuals. When comparing severe vs. mild influenza A cases, only one SNP was close to the conventional p = 5×10-8. This SNP, rs28454025, sits in an intron of the GSK233 gene, which is involved in a neural development, but seems not to have any connections with immunological or inflammatory functions. Indirectly, a previous association reported with CD55 was replicated. Although sample sizes are low, we show that the statistical power in our design was sufficient to detect highly-penetrant, quasi-Mendelian genetic factors. Hence, and assuming that rs28454025 is likely to be a false positive, no major genetic factor was detected that could explain poor influenza A course.

Parkinson's disease (PD) can be divided into familial (Mendelian) and sporadic forms. A number of causal genes have been discovered for the Mendelian form, which constitutes 10-20% of the total cases. Genome-wide association studies have successfully uncovered a number of susceptibility loci for sporadic cases but those only explain a small fraction (6-7%) of PD heritability. It has been observed that some genes that confer susceptibility to PD through common risk variants also contain rare causing mutations for the Mendelian forms of the disease. These results suggest a possible functional link between Mendelian and sporadic PD and led us to investigate the role that rare and low-frequency variants could have on the sporadic form. Through a targeting approach, we have resequenced at 49× coverage the exons and regulatory regions of 38 genes (including Mendelian and susceptibility PD genes) in 249 sporadic PD patients and 145 unrelated controls of European origin. Unlike susceptibility genes, Mendelian genes show a clear general enrichment of rare functional variants in PD cases, observed directly as well as with Tajima's D statistic and several collapsing methods. Our findings suggest that rare variation on PD Mendelian genes may have a role in the sporadic forms of the disease., (© The Author 2014. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

DNA methylation is one of the main epigenetic mechanisms that can regulate gene expression and is an important means for creating phenotypic variation. In the present study, we performed methylation profiling of 2 candidate genes for personality traits, namely DRD4 and SERT, in the great tit Parus major to ascertain whether personality traits and behavior within different habitats have evolved with the aid of epigenetic variation. We applied bisulphite PCR and strand-specific sequencing to determine the methylation profile of the CpG dinucleotides in the DRD4 and SERT promoters and also in the CpG island overlapping DRD4 exon 3. Furthermore, we performed pyrosequencing to quantify the total methylation levels at each CpG location. Our results indicated that methylation was ∼1-4% higher in urban than in forest birds, for all loci and tissues analyzed, suggesting that this epigenetic modification is influenced by environmental conditions. Screening of genomic DNA sequence revealed that the SERT promoter is CpG poor region. The methylation at a single CpG dinucleotide located 288 bp from the transcription start site was related to exploration score in urban birds. In addition, the genotypes of the SERT polymorphism SNP234 located within the minimal promoter were significantly correlated with novelty seeking behavior in captivity, with the allele increasing this behavior being more frequent in urban birds. As a conclusion, it seems that both genetic and methylation variability of the SERT gene have an important role in shaping personality traits in great tits, whereas genetic and methylation variation at the DRD4 gene is not strongly involved in behavior and personality traits.

The large amount of DNA needed to prepare a library in next generation sequencing protocols hinders direct sequencing of small DNA samples. This limitation is usually overcome by the enrichment of such samples with whole genome amplification (WGA), mostly by multiple displacement amplification (MDA) based on φ29 polymerase. However, this technique can be biased by the GC content of the sample and is prone to the development of chimeras as well as contamination during enrichment, which contributes to undesired noise during sequence data analysis, and also hampers the proper functional and/or taxonomic assignments. An alternative to MDA is direct DNA sequencing (DS), which represents the theoretical gold standard in genome sequencing. In this work, we explore the possibility of sequencing the genome of Escherichia coli fs 24 from the minimum number of DNA molecules required for pyrosequencing, according to the notion of one-bead-one-molecule. Using an optimized protocol for DS, we constructed a shotgun library containing the minimum number of DNA molecules needed to fill a selected region of a picotiterplate. We gathered most of the reference genome extension with uniform coverage. We compared the DS method with MDA applied to the same amount of starting DNA. As expected, MDA yielded a sparse and biased read distribution, with a very high amount of unassigned and unspecific DNA amplifications. The optimized DS protocol allows unbiased sequencing to be performed from samples with a very small amount of DNA.

Background: Metagenomics is the genomic study of uncultured environmental samples, which has been greatly facilitated by the advent of shotgun-sequencing technologies. One of the main focuses of metagenomics is the discovery of previously uncultured microorganisms, which makes the assignment of sequences to a particular taxon a challenge and a crucial step. Recently, several methods have been developed to perform this task, based on different methodologies such as sequence composition or sequence similarity. The sequence composition methods have the ability to completely assign the whole dataset. However, their use in metagenomics and the study of their performance with real data is limited. In this work, we assess the consistency of three different methods (BLAST + Lowest Common Ancestor, Phymm, and Naïve Bayesian Classifier) in assigning real and simulated sequence reads., Results: Both in real and in simulated data, BLAST + Lowest Common Ancestor (BLAST + LCA), Phymm, and Naïve Bayesian Classifier consistently assign a larger number of reads in higher taxonomic levels than in lower levels. However, discrepancies increase at lower taxonomic levels. In simulated data, consistent assignments between all three methods showed greater precision than assignments based on Phymm or Bayesian Classifier alone, since the BLAST + LCA algorithm performed best. In addition, assignment consistency in real data increased with sequence read length, in agreement with previously published simulation results., Conclusions: The use and combination of different approaches is advisable to assign metagenomic reads. Although the sensitivity could be reduced, the reliability can be increased by using the reads consistently assigned to the same taxa by, at least, two methods, and by training the programs using all available information.

In the last decade, an extensive effort has been made to characterize the human microbiota, due to its clinical and economic interests. However, a metagenomic approach to the skin microbiota is hampered by the high proportion of host DNA that is recovered. In contrast with the burgeoning field of gut metagenomics, skin metagenomics has been hindered by the absence of an efficient method to avoid sequencing the host DNA. We present here a method for recovering microbial DNA from skin samples, based on a combination of molecular techniques. We have applied this method to mouse skin, and have validated it by standard, quantitative PCR and amplicon sequencing of 16S rRNA. The taxonomic diversity recovered was not altered by this new method, as proved by comparing the phylogenetic structure revealed by 16S rRNA sequencing in untreated vs. treated samples. As proof of concept, we also present the first two mouse skin metagenomes, which allowed discovering new taxa (not only prokaryotes but also viruses and eukaryots) not reachable by 16S rRNA sequencing, as well as to characterize the skin microbiome functional landscape. Our method paves the way for the development of skin metagenomics, which will allow a much deeper knowledge of the skin microbiome and its relationship with the host, both in a healthy state and in relation to disease.

A mummified head was identified in 2010 as belonging to Henri IV, King of France. A putative blood sample from the King Louis XVI preserved into a pyrographically decorated gourd was analyzed in 2011. Both kings are in a direct male-line descent, separated by seven generations. We have retrieved the hypervariable region 1 of the mitochondrial DNA as well as a partial Y-chromosome profile from Henri IV. Five STR loci match the alleles found in Louis XVI, while another locus shows an allele that is just one mutation step apart. Taking into consideration that the partial Y-chromosome profile is extremely rare in modern human databases, we concluded that both males could be paternally related. The likelihood ratio of the two samples belonging to males separated by seven generations (as opposed to unrelated males) was estimated as 246.3, with a 95% confidence interval between 44.2 and 9729. Historically speaking, this forensic DNA data would confirm the identity of the previous Louis XVI sample, and give another positive argument for the authenticity of the head of Henri IV., (Copyright © 2012 Elsevier Ireland Ltd. All rights reserved.)

Host-commensal relationships in the skin are a complex system governed by variables related to the host, the bacteria and the environment. A disruption of this system may lead to new steady states, which, in turn, may lead to disease. We have studied one such disruption by characterizing the skin microbiota in healthy and immunodepressed (ID) mice. A detailed anatomopathological study failed to reveal any difference between the skin of healthy and ID mice. We sequenced the 16S rDNA V1-V2 gene region to saturation in 10 healthy and 10 ID 8 week-old mice, and found than all of the healthy and two of the ID mice had bacterial communities that were similar in composition to that of human skin, although, presumably because of the uniform raising conditions, less interindividual variation was found in mice. However, eight ID mice showed microbiota dominated by Staphylococcus epidermidis. Quantitative PCR amplification of 16S rDNA gene and of the Staphylococcus-specific TstaG region confirmed the previous results and indicated that the quantitative levels of Staphylococcus were similar in both groups while the total number of 16S copies was greater in the healthy mice. Thus, it is possible that, under long-term immunodeficiency, which removes the acquired but not the native immune system, S.epidermidis may inhibit the growth of other bacteria but does not cause a pathogenic state., (© 2012 Society for Applied Microbiology and Blackwell Publishing Ltd.)

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., (Copyright © 2012 Wiley Periodicals, Inc.)

Background: IgA nephropathy (IgAN) is not generally considered a hereditary disease, even though extensive evidence suggests an undefined genetic influence. Linkage analysis identified a number of genome regions that could contain variations linked to IgAN., Methods: In this case-control association study, genes possibly involved in the development of IgAN were investigated. DNA samples from 460 North Italian patients with IgAN and 444 controls were collected. Candidate genes were selected based on their possible functional involvement (6 genes) or because of their location within linkage-identified genomic regions IGAN2 and IGAN3 (5 and 13 genes within chromosome 4q26-31 and 17q12-22, respectively). One hundred and ninety-two tag and functional single-nucleotide polymorphisms (SNPs) were typed with Veracode GoldenGate technology (Illumina)., Results: C1GALT1 showed an association with IgAN (rs1008898: P = 0.0019 and rs7790522: P = 0.0049). Associations were found when the population was stratified by gender (C1GALT1, CD300LG, GRN, ITGA2B, ITGB3 in males and C1GALT1, TRPC3, B4GALNT2 in females) and by age (TLR4, ITGB3 in patients aged <27 years). Furthermore, rs7873784 in TLR4 showed an association with proteinuria (G allele: P = 0.0091; GG genotype: P = 0.0077)., Conclusions: Age and gender are likely to evidence distinct immunological and inflammatory reactions leading to individual susceptibility to IgAN. Overall, a genetic predisposition to sporadic IgAN was found. We might hypothesize that C1GALT1 and TLR4 polymorphisms influence the risk to develop IgAN and proteinuria, respectively.

We have analyzed human genetic diversity in 33 Old World populations including 23 populations obtained through Genographic Project studies. A set of 1,536 SNPs in five X chromosome regions were genotyped in 1,288 individuals (mostly males). We use a novel analysis employing subARG network construction with recombining chromosomal segments. Here, a subARG is constructed independently for each of five gene-free regions across the X chromosome, and the results are aggregated across them. For PCA, MDS and ancestry inference with STRUCTURE, the subARG is processed to obtain feature vectors of samples and pairwise distances between samples. The observed population structure, estimated from the five short X chromosomal segments, supports genome-wide frequency-based analyses: African populations show higher genetic diversity, and the general trend of shared variation is seen across the globe from Africa through Middle East, Europe, Central Asia, Southeast Asia, and East Asia in broad patterns. The recombinational analysis was also compared with established methods based on SNPs and haplotypes. For haplotypes, we also employed a fixed-length approach based on information-content optimization. Our recombinational analysis suggested a southern migration route out of Africa, and it also supports a single, rapid human expansion from Africa to East Asia through South Asia.

According to most historians, Christopher Columbus was born in Genoa, Italy. However, based on some key facts in the discoverer's biography, as well as in the linguistic analysis of his texts, some historians and linguists believe that Columbus could have been of Catalan origin. A Ligurian Columbus would have carried the Colombo surname, whereas he would have been called Colom if he were Catalan. In order to test whether it would be possible to discriminate between a Ligurian or a Catalan origin were Columbus' Y-chromosome haplotype to be retrieved, we genotyped 17 Y-chromosome STRs in 238 Spanish (from Catalonia, Valencia, and the Balearic Islands) and French Colom men, and 114 North Italian Colombo (from Liguria, Lombardy, and Piedmont). The Italian samples and, in particular, the Lombard Colombos were genetically as diverse as the general population, and we found little evidence of clusters of haplotypes that could indicate descent from a single founder. Colombo is actually the most frequent surname in Lombardy, where foundlings and orphans used to be given the surname Colombo. By contrast, Y-chromosome diversity was reduced in the Iberian Colom, where most of the men had Y chromosomes belonging to a few lineages. This implies that a positive identification would be more likely if Columbus were of Catalan descent. In this study, we have shown the diverse dynamics of two surnames linked by their etymology, in what is, to the best of our knowledge, the first genetic analysis of a surname in Southern Europe.

Vlad III The Impaler, also known as Dracula, descended from the dynasty of Basarab, the first rulers of independent Wallachia, in present Romania. Whether this dynasty is of Cuman (an admixed Turkic people that reached Wallachia from the East in the 11(th) century) or of local Romanian (Vlach) origin is debated among historians. Earlier studies have demonstrated the value of investigating the Y chromosome of men bearing a historical name, in order to identify their genetic origin. We sampled 29 Romanian men carrying the surname Basarab, in addition to four Romanian populations (from counties Dolj, N = 38; Mehedinti, N = 11; Cluj, N = 50; and Brasov, N = 50), and compared the data with the surrounding populations. We typed 131 SNPs and 19 STRs in the non-recombinant part of the Y-chromosome in all the individuals. We computed a PCA to situate the Basarab individuals in the context of Romania and its neighboring populations. Different Y-chromosome haplogroups were found within the individuals bearing the Basarab name. All haplogroups are common in Romania and other Central and Eastern European populations. In a PCA, the Basarab group clusters within other Romanian populations. We found several clusters of Basarab individuals having a common ancestor within the period of the last 600 years. The diversity of haplogroups found shows that not all individuals carrying the surname Basarab can be direct biological descendants of the Basarab dynasty. The absence of Eastern Asian lineages in the Basarab men can be interpreted as a lack of evidence for a Cuman origin of the Basarab dynasty, although it cannot be positively ruled out. It can be therefore concluded that the Basarab dynasty was successful in spreading its name beyond the spread of its genes.

A text on a pyrographically decorated gourd dated to 1793 explains that it contains a handkerchief dipped with the blood of Louis XVI, king of France, after his execution. Biochemical analyses confirmed that the material contained within the gourd was blood. The mitochondrial DNA (mtDNA) hypervariable region 1 (HVR1) and 2 (HVR2), the Y-chromosome STR profile, some autosomal STR markers and a SNP in HERC2 gene associated to blue eyes, were retrieved, and some results independently replicated in two different laboratories. The uncommon mtDNA sequence retrieved can be attributed to a N1b haplotype, while the novel Y-chromosome haplotype belongs to haplogroup G2a. The HERC2 gene showed that the subject analyzed was a heterozygote, which is compatible with a blue-eyed person, as king Louis XVI was. To confirm the identity of the subject, an analysis of the dried heart of his son, Louis XVII, could be undertaken., (Copyright © 2010 Elsevier Ireland Ltd. All rights reserved.)

Summary: Given a set of extant haplotypes IRiS first detects high confidence recombination events in their shared genealogy. Next using the local sequence topology defined by each detected event, it integrates these recombinations into an ancestral recombination graph. While the current system has been calibrated for human population data, it is easily extendible to other species as well., Availability: IRiS (Identification of Recombinations in Sequences) binary files are available for non-commercial use in both Linux and Microsoft Windows, 32 and 64 bit environments from https://researcher.ibm.com/researcher/view&#95;project.php?id = 2303, Contact: parida@us.ibm.com.

Background: Pathogens have represented an important selective force during the adaptation of modern human populations to changing social and other environmental conditions. The evolution of the immune system has therefore been influenced by these pressures. Genomic scans have revealed that immune system is one of the functions enriched with genes under adaptive selection., Results: Here, we describe how the innate immune system has responded to these challenges, through the analysis of resequencing data for 132 innate immunity genes in two human populations. Results are interpreted in the context of the functional and interaction networks defined by these genes. Nucleotide diversity is lower in the adaptors and modulators functional classes, and is negatively correlated with the centrality of the proteins within the interaction network. We also produced a list of candidate genes under positive or balancing selection in each population detected by neutrality tests and showed that some functional classes are preferential targets for selection., Conclusions: We found evidence that the role of each gene in the network conditions the capacity to evolve or their evolvability: genes at the core of the network are more constrained, while adaptation mostly occurred at particular positions at the network edges. Interestingly, the functional classes containing most of the genes with signatures of balancing selection are involved in autoinflammatory and autoimmune diseases, suggesting a counterbalance between the beneficial and deleterious effects of the immune response.

Human population movements in North Africa have been mostly restricted to an east-west direction due to the geographical barriers imposed by the Sahara Desert and the Mediterranean Sea. Although these barriers have not completely impeded human migrations, genetic studies have shown that an east-west genetic gradient exists. However, the lack of genetic information of certain geographical areas and the focus of some studies in parts of the North African landscape have limited the global view of the genetic pool of North African populations. To provide a global view of the North African genetic landscape and population structure, we have analyzed ∼2,300 North African mitochondrial DNA lineages (including 269 new sequences from Libya, in the first mtDNA study of the general Libyan population). Our results show a clinal distribution of certain haplogroups, some of them more frequent in Western (H, HV0, L1b, L3b, U6) or Eastern populations (L0a, R0a, N1b, I, J) that might be the result of human migrations from the Middle East, sub-Saharan Africa, and Europe. Despite this clinal pattern, a genetic discontinuity is found in the Libyan/Egyptian border, suggesting a differential gene flow in the Nile River Valley. Finally, frequency of the post-LGM subclades H1 and H3 is predominant in Libya within the H sequences, highlighting the magnitude of the LGM expansion in North Africa., (Copyright © 2011 Wiley-Liss, Inc.)

Recombination varies greatly among species, as illustrated by the poor conservation of the recombination landscape between humans and chimpanzees. Thus, shorter evolutionary time frames are needed to understand the evolution of recombination. Here, we analyze its recent evolution in humans. We calculated the recombination rates between adjacent pairs of 636,933 common single-nucleotide polymorphism loci in 28 worldwide human populations and analyzed them in relation to genetic distances between populations. We found a strong and highly significant correlation between similarity in the recombination rates corrected for effective population size and genetic differentiation between populations. This correlation is observed at the genome-wide level, but also for each chromosome and when genetic distances and recombination similarities are calculated independently from different parts of the genome. Moreover, and more relevant, this relationship is robustly maintained when considering presence/absence of recombination hotspots. Simulations show that this correlation cannot be explained by biases in the inference of recombination rates caused by haplotype sharing among similar populations. This result indicates a rapid pace of evolution of recombination, within the time span of differentiation of modern humans.

Pygmy populations are among the few hunter-gatherers currently living in sub-Saharan Africa and are mainly represented by two groups, Eastern and Western, according to their current geographical distribution. They are scattered across the Central African belt and surrounded by Bantu-speaking farmers, with whom they have complex social and economic interactions. To investigate the demographic history of Pygmy groups, a population approach was applied to the analysis of 205 complete mitochondrial DNA (mtDNA) sequences from ten central African populations. No sharing of maternal lineages was observed between the two Pygmy groups, with haplogroup L1c being characteristic of the Western group but most of Eastern Pygmy lineages falling into subclades of L0a, L2a, and L5. Demographic inferences based on Bayesian coalescent simulations point to an early split among the maternal ancestors of Pygmies and those of Bantu-speaking farmers (∼ 70,000 years ago [ya]). Evidence for population growth in the ancestors of Bantu-speaking farmers has been observed, starting ∼ 65,000 ya, well before the diffusion of Bantu languages. Subsequently, the effective population size of the ancestors of Pygmies remained constant over time and ∼ 27,000 ya, coincident with the Last Glacial Maximum, Eastern and Western Pygmies diverged, with evidence of subsequent migration only among the Western group and the Bantu-speaking farmers. Western Pygmies show signs of a recent bottleneck 4,000-650 ya, coincident with the diffusion of Bantu languages, whereas Eastern Pygmies seem to have experienced a more ancient decrease in population size (20,000-4,000 ya). In conclusion, the results of this first attempt at analyzing complete mtDNA sequences at the population level in sub-Saharan Africa not only support previous findings but also offer new insights into the demographic history of Pygmy populations, shedding new light on the ancient peopling of the African continent.

Background: Searching for associations between genetic variants and complex diseases has been a very active area of research for over two decades. More than 51,000 potential associations have been studied and published, a figure that keeps increasing, especially with the recent explosion of array-based Genome-Wide Association Studies. Even if the number of true associations described so far is high, many of the putative risk variants detected so far have failed to be consistently replicated and are widely considered false positives. Here, we focus on the world-wide patterns of replicability of published association studies., Results: We report three main findings. First, contrary to previous results, genes associated to complex diseases present lower degrees of genetic differentiation among human populations than average genome-wide levels. Second, also contrary to previous results, the differences in replicability of disease associated-loci between Europeans and East Asians are highly correlated with genetic differentiation between these populations. Finally, highly replicated genes present increased levels of high-frequency derived alleles in European and Asian populations when compared to African populations., Conclusions: Our findings highlight the heterogeneous nature of the genetic etiology of complex disease, confirm the importance of the recent evolutionary history of our species in current patterns of disease susceptibility and could cast doubts on the status as false positives of some associations that have failed to replicate across populations.

Previous genetic, anthropological and linguistic studies have shown that Roma (Gypsies) constitute a founder population dispersed throughout Europe whose origins might be traced to the Indian subcontinent. Linguistic and anthropological evidence point to Indo-Aryan ethnic groups from North-western India as the ancestral parental population of Roma. Recently, a strong genetic hint supporting this theory came from a study of a private mutation causing primary congenital glaucoma. In the present study, complete mitochondrial control sequences of Iberian Roma and previously published maternal lineages of other European Roma were analyzed in order to establish the genetic affinities among Roma groups, determine the degree of admixture with neighbouring populations, infer the migration routes followed since the first arrival to Europe, and survey the origin of Roma within the Indian subcontinent. Our results show that the maternal lineage composition in the Roma groups follows a pattern of different migration routes, with several founder effects, and low effective population sizes along their dispersal. Our data allowed the confirmation of a North/West migration route shared by Polish, Lithuanian and Iberian Roma. Additionally, eleven Roma founder lineages were identified and degrees of admixture with host populations were estimated. Finally, the comparison with an extensive database of Indian sequences allowed us to identify the Punjab state, in North-western India, as the putative ancestral homeland of the European Roma, in agreement with previous linguistic and anthropological studies.

Background: Despite the successful retrieval of genomes from past remains, the prospects for human palaeogenomics remain unclear because of the difficulty of distinguishing contaminant from endogenous DNA sequences. Previous sequence data generated on high-throughput sequencing platforms indicate that fragmentation of ancient DNA sequences is a characteristic trait primarily arising due to depurination processes that create abasic sites leading to DNA breaks. METHODOLOGY/PRINCIPALS FINDINGS: To investigate whether this pattern is present in ancient remains from a temperate environment, we have 454-FLX pyrosequenced different samples dated between 5,500 and 49,000 years ago: a bone from an extinct goat (Myotragus balearicus) that was treated with a depurinating agent (bleach), an Iberian lynx bone not subjected to any treatment, a human Neolithic sample from Barcelona (Spain), and a Neandertal sample from the El Sidrón site (Asturias, Spain). The efficiency of retrieval of endogenous sequences is below 1% in all cases. We have used the non-human samples to identify human sequences (0.35 and 1.4%, respectively), that we positively know are contaminants., Conclusions: We observed that bleach treatment appears to create a depurination-associated fragmentation pattern in resulting contaminant sequences that is indistinguishable from previously described endogenous sequences. Furthermore, the nucleotide composition pattern observed in 5' and 3' ends of contaminant sequences is much more complex than the flat pattern previously described in some Neandertal contaminants. Although much research on samples with known contaminant histories is needed, our results suggest that endogenous and contaminant sequences cannot be distinguished by the fragmentation pattern alone.

Studies of large sets of single nucleotide polymorphism (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 2841 SNPs in 12 sub-Saharan African populations, including a previously unsampled region of southeastern Africa (Mozambique). We show that robust results in a world-wide perspective can be obtained when analyzing only 1000 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 southeastern Bantu population from Mozambique, which suggests an assimilation of a pre-Bantu substrate by Bantu speakers in the region.

Recombination is one of the main forces shaping genome diversity, but the information it generates is often overlooked. A recombination event creates a junction between two parental sequences that may be transmitted to the subsequent generations. Just like mutations, these junctions carry evidence of the shared past of the sequences. We present the IRiS algorithm, which detects past recombination events from extant sequences and specifies the place of each recombination and which are the recombinants sequences. We have validated and calibrated IRiS for the human genome using coalescent simulations replicating standard human demographic history and a variable recombination rate model, and we have fine-tuned IRiS parameters to simultaneously optimize for false discovery rate, sensitivity, and accuracy in placing the recombination events in the sequence. Newer recombinations overwrite traces of past ones and our results indicate more recent recombinations are detected by IRiS with greater sensitivity. IRiS analysis of the MS32 region, previously studied using sperm typing, showed good concordance with estimated recombination rates. We also applied IRiS to haplotypes for 18 X-chromosome regions in HapMap Phase 3 populations. Recombination events detected for each individual were recoded as binary allelic states and combined into recotypes. Principal component analysis and multidimensional scaling based on recotypes reproduced the relationships between the eleven HapMap Phase III populations that can be expected from known human population history, thus further validating IRiS. We believe that our new method will contribute to the study of the distribution of recombination events across the genomes and, for the first time, it will allow the use of recombination as genetic marker to study human genetic variation.

Basques are a cultural isolate, and, according to mainly allele frequencies of classical polymorphisms, also a genetic isolate. We investigated the differentiation of Spanish Basques from the rest of Iberian populations by means of a dense, genome-wide SNP array. We found that F (ST) distances between Spanish Basques and other populations were similar to those between pairs of non-Basque populations. The same result is found in a PCA of individuals, showing a general distinction between Iberians and other South Europeans independently of being Basques. Pathogen-mediated natural selection may be responsible for the high differentiation previously reported for Basques at very specific genes such as ABO, RH, and HLA. Thus, Basques cannot be considered a genetic outlier under a general genome scope and interpretations on their origin may have to be revised.

The level of Factor XII (FXII) is an important phenotype that exhibits a high genetic component and is associated with thrombotic disease. In a genome-wide linkage scan, we demonstrated that the F12 gene represents a quantitative trait locus (QTL) that influences FXII levels. The current study investigated the genetic architecture of the F12 gene to locate polymorphism(s) responsible for the variation of FXII levels. Re-sequencing of the F12 gene in 40 unrelated individuals (selected from the tails of normal distribution of FXII levels) identified 26 polymorphisms which were genotyped in 398 individuals belonging to 21 families from the GAIT Project. By a measured genotype association analysis, eight of 26 SNPs showed significant P-values less than 10(-5) (after multiple test correction) with FXII levels. In addition, the Bayesian Quantitative Trait Nucleotide method, which infers those polymorphisms most likely to have a direct influence on the trait under study, provided evidence that only rs1801020 variation accounted for the variance attributed to this QTL. Moreover, we have analyzed the evolutionary processes that produced the variation in F12 gene and concluded that is evolutionarily neutral and that the T allele of the rs1801020 appeared approximately 100 000 years ago and spread to most human populations rising to high frequencies by genetic drift. Our study provides a template for future genetic studies of human quantitative traits, as we move beyond QTL localization to the polymorphisms responsible for the variation of important biomedical phenotypes.

The Basques are a culturally isolated population, living across the western border between France and Spain and speaking a non-Indo-European language. They show outlier allele frequencies in the ABO, RH, and HLA loci. To test whether Basques are a genetic isolate with the features that would make them good candidates in genetic association studies, we genotyped 123 SNPs in a 1-Mb region in chromosome 22 in Basque samples from France and Spain, as well as in samples from northern and southern Spain, and in three North African samples. Both Basque samples showed similar levels of heterozygosity to the other populations, and the decay of linkage disequilibrium with physical distance was not different between Basques and non-Basques. Thus, Basques do not show the genetic properties expected in population isolates.

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 the American continent. The Mesoamerican isthmus has constituted an important geographic barrier that has severely restricted gene flow between North and South America in pre-historical times. Although the Native American component has been already described in admixed Mexican populations, few studies have been carried out in native Mexican populations. In this study, we present mitochondrial DNA (mtDNA) sequence data for the first hypervariable region (HVR-I) in 477 unrelated individuals belonging to 11 different native populations from Mexico. Almost all of the Native Mexican mtDNAs could be classified into the four pan-Amerindian haplogroups (A2, B2, C1, and D1); only two 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, Native Mexican populations show a high degree of heterogeneity in their patterns of haplogroup frequencies. Haplogroup X2a 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 in 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.

Coding variants in tumor necrosis factor receptor superfamily member 13B (TNFRSF13B) have been implicated in common variable immunodeficiency (CVID), but the functional effects of such mutations in relation to the development of the disease have not been entirely established. To examine the potential contribution of TNFRSF13B variants to CVID, we have applied an evolutionary approach by sequencing its coding region in 451 individuals belonging to 26 worldwide populations, in addition to controls, patients with CVID and selective IgA deficiency (IgAD) from Italy. The low level of geographical structure for the observed genetic diversity and the several neutrality tests performed confirm the absence of recent population-specific selective pressures, suggesting that TNFRSF13B may be involved also in innate immune functions, rather than in adaptive immunity only. A slight excess of rare derived alleles was found in patients with CVID, and thus some of these variants may contribute to the disease, implying that CVID probably fits the rare variants rather than the common disease/common variant paradigm. This also confirms the previous suggestion that TNFRSF13B defects alone do not cause CVID and that such an extremely heterogeneous immunodeficiency might be more likely related to additional, still unknown environmental and genetic factors.

Background: It is well known that the pattern of linkage disequilibrium varies between human populations, with remarkable geographical stratification. Indirect association studies routinely exploit linkage disequilibrium around genes, particularly in isolated populations where it is assumed to be higher. Here, we explore both the amount and the decay of linkage disequilibrium with physical distance along 211 gene regions, most of them related to complex diseases, across 39 HGDP-CEPH population samples, focusing particularly on the populations defined as isolates. Within each gene region and population we use r2 between all possible single nucleotide polymorphism (SNP) pairs as a measure of linkage disequilibrium and focus on the proportion of SNP pairs with r2 greater than 0.8., Results: Although the average r2 was found to be significantly different both between and within continental regions, a much higher proportion of r2 variance could be attributed to differences between continental regions (2.8% vs. 0.5%, respectively). Similarly, while the proportion of SNP pairs with r2 > 0.8 was significantly different across continents for all distance classes, it was generally much more homogenous within continents, except in the case of Africa and the Americas. The only isolated populations with consistently higher LD in all distance classes with respect to their continent are the Kalash (Central South Asia) and the Surui (America). Moreover, isolated populations showed only slightly higher proportions of SNP pairs with r2 > 0.8 per gene region than non-isolated populations in the same continent. Thus, the number of SNPs in isolated populations that need to be genotyped may be only slightly less than in non-isolates., Conclusion: The "isolated population" label by itself does not guarantee a greater genotyping efficiency in association studies, and properties other than increased linkage disequilibrium may make these populations interesting in genetic epidemiology.

The GT6 glycosyltransferases gene family, that includes the ABO blood group, shows a complex evolution pattern, with multiple events of gain and loss in different mammal species. In humans the ABO gene is considered the sole functional member although the O allele is null and is fixed in certain populations. Here, we analyze the human GT6 pseudogene sequences (Forssman, IGB3, GGTA1, GT6m5, GT6m6, and GT6m7) from an evolutionary perspective, by the study of (i) their diversity levels in populations through the resequencing analysis of European and African individuals; (ii) the interpopulation differentiation, with genotyping data from a survey of populations covering most of human genetic diversity; and (iii) the interespecific divergence, by the comparison of the human and some other primate species sequences. Since pseudogenes are expected to evolve under neutrality, they should show an evolutionary pattern different to that of functional sequences, with higher levels of diversity as well as a ratio of nonsynonymous to synonymous changes close to 1. We describe some departures from these expectations, including selection for inactivation in IGB3, GGTA1, and the interesting case of FS (Forssman) with a probable shift of its initial function in the primate lineage, which put it apart from a pure neutral pseudogene. These results suggest that some of these GT6 human pseudogenes may still be functional and retain some valuable unknown function in humans, in some case even at the protein level. The evolutionary analysis of all members of the GT6 family in humans allows an insight into their functional history, a process likely due to the interaction of the host glycans that they synthesize with pathogens; the past process that can be unraveled through the footprints left by natural selection in the extant genome variation.

Background: We address the problem of studying recombinational variations in (human) populations. In this paper, our focus is on one computational aspect of the general task: Given two networks G1 and G2, with both mutation and recombination events, defined on overlapping sets of extant units the objective is to compute a consensus network G3 with minimum number of additional recombinations. We describe a polynomial time algorithm with a guarantee that the number of computed new recombination events is within = sz(G1, G2) (function sz is a well-behaved function of the sizes and topologies of G1 and G2) of the optimal number of recombinations. To date, this is the best known result for a network consensus problem., Results: Although the network consensus problem can be applied to a variety of domains, here we focus on structure of human populations. With our preliminary analysis on a segment of the human Chromosome X data we are able to infer ancient recombinations, population-specific recombinations and more, which also support the widely accepted 'Out of Africa' model. These results have been verified independently using traditional manual procedures. To the best of our knowledge, this is the first recombinations-based characterization of human populations., Conclusion: We show that our mathematical model identifies recombination spots in the individual haplotypes; the aggregate of these spots over a set of haplotypes defines a recombinational landscape that has enough signal to detect continental as well as population divide based on a short segment of Chromosome X. In particular, we are able to infer ancient recombinations, population-specific recombinations and more, which also support the widely accepted 'Out of Africa' model. The agreement with mutation-based analysis can be viewed as an indirect validation of our results and the model. Since the model in principle gives us more information embedded in the networks, in our future work, we plan to investigate more non-traditional questions via these structures computed by our methodology.

Most studies of European genetic diversity have focused on large-scale variation and interpretations based on events in prehistory, but migrations and invasions in historical times could also have had profound effects on the genetic landscape. The Iberian Peninsula provides a suitable region for examination of the demographic impact of such recent events, because its complex recent history has involved the long-term residence of two very different populations with distinct geographical origins and their own particular cultural and religious characteristics-North African Muslims and Sephardic Jews. To address this issue, we analyzed Y chromosome haplotypes, which provide the necessary phylogeographic resolution, in 1140 males from the Iberian Peninsula and Balearic Islands. Admixture analysis based on binary and Y-STR haplotypes indicates a high mean proportion of ancestry from North African (10.6%) and Sephardic Jewish (19.8%) sources. Despite alternative possible sources for lineages ascribed a Sephardic Jewish origin, these proportions attest to a high level of religious conversion (whether voluntary or enforced), driven by historical episodes of social and religious intolerance, that ultimately led to the integration of descendants. In agreement with the historical record, analysis of haplotype sharing and diversity within specific haplogroups suggests that the Sephardic Jewish component is the more ancient. The geographical distribution of North African ancestry in the peninsula does not reflect the initial colonization and subsequent withdrawal and is likely to result from later enforced population movement-more marked in some regions than in others-plus the effects of genetic drift.

Unlabelled: We have developed a software that applies ascertainment bias on simulated DNA sequences and calculates F(ST) on them, so they can be used to generate neutral distributions that are appropriate to test whether the genetic differentiation of a particular gene between populations is compatible with neutral evolution, or, on the contrary, suggests local adaptation by natural selection., Availability: FABSIM is available from http://www.snpator.com/public/downloads/aRamirez/FABSIM/.

Traditionally nonrecombinant genome, i.e., mtDNA or Y chromosome, has been used for phylogeography, notably for ease of analysis. The topology of the phylogeny structure in this case is an acyclic graph, which is often a tree, is easy to comprehend and is somewhat easy to infer. However, recombination is an undeniable genetic fact for most part of the genome. Driven by the need for a more complete analysis, we address the problem of estimating the ancestral recombination graph (ARG) from a collection of extant sequences. We exploit the coherence that is observed in the human haplotypes as patterns and present a network model of patterns to reconstruct the ARG. We test our model on simulations that closely mimic the observed haplotypes and observe promising results.

The ABO polymorphism has long been suspected to be under balancing selection. To explore this possibility, we analyzed two datasets: (1) a set of 94 23-Kb sequences in European- and African-Americans produced by the Seattle SNPs project, and (2) a set of 814 2-Kb sequences in O alleles from seven worldwide populations. A phylogenetic analysis of the Seattle sequences showed a complex pattern in which the action of recombination and gene conversion are evident, and in which four main lineages could be individuated. The sequence patterns could be linked to the expected blood group phenotype; in particular, the main mutation giving rise to the null O allele is likely to have appeared at least three times in human evolution, giving rise to allele lineages O02, O01, and O09. However, the genealogy changes along the gene and variations of both numbers of branches and of their time depth were observed, which could result from a combined action of recombination and selection. Several neutrality tests clearly demonstrated deviations compatible with balancing selection, peaking at several locations along the gene. The time depth of the genealogy was also incompatible with neutral evolution, particularly in the region from exons 6 to 7, which codes for most of the catalytic domain.

Background: Before the arrival of Europeans to Cuba, the island was inhabited by two Native American groups, the Tainos and the Ciboneys. Most of the present archaeological, linguistic and ancient DNA evidence indicates a South American origin for these populations. In colonial times, Cuban Native American people were replaced by European settlers and slaves from Africa. It is still unknown however, to what extent their genetic pool intermingled with and was 'diluted' by the arrival of newcomers. In order to investigate the demographic processes that gave rise to the current Cuban population, we analyzed the hypervariable region I (HVS-I) and five single nucleotide polymorphisms (SNPs) in the mitochondrial DNA (mtDNA) coding region in 245 individuals, and 40 Y-chromosome SNPs in 132 male individuals., Results: The Native American contribution to present-day Cubans accounted for 33% of the maternal lineages, whereas Africa and Eurasia contributed 45% and 22% of the lineages, respectively. This Native American substrate in Cuba cannot be traced back to a single origin within the American continent, as previously suggested by ancient DNA analyses. Strikingly, no Native American lineages were found for the Y-chromosome, for which the Eurasian and African contributions were around 80% and 20%, respectively., Conclusion: While the ancestral Native American substrate is still appreciable in the maternal lineages, the extensive process of population admixture in Cuba has left no trace of the paternal Native American lineages, mirroring the strong sexual bias in the admixture processes taking place during colonial times.

The evolutionarily recent geographic expansion of humans, and the even more recent development of large, relatively dense human settlements, has exposed our species to new pathogenic environments. Potentially lethal pathogens are likely to have exerted important selective pressures on our genome, so immunity genes can be expected to show molecular signatures of the adaptation of human populations to these recent conditions. While genes related to the acquired immunity system have indeed been reported to show traces of local adaptation, little is known about the response of the innate immunity system. In this study, we analyze the variability patterns in different human populations of fifteen genes related to innate immunity. We have used both single nucleotide polymorphism and sequence data, and through the analysis of interpopulation differentiation, the linkage disequilibrium pattern, and intrapopulation diversity, we have discovered some signatures of positive and especially balancing selection in these genes, thus confirming the importance of the immune system genetic plasticity in the evolutionary adaptive process. Interestingly, the strongest evidence is found in three TLR genes and CD14. These innate immunity genes play a pivotal role, being involved in the primary recognition of pathogens. In general, more evidences of selection appear in the European populations, in some case possibly related to severe population specific pressures. However, we also describe evidence from African populations, which may reflect parallel or long-term selective forces acting in different geographic areas.

Unlabelled: Single nucleotide polymorphisms (SNPs) are the most widely used marker in studies to assess associations between genetic variants and complex traits or diseases. They are also becoming increasingly important in the study of the evolution and history of humans and other species. The analysis and processing of SNPs obtained thanks to high-throughput technologies imply the time consuming and costly use of different, complex and usually format-incompatible software. SNPator is a user-friendly web-based SNP data analysis suite that integrates, among many other algorithms, the most common steps of a SNP association study. It frees the user from the need to have large computer facilities and an in depth knowledge of genetic software installation and management. Genotype data is directly read from the output files of the usual genotyping platforms. Phenotypic data on the samples can also be easily uploaded. Many different quality control and analysis procedures can be performed either by using built-in SNPator algorithms or by calling standard genetic software., Availability: Access is granted from the SNPator webpage http://www.snpator.org.

In repeat expansion disorders, the size of pathological alleles is the most relevant factor accounting for the disease severity and age-at-onset, emphasizing the clinical significance of their underlying intergenerational instability. In one of these diseases, Machado-Joseph disease (MJD), the sex of transmitting progenitor and the C(987)GG/G(987)GG polymorphism are the best studied factors acting on intergenerational instability of expanded alleles. Here, we assessed the influence of other cis and inter-allelic acting factors, at the ATXN3 locus, through the analysis of MJD lineages, flanking STR-based haplotypes, the initial repeat size and parental age. A total of 100 transmissions of the expanded MJD allele were analyzed according to the sex of the transmitting parent. We have shown that independent origin mutations (identified by intragenic SNP-based haplotypes) behave differently, as the status of instability (contraction, no change or further expansion) is concerned. Indeed, 72% of expansions were associated to the worldwide spread TTACAC lineage, whereas the GTGGCA displayed 75% of all contractions observed. The analysis of flanking recombinant haplotypes did not suggest any further distant cis elements acting up- or downstream the ATXN3 locus. Considering the increased amplitude of expansions seen in older transmitting fathers, a repair-based mechanism may be suggested for the meiotic instability at this locus; furthermore, the lack of correlation between the initial repeat size and degree of instability did not support a replication-based mechanism. In summary, our findings point to different mechanisms of instability underlying male and female meioses, as well as contraction and expansion processes in MJD., ((c) 2007 Wiley-Liss, Inc.)

La Réunion, one of the three Mascarene islands located in the Indian Ocean, remained devoid of inhabitants until it was first colonized by the French in the middle of the 17th century. The continuous flow of foreign-born slaves and immigrant workers from Africa, India, Southeast Asia, and China to work on coffee and sugar cane plantations led to the island becoming a melting pot of people of multiple ethnic origins. To establish the impact of the different incoming ethnic groups on the present Reunionese gene pool, we have sequenced both hypervariable regions I and II of the mitochondrial DNA molecule, the 9 bp COII/tRNA(Lys) deletion, and four SNPs located in the coding region in a total of 41 samples of the general population, and a further 18 STRs and 35 SNPs on the Y chromosome in 26 of these samples. Our results show that there was a strong sexual bias (asymmetrical gene flow) in the peopling of La Réunion, where admixture events were mainly between male settlers and females from the incoming slave groups. Most of the Y-chromosome gene pool is of European/Middle Eastern ancestry (85%), whereas the mtDNA gene pool is mainly of Indian and East Asian ancestry (70%). The absence of genetic diversity within these two major components of the mtDNA gene pool suggests these populations may have witnessed strong founder effects during the colonization process., ((c) 2008 Wiley-Liss, Inc.)

Several tests have been proposed to detect departures of nucleotide variability patterns from neutral expectations. However, very different kinds of evolutionary processes, such as selective events or demographic changes, can produce similar deviations from these tests, thus making interpretation difficult when a significant departure of neutrality is detected. Here we study the effects of demography and recombination upon neutrality tests by analyzing their power under sudden population expansions, sudden contractions, and bottlenecks. We evaluate tests based on the frequency spectrum of mutations and the distribution of haplotypes and explore the consequences of using incorrect estimates of the rates of recombination when testing for neutrality. We show that tests that rely on haplotype frequencies-especially Fs and ZnS, which are based, respectively, on the number of different haplotypes and on the r2 values between all pairs of polymorphic sites-are the most powerful for detecting expansions on nonrecombining genomic regions. Nevertheless, they are strongly affected by misestimations of recombination, so they should not be used when recombination levels are unknown. Instead, class I tests, particularly Tajima's D or R2, are recommended.

The European Gypsies, commonly referred to as Roma, are represented by a vast number of groups spread across many countries. Although sharing a common origin, the Gypsy groups are highly heterogeneous as a consequence of genetic drift and different levels of admixture with surrounding populations. With this study we aimed at contributing to the knowledge of the Roma history by studying 17 Y-STR and 34 Y-SNP loci in a sample of 126 Portuguese Gypsies. Distinct genetic hallmarks of their past and migration route were detected, namely: an ancestral component, shared by all Roma groups, that reflects their origin in India (H1a-M82; approximately 17%); an influence from their long permanence in the Balkans/Middle-East region (J2a1b-M67, J2a1b1-M92, I-M170, Q-M242; approximately 31%); traces of contacts with European populations preceding the entrance in the Iberian Peninsula (R1b1c-M269, J2b1a-M241; approximately 10%); and a high proportion of admixture with the non-Gypsy population from Iberia (R1b1c-M269, R1-M173/del.M269, J2a-M410, I1b1b-M26, E3b1b-M81; approximately 37%). Among the Portuguese Gypsies the proportion of introgression from host populations is higher than observed in other groups, a fact which is somewhat unexpected since the arrival of the Roma to Portugal is documented to be more recent than in Central or East Europe.

The human olfactory receptor (OR) repertoire is reduced in comparison to other mammals and to other nonhuman primates. Nonetheless, this olfactory decline opens an opportunity for evolutionary innovation and improvement. In the present study, we focus on an OR gene, OR5I1, which had previously been shown to present an excess of amino acid replacement substitutions between humans and chimpanzees. We analyze the genetic variation in OR5I1 in a large worldwide human panel and find an excess of derived alleles segregating at relatively high frequencies in all populations. Additional evidence for selection includes departures from neutrality in allele frequency spectra tests but no unusually extended haplotype structure. Moreover, molecular structural inference suggests that one of the nonsynonymous polymorphisms defining the presumably adaptive protein form of OR5I1 may alter the functional binding properties of the OR. These results are compatible with positive selection having modeled the pattern of variation found in the OR5I1 gene and with a relatively ancient, mild selective sweep predating the "Out of Africa" expansion of modern humans.