The population history of northeastern Siberia since the Pleistocene.

Northeastern Siberia has been inhabited by humans for more than 40,000 years but its deep population history remains poorly understood. Here we investigate the late Pleistocene population history of northeastern Siberia through analyses of 34 newly recovered ancient genomes that date to between 31,000 and 600 years ago. We document complex population dynamics during this period, including at least three major migration events: an initial peopling by a previously unknown Palaeolithic population of 'Ancient North Siberians' who are distantly related to early West Eurasian hunter-gatherers; the arrival of East Asian-related peoples, which gave rise to 'Ancient Palaeo-Siberians' who are closely related to contemporary communities from far-northeastern Siberia (such as the Koryaks), as well as Native Americans; and a Holocene migration of other East Asian-related peoples, who we name 'Neo-Siberians', and from whom many contemporary Siberians are descended. Each of these population expansions largely replaced the earlier inhabitants, and ultimately generated the mosaic genetic make-up of contemporary peoples who inhabit a vast area across northern Eurasia and the Americas.

Unmasking the tissue-resident eukaryotic DNA virome in humans.

Little is known on the landscape of viruses that reside within our cells, nor on the interplay with the host imperative for their persistence. Yet, a lifetime of interactions conceivably have an imprint on our physiology and immune phenotype. In this work, we revealed the genetic make-up and unique composition of the known eukaryotic human DNA virome in nine organs (colon, liver, lung, heart, brain, kidney, skin, blood, hair) of 31 Finnish individuals. By integration of quantitative (qPCR) and qualitative (hybrid-capture sequencing) analysis, we identified the DNAs of 17 species, primarily herpes-, parvo-, papilloma- and anello-viruses (>80% prevalence), typically persisting in low copies (mean 540 copies/ million cells). We assembled in total 70 viral genomes (>90% breadth coverage), distinct in each of the individuals, and identified high sequence homology across the organs. Moreover, we detected variations in virome composition in two individuals with underlying malignant conditions. Our findings reveal unprecedented prevalences of viral DNAs in human organs and provide a fundamental ground for the investigation of disease correlates. Our results from post-mortem tissues call for investigation of the crosstalk between human DNA viruses, the host, and other microbes, as it predictably has a significant impact on our health.

Current status of undergraduate teaching in forensic & legal medicine in Europe.

The European Council of Legal Medicine (ECLM) is the body established in 1992 to represent practitioners forensic & legal medicine and is composed of delegates of the countries of the European Union (EU) and from other countries which form part of Europe to a current total of 34 member countries. The aims of this study were to determine the current status of undergraduate forensic & legal medicine teaching in the curriculum of medical studies in ECLM countries and to use the results of this study to determine whether it would be appropriate to develop new guidelines and standards for harmonising the content of undergraduate forensic medicine training across ECLM member countries. A detailed questionnaire was sent to all individuals or organisations listed on the ECLM contact database. Responses were received from 21 of 33 countries on the database. These responses showed considerable emphasis on undergraduate teaching of forensic medicine in all countries with the exception of Belgium and the United Kingdom. There was great general consistency in the subjects taught. The data from this survey provide a baseline which should assist in developing a strategy to harmonise forensic & legal medicine undergraduate training in member countries of the ECLM. The ECLM is now in a good position to establish a pan-European working group to coordinate a consensus document identifying an appropriate and modern core undergraduate forensic medicine curriculum that can be presented to the medical education authorities in each country, and which can be adapted for local requirements, based on available personnel, the forensic medicine structure in the country, and most importantly, the needs of the local population.

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.

A cost-benefit analysis for use of large SNP panels and high throughput typing for forensic investigative genetic genealogy.

