Y-chromosome variation in a Norwegian population sample

Microsatellites' mutation modeling through the analysis of the Y-chromosomal transmission: Results of a GHEP-ISFG collaborative study

The Spanish and Portuguese Speaking Working Group of the International Society for Forensic Genetics (GHEP-ISFG) organized a collaborative study on mutations of Y-chromosomal short tandem repeats (Y-STRs). New data from 2225 father-son duos and data from 44 previously published reports, corresponding to 25,729 duos, were collected and analyzed. Marker-specific mutation rates were estimated for 33 Y-STRs. Although highly dependent on the analyzed marker, mutations compatible with the gain or loss of a single repeat were 23.2 times more likely than those involving a greater number of repeats. Longer alleles (relatively to the modal one) showed to be nearly twice more mutable than the shorter ones. Within the subset of longer alleles, the loss of repeats showed to be nearly twice more likely than the gain. Conversely, shorter alleles showed a symmetrical trend, with repeat gains being twofold more frequent than reductions. A positive correlation between the paternal age and the mutation rate was observed, strengthening previous findings. The results of a machine learning approach, via logistic regression analyses, allowed the establishment of algebraic formulas for estimating the probability of mutation depending on paternal age and allele length for DYS389I, DYS393 and DYS627. Algebraic formulas could also be established considering only the allele length as predictor for DYS19, DYS389I, DYS389II-I, DYS390, DYS391, DYS393, DYS437, DYS439, DYS449, DYS456, DYS458, DYS460, DYS481, DYS518, DYS533, DYS576, DYS626 and DYS627 loci. For the remaining Y-STRs, a lack of statistical significance was observed, probably as a consequence of the small effective size of the subsets available, a common difficulty in the modeling of rare events as is the case of mutations. The amount of data used in the different analyses varied widely, depending on how the data were reported in the publications analyzed. This shows a regrettable waste of produced data, due to inadequate communication of the results, supporting an urgent need of publication guidelines for mutation studies.

The sequence of the repetitive motif influences the frequency of multistep mutations in Short Tandem Repeats

Microsatellites, or Short Tandem Repeats (STRs), are subject to frequent length mutations that involve the loss or gain of an integer number of repeats. This work aimed to investigate the correlation between STRs' specific repetitive motif composition and mutational dynamics, specifically the occurrence of single- or multistep mutations. Allelic transmission data, comprising 323,818 allele transfers and 1,297 mutations, were gathered for 35 Y-chromosomal STRs with simple structure. Six structure groups were established: ATT, CTT, TCTA/GATA, GAAA/CTTT, CTTTT, and AGAGAT, according to the repetitive motif present in the DNA leading strand of the markers. Results show that the occurrence of multistep mutations varies significantly among groups of markers defined by the repetitive motif. The group of markers with the highest frequency of multistep mutations was the one with repetitive motif CTTTT (25% of the detected mutations) and the lowest frequency corresponding to the group with repetitive motifs TCTA/GATA (0.93%). Statistically significant differences (α = 0.05) were found between groups with repetitive motifs with different lengths, as is the case of TCTA/GATA and ATT (p = 0.0168), CTT (p < 0.0001) and CTTTT (p < 0.0001), as well as between GAAA/CTTT and CTTTT (p = 0.0102). The same occurred between the two tetrameric groups GAAA/CTTT and TCTA/GATA (p < 0.0001) - the first showing 5.7 times more multistep mutations than the second. When considering the number of repeats of the mutated paternal alleles, statistically significant differences were found for alleles with 10 or 12 repeats, between GATA and ATT structure groups. These results, which demonstrate the heterogeneity of mutational dynamics across repeat motifs, have implications in the fields of population genetics, epidemiology, or phylogeography, and whenever STR mutation models are used in evolutionary studies in general.

Joint Genetic Analyses of Mitochondrial and Y-Chromosome Molecular Markers for a Population from Northwest China

The genetic markers on mitochondria DNA (mtDNA) and Y-chromosome can be applied as a powerful tool in population genetics. We present a study to reveal the genetic background of Kyrgyz group, a Chinese ethnic group living in northwest China, and genetic polymorphisms of 60 loci on maternal inherited mtDNA and 24 loci on paternal inherited Y-chromosome short tandem repeats (Y-STRs) were investigated. The relationship between the two systems was tested, and the result indicated that they were statistically independent from each other. The genetic distances between Kyrgyz group and 11 reference populations for mtDNA, and 13 reference populations for Y-STRs were also calculated, respectively. The present results demonstrated that the Kyrgyz group was genetically closer to East Asian populations than European populations based on the mtDNA loci but the other way around for the Y-STRs. The genetic analyses could largely strengthen the understanding for the genetic background of the Kyrgyz group.

