The genome-wide structure of the Jewish people

Founder mutation for Huntington disease in Caucasus Jews

Huntington disease (HD), an autosomal dominant disorder involving HTT, is characterized by chorea, psychiatric illness and cognitive decline. Diagnosis and age of onset depend on the degree of expansion of the trinucleotide CAG repeat within the gene. The prevalence of HD is known for Europeans but has not been studied in the Israeli population. Between 2006 and 2011 we diagnosed in our adult genetics clinic ten HD probands, nine of whom were Caucasus Jews (CJ) (Azerbaijani), and one Ashkenazi Jewish. We performed haplotype analysis to look for evidence of a founder mutation, and found that of the nine CJ, eight shared the same haplotype that was compatible with the A1 haplogroup. We calculated the coalescence age of the mutation to be between 80 and 150 years. Ninety percent of our HD patients are CJ, as are 27% of the HD patients in Israel, although the CJ comprise only 1.4% of the Israeli population. Our findings suggest a higher prevalence of HD among CJ compared to the general Israeli population and are consistent with a recent founder mutation. We recommend a higher degree of suspicion for HD in CJ with subtle clinical findings.

Next Generation Statistical Genetics: Modeling, Penalization, and Optimization in High-Dimensional Data

Statistical genetics is undergoing the same transition to big data that all branches of applied statistics are experiencing. With the advent of inexpensive DNA sequencing, the transition is only accelerating. This brief review highlights some modern techniques with recent successes in statistical genetics. These include: (a) lasso penalized regression and association mapping, (b) ethnic admixture estimation, (c) matrix completion for genotype and sequence data, (d) the fused lasso and copy number variation, (e) haplotyping, (f) estimation of relatedness, (g) variance components models, and (h) rare variant testing. For more than a century, genetics has been both a driver and beneficiary of statistical theory and practice. This symbiotic relationship will persist for the foreseeable future.

The GenoChip: a new tool for genetic anthropology

The Genographic Project is an international effort aimed at charting human migratory history. The project is nonprofit and nonmedical, and, through its Legacy Fund, supports locally led efforts to preserve indigenous and traditional cultures. Although the first phase of the project was focused on uniparentally inherited markers on the Y-chromosome and mitochondrial DNA (mtDNA), the current phase focuses on markers from across the entire genome to obtain a more complete understanding of human genetic variation. Although many commercial arrays exist for genome-wide single-nucleotide polymorphism (SNP) genotyping, they were designed for medical genetic studies and contain medically related markers that are inappropriate for global population genetic studies. GenoChip, the Genographic Project’s new genotyping array, was designed to resolve these issues and enable higher resolution research into outstanding questions in genetic anthropology. The GenoChip includes ancestry informative markers obtained for over 450 human populations, an ancient human (Saqqaq), and two archaic hominins (Neanderthal and Denisovan) and was designed to identify all known Y-chromosome and mtDNA haplogroups. The chip was carefully vetted to avoid inclusion of medically relevant markers. To demonstrate its capabilities, we compared the F_ST distributions of GenoChip SNPs to those of two commercial arrays. Although all arrays yielded similarly shaped (inverse J) F_ST distributions, the GenoChip autosomal and X-chromosomal distributions had the highest mean F_ST, attesting to its ability to discern subpopulations. The chip performances are illustrated in a principal component analysis for 14 worldwide populations. In summary, the GenoChip is a dedicated genotyping platform for genetic anthropology. With an unprecedented number of approximately 12,000 Y-chromosomal and approximately 3,300 mtDNA SNPs and over 130,000 autosomal and X-chromosomal SNPs without any known health, medical, or phenotypic relevance, the GenoChip is a useful tool for genetic anthropology and population genetics.

The light skin allele of SLC24A5 in South Asians and Europeans shares identity by descent

Skin pigmentation is one of the most variable phenotypic traits in humans. A non-synonymous substitution (rs1426654) in the third exon of SLC24A5 accounts for lighter skin in Europeans but not in East Asians. A previous genome-wide association study carried out in a heterogeneous sample of UK immigrants of South Asian descent suggested that this gene also contributes significantly to skin pigmentation variation among South Asians. In the present study, we have quantitatively assessed skin pigmentation for a largely homogeneous cohort of 1228 individuals from the Southern region of the Indian subcontinent. Our data confirm significant association of rs1426654 SNP with skin pigmentation, explaining about 27% of total phenotypic variation in the cohort studied. Our extensive survey of the polymorphism in 1573 individuals from 54 ethnic populations across the Indian subcontinent reveals wide presence of the derived-A allele, although the frequencies vary substantially among populations. We also show that the geospatial pattern of this allele is complex, but most importantly, reflects strong influence of language, geography and demographic history of the populations. Sequencing 11.74 kb of SLC24A5 in 95 individuals worldwide reveals that the rs1426654-A alleles in South Asian and West Eurasian populations are monophyletic and occur on the background of a common haplotype that is characterized by low genetic diversity. We date the coalescence of the light skin associated allele at 22–28 KYA. Both our sequence and genome-wide genotype data confirm that this gene has been a target for positive selection among Europeans. However, the latter also shows additional evidence of selection in populations of the Middle East, Central Asia, Pakistan and North India but not in South India.

Human skin color is one of the most visible aspects of human diversity. The genetic basis of pigmentation in Europeans has been understood to some extent, but our knowledge about South Asians has been restricted to a handful of studies. It has been suggested that a single nucleotide difference in SLC24A5 accounts for 25–38% European-African pigmentation differences and correlates with lighter skin. This genetic variant has also been associated with skin color variation among South Asians living in the UK. Here, we report a study based on a homogenous cohort of South India. Our results confirm that SLC24A5 plays a key role in pigmentation diversity of South Asians. Country-wide screening of the variant reveals that the light skin associated allele is widespread in the Indian subcontinent and its complex patterning is shaped by a combination of processes involving selection and demographic history of the populations. By studying the variation of SLC24A5 sequences among a diverse set of individuals, we show that the light skin associated allele in South Asians is identical by descent to that found in Europeans. Our study also provides new insights into positive selection acting on the gene and the evolutionary history of light skin in humans.

Molecular phylogeography of a human autosomal skin color locus under natural selection

Divergent natural selection caused by differences in solar exposure has resulted in distinctive variations in skin color between human populations. The derived light skin color allele of the SLC24A5 gene, A111T, predominates in populations of Western Eurasian ancestry. To gain insight into when and where this mutation arose, we defined common haplotypes in the genomic region around SLC24A5 across diverse human populations and deduced phylogenetic relationships between them. Virtually all chromosomes carrying the A111T allele share a single 78-kb haplotype that we call C11, indicating that all instances of this mutation in human populations share a common origin. The C11 haplotype was most likely created by a crossover between two haplotypes, followed by the A111T mutation. The two parental precursor haplotypes are found from East Asia to the Americas but are nearly absent in Africa. The distributions of C11 and its parental haplotypes make it most likely that these two last steps occurred between the Middle East and the Indian subcontinent, with the A111T mutation occurring after the split between the ancestors of Europeans and East Asians.

The Nature/Culture of Genetic Facts

This review aims to explore the relationship between anthropology and genetics, an intellectual zone that has been occupied in different ways over the past century. One way to think about it is to contrast a classical “anthropological genetics” ( Roberts 1965 ), that is to say, a genetics that presumably informs anthropological issues or questions, with a “genomic anthropology” ( Pálsson 2008 ), that is to say, an anthropology that complements and relativizes modern genomics (on the model of, say, medical anthropology and legal anthropology). 1 This review argues that a principal contribution of anthropology to the study of human heredity lies in the ontology of genetic facts. For anthropology, genetic facts are not natural, with meanings inscribed on them, but are instead natural/cultural: The natural facts have cultural information (values, ideologies, meanings) integrated into them, not layered on them. To understand genetic facts involves confronting their production, which has classically been restricted to questions of methodology but which may be conceptualized more broadly. This review is not intended as a critique of the field of anthropological genetics, but as a reformulation of its central objects of study. I argue for reconceptualizing the ontology of scientific facts in anthropological genetics, not as (value-neutral) biological facts situated in a cultural context, but instead as inherently biocultural facts.

A general linear model-based approach for inferring selection to climate

Background

Many efforts have been made to detect signatures of positive selection in the human genome, especially those associated with expansion from Africa and subsequent colonization of all other continents. However, most approaches have not directly probed the relationship between the environment and patterns of variation among humans. We have designed a method to identify regions of the genome under selection based on Mantel tests conducted within a general linear model framework, which we call MAntel-GLM to Infer Clinal Selection (MAGICS). MAGICS explicitly incorporates population-specific and genome-wide patterns of background variation as well as information from environmental values to provide an improved picture of selection and its underlying causes in human populations.

Results

Our results significantly overlap with those obtained by other published methodologies, but MAGICS has several advantages. These include improvements that: limit false positives by reducing the number of independent tests conducted and by correcting for geographic distance, which we found to be a major contributor to selection signals; yield absolute rather than relative estimates of significance; identify specific geographic regions linked most strongly to particular signals of selection; and detect recent balancing as well as directional selection.

Conclusions

We find evidence of selection associated with climate (P < 10^-5) in 354 genes, and among these observe a highly significant enrichment for directional positive selection. Two of our strongest 'hits’, however, ADRA2A and ADRA2C, implicated in vasoconstriction in response to cold and pain stimuli, show evidence of balancing selection. Our results clearly demonstrate evidence of climate-related signals of directional and balancing selection.

Global patterns of diversity and selection in human tyrosinase gene

Global variation in skin pigmentation is one of the most striking examples of environmental adaptation in humans. More than two hundred loci have been identified as candidate genes in model organisms and a few tens of these have been found to be significantly associated with human skin pigmentation in genome-wide association studies. However, the evolutionary history of different pigmentation genes is rather complex: some loci have been subjected to strong positive selection, while others evolved under the relaxation of functional constraints in low UV environment. Here we report the results of a global study of the human tyrosinase gene, which is one of the key enzymes in melanin production, to assess the role of its variation in the evolution of skin pigmentation differences among human populations. We observe a higher rate of non-synonymous polymorphisms in the European sample consistent with the relaxation of selective constraints. A similar pattern was previously observed in the MC1R gene and concurs with UV radiation-driven model of skin color evolution by which mutations leading to lower melanin levels and decreased photoprotection are subject to purifying selection at low latitudes while being tolerated or even favored at higher latitudes because they facilitate UV-dependent vitamin D production. Our coalescent date estimates suggest that the non-synonymous variants, which are frequent in Europe and North Africa, are recent and have emerged after the separation of East and West Eurasian populations.

The effect of very high versus very low sustained loading on the lower back and knees in middle life

To evaluate the effect of the extremes of long term high and low physical activities on musculoskeletal heath in middle age, a historical cohort study was performed. The MRI knee and back findings of 25 randomly selected subjects who were inducted into the armed forces in 1983 and served at least 3 years as elite infantry soldiers were compared 25 years later, with 20 randomly selected subjects who were deferred from army service for full time religious studies at the same time. Both cohorts were from the same common genome. The two primary outcome measures were degenerative lumbar disc disease evaluated by the Pfirrmann score and degenerative knee changes evaluated by the WORMS score. At the 25-year follow up, the mean Pfirrmann score (8.6) for the L1 to S1 level of the elite infantry group was significantly higher than that of the sedentary group (6.7), (P = 0.003). There was no statistically significant difference between the WORMS knee scores between the two cohorts (P = 0.7). In spite of the much greater musculoskeletal loading history of the elite infantry cohort, only their lumbar spines but not their knees showed increased degenerative changes at middle age by MRI criteria.

Genetic evidence for recent population mixture in India

Most Indian groups descend from a mixture of two genetically divergent populations: Ancestral North Indians (ANI) related to Central Asians, Middle Easterners, Caucasians, and Europeans; and Ancestral South Indians (ASI) not closely related to groups outside the subcontinent. The date of mixture is unknown but has implications for understanding Indian history. We report genome-wide data from 73 groups from the Indian subcontinent and analyze linkage disequilibrium to estimate ANI-ASI mixture dates ranging from about 1,900 to 4,200 years ago. In a subset of groups, 100% of the mixture is consistent with having occurred during this period. These results show that India experienced a demographic transformation several thousand years ago, from a region in which major population mixture was common to one in which mixture even between closely related groups became rare because of a shift to endogamy.

Biological races in humans

Races may exist in humans in a cultural sense, but biological concepts of race are needed to access their reality in a non-species-specific manner and to see if cultural categories correspond to biological categories within humans. Modern biological concepts of race can be implemented objectively with molecular genetic data through hypothesis-testing. Genetic data sets are used to see if biological races exist in humans and in our closest evolutionary relative, the chimpanzee. Using the two most commonly used biological concepts of race, chimpanzees are indeed subdivided into races but humans are not. Adaptive traits, such as skin color, have frequently been used to define races in humans, but such adaptive traits reflect the underlying environmental factor to which they are adaptive and not overall genetic differentiation, and different adaptive traits define discordant groups. There are no objective criteria for choosing one adaptive trait over another to define race. As a consequence, adaptive traits do not define races in humans. Much of the recent scientific literature on human evolution portrays human populations as separate branches on an evolutionary tree. A tree-like structure among humans has been falsified whenever tested, so this practice is scientifically indefensible. It is also socially irresponsible as these pictorial representations of human evolution have more impact on the general public than nuanced phrases in the text of a scientific paper. Humans have much genetic diversity, but the vast majority of this diversity reflects individual uniqueness and not race.

Mitogenomes from two uncommon haplogroups mark late glacial/postglacial expansions from the near east and neolithic dispersals within Europe

The current human mitochondrial (mtDNA) phylogeny does not equally represent all human populations but is biased in favour of representatives originally from north and central Europe. This especially affects the phylogeny of some uncommon West Eurasian haplogroups, including I and W, whose southern European and Near Eastern components are very poorly represented, suggesting that extensive hidden phylogenetic substructure remains to be uncovered. This study expanded and re-analysed the available datasets of I and W complete mtDNA genomes, reaching a comprehensive 419 mitogenomes, and searched for precise correlations between the ages and geographical distributions of their numerous newly identified subclades with events of human dispersal which contributed to the genetic formation of modern Europeans. Our results showed that haplogroups I (within N1a1b) and W originated in the Near East during the Last Glacial Maximum or pre-warming period (the period of gradual warming between the end of the LGM, ∼19 ky ago, and the beginning of the first main warming phase, ∼15 ky ago) and, like the much more common haplogroups J and T, may have been involved in Late Glacial expansions starting from the Near East. Thus our data contribute to a better definition of the Late and postglacial re-peopling of Europe, providing further evidence for the scenario that major population expansions started after the Last Glacial Maximum but before Neolithic times, but also evidencing traces of diffusion events in several I and W subclades dating to the European Neolithic and restricted to Europe.

Neanderthal and Denisova genetic affinities with contemporary humans: introgression versus common ancestral polymorphisms

Analyses of the genetic relationships among modern humans, Neanderthals and Denisovans have suggested that 1-4% of the non-Sub-Saharan African gene pool may be Neanderthal derived, while 6-8% of the Melanesian gene pool may be the product of admixture between the Denisovans and the direct ancestors of Melanesians. In the present study, we analyzed single nucleotide polymorphism (SNP) diversity among a worldwide collection of contemporary human populations with respect to the genetic constitution of these two archaic hominins and Pan troglodytes (chimpanzee). We partitioned SNPs into subsets, including those that are derived in both archaic lineages, those that are ancestral in both archaic lineages and those that are only derived in one archaic lineage. By doing this, we have conducted separate examinations of subsets of mutations with higher probabilities of divergent phylogenetic origins. While previous investigations have excluded SNPs from common ancestors in principal component analyses, we included common ancestral SNPs in our analyses to visualize the relative placement of the Neanderthal and Denisova among human populations. To assess the genetic similarities among the various hominin lineages, we performed genetic structure analyses to provide a comparison of genetic patterns found within contemporary human genomes that may have archaic or common ancestral roots. Our results indicate that 3.6% of the Neanderthal genome is shared with roughly 65.4% of the average European gene pool, which clinally diminishes with distance from Europe. Our results suggest that Neanderthal genetic associations with contemporary non-Sub-Saharan African populations, as well as the genetic affinities observed between Denisovans and Melanesians most likely result from the retention of ancient mutations in these populations.

CYP2B6 SNPs are associated with methadone dose required for effective treatment of opioid addiction

Adequate methadone dosing in methadone maintenance treatment (MMT) for opioid addiction is critical for therapeutic success. One of the challenges in dose determination is the inter-individual variability in dose-response. Methadone metabolism is attributed primarily to cytochrome P450 enzymes CYP3A4, CYP2B6 and CYP2D6. The CYP2B6*6 allele [single nucleotide polymorphisms (SNPs) 785A>G (rs2279343) and 516G>T (rs3745274)] was associated with slow methadone metabolism. To explore the effects of CYP2B6*6 allele on methadone dose requirement, it was genotyped in a well-characterized sample of 74 Israeli former heroin addicts in MMT. The sample is primarily of Middle Eastern/European ancestry, based on ancestry informative markers (AIMs). Only patients with no major co-medication that may affect methadone metabolism were included. The stabilizing daily methadone dose in this sample ranges between 13 and 260mg (mean 140±52mg). The mean methadone doses required by subjects homozygous for the variant alleles of the CYP2B6 SNPs 785A>G and 516G>T (88, 96mg, respectively) were significantly lower than those of the heterozygotes (133, 129mg, respectively) and the non-carriers (150, 151mg, respectively) (nominal P=0.012, 0.048, respectively). The results remain significant after controlling for age, sex and the ABCB1 SNP 1236C>T (rs1128503), which was previously shown to be associated with high methadone dose requirement in this population (P=0.006, 0.030, respectively). An additional 77 CYP2B6, CYP3A4 and CYP2D6 SNPs were genotyped. Of these, 24 SNPs were polymorphic and none showed significant association with methadone dose. Further studies are necessary to replicate these preliminary findings in additional subjects and other populations.

Variation at diabetes- and obesity-associated Loci may mirror neutral patterns of human population diversity and diabetes prevalence in India

South Asian populations harbor a high degree of genetic diversity, due in part to demographic history. Two studies on genome-wide variation in Indian populations have shown that most Indian populations show varying degrees of admixture between ancestral north Indian and ancestral south Indian components. As a result of this structure, genetic variation in India appears to follow a geographic cline. Similarly, Indian populations seem to show detectable differences in diabetes and obesity prevalence between different geographic regions of the country. We tested the hypothesis that genetic variation at diabetes- and obesity-associated loci may be potentially related to different genetic ancestries. We genotyped 2977 individuals from 61 populations across India for 18 SNPs in genes implicated in T2D and obesity. We examined patterns of variation in allele frequency across different geographical gradients and considered state of origin and language affiliation. Our results show that most of the 18 SNPs show no significant correlation with latitude, the geographic cline reported in previous studies, or by language family. Exceptions include KCNQ1 with latitude and THADA and JAK1 with language, which suggests that genetic variation at previously ascertained diabetes-associated loci may only partly mirror geographic patterns of genome-wide diversity in Indian populations.

Autosomal and uniparental portraits of the native populations of Sakha (Yakutia): implications for the peopling of Northeast Eurasia

BACKGROUND

Sakha--an area connecting South and Northeast Siberia--is significant for understanding the history of peopling of Northeast Eurasia and the Americas. Previous studies have shown a genetic contiguity between Siberia and East Asia and the key role of South Siberia in the colonization of Siberia.

RESULTS

We report the results of a high-resolution phylogenetic analysis of 701 mtDNAs and 318 Y chromosomes from five native populations of Sakha (Yakuts, Evenks, Evens, Yukaghirs and Dolgans) and of the analysis of more than 500,000 autosomal SNPs of 758 individuals from 55 populations, including 40 previously unpublished samples from Siberia. Phylogenetically terminal clades of East Asian mtDNA haplogroups C and D and Y-chromosome haplogroups N1c, N1b and C3, constituting the core of the gene pool of the native populations from Sakha, connect Sakha and South Siberia. Analysis of autosomal SNP data confirms the genetic continuity between Sakha and South Siberia. Maternal lineages D5a2a2, C4a1c, C4a2, C5b1b and the Yakut-specific STR sub-clade of Y-chromosome haplogroup N1c can be linked to a migration of Yakut ancestors, while the paternal lineage C3c was most likely carried to Sakha by the expansion of the Tungusic people. MtDNA haplogroups Z1a1b and Z1a3, present in Yukaghirs, Evens and Dolgans, show traces of different and probably more ancient migration(s). Analysis of both haploid loci and autosomal SNP data revealed only minor genetic components shared between Sakha and the extreme Northeast Siberia. Although the major part of West Eurasian maternal and paternal lineages in Sakha could originate from recent admixture with East Europeans, mtDNA haplogroups H8, H20a and HV1a1a, as well as Y-chromosome haplogroup J, more probably reflect an ancient gene flow from West Eurasia through Central Asia and South Siberia.

CONCLUSIONS

Our high-resolution phylogenetic dissection of mtDNA and Y-chromosome haplogroups as well as analysis of autosomal SNP data suggests that Sakha was colonized by repeated expansions from South Siberia with minor gene flow from the Lower Amur/Southern Okhotsk region and/or Kamchatka. The minor West Eurasian component in Sakha attests to both recent and ongoing admixture with East Europeans and an ancient gene flow from West Eurasia.

The missing link of Jewish European ancestry: contrasting the Rhineland and the Khazarian hypotheses

The question of Jewish ancestry has been the subject of controversy for over two centuries and has yet to be resolved. The “Rhineland hypothesis” depicts Eastern European Jews as a “population isolate” that emerged from a small group of German Jews who migrated eastward and expanded rapidly. Alternatively, the “Khazarian hypothesis” suggests that Eastern European Jews descended from the Khazars, an amalgam of Turkic clans that settled the Caucasus in the early centuries CE and converted to Judaism in the 8th century. Mesopotamian and Greco–Roman Jews continuously reinforced the Judaized empire until the 13th century. Following the collapse of their empire, the Judeo–Khazars fled to Eastern Europe. The rise of European Jewry is therefore explained by the contribution of the Judeo–Khazars. Thus far, however, the Khazars’ contribution has been estimated only empirically, as the absence of genome-wide data from Caucasus populations precluded testing the Khazarian hypothesis. Recent sequencing of modern Caucasus populations prompted us to revisit the Khazarian hypothesis and compare it with the Rhineland hypothesis. We applied a wide range of population genetic analyses to compare these two hypotheses. Our findings support the Khazarian hypothesis and portray the European Jewish genome as a mosaic of Near Eastern-Caucasus, European, and Semitic ancestries, thereby consolidating previous contradictory reports of Jewish ancestry. We further describe a major difference among Caucasus populations explained by the early presence of Judeans in the Southern and Central Caucasus. Our results have important implications for the demographic forces that shaped the genetic diversity in the Caucasus and for medical studies.

Population structure in a comprehensive genomic data set on human microsatellite variation

Over the past two decades, microsatellite genotypes have provided the data for landmark studies of human population-genetic variation. However, the various microsatellite data sets have been prepared with different procedures and sets of markers, so that it has been difficult to synthesize available data for a comprehensive analysis. Here, we combine eight human population-genetic data sets at the 645 microsatellite loci they share in common, accounting for procedural differences in the production of the different data sets, to assemble a single data set containing 5795 individuals from 267 worldwide populations. We perform a systematic analysis of genetic relatedness, detecting 240 intra-population and 92 inter-population pairs of previously unidentified close relatives and proposing standardized subsets of unrelated individuals for use in future studies. We then augment the human data with a data set of 84 chimpanzees at the 246 loci they share in common with the human samples. Multidimensional scaling and neighbor-joining analyses of these data sets offer new insights into the structure of human populations and enable a comparison of genetic variation patterns in chimpanzees with those in humans. Our combined data sets are the largest of their kind reported to date and provide a resource for use in human population-genetic studies.

High mammographic density in women of Ashkenazi Jewish descent

Introduction

Percent mammographic density (PMD) adjusted for age and body mass index is one of the strongest risk factors for breast cancer and is known to be approximately 60% heritable. Here we report a finding of an association between genetic ancestry and adjusted PMD.

Methods

We selected self-identified Caucasian women in the California Pacific Medical Center Research Institute Cohort whose screening mammograms placed them in the top or bottom quintiles of age-adjusted and body mass index-adjusted PMD. Our final dataset included 474 women with the highest adjusted PMD and 469 with the lowest genotyped on the Illumina 1 M platform. Principal component analysis (PCA) and identity-by-descent analyses allowed us to infer the women's genetic ancestry and correlate it with adjusted PMD.

Results

Women of Ashkenazi Jewish ancestry, as defined by the first principal component of PCA and identity-by-descent analyses, represented approximately 15% of the sample. Ashkenazi Jewish ancestry, defined by the first principal component of PCA, was associated with higher adjusted PMD (P = 0.004). Using multivariate regression to adjust for epidemiologic factors associated with PMD, including age at parity and use of postmenopausal hormone therapy, did not attenuate the association.

Conclusions

Women of Ashkenazi Jewish ancestry, based on genetic analysis, are more likely to have high age-adjusted and body mass index-adjusted PMD. Ashkenazi Jews may have a unique set of genetic variants or environmental risk factors that increase mammographic density.

The KCNJ8-S422L variant previously associated with J-wave syndromes is found at an increased frequency in Ashkenazi Jews

J-wave syndromes have been associated with increased risk of ventricular fibrillation and sudden cardiac death. Previous studies have identified the KCNJ8-S422L variant in heterozygous form in individuals with J-wave syndromes. Its absence in over 1500 controls, coupled with in vitro analysis, have led to the conclusion that S422L is pathogenic. We previously performed whole-genome sequencing in a family quartet of Ashkenazi Jewish decent with no history of J-wave syndrome. Re-examination of these data reveals that both parents are heterozygous for the S422L variant, while the 12-year old son carries two copies--thus representing the first reported case of a S422L homozygote. In order to examine whether the S422L mutation might segregate at appreciable frequencies in specific populations, we genotyped the variant in a panel consisting of 722 individuals from 22 European, Middle Eastern non-Jewish, Ashkenazi Jewish, and non-Ashkenazi Jewish populations. We found that the S422L allele was at a significantly higher frequency in Ashkenazi Jews (~4%) compared with other populations in our survey, which have frequencies <0.25%. We also performed ECGs in both male members of the family quartet. The homozygous boy demonstrated no clinically significant ECG abnormalities, while the heterozygous father presented with a subtle J-wave point elevation. Our results suggest that either (a) previous studies implicating S422L as pathogenic for J-wave syndromes failed to appropriately account for European population structure and the variant is likely benign, or (b) Ashkenazi Jews may be at significantly increased risk of J-wave syndromes and ultimately sudden cardiac death.

Characterization of SNPs associated with prostate cancer in men of Ashkenazic descent from the set of GWAS identified SNPs: impact of cancer family history and cumulative SNP risk prediction

BACKGROUND

Genome-wide association studies (GWAS) have identified multiple SNPs associated with prostate cancer (PrCa). Population isolates may have different sets of risk alleles for PrCa constituting unique population and individual risk profiles.

METHODS

To test this hypothesis, associations between 31 GWAS SNPs of PrCa were examined among 979 PrCa cases and 1,251 controls of Ashkenazic descent using logistic regression. We also investigated risks by age at diagnosis, pathological features of PrCa, and family history of cancer. Moreover, we examined associations between cumulative number of risk alleles and PrCa and assessed the utility of risk alleles in PrCa risk prediction by comparing the area under the curve (AUC) for different logistic models.

RESULTS

Of the 31 genotyped SNPs, 8 were associated with PrCa at p ≤ 0.002 (corrected p-value threshold) with odds ratios (ORs) ranging from 1.22 to 1.42 per risk allele. Four SNPs were associated with aggressive PrCa, while three other SNPs showed potential interactions for PrCa by family history of PrCa (rs8102476; 19q13), lung cancer (rs17021918; 4q22), and breast cancer (rs10896449; 11q13). Men in the highest vs. lowest quartile of cumulative number of risk alleles had ORs of 3.70 (95% CI 2.76-4.97); 3.76 (95% CI 2.57-5.50), and 5.20 (95% CI 2.94-9.19) for overall PrCa, aggressive cancer and younger age at diagnosis, respectively. The addition of cumulative risk alleles to the model containing age at diagnosis and family history of PrCa yielded a slightly higher AUC (0.69 vs. 0.64).

CONCLUSION

These data define a set of risk alleles associated with PrCa in men of Ashkenazic descent and indicate possible genetic differences for PrCa between populations of European and Ashkenazic ancestry. Use of genetic markers might provide an opportunity to identify men at highest risk for younger age of onset PrCa; however, their clinical utility in identifying men at highest risk for aggressive cancer remains limited.

Complex genetic origin of Indian populations and its implications

Indian populations are classified into various caste, tribe and religious groups, which altogether makes them very unique compared to rest of the world. The long-term firm socio-religious boundaries and the strict endogamy practices along with the evolutionary forces have further supplemented the existing high-level diversity. As a result, drawing definite conclusions on its overall origin, affinity, health and disease conditions become even more sophisticated than was thought earlier. In spite of these challenges, researchers have undertaken tireless and extensive investigations using various genetic markers to estimate genetic variation and its implication in health and diseases. We have demonstrated that the Indian populations are the descendents of the very first modern humans, who ventured the journey of out-of-Africa about 65,000 years ago. The recent gene flow from east and west Eurasia is also evident. Thus, this review attempts to summarize the unique genetic variation among Indian populations as evident from our extensive study among approximately 20,000 samples across India.

A genome-wide analysis of populations from European Russia reveals a new pole of genetic diversity in northern Europe

Several studies examined the fine-scale structure of human genetic variation in Europe. However, the European sets analyzed represent mainly northern, western, central, and southern Europe. Here, we report an analysis of approximately 166,000 single nucleotide polymorphisms in populations from eastern (northeastern) Europe: four Russian populations from European Russia, and three populations from the northernmost Finno-Ugric ethnicities (Veps and two contrast groups of Komi people). These were compared with several reference European samples, including Finns, Estonians, Latvians, Poles, Czechs, Germans, and Italians. The results obtained demonstrated genetic heterogeneity of populations living in the region studied. Russians from the central part of European Russia (Tver, Murom, and Kursk) exhibited similarities with populations from central-eastern Europe, and were distant from Russian sample from the northern Russia (Mezen district, Archangelsk region). Komi samples, especially Izhemski Komi, were significantly different from all other populations studied. These can be considered as a second pole of genetic diversity in northern Europe (in addition to the pole, occupied by Finns), as they had a distinct ancestry component. Russians from Mezen and the Finnic-speaking Veps were positioned between the two poles, but differed from each other in the proportions of Komi and Finnic ancestries. In general, our data provides a more complete genetic map of Europe accounting for the diversity in its most eastern (northeastern) populations.

Genome-wide diversity in the levant reveals recent structuring by culture

The Levant is a region in the Near East with an impressive record of continuous human existence and major cultural developments since the Paleolithic period. Genetic and archeological studies present solid evidence placing the Middle East and the Arabian Peninsula as the first stepping-stone outside Africa. There is, however, little understanding of demographic changes in the Middle East, particularly the Levant, after the first Out-of-Africa expansion and how the Levantine peoples relate genetically to each other and to their neighbors. In this study we analyze more than 500,000 genome-wide SNPs in 1,341 new samples from the Levant and compare them to samples from 48 populations worldwide. Our results show recent genetic stratifications in the Levant are driven by the religious affiliations of the populations within the region. Cultural changes within the last two millennia appear to have facilitated/maintained admixture between culturally similar populations from the Levant, Arabian Peninsula, and Africa. The same cultural changes seem to have resulted in genetic isolation of other groups by limiting admixture with culturally different neighboring populations. Consequently, Levant populations today fall into two main groups: one sharing more genetic characteristics with modern-day Europeans and Central Asians, and the other with closer genetic affinities to other Middle Easterners and Africans. Finally, we identify a putative Levantine ancestral component that diverged from other Middle Easterners ∼23,700–15,500 years ago during the last glacial period, and diverged from Europeans ∼15,900–9,100 years ago between the last glacial warming and the start of the Neolithic.

Population stratification caused by nonrandom mating between groups of the same species is often due to geographical distances leading to physical separation followed by genetic drift of allele frequencies in each group. In humans, population structures are also often driven by geographical barriers or distances; however, humans might also be structured by abstract factors such as culture, a consequence of their reasoning and self-awareness. Religion in particular, is one of the unusual conceptual factors that can drive human population structures. This study explores the Levant, a region flanked by the Middle East and Europe, where individual and population relationships are still strongly influenced by religion. We show that religious affiliation had a strong impact on the genomes of the Levantines. In particular, conversion of the region's populations to Islam appears to have introduced major rearrangements in populations' relations through admixture with culturally similar but geographically remote populations, leading to genetic similarities between remarkably distant populations like Jordanians, Moroccans, and Yemenis. Conversely, other populations, like Christians and Druze, became genetically isolated in the new cultural environment. We reconstructed the genetic structure of the Levantines and found that a pre-Islamic expansion Levant was more genetically similar to Europeans than to Middle Easterners.

The variance of identity-by-descent sharing in the Wright-Fisher model

Widespread sharing of long, identical-by-descent (IBD) genetic segments is a hallmark of populations that have experienced recent genetic drift. Detection of these IBD segments has recently become feasible, enabling a wide range of applications from phasing and imputation to demographic inference. Here, we study the distribution of IBD sharing in the Wright-Fisher model. Specifically, using coalescent theory, we calculate the variance of the total sharing between random pairs of individuals. We then investigate the cohort-averaged sharing: the average total sharing between one individual and the rest of the cohort. We find that for large cohorts, the cohort-averaged sharing is distributed approximately normally. Surprisingly, the variance of this distribution does not vanish even for large cohorts, implying the existence of "hypersharing" individuals. The presence of such individuals has consequences for the design of sequencing studies, since, if they are selected for whole-genome sequencing, a larger fraction of the cohort can be subsequently imputed. We calculate the expected gain in power of imputation by IBD and subsequently in power to detect an association, when individuals are either randomly selected or specifically chosen to be the hypersharing individuals. Using our framework, we also compute the variance of an estimator of the population size that is based on the mean IBD sharing and the variance in the sharing between inbred siblings. Finally, we study IBD sharing in an admixture pulse model and show that in the Ashkenazi Jewish population the admixture fraction is correlated with the cohort-averaged sharing.

Genetic characterization of northeastern Italian population isolates in the context of broader European genetic diversity

Population genetic studies on European populations have highlighted Italy as one of genetically most diverse regions. This is possibly due to the country's complex demographic history and large variability in terrain throughout the territory. This is the reason why Italy is enriched for population isolates, Sardinia being the best-known example. As the population isolates have a great potential in disease-causing genetic variants identification, we aimed to genetically characterize a region from northeastern Italy, which is known for isolated communities. Total of 1310 samples, collected from six geographically isolated villages, were genotyped at >145 000 single-nucleotide polymorphism positions. Newly genotyped data were analyzed jointly with the available genome-wide data sets of individuals of European descent, including several population isolates. Despite the linguistic differences and geographical isolation the village populations still show the greatest genetic similarity to other Italian samples. The genetic isolation and small effective population size of the village populations is manifested by higher levels of genomic homozygosity and elevated linkage disequilibrium. These estimates become even more striking when the detected substructure is taken into account. The observed level of genetic isolation in Friuli-Venezia Giulia region is more extreme according to several measures of isolation compared with Sardinians, French Basques and northern Finns, thus proving the status of an isolate.

Identification of six novel P450 oxidoreductase missense variants in Ashkenazi and Moroccan Jewish populations

BACKGROUND

The enzyme NADPH-P450 oxidoreductase (POR) is the main electron donor to all microsomal CYPs. The possible contribution of common POR variants to inter- and intra-individual variability in drug metabolism is of great pharmacogenetic interest.

AIM

To search for POR polymorphic alleles and estimate their frequencies in a Jewish population.

MATERIALS & METHODS

We analyzed the POR gene in 301 Ashkenazi and Moroccan Jews.

RESULTS

A total of 30 POR SNPs were identified, nine in the noncoding regions and 21 in the protein-coding regions (ten synonymous, 11 missense). Six of these missense variants are previously undescribed (S102P, V164M, V191M, D344N, E398A and D648N).

CONCLUSION

The data collected in this study on missense POR SNPs, interpreted in light of the crystallographic structure of human POR, indicate that some POR missense variants may be potential biomarkers for future POR pharmacogenetic screening.

Reconstructing the population history of European Romani from genome-wide data

The Romani, the largest European minority group with approximately 11 million people, constitute a mosaic of languages, religions, and lifestyles while sharing a distinct social heritage. Linguistic and genetic studies have located the Romani origins in the Indian subcontinent. However, a genome-wide perspective on Romani origins and population substructure, as well as a detailed reconstruction of their demographic history, has yet to be provided. Our analyses based on genome-wide data from 13 Romani groups collected across Europe suggest that the Romani diaspora constitutes a single initial founder population that originated in north/northwestern India ~1.5 thousand years ago (kya). Our results further indicate that after a rapid migration with moderate gene flow from the Near or Middle East, the European spread of the Romani people was via the Balkans starting ~0.9 kya. The strong population substructure and high levels of homozygosity we found in the European Romani are in line with genetic isolation as well as differential gene flow in time and space with non-Romani Europeans. Overall, our genome-wide study sheds new light on the origins and demographic history of European Romani.

Divorcing the Late Upper Palaeolithic demographic histories of mtDNA haplogroups M1 and U6 in Africa

Background

A Southwest Asian origin and dispersal to North Africa in the Early Upper Palaeolithic era has been inferred in previous studies for mtDNA haplogroups M1 and U6. Both haplogroups have been proposed to show similar geographic patterns and shared demographic histories.

Results

We report here 24 M1 and 33 U6 new complete mtDNA sequences that allow us to refine the existing phylogeny of these haplogroups. The resulting phylogenetic information was used to genotype a further 131 M1 and 91 U6 samples to determine the geographic spread of their sub-clades. No southwest Asian specific clades for M1 or U6 were discovered. U6 and M1 frequencies in North Africa, the Middle East and Europe do not follow similar patterns, and their sub-clade divisions do not appear to be compatible with their shared history reaching back to the Early Upper Palaeolithic. The Bayesian Skyline Plots testify to non-overlapping phases of expansion, and the haplogroups’ phylogenies suggest that there are U6 sub-clades that expanded earlier than those in M1. Some M1 and U6 sub-clades could be linked with certain events. For example, U6a1 and M1b, with their coalescent ages of ~20,000–22,000 years ago and earliest inferred expansion in northwest Africa, could coincide with the flourishing of the Iberomaurusian industry, whilst U6b and M1b1 appeared at the time of the Capsian culture.

Conclusions

Our high-resolution phylogenetic dissection of both haplogroups and coalescent time assessments suggest that the extant main branching pattern of both haplogroups arose and diversified in the mid-later Upper Palaeolithic, with some sub-clades concomitantly with the expansion of the Iberomaurusian industry. Carriers of these maternal lineages have been later absorbed into and diversified further during the spread of Afro-Asiatic languages in North and East Africa.

Empirical distributions of F(ST) from large-scale human polymorphism data

Studies of the apportionment of human genetic variation have long established that most human variation is within population groups and that the additional variation between population groups is small but greatest when comparing different continental populations. These studies often used Wright's F(ST) that apportions the standardized variance in allele frequencies within and between population groups. Because local adaptations increase population differentiation, high-F(ST) may be found at closely linked loci under selection and used to identify genes undergoing directional or heterotic selection. We re-examined these processes using HapMap data. We analyzed 3 million SNPs on 602 samples from eight worldwide populations and a consensus subset of 1 million SNPs found in all populations. We identified four major features of the data: First, a hierarchically F(ST) analysis showed that only a paucity (12%) of the total genetic variation is distributed between continental populations and even a lesser genetic variation (1%) is found between intra-continental populations. Second, the global F(ST) distribution closely follows an exponential distribution. Third, although the overall F(ST) distribution is similarly shaped (inverse J), F(ST) distributions varies markedly by allele frequency when divided into non-overlapping groups by allele frequency range. Because the mean allele frequency is a crude indicator of allele age, these distributions mark the time-dependent change in genetic differentiation. Finally, the change in mean-F(ST) of these groups is linear in allele frequency. These results suggest that investigating the extremes of the F(ST) distribution for each allele frequency group is more efficient for detecting selection. Consequently, we demonstrate that such extreme SNPs are more clustered along the chromosomes than expected from linkage disequilibrium for each allele frequency group. These genomic regions are therefore likely candidates for natural selection.

Genome-wide scan with nearly 700,000 SNPs in two Sardinian sub-populations suggests some regions as candidate targets for positive selection

This paper explores the genetic structure and signatures of natural selection in different sub-populations from the Island of Sardinia, exploiting information from nearly 700,000 autosomal SNPs genotyped with the Affymetrix Genome-Wide Human SNP 6.0 Array. The genetic structure of the Sardinian population and its position within the context of other Mediterranean and European human groups were investigated in depth by comparing our data with publicly available data sets. Principal components and admixture analyses suggest a clustering of the examined samples in two significantly differentiated sub-populations (Ogliastra and Southern Sardinia), as confirmed by AMOVA (F(ST)=0.011; P<0.001). Differentiation of these sub-populations was still evident when they were pooled together with supplementary Sardinian samples from HGDP and compared with several other European, North-African and Near Eastern populations, confirming the uniqueness of the Sardinian genetic background. Moreover, by applying several statistical approaches aimed at assessing differences at the SNP level, the highest differentiated genomic regions between Ogliastra and Southern Sardinia were thus investigated via an extended haplotype homozygosity (EHH)-based test to point out potential selective sweeps. Using this approach, 40 genomic regions were detected, with significant differences between Ogliastra and Southern Sardinia. These regions were subsequently investigated using a long-range haplotype test, which found significant REHH values for SNPs rs11070188 and rs11070192 in the Ogliastra sub-population. In the light of these results and the overlap of the different computed statistics, the region encompassing these loci can be considered a strong candidate to have undergone selective pressure in Ogliastra.

NetView: a high-definition network-visualization approach to detect fine-scale population structures from genome-wide patterns of variation

High-throughput sequencing and single nucleotide polymorphism (SNP) genotyping can be used to infer complex population structures. Fine-scale population structure analysis tracing individual ancestry remains one of the major challenges. Based on network theory and recent advances in SNP chip technology, we investigated an unsupervised network clustering method called Super Paramagnetic Clustering (Spc). When applied to whole-genome marker data it identifies the natural divisions of groups of individuals into population clusters without use of prior ancestry information. Furthermore, we optimised an analysis pipeline called NetView, a high-definition network visualization, starting with computation of genetic distance, followed clustering using Spc and finally visualization of clusters with Cytoscape. We compared NetView against commonly used methodologies including Principal Component Analyses (PCA) and a model-based algorithm, Admixture, on whole-genome-wide SNP data derived from three previously described data sets: simulated (2.5 million SNPs, 5 populations), human (1.4 million SNPs, 11 populations) and cattle (32,653 SNPs, 19 populations). We demonstrate that individuals can be effectively allocated to their correct population whilst simultaneously revealing fine-scale structure within the populations. Analyzing the human HapMap populations, we identified unexpected genetic relatedness among individuals, and population stratification within the Indian, African and Mexican samples. In the cattle data set, we correctly assigned all individuals to their respective breeds and detected fine-scale population sub-structures reflecting different sample origins and phenotypes. The NetView pipeline is computationally extremely efficient and can be easily applied on large-scale genome-wide data sets to assign individuals to particular populations and to reproduce fine-scale population structures without prior knowledge of individual ancestry. NetView can be used on any data from which a genetic relationship/distance between individuals can be calculated.

The population genetics of the Jewish people

Adherents to the Jewish faith have resided in numerous geographic locations over the course of three millennia. Progressively more detailed population genetic analysis carried out independently by multiple research groups over the past two decades has revealed a pattern for the population genetic architecture of contemporary Jews descendant from globally dispersed Diaspora communities. This pattern is consistent with a major, but variable component of shared Near East ancestry, together with variable degrees of admixture and introgression from the corresponding host Diaspora populations. By combining analysis of monoallelic markers with recent genome-wide variation analysis of simple tandem repeats, copy number variations, and single-nucleotide polymorphisms at high density, it has been possible to determine the relative contribution of sex-specific migration and introgression to map founder events and to suggest demographic histories corresponding to western and eastern Diaspora migrations, as well as subsequent microevolutionary events. These patterns have been congruous with the inferences of many, but not of all historians using more traditional tools such as archeology, archival records, linguistics, comparative analysis of religious narrative, liturgy and practices. Importantly, the population genetic architecture of Jews helps to explain the observed patterns of health and disease-relevant mutations and phenotypes which continue to be carefully studied and catalogued, and represent an important resource for human medical genetics research. The current review attempts to provide a succinct update of the more recent developments in a historical and human health context.

Genomic view on the peopling of India

India is known for its vast human diversity, consisting of more than four and a half thousand anthropologically well-defined populations. Each population differs in terms of language, culture, physical features and, most importantly, genetic architecture. The size of populations varies from a few hundred to millions. Based on the social structure, Indians are classified into various caste, tribe and religious groups. These social classifications are very rigid and have remained undisturbed by emerging urbanisation and cultural changes. The variable social customs, strict endogamy marriage practices, long-term isolation and evolutionary forces have added immensely to the diversification of the Indian populations. These factors have also led to these populations acquiring a set of Indian-specific genetic variations responsible for various diseases in India. Interestingly, most of these variations are absent outside the Indian subcontinent. Thus, this review is focused on the peopling of India, the caste system, marriage practice and the resulting health and forensic implications.

North African Jewish and non-Jewish populations form distinctive, orthogonal clusters

North African Jews constitute the second largest Jewish Diaspora group. However, their relatedness to each other; to European, Middle Eastern, and other Jewish Diaspora groups; and to their former North African non-Jewish neighbors has not been well defined. Here, genome-wide analysis of five North African Jewish groups (Moroccan, Algerian, Tunisian, Djerban, and Libyan) and comparison with other Jewish and non-Jewish groups demonstrated distinctive North African Jewish population clusters with proximity to other Jewish populations and variable degrees of Middle Eastern, European, and North African admixture. Two major subgroups were identified by principal component, neighbor joining tree, and identity-by-descent analysis-Moroccan/Algerian and Djerban/Libyan-that varied in their degree of European admixture. These populations showed a high degree of endogamy and were part of a larger Ashkenazi and Sephardic Jewish group. By principal component analysis, these North African groups were orthogonal to contemporary populations from North and South Morocco, Western Sahara, Tunisia, Libya, and Egypt. Thus, this study is compatible with the history of North African Jews-founding during Classical Antiquity with proselytism of local populations, followed by genetic isolation with the rise of Christianity and then Islam, and admixture following the emigration of Sephardic Jews during the Inquisition.

Methods for detecting and correcting for population stratification

Population stratification (PS) is a primary consideration in studies of the genetic determinants of human traits. Failure to control for it may lead to confounding, causing a study to fail for lack of significant results or resources to be wasted following false-positive signals. Here we review historical and current approaches for addressing PS when performing genetic association studies in human populations. We describe methods for detecting the presence of PS including global and local ancestry methods. We also describe approaches for accounting for PS when calculating association statistics, such that measures of association are not confounded. This is particularly important since many traits are being examined in minority populations, populations that may inherently feature PS.

Ethiopian genetic diversity reveals linguistic stratification and complex influences on the Ethiopian gene pool

Humans and their ancestors have traversed the Ethiopian landscape for millions of years, and present-day Ethiopians show great cultural, linguistic, and historical diversity, which makes them essential for understanding African variability and human origins. We genotyped 235 individuals from ten Ethiopian and two neighboring (South Sudanese and Somali) populations on an Illumina Omni 1M chip. Genotypes were compared with published data from several African and non-African populations. Principal-component and STRUCTURE-like analyses confirmed substantial genetic diversity both within and between populations, and revealed a match between genetic data and linguistic affiliation. Using comparisons with African and non-African reference samples in 40-SNP genomic windows, we identified "African" and "non-African" haplotypic components for each Ethiopian individual. The non-African component, which includes the SLC24A5 allele associated with light skin pigmentation in Europeans, may represent gene flow into Africa, which we estimate to have occurred ~3 thousand years ago (kya). The non-African component was found to be more similar to populations inhabiting the Levant rather than the Arabian Peninsula, but the principal route for the expansion out of Africa ~60 kya remains unresolved. Linkage-disequilibrium decay with genomic distance was less rapid in both the whole genome and the African component than in southern African samples, suggesting a less ancient history for Ethiopian populations.

Distinguishing the co-ancestries of haplogroup G Y-chromosomes in the populations of Europe and the Caucasus

Haplogroup G, together with J2 clades, has been associated with the spread of agriculture, especially in the European context. However, interpretations based on simple haplogroup frequency clines do not recognize underlying patterns of genetic diversification. Although progress has been recently made in resolving the haplogroup G phylogeny, a comprehensive survey of the geographic distribution patterns of the significant sub-clades of this haplogroup has not been conducted yet. Here we present the haplogroup frequency distribution and STR variation of 16 informative G sub-clades by evaluating 1472 haplogroup G chromosomes belonging to 98 populations ranging from Europe to Pakistan. Although no basal G-M201* chromosomes were detected in our data set, the homeland of this haplogroup has been estimated to be somewhere nearby eastern Anatolia, Armenia or western Iran, the only areas characterized by the co-presence of deep basal branches as well as the occurrence of high sub-haplogroup diversity. The P303 SNP defines the most frequent and widespread G sub-haplogroup. However, its sub-clades have more localized distribution with the U1-defined branch largely restricted to Near/Middle Eastern and the Caucasus, whereas L497 lineages essentially occur in Europe where they likely originated. In contrast, the only U1 representative in Europe is the G-M527 lineage whose distribution pattern is consistent with regions of Greek colonization. No clinal patterns were detected suggesting that the distributions are rather indicative of isolation by distance and demographic complexities.

DOCK4 and CEACAM21 as novel schizophrenia candidate genes in the Jewish population

It is well accepted that schizophrenia has a strong genetic component. Several genome-wide association studies (GWASs) of schizophrenia have been published in recent years; most of them population based with a case-control design. Nevertheless, identifying the specific genetic variants which contribute to susceptibility to the disorder remains a challenging task. A family-based GWAS strategy may be helpful in the identification of schizophrenia susceptibility genes since it is protected against population stratification, enables better accounting for genotyping errors and is more sensitive for identification of rare variants which have a very low frequency in the general population. In this project we implemented a family-based GWAS of schizophrenia in a sample of 107 Jewish-Israeli families. We found one genome-wide significant association in the intron of the DOCK4 gene (rs2074127, p value=1.134×10⁻⁷) and six additional nominally significant association signals with p<1×10⁻⁵. One of the top single nucleotide polymorphisms (p<1×10⁻⁵) which is located in the predicted intron of the CEACAM21 gene was significantly replicated in independent family-based sample of Arab-Israeli origin (rs4803480: p value=0.002; combined p value=9.61×10⁻⁸), surviving correction for multiple testing. Both DOCK4 and CEACAM21 are biologically reasonable candidate genes for schizophrenia although generalizability of the association of DOCK4 with schizophrenia should be investigated in further studies. In addition, gene-wide significant associations were found within three schizophrenia candidate genes: PGBD1, RELN and PRODH, replicating previously reported associations. By application of a family-based strategy to GWAS, our study revealed new schizophrenia susceptibility loci in the Jewish-Israeli population.

Single-cell exome sequencing and monoclonal evolution of a JAK2-negative myeloproliferative neoplasm

Tumor heterogeneity presents a challenge for inferring clonal evolution and driver gene identification. Here, we describe a method for analyzing the cancer genome at a single-cell nucleotide level. To perform our analyses, we first devised and validated a high-throughput whole-genome single-cell sequencing method using two lymphoblastoid cell line single cells. We then carried out whole-exome single-cell sequencing of 90 cells from a JAK2-negative myeloproliferative neoplasm patient. The sequencing data from 58 cells passed our quality control criteria, and these data indicated that this neoplasm represented a monoclonal evolution. We further identified essential thrombocythemia (ET)-related candidate mutations such as SESN2 and NTRK1, which may be involved in neoplasm progression. This pilot study allowed the initial characterization of the disease-related genetic architecture at the single-cell nucleotide level. Further, we established a single-cell sequencing method that opens the way for detailed analyses of a variety of tumor types, including those with high genetic complex between patients.

The effect of high versus low loading on bone strength in middle life

While bone mass and geometry are largely genetically determined, mechanical loading is considered to be an important additional determinant. This study investigates to what extent very high mechanical loading begun at a young age and sustained afterward can affect tibia bone mass and geometry in middle age. Cohorts from a common ethnic background, with a history of very high and very low tibia bone loading based on an assessment of their activities according their strain levels were compared. The study hypothesis was that the tibia bone density and geometric strength parameters would be greater in the high bone loading cohort. Subjects from a group of elite infantry recruits who sustained a 31% incidence of stress fractures during their basic training in 1983, were reviewed 25 years later. The tibia bone strength of 25 of these soldiers, 11 of whom had sustained stress fractures, was compared to a group of 20 subjects who received exemption from military service in 1982-5 because they were religious scholars and who continued these studies afterwards. Anthropometric measurements were made. The bone density and geometric strength of the tibia was assessed by quantitative computerized tomography (QCT). The average daily dietary intake and metabolic expenditure of subjects were assessed by questionnaires. At the 25 year follow-up soldiers were on an average 3 cm taller than the religious scholars (p=0.02) and had lower abdominal girths (p=0.03). There was no difference in the tibia cortical density between cohorts in spite of the fact that the religious scholars had lower daily calcium intakes (p=0.02). Soldiers had stronger tibias based on geometric engineering criteria. The mean area moments of inertia (p=0.02, p=0.04) and polar moments of inertia (p=0.02) were 16% larger in the soldier cohort. By multivariate regression analysis greater height, weight and daily energy expenditure were related to larger bone geometric strength parameters. According to semipartial eta-square analysis, between 39% to 45% of the variance in the area moments of inertia between the cohorts was attributable to these three parameters. The religious scholars burned less calories daily, principally because they did no sport activity (p=0.001). There was no difference in tibia bone strength parameters between soldiers who did and did not sustain stress fractures in their 1983 basic training. In conclusion, in a middle age population with a common ethnic origin, the high bone loading cohort had stronger tibias than the low bone loading cohort based on larger geometric strength properties and not because of higher cortical density. In spite of being at the extremes of the bone loading spectra, the tibia area moment of inertia of the two cohorts in this study differed by only 16%, with part of this difference attributable to factors other than bone loading. We do not know for sure if the difference in the geometric properties is related to high bone loading or whether people with stronger bones are more likely to engage in high bone loading. Healthy male subjects who sustained stress fractures at a young age do not have weaker tibias at middle age according to QCT measurements.

Association of the ZFPM2 gene with antipsychotic-induced parkinsonism in schizophrenia patients

RATIONALE

Antipsychotic-induced parkinsonism (AIP) is a severe adverse affect of antipsychotic drug treatment. Recently, our group performed a genome-wide association study (GWAS) for AIP severity, and identified several potential AIP risk variants.

OBJECTIVES

The aim of this study was to validate our original AIP-GWAS susceptibility variants and to understand their possible function.

METHODS

We conducted a validation study of 15 single-nucleotide polymorphisms (SNPs) in an independent sample of 178 US schizophrenia patients treated for at least a month with typical or atypical antipsychotics. Then, a sample of 49 Jewish Israeli Parkinson's disease (PD) patients with available neuroimaging ([(123)I]-FP-CIT-SPECT) data was analyzed, to study association of confirmed AIP SNPs with level of dopaminergic deficits in the putamen.

RESULTS

Using logistic regression and controlling for possible confounders, we found nominal association of the intronic SNP, rs12678719, in the Zinc Finger Protein Multitype 2 (ZFPM2) gene with AIP (62 affected/116 unaffected), in the whole sample (p = 0.009; P = 5.97 × 10(-5) in the GWAS), and in the African American sub-sample (N = 111; p = 0.002). The same rs12678719-G AIP susceptibility allele was associated with lower levels of dopaminergic neuron related ligand binding in the contralateral putamen of PD patients (p = 0.026).

CONCLUSIONS

Our preliminary findings support association of the ZFPM2 SNP, rs12678719, with AIP. At the functional level, this variant is associated with deficits in the nigrostriatal pathway in PD patients that may be related to latent subclinical deficits among AIP-prone individuals with schizophrenia. Further validation studies in additional populations are required.