Search

The strategic location of North Africa has led to cultural and demographic shifts, shaping its genetic structure. Historical migrations brought different genetic components that are evident in present-day North African genomes, along with autochthonous components. The Imazighen (plural of Amazigh) are believed to be the descendants of autochthonous North Africans and speak various Amazigh languages, which belong to the Afro-Asiatic language family. However, the arrival of different human groups, especially during the Arab conquest, caused cultural and linguistic changes in local populations, increasing their heterogeneity. We aim to characterize the genetic structure of the region, using the largest Amazigh dataset to date and other reference samples. Our findings indicate microgeographical genetic heterogeneity among Amazigh populations, modeled by various admixture waves and different effective population sizes. A first admixture wave is detected group-wide around the twelfth century, whereas a second wave appears in some Amazigh groups around the nineteenth century. These events involved populations with higher genetic ancestry from south of the Sahara compared to the current North Africans. A plausible explanation would be the historical trans-Saharan slave trade, which lasted from the Roman times to the nineteenth century. Furthermore, our investigation shows that assortative mating in North Africa has been rare., (© 2024. The Author(s).)

The strategic location of North Africa has made the region the core of a wide range of human demographic events, including migrations, bottlenecks, and admixture processes. This has led to a complex and heterogeneous genetic and cultural landscape, which remains poorly studied compared to other world regions. Whole-exome sequencing is particularly relevant to determine the effects of these demographic events on current-day North Africans' genomes, since it allows to focus on those parts of the genome that are more likely to have direct biomedical consequences. Whole-exome sequencing can also be used to assess the effect of recent demography in functional genetic variation and the efficacy of natural selection, a long-lasting debate. In the present work, we use newly generated whole-exome sequencing and genome-wide array genotypes to investigate the effect of demography in functional variation in 7 North African populations, considering both cultural and demographic differences and with a special focus on Amazigh (plur. Imazighen) groups. We detect genetic differences among populations related to their degree of isolation and the presence of bottlenecks in their recent history. We find differences in the functional part of the genome that suggest a relaxation of purifying selection in the more isolated groups, allowing for an increase of putatively damaging variation. Our results also show a shift in mutational load coinciding with major demographic events in the region and reveal differences within and between cultural and geographic groups., Competing Interests: Conflict of interest statement. None declared., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

We provide the first whole genome Copy Number Variant (CNV) study addressing Roma, along with reference populations from South Asia, the Middle East and Europe. Using CNV calling software for short-read sequence data, we identified 3171 deletions and 489 duplications. Taking into account the known population history of the Roma, as inferred from whole genome nucleotide variation, we could discern how this history has shaped CNV variation. As expected, patterns of deletion variation, but not duplication, in the Roma followed those obtained from single nucleotide polymorphisms (SNPs). Reduced effective population size resulting in slightly relaxed natural selection may explain our observation of an increase in intronic (but not exonic) deletions within Loss of Function (LoF)-intolerant genes. Over-representation analysis for LoF-intolerant gene sets hosting intronic deletions highlights a substantial accumulation of shared biological processes in Roma, intriguingly related to signaling, nervous system and development features, which may be related to the known profile of private disease in the population. Finally, we show the link between deletions and known trait-related SNPs reported in the genome-wide association study (GWAS) catalog, which exhibited even frequency distributions among the studied populations. This suggests that, in general human populations, the strong association between deletions and SNPs associated to biomedical conditions and traits could be widespread across continental populations, reflecting a common background of potentially disease/trait-related CNVs., (© 2023. The Author(s).)

The use of mitochondrial DNA (mtDNA) in the field of forensic genetics is widely spread mainly due to its advantages when identifying highly degraded samples. In this sense, massive parallel sequencing has made the analysis of the whole mitogenome more accessible, noticeably increasing the informativeness of mtDNA haplotypes. The civil war (1980-1992) in El Salvador caused many deaths and disappearances (including children) all across the country and the economic and social instability after the war forced many people to emigration. For this reason, different organizations have collected DNA samples from relatives with the aim of identifying missing people. Thus, we present a dataset containing 334 complete mitogenomes from the Salvadoran general population. To the best of our knowledge, this is the first publication of a nationwide forensic-quality complete mitogenome database of any Latin American country. We found 293 different haplotypes, with a random match probability of 0.0041 and 26.6 mean pairwise differences, which is similar to other Latin American populations, and which represent a marked improvement from the values obtained with just control region sequences. These haplotypes belong to 54 different haplogroups, being 91% of them of Native American origin. Over a third (35.9%) of the individuals carried at least a heteroplasmic site (excluding length heteroplasmies). Ultimately, the present database aims to represent mtDNA haplotype diversity in the general Salvadoran populations as a basis for the identification of people that disappeared during or after the civil war., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Authors. Published by Elsevier B.V. All rights reserved.)

Despite being enclosed between the Mediterranean Sea and the Sahara Desert, North Africa has been the scenario of multiple human migrations that have shaped the genetic structure of its present-day populations. Despite its richness, North Africa remains underrepresented in genomic studies. To overcome this, we have sequenced and analyzed 264 mitogenomes from the Algerian Chaoui-speaking Imazighen (a.k.a. Berbers) living in the Aurès region. The maternal genetic composition of the Aurès is similar to Arab populations in the region, dominated by West Eurasian lineages with a moderate presence of M1/U6 North African and L sub-Saharan lineages. When focusing on the time and geographic origin of the North African specific clades within the non-autochthonous haplogroups, different geographical neighboring regions contributed to the North African maternal gene pool during time periods that could be attributed to previously suggested admixture events in the region, since Paleolithic times to recent historical movements such as the Arabization. We have also observed the role of North Africa as a source of geneflow mainly in Southern European regions since Neolithic times. Finally, the present work constitutes an effort to increase the representation of North African populations in genetic databases, which is key to understand their history., (© 2023. The Author(s).)