Genomic and Strontium Isotope Variation Reveal Immigration Patterns in a Viking Age Town

The impact of human mobility on the northern European urban populations during the Viking and Early Middle Ages and its repercussions in Scandinavia itself are still largely unexplored. Our study of the demographics in the final phase of the Viking era is the first comprehensive multidisciplinary investigation that includes genetics, isotopes, archaeology, and osteology on a larger scale. This early Christian dataset is particularly important as the earlier common pagan burial tradition during the Iron Age was cremation, hindering large-scale DNA analyses. We present genome-wide sequence data from 23 individuals from the 10^th to 12^th century Swedish town of Sigtuna. The data revealed high genetic diversity among the early urban residents. The observed variation exceeds the genetic diversity in distinct modern-day and Iron Age groups of central and northern Europe. Strontium isotope data suggest mixed local and non-local origin of the townspeople. Our results uncover the social system underlying the urbanization process of the Viking World of which mobility was an intricate part and was comparable between males and females. The inhabitants of Sigtuna were heterogeneous in their genetic affinities, probably reflecting both close and distant connections through an established network, confirming that early urbanization processes in northern Europe were driven by migration.

Empirical evaluation reveals best fit of a logistic mutation model for human Y-chromosomal microsatellites

The rate of microsatellite mutation is dependent upon both the allele length and the repeat motif, but the exact nature of this relationship is still unknown. We analyzed data on the inheritance of human Y-chromosomal microsatellites in father-son duos, taken from 24 published reports and comprising 15,285 directly observable meioses. At the six microsatellites analyzed (DYS19, DYS389I, DYS390, DYS391, DYS392, and DYS393), a total of 162 mutations were observed. For each locus, we employed a maximum-likelihood approach to evaluate one of several single-step mutation models on the basis of the data. For five of the six loci considered, a novel logistic mutation model was found to provide the best fit according to Akaike's information criterion. This implies that the mutation probability at the loci increases (nonlinearly) with allele length at a rate that differs between upward and downward mutations. For DYS392, the best fit was provided by a linear model in which upward and downward mutation probabilities increase equally with allele length. This is the first study to empirically compare different microsatellite mutation models in a locus-specific fashion.

PCR-free digital minisatellite tandem repeat genotyping

We demonstrated a proof-of-concept for novel minisatellite tandem repeat typing, called PCR-free digital VNTR (variable number tandem repeat) typing, which is composed of three steps: a ligation reaction instead of PCR thermal cycling, magnetic bead-based solid-phase capture for purification, and an elongated sample stacking microcapillary electrophoresis (μCE) for sensitive digital coding of repeat number. We designed a 16-bp fluorescently labeled ligation probe which is complementary to a repeat unit of a biotinylated synthetic template mimicking the human D1S80 VNTR locus and is randomly hybridized with the minisatellite tandem repeats. A quick isothermal ligation reaction was followed to link the adjacent ligation probes on the DNA templates, and then the ligated products were purified by streptavidin-coated magnetic beads. After a denaturing step, a large amount of ligated products whose size difference was equivalent to the repeat unit were released and recovered. Through the elongated sample stacking μCE separation on a microdevice, the fluorescence signal of the ligated products was generated in the electropherogram and the peak number was directly counted which was exactly matched with the repeat number of VNTR locus. We could successfully identify the minisatellite tandem repeat number with only 5 fmol of DNA template in 30 min.

Investigations of the Y chromosome, male founder structure and YSTR mutation rates in the Old Order Amish

OBJECTIVES

Using Y chromosome short tandem repeat (YSTR) genotypes, (1) evaluate the accuracy and completeness of the Lancaster County Old Order Amish (OOA) genealogical records and (2) estimate YSTR mutation rates.

METHODS

Nine YSTR markers were genotyped in 739 Old Order Amish males who participated in several ongoing genetic studies of complex traits and could be connected into one of 28 all-male lineage pedigrees constructed using the Anabaptist Genealogy Database and the query software Ped-Hunter. A putative founder YSTR haplotype was constructed for each pedigree, and observed and inferred father-son transmissions were used to estimate YSTR mutation rates.

RESULTS

We inferred 27 distinct founder Y chromosome haplotypes in the 28 male lineages, which encompassed 27 surnames accounting for 98% of Lancaster OOA households. Nearly all deviations from founder haplotypes were consistent with mutation events rather than errors. The estimated marker-specific mutation rates ranged from 0 to 1.09% (average 0.33% using up to 283 observed meioses only and 0.28% using up to 1,232 observed and inferred meioses combined).

