Characteristics and frequency of germline mutations at microsatellite loci from the human Y chromosome, as revealed by direct observation in father/son pairs

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.

Improving the strategy to identify historical military remains: a literature review and Y-STR meta-analysis

The identification of historical military remains by Unrecovered War Casualties-Army (UWC-A) currently relies on Y-chromosome Short Tandem Repeat (Y-STR) testing when maternal relatives are not available, or when a mitochondrial DNA match does not provide sufficient certainty of identification. However, common Y-STR profiles (using Yfiler™) between sets of remains or families often prevent identification. To resolve these cases, an investigation of additional Y-DNA markers is needed for their potential inclusion into the DNA identification strategy. The number of genetic transmissions between missing soldiers and their living relatives needs to be considered to avoid false exclusions between paternal relatives. Analysis of 236 World War I/II (WWI/II) era pairs of relatives identified up to seven genetic transmissions between WWII soldiers and their living relatives, and nine for WWI. Previous Y-STR meta-analyses were published approximately 10 years ago when rapidly mutating markers were relatively new. This paper reports a contemporary literature review and meta-analysis of 35 studies (which includes 23 studies not previously used in meta-analysis) and 23 commonly used Y-STR's mutation rates to inform the inclusion of additional loci to UWC-A's DNA identification strategy. Meta-analysis found mutation data for a given Y-STR locus could be pooled between studies and that the mutation rates were significantly different between some loci (at P < 0.05). Based on this meta-analysis, we have identified two additional markers from PowerPlex® Y23 for potential inclusion in UWC-A's identification strategy. Further avenues for potential experimental exploration are discussed.

Key points

From 236 UWC-A pairs of relatives, we observed up to nine genetic transmissions between WWI soldiers and their living relatives, and seven for WWII.MedCalc® software for meta-analysis utilizing the Freeman-Tukey transformation was run, which analysed 35 published studies and 23 commonly used loci. Previous Y-STR mutation rate meta-analyses are now 10 years old; this paper includes 23 studies that were not included in previous meta-analyses.Through meta-analysis, we identify two markers from PowerPlex® Y23 for potential inclusion in UWC-A's historical remains identification strategy (alongside Yfiler™). We discuss potential next steps for experimental exploration of additional Y-DNA markers.

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.

Analysis of mutation rates and haplotypes of 23 Y-chromosomal STRs in Korean father-son pairs

Y-chromosomal short tandem repeats (Y-STRs) have been widely used in forensic genetics, and accurate knowledge of mutation rates at Y-STR loci is essential in kinship analysis. The main aim of this study was to estimate Y-STR mutation rates in Korean males. To obtain locus-specific mutations and haplotypes at 23 Y-STRs, we analyzed samples from 620 Korean father-son pairs. In addition, we also analyzed 476 unrelated individuals using the PowerPlex® Y23 System, with the aim of augmenting the available data for the Korean population. The PowerPlex® Y23 system facilitates analysis of the 23 Y-STR loci (DYS576, DYS570, DYS458, DYS635, DYS389 II, DYS549, DYS385, DYS481, DYS439, DYS456, DYS389 I, DYS19, DYS393, DYS391, DYS533, DYS437, DYS390, Y GATA H4, DYS448, DYS438, DYS392, and DYS643). Locus-specific mutation rate estimates varied from 0.00 to 8.06 × 10^-3 per generation, with an average mutation rate of 2.17 × 10^-3 (95% CI, 1.5-3.1 × 10^-3). To obtain comprehensive genetic values for the Korean population, we combined data obtained in this study with previously reported values, thereby enabling us to estimate the locus-specific mutation rates regarding 22,711 allele transmissions. By combining these data, we obtained an overall average mutation rate of 2.91 × 10^-3 (95% CI, 2.3-3.7 × 10^-3). In addition, among the 476 unrelated Korean males, we detected 467 different haplotypes, with an overall haplotype diversity value of 0.9999. By extracting haplotypes of Y-STRs described in previous literature on 23 Y-STR reported in Korea, we obtained gene diversities for 1133 Korean individuals. We believe that the values and characteristics of the 23 Y-STRs analyzed in this study will contribute to establishing criteria for forensic genetic interpretation, including kinship analysis.

Revisiting informed consent in forensic genomics in light of current technologies and the times