Next-generation sequencing (NGS), also known as massively sequencing, enables large dense SNP panel analyses which generate the genetic component of forensic investigative genetic genealogy (FIGG). While the costs of implementing large SNP panel analyses into the laboratory system may seem high and daunting, the benefits of the technology may more than justify the investment. To determine if an infrastructural investment in public laboratories and using large SNP panel analyses would reap substantial benefits to society, a cost-benefit analysis (CBA) was performed. This CBA applied the logic that an increase of DNA profile uploads to a DNA database due to a sheer increase in number of markers and a greater sensitivity of detection afforded with NGS and a higher hit/association rate due to large SNP/kinship resolution and genealogy will increase investigative leads, will be more effective for identifying recidivists which in turn reduces future victims of crime, and will bring greater safety and security to communities. Analyses were performed for worst case/best case scenarios as well as by simulation sampling the range spaces with multiple input values simultaneously to generate best estimate summary statistics. This study shows that the benefits, both tangible and intangible, over the lifetime of an advanced database system would be huge and can be projected to be for less than $1 billion per year (over a 10-year period) investment can reap on average > $4.8 billion in tangible and intangible cost-benefits per year. More importantly, on average > 50,000 individuals need not become victims if FIGG were employed, assuming investigative associations generated were acted upon. The benefit to society is immense making the laboratory investment a nominal cost. The benefits likely are underestimated herein. There is latitude in the estimated costs, and even if they were doubled or tripled, there would still be substantial benefits gained with a FIGG-based approach. While the data used in this CBA are US centric (primarily because data were readily accessible), the model is generalizable and could be used by other jurisdictions to perform relevant and representative CBAs.

Forensic neuropathology in the past decade: a scoping literature review.

While there has been notable research activity in the field of clinical neuropathology over the recent years, forensic approaches have been less frequent. This scoping literature review explored original research on forensic neuropathology over the past decade (January 1, 2010, until February 12, 2022) using the MEDLINE database. The aims were to (1) analyze the volume of research on the topic, (2) describe meta-level attributes and sample characteristics, and (3) summarize key research themes and methods. Of 5053 initial hits, 2864 fell within the target timeframe, and 122 were included in the review. Only 3-17 articles were published per year globally. Most articles originated from the Europe (39.3%) and Asia (36.1%) and were published in forensic journals (57.4%). A median sample included 57 subjects aged between 16 and 80 years. The most common research theme was traumatic intracranial injury (24.6%), followed by anatomy (12.3%) and substance abuse (11.5%). Key methods included immunotechniques (31.1%) and macroscopic observation (21.3%). Although a number of novel findings were reported, most were of preliminary nature and will require further validation. In order to reach breakthroughs and validate novel tools for routine use, more research input is urged from researchers across the world. It would be necessary to ensure appropriate sample sizes and make use of control groups.

A prospective cost-benefit analysis for nylon 4N6FLOQSwabs®: example of the process and potential benefits.

Laboratories and their criminal justice systems are confronted with challenges for implementing new technologies, practices, and policies even when there appears to be demonstrative benefits to operational performance. Impacting decisions are the often higher costs associated with, for example, new technologies, limited current budgets, and making hard decisions on what to sacrifice to take on the seemingly better approach. A prospective cost-benefit analysis (CBA) could help an agency better formulate its strategies and plans and more importantly delineate how a relatively small increase to take on, for example, a new technology can have large impact on the system (e.g., the agency, other agencies, victims and families, and taxpayers). To demonstrate the process and potential value a CBA was performed on the use of an alternate and more expensive swab with reported better DNA yield and being certified human DNA free (i.e., nylon 4N6FLOQSwabs®), versus the traditional less costly swab (i.e., cotton swab). Assumptions are described, potential underestimates and overestimates noted, different values applied (for low and modest to high), and potential benefits (monetary and qualitative) presented. The overall outcome is that the cost of using the more expensive technology pales compared with the potential tangible and intangible benefits. This approach could be a guide for laboratories (and associated criminal justice systems) worldwide to support increased funding, although the costs and benefits may vary locally and for different technologies, practices, and policies. With well-developed CBAs, goals of providing the best services to support the criminal justice system and society can be attained.

A genome-wide association study of tramadol metabolism from post-mortem samples.

Phase I tramadol metabolism requires cytochrome p450 family 2, subfamily D, polypeptide 6 (CYP2D6) to form O-desmethyltramadol (M1). CYP2D6 genetic variants may infer metabolizer phenotype; however, drug ADME (absorption, distribution, metabolism, and excretion) and response depend on protein pathway(s), not CYP2D6 alone. There is a paucity of data regarding the contribution of trans-acting proteins to idiosyncratic phenotypes following drug exposure. A genome-wide association study identified five markers (rs79983226/kgp11274252, rs9384825, rs62435418/kgp10370907, rs72732317/kgp3743668, and rs184199168/exm1592932) associated with the conversion of tramadol to M1 (M1:T). These SNPs reside within five genes previously implicated with adverse reactions. Analysis of accompanying toxicological meta-data revealed a significant positive linear relationship between M1:T and degree of sample polypharmacy. Taken together, these data identify candidate loci for potential clinical inferences of phenotype following exposure to tramadol and highlight sample polypharmacy as a possible diagnostic covariate in post-mortem genetic studies.

