Vaidutis Kučinskas

Y-Chromosomal Diversity in Europe Is Clinal and Influenced Primarily by Geography, Rather than by Language

Clinal patterns of autosomal genetic diversity within Europe have been interpreted in previous studies in terms of a Neolithic demic diffusion model for the spread of agriculture; in contrast, studies using mtDNA have traced many founding lineages to the Paleolithic and have not shown strongly clinal variation. We have used 11 human Y-chromosomal biallelic polymorphisms, defining 10 haplogroups, to analyze a sample of 3,616 Y chromosomes belonging to 47 European and circum-European populations. Patterns of geographic differentiation are highly nonrandom, and, when they are assessed using spatial autocorrelation analysis, they show significant clines for five of six haplogroups analyzed. Clines for two haplogroups, representing 45% of the chromosomes, are continentwide and consistent with the demic diffusion hypothesis. Clines for three other haplogroups each have different foci and are more regionally restricted and are likely to reflect distinct population movements, including one from north of the Black Sea. Principal-components analysis suggests that populations are related primarily on the basis of geography, rather than on the basis of linguistic affinity. This is confirmed in Mantel tests, which show a strong and highly significant partial correlation between genetics and geography but a low, nonsignificant partial correlation between genetics and language. Genetic-barrier analysis also indicates the primacy of geography in the shaping of patterns of variation. These patterns retain a strong signal of expansion from the Near East but also suggest that the demographic history of Europe has been complex and influenced by other major population movements, as well as by linguistic and geographic heterogeneities and the effects of drift.

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 (λ) (ranging from 1.00 to 4.21), fixation index (Fst) (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).

Mutation history of the Roma/Gypsies

The 8-10 million European Roma/Gypsies are a founder population of common origins that has subsequently split into multiple socially divergent and geographically dispersed Gypsy groups. Unlike other founder populations, whose genealogy has been extensively documented, the demographic history of the Gypsies is not fully understood and, given the lack of written records, has to be inferred from current genetic data. In this study, we have used five disease loci harboring private Gypsy mutations to examine some missing historical parameters and current structure. We analyzed the frequency distribution of the five mutations in 832-1,363 unrelated controls, representing 14 Gypsy populations, and the diversification of chromosomal haplotypes in 501 members of affected families. Sharing of mutations and high carrier rates supported a strong founder effect, and the identity of the congenital myasthenia 1267delG mutation in Gypsy and Indian/Pakistani chromosomes provided the best evidence yet of the Indian origins of the Gypsies. However, dramatic differences in mutation frequencies and haplotype divergence and very limited haplotype sharing pointed to strong internal differentiation and characterized the Gypsies as a founder population comprising multiple subisolates. Using disease haplotype coalescence times at the different loci, we estimated that the entire Gypsy population was founded approximately 32-40 generations ago, with secondary and tertiary founder events occurring approximately 16-25 generations ago. The existence of multiple subisolates, with endogamy maintained to the present day, suggests a general approach to complex disorders in which initial gene mapping could be performed in large families from a single Gypsy group, whereas fine mapping would rely on the informed sampling of the divergent subisolates and searching for the shared genomic region that displays the strongest linkage disequilibrium with the disease.

Duplications in addition to terminal deletions are present in a proportion of ring chromosomes. Clues to the mechanisms of formation

Background and methods: Ring chromosomes are often associated with abnormal phenotypes because of loss of genomic material at one or both ends. In some cases no deletion has been detected and the abnormal phenotype has been attributed to mitotic ring instability. We investigated 33 different ring chromosomes in patients with phenotypic abnormalities by array based comparative genomic hybridisation (CGH) and fluorescence in situ hybridisation (FISH). Results: In seven cases we found not only the expected terminal deletion but also a contiguous duplication. FISH analysis in some of these cases demonstrated that the duplication was inverted. Thus these ring chromosomes derived through a classical inv dup del rearrangement consisting of a deletion and an inverted duplication. Discussion: Inv dup del rearrangements have been reported for several chromosomes, but hardly ever in ring chromosomes. Our findings highlight a new mechanism for the formation of some ring chromosomes and show that inv dup del rearrangements may be stabilised not only through telomere healing and telomere capture but also through circularisation. This type of mechanism must be kept in mind when evaluating possible genotype–phenotype correlations in ring chromosomes since in these cases: (1) the deletion may be larger or smaller than first estimated based on the size of the ring, with a different impact on the phenotype; and (2) the associated duplication will in general cause further phenotypic anomalies and might confuse the genotype–phenotype correlation. Moreover, these findings explain some phenotypic peculiarities which previously were attributed to a wide phenotypic variation or hidden mosaicism related to the instability of the ring.

Reconstructing the population history of European Romani from genome-wide data

The Romani, the largest European minority group with approximately 11 million people, constitute a mosaic of languages, religions, and lifestyles while sharing a distinct social heritage. Linguistic and genetic studies have located the Romani origins in the Indian subcontinent. However, a genome-wide perspective on Romani origins and population substructure, as well as a detailed reconstruction of their demographic history, has yet to be provided. Our analyses based on genome-wide data from 13 Romani groups collected across Europe suggest that the Romani diaspora constitutes a single initial founder population that originated in north/northwestern India ~1.5 thousand years ago (kya). Our results further indicate that after a rapid migration with moderate gene flow from the Near or Middle East, the European spread of the Romani people was via the Balkans starting ~0.9 kya. The strong population substructure and high levels of homozygosity we found in the European Romani are in line with genetic isolation as well as differential gene flow in time and space with non-Romani Europeans. Overall, our genome-wide study sheds new light on the origins and demographic history of European Romani.

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 >145 000 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.

Y Chromosome and Mitochondrial DNA Variation in Lithuanians

The genetic composition of the Lithuanian population was investigated by analysing mitochondrial DNA hypervariable region 1, RFLP polymorphisms and Y chromosomal biallelic and STR markers in six ethnolinguistic groups of Lithuanians, to address questions about the origin and genetic structure of the present day population. There were no significant genetic differences among ethnolinguistic groups, and an analysis of molecular variance confirmed the homogeneity of the Lithuanian population. MtDNA diversity revealed that Lithuanians are close to both Slavic (Indo‐European) and Finno‐Ugric speaking populations of Northern and Eastern Europe. Y‐chromosome SNP haplogroup analysis showed Lithuanians to be closest to Latvians and Estonians. Significant differences between Lithuanian and Estonian Y chromosome STR haplotypes suggested that these populations have had different demographic histories. We suggest that the observed pattern of Y chromosome diversity in Lithuanians may be explained by a population bottleneck associated with Indo‐European contact. Different Y chromosome STR distributions in Lithuanians and Estonians might be explained by different origins or, alternatively, be the result of some period of isolation and genetic drift after the population split.

Replication of Novel Susceptibility Locus for Nonsyndromic Cleft Lip With or Without Cleft Palate on Chromosome 8q24 in Estonian and Lithuanian Patients

Cleft lip with or without cleft palate (CL/P) is one of the mostcommonhumanbirthdefects.CL/Pmayoccuraspartofacomplexmalformation syndrome or, as is more commonly the case, as anonsyndromic malformation. In European populations non-syndromic cleft lip with or without cleft palate (NSCL/P) has abirth prevalence of 1/700–1/1,000. NSCL/P has a multifactorialetiology resulting from interactions between multiple genetic andenvironmentalfactors.Nonsyndromiccleftpalateonly(NSCPO)isthe second most common form of nonsyndromic clefting, with aprevalence of 1/2,000 in most European populations. AlthoughNSCL/P and NSCPO are generally considered to be develop-mentally and genetically distinct entities, there is evidence thatthey may have an etiological overlap [Kondo et al., 2002; Sivertsenet al., 2008].The identification of susceptibility genes for NSCL/P has beenthe subject of extensive research. Of the many candidate genesinvestigated, however, only the IRF6 gene has shown a convincingdegreeofconsistencyacrossstudies[Zuccheroetal.,2004;Rahimovet al., 2008]. The first genome-wide association scan (GWAS)recently reported a new major susceptibility locus at chromosome8q24.21,with the strongestassociationsignal having been observedfor rs987525 (P¼3.34 10

Genetic Heritage of the Balto-Slavic Speaking Populations: A Synthesis of Autosomal, Mitochondrial and Y-Chromosomal Data

The Slavic branch of the Balto-Slavic sub-family of Indo-European languages underwent rapid divergence as a result of the spatial expansion of its speakers from Central-East Europe, in early medieval times. This expansion–mainly to East Europe and the northern Balkans–resulted in the incorporation of genetic components from numerous autochthonous populations into the Slavic gene pools. Here, we characterize genetic variation in all extant ethnic groups speaking Balto-Slavic languages by analyzing mitochondrial DNA (n = 6,876), Y-chromosomes (n = 6,079) and genome-wide SNP profiles (n = 296), within the context of other European populations. We also reassess the phylogeny of Slavic languages within the Balto-Slavic branch of Indo-European. We find that genetic distances among Balto-Slavic populations, based on autosomal and Y-chromosomal loci, show a high correlation (0.9) both with each other and with geography, but a slightly lower correlation (0.7) with mitochondrial DNA and linguistic affiliation. The data suggest that genetic diversity of the present-day Slavs was predominantly shaped in situ, and we detect two different substrata: ‘central-east European’ for West and East Slavs, and ‘south-east European’ for South Slavs. A pattern of distribution of segments identical by descent between groups of East-West and South Slavs suggests shared ancestry or a modest gene flow between those two groups, which might derive from the historic spread of Slavic people.

Variation in FGF1, FOXE1, and TIMP2genes is associated with nonsyndromic cleft lip with or without cleft palate

BACKGROUND Nonsyndromic cleft lip with or without cleft palate (CL/P) is a common complex birth defect caused by the interaction between multiple genes and environmental factors. METHODS Five hundred and eighty-seven single nucleotide polymorphisms in 40 candidate genes related to orofacial clefting were tested for association with CL/P in a clefting sample composed of 300 patients and 606 controls from Estonian, Latvian, and Lithuanian populations. RESULTS In case-control comparisons, the minor alleles of FGF1 rs34010 (p = 4.56 × 10(-4) ), WNT9B rs4968282 (p = 0.0013), and FOXE1 rs7860144 (p = 0.0021) were associated with a decreased risk of CL/P. Multiple haplotypes in FGF1, FOXE1, and TIMP2 and haplotypes in WNT9B, PVRL2, and LHX8 were associated with CL/P. The strongest association was found for protective haplotype rs250092/rs34010 GT in the FGF1 gene (p = 5.01 × 10(-4) ). The strongest epistatic interaction was observed between the COL2A1 and WNT3 genes. CONCLUSIONS Our results provide for the first time evidence implicating FGF1 in the occurrence of CL/P, and support TIMP2 and WNT9B as novel loci predisposing to CL/P. We have also replicated recently reported significant associations between variants in or near FOXE1 and CL/P. It is likely that variation in FOXE1, TIMP2, and the FGF and Wnt signaling pathway genes confers susceptibility to nonsyndromic CL/P in Northeastern European populations.