The Y chromosome can yield a unique perspective into the study of human demographic history. However, due to the repetitive nature of part of its sequence, only a small set of regions are suitable for variant calling and discovery from short-read sequencing data. These regions combined represent 8.9 Mbp or 0.14% of a diploid human genome. Consequently, investing in whole-genome sequencing to resolve Y-chromosome questions is poorly efficient. Here we use, as an alternative, target enrichment technology to greatly increase sequencing effectiveness, validating and applying the technique to 181 males, for 162 of whom we obtained a positive result. Additionally, 75 samples sequenced for the whole genome were also included, for a total sample size of 237. These samples were chosen for their Y chromosome haplogroup: R1b-DF27. In the context of European populations, and particularly in Iberia, this haplogroup stands out for its high frequency and its demographic history. Current evidence indicates that the diffusion of this haplogroup is related to the population movements that mark the cultural Bronze Age transition, making it remarkably interesting for population geneticists. The results of this study show the effects of the rapid radiation of the haplogroup in Spain, as even with the higher discriminating power of whole sequences, most haplotypes still fall within the R1b-DF27* paragroup rather than in the main derived branches. However, we were able to refine the ISOGG 2019-2020 phylogeny, and its two main subbranches, namely L176.2 and Z272, which present geographical differentiation between the Atlantic and Mediterranean coasts of Iberia., (© 2022. The Author(s).)

The Roma are a group of populations with a common origin that share the Romani identity and cultural heritage. Their genetic history has been inferred through multiple studies based on uniparental and autosomal markers, and current genomic data have provided novel insights into their genetic background. This review was prompted by two factors: (i) new developments to estimate the genetic structure of the Roma at a fine-scale resolution have precisely identified the ancestral components and traced migrations that were previously documented only in historical sources, clarifying and solving debates on the origins and the diaspora of the Roma; (ii) while there has been an effort to review the health determinants of the Roma, the increasing literature on their population genetics has not been subjected to a dedicated review in the last two decades. We believe that a summary on the state of the art will benefit both the public and scholars that are approaching the subject.

The Roma are the largest ethnic minority in Europe. With a Northwestern Indian origin around ~ 1.5 kya, they travelled throughout West Asia until their arrival in Europe around the eleventh century CE. Their diaspora through Europe is characterized by population bottlenecks and founder events which have contributed to their present day genetic and cultural diversity. In our study, we focus on the effects of founder effects in the mitochondrial DNA (mtDNA) pool of Iberian Roma by producing and analyzing 144 novel whole mtDNA sequences of Iberian Roma. Over 60% of their mtDNA pool is composed by founder lineages of South Asian origin or acquired by gene flow during their diaspora in the Middle East or locally in Europe in Europe. The TMRCA of these lineages predates the historical record of the Roma arrival in Spain. The abundance of founder lineages is in contrast with ~ 0.7% of autochthonous founder lineages present in the non-Roma Iberian population. Within those founder lineages, we found a substantial amount of South Asian M5a1b1a1 haplotypes and high frequencies of West Eurasian founder lineages (U3b1c, J2b1c, J1c1b, J1b3a, H88, among others), which we characterized phylogenetically and put in phylogeographical context. Besides, we found no evidence of genetic substructure of Roma within the Iberian Peninsula. These results show the magnitude of founder effects in the Iberian Roma and further explain the Roma history and genetic diversity from a matrilineal point of view., (© 2022. The Author(s).)

Genomic reference databases of residing populations are available in different countries and regions. Since they represent the whole genetic diversity of a geographical region, they have wide applications, from biomedical studies to forensic identifications. Uniparentally transmitted portions of the genome specifically are highly suitable for kinship analyses, mixed DNA cases and geographical ancestry inferences. We have sampled 808 individuals currently residing in Catalonia within the GCAT cohort, from which we have generated 808 high-quality whole mitochondrial DNA (mtDNA) genomes and 399 sequences of the male-specific part of the Y chromosome (MSY). We observe higher genetic diversity than in classical population genetics datasets. We test the robustness of whole sequences for unequivocal identifications, and we found that they have higher resolution than mitochondrial control region and Y chromosome short tandem repeats (Y-STRs), and that most of the variants they present are at low frequencies, increasing the discrimination capacity between individuals. These results confirm the forensic applicability of whole uniparental sequences and provide one of the largest high-quality reference datasets ever published., (Copyright © 2022 The Authors. Published by Elsevier B.V. All rights reserved.)

Background: STR allele frequency databases from populations are necessary to take full advantage of the increased power of discrimination offered by massively parallel sequencing (MPS) platforms., Material and Methods: For this reason, we sequenced 58 STRs (aSTRs, X-STRs, and Y-STRs) and 94 identity informative SNPs (iiSNPs) on 105 Mestizo (admixed) individuals from Monterrey City (Northeast, Mexico), with the Primer Set-A of the ForenSeq™ DNA Signature Prep Kit., Results: Most of the STR markers were in Hardy Weinberg equilibrium, with a few exceptions. We found 346 different length-based alleles for these 58 STRs; nevertheless, they became 528 alleles when the sequence was assessed. The combined power of discrimination from autosomal STRs (aSTRs) was -virtually- 100% in both length and sequence-based alleles, while the power of exclusion was 99.9999999976065 and 99.9999999999494%, respectively. Haplotypes based on X-STRs and Y-STRs showed 100% of discriminatory capacity., Conclusions: These results provide -for the first time- forensic genomic population data from Mexico necessary for interpretation in kinship and criminal analyses., (© 2022. The Author(s), under exclusive licence to Springer Nature B.V.)