Informed consent is based on basic ethical principles that should be considered when conducting biomedical and behavioral research involving human subjects. These principles-respect, beneficence, and justice-form the foundations of informed consent which in itself is grounded on three fundamental elements: information, comprehension, and voluntary participation. While informed consent has focused on human subjects and research, the practice has been adopted willingly in the forensic science arena primarily to acquire reference samples from family members to assist in identifying missing persons. With advances in molecular biology technologies, data mining, and access to metadata, it is important to assess whether the past informed consent process and in particular associated risks are concomitant with these increased capabilities. Given the state-of-the-art, areas in which informed consent may need to be modified and augmented are as follows: reference samples from family members in missing persons or unidentified human remains cases; targeted analysis of an individual(s) during forensic genetic genealogy cases to reduce an investigative burden; donors who provide their samples for validation studies (to include population studies and entry into databases that would be applied to forensic statistical calculations) to support implementation of procedures and operations of the forensic laboratory; family members that may contribute samples or obtain genetic information from a molecular autopsy; and use of medical and other acquired samples that could be informative for identification purposes. The informed consent process should cover (1) purpose for collection of samples; (2) process to analyze the samples (to include type of data); (3) benefits (to donor, target, family, community, etc. as applicable); (4) risks (to donor, target, family, community, etc. as applicable); (5) access to data/reports by the donor; (6) sample disposition; (7) removal of data process (i.e., expungement); (8) process to ask questions/assessment of comprehension; (9) follow-up processes; and (10) voluntary, signed, and dated consent. Issues surrounding these topics are discussed with an emphasis on addressing risk factors. Addressing informed consent will allow human subjects to make decisions voluntarily and with autonomy as well as secure the use of samples for intended use.

The mutational dynamics of short tandem repeats in large, multigenerational families

BACKGROUND

Short tandem repeats (STRs) compose approximately 3% of the genome, and mutations at STR loci have been linked to dozens of human diseases including amyotrophic lateral sclerosis, Friedreich ataxia, Huntington disease, and fragile X syndrome. Improving our understanding of these mutations would increase our knowledge of the mutational dynamics of the genome and may uncover additional loci that contribute to disease. To estimate the genome-wide pattern of mutations at STR loci, we analyze blood-derived whole-genome sequencing data for 544 individuals from 29 three-generation CEPH pedigrees. These pedigrees contain both sets of grandparents, the parents, and an average of 9 grandchildren per family.

RESULTS

We use HipSTR to identify de novo STR mutations in the 2nd generation of these pedigrees and require transmission to the third generation for validation. Analyzing approximately 1.6 million STR loci, we estimate the empirical de novo STR mutation rate to be 5.24 × 10^-5 mutations per locus per generation. Perfect repeats mutate about 2 × more often than imperfect repeats. De novo STRs are significantly enriched in Alu elements.

CONCLUSIONS

Approximately 30% of new STR mutations occur within Alu elements, which compose only 11% of the genome, but only 10% are found in LINE-1 insertions, which compose 17% of the genome. Phasing these mutations to the parent of origin shows that parental transmission biases vary among families. We estimate the average number of de novo genome-wide STR mutations per individual to be approximately 85, which is similar to the average number of observed de novo single nucleotide variants.

Assessing the factors influencing the performance of machine learning for classifying haplogroups from Y-STR haplotypes

Two distinct genetic markers, single nucleotide polymorphisms (Y-SNPs) and short tandem repeats (Y-STRs), exist simultaneously in the non-recombining portion of the Y chromosome. Because of their different rates of mutation, Y-STRs and Y-SNPs play distinct roles in forensic and evolutionary genetics. Current approaches to infer haplogroup status rely on genotyping lots of Y-SNP loci. Given the relationship between haplotype and haplogroup of a Y chromosome, a cost-effective strategy of Y-STRs typing had an advantage in haplogroup prediction. Many machine learning algorithms have sprung up for assigning a Y-STR haplotype to a haplogroup. However, a series of issues must be solved before the using of machine learning method in practice. Thus, the k-nearest neighbor (kNN) classifier was built respectively based on different situations in this study. We assessed different factors which may influence the performance of the kNN prediction model for classifying haplogroups. The training set was based on a diverse ground-truth data set comprising Y-STR haplotypes and corresponding Y-SNP haplogroups. Our results showed that combining different levels of haplogroups into the observations or transracial prediction was impractical. Moreover, using more slow mutation Y-STR loci in the category is good for promoting classification accuracy. The preconditions for an effective and accurate haplogroup assignment by the kNN classifier were revealed.

Resolving the recombination pattern of 38 X-STRs from Chinese Han three-generation pedigrees

X-chromosomal markers have been proved as a useful tool for solving complex kinship cases due to its sex-linked inheriting feature. Among these markers, tightly linked X-STR clusters forming haplotypes are highly informative. The analysis of the haplotypes requires determination of linkage disequilibrium. In this study, genetic linkage, recombination fractions and mutation rates of 38 X-STR loci in 177 three-generation pedigrees were investigated. Genetic linkage analysis and calculation of recombination fractions were performed within each pair of markers and clusters. Then mutation rates were calculated. The results showed that, a) 22 recombination events happened within the tightly linked X-STR clusters, which span<1.0 Mb; b) significantly linked marker pairs were observed with the LOD (logarithm of the odds) scores > 2.0 (2.0104 to 54.8316); c) the average mutation rate of the 38 X-STR loci was 1.32 × 10^-3 per meiosis in the Chinese Han population, with DXS10135 and DXS8377 presenting notably high mutation rate (6.5 × 10^-3). Our results confirmed that meiotic recombination was not a simple function of physical distance, so that whether recombination occurred at the closely clustered X-STRs or not should be assumed cautiously considering the stability of haplotypes in inheritance process for kinship analysis. This study supplemented the existing database and laid an experimental foundation for the future study on genetic characteristics, recombination, and mutation of the X-STRs.