Completed suicides of citalopram users-the role of CYP genotypes and adverse drug interactions.

Depression is known to be a risk factor for suicide. Currently, the most used antidepressants are selective serotonin reuptake inhibitors (SSRIs). Not all users, however, benefit from them. In such cases, treatment failure can be explained in part by genetic differences. In this study, we investigated the role of pharmacogenetic factors in citalopram-positive completed suicides (n = 349). Since citalopram is metabolized by CYP2C19 and CYP2D6 enzymes, the study population was genotyped for clinically relevant CYP2C19 and CYP2D6 polymorphisms and CYP2D6 copy number variation. To assess genetic differences between suicide cases and Finns in general, Finnish population samples (n = 855) were used as controls. Also, the role of drug interactions among suicide cases was evaluated. We found enrichment of a combined group of genetically predicted poor and ultrarapid metabolizer phenotypes (gMPs) of CYP2C19 among suicide victims compared to controls 0.356 [0.31-0.41] vs. 0.265 [0.24-0.30] (p = 0.0065). In CYP2D6 gMPs, there was no difference between cases and controls when the study population was analyzed as a whole. However, there were significantly more poor metabolizers among females who committed suicide by poisoning compared to female controls. In 8% of all drug poisoning deaths, lifetime drug-drug interaction was evaluated having a contribution to the fatal outcome. From clinical perspective, pharmacogenetic testing prior to initiation of SSRI drug could be beneficial. It may also be useful in medico-legal settings as it may elucidate obscure poisoning cases. Also, the possibility of unintentional drug interactions should be taken into account in drug poisoning deaths.

Supervised Classification of CYP2D6 Genotype and Metabolizer Phenotype With Postmortem Tramadol-Exposed Finns.

Cytochrome p450 family 2, subfamily D, polypeptide 6 (CYP2D6) may be used to infer the metabolizer phenotype (MP) of an individual as poor, intermediate, extensive/normal, or ultrarapid. Metabolizer phenotypes may suggest idiosyncratic drug responses as contributing factors to cause and/or manner of death in postmortem investigations. Application of CYP2D6 has used long-range amplification of the locus and restriction enzyme digestion to detect single-nucleotide variants (SNVs) associated with MPs. This process can be cumbersome and requires knowledge of genotype phase. Phase may be achieved using long-read DNA sequencing and/or computational methods; however, both can be error prone, which may make it difficult or impractical for implementation into medicolegal practice. CYP2D6 was interrogated in postmortem autopsied Finns using supervised machine learning and feature selection to identify SNVs indicative of MP and/or rate of tramadol O-demethylation (T:M1). A subset of 18 CYP2D6 SNVs could predict MP/T:M1 with up to 96.3% accuracy given phased data. These data indicate that phase contributes to classification accuracy when using CYP2D6 data. Of these 18 SNVs, 3 are novel loci putatively associated with T:M1. These findings may enable design of small multiplexes for easy forensic application of MP prediction when cause and/or manner of death is unknown.

Post-mortem analysis of suicide victims shows ABCB1 haplotype 1236T-2677T-3435T as a candidate predisposing factor behind adverse drug reactions in females.