The coping style of an individual in relation to potentially dangerous situations has been suggested to be inherited in a polygenic fashion, SERT being one of the candidate genes. In this study, we assessed in free-living great tits (Parus major) the association between SNP290 in the SERT promoter and three standard fear-related behaviors: the response of the birds to a black-and-white flag fixed to the top of the nest box, distress calling rate of the birds in the hand once captured and the hissing call of incubating females when approached by a predator. We found a strong association between SNP290 polymorphism and the three risk-taking behaviors, with birds with genotype CT entering the nest box with the flag faster and displaying more distress calls and fewer hissing calls. CT birds could therefore be described as more proactive than CC individuals. These results also suggest that hissing behavior should be regarded as a fear-induced shy behavior, and confirm that SERT has an important function in relation to risk aversion behaviors and coping style., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2022. Published by The Company of Biologists Ltd.)

We have genotyped the 58 STRs (27 autosomal, 24 Y-STRs and 7 X-STRs) and 94 autosomal SNPs in Illumina ForenSeq™ Primer Mix A in a sample of 248 men and 143 women from El Salvador, Central America. Regional division (Centro, Oriente, Occidente) showed in almost all cases F ST values not significantly different from 0, and further analyses were applied only to the undivided, country-wide population. The overall random match probability (RMP) decreased from 6.79 × 10 -31 in length-based genotypes in the 27 autosomal STRs to 1.47 × 10 -34 in repeat-sequence based genotypes. Combining the autosomal loci in this set, RMP reaches 2.97 × 10 -70 . In a population genetic analysis, El Salvador showed the lowest F ST values with US Hispanics both for autosomal and X-STRs; however, it was much closer to Native Americans for the latter than for the former, in accordance with the well-known gender-biased admixture that created most Latin American populations., (Copyright © 2021 Elsevier B.V. All rights reserved.)

Human populations are genetically affected by their demographic history, which shapes the distribution of their functional genomic variation. However, the genetic impact of recent demography is debated. This issue has been studied in different populations, but never in North Africans, despite their relevant cultural and demographic diversity. In this study we address the question by analyzing new whole-exome sequences from two culturally different Tunisian populations, an isolated Amazigh population and a close non-isolated Arab-speaking population, focusing on the distribution of functional variation. Both populations present clear differences in their variant frequency distribution, in general and for putatively damaging variation. This suggests a relevant effect in the Amazigh population of genetic isolation, drift, and inbreeding, pointing to relaxed purifying selection. We also discover the enrichment in Imazighen of variation associated to specific diseases or phenotypic traits, but the scarce genetic and biomedical data in the region limits further interpretation. Our results show the genomic impact of recent demography and reveal a clear genetic differentiation probably related to culture. These findings highlight the importance of considering cultural and demographic heterogeneity within North Africa when defining population groups, and the need for more data to improve knowledge on the region's health and disease landscape., (© 2021. The Author(s).)

The Spanish and Portuguese-Speaking Working Group of the International Society for Forensic Genetics (GHEP-ISFG) has organized a second collaborative exercise on a simulated case of Disaster Victim Identification (DVI), with the participation of eighteen laboratories. The exercise focused on the analysis of a simulated plane crash case of medium-size resulting in 66 victims with varying degrees of fragmentation of the bodies (with commingled remains). As an additional difficulty, this second exercise included 21 related victims belonging to 6 families among the 66 missings to be identified. A total number of 228 post-mortem samples were represented with aSTR and mtDNA profiles, with a proportion of partial aSTR profiles simulating charred remains. To perform the exercise, participants were provided with aSTR and mtDNA data of 51 reference pedigrees -some of which deficient-including 128 donors for identification purposes. The exercise consisted firstly in the comparison of the post-mortem genetic profiles in order to re-associate fragmented remains to the same individual and secondly in the identification of the re-associated remains by comparing aSTR and mtDNA profiles with reference pedigrees using pre-established thresholds to report a positive identification. Regarding the results of the post-mortem samples re-associations, only a small number of discrepancies among participants were detected, all of which were from just a few labs. However, in the identification process by kinship analysis with family references, there were more discrepancies in comparison to the correct results. The identification results of single victims yielded fewer problems than the identification of multiple related victims within the same family groups. Several reasons for the discrepant results were detected: a) the identity/non-identity hypotheses were sometimes wrongly expressed in the likelihood ratio calculations, b) some laboratories failed to use all family references to report the DNA match, c) In families with several related victims, some laboratories firstly identified some victims and then unnecessarily used their genetic information to identify the remaining victims within the family, d) some laboratories did not correctly use "prior odds" values for the Bayesian treatment of the episode for both post-mortem/post-mortem re-associations as well as the ante-mortem/post-mortem comparisons to evaluate the probability of identity. For some of the above reasons, certain laboratories failed to identify some victims. This simulated "DNA-led" identification exercise may help forensic genetic laboratories to gain experience and expertize for DVI or MPI in using genetic data and comparing their own results with the ones in this collaborative exercise., (Copyright © 2021 Elsevier B.V. All rights reserved.)