Short Tandem Repeat Variation in the CNR1 Gene Associated With Analgesic Requirements of Opioids in Postoperative Pain Management

Short tandem repeats (STRs) and variable number of tandem repeats (VNTRs) that have been identified at approximately 0.7 and 0.5 million loci in the human genome, respectively, are highly multi-allelic variations rather than single-nucleotide polymorphisms. The number of repeats of more than a few thousand STRs was associated with the expression of nearby genes, indicating that STRs are influential genetic variations in human traits. Analgesics act on the central nervous system via their intrinsic receptors to produce analgesic effects. In the present study, we focused on STRs and VNTRs in the CNR1, GRIN2A, PENK, and PDYN genes and analyzed two peripheral pain sensation-related traits and seven analgesia-related traits in postoperative pain management. A total of 192 volunteers who underwent the peripheral pain sensation tests and 139 and 252 patients who underwent open abdominal and orthognathic cosmetic surgeries, respectively, were included in the study. None of the four STRs or VNTRs were associated with peripheral pain sensation. Short tandem repeats in the CNR1, GRIN2A, and PENK genes were associated with the frequency of fentanyl use, fentanyl dose, and visual analog scale pain scores 3 h after orthognathic cosmetic surgery (Spearman’s rank correlation coefficient ρ = 0.199, p = 0.002, ρ = 0.174, p = 0.006, and ρ = 0.135, p = 0.033, respectively), analgesic dose, including epidural analgesics after open abdominal surgery (ρ = −0.200, p = 0.018), and visual analog scale pain scores 24 h after orthognathic cosmetic surgery (ρ = 0.143, p = 0.023), respectively. The associations between STRs in the CNR1 gene and the frequency of fentanyl use and fentanyl dose after orthognathic cosmetic surgery were confirmed by Holm’s multiple-testing correction. These findings indicate that STRs in the CNR1 gene influence analgesia in the orofacial region.

Cloning and Selfing Affect Population Genetic Variation in Simulations of Outcrossing, Sexual Sea Stars

AbstractMany sea stars are well known for facultative or obligate asexual reproduction in both the adult and larval life-cycle stages. Some species and lineages are also capable of facultative or obligate hermaphroditic reproduction with self-fertilization. However, models of population genetic variation and empirical analyses of genetic data typically assume only sexual reproduction and outcrossing. A recent reanalysis of previously published empirical data (microsatellite genotypes) from two studies of one of the most well-known sea star species (the crown-of-thorns sea star; Acanthaster sp.) concluded that cloning and self-fertilization in that species are rare and contribute little to patterns of population genetic variation. Here we reconsider that conclusion by simulating the contribution of cloning and selfing to genetic variation in a series of models of sea star demography. Simulated variation in two simple models (analogous to previous analyses of empirical data) was consistent with high rates of cloning or selfing or both. More realistic scenarios that characterize population flux in sea stars of ecological significance, including outbreaks of crown-of-thorns sea stars that devastate coral reefs, invasions by Asterias amurensis, and epizootics of sea star wasting disease that kill Pisaster ochraceus, also showed significant but smaller effects of cloning and selfing on variation within subpopulations and differentiation between subpopulations. Future models or analyses of genetic variation in similar study systems might benefit from simulation modeling to characterize possible contributions of cloning or selfing to genetic variation in population samples or to understand the limits on inferring the effects of cloning or selfing in nature.

Uniparental Lineages from the Oldest Indigenous Population of Ecuador: The Tsachilas

Together with Cayapas, the Tsachilas constitute the oldest population in the country of Ecuador and, according to some historians, they are the last descendants of the ancient Yumbos. Several anthropological issues underlie the interest towards this peculiar population: the uncertainty of their origin, their belonging to the Barbacoan linguistic family, which is still at the center of an intense linguistic debate, and the relations of their Yumbo ancestors with the Inca invaders who occupied their ancient territory. Our contribution to the knowledge of their complex past was the reconstruction of their genetic maternal and paternal inheritance through the sequencing of 70 entire mitochondrial genomes and the characterization of the non-recombinant region of the Y chromosome in 26 males. For both markers, we built comprehensive datasets of various populations from the surrounding geographical area, northwestern South America, NW, with a known linguistic affiliation, and we could then compare our sample against the overall variability to infer relationships with other Barbacoan people and with other NW natives. We found contrasting patterns of genetic diversity for the two markers, but generally, our results indicated a possible common origin between the Tsachilas, the Chachi, and other Ecuadorian and Colombian Barbacoans and are suggestive of an interesting ancient linkage to the Inca invaders in Yumbo country.

Improved Models of Coalescence Ages of Y-DNA Haplogroups

Databases of commercial DNA-testing companies now contain more customers with sequenced DNA than any completed academic study, leading to growing interest from academic and forensic entities. An important result for both these entities and the test takers themselves is how closely two individuals are related in time, as calculated through one or more molecular clocks. For Y-DNA, existing interpretations of these clocks are insufficiently accurate to usefully measure relatedness in historic times. In this article, I update the methods used to calculate coalescence ages (times to most-recent common ancestor, or TMRCAs) using a new, probabilistic statistical model that includes Y-SNP, Y-STR and ancilliary historical data, and provide examples of its use.