BACKGROUND: Genetic variation in efflux transporter, permeability glycoprotein (P-gp), has recently been associated with completed violent suicides and also violent suicide attempts. As depression is known to be a risk factor for suicide and many antidepressants are P-gp substrates, it has been speculated that inadequate antidepressant treatment response or adverse side effects could be involved. OBJECTIVES: The aim of this study was to investigate whether there is an association between the P-gp coding ABCB1 gene and completed suicides in citalopram users. Also, the effect of sex and suicide method used (violent vs. non-violent) was evaluated. MATERIALS AND METHODS: All cases included in the study population, 349 completed suicide victims and 284 controls, were shown to be positive for antidepressant citalopram in a post-mortem toxicological drug screen. ABCB1 1236C>T, 2677G>T/A and 3435C>T polymorphisms were determined by TaqMan genotyping assays. Haplotypes were constructed from genotype data using the PHASE software. The association between the manner of death and the ABCB1 haplotype was tested with logistic regression analysis. RESULTS: No statistically significant differences were observed in the ABCB1 allele or genotype frequencies between the suicide and control groups. However, the ABCB1 1236T-2677T-3435T haplotype was associated with completed suicides of female citalopram users (odds ratio: 2.23; 95% confidence interval: 1.22-4.07; P=0.009). After stratification by the method used for suicide, the association emerged in fatal intoxications (odds ratio: 2.51; 95% confidence interval: 1.29-4.87; P=0.007). In other groups, no statistically significant associations were observed. CONCLUSION: Our results suggest that female citalopram users with ABCB1 1236T-2677T-3435T are more vulnerable to adverse effects of the drugs as this haplotype was enriched in non-violent suicides of female citalopram users. Even though the biological mechanism behind this observation is unknown, the results provide another example of the importance of sex-based segregation in pharmacogenetics studies.

Monocytes accumulate in the airways of children with fatal asthma.

BACKGROUND: Activated T helper type 2 (Th2) cells are believed to play a pivotal role in allergic airway inflammation, but which cells attract and activate Th2 cells locally have not been fully determined. Recently, it was shown in an experimental human model of allergic rhinitis (AR) that activated monocytes rapidly accumulate in the nasal mucosa after local allergen challenge, where they promote recruitment of Th2 cells and eosinophils. OBJECTIVE: To investigate whether monocytes are recruited to the lungs in paediatric asthma. METHODS: Tissue samples obtained from children and adolescents with fatal asthma attack (n = 12), age-matched non-atopic controls (n = 9) and allergen-challenged AR patients (n = 8) were subjected to in situ immunostaining. RESULTS: Monocytes, identified as CD68+S100A8/A9+ cells, were significantly increased in the lower airway mucosa and in the alveoli of fatal asthma patients compared with control individuals. Interestingly, cellular aggregates containing CD68+S100A8/A9+ monocytes obstructing the lumen of bronchioles were found in asthmatics (8 out of 12) but not in controls. Analysing tissue specimens from challenged AR patients, we confirmed that co-staining with CD68 and S100A8/A9 was a valid method to identify recently recruited monocytes. We also showed that the vast majority of accumulating monocytes both in the lungs and in the nasal mucosa expressed matrix metalloproteinase 10, suggesting that this protein may be involved in their migration within the tissue. CONCLUSIONS AND CLINICAL RELEVANCE: Monocytes accumulated in the lungs of children and adolescents with fatal asthma attack. This finding strongly suggests that monocytes are directly involved in the immunopathology of asthma and that these pro-inflammatory cells are potential targets for therapy.

Challenges in investigation of diabetes-related aviation fatalities-an analysis of 1491 subsequent aviation fatalities in USA during 2011-2016.

Diabetes mellitus (DM) could cause pilot incapacitation and result in aviation fatalities. The mechanisms could be directly as a consequence of acute hypoglycemia/subacute diabetic ketoacidosis (DKA) or indirectly as an acute cardiovascular event by contributing to the development of atherosclerosis in coronary or carotid and cerebral arteries. In this study, DM-related fatal flight accidents in the US National Transport Bureau's database between years 2011-2016 were analyzed with special emphasis on postmortem (PM) glucose levels and correlation of toxicological reports with anamnestic information on DM. Additionally, autopsy results on coronary arteries were reviewed. In 43 out of 1491 (~ 3%) fatal accidents pilots had DM. Postmortem glucose or glycated hemoglobin percentage (Hb1Ac) was measured in 12 of the 43 cases; while antidiabetic medication was found in 14 of the cases (only two of the cases had both glucose measurements and medication). With the increasing prevalence of DM, a possibility of pilot incapacitation due to DM or complications of DM should be actively studied, even if no anamnestic information of DM was available. While PM hypoglycemia is difficult to assess, we propose a systematic investigation based on measurement of glucose, Hb1Ac%, and ketone bodies, and documentation of atherosclerotic lesions in major arteries to identify or rule out DM as a cause of pilot incapacitation.