Demographic history plays a major role in shaping the distribution of genomic variation. Yet the interaction between different demographic forces and their effects in the genomes is not fully resolved in human populations. Here, we focus on the Roma population, the largest transnational ethnic minority in Europe. They have a South Asian origin and their demographic history is characterized by recent dispersals, multiple founder events, and extensive gene flow from non-Roma groups. Through the analyses of new high-coverage whole exome sequences and genome-wide array data for 89 Iberian Roma individuals together with forward simulations, we show that founder effects have reduced their genetic diversity and proportion of rare variants, gene flow has counteracted the increase in mutational load, runs of homozygosity show ancestry-specific patterns of accumulation of deleterious homozygotes, and selection signals primarily derive from preadmixture adaptation in the Roma population sources. The present study shows how two demographic forces, bottlenecks and admixture, act in opposite directions and have long-term balancing effects on the Roma genomes. Understanding how demography and gene flow shape the genome of an admixed population provides an opportunity to elucidate how genomic variation is modeled in human populations., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Genetic patterns of inter-population variation are a result of different demographic and adaptive histories, which gradually shape the frequency distribution of the variants. However, the study of clinically relevant mutations has a Eurocentric bias. The Romani, the largest transnational minority ethnic group in Europe, originated in South Asia and received extensive gene flow from West Eurasia. Most medical genetic studies have only explored founder mutations related to Mendelian disorders in this population. Here we analyze exome sequences and genome-wide array data of 89 healthy Spanish Roma individuals to study complex traits and disease. We apply a different framework and focus on variants with both increased and decreased allele frequencies, taking into account their local ancestry. We report several OMIM traits enriched for genes with deleterious variants showing increased frequencies in Roma or in non-Roma (e.g., obesity is enriched in Roma, with an associated variant linked to South Asian ancestry; while non-insulin dependent diabetes is enriched in non-Roma Europeans). In addition, previously reported pathogenic variants also show differences among populations, where some variants segregating at low frequency in non-Roma are virtually absent in the Roma. Lastly, we describe frequency changes in drug-response variation, where many of the variants increased in Roma are clinically associated with metabolic and cardiovascular-related drugs. These results suggest that clinically relevant variation in Roma cannot only be characterized in terms of founder mutations. Instead, we observe frequency differences compared to non-Roma: some variants are absent, while other have drifted to higher frequencies. As a result of the admixture events, these clinically damaging variants can be traced back to both European and South Asian-related ancestries. This can be attributed to a different prevalence of some genetic disorders or to the fact that genetic susceptibility variants are mostly studied in populations of European descent, and can differ in individuals with different ancestries., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Font-Porterias, Giménez, Carballo-Mesa, Calafell and Comas.)

Basques have historically lived along the Western Pyrenees, in the Franco-Cantabrian region, straddling the current Spanish and French territories. Over the last decades, they have been the focus of intense research due to their singular cultural and biological traits that, with high controversy, placed them as a heterogeneous, isolated, and unique population. Their non-Indo-European language, Euskara, is thought to be a major factor shaping the genetic landscape of the Basques. Yet there is still a lively debate about their history and assumed singularity due to the limitations of previous studies. Here, we analyze genome-wide data of Basque and surrounding groups that do not speak Euskara at a micro-geographical level. A total of ∼629,000 genome-wide variants were analyzed in 1,970 modern and ancient samples, including 190 new individuals from 18 sampling locations in the Basque area. For the first time, local- and wide-scale analyses from genome-wide data have been performed covering the whole Franco-Cantabrian region, combining allele frequency and haplotype-based methods. Our results show a clear differentiation of Basques from the surrounding populations, with the non-Euskara-speaking Franco-Cantabrians located in an intermediate position. Moreover, a sharp genetic heterogeneity within Basques is observed with significant correlation with geography. Finally, the detected Basque differentiation cannot be attributed to an external origin compared to other Iberian and surrounding populations. Instead, we show that such differentiation results from genetic continuity since the Iron Age, characterized by periods of isolation and lack of recent gene flow that might have been reinforced by the language barrier., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Currently, there are 18 different religious communities living in Lebanon. While evolving primarily within Lebanon, these communities show a level of local isolation as demonstrated previously from their Y-haplogroup distributions. In order to trace the origins and migratory patterns that may have led to the genetic isolation and autosomal clustering in some of these communities we analyzed Y-chromosome STR and SNP sample data from 6327 individuals, in addition to whole genome autosomal sample data from 609 individuals, from Mount Lebanon and other surrounding communities. We observed Y chromosome L1b Levantine STR branching that occurred around 5000 years ago. Autosomal DNA analyses suggest that the North Lebanese Mountain Maronite community possesses an ancestral Fertile Crescent genetic component distinct from other populations in the region. We suggest that the Levantine L1b group split from the Caucasus ancestral group around 7300 years ago and migrated to the Levant. This event was distinct from the earlier expansions from the Caucasus region that contributed to the wider Levantine populations. Differential cultural adaption by populations from the North Lebanese Mountains are clearly aligned with the L1b haplotype STR haplogroup clusters, indicating pre-existing and persistent cultural barriers marked by the transmission of L1b lineages. Our findings highlight the value of uniparental haplogroups and STR haplotype data for elucidating biosocial events among these populations.

The Roma Diaspora-traditionally known as Gypsies-remains among the least explored population migratory events in historical times. It involved the migration of Roma ancestors out-of-India through the plateaus of Western Asia ultimately reaching Europe. The demographic effects of the Diaspora-bottlenecks, endogamy, and gene flow-might have left marked molecular traces in the Roma genomes. Here, we analyze the whole-genome sequence of 46 Roma individuals pertaining to four migrant groups in six European countries. Our analyses revealed a strong, early founder effect followed by a drastic reduction of ∼44% in effective population size. The Roma common ancestors split from the Punjabi population, from Northwest India, some generations before the Diaspora started, <2,000 years ago. The initial bottleneck and subsequent endogamy are revealed by the occurrence of extensive runs of homozygosity and identity-by-descent segments in all Roma populations. Furthermore, we provide evidence of gene flow from Armenian and Anatolian groups in present-day Roma, although the primary contribution to Roma gene pool comes from non-Roma Europeans, which accounts for >50% of their genomes. The linguistic and historical differentiation of Roma in migrant groups is confirmed by the differential proportion, but not a differential source, of European admixture in the Roma groups, which shows a westward cline. In the present study, we found that despite the strong admixture Roma had in their diaspora, the signature of the initial bottleneck and the subsequent endogamy is still present in Roma genomes., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