Next generation sequencing of Y-STRs in father-son pairs and comparison with traditional capillary electrophoresis

To evaluate the promising advantages of massively parallel sequencing (MPS) in our casework, we analysed a total of 33 Y-chromosomal short tandem repeats (Y-STRs) with traditional capillary electrophoresis (CE) and 25 Y-STRs using the newer MPS technology. We studied the outcome of both technologies in 64 father-son pairs using stock and custom-designed kits. Current MPS technology confirmed the 13 mutational events observed with CE and improved our understanding of the complex nature of STR mutations. By detecting isometric sequence variants between unrelated males, we show that sequencing Y-STRs using MPS can boost discrimination power.

Current genetic conservation of Chinese indigenous horses revealed with Y-chromosomal and mitochondrial DNA polymorphisms

To investigate the genetic diversity of Chinese indigenous horses and determine the genetic status of extant horse breeds, novel Y chromosomal microsatellite markers and known Y chromosomal SNPs and mtDNA loop sequences, were employed to study the genetic diversity levels of 13 Chinese indigenous horse populations and four introduced breeds. Sixteen Y-chromosomal microsatellite markers, including seven newly identified loci, were used in the genotyping. The results showed that 4 out of the 16 loci were highly polymorphic in Chinese indigenous horse populations, in which the polymorphisms of 3 loci, ECAYP12, ECAYP13, and ECAYCAU3, were first reported in the present study. The polymorphic Y chromosomal microsatellite markers result in 19 haplotypes in the studied horses and formed 24 paternal lines when merged with the 14 Y chromosomal SNPs reported previously. The haplotypes CHT18 and SS24 harboring AMELY gene mutation were the ancestral haplotypes, and other haplotypes were derived from them by one or more mutation steps. The horse populations in mountainous and remote areas of southwestern China have the most ancient paternal lines, which suggests that ancient paternal lines preserved in local populations attributed to less human interventions. Our results also showed that the northern local breeds had higher mtDNA diversity than the southern ones in China. The frequency of haplogroup B, F, and G of mtDNA in Chinese indigenous horses has declined in recent years, and some breeds are in endangered status mainly due to small population sizes. Urgent actions should be taken to conserve the genetic diversity of the indigenous horse populations, especially the rare paternal lines. Our findings help to elucidate the genetic diversity and evolutionary history of Chinese domestic horses, which will facilitate the conservation of the indigenous horses in the future.

Application of CPI cutoff value based on parentage testing of duos and trios typed by four autosomal kits

In this study, we analyzed the application of four autosomal kits and the sensitivity of the combined paternity index (CPI) cutoff value (CPI≥10000) in parentage testing. First, 1442 real trios and 803 real duos were tested using the Goldeneye 25A kit. The Goldeneye 25A kit covers the autosomal short tandem repeat (STR) loci of the other three kits, so we calculated the CPI value of every case for the four kits. Second, three complex close relative kinship cases were also analyzed to evaluate the application of the CPI cutoff value. The CPI values of all trio cases were higher than 10000 using the four kits; the CPI values of all duo cases were higher than 10000 using the Goldeneye 25A kit; and the CPI values of a portion of the duo cases were lower than 10000 using the other three kits. In the three complex close relative cases, the alleged father or mother was not excluded using 40 autosomal STRs. Adding X chromosome short tandem repeats (X-STR) and samples of biological fathers or mothers, the conclusions were confirmed. The four kits were adequate to draw conclusions in the trio cases; the Goldeneye 25A Kit was adequate to draw conclusions in the duo cases; and the other three kits were not sufficient for a portion of the duo cases. The CPI cutoff value was sensitive for real trio and duo cases. In complex close relative kinship cases, high CPI values may result in false conclusions.

Forensic characteristics and genetic analysis of both 27 Y-STRs and 143 Y-SNPs in Eastern Han Chinese population

Y-chromosome short tandem repeat (Y-STR) and Y-chromosome single nucleotide polymorphism (Y-SNP) frequency distributions provide resources for assessment of male population stratification among world-wide populations. Currently, the Y-STR Haplotype Reference Database (YHRD) contains numerous Y-chromosome haplotype profiles from various populations and countries around the world. However, for many of the recently discovered and already phylogenetically mapped Y-SNPs, the population data are scarce. Herein, the typing of 27 Y-STRs (Yfiler Plus) and 143 Y-SNPs (self-designed Y-SNP panel) was performed on 1269 unrelated males from 11 Han Chinese populations. Haplogroup O-M175 was the most predominant haplogroup in our Han Chinese data, ranging from 67.34% (Henan Han) to 93.16% (Guangdong Han). The highest haplogroup diversity (0.967056) was observed in Heilongjiang Han, with a discrimination capacity (DC) value of 0.3723. The number of alleles at single-copy loci varied from 2 for DYS391 (Guangdong Han) to 16 for DYS518 (Henan Han). For the majority of the populations (8/11), both the haplotype diversity and DC values are 1.0000. Furthermore, genetic differentiations were observed between Northern and Southern Han Chinese. These genetic differences were mainly reflected in haplogroup distribution and frequency, and they were confirmed by statistical analysis.