Improved Y-STR typing for disaster victim identification, missing persons investigations, and historical human skeletal remains.

Bones are a valuable source of DNA in forensic, anthropological, and archaeological investigations. There are a number of scenarios in which the only samples available for testing are highly degraded and/or skeletonized. Often it is necessary to perform more than one type of marker analysis on such samples in order to compile sufficient data for identification. Lineage markers, such as Y-STRs and mitochondrial DNA (mtDNA), represent important systems to complement autosomal DNA markers and anthropological metadata in making associations between unidentified remains and living relatives or for characterization of the remains for historical and archaeological studies. In this comparative study, Y-STR typing with both Yfiler™ and Yfiler™ Plus (Thermo Fisher Scientific, Waltham, MA, USA) was performed on a variety of human skeletal remains, including samples from the American Civil War (1861-1865), the late nineteenth century gold rush era in Deadwood, SD, USA (1874-1877), the Seven Years' War (1756-1763), a seventeenth-century archaeological site in Raspenava, Bohemia (Czech Republic), and World War II (1939-1945). The skeletal remains used for this study were recovered from a wide range of environmental conditions and were extracted using several common methods. Regardless of the DNA extraction method used and the age/condition of the remains, 22 out of 24 bone samples yielded a greater number of alleles using the Yfiler™ Plus kit compared to the Yfiler™ kit using the same quantity of input DNA. There was no discernable correlation with the degradation index values for these samples. Overall, the efficacy of the Yfiler™ Plus assay was demonstrated on degraded DNA from skeletal remains. Yfiler™ Plus increases the discriminatory power over the previous generation multiplex due to the larger set of Y-STR markers available for analysis and buffer modifications with the newer version kit. Increased haplotype resolution is provided to infer or refute putative genetic relationships.

Full-gene haplotypes refine CYP2D6 metabolizer phenotype inferences.

CYP2D6 is a critical pharmacogenetic target, and polymorphisms in the gene region are commonly used to infer enzyme activity score and predict resulting metabolizer phenotype: poor, intermediate, extensive/normal, or ultrarapid which can be useful in determining cause and/or manner of death in some autopsies. Current genotyping approaches are incapable of identifying novel and/or rare variants, so CYP2D6 star allele definitions are limited to polymorphisms known a priori. While useful for most predictions, recent studies using massively parallel sequencing data have identified additional polymorphisms in CYP2D6 that are predicted to alter enzyme function but are not considered in current star allele nomenclature. The 1000 Genomes Project data were used to produce full-gene haplotypes, describe their distribution in super-populations, and predict enzyme activity scores. Full-gene haplotypes generated lower activity scores than current approaches due to inclusion of additional damaging polymorphisms in the star allele. These findings are critical for clinical implementation of metabolizer phenotype prediction because a fraction of the population may be incorrectly considered normal metabolizers but actually may be poor or intermediate metabolizers.

Immunohistology and remodeling in fatal pediatric and adolescent asthma.

BACKGROUND: Thickening of reticular basement membrane, increased airway smooth muscle mass and eosinophilic inflammation are found in adult fatal asthma. At the present study the histopathology of fatal paediatric and adolescent asthma is evaluated. METHODS: Post-mortem lung autopsies from 12 fatal asthma cases and 8 non-asthmatic control subjects were examined. Thickness of reticular basement membrane (RBM) and percentage of airway smooth muscle (ASM%) mass area were measured and inflammatory cells were counted. Patient records were reviewed for clinical history. RESULTS: The age range of the cases was from 0.9 to 19.5 years, eight were males and five had received inhaled corticosteroids. Thickened RBM was detected in majority of the cases without any correlation to treatment delay, age at onset of symptoms or diagnosis. In the large airways ASM was clearly increased in one third of the cases whereas the median ASM% did not differ from that in healthy controls (14.0% vs. 14.0%). In small airways no increase of ASM was found, instead mucous plugs were seen in fatal asthma. The number of eosinophils, plasmacytoid dendritic cells, macrophages, and B-cells were significantly increased in fatal asthma cases compared with controls and the two latter correlated with the length of the fatal exacerbation. CONCLUSIONS: The findings highlight the strong presence of eosinophils and mucous plugs even in small airways in children and adolescents with fatal asthma. Thickened RBM was obvious in majority of the patients. Contrary to our hypothesis, increased ASM% was detected in only one third of the patients.