The Roma population is a European ethnic minority characterized by recent and multiple dispersals and founder effects. After their origin in South Asia around 1,500 years ago, they migrated West. In Europe, they diverged into ethnolinguistically distinct migrant groups that spread across the continent. Previous genetic studies based on genome-wide data and uniparental markers detected Roma founder events and West-Eurasian gene flow. However, to the best of our knowledge, it has not been assessed whether these demographic processes have equally affected both sexes in the population. The present study uses the largest and most comprehensive dataset of complete mitochondrial and Y chromosome Roma sequences to unravel the sex-biased patterns that have shaped their genetic history. The results show that the Roma maternal genetic pool carries a higher lineage diversity from South Asia, as opposed to a single paternal South Asian lineage. Nonetheless, the European gene flow events mainly occurred through the maternal lineages; however, a signal of this gene flow is also traceable in the paternal lineages. We also detect a higher female migration rate among European Roma groups. Altogether, these results suggest that sociocultural factors influenced the emergence of sex-biased genetic patterns at global and local scales in the Roma population through time.

Unlike other European countries, the human population genetics and demographic history of Metropolitan France is surprisingly understudied. In this work, we combined newly genotyped samples from various zones in France with publicly available data and applied both allele frequency and haplotype-based methods to describe the internal structure of this country, using genome-wide single nucleotide polymorphism (SNP) array genotypes. We found out that French Basques, already known for their linguistic uniqueness, are genetically distinct from all other groups and that the populations from southwest France (namely the Gascony region) share a large proportion of their ancestry with Basques. Otherwise, the genetic makeup of the French population is relatively homogeneous and mostly related to Southern and Central European groups. However, a fine-grained, haplotype-based analysis revealed that Bretons slightly separated from the rest of the groups, due mostly to gene flow from the British Isles in a time frame that coincides both historically attested Celtic population movements to this area between the 3th and the ninth centuries CE, but also with a more ancient genetic continuity between Brittany and the British Isles related to the shared drift with hunter-gatherer populations. Haplotype-based methods also unveiled subtle internal structures and connections with the surrounding modern populations, particularly in the periphery of the country.

Paternity testing using genetic markers has shown that extra-pair paternity (EPP) is common in many pair-bonded species [1, 2]. Evolutionary theory and empirical data show that extra-pair copulations can increase the fitness of males as well as females [3, 4]. This can carry a significant fitness cost for the social father, who then invests in rearing offspring that biologically are not his own [5]. In human populations, the incidence and correlates of extra-pair paternity remain highly contentious [2, 6, 7]. Here, we use a population-level genetic genealogy approach [6, 8] to reconstruct spatiotemporal patterns in human EPP rates. Using patrilineal genealogies from the Low Countries spanning a period of over 500 years and Y chromosome genotyping of living descendants, our analysis reveals that historical EPP rates, while low overall, were strongly impacted by socioeconomic and demographic factors. Specifically, we observe that estimated EPP rates among married couples varied by more than an order of magnitude, from 0.4% to 5.9%, and peaked among families with a low socioeconomic background living in densely populated cities of the late 19 th century. Our results support theoretical predictions that social context can strongly affect the outcomes of sexual conflict in human populations by modulating the incentives and opportunities for engaging in extra-pair relationships [9-11]. These findings show how contemporary genetic data combined with in-depth genealogies open up a new window on the sexual behavior of our ancestors., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

The Roma population is the largest transnational ethnic minority in Europe, characterized by a linguistic, cultural and historical heterogeneity. Comparative linguistics and genetic studies have placed the origin of European Roma in the Northwest of India. After their migration across Persia, they entered into the Balkan Peninsula, from where they spread into Europe, arriving in the Iberian Peninsula in the 15th century. Their particular demographic history has genetic implications linked to rare and common diseases. However, the South Asian source of the proto-Roma remains still untargeted and the West Eurasian Roma component has not been yet deeply characterized. Here, in order to describe both the South Asian and West Eurasian ancestries, we analyze previously published genome-wide data of 152 European Roma and 34 new Iberian Roma samples at a fine-scale and haplotype-based level, with special focus on the Iberian Roma genetic substructure. Our results suggest that the putative origin of the proto-Roma involves a Punjabi group with low levels of West Eurasian ancestry. In addition, we have identified a complex West Eurasian component (around 65%) in the Roma, as a result of the admixture events occurred with non-proto-Roma populations between 1270-1580. Particularly, we have detected the Balkan genetic footprint in all European Roma, and the Baltic and Iberian components in the Northern and Western Roma groups, respectively. Finally, our results show genetic substructure within the Iberian Roma, with different levels of West Eurasian admixture, as a result of the complex historical events occurred in the Peninsula., Competing Interests: The authors have declared that no competing interests exist.

The analysis of geographically specific regions and the characterization of fine-scale patterns of genetic diversity may facilitate a much better understanding of the microevolutionary processes affecting local human populations. Here we generated genome-wide high-density SNP genotype data in 425 individuals from six geographical regions in Lithuania and combined our dataset with available ancient and modern data to explore genetic population structure, ancestry components and signatures of natural positive selection in the Lithuanian population. Our results show that Lithuanians are a homogenous population, genetically differentiated from neighbouring populations but within the general expected European context. Moreover, we not only confirm that Lithuanians preserve one of the highest proportions of western, Scandinavian and eastern hunter-gather ancestry components found in European populations but also that of an steppe Early to Middle Bronze Age pastoralists, which together configure the genetic distinctiveness of the Lithuanian population. Finally, among the top signatures of positive selection detected in Lithuanians, we identified several candidate genes related with diet (PNLIP, PPARD), pigmentation (SLC24A5, TYRP1 and PPARD) and the immune response (BRD2, HLA-DOA, IL26 and IL22).

In this study, we seek to understand and to correlate the genetic patterns observed in the population of the island of Ibiza in the Western Mediterranean basin with past events. Genome-wide genotypes of 189 samples representing 13 of 17 regions in Spain have been analyzed, in addition to 105 samples from the Levant, 157 samples from North Africa, and one ancient sample from the Phoenician Cas Molí site in Ibiza. Before the Catalans conquered the island in 1235 CE, Ibiza (Eivissa) had already been influenced by several cultures, starting with the Phoenicians, then the Carthaginians, followed by the Umayyads. The impact of these various cultures on the genetic structure of the islanders is still unexplored. Our results show a clear distinction between Ibiza and the rest of Spain. To investigate whether this was due to the Phoenician colonization or to more recent events, we compared the genomes of current Ibizans to that of an ancient Phoenician sample from Ibiza and to both modern Levantine and North African genomes. We did not identify any trace of Phoenician ancestry in the current Ibizans. Interestingly, the analysis of runs of homozygosity and changes in the effective population size through time support the idea that drift has shaped the genetic structure of current Ibizans. In addition to the small carrying capacity of the island, Ibiza experienced a series of dramatic demographic changes due to several instances of famine, war, malaria and plague that could have significantly contributed to its current genetic differentiation.

Ibiza was permanently settled around the 7 th century BCE by founders arriving from west Phoenicia. The founding population grew significantly and reached its height during the 4 th century BCE. We obtained nine complete mitochondrial genomes from skeletal remains from two Punic necropoli in Ibiza and a Bronze Age site from Formentara. We also obtained low coverage (0.47X average depth) of the genome of one individual, directly dated to 361-178 cal BCE, from the Cas Molí site on Ibiza. We analysed and compared ancient DNA results with 18 new mitochondrial genomes from modern Ibizans to determine the ancestry of the founders of Ibiza. The mitochondrial results indicate a predominantly recent European maternal ancestry for the current Ibizan population while the whole genome data suggest a significant Eastern Mediterranean component. Our mitochondrial results suggest a genetic discontinuity between the early Phoenician settlers and the island's modern inhabitants. Our data, while limited, suggest that the Eastern or North African influence in the Punic population of Ibiza was primarily male dominated.

Basques show specific cultural, demographic, and genetic characteristics that have placed them as an isolated and unique population within Europe, such as their non-Indo-European language, Euskara. They have historically lived along the Western Pyrenees, between Spain and France, in one of the most important European glacial refugia during the Last Glacial Maximum. The most striking genetic characteristic is their highest frequency of the RhD blood group negative allele, a variant related to the hemolytic disease of the newborn. Both demographic and adaptive processes have been suggested as possible causes of the high frequency of RhD negative in Basques, but neither hypothesis has been clearly demonstrated. While previous studies on the Rh system in Basques have been mostly focused on serological and genotyping diversity, in this work we analyze genotyping and next generation sequencing data in order to provide a general framework of the genetic scenario of the system in Basques. In particular, we genotyped the most relevant variants of the system (D/d, E/e, and C/c), and sequenced three ~6 kb flanking regions surrounding the Rh genes in Basques and also in other populations for comparison. Our results are in agreement with previous studies, with Basques presenting the highest frequency of the RHD deletion (47.2%). Haplotype analyses of D/d, E/e, and C/c variants confirmed an association between the RhC allele, previously suggested to be under positive selection, and the RhD positive variant in non-sub-Saharan populations, including Basques. We also found extreme differentiation for the C/c variant when comparing sub-Saharan to non-sub-Saharan populations.

Objectives: War atrocities committed by the Spanish army in the Low Countries during the 16th century are so ingrained in the collective memory of Belgian and Dutch societies that they generally assume a signature of this history to be present in their genetic ancestry. Historians claim this assumption is a consequence of the so-called "Black Legend" and negative propaganda portraying and remembering Spanish soldiers as extreme sexual aggressors. The impact of the presence of Spaniards during the Dutch Revolt on the genetic variation in the Low Countries has been verified in this study., Materials and Methods: A recent population genetic analysis of Iberian-associated Y-chromosomal variation among Europe is enlarged with representative samples of Dutch (N = 250) and Flemish (N = 1,087) males. Frequencies of these variants are also compared between donors whose oldest reported paternal ancestors lived in-nowadays Flemish-cities affected by so-called Spanish Furies (N = 116) versus other patrilineages in current Flemish territory (N = 971)., Results: The frequencies of Y-chromosomal markers Z195 and SRY2627 decline steeply going north from Spain and the data for the Flemish and Dutch populations fits within this pattern. No trend of higher frequencies of these variants has been found within the well-ascertained samples associated with Spanish Fury cities., Discussion: Although sexual aggression did occur in the 16th century, these activities did not leave a traceable "Spanish" genetic signature in the autochthonous genome of the Low Countries. Our results support the view that the 'Black Legend' and historical propaganda on sexual aggression have nurtured today's incorrect assumptions regarding genetic ancestry., (© 2018 Wiley Periodicals, Inc.)

Background: The genetic composition of human North African populations is an amalgam of different ancestral components coming from the Middle East, Europe, south-Saharan Africa and autochthonous to North Africa. This complex genetic pattern is the result of migrations and admixtures in the region since Palaeolithic times., Aims: The objective of the present study is to refine knowledge of the population history of North African populations through the analysis of complete mitochondrial sequences., Subjects and Methods: This study has sequenced complete mitochondrial DNAs (mtDNAs) in several North African and neighbouring individuals., Results: The mtDNA haplogroup classification and phylogeny shows a high genetic diversity in the region as a result of continuous admixture. The phylogenetic analysis allowed us to identify a new haplogroup characterised by positions 10 101 C and 146 C (H1v2), a sub-branch of H1v, which is restricted to North Africa and whose origins are estimated as ∼4000 years ago., Conclusions: The analysis of the complete mtDNA genome has allowed for the identification of a North African sub-lineage that might be ignored by the analysis of partial mtDNA control region sequences, highlighting the phylogeographic relevance of mtDNA complete sequence analysis.

E-M183 (E-M81) is the most frequent paternal lineage in North Africa and thus it must be considered to explore past historical and demographical processes. Here, by using whole Y chromosome sequences from 32 North African individuals, we have identified five new branches within E-M183. The validation of these variants in more than 200 North African samples, from which we also have information of 13 Y-STRs, has revealed a strong resemblance among E-M183 Y-STR haplotypes that pointed to a rapid expansion of this haplogroup. Moreover, for the first time, by using both SNP and STR data, we have provided updated estimates of the times-to-the-most-recent-common-ancestor (TMRCA) for E-M183, which evidenced an extremely recent origin of this haplogroup (2,000-3,000 ya). Our results also showed a lack of population structure within the E-M183 branch, which could be explained by the recent and rapid expansion of this haplogroup. In spite of a reduction in STR heterozygosity towards the West, which would point to an origin in the Near East, ancient DNA evidence together with our TMRCA estimates point to a local origin of E-M183 in NW Africa.

We have genotyped the 58 STRs (27 autosomal, 24 Y-STRs and 7 X-STRs) and 94 autosomal SNPs in Illumina ForenSeq™ Primer Mix A in 88 Spanish Roma (Gypsy) samples and 143 Catalans. Since this platform is based in massive parallel sequencing, we have used simple R scripts to uncover the sequence variation in the repeat region. Thus, we have found, across 58 STRs, 541 length-based alleles, which, after considering repeat-sequence variation, became 804 different alleles. All loci in both populations were in Hardy-Weinberg equilibrium. F ST between both populations was 0.0178 for autosomal SNPs, 0.0146 for autosomal STRs, 0.0101 for X-STRs and 0.1866 for Y-STRs. Combined a priori statistics showed quite large; for instance, pooling all the autosomal loci, the a priori probabilities of discriminating a suspect become 1-(2.3×10 -70 ) and 1-(5.9×10 -73 ), for Roma and Catalans respectively, and the chances of excluding a false father in a trio are 1-(2.6×10 -20 ) and 1-(2.0×10 -21 )., (Copyright © 2017 Elsevier B.V. All rights reserved.)

Haplogroup R1b-M269 comprises most Western European Y chromosomes; of its main branches, R1b-DF27 is by far the least known, and it appears to be highly prevalent only in Iberia. We have genotyped 1072 R1b-DF27 chromosomes for six additional SNPs and 17 Y-STRs in population samples from Spain, Portugal and France in order to further characterize this lineage and, in particular, to ascertain the time and place where it originated, as well as its subsequent dynamics. We found that R1b-DF27 is present in frequencies ~40% in Iberian populations and up to 70% in Basques, but it drops quickly to 6-20% in France. Overall, the age of R1b-DF27 is estimated at ~4,200 years ago, at the transition between the Neolithic and the Bronze Age, when the Y chromosome landscape of W Europe was thoroughly remodeled. In spite of its high frequency in Basques, Y-STR internal diversity of R1b-DF27 is lower there, and results in more recent age estimates; NE Iberia is the most likely place of origin of DF27. Subhaplogroup frequencies within R1b-DF27 are geographically structured, and show domains that are reminiscent of the pre-Roman Celtic/Iberian division, or of the medieval Christian kingdoms.

Background: The co-occurrence of the c.709-1G>A GRN mutation and the p.A152T MAPT variant has been identified in 18 Basque families affected by frontotemporal dementia (FTD). We aimed to investigate the influence of the p.A152T MAPT variant on the clinical and neuropathological features of these Basque GRN families., Methods and Findings: We compared clinical characteristics of 14 patients who carried the c.709-1G>A GRN mutation (GRN+/A152T-) with 21 patients who carried both the c.709-1G>A GRN mutation and the p.A152T MAPT variant (GRN+/A152T+). Neuropsychological data (n = 17) and plasma progranulin levels (n = 23) were compared between groups, and 7 subjects underwent neuropathological studies. We genotyped six short tandem repeat markers in the two largest families. By the analysis of linkage disequilibrium decay in the haplotype block we estimated the time when the first ancestor to carry both genetic variants emerged. GRN+/A152T+ and GRN+/A152T- patients shared similar clinical and neuropsychological features and plasma progranulin levels. All were diagnosed with an FTD disorder, including behavioral variant FTD or non fluent / agrammatic variant primary progressive aphasia, and shared a similar pattern of neuropsychological deficits, predominantly in executive function, memory, and language. All seven participants with available brain autopsies (6 GRN+/A152T+, 1 GRN+/A152T-) showed frontotemporal lobar degeneration with TDP-43 inclusions (type A classification), which is characteristic of GRN carriers. Additionally, all seven showed mild to moderate tau inclusion burden: five cases lacked β-amyloid pathology and two cases had Alzheimer's pathology. The co-occurrence of both genes within one individual is recent, with the birth of the first GRN+/A152T+ individual estimated to be within the last 50 generations (95% probability)., Conclusions: In our sample, the p.A152T MAPT variant does not appear to show a discernible influence on the clinical phenotype of GRN carriers. Whether p.A152T confers a greater than expected propensity for tau pathology in these GRN carriers remains an open question.

The Y chromosome is currently by far the most popular marker in genetic genealogy that combines genetic data and family history. This popularity is based on its haploid character and its close association with the patrilineage and paternal inherited surname. Other markers have not been found (yet) to overrule this status due to the low sensitivity and precision of autosomal DNA for genetic genealogical applications, given the vagaries of recombination, and the lower capacities of mitochondrial DNA combined with an in general much lower interest in maternal lineages. The current knowledge about the Y chromosome and the availability of markers with divergent mutation rates make it possible to answer questions on relatedness levels which differ in time depth; from the individual and familial level to the surnames, clan and population level. The use of the Y chromosome in genetic genealogy has led to applications in several well-established research disciplines; namely in, e.g., family history, demography, anthropology, forensic sciences, population genetics and sex chromosome evolution. The information obtained from analysing this chromosome is not only interesting for academic scientists but also for the huge and lively community of amateur genealogists and citizen-scientists, fascinated in analysing their own genealogy or surname. This popularity, however, has also some drawbacks, mainly for privacy reasons related to the DNA donor, his close family and far-related namesakes. In this review paper we argue why Y-chromosomal analysis and its genetic genealogical applications will still perform an important role in future interdisciplinary research.

We present 42 new Y-chromosomal sequences from diverse Indian tribal and non-tribal populations, including the Jarawa and Onge from the Andaman Islands, which are analysed within a calibrated Y-chromosomal phylogeny incorporating South Asian (in total 305 individuals) and worldwide (in total 1286 individuals) data from the 1000 Genomes Project. In contrast to the more ancient ancestry in the South than in the North that has been claimed, we detected very similar coalescence times within Northern and Southern non-tribal Indian populations. A closest neighbour analysis in the phylogeny showed that Indian populations have an affinity towards Southern European populations and that the time of divergence from these populations substantially predated the Indo-European migration into India, probably reflecting ancient shared ancestry rather than the Indo-European migration, which had little effect on Indian male lineages. Among the tribal populations, the Birhor (Austro-Asiatic-speaking) and Irula (Dravidian-speaking) are the nearest neighbours of South Asian non-tribal populations, with a common origin in the last few millennia. In contrast, the Riang (Tibeto-Burman-speaking) and Andamanese have their nearest neighbour lineages in East Asia. The Jarawa and Onge shared haplogroup D lineages with each other within the last ~7000 years, but had diverged from Japanese haplogroup D Y-chromosomes ~53000 years ago, most likely by a split from a shared ancestral population. This analysis suggests that Indian populations have complex ancestry which cannot be explained by a single expansion model.

The European paternal lineage R-DF27 has been proposed as a haplogroup of Iberian origin due to its maximum frequencies in the Iberian Peninsula. In this study, the distribution and structure of DF27 were characterized in 591 unrelated male individuals from four key populations of the north area of the Iberian Peninsula through the analysis of 12 Y-SNPs that define DF27 main sublineages. Additionally, Y-SNP allele frequencies were also gathered from the reference populations in the 1000 Genomes Project to compare and obtain a better landscape of the distribution of DF27. Our results reveal frequencies over 35% of DF27 haplogroup in the four North Iberian populations analyzed and high frequencies for its subhaplogroups. Considering the low frequency of DF27 and its sublineages in most populations outside of the Iberian Peninsula, this haplogroup seems to have geographical significance; thus, indicating a possible Iberian patrilineal origin of vestiges bearing this haplogroup. The dataset presented here contributes with new data to better understand the complex genetic variability of the Y chromosome in the Iberian Peninsula, that can be applied in Forensic Genetics., (Copyright © 2016 Elsevier Ireland Ltd. All rights reserved.)

Next-generation sequencing technologies have opened new opportunities in forensic genetics. Here, we assess the applicability and performance of the MiSeq FGx™ & ForenSeq™ DNA Signature Prep Kit (Illumina) for the identification of individuals from the mass graves of the Spanish Civil War (1936-1939). The main limitations for individual identification are the low number of possible first-degree living relatives and the high levels of DNA degradation reported in previous studies. Massively parallel sequencing technologies enabling the analysis of hundreds of regions and prioritizing short length amplicons constitute a promising tool for this kind of approaches. In this study, we first explore the power of this new technology to detect first- and second-degree kinship given different scenarios of DNA degradation. Second, we specifically assess its performance in a set of low DNA input samples previously analyzed with CE technologies. We conclude that this methodology will allow identification of up to second-degree relatives, even in situations with low sequencing performance and important levels of allele drop-out; it is thus a technology that resolves previous drawbacks and that will allow a successful approximation to the identification of remains., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Essential trace elements possess vital functions at molecular, cellular, and physiological levels in health and disease, and they are tightly regulated in the human body. In order to assess variability and potential adaptive evolution of trace element homeostasis, we quantified 18 trace elements in 150 liver samples, together with the expression levels of 90 genes and abundances of 40 proteins involved in their homeostasis. Additionally, we genotyped 169 single nucleotide polymorphism (SNPs) in the same sample set. We detected significant associations for 8 protein quantitative trait loci (pQTL), 10 expression quantitative trait loci (eQTLs), and 15 micronutrient quantitative trait loci (nutriQTL). Six of these exceeded the false discovery rate cutoff and were related to essential trace elements: 1) one pQTL for GPX2 (rs10133290); 2) two previously described eQTLs for HFE (rs12346) and SELO (rs4838862) expression; and 3) three nutriQTLs: The pathogenic C282Y mutation at HFE affecting iron (rs1800562), and two SNPs within several clustered metallothionein genes determining selenium concentration (rs1811322 and rs904773). Within the complete set of significant QTLs (which involved 30 SNPs and 20 gene regions), we identified 12 SNPs with extreme patterns of population differentiation (FST values in the top 5% percentile in at least one HapMap population pair) and significant evidence for selective sweeps involving QTLs at GPX1, SELENBP1, GPX3, SLC30A9, and SLC39A8. Overall, this detailed study of various molecular phenotypes illustrates the role of regulatory variants in explaining differences in trace element homeostasis among populations and in the human adaptive response to environmental pressures related to micronutrients., (© The Author 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

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.