Suppressed Recombination of Sex Chromosomes Is Not Caused by Chromosomal Reciprocal Translocation in Spiny Frog (Quasipaa boulengeri)

Chromosome rearrangements (CRs) are perceived to be related to sex chromosome evolution, but it is a matter of controversy whether CRs are the initial causative mechanism of suppressed recombination for sex differentiation. The early stages of sex chromosome evolution in amphibians may represent intermediate states of differentiation, and if so, they potentially shed light on the ultimate cause of suppressed recombination and the role of CRs in sex chromosome differentiation. In this paper, we showed that sex determination differs among 16 populations of spiny frog (Quasipaa boulengeri), in which individuals have normal and rearranged chromosomes caused by reciprocal translocation. In eastern areas, without translocation, genetic differentiation between sexes was relatively low, suggesting unrestricted recombination. In comparison, in western populations that have both normal and translocated chromosomes, a male-heterogametic system and lack of X-Y recombination were identified by male-specific alleles and heterozygote excess. However, such genetic differentiation between sexes in western populations was not directly related to karyotypes, as it was found in individuals with both normal and translocated karyotypes. In the western Sichuan Basin, male-specific and translocation-specific allelic frequency distributions suggested that recombination of sex-differentiation ceased in all populations, but recombination suppression caused by translocation did not exist in some populations. Combined with phylogenetic inference, this indicated that the establishment of sex-linkage had taken place independently of reciprocal translocation, and translocation was not the ultimate cause of sex chromosome differentiation. Furthermore, comparison of the genetic diversity of alleles on Y chromosomes, X chromosomes, and autosomes in western populations showed a reduction of effective population size on sex chromosomes, which may be caused by reciprocal translocation. It indicates that, although it is not the ultimate cause of recombination suppression, reciprocal translocation may enhance sex chromosome differentiation.

Haplotypic polymorphisms and mutation rate estimates of 22 Y-chromosome STRs in the Northern Chinese Han father-son pairs

Y chromosome Short tandem repeats (Y-STRs) analysis has been widely used in forensic identification, kinship testing, and population evolution. An accurate understanding of haplotype and mutation rate will benefit these applications. In this work, we analyzed 1123 male samples from Northern Chinese Han population which including 578 DNA-confirmed father-son pairs at 22 Y-STRs loci. A total of 537 haplotypes were observed and the overall haplotype diversity was calculated as 1.0000 ± 0.0001. Except that only two haplotypes were observed twice, all the rest of the 535 were unique. Furthermore, totally 47 mutations were observed during 13,872 paternal meiosis. The mutation rate for each locus estimates ranged from 0.0 to 15.6 × 10-3 with an average mutation rate 3.4 × 10-3 (95% CI 2.5-4.5 × 10-3). Among the 22 loci, DYS449, DYS389 II and DYS458 are the most prone to mutations. This study adds to the growing data on Y-STR haplotype diversity and mutation rates and could be very useful for population and forensic genetics.

Expanded Croatian 12 X-STR loci database with an overview of anomalous profiles