Glucose metabolism in completed suicide: a forensic-pathological pilot study.

AIM: To determine whether antemortem blood levels of glycated hemoglobin (HbA1c) and glucose predict completed suicide and, by extension, whether markers of glucose metabolism might be associated with a prosuicidal trait or state. METHOD: From consecutively performed autopsies, samples of blood and vitreous humor from 17 suicide victims and 27 non-suicide controls were compared with regard to levels of glucose, lactate, and HbA1c. RESULTS: Mean HbA1c was higher, and mean estimated blood glucose lower, among suicide victims, although tests revealed no significant differences (P=0.171 and P=0.395, respectively). HbA1c levels exceeding 48.0 mmol/mol, which were indicative of persistent hyperglycemia, were twice as common in suicide victims (59% vs 30%; P=0.068). CONCLUSION: The finding of this pilot study suggest that deranged glucose metabolism may reflect biological events antecedent to, or concomitant with, completed suicide, with the following clinical implications: recurring hyperglycemia due to defective glucose transport, which may give rise to depression and suicidal ideation, and elevated HbA1c levels, which may represent an assayable correlate to neurobiological conditions predisposing to suicide.

Serotonergic 5HTTLPR/rs25531 s-allele homozygosity associates with violent suicides in male citalopram users.

Depressive disorders are involved as a background factor in over 50% of suicide cases. The most widely used antidepressants today are serotonin selective reuptake inhibitors (SSRIs). However, not all users benefit from SSRI medication. Although the overall number of suicides in Finland have decreased notably during the last decade, the annual rate is still relatively high, particularly in male population. In this study, we tested the hypothesis that the genetic variants associated with decreased citalopram efficiency, 5HTTLPR/rs25531, and increased impulsive behavior, MAOA-uVNTR and HTR2B Q20*, are more frequent among citalopram users committing suicide than among the citalopram users in general. Also the effect of alcohol was evaluated. The study population comprised 349 suicide victims (184 males and 165 females). Based on the suicide method used, cases were divided into two groups; violent (88 males and 49 females) and non-violent (96 males and 116 females). The control group (284; 159 males and 125 females) consisted of citalopram users who died of causes other than suicide. We found that male citalopram users with low functioning s/s genotype of 5HTTLPR/rs25531 were in increased risk to commit violent suicide (OR 2.50, 95%CI 1.15-5.42, p = 0.020). Surprisingly, high blood alcohol concentration was observed to be a risk factor only in non-violent suicides (both males and females), but not in violent ones. No association between suicides and MAOA-uVNTR and HTR2B Q20*, which have been previously connected to violent and impulsive behavior, was detected.

The Y-chromosome tree bursts into leaf: 13,000 high-confidence SNPs covering the majority of known clades.

Many studies of human populations have used the male-specific region of the Y chromosome (MSY) as a marker, but MSY sequence variants have traditionally been subject to ascertainment bias. Also, dating of haplogroups has relied on Y-specific short tandem repeats (STRs), involving problems of mutation rate choice, and possible long-term mutation saturation. Next-generation sequencing can ascertain single nucleotide polymorphisms (SNPs) in an unbiased way, leading to phylogenies in which branch-lengths are proportional to time, and allowing the times-to-most-recent-common-ancestor (TMRCAs) of nodes to be estimated directly. Here we describe the sequencing of 3.7 Mb of MSY in each of 448 human males at a mean coverage of 51×, yielding 13,261 high-confidence SNPs, 65.9% of which are previously unreported. The resulting phylogeny covers the majority of the known clades, provides date estimates of nodes, and constitutes a robust evolutionary framework for analyzing the history of other classes of mutation. Different clades within the tree show subtle but significant differences in branch lengths to the root. We also apply a set of 23 Y-STRs to the same samples, allowing SNP- and STR-based diversity and TMRCA estimates to be systematically compared. Ongoing purifying selection is suggested by our analysis of the phylogenetic distribution of nonsynonymous variants in 15 MSY single-copy genes.