CONCLUSIONS

These data confirm the accuracy and completeness of the male lineage portion of the Anabaptist Genealogy Database and contribute mutation rate estimates for several commonly used Y chromosome STR markers.

Regional differences among the Finns: A Y-chromosomal perspective

Twenty-two Y-chromosomal markers, consisting of fourteen biallelic markers (YAP/DYS287, M170, M253, P37, M223, 12f2, M9, P43, Tat, 92R7, P36, SRY-1532, M17, P25) and eight STRs (DYS19, DYS385a/b, DYS388, DYS389I/II, DYS390, DYS391, DYS392, DYS393), were analyzed in 536 unrelated Finnish males from eastern and western subpopulations of Finland. The aim of the study was to analyze regional differences in genetic variation within the country, and to analyze the population history of the Finns. Our results gave further support to the existence of a sharp genetic border between eastern and western Finns so far observed exclusively in Y-chromosomal variation. Both biallelic haplogroup and STR haplotype networks showed bifurcated structures, and similar clustering was evident in haplogroup and haplotype frequencies and genetic distances. These results suggest that the western and eastern parts of the country have been subject to partly different population histories, which is also supported by earlier archaeological, historical and genetic data. It seems probable that early migrations from Finno-Ugric sources affected the whole country, whereas subsequent migrations from Scandinavia had an impact mainly on the western parts of the country. The contacts between Finland and neighboring Finno-Ugric, Scandinavian and Baltic regions are evident. However, there is no support for recent migrations from Siberia and Central Europe. Our results emphasize the importance of incorporating Y-chromosomal data to reveal the population substructure which is often left undetected in mitochondrial DNA variation. Early assumptions of the homogeneity of the isolated Finnish population have now proven to be false, which may also have implications for future association studies.

A study of mutation rates and the characterisation of intermediate, null and duplicated alleles for 13 Y chromosome STRs

Previously reported Y chromosome STR haplotype databases for three UK population groups, plus additionally analysed samples, have been scrutinised for the presence of non-standard (intermediate, null and duplicated) alleles. These alleles have been characterised by sequencing, some showing changes in the repeat structure, and the frequencies reported. Mutation rates for each of the 13 STRs have been calculated when analysis of father-son pairs has been possible. An example illustrating the use of non-standard alleles in a large family tree is outlined.

Y-chromosomal STR haplotypes in a Northeast Italian population sample using 17plex loci PCR assay

One hundred fifty-five unrelated, autochthonous healthy males from Northeast Italy were typed for the 17 Y-chromosome short tandem repeat (STR) (Y-STR) loci DYS456, DYS389I, DYS390, DYS389II, DYS458, DYS19, DYS385, DYS393, DYS391, DYS439, DYS635, DYS392, Y GATA H4, DYS437, DYS438, DYS448 using the AmpFLSTR Yfiler polymerase chain reaction amplification kit. A total of 153 different haplotypes were observed, and among these, 151 were unique, while 2 were found two times. The overall haplotype diversity was 0.9997. Furthermore, 50 father-son pairs, previously confirmed by autosomal STR analysis, were typed using the same set of 17 Y-STR loci, and, among 850 allele transfers, three mutation events were identified, giving an average mutation rate of 3.53 x 10(-3) per locus per generation (95% confidence interval 0.73-1.03).

Y-STR haplotypes in Central Anatolia region of Turkey

Allele frequencies and haplotypes of the 11 Y-chromosome STRs loci, namely DYS19, DYS385a/b, DYS389I/II, DYS390, DYS391, DYS392, DYS393, DYS438, and DYS439 were determined in a sample of 113 unrelated males from the Central Anatolia region of Turkey. In the 113 samples 106 different haplotypes were encountered, of which 100 were observed only once. The overall haplotype diversity was 0.9987. In the study, a duplication at locus DYS19 and locus DYS393 was observed. The results demonstrate that these loci will be very useful for human identification in forensic cases and paternity tests in the Central Anatolia region.

Y-chromosomal microsatellite mutation rates: differences in mutation rate between and within loci

Precise estimates of mutation rates at Y-chromosomal microsatellite STR (short tandem repeat) loci make an important basis for paternity diagnostics and dating of Y chromosome lineage origins. There are indications of considerable locus mutation rate variability between (inter-) and within (intra-) loci. We have studied nine Y-STR loci-DYS19, DYS389I/II, DYS390, DYS391, DYS392, DYS393, DYS385, and DYS388-in 1,766 father-son pairs of confirmed paternity (a total of 15,894 meioses). Five biallelic markers were also analyzed in the fathers-Tat, YAP, 12f2, SRY1532, and 92R7-defining haplogroups 1, 2, 3, 4, 9, and 16, respectively. A total of 36 fragment length mutations were observed: 24 gains (22 single-step, two double-step) and 12 single-step losses. Thus, there was a significant surplus of gains (p=0.045). Overall, the mutation rate was positively correlated to STR repeat length and there was a significant relative excess of losses in long alleles and gains in short alleles (p=0.043). In contrast to the situation in autosomal STR loci and in MSY-1, no noteworthy correlation between mutation rate and the father's age at the child's birth was observed. We observed significant interlocus differences in Y-STR mutation rates (p<0.01). The number of observed mutations ranged from zero in DYS392 to eight in DYS391 and DYS390. We have also demonstrated obvious differences in mutation rates between the haplogroups studied (p=0.024), a phenomenon that is a reflection of the dependence of mutation rate on allele size. Our study has thus demonstrated the necessity of not only locus-specific, but even allele-specific, mutation rate estimates for forensic and population genetic purposes, and provides a considerable basis for such estimates.

Genomics of the human Y-chromosome. 1. Association with male infertility

The human Y chromosome contains over 60 million nucleotides, but least number of genes compared to any other chromosome and acts as a genetic determinant of the male characteristic features. The male specific region, MSY, comprising 95% of the Y chromosome represents a mosaic of heterochromatic and three classes of euchromatic (X-transposed, X-degenerate and ampliconic) sequences. Thus far, 156 transcription units, 78 protein-coding genes and 27 distinct proteins of the Y chromosome have been identified. The MSY euchromatic sequences show frequent gene conversion. Of the eight massive palindromes identified on the human Y chromosome, six harbor vital testis specific genes. The human male infertility has been attributed to mutations in the genes on Y chromosome and autosomes and failures of several physical and physiological attributes including paracrine controls. In addition, deletion of any one or all the three azoospermia (AZFa, AZFb or AZFc) factor(s) and some still unidentified regulatory elements located elsewhere in the genome result in infertility. Characterization of palindromic complexes on the long arm of Y chromosome encompassing AZFb and AZFc regions and identification of HERV15 class of endogenous retroviruses close to AZFa region have facilitated our understanding on the organization of azoospermia factors. Considerable overlap of the AZFb and AZFc regions encompassing a number of genes and transcripts has been shown to exist. However, barring details on AZF, information on the exact number of genes or the types of mutations prevalent in the infertile male is not available. Similarly, roles of sizable body of repetitive DNA present in close association with transcribing sequences on the Y chromosome are yet not clear. In a clinical setting with known cases of infertility, systematic search for loss or gain of these repeat elements would help understand their biological role(s). We present a brief overview on the genetic complexity of the human Y chromosome in the context of human male infertility.

Haploid allele mapping of Y-chromosome minisatellite, MSY1 (DYF155S1), to a Japanese population

The present study analyses the human Y-chromosome minisatellite locus, MSY1 (DYF155S1), in 205 Japanese males of 191 pedigrees using the minisatellite variant repeat (MVR) mapping system. The internal haploid structures of the detected alleles considerably varied and consisted of three major repeat units: types 2, 3 and 4. A comparison of the haploid profiles of the MVR codes identified 185 distinct alleles, of which only five were shared. We did not detect a type 1 repeat unit, and variations were frequent at the 5' end of the minisatellite locus. Within an analysis of 24 paternally linked DNA samples donated by ten families, no mutational events were identified even over two generation gaps. Furthermore, we applied this mapping system to a paternity test in which the alleged father was missing.

High level of male-biased Scandinavian admixture in Greenlandic Inuit shown by Y-chromosomal analysis

We have used binary markers and microsatellites on the Y chromosome to analyse diversity in a sample of Greenlandic Inuit males. This sample contains Y chromosomes typical of those found in European populations. Because the Y chromosome has a unique and robust phylogeny of a time depth that precedes the split between European and Native American populations, it is possible to assign chromosomes in an admixed population to either continental source. On this basis, 58+/-6% of these Y chromosomes have been assigned to a European origin. The high proportion of European Y chromosomes contrasts with a complete absence of European mitochondrial DNA and indicates strongly male-biased European admixture into Inuit. Comparison of the European component of Inuit Y chromosomes with European population data suggests that they have their origins in Scandinavia. There are two potential source populations: Norse settlers from Iceland, who may have been assimilated 500 years ago, and the Danish-Norwegian colonists of the eighteenth century. Insufficient differentiation between modern Icelandic and Danish Y chromosomes means that a choice between these cannot be made on the basis of diversity analysis. However, the extreme sex bias in the admixture makes the later event more likely as the source.

Different genetic components in the Norwegian population revealed by the analysis of mtDNA and Y chromosome polymorphisms

The genetic composition of the Norwegian population was investigated by analysing polymorphisms associated with both the mitochondrial DNA (mtDNA) and Y chromosome loci in a sample of 74 Norwegian males. The combination of their uniparental mode of inheritance and the absence of recombination make these haplotypic stretches of DNA the tools of choice in evaluating the different components of a population's gene pool. The sequencing of the Dloop and two diagnostic RFLPs (AluI 7025 and HinfI at 12 308) allowed us to classify the mtDNA molecules in 10 previously described groups. As for the Y chromosome the combination of binary markers and microsatellites allowed us to compare our results to those obtained elsewhere in Europe. Both mtDNA and Y chromosome polymorphisms showed a noticeable genetic affinity between Norwegians and central Europeans, especially Germans. When the phylogeographic analysis of the Y chromosome haplotypes was attempted some interesting clues on the peopling of Norway emerged. Although Y chromosome binary and microsatellite data indicate that 80% of the haplotypes are closely related to Central and western Europeans, the remainder share a unique binary marker (M17) common in eastern Europeans with informative microsatellite haplotypes suggesting a different demographic history. Other minor genetic influences on the Norwegian population from Uralic speakers and Mediterranean populations were also highlighted.

Mutation at minisatellite locus DYF155S1: allele length mutation rate is affected by age of progenitor

A father/son material consisting of 1071 pairs was screened for de novo allele length mutation in locus DYF155S1. Six hundred of these pairs were also analyzed in locus DYF155S1 to detect de novo mutations in the minisatellite variant repeat (MVR)-code not resulting in a length change ("boundary switch" mutations). A modified MVR-polymerase chain reaction (PCR) method was used for this purpose. Twenty-seven de novo allele length mutations and eight "boundary switch" mutations were detected indicating mutation frequencies of approximately 2.5% and 1.3%, respectively. The combined mutation rate for MVR-code mutation is approximately 3.8%. There is a significant increase in mutation rate with paternal age (p = 0.049) in allele length mutations. In the present material, the mutation rate in the oldest age group is three times that of the youngest age group. A similar age relationship is not observed in "boundary switch" mutations. A comparison between progenitors and the other fathers in the material revealed no obvious association between mutation rate and allele length or modular structure (variation in repeat sequence). More than 75% of the length mutations involved the gain or loss of one repeat only. This finding as well as the observed paternal age influence on mutation rate, suggests replication slippage to be the major mutation mechanism in length mutations. However, in one particular case, an allele length mutant revealed rearrangements with direct duplication of repeats at distant sites within the repeat array, and with both loss and gain of repeats. Such complex structural changes could indicate that some of the mutants might arise from sister chromatide exchange. The mutation rate of "boundary switch" mutations is by far higher than would be expected if these mutations are two independent one-step allele length mutations. A different age distribution of "boundary switch" mutations than of allele length mutations also argue against such a hypothesis. Together this could indicate that "boundary switches" are products of another mutation mechanism than the one-step allele length mutations.

DXYS267: DYS393 and its X chromosome counterpart

The GATA repeat DYS393 was reported in 1987 among other Y-specific short tandem repeats. It has since been used for forensic and evolutionary studies. We decided to test its Y-specificity when we found that female DNA gave amplicons, in agreement with recent GDB-recorded experiences on radiation hybrids. Parent-child triplets revealed that heterozygous daughters always carried the same paternally derived amplicon which, however, was not amplified in their fathers' DNAs. The X-assignment was verified in larger families. A half-new primer set with a new reverse DYS393 primer, outside the old one, resulted in X amplicons in females as well as Y and X amplicons in males. This new primer set defines the new DXYS267 (GDB Data Curation). DNA-sequencing revealed four base pair differences between the Y- and the X-sequences. Two are within the reverse primer site sequence, thus probably causing preferential hybridization to the Y sequence when using the conventional primers. The two others are within the repeat array, giving the regular repeat GATA in the Y-sequence, and TATA and GACA, respectively, in the X-sequence. Allele frequency distribution in DYS393 was studied in 300 unrelated Norwegian males, allele distribution in the X-locus in 48 Norwegian women. Even if allele repeat numbers are overlapping between the loci, leading to identical fragment lengths, the allele distribution is different between DYS393 and the X-chromosome locus. The differences between the two homologous loci on the Y and X indicate a considerable lap of time since common ancestry. To avoid co-amplification of the X-locus in DYS393 typing, primer A was elongated to include one of the sequence differences between the two loci. This to a considerable extent improved the specificity of the DYS393 primers.