In order to implement X-chromosome short tandem repeat (X-STR) typing into routine forensic practice, reference database of a given population should be established. Therefore we extended already published data with additional 397 blood samples from unrelated Croatian citizens, and analyzed the total of 995 samples (549 male and 446 female) typed by Investigator^® Argus X-12 Kit. To test genetic homogeneity of consecutively processed five historic-cultural regions covering the entire national territory, we calculated pairwise Fst genetic distances between regions based on allele and full haplotype frequencies. Since the comparison did not yield any statistically significant difference, we integrated STR profile information from all regions and used the whole data set to calculate forensic parameters. The most informative marker is DXS10135 (polymorphism information content (PIC = 0.929) and the most informative linkage group (LG) is LG1 (PIC = 0.996). We confirmed linkage disequilibrium (LD) for seven marker pairs belonging to LG2, LG3 and LG4. By including LD information, we calculated cumulative power of discrimination that amounted to 0.999999999997 in females and 0.999999005 in males. We also compared Croatia with 13 European populations based on haplotype frequencies and detected no statistically significant Fst values after Bonferroni correction in any LG. Multi-dimensional scaling plot revealed tight grouping of four Croatian regions amongst populations of southern, central and northern Europe, with the exception of northern Croatia. In this study we gave the first extensive overview of aberrant profiles encountered during Investigator^® Argus X-12 typing. We found ten profiles consistent with single locus duplication followed by tetranucleotide tract length polymorphism. Locus DXS10079 is by far the most frequently affected one, presumably mutated in eight samples. We also found four profiles consistent with X-chromosome aneuploidy (three profiles with XXX pattern and one profile with XXY pattern). In conclusion, we established integral forensic Croatian X-chromosome database, proved forensic pertinence of Investigator^® Argus X-12 Kit for the entire Croatian population and identified locus DXS10079 as a potential duplication hotspot.

Prostate cancer susceptibility and growth linked to Y chromosome genes

The role of Y chromosome in prostate cancer progression and incidence is not well known. Among the 46 chromosomes, Y chromosome determines the male gender. The Y chromosome is smaller than the X chromosome and contains only 458 genes compared to over 2000 genes found in the X chromosome. The Y chromosome is prone to high mutation rates, created exclusively in sperm cells due to the highly oxidative environment of the testis. Y chromosome harbors epigenetic information, which affects the expression of genes associated with the incidence and progression of prostate cancer. In this review, we focus on Y chromosome related genetic abnormalities, likely to be involved in the development and progression of prostate cancer.

Modular knowledge systems accelerate human migration in asymmetric random environments

Migration is a key mechanism for expansion of communities. In spatially heterogeneous environments, rapidly gaining knowledge about the local environment is key to the evolutionary success of a migrating population. For historical human migration, environmental heterogeneity was naturally asymmetric in the north-south (NS) and east-west (EW) directions. We here consider the human migration process in the Americas, modelled as random, asymmetric, modularly correlated environments. Knowledge about the environments determines the fitness of each individual. We present a phase diagram for asymmetry of migration as a function of carrying capacity and fitness threshold. We find that the speed of migration is proportional to the inverse complement of the spatial environmental gradient, and in particular, we find that NS migration rates are lower than EW migration rates when the environmental gradient is higher in the NS direction. Communication of knowledge between individuals can help to spread beneficial knowledge within the population. The speed of migration increases when communication transmits pieces of knowledge that contribute in a modular way to the fitness of individuals. The results for the dependence of migration rate on asymmetry and modularity are consistent with existing archaeological observations. The results for asymmetry of genetic divergence are consistent with patterns of human gene flow.

The Genetic Legacy of Zoroastrianism in Iran and India: Insights into Population Structure, Gene Flow, and Selection

Zoroastrianism is one of the oldest extant religions in the world, originating in Persia (present-day Iran) during the second millennium BCE. Historical records indicate that migrants from Persia brought Zoroastrianism to India, but there is debate over the timing of these migrations. Here we present genome-wide autosomal, Y chromosome, and mitochondrial DNA data from Iranian and Indian Zoroastrians and neighboring modern-day Indian and Iranian populations and conduct a comprehensive genome-wide genetic analysis in these groups. Using powerful haplotype-based techniques, we find that Zoroastrians in Iran and India have increased genetic homogeneity relative to other sampled groups in their respective countries, consistent with their current practices of endogamy. Despite this, we infer that Indian Zoroastrians (Parsis) intermixed with local groups sometime after their arrival in India, dating this mixture to 690–1390 CE and providing strong evidence that Iranian Zoroastrian ancestry was maintained primarily through the male line. By making use of the rich information in DNA from ancient human remains, we also highlight admixture in the ancestors of Iranian Zoroastrians dated to 570 BCE–746 CE, older than admixture seen in any other sampled Iranian group, consistent with a long-standing isolation of Zoroastrians from outside groups. Finally, we report results, and challenges, from a genome-wide scan to identify genomic regions showing signatures of positive selection in present-day Zoroastrians that might correlate to the prevalence of particular diseases among these communities.

East of the Andes: The genetic profile of the Peruvian Amazon populations

OBJECTIVES

Assuming that the differences between the Andes and the Amazon rainforest at environmental and historical levels have influenced the distribution patterns of genes, languages, and cultures, the maternal and paternal genetic reconstruction of the Peruvian Amazon populations was used to test the relationships within and between these two extreme environments.

MATERIALS AND METHODS

We analyzed four Peruvian Amazon communities (Ashaninka, Huambisa, Cashibo, and Shipibo) for both Y chromosome (17 STRs and 8 SNPs) and mtDNA data (control region sequences, two diagnostic sites of the coding region, and one INDEL), and we studied their variability against the rest of South America.

RESULTS

We detected a high degree of genetic diversity in the Peruvian Amazon people, both for mtDNA than for Y chromosome, excepting for Cashibo people, who seem to have had no exchanges with their neighbors, in contrast with the others communities. The genetic structure follows the divide between the Andes and the Amazon, but we found a certain degree of gene flow between these two environments, as particularly emerged with the Y chromosome descent cluster's (DCs) analysis.

DISCUSSION

The Peruvian Amazon is home to an array of populations with differential rates of genetic exchanges with their neighbors and with the Andean people, depending on their peculiar demographic histories. We highlighted some successful Y chromosome lineages expansions originated in Peru during the pre-Columbian history which involved both Andeans and Amazon Arawak people, showing that at least a part of the Amazon rainforest did not remain isolated from those exchanges.

Population-Scale Sequencing Data Enable Precise Estimates of Y-STR Mutation Rates

Short tandem repeats (STRs) are mutation-prone loci that span nearly 1% of the human genome. Previous studies have estimated the mutation rates of highly polymorphic STRs by using capillary electrophoresis and pedigree-based designs. Although this work has provided insights into the mutational dynamics of highly mutable STRs, the mutation rates of most others remain unknown. Here, we harnessed whole-genome sequencing data to estimate the mutation rates of Y chromosome STRs (Y-STRs) with 2-6 bp repeat units that are accessible to Illumina sequencing. We genotyped 4,500 Y-STRs by using data from the 1000 Genomes Project and the Simons Genome Diversity Project. Next, we developed MUTEA, an algorithm that infers STR mutation rates from population-scale data by using a high-resolution SNP-based phylogeny. After extensive intrinsic and extrinsic validations, we harnessed MUTEA to derive mutation-rate estimates for 702 polymorphic STRs by tracing each locus over 222,000 meioses, resulting in the largest collection of Y-STR mutation rates to date. Using our estimates, we identified determinants of STR mutation rates and built a model to predict rates for STRs across the genome. These predictions indicate that the load of de novo STR mutations is at least 75 mutations per generation, rivaling the load of all other known variant types. Finally, we identified Y-STRs with potential applications in forensics and genetic genealogy, assessed the ability to differentiate between the Y chromosomes of father-son pairs, and imputed Y-STR genotypes.

Analysis of genetic admixture in Uyghur using the 26 Y-STR loci system

The Uyghur population has experienced extensive interaction with European and Eastern Asian populations historically. A set of high-resolution genetic markers could be useful to infer the genetic relationships between the Uyghur population and European and Asian populations. In this study we typed 100 unrelated Uyghur males living in southern Xinjiang at 26 Y-STR loci. Using the high-resolution 26 Y-STR loci system, we investigated genetic and phylogenetic relationship between the Uyghur population and 23 reference European or Asian populations. We found that the Uyghur population exhibited a genetic admixture of Eastern Asian and European populations, and had a slightly closer relationship with the selected European populations than the Eastern Asian populations. We also demonstrated that the 26 Y-STR loci system was potentially useful in forensic sciences because it has a large power of discrimination and rarely exhibits common haplotypes. However, ancestry inference of Uyghur samples could be challenging due to the admixed nature of the population.

Evolution of Nine Microsatellite Loci in the Fungus Fusarium oxysporum

The evolution of nine microsatellites and one minisatellite was investigated in the fungus Fusarium oxysporum and sister taxa Fusarium redolens and Fusarium verticillioides. Compared to other organisms, fungi have been reported to contain fewer and less polymorphic microsatellites. Mutational patterns over evolutionary time were studied for these ten loci by mapping changes in core repeat numbers onto a phylogeny based on the sequence of the conserved translation elongation factor 1-α gene. The patterns of microsatellite formation, expansion, and interruption by base substitutions were followed across the phylogeny, showing that these loci are evolving in a manner similar to that of microsatellites in other eukaryotes. Most mutations could be fit to a stepwise mutation model, but a few appear to have involved multiple repeat units. No evidence of gene conversion was seen at the minisatellite locus, which may also be mutating by replication slippage. Some homoplastic numbers of repeat units were observed for these loci, and polymorphisms in the regions flanking the microsatellites may provide better genetic markers for population genetics studies of these species.

Y-Chromosomal Lineages of Latvians in the Context of the Genetic Variation of the Eastern-Baltic Region

Variations of the nonrecombining Y-chromosomal region were investigated in 159 unrelated Baltic-speaking ethnic Latvians from four different geographic regions, using 28 biallelic markers and 12 short tandem repeats. Eleven different haplogroups (hgs) were detected in a regionally homogeneous Latvian population, among which N1c, R1a, and I1 cover more than 85% of its paternal lineages. When compared its closest geographic neighbors, the composition of the Latvian Y-chromosomal gene pool was found to be very similar to those of Lithuanians and Estonians. Despite the comparable frequency distribution of hg N1c in Latvians and Lithuanians with the Finno-Ugric-speaking populations from the Eastern coast of the Baltic Sea, the observed differences in allelic variances of N1c haplotypes between these two groups are in concordance with the previously stated hypothesis of different dispersal ways of this lineage in the region. More than a third of Latvian paternal lineages belong specifically to a recently defined R1a-M558 hg, indicating an influence from a common source within Eastern Slavic populations on the formation of the present-day Latvian Y-chromosome gene pool.

Y-chromosome descent clusters and male differential reproductive success: young lineage expansions dominate Asian pastoral nomadic populations

High-frequency microsatellite haplotypes of the male-specific Y-chromosome can signal past episodes of high reproductive success of particular men and their patrilineal descendants. Previously, two examples of such successful Y-lineages have been described in Asia, both associated with Altaic-speaking pastoral nomadic societies, and putatively linked to dynasties descending, respectively, from Genghis Khan and Giocangga. Here we surveyed a total of 5321 Y-chromosomes from 127 Asian populations, including novel Y-SNP and microsatellite data on 461 Central Asian males, to ask whether additional lineage expansions could be identified. Based on the most frequent eight-microsatellite haplotypes, we objectively defined 11 descent clusters (DCs), each within a specific haplogroup, that represent likely past instances of high male reproductive success, including the two previously identified cases. Analysis of the geographical patterns and ages of these DCs and their associated cultural characteristics showed that the most successful lineages are found both among sedentary agriculturalists and pastoral nomads, and expanded between 2100 BCE and 1100 CE. However, those with recent origins in the historical period are almost exclusively found in Altaic-speaking pastoral nomadic populations, which may reflect a shift in political organisation in pastoralist economies and a greater ease of transmission of Y-chromosomes through time and space facilitated by the use of horses.

Gene conversion violates the stepwise mutation model for microsatellites in y-chromosomal palindromic repeats

The male-specific region of the human Y chromosome (MSY) contains eight large inverted repeats (palindromes), in which high-sequence similarity between repeat arms is maintained by gene conversion. These palindromes also harbor microsatellites, considered to evolve via a stepwise mutation model (SMM). Here, we ask whether gene conversion between palindrome microsatellites contributes to their mutational dynamics. First, we study the duplicated tetranucleotide microsatellite DYS385a,b lying in palindrome P4. We show, by comparing observed data with simulated data under a SMM within haplogroups, that observed heteroallelic combinations in which the modal repeat number difference between copies was large, can give rise to homoallelic combinations with zero-repeats difference, equivalent to many single-step mutations. These are unlikely to be generated under a strict SMM, suggesting the action of gene conversion. Second, we show that the intercopy repeat number difference for a large set of duplicated microsatellites in all palindromes in the MSY reference sequence is significantly reduced compared with that for nonpalindrome-duplicated microsatellites, suggesting that the former are characterized by unusual evolutionary dynamics. These observations indicate that gene conversion violates the SMM for microsatellites in palindromes, homogenizing copies within individual Y chromosomes, but increasing overall haplotype diversity among chromosomes within related groups.

Usefulness of SNPs as Supplementary Markers in a Paternity Case with 3 Genetic Incompatibilities at Autosomal and Y Chromosomal Loci

BACKGROUND

In kinship testing, investigation of 15 short tandem repeats (STRs) usually provides decisive genetic information for resolving relationship cases. However, in complex deficiency cases, in cases with more than 2 mutations at different STR loci or when close (untested) relatives of the alleged father are suggested to be the biological father of the child, STR typing alone may not be sufficient. In these cases, the application of supplementary markers such as single nucleotide polymorphisms (SNPs) is recommended.

METHODS

We describe a paternity case with 3 genetic incompatibilities (Penta D, VWA, and DYS385) between the alleged father and the child after analyzing 23 autosomal and 16 Y chromosomal STR loci. The question arose as to whether the alleged father could be excluded and a related person could be the biological father of the child, or whether the observed genetic incompatibilities were mutations. Interestingly, the 2 excluded full brothers of the alleged father possessed identical genetic incompatibilities at locus VWA and DYS385 as the alleged father.

RESULTS AND CONCLUSIONS

Additional performance of a 50-plex SNP assay demonstrated that the observed mismatches were indeed mutations and the alleged father was the biological father of the child. The results show the usefulness of SNPs as supplementary markers in relationship testing when STR analyses show ambiguous results.

Mitochondrial DNA variation, genetic structure and demographic history of Iranian populations

In order to survey the evolutionary history and impact of historical events on the genetic structure of Iranian people, the HV2 region of 141 mtDNA sequences related to six Iranian populations were analyzed. Slight and non-significant F_ST distances among the Central-western Persian speaking populations of Iran testify to the common origin of these populations from one proto-population. Mismatch distribution suggests that this proto-Iranian population started to colonize Iran about 30000 years ago which is almost consistent with the timing of arrival and colonization of western Asia by the anatomically modern human. Star-like haplotype network structures, significant and negative Tajima's D (D=-2.08, P<0.05) and unimodal mismatch distributions support the genetic effects of this expansion. Iranian populations presented mtDNA lineages that clearly belong to the European gene pool (i.e. H and U), while the Mashhad population was characterized by the presence of eastern and central Asian mtDNA lineages (i.e. M, B and D). Furthermore, the low diversity (h=0.428) observed in Mashhad may indicated the presence of inbreeding, drift or bottleneck events. The application of Monmonier's maximum differences algorithm revealed a geographic zone of genetic discontinuity between the Arab people of Khuzestan and rest of Iranian populations. Geographical factors, in cooperation with cultural/linguistic differences, are the main reasons for this differentiation. The lack of a sharp geographical or ethno-linguistic structure for mtDNA HV2 sequence diversity was statistically supported by AMOVA and Mantel (r=0.19, P<0.05) tests.

Y-chromosomal analysis identifies the skeletal remains of Swiss national hero Jörg Jenatsch (1596-1639)

Jörg Jenatsch was a Swiss defender of independence and a fighter for liberty in the 17th century. With the help of three living male members of the Jenatsch family, we successfully identified a skeleton exhumed from Chur cathedral as the remains of Jörg Jenatsch. Our conclusion was based upon complete Y-STR and Y-SNP profiles that could be generated by replicate analyses of a bone sample available to us. The skeleton and the three living family members carried the same Y-SNP haplogroup, but were discordant at three of 23 Y-STR loci. This notwithstanding, conservative biostatistical evaluation of the data suggests that the Chur skeleton is at least 20 times more likely than not to be Jörg Jenatsch.