Rafał Płoski

Correlation between Genetic and Geographic Structure in Europe

Understanding the genetic structure of the European population is important, not only from a historical perspective, but also for the appropriate design and interpretation of genetic epidemiological studies. Previous population genetic analyses with autosomal markers in Europe either had a wide geographic but narrow genomic coverage [1, 2], or vice versa [3-6]. We therefore investigated Affymetrix GeneChip 500K genotype data from 2,514 individuals belonging to 23 different subpopulations, widely spread over Europe. Although we found only a low level of genetic differentiation between subpopulations, the existing differences were characterized by a strong continent-wide correlation between geographic and genetic distance. Furthermore, mean heterozygosity was larger, and mean linkage disequilibrium smaller, in southern as compared to northern Europe. Both parameters clearly showed a clinal distribution that provided evidence for a spatial continuity of genetic diversity in Europe. Our comprehensive genetic data are thus compatible with expectations based upon European population history, including the hypotheses of a south-north expansion and/or a larger effective population size in southern than in northern Europe. By including the widely used CEPH from Utah (CEU) samples into our analysis, we could show that these individuals represent northern and western Europeans reasonably well, thereby confirming their assumed regional ancestry.

GJB2 Mutations and Degree of Hearing Loss: A Multicenter Study

Hearing impairment (HI) affects 1 in 650 newborns, which makes it the most common congenital sensory impairment. Despite extraordinary genetic heterogeneity, mutations in one gene, GJB2, which encodes the connexin 26 protein and is involved in inner ear homeostasis, are found in up to 50% of patients with autosomal recessive nonsyndromic hearing loss. Because of the high frequency of GJB2 mutations, mutation analysis of this gene is widely available as a diagnostic test. In this study, we assessed the association between genotype and degree of hearing loss in persons with HI and biallelic GJB2 mutations. We performed cross-sectional analyses of GJB2 genotype and audiometric data from 1,531 persons, from 16 different countries, with autosomal recessive, mild-to-profound nonsyndromic HI. The median age of all participants was 8 years; 90% of persons were within the age range of 0-26 years. Of the 83 different mutations identified, 47 were classified as nontruncating, and 36 as truncating. A total of 153 different genotypes were found, of which 56 were homozygous truncating (T/T), 30 were homozygous nontruncating (NT/NT), and 67 were compound heterozygous truncating/nontruncating (T/NT). The degree of HI associated with biallelic truncating mutations was significantly more severe than the HI associated with biallelic nontruncating mutations (P<.0001). The HI of 48 different genotypes was less severe than that of 35delG homozygotes. Several common mutations (M34T, V37I, and L90P) were associated with mild-to-moderate HI (median 25-40 dB). Two genotypes--35delG/R143W (median 105 dB) and 35delG/dela(GJB6-D13S1830) (median 108 dB)--had significantly more-severe HI than that of 35delG homozygotes.

Mutability of Y-chromosomal microsatellites: rates, characteristics, molecular bases, and forensic implications.

Nonrecombining Y-chromosomal microsatellites (Y-STRs) are widely used to infer population histories, discover genealogical relationships, and identify males for criminal justice purposes. Although a key requirement for their application is reliable mutability knowledge, empirical data are only available for a small number of Y-STRs thus far. To rectify this, we analyzed a large number of 186 Y-STR markers in nearly 2000 DNA-confirmed father-son pairs, covering an overall number of 352,999 meiotic transfers. Following confirmation by DNA sequence analysis, the retrieved mutation data were modeled via a Bayesian approach, resulting in mutation rates from 3.78 × 10(-4) (95% credible interval [CI], 1.38 × 10(-5) - 2.02 × 10(-3)) to 7.44 × 10(-2) (95% CI, 6.51 × 10(-2) - 9.09 × 10(-2)) per marker per generation. With the 924 mutations at 120 Y-STR markers, a nonsignificant excess of repeat losses versus gains (1.16:1), as well as a strong and significant excess of single-repeat versus multirepeat changes (25.23:1), was observed. Although the total repeat number influenced Y-STR locus mutability most strongly, repeat complexity, the length in base pairs of the repeated motif, and the fathers age also contributed to Y-STR mutability. To exemplify how to practically utilize this knowledge, we analyzed the 13 most mutable Y-STRs in an independent sample set and empirically proved their suitability for distinguishing close and distantly related males. This finding is expected to revolutionize Y-chromosomal applications in forensic biology, from previous male lineage differentiation toward future male individual identification.

Lymphoid tyrosine phosphatase (PTPN22/LYP) variant and Graves' disease in a Polish population: association and gene dose-dependent correlation with age of onset.

Objective  Susceptibility to Graves’ disease (GD) is to a significant extent determined by genetic factors of which the best known are those associated with the HLA and the CTLA4 locus. Recently, two studies on British Caucasians reported that a single nucleotide polymorphism, 1858 C>T in PTPN22, encoding Arg620Trp in the lymphoid protein tyrosine phosphatase (LYP), which is a negative regulator of T‐cell activation, increases the risk of GD. The purpose of our study was to investigate whether the PTPN22‘T’ allele is associated with GD and/or its subsets, defined by clinical or genetic parameters, in a Polish population.

Comprehensive mutation analysis of 17 Y-chromosomal short tandem repeat polymorphisms included in the AmpFlSTR® Yfiler® PCR amplification kit

The Y-chromosomal short tandem repeat (Y-STR) polymorphisms included in the AmpFlSTR® Yfiler® polymerase chain reaction amplification kit have become widely used for forensic and evolutionary applications where a reliable knowledge on mutation properties is necessary for correct data interpretation. Therefore, we investigated the 17 Yfiler Y-STRs in 1,730–1,764 DNA-confirmed father–son pairs per locus and found 84 sequence-confirmed mutations among the 29,792 meiotic transfers covered. Of the 84 mutations, 83 (98.8%) were single-repeat changes and one (1.2%) was a double-repeat change (ratio, 1:0.01), as well as 43 (51.2%) were repeat gains and 41 (48.8%) repeat losses (ratio, 1:0.95). Medians from Bayesian estimation of locus-specific mutation rates ranged from 0.0003 for DYS448 to 0.0074 for DYS458, with a median rate across all 17 Y-STRs of 0.0025. The mean age (at the time of son’s birth) of fathers with mutations was with 34.40 (±11.63) years higher than that of fathers without ones at 30.32 (±10.22) years, a difference that is highly statistically significant (p < 0.001). A Poisson-based modeling revealed that the Y-STR mutation rate increased with increasing father’s age on a statistically significant level (α = 0.0294, 2.5% quantile = 0.0001). From combining our data with those previously published, considering all together 135,212 meiotic events and 331 mutations, we conclude for the Yfiler Y-STRs that (1) none had a mutation rate of >1%, 12 had mutation rates of >0.1% and four of <0.1%, (2) single-repeat changes were strongly favored over multiple-repeat ones for all loci but 1 and (3) considerable variation existed among loci in the ratio of repeat gains versus losses. Our finding of three Y-STR mutations in one father–son pair (and two pairs with two mutations each) has consequences for determining the threshold of allelic differences to conclude exclusion constellations in future applications of Y-STRs in paternity testing and pedigree analyses.

Significant genetic differentiation between Poland and Germany follows present-day political borders, as revealed by Y-chromosome analysis

To test for human population substructure and to investigate human population history we have analysed Y-chromosome diversity using seven microsatellites (Y-STRs) and ten binary markers (Y-SNPs) in samples from eight regionally distributed populations from Poland (n=913) and 11 from Germany (n=1,215). Based on data from both Y-chromosome marker systems, which we found to be highly correlated (r=0.96), and using spatial analysis of the molecular variance (SAMOVA), we revealed statistically significant support for two groups of populations: (1) all Polish populations and (2) all German populations. By means of analysis of the molecular variance (AMOVA) we observed a large and statistically significant proportion of 14% (for Y-SNPs) and 15% (for Y-STRs) of the respective total genetic variation being explained between both countries. The same population differentiation was detected using Monmonier’s algorithm, with a resulting genetic border between Poland and Germany that closely resembles the course of the political border between both countries. The observed genetic differentiation was mainly, but not exclusively, due to the frequency distribution of two Y-SNP haplogroups and their associated Y-STR haplotypes: R1a1*, most frequent in Poland, and R1*(xR1a1), most frequent in Germany. We suggest here that the pronounced population differentiation between the two geographically neighbouring countries, Poland and Germany, is the consequence of very recent events in human population history, namely the forced human resettlement of many millions of Germans and Poles during and, especially, shortly after World War II. In addition, our findings have consequences for the forensic application of Y-chromosome markers, strongly supporting the implementation of population substructure into forensic Y chromosome databases, and also for genetic association studies.

Homogeneity and distinctiveness of Polish paternal lineages revealed by Y chromosome microsatellite haplotype analysis

Abstract. Different regional populations from Poland were studied in order to assess the genetic heterogeneity within Poland, investigate the genetic relationships with other European populations and provide a population-specific reference database for anthropological and forensic studies. Nine Y-chromosomal microsatellites were analysed in a total of 919 unrelated males from six regions of Poland and in 1,273 male individuals from nine other European populations. AMOVA revealed that all of the molecular variation in the Polish dataset is due to variation within populations, and no variation was detected among populations of different regions of Poland. However, in the non-Polish European dataset 9.3% (P<0.0001) of the total variation was due to differences among populations. Consequently, differences in RST-values between all possible pairs of Polish populations were not statistically significant, whereas significant differences were observed in nearly all comparisons of Polish and non-Polish European populations. Phylogenetic analyses demonstrated tight clustering of Polish populations separated from non-Polish groups. Population clustering based on Y-STR haplotypes generally correlates well with the geography and history of the region. Thus, our data are consistent with the assumption of homogeneity of present-day paternal lineages within Poland and their distinctiveness from other parts of Europe, at least in respect to their Y-STR haplotypes. Electronic supplementary material to this paper can be obtained by using the Springer LINK server located at http://dx.doi.org/10.1007/s00439-002-0728-0.

M34T and V37I mutations in GJB2 associated hearing impairment: Evidence for pathogenicity and reduced penetrance

Despite research the role of the M34T and V37I variants of GJB2 in causing hearing impairment (HI) remains controversial. Our purpose was to test a hypothesis that M34T and V37I are pathogenic but have distinct features resulting in a reduced penetrance. We screened for known GJB2/GJB6 mutations 233 Polish consecutive unrelated subjects with non‐syndromic, sensorineural HI who were previously found to carry 35delG mutation on one chromosome. The most frequent mutations were also analyzed in ∼1,000 controls. We found that M34T and V37I were significantly (P ≪ 10−6) overrepresented among patients, but their penetrance was estimated as 1/10 relative to mutations of undisputed pathogenicity. This finding apparently could not be explained by low degree of HI associated with M34T and V37I since another mutation causing comparably mild HI (L90P) did not have reduced penetrance. Subsequent analyses showed that the patients with M34T/35delG and V37I/35delG had significantly later onset of HI than patients with other genotypes (P < 10−6) including the L90P/35delG (P = 0.006). Also, among these patients (but not others) a strong correlation between the degree of HI and its duration was found (r = 0.79, P < 10−5). We tentatively suggest that M34T and V37I might cause mild HI characterized by relatively late onset and progression.

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.

Examination of DNA methylation status of the ELOVL2 marker may be useful for human age prediction in forensic science

Age estimation in forensic investigations may complement the prediction of externally visible characteristics and the inference of biogeographical ancestry, thus allowing a better description of an unknown individual. Multiple CpG sites that show linear correlation between age and degree of DNA methylation have been identified in the human genome, providing a selection of candidates for age prediction. In this study, we optimized an assay based on bisulfite conversion and pyrosequencing of 7 CpG sites located in the ELOVL2 gene. Examination of 303 blood samples collected from individuals aged 2-75 years allowed selection of the most informative site, explaining 83% of variation in age. The final linear regression model included two CpG sites in ELOVL2 and enabled age prediction with R(2)=0.859, prediction error=6.85 and mean absolute deviation MAD=5.03. Examination of a testing set of 124 blood samples (MAD=5.75) showed that 68.5% of samples were correctly predicted, assuming that chronological and predicted ages matched ± 7 years. It was found that the ELOVL2 methylation status in bloodstains had not changed significantly after 4 weeks of storage in room temperature conditions. Analysis of 45 bloodstains deposited on tissue paper after 5, 10 and 15 years of storage in room conditions indicated that although a gradual decrease of positive PCR results was observed, the general age prediction success rate remained similar and equaled 60-78%. The obtained results show that the ELOVL2 locus provides a very good source of information about human chronological age based on analysis of blood, including bloodstains, and it may constitute a powerful and reliable predictor in future forensic age estimation models.