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

An African American paternal lineage adds an extremely ancient root to the human Y chromosome phylogenetic tree

We report the discovery of an African American Y chromosome that carries the ancestral state of all SNPs that defined the basal portion of the Y chromosome phylogenetic tree. We sequenced ∼240 kb of this chromosome to identify private, derived mutations on this lineage, which we named A00. We then estimated the time to the most recent common ancestor (TMRCA) for the Y tree as 338 thousand years ago (kya) (95% confidence interval = 237-581 kya). Remarkably, this exceeds current estimates of the mtDNA TMRCA, as well as those of the age of the oldest anatomically modern human fossils. The extremely ancient age combined with the rarity of the A00 lineage, which we also find at very low frequency in central Africa, point to the importance of considering more complex models for the origin of Y chromosome diversity. These models include ancient population structure and the possibility of archaic introgression of Y chromosomes into anatomically modern humans. The A00 lineage was discovered in a large database of consumer samples of African Americans and has not been identified in traditional hunter-gatherer populations from sub-Saharan Africa. This underscores how the stochastic nature of the genealogical process can affect inference from a single locus and warrants caution during the interpretation of the geographic location of divergent branches of the Y chromosome phylogenetic tree for the elucidation of human origins.

Detecting genetic isolation in human populations: a study of European language minorities

The identification of isolation signatures is fundamental to better understand the genetic structure of human populations and to test the relations between cultural factors and genetic variation. However, with current approaches, it is not possible to distinguish between the consequences of long-term isolation and the effects of reduced sample size, selection and differential gene flow. To overcome these limitations, we have integrated the analysis of classical genetic diversity measures with a Bayesian method to estimate gene flow and have carried out simulations based on the coalescent. Combining these approaches, we first tested whether the relatively short history of cultural and geographical isolation of four “linguistic islands” of the Eastern Alps (Lessinia, Sauris, Sappada and Timau) had left detectable signatures in their genetic structure. We then compared our findings to previous studies of European population isolates. Finally, we explored the importance of demographic and cultural factors in shaping genetic diversity among the groups under study. A combination of small initial effective size and continued genetic isolation from surrounding populations seems to provide a coherent explanation for the diversity observed among Sauris, Sappada and Timau, which was found to be substantially greater than in other groups of European isolated populations. Simulations of micro-evolutionary scenarios indicate that ethnicity might have been important in increasing genetic diversity among these culturally related and spatially close populations.

Y-linked microsatellites in Amazonian Amerindians applied to ancestry estimates in Brazilian Afro-derived populations

OBJECTIVES

Proper ancestral populations are required to determine accurate ancestry estimates for Afro-derived Brazilian populations. Herein, we have genotyped Y-STRs in Amazonian Amerindians to determine the ancestral contribution in quilombo remnant communities.

METHODS

The frequencies for five Y-chromosome linked microsatellites (DYS19, DYS390, DYS391, DYS392, and DYS393) were characterized in four Amerindian tribes from Brazilian Amazon (Tikúna, Baníwa, Kashinawa, and Kanamarí), and in four quilombo remnants (Mimbó, Sítio Velho, Gaucinha, and São Gonçalo) and two urban populations (Teresina and Jequié) from Northeastern Brazil. We then estimated the male genetic ancestry in each admixed population. Moreover, we performed analysis of molecular variance (AMOVA), FST , haplotype diversity, and principal component analysis.

RESULTS

Lower haplotype diversity (h) values were observed for Tikúna compared with other tribes. Quilombo remnants exhibited higher h levels ranging from 0.893 ± 0.027 in Sítio Velho to 0.963 ± 0.033 in São Gonçalo. African ancestry estimates ranged from 0.529 ± 0.027 in Mimbó to 0.602 ± 0.086 in Sítio Velho. Conversely, European contribution was 0.795 ± 0.045 in Teresina and 0.826 ± 0.040 in Jequié.

CONCLUSIONS

FST and principal component analysis indicate homogeneity in the male genetic constitution among the quilombo remnants analyzed. Data on Amerindians allowed accurate ancestry estimates, which indicated a higher African contribution, followed by a considerable European contribution for these quilombo remnants.

Towards improvements in the estimation of the coalescent: implications for the most effective use of Y chromosome short tandem repeat mutation rates

Over the past two decades, many short tandem repeat (STR) microsatellite loci on the human Y chromosome have been identified together with mutation rate estimates for the individual loci. These have been used to estimate the coalescent age, or the time to the most recent common ancestor (TMRCA) expressed in generations, in conjunction with the average square difference measure (ASD), an unbiased point estimator of TMRCA based upon the average within-locus allele variance between haplotypes. The ASD estimator, in turn, depends on accurate mutation rate estimates to be able to produce good approximations of the coalescent age of a sample. Here, a comparison is made between three published sets of per locus mutation rate estimates as they are applied to the calculation of the coalescent age for real and simulated population samples. A novel evaluation method is developed for estimating the degree of conformity of any Y chromosome STR locus of interest to the strict stepwise mutation model and specific recommendations are made regarding the suitability of thirty-two commonly used Y-STR loci for the purpose of estimating the coalescent. The use of the geometric mean for averaging ASD and across loci is shown to improve the consistency of the resulting estimates, with decreased sensitivity to outliers and to the number of STR loci compared or the particular set of mutation rates selected.

A direct characterization of human mutation based on microsatellites

Mutations are the raw material of evolution, but have been difficult to study directly. We report the largest study of new mutations to date: 2,058 germline changes discovered by analyzing 85,289 Icelanders at 2,477 microsatellites. The paternal-to-maternal mutation rate ratio is 3.3, and the rate in fathers doubles from age 20 to 58 whereas there is no association with age in mothers. Longer microsatellite alleles are more mutagenic and tend to decrease in length, whereas the opposite is seen for shorter alleles. We use these empirical observations to build a model that we apply to individuals for whom we have both genome sequence and microsatellite data, allowing us to estimate key parameters of evolution without calibration to the fossil record. We infer that the sequence mutation rate is 1.4–2.3×10⁻⁸ per base pair per generation (90% credible interval), and that human-chimpanzee speciation occurred 3.7–6.6 million years ago.

In the name of the migrant father--analysis of surname origins identifies genetic admixture events undetectable from genealogical records

Patrilineal heritable surnames are widely used to select autochthonous participants for studies on small-scale population genetic patterns owing to the unique link between the surname and a genetic marker, the Y-chromosome (Y-chr). Today, the question arises as to whether the surname origin will be informative on top of in-depth genealogical pedigrees. Admixture events that happened in the period after giving heritable surnames but before the start of genealogical records may be informative about the additional value of the surname origin. In this context, an interesting historical event is the demic migration from French-speaking regions in Northern France to the depopulated and Dutch-speaking region Flanders at the end of the sixteenth century. Y-chr subhaplogroups of individuals with a French/Roman surname that could be associated with this migration event were compared with those of a group with autochthonous Flemish surnames. Although these groups could not be differentiated based on in-depth genealogical data, they were significantly genetically different from each other. Moreover, the observed genetic divergence was related to the differences in the distributions of main Y-subhaplogroups between contemporary populations from Northern France and Flanders. Therefore, these results indicate that the surname origin can be an important feature on top of in-depth genealogical results to select autochthonous participants for a regional population genetic study based on Y-chromosomes.

Pasture names with Romance and Slavic roots facilitate dissection of Y chromosome variation in an exclusively German-speaking alpine region

The small alpine district of East Tyrol (Austria) has an exceptional demographic history. It was contemporaneously inhabited by members of the Romance, the Slavic and the Germanic language groups for centuries. Since the Late Middle Ages, however, the population of the principally agrarian-oriented area is solely Germanic speaking. Historic facts about East Tyrol's colonization are rare, but spatial density-distribution analysis based on the etymology of place-names has facilitated accurate spatial mapping of the various language groups' former settlement regions. To test for present-day Y chromosome population substructure, molecular genetic data were compared to the information attained by the linguistic analysis of pasture names. The linguistic data were used for subdividing East Tyrol into two regions of former Romance (A) and Slavic (B) settlement. Samples from 270 East Tyrolean men were genotyped for 17 Y-chromosomal microsatellites (Y-STRs) and 27 single nucleotide polymorphisms (Y-SNPs). Analysis of the probands' surnames revealed no evidence for spatial genetic structuring. Also, spatial autocorrelation analysis did not indicate significant correlation between genetic (Y-STR haplotypes) and geographic distance. Haplogroup R-M17 chromosomes, however, were absent in region A, but constituted one of the most frequent haplogroups in region B. The R-M343 (R1b) clade showed a marked and complementary frequency distribution pattern in these two regions. To further test East Tyrol's modern Y-chromosomal landscape for geographic patterning attributable to the early history of settlement in this alpine area, principal coordinates analysis was performed. The Y-STR haplotypes from region A clearly clustered with those of Romance reference populations and the samples from region B matched best with Germanic speaking reference populations. The combined use of onomastic and molecular genetic data revealed and mapped the marked structuring of the distribution of Y chromosomes in an alpine region that has been culturally homogeneous for centuries.

Disclosing the genetic structure of Brazil through analysis of male lineages with highly discriminating haplotypes

In a large variety of genetic studies, probabilistic inferences are made based on information available in population databases. The accuracy of the estimates based on population samples are highly dependent on the number of chromosomes being analyzed as well as the correct representation of the reference population. For frequency calculations the size of a database is especially critical for haploid markers, and for countries with complex admixture histories it is important to assess possible substructure effects that can influence the coverage of the database. Aiming to establish a representative Brazilian population database for haplotypes based on 23 Y chromosome STRs, more than 2,500 Y chromosomes belonging to Brazilian, European and African populations were analyzed. No matter the differences in the colonization history of the five geopolitical regions that currently exist in Brazil, for the Y chromosome haplotypes of the 23 studied Y-STRs, a lack of genetic heterogeneity was found, together with a predominance of European male lineages in all regions of the country. Therefore, if we do not consider the diverse Native American or Afro-descendent isolates, which are spread through the country, a single Y chromosome haplotype frequency database will adequately represent the urban populations in Brazil. In comparison to the most commonly studied group of 17 Y-STRs, the 23 markers included in this work allowed a high discrimination capacity between haplotypes from non-related individuals within a population and also increased the capacity to discriminate between paternal relatives. Nevertheless, the expected haplotype mutation rate is still not enough to distinguish the Y chromosome profiles of paternally related individuals. Indeed, even for rapidly mutating Y-STRs, a very large number of markers will be necessary to differentiate male lineages from paternal relatives.

The impact of recent events on human genetic diversity

The historical record tells us stories of migrations, population expansions and colonization events in the last few thousand years, but what was their demographic impact? Genetics can throw light on this issue, and has mostly done so through the maternally inherited mitochondrial DNA (mtDNA) and the male-specific Y chromosome. However, there are a number of problems, including marker ascertainment bias, possible influences of natural selection, and the obscuring layers of the palimpsest of historical and prehistorical events. Y-chromosomal lineages are particularly affected by genetic drift, which can be accentuated by recent social selection. A diversity of approaches to expansions in Europe is yielding insights into the histories of Phoenicians, Roma, Anglo-Saxons and Vikings, and new methods for producing and analysing genome-wide data hold much promise. The field would benefit from more consensus on appropriate methods, and better communication between geneticists and experts in other disciplines, such as history, archaeology and linguistics.

Kinship Analysis with Diallelic SNPs - Experiences with the SNPforID Multiplex in an ISO17025 Accreditated Laboratory

BACKGROUND: The mutation rate of single nucleotide polymorphisms (SNPs) is estimated to be 100,000 times lower than that of short tandem repeats (STRs), which makes SNPs very suitable for relationship testing. The SNPforID multiplex assay was the first SNP typing assay that was a real alternative to the commonly used STR kits in kinship and crime case work and the first SNP assay to be validated in a forensic laboratory accredited according to the ISO17025 standard. METHODS: A total of 54 crime case samples were typed with the SNPforID multiplex assay. 30 samples from relationship cases were sequenced in selected SNP loci. RESULTS: It was demonstrated that mixtures were easily detected with the SNPforID assay by analyzing the signal strengths of the detected alleles. Unusual imbalances in signal strengths that were observed in a few individuals could be explained by unexpected SNPs in one of the primer binding sites. A complicated relationship case with four closely related individuals is presented. CONCLUSION: Mixtures can be detected with bi-allelic SNPs. The SNPforID assay is a very useful supplement to the STR kits in relationship testing.

lobSTR: A short tandem repeat profiler for personal genomes

Short tandem repeats (STRs) have a wide range of applications, including medical genetics, forensics, and genetic genealogy. High-throughput sequencing (HTS) has the potential to profile hundreds of thousands of STR loci. However, mainstream bioinformatics pipelines are inadequate for the task. These pipelines treat STR mapping as gapped alignment, which results in cumbersome processing times and a biased sampling of STR alleles. Here, we present lobSTR, a novel method for profiling STRs in personal genomes. lobSTR harnesses concepts from signal processing and statistical learning to avoid gapped alignment and to address the specific noise patterns in STR calling. The speed and reliability of lobSTR exceed the performance of current mainstream algorithms for STR profiling. We validated lobSTR's accuracy by measuring its consistency in calling STRs from whole-genome sequencing of two biological replicates from the same individual, by tracing Mendelian inheritance patterns in STR alleles in whole-genome sequencing of a HapMap trio, and by comparing lobSTR results to traditional molecular techniques. Encouraged by the speed and accuracy of lobSTR, we used the algorithm to conduct a comprehensive survey of STR variations in a deeply sequenced personal genome. We traced the mutation dynamics of close to 100,000 STR loci and observed more than 50,000 STR variations in a single genome. lobSTR's implementation is an end-to-end solution. The package accepts raw sequencing reads and provides the user with the genotyping results. It is written in C/C++, includes multi-threading capabilities, and is compatible with the BAM format.

Y chromosome lineages in men of west African descent

The early African experience in the Americas is marked by the transatlantic slave trade from ∼1619 to 1850 and the rise of the plantation system. The origins of enslaved Africans were largely dependent on European preferences as well as the availability of potential laborers within Africa. Rice production was a key industry of many colonial South Carolina low country plantations. Accordingly, rice plantations owners within South Carolina often requested enslaved Africans from the so-called "Grain Coast" of western Africa (Senegal to Sierra Leone). Studies on the African origins of the enslaved within other regions of the Americas have been limited. To address the issue of origins of people of African descent within the Americas and understand more about the genetic heterogeneity present within Africa and the African Diaspora, we typed Y chromosome specific markers in 1,319 men consisting of 508 west and central Africans (from 12 populations), 188 Caribbeans (from 2 islands), 532 African Americans (AAs from Washington, DC and Columbia, SC), and 91 European Americans. Principal component and admixture analyses provide support for significant Grain Coast ancestry among African American men in South Carolina. AA men from DC and the Caribbean showed a closer affinity to populations from the Bight of Biafra. Furthermore, 30-40% of the paternal lineages in African descent populations in the Americas are of European ancestry. Diverse west African ancestries and sex-biased gene flow from EAs has contributed greatly to the genetic heterogeneity of African populations throughout the Americas and has significant implications for gene mapping efforts in these populations.

The dual origin of Tati-speakers from Dagestan as written in the genealogy of uniparental variants

OBJECTIVES

Tat language is classified in an Iranian subbranch of the Indo-European family. It is spoken in the Caucasus and in the West Caspian region by populations with heterogeneous cultural traditions and religion whose ancestry is unknown. The aim of this study is to get a first insight about the genetic history of this peculiar linguistic group.

METHODS

We investigated the uniparental gene pools, defined by NRY and mtDNA high-resolution markers, in two Tati-speaking communities from Dagestan: Mountain Jews or Juhur, who speak the Judeo-Tat dialect, and the Tats, who speak the Muslim-Tat dialect. The samples have been collected in monoethnic rural villages and selected on the basis of genealogical relationships. A novel approach aimed at resolving cryptic cases in the recent history of human populations, which combines the properties of uniparental genetic markers with the potential of "forward-in-time" computer simulations, is presented.

RESULTS

Judeo-Tats emerged as a group with tight matrilineal genetic legacy who separated early from other Jewish communities. Tats exhibited genetic signals of a much longer in situ evolution, which appear as substantially unlinked with other Indo-Iranian enclaves in the Caucasus.

CONCLUSIONS

The independent demographic histories of the two samples, with mutually reversed profiles at paternally and maternally transmitted genetic systems, suggest that geographic proximity and linguistic assimilation of Tati-speakers from Dagestan do not reflect a common ancestry.

Next-generation sequencing technologies and applications for human genetic history and forensics

Rapid advances in the development of sequencing technologies in recent years have enabled an increasing number of applications in biology and medicine. Here, we review key technical aspects of the preparation of DNA templates for sequencing, the biochemical reaction principles and assay formats underlying next-generation sequencing systems, methods for imaging and base calling, quality control, and bioinformatic approaches for sequence alignment, variant calling and assembly. We also discuss some of the most important advances that the new sequencing technologies have brought to the fields of human population genetics, human genetic history and forensic genetics.

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.

Time-dependent rates of molecular evolution

For over half a century, it has been known that the rate of morphological evolution appears to vary with the time frame of measurement. Rates of microevolutionary change, measured between successive generations, were found to be far higher than rates of macroevolutionary change inferred from the fossil record. More recently, it has been suggested that rates of molecular evolution are also time dependent, with the estimated rate depending on the timescale of measurement. This followed surprising observations that estimates of mutation rates, obtained in studies of pedigrees and laboratory mutation-accumulation lines, exceeded long-term substitution rates by an order of magnitude or more. Although a range of studies have provided evidence for such a pattern, the hypothesis remains relatively contentious. Furthermore, there is ongoing discussion about the factors that can cause molecular rate estimates to be dependent on time. Here we present an overview of our current understanding of time-dependent rates. We provide a summary of the evidence for time-dependent rates in animals, bacteria and viruses. We review the various biological and methodological factors that can cause rates to be time dependent, including the effects of natural selection, calibration errors, model misspecification and other artefacts. We also describe the challenges in calibrating estimates of molecular rates, particularly on the intermediate timescales that are critical for an accurate characterization of time-dependent rates. This has important consequences for the use of molecular-clock methods to estimate timescales of recent evolutionary events.

Improving human forensics through advances in genetics, genomics and molecular biology

Forensic DNA profiling currently allows the identification of persons already known to investigating authorities. Recent advances have produced new types of genetic markers with the potential to overcome some important limitations of current DNA profiling methods. Moreover, other developments are enabling completely new kinds of forensically relevant information to be extracted from biological samples. These include new molecular approaches for finding individuals previously unknown to investigators, and new molecular methods to support links between forensic sample donors and criminal acts. Such advances in genetics, genomics and molecular biology are likely to improve human forensic case work in the near future.