The peopling of Europe and the cautionary tale of Y chromosome lineage R-M269.

Recently, the debate on the origins of the major European Y chromosome haplogroup R1b1b2-M269 has reignited, and opinion has moved away from Palaeolithic origins to the notion of a younger Neolithic spread of these chromosomes from the Near East. Here, we address this debate by investigating frequency patterns and diversity in the largest collection of R1b1b2-M269 chromosomes yet assembled. Our analysis reveals no geographical trends in diversity, in contradiction to expectation under the Neolithic hypothesis, and suggests an alternative explanation for the apparent cline in diversity recently described. We further investigate the young, STR-based time to the most recent common ancestor estimates proposed so far for R-M269-related lineages and find evidence for an appreciable effect of microsatellite choice on age estimates. As a consequence, the existing data and tools are insufficient to make credible estimates for the age of this haplogroup, and conclusions about the timing of its origin and dispersal should be viewed with a large degree of caution.

Bulgarians vs the other European populations: a mitochondrial DNA perspective.

To define the matrilineal relationships between Bulgarians and other European populations, we have evaluated the mitochondrial DNA (mtDNA) variation in a sample of 855 Bulgarian subjects from the mtDNA perspective. The molecular survey was performed by sequencing ∼750 bp of the control region, which resulted in 557 different haplotypes, and by a subsequent restriction fragment length polymorphism analysis to confirm haplogroup/subhaplogroup affiliation. The classification was carried out according to the most updated criteria as reported by van Oven and Kayser (Hum Mutat 30:386-394, 2009), allowing the identification of 45 mitochondrial clades. The observed pattern of mtDNA variation indicates that the Bulgarian mitochondrial pool is geographically homogeneous across the country, and that is characterized by an overall extremely high frequency of western Eurasian lineages. In the principal component analysis, Bulgarians locate in an intermediate position between Eastern European and Mediterranean populations, which is in agreement with historical events. Thus, while the Mediterranean legacy could be attributed to the Thracians, indigenous people that firstly inhabited the Balkans, the Eastern contribution is likely due to the Proto-Bulgarians originating from the Middle East and to the Slavs migrating from northeast Europe.

Case-control association study of 59 candidate genes reveals the DRD2 SNP rs6277 (C957T) as the only susceptibility factor for schizophrenia in the Bulgarian population.

The development of molecular psychiatry in the last few decades identified a number of candidate genes that could be associated with schizophrenia. A great number of studies often result with controversial and non-conclusive outputs. However, it was determined that each of the implicated candidates would independently have a minor effect on the susceptibility to that disease. Herein we report results from our replication study for association using 255 Bulgarian patients with schizophrenia and schizoaffective disorder and 556 Bulgarian healthy controls. We have selected from the literatures 202 single nucleotide polymorphisms (SNPs) in 59 candidate genes, which previously were implicated in disease susceptibility, and we have genotyped them. Of the 183 SNPs successfully genotyped, only 1 SNP, rs6277 (C957T) in the DRD2 gene (P=0.0010, odds ratio=1.76), was considered to be significantly associated with schizophrenia after the replication study using independent sample sets. Our findings support one of the most widely considered hypotheses for schizophrenia etiology, the dopaminergic hypothesis.

Y-chromosome diversity in modern Bulgarians: new clues about their ancestry.

To better define the structure and origin of the Bulgarian paternal gene pool, we have examined the Y-chromosome variation in 808 Bulgarian males. The analysis was performed by high-resolution genotyping of biallelic markers and by analyzing the STR variation within the most informative haplogroups. We found that the Y-chromosome gene pool in modern Bulgarians is primarily represented by Western Eurasian haplogroups with ∼ 40% belonging to haplogroups E-V13 and I-M423, and 20% to R-M17. Haplogroups common in the Middle East (J and G) and in South Western Asia (R-L23*) occur at frequencies of 19% and 5%, respectively. Haplogroups C, N and Q, distinctive for Altaic and Central Asian Turkic-speaking populations, occur at the negligible frequency of only 1.5%. Principal Component analyses group Bulgarians with European populations, apart from Central Asian Turkic-speaking groups and South Western Asia populations. Within the country, the genetic variation is structured in Western, Central and Eastern Bulgaria indicating that the Balkan Mountains have been permeable to human movements. The lineage analysis provided the following interesting results: (i) R-L23* is present in Eastern Bulgaria since the post glacial period; (ii) haplogroup E-V13 has a Mesolithic age in Bulgaria from where it expanded after the arrival of farming; (iii) haplogroup J-M241 probably reflects the Neolithic westward expansion of farmers from the earliest sites along the Black Sea. On the whole, in light of the most recent historical studies, which indicate a substantial proto-Bulgarian input to the contemporary Bulgarian people, our data suggest that a common paternal ancestry between the proto-Bulgarians and the Altaic and Central Asian Turkic-speaking populations either did not exist or was negligible.

Whole-genome-wide association study in the Bulgarian population reveals HHAT as schizophrenia susceptibility gene.

OBJECTIVE: Schizophrenia, the most common major psychiatric disorder (or group of disorders), entails severe decline of higher functions, principally with alterations in cognitive functioning and reality perception. Both genetic and environmental factors are involved in its pathogenesis; however, its genetic background still needs to be clarified. The objective of the study was to reveal genetic markers associated with schizophrenia in the Bulgarian population. METHODS: We have conducted a genome-wide association study using 554 496 single nucleotide polymorphisms (SNPs) in 188 affected and 376 unaffected Bulgarian individuals. Subsequently, the 100 candidate SNPs that revealed the smallest P-values were further evaluated in an additional set of 99 case and 328 control samples. RESULTS: We found a significant association between schizophrenia and the intronic SNP rs7527939 in the HHAT gene (P-value of 6.49×10 with an odds ratio of 2.63, 95% confidence interval of 1.89-3.66). We also genotyped additional SNPs within a 58-kb linkage disequilibrium block surrounding the landmark SNP. CONCLUSION: We suggest rs7527939 to be the strongest indicator of susceptibility to schizophrenia in the Bulgarian population within the HHAT locus.

Genetic characterization of northeastern Italian population isolates in the context of broader European genetic diversity.

Population genetic studies on European populations have highlighted Italy as one of genetically most diverse regions. This is possibly due to the country's complex demographic history and large variability in terrain throughout the territory. This is the reason why Italy is enriched for population isolates, Sardinia being the best-known example. As the population isolates have a great potential in disease-causing genetic variants identification, we aimed to genetically characterize a region from northeastern Italy, which is known for isolated communities. Total of 1310 samples, collected from six geographically isolated villages, were genotyped at >145000 single-nucleotide polymorphism positions. Newly genotyped data were analyzed jointly with the available genome-wide data sets of individuals of European descent, including several population isolates. Despite the linguistic differences and geographical isolation the village populations still show the greatest genetic similarity to other Italian samples. The genetic isolation and small effective population size of the village populations is manifested by higher levels of genomic homozygosity and elevated linkage disequilibrium. These estimates become even more striking when the detected substructure is taken into account. The observed level of genetic isolation in Friuli-Venezia Giulia region is more extreme according to several measures of isolation compared with Sardinians, French Basques and northern Finns, thus proving the status of an isolate.

Distinguishing the co-ancestries of haplogroup G Y-chromosomes in the populations of Europe and the Caucasus.

Haplogroup G, together with J2 clades, has been associated with the spread of agriculture, especially in the European context. However, interpretations based on simple haplogroup frequency clines do not recognize underlying patterns of genetic diversification. Although progress has been recently made in resolving the haplogroup G phylogeny, a comprehensive survey of the geographic distribution patterns of the significant sub-clades of this haplogroup has not been conducted yet. Here we present the haplogroup frequency distribution and STR variation of 16 informative G sub-clades by evaluating 1472 haplogroup G chromosomes belonging to 98 populations ranging from Europe to Pakistan. Although no basal G-M201* chromosomes were detected in our data set, the homeland of this haplogroup has been estimated to be somewhere nearby eastern Anatolia, Armenia or western Iran, the only areas characterized by the co-presence of deep basal branches as well as the occurrence of high sub-haplogroup diversity. The P303 SNP defines the most frequent and widespread G sub-haplogroup. However, its sub-clades have more localized distribution with the U1-defined branch largely restricted to Near/Middle Eastern and the Caucasus, whereas L497 lineages essentially occur in Europe where they likely originated. In contrast, the only U1 representative in Europe is the G-M527 lineage whose distribution pattern is consistent with regions of Greek colonization. No clinal patterns were detected suggesting that the distributions are rather indicative of isolation by distance and demographic complexities.

Case-control association study of 65 candidate genes revealed a possible association of a SNP of HTR5A to be a factor susceptible to bipolar disease in Bulgarian population.

BACKGROUND: Bipolar affective disorder (BAD) is a psychiatric illness characterized by episodes of mania and depression. Although the etiology is not clear, epidemiological studies suggest it is a result of an interaction of genetic and environmental factors. Despite of enormous efforts and abundant studies conducted, none has yet been identified definitively a gene susceptible to bipolar disorder. METHODS: Ninety-four Bulgarian patients diagnosed with bipolar disorder and 184 Bulgarian healthy individuals, were used for genotyping of 191 single nucleotide polymorphisms (SNPs) by TaqMan and/or Invader assays. Seventeen SNPs that revealed P value less than 0.05 in the first screening were genotyped using an additional independent set of samples, consisting of 78 BAD cases and 372 controls. RESULTS: After applying the Bonferonni correction on genotyping results of 172 cases and 556 controls, only one SNP, rs1800883, in the HTR5A gene revealed a significant level of P value (P=0.000097; odds ratio=1.80 (95%CI, 1.27-2.54); corrected P=0.017). CONCLUSIONS: Our findings suggest that HTR5A gene could play an important role in the pathogenesis of bipolar disorder in our population. However these findings should be viewed with caution and replication studies in other populations are necessary in support of these findings.

Genetic structure of Europeans: a view from the North-East.

Using principal component (PC) analysis, we studied the genetic constitution of 3,112 individuals from Europe as portrayed by more than 270,000 single nucleotide polymorphisms (SNPs) genotyped with the Illumina Infinium platform. In cohorts where the sample size was >100, one hundred randomly chosen samples were used for analysis to minimize the sample size effect, resulting in a total of 1,564 samples. This analysis revealed that the genetic structure of the European population correlates closely with geography. The first two PCs highlight the genetic diversity corresponding to the northwest to southeast gradient and position the populations according to their approximate geographic origin. The resulting genetic map forms a triangular structure with a) Finland, b) the Baltic region, Poland and Western Russia, and c) Italy as its vertexes, and with d) Central- and Western Europe in its centre. Inter- and intra- population genetic differences were quantified by the inflation factor lambda (lambda) (ranging from 1.00 to 4.21), fixation index (F(st)) (ranging from 0.000 to 0.023), and by the number of markers exhibiting significant allele frequency differences in pair-wise population comparisons. The estimated lambda was used to assess the real diminishing impact to association statistics when two distinct populations are merged directly in an analysis. When the PC analysis was confined to the 1,019 Estonian individuals (0.1% of the Estonian population), a fine structure emerged that correlated with the geography of individual counties. With at least two cohorts available from several countries, genetic substructures were investigated in Czech, Finnish, German, Estonian and Italian populations. Together with previously published data, our results allow the creation of a comprehensive European genetic map that will greatly facilitate inter-population genetic studies including genome wide association studies (GWAS).