European population genetic substructure: further definition of ancestry informative markers for distinguishing among diverse European ethnic groups

Principal Component Analyses (PCA)-based findings in population genetic studies are highly biased and must be reevaluated

Principal Component Analysis (PCA) is a multivariate analysis that reduces the complexity of datasets while preserving data covariance. The outcome can be visualized on colorful scatterplots, ideally with only a minimal loss of information. PCA applications, implemented in well-cited packages like EIGENSOFT and PLINK, are extensively used as the foremost analyses in population genetics and related fields (e.g., animal and plant or medical genetics). PCA outcomes are used to shape study design, identify, and characterize individuals and populations, and draw historical and ethnobiological conclusions on origins, evolution, dispersion, and relatedness. The replicability crisis in science has prompted us to evaluate whether PCA results are reliable, robust, and replicable. We analyzed twelve common test cases using an intuitive color-based model alongside human population data. We demonstrate that PCA results can be artifacts of the data and can be easily manipulated to generate desired outcomes. PCA adjustment also yielded unfavorable outcomes in association studies. PCA results may not be reliable, robust, or replicable as the field assumes. Our findings raise concerns about the validity of results reported in the population genetics literature and related fields that place a disproportionate reliance upon PCA outcomes and the insights derived from them. We conclude that PCA may have a biasing role in genetic investigations and that 32,000-216,000 genetic studies should be reevaluated. An alternative mixed-admixture population genetic model is discussed.

Genetic Diversity on a Rare Terrestrial Orchid, Habenaria linearifolia in South Korea: Implications for Conservation Offered by Genome-Wide Single Nucleotide Polymorphisms

Monitoring intraspecific diversity offers invaluable insights on conservation practices as the variation is the product of species evolution. Accordingly, the role of population genetic diversity has drawn great attention over the last century responding to the biodiversity loss induced by a series of anthropogenic changes. Orchids are one of the most diverse, yet ironically most rapidly disappearing plant groups due to the specialized habitat preferences. Thus, population-level genetic diversity studies may offer a powerful tool for orchid conservation programs. Using the 3 restriction site-associated DNA (3RAD) approach, 2,734 genome-wide single nucleotide polymorphisms (SNPs) were isolated. With the 2,734 SNPs, we investigated genetic diversity and population structure on 72 individuals of Habenaria linearifolia and Habenaria cruciformis in South Korea. Overall, the genetic diversity was well maintained in South Korean Habenaria, but high F _ST values were estimated suggesting large population diversification with limited gene flow. Bayesian assignment analysis revealed a morphologically cryptic diversity pattern in Jeju Island populations, which might serve as an evolutionarily significant unit.

Genetic differentiation and signatures of local adaptation revealed by RADseq for a highly dispersive mud crab Scylla olivacea (Herbst, 1796) in the Sulu Sea

Connectivity of marine populations is shaped by complex interactions between biological and physical processes across the seascape. The influence of environmental features on the genetic structure of populations has key implications for the dynamics and persistence of populations, and an understanding of spatial scales and patterns of connectivity is crucial for management and conservation. This study employed a seascape genomics approach combining larval dispersal modeling and population genomic analysis using single nucleotide polymorphisms (SNPs) obtained from RADseq to examine environmental factors influencing patterns of genetic structure and connectivity for a highly dispersive mud crab Scylla olivacea (Herbst, 1796) in the Sulu Sea. Dispersal simulations reveal widespread but asymmetric larval dispersal influenced by persistent southward and westward surface circulation features in the Sulu Sea. Despite potential for widespread dispersal across the Sulu Sea, significant genetic differentiation was detected among eight populations based on 1,655 SNPs (FST  = 0.0057, p < .001) and a subset of 1,643 putatively neutral SNP markers (FST  = 0.0042, p < .001). Oceanography influences genetic structure, with redundancy analysis (RDA) indicating significant contribution of asymmetric ocean currents to neutral genetic variation ( R adj 2  = 0.133, p = .035). Genetic structure may also reflect demographic factors, with divergent populations characterized by low effective population sizes (N e < 50). Pronounced latitudinal genetic structure was recovered for loci putatively under selection (FST  = 0.2390, p < .001), significantly correlated with sea surface temperature variabilities during peak spawning months for S. olivacea ( R adj 2  = 0.692-0.763; p < .050), suggesting putative signatures of selection and local adaptation to thermal clines. While oceanography and dispersal ability likely shape patterns of gene flow and genetic structure of S. olivacea across the Sulu Sea, the impacts of genetic drift and natural selection influenced by sea surface temperature also appear as likely drivers of population genetic structure. This study contributes to the growing body of literature documenting population genetic structure and local adaptation for highly dispersive marine species, and provides information useful for spatial management of the fishery resource.

Varying Effects of APOE Alleles on Extreme Longevity in European Ethnicities

APOE is a well-studied gene with multiple effects on aging and longevity. The gene has three alleles: e2, e3, and e4, whose frequencies vary by ethnicity. While the e2 is associated with healthy cognitive aging, the e4 allele is associated with Alzheimer's disease and early mortality and therefore its prevalence among people with extreme longevity (EL) is low. Using the PopCluster algorithm, we identified several ethnically different clusters in which the effect of the e2 and e4 alleles on EL changed substantially. For example, PopCluster discovered a large group of 1,309 subjects enriched of Southern Italian genetic ancestry with weaker protective effect of e2 (odds ratio [OR] = 1.27, p = .14) and weaker damaging effect of e4 (OR = 0.82, p = .31) on the phenotype of EL compared to other European ethnicities. Further analysis of this cluster suggests that the odds for EL in carriers of the e4 allele with Southern Italian genetic ancestry differ depending on whether they live in the United States (OR = 0.29, p = .009) or Italy (OR = 1.21, p = .38). PopCluster also found clusters enriched of subjects with Danish ancestry with varying effect of e2 on EL. The country of residence (Denmark or United States) appears to change the odds for EL in the e2 carriers.

The FTO Gene and Measured Food Intake in 5- to 10-Year-Old Children Without Obesity

OBJECTIVE

Genetic variation in the first intron of FTO (e.g., single-nucleotide polymorphism [SNP] rs9939609) is strongly associated with adiposity. This effect is thought to be mediated (at least in part) via increasing caloric intake, although the precise molecular genetic mechanisms are not fully understood. Prior pediatric studies of FTO have included youth with overweight and obesity; however, they have not informed whether a genotypic effect on ingestive behavior is present prior to obesity onset. Therefore, this study investigated the association between FTO and caloric intake in children aged 5 to 10 years without obesity (adiposity ≤ 95th percentile).

METHODS

A total of 122 children were genotyped for rs9939609 and ate ad libitum from a laboratory lunch buffet following a standardized breakfast. Linear regressions, adjusting for body mass, were used to examine the association between FTO "dose" (number of copies of SNP rs9939609) and intake variables.

RESULTS

There was a significant association between FTO and total intake. Each risk allele predicted an additional 64 calories, accounting for 3% of the variance. There were no associations between FTO and macronutrient preference, energy density, or diet variety. Results were influenced by race.

CONCLUSIONS

Results corroborate and extend prior work by showing a dose-dependent effect on food intake in children without obesity.

Selection of Optimal Ancestry Informative Markers for Classification and Ancestry Proportion Estimation in Pigs

Using small sets of ancestry informative markers (AIMs) constitutes a cost-effective method to accurately estimate the ancestry proportions of individuals. This study aimed to generate a small and effective number of AIMs from ∼60 K single nucleotide polymorphism (SNP) data of porcine and estimate three ancestry proportions [East China pig (ECHP), South China pig (SCHP), and European commercial pig (EUCP)] from Asian breeds and European domestic breeds. A total of 186 samples of 10 pure breeds were divided into three groups: ECHP, SCHP, and EUCP. Using these samples and a one-vs.-rest SVM classifier, we found that using only seven AIMs could completely separate the three groups. Subsequently, we utilized supervised ADMIXTURE to calculate ancestry proportions and found that the 129 AIMs performed well on ancestry estimates when pseudo admixed individuals were used. Furthermore, another 969 samples of 61 populations were applied to evaluate the performance of the 129 AIMs. We also observed that the 129 AIMs were highly correlated with estimates using ∼60 K SNP data for three ancestry components: ECHP (Pearson correlation coefficient (r) = 0.94), SCHP (r = 0.94), and EUCP (r = 0.99). Our results provided an example of using a small number of pig AIMs for classifications and estimating ancestry proportions with high accuracy and in a cost-effective manner.

A single nucleotide polymorphism panel for individual identification and ancestry assignment in Caucasians and four East and Southeast Asian populations using a machine learning classifier

Single nucleotide polymorphism (SNP) profiling is an effective means of individual identification and ancestry inferences in forensic genetics. This study established a SNP panel for the simultaneous individual identification and ancestry assignment of Caucasian and four East and Southeast Asian populations. We analyzed 220 SNPs (125 autosomal, 17 X-chromosomal, 30 Y-chromosomal, and 48 mitochondrial SNPs) of the DNA samples from 563 unrelated individuals of five populations (89 Caucasian, 234 Taiwanese Han, 90 Filipino, 79 Indonesian and 71 Vietnamese) and 18 degraded DNA samples. Informativeness for assignment (In) was used to select ancestry informative SNPs (AISNPs). A machine learning classifier, support vector machine (SVM), was used for ancestry assignment. Of the 220 SNPs, 62 were individual identification SNPs (IISNPs) (51 autosomal and 11 X-chromosomal SNPs) and 191 were AISNPs (100 autosomal, 13 X-chromosomal, 30 Y-chromosomal, and 48 mitochondrial SNPs). The 51 autosomal IISNPs offered cumulative random match probabilities (cRMPs) ranging from 1.56 × 10^-21 to 3.16 × 10^-22 among these five populations. Using AISNPs with the SVM, the overall accuracy rate of ancestry inference achieved in the testing dataset between Caucasian, Taiwanese Han, and Filipino populations was 88.9%, whereas it was 70.0% between Caucasians and each of the four East and Southeast Asian populations. For the 18 degraded DNA samples with incomplete profiling, the accuracy rate of ancestry assignment was 94.4%. We have developed a 220-SNP panel for simultaneous individual identification and ethnic origin differentiation between Caucasian and the four East and Southeast Asian populations. This SNP panel may assist with DNA analysis of forensic casework.

Pancreatic adenocarcinoma, chronic pancreatitis, and MODY-8 diabetes: is bile salt-dependent lipase (or carboxyl ester lipase) at the crossroads of pancreatic pathologies?

Pancreatic adenocarcinomas and diabetes mellitus are responsible for the deaths of around two million people each year worldwide. Patients with chronic pancreatitis do not die directly of this disease, except where the pathology is hereditary. Much current literature supports the involvement of bile salt-dependent lipase (BSDL), also known as carboxyl ester lipase (CEL), in the pathophysiology of these pancreatic diseases. The purpose of this review is to shed light on connections between chronic pancreatitis, diabetes, and pancreatic adenocarcinomas by gaining an insight into BSDL and its variants. This enzyme is normally secreted by the exocrine pancreas, and is diverted within the intestinal lumen to participate in the hydrolysis of dietary lipids. However, BSDL is also expressed by other cells and tissues, where it participates in lipid homeostasis. Variants of BSDL resulting from germline and/or somatic mutations (nucleotide insertion/deletion or nonallelic homologous recombination) are expressed in the pancreas of patients with pancreatic pathologies such as chronic pancreatitis, MODY-8, and pancreatic adenocarcinomas. We discuss the possible link between the expression of BSDL variants and these dramatic pancreatic pathologies, putting forward the suggestion that BSDL and its variants are implicated in the cell lipid metabolism/reprogramming that leads to the dyslipidemia observed in chronic pancreatitis, MODY-8, and pancreatic adenocarcinomas. We also propose potential strategies for translation to therapeutic applications.

Population genomic analysis suggests strong influence of river network on spatial distribution of genetic variation in invasive saltcedar across the southwestern United States

Understanding the complex influences of landscape and anthropogenic elements that shape the population genetic structure of invasive species provides insight into patterns of colonization and spread. The application of landscape genomics techniques to these questions may offer detailed, previously undocumented insights into factors influencing species invasions. We investigated the spatial pattern of genetic variation and the influences of landscape factors on population similarity in an invasive riparian shrub, saltcedar (Tamarix L.) by analysing 1,997 genomewide SNP markers for 259 individuals from 25 populations collected throughout the southwestern United States. Our results revealed a broad-scale spatial genetic differentiation of saltcedar populations between the Colorado and Rio Grande river basins and identified potential barriers to population similarity along both river systems. River pathways most strongly contributed to population similarity. In contrast, low temperature and dams likely served as barriers to population similarity. We hypothesize that large-scale geographic patterns in genetic diversity resulted from a combination of early introductions from distinct populations, the subsequent influence of natural selection, dispersal barriers and founder effects during range expansion.

Evaluating genetic ancestry and self-reported ethnicity in the context of carrier screening

Background

Current professional society guidelines recommend genetic carrier screening be offered on the basis of ethnicity, or when using expanded carrier screening panels, they recommend to compute residual risk based on ethnicity. We investigated the reliability of self-reported ethnicity in 9138 subjects referred to carrier screening. Self-reported ethnicity gathered from test requisition forms and during post-test genetic counseling, and genetic ancestry predicted by a statistical model, were compared for concordance.

Results

We identified several discrepancies between the two sources of self-reported ethnicity and genetic ancestry. Only 30.3% of individuals who indicated Mediterranean ancestry during consultation self-reported this on requisition forms. Additionally, the proportion of individuals who reported Southeast Asian but were estimated to have a different genetic ancestry was found to depend on the source of self-report. Finally, individuals who reported Latin American demonstrated a high degree of ancestral admixture. As a result, carrier rates and residual risks provided for patient decision-making are impacted if using self-reported ethnicity.

Conclusion

Our analysis highlights the unreliability of ethnicity classification based on patient self-reports. We recommend the routine use of pan-ethnic carrier screening panels in reproductive medicine. Furthermore, the use of an ancestry model would allow better estimation of carrier rates and residual risks.

Electronic supplementary material

The online version of this article (10.1186/s12863-017-0570-y) contains supplementary material, which is available to authorized users.

Assortative Mating by Ethnicity in Longevous Families

Recent work shows strong evidence of ancestry-based assortative mating in spouse pairs of the older generation of the Framingham Heart Study. Here, we extend this analysis to two studies of human longevity: the Long Life Family Study (LLFS), and the New England Centenarian Study (NECS). In the LLFS, we identified 890 spouse pairs spanning two generations, while in the NECS we used data from 102 spouse pairs including offspring of centenarians. We used principal components of genome-wide genotype data to demonstrate strong evidence of ancestry-based assortative mating in spouse pairs of the older generation and also confirm the decreasing trend of endogamy in more recent generations. These findings in studies of human longevity suggest that spouses marrying into longevous families may not be powerful controls for genetic association studies, and that there may be important ethnicity-specific, genetic influences and/or gene–environment interactions that influence extreme survival in old generations. In addition, the decreasing trend of genetic similarity of more recent generations might have ramifications for the incidence of homozygous rare variants necessary for survival to the most extreme ages.

Ancestry inference using principal component analysis and spatial analysis: a distance-based analysis to account for population substructure

Background

Accurate inference of genetic ancestry is of fundamental interest to many biomedical, forensic, and anthropological research areas. Genetic ancestry memberships may relate to genetic disease risks. In a genome association study, failing to account for differences in genetic ancestry between cases and controls may also lead to false-positive results. Although a number of strategies for inferring and taking into account the confounding effects of genetic ancestry are available, applying them to large studies (tens thousands samples) is challenging. The goal of this study is to develop an approach for inferring genetic ancestry of samples with unknown ancestry among closely related populations and to provide accurate estimates of ancestry for application to large-scale studies.

Methods

In this study we developed a novel distance-based approach, Ancestry Inference using Principal component analysis and Spatial analysis (AIPS) that incorporates an Inverse Distance Weighted (IDW) interpolation method from spatial analysis to assign individuals to population memberships.

Results

We demonstrate the benefits of AIPS in analyzing population substructure, specifically related to the four most commonly used tools EIGENSTRAT, STRUCTURE, fastSTRUCTURE, and ADMIXTURE using genotype data from various intra-European panels and European-Americans. While the aforementioned commonly used tools performed poorly in inferring ancestry from a large number of subpopulations, AIPS accurately distinguished variations between and within subpopulations.

Conclusions

Our results show that AIPS can be applied to large-scale data sets to discriminate the modest variability among intra-continental populations as well as for characterizing inter-continental variation. The method we developed will protect against spurious associations when mapping the genetic basis of a disease. Our approach is more accurate and computationally efficient method for inferring genetic ancestry in the large-scale genetic studies.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4166-8) contains supplementary material, which is available to authorized users.

A GHEP-ISFG collaborative study on the genetic variation of 38 autosomal indels for human identification in different continental populations

A collaborative effort was carried out by the Spanish and Portuguese Speaking Working Group of the International Society for Forensic Genetics (GHEP-ISFG) to promote knowledge exchange between associate laboratories interested in the implementation of indel-based methodologies and build allele frequency databases of 38 indels for forensic applications. These databases include populations from different countries that are relevant for identification and kinship investigations undertaken by the participating laboratories. Before compiling population data, participants were asked to type the 38 indels in blind samples from annual GHEP-ISFG proficiency tests, using an amplification protocol previously described. Only laboratories that reported correct results contributed with population data to this study. A total of 5839 samples were genotyped from 45 different populations from Africa, America, East Asia, Europe and Middle East. Population differentiation analysis showed significant differences between most populations studied from Africa and America, as well as between two Asian populations from China and East Timor. Low F_ST values were detected among most European populations. Overall diversities and parameters of forensic efficiency were high in populations from all continents.

What Ancestry Can Tell Us About the Genetic Origins of Inter-Ethnic Differences in Asthma Expression

Differences in asthma prevalence have been described across different populations, suggesting that genetic ancestry can play an important role in this disease. In fact, several studies have demonstrated an association between African ancestry with increased asthma susceptibility and severity, higher immunoglobulin E levels, and lower lung function. In contrast, Native American ancestry has been shown to have a protective role for this disease. Genome-wide association studies have allowed the identification of population-specific genetic variants with varying allele frequency among populations. Additionally, the correlation of genetic ancestry at the chromosomal level with asthma and related traits by means of admixture mapping has revealed regions of the genome where ancestry is correlated with the disease. In this review, we discuss the evidence supporting the association of genetic ancestry with asthma susceptibility and asthma-related traits, and highlight the regions of the genome harboring ancestry-specific genetic risk factors.

Forensic genetic analyses in isolated populations with examples of central European Valachs and Roma

Isolated populations present a constant threat to the correctness of forensic genetic casework. In this review article we present several examples of how analyzing samples from isolated populations can bias the results of the forensic statistics and analyses. We select our examples from isolated populations from central and southeastern Europe, namely the Valachs and the European Roma. We also provide the reader with general strategies and principles to improve the laboratory practice (best practice) and reporting of samples from supposedly isolated populations. These include reporting the precise population data used for computing the forensic statistics, using the appropriate θ correction factor for calculating allele frequencies, typing ancestry informative markers in samples of unknown or uncertain ethnicity and establishing ethnic-specific forensic databases.

Structured mating: Patterns and implications

Genetic similarity of spouses can reflect factors influencing mate choice, such as physical/behavioral characteristics, and patterns of social endogamy. Spouse correlations for both genetic ancestry and measured traits may impact genotype distributions (Hardy Weinberg and linkage equilibrium), and therefore genetic association studies. Here we evaluate white spouse-pairs from the Framingham Heart Study (FHS) original and offspring cohorts (N = 124 and 755, respectively) to explore spousal genetic similarity and its consequences. Two principal components (PCs) of the genome-wide association (GWA) data were identified, with the first (PC1) delineating clines of Northern/Western to Southern European ancestry and the second (PC2) delineating clines of Ashkenazi Jewish ancestry. In the original (older) cohort, there was a striking positive correlation between the spouses in PC1 (r = 0.73, P = 3x10^-22) and also for PC2 (r = 0.80, P = 7x10^-29). In the offspring cohort, the spouse correlations were lower but still highly significant for PC1 (r = 0.38, P = 7x10^-28) and for PC2 (r = 0.45, P = 2x10^-39). We observed significant Hardy-Weinberg disequilibrium for single nucleotide polymorphisms (SNPs) loading heavily on PC1 and PC2 across 3 generations, and also significant linkage disequilibrium between unlinked SNPs; both decreased with time, consistent with reduced ancestral endogamy over generations and congruent with theoretical calculations. Ignoring ancestry, estimates of spouse kinship have a mean significantly greater than 0, and more so in the earlier generations. Adjusting kinship estimates for genetic ancestry through the use of PCs led to a mean spouse kinship not different from 0, demonstrating that spouse genetic similarity could be fully attributed to ancestral assortative mating. These findings also have significance for studies of heritability that are based on distantly related individuals (kinship less than 0.05), as we also demonstrate the poor correlation of kinship estimates in that range when ancestry is or is not taken into account.

We analyzed three generations of whites from the Framingham Heart Study (FHS) using genome-wide genotype data to characterize their genetic ancestry. By examination of spouse-pairs, we observed that individuals of Northern/Western European, Southern European and Ashkenazi ancestry preferentially chose spouses of the same ancestry, however, the degree of endogamy decreased in each successive generation, especially between Northern/Western and Southern Europeans. We then showed that the mating pattern results in Hardy-Weinberg disequilibrium (HWD) at ancestrally-informative SNPs, and also results in linkage disequilibrium (LD) between unlinked loci. The HWD and LD decrease as theoretically expected with the decrease in endogamy noted in each generation. In the FHS sample, spouse genetic similarity can be explained by ancestry-related assortative mating.

TMEM187-IRAK1 Polymorphisms Associated with Rheumatoid Arthritis Susceptibility in Tunisian and French Female Populations: Influence of Geographic Origin

Polymorphisms have been identified in the Xq28 locus as risk loci for rheumatoid arthritis (RA). Here, we investigated the association between three polymorphisms in the Xq28 region containing TMEM187 and IRAK1 (rs13397, rs1059703, and rs1059702) in two unstudied populations: Tunisian and French. The rs13397 G and rs1059703 T major alleles were significantly increased in RA patients (n = 408) compared with age-matched controls (n = 471) in both Tunisian and French women. These results were confirmed by a meta-analysis replication study including two independent Greek and Korean cohorts. The rs1059702 C major allele was significantly associated with RA, only with French women. In the French population, the GTC haplotype displayed a protective effect against RA, while the ATC, GCC, and GTT haplotypes conferred significant risk for RA. No association for these haplotypes was found in the Tunisian population. Our results replicated for the first time the association of the three Xq28 polymorphisms with RA risk in Tunisian and French populations and suggested that RA susceptibility is associated with TMEM187-IRAK1 polymorphisms in women. Our data further support the involvement of X chromosome in RA susceptibility and evidence ethnicities differences that might be explained by differences in the frequencies of SE HLA-DRB1 alleles between both populations.

Expression of truncated bile salt-dependent lipase variant in pancreatic pre-neoplastic lesions

Pancreatic adenocarcinoma (PDAC) is a dismal disease. The lack of specific symptoms still leads to a delay in diagnosis followed by death within months for most patients. Exon 11 of the bile salt-dependent lipase (BSDL) gene encoding variable number of tandem repeated (VNTR) sequences has been involved in pancreatic pathologies. We hypothesized that BSDL VNTR sequences may be mutated in PDAC. The amplification of BSDL VNTR from RNA extracted from pancreatic SOJ-6 cells allowed us to identify a BSDL amplicon in which a cytosine residue is inserted in a VNTR sequence. This insertion gives rise to a premature stop codon, resulting in a truncated protein and to a modification of the C-terminal amino-acid sequence; that is PRAAHG instead of PAVIRF. We produced antibodies directed against these sequences and examined pancreatic tissues from patients with PDAC and PanIN. Albeit all tissues were positive to anti-PAVIRF antibodies, 72.2% of patient tissues gave positive reaction with anti-PRAAHG antibodies, particularly in dysplastic areas of the tumor. Neoplastic cells with ductal differentiation were not reactive to anti-PRAAHG antibodies. Some 70% of PanIN tissues were also reactive to anti-PRAAHG antibodies, suggesting that the C insertion occurs early during pancreatic carcinogenesis. Data suggest that anti-PRAAHG antibodies were uniquely reactive with a short isoform of BSDL specifically expressed in pre-neoplastic lesions of the pancreas. The detection of truncated BSDL reactive to antibodies against the PRAAHG C-terminal sequence in pancreatic juice or in pancreatic biopsies may be a new tool in the early diagnosis of PDAC.

Localizing Ashkenazic Jews to Primeval Villages in the Ancient Iranian Lands of Ashkenaz

The Yiddish language is over 1,000 years old and incorporates German, Slavic, and Hebrew elements. The prevalent view claims Yiddish has a German origin, whereas the opposing view posits a Slavic origin with strong Iranian and weak Turkic substrata. One of the major difficulties in deciding between these hypotheses is the unknown geographical origin of Yiddish speaking Ashkenazic Jews (AJs). An analysis of 393 Ashkenazic, Iranian, and mountain Jews and over 600 non-Jewish genomes demonstrated that Greeks, Romans, Iranians, and Turks exhibit the highest genetic similarity with AJs. The Geographic Population Structure analysis localized most AJs along major primeval trade routes in northeastern Turkey adjacent to primeval villages with names that may be derived from "Ashkenaz." Iranian and mountain Jews were localized along trade routes on the Turkey's eastern border. Loss of maternal haplogroups was evident in non-Yiddish speaking AJs. Our results suggest that AJs originated from a Slavo-Iranian confederation, which the Jews call "Ashkenazic" (i.e., "Scythian"), though these Jews probably spoke Persian and/or Ossete. This is compatible with linguistic evidence suggesting that Yiddish is a Slavic language created by Irano-Turko-Slavic Jewish merchants along the Silk Roads as a cryptic trade language, spoken only by its originators to gain an advantage in trade. Later, in the 9th century, Yiddish underwent relexification by adopting a new vocabulary that consists of a minority of German and Hebrew and a majority of newly coined Germanoid and Hebroid elements that replaced most of the original Eastern Slavic and Sorbian vocabularies, while keeping the original grammars intact.

Glutamatergic and GABAergic susceptibility loci for heroin and cocaine addiction in subjects of African and European ancestry

BACKGROUND

Drug addiction, a leading health problem, is a chronic brain disease with a significant genetic component. Animal models and clinical studies established the involvement of glutamate and GABA neurotransmission in drug addiction. This study was designed to assess if 258 variants in 27 genes of these systems contribute to the vulnerability to develop drug addiction.

METHODS

Four independent analyses were conducted in a sample of 1860 subjects divided according to drug of abuse (heroin or cocaine) and ancestry (African and European).

RESULTS

A total of 11 SNPs in eight genes showed nominally significant associations (P<0.01) with heroin and/or cocaine addiction in one or both ancestral groups but the associations did not survive correction for multiple testing. Of these SNPs, the GAD1 upstream SNP rs1978340 is potentially functional as it was shown to affect GABA concentrations in the cingulate cortex. In addition, SNPs GABRB3 rs7165224; DBI rs12613135; GAD1 SNPs rs2058725, rs1978340, rs2241164; and GRIN2A rs1650420 were previously reported in associations with drug addiction or related phenotypes.

CONCLUSIONS

The study supports the involvement of genetic variation in the glutamatergic and GABAergic systems in drug addiction with partial overlap in susceptibility loci between cocaine and heroin addiction.

Dopaminergic pathway polymorphisms and heroin addiction: further support for association of CSNK1E variants

BACKGROUND & AIM

The dopaminergic pathways have been implicated in the etiology of drug addictions. The aim of this study was to determine if variants in dopaminergic genes are associated with heroin addiction.

MATERIALS & METHODS

The study includes 828 former heroin addicts and 232 healthy controls, of predominantly European ancestry. Ninety seven SNPs (13 genes) were analyzed.

RESULTS

Nine nominally significant associations were observed at CSNK1E, ANKK1, DRD2 and DRD3.

CONCLUSION

The results support our previous report of association of CSNK1E SNP rs1534891 with protection from heroin addiction. CSNK1E interacts with circadian rhythms and DARPP-32 and has been implicated in negative regulation of sensitivity to opioids in rodents. It may be a target for drug addiction treatment. Original submitted 8 August 2014; Revision submitted 8 October 2014.

Susceptibility loci for heroin and cocaine addiction in the serotonergic and adrenergic pathways in populations of different ancestry

BACKGROUND

Drug addiction is influenced by genetic factors.

AIM

To determine if genetic variants in the serotonergic and adrenergic pathways are associated with heroin and/or cocaine addiction.

SUBJECTS & METHODS

The study examined 140 polymorphisms in 19 genes in 1855 subjects with predominantly European or African ancestries.

RESULTS

A total of 38 polymorphisms (13 genes) showed nominal associations, including novel associations in S100A10 (p11) and SLC18A2 (VMAT2). The association of HTR3B SNP rs11606194 with heroin addiction in the European ancestry subgroup remained significant after correction for multiple testing (p(corrected) = 0.04).

CONCLUSION

The study strengthens our previous findings of association of polymorphisms in HTR3A, HTR3B and ADRA1A. The study suggests partial overlap in genetic susceptibility between populations of different ancestry and between heroin and cocaine addiction.

Classification, epidemiology and clinical subgrouping of antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis

It is now 25 years since the first European studies on vasculitis--the anti-neutrophil cytoplasmic antibody (ANCA) standardization project. Over that period of time, there have been major developments in the classification of the vasculitides, which has permitted the conduct of high-quality epidemiology studies. Studying the epidemiology of rare diseases such as the ANCA-associated vasculitides (AAV) poses considerable challenges to epidemiologists. The first is the need for a clear definition of a case with good differentiation from similar disorders. The second is case capture. The vasculitides are rare, and therefore, a large population is required to determine the incidence and prevalence, and this poses questions of feasibility. A large population increases the risk of incomplete case detection but permits a reasonable number of cases to be collected in a practicable time frame, whereas a smaller population requires a much longer time frame to collect the necessary cases, which may also not be feasible. Statistical methods of capture-recapture analysis enable estimates to be made of the number of missing cases. The third is case ascertainment. The AAV are virtually always managed in secondary care, and therefore, hospital-based case ascertainment may be appropriate. Fourthly, the rarity of the conditions makes prospective case-control studies investigating risk factors difficult to conduct because the population size required to achieve statistical confidence is in excess of that which is readily available. Thus, much of the data on risk factors are derived from retrospective studies with inherent potential bias.

Immunogenetic influences on acquisition of HIV-1 infection: consensus findings from two African cohorts point to an enhancer element in IL19 (1q32.2)

Numerous reports have suggested that immunogenetic factors may influence HIV-1 acquisition, yet replicated findings that translate between study cohorts remain elusive. Our work aimed to test several hypotheses about genetic variants within the IL10-IL24 gene cluster that encodes interleukin (IL)-10, IL-19, IL-20, and IL-24. In aggregated data from 515 Rwandans and 762 Zambians with up to 12 years of follow-up, 190 single nucleotide polymorphisms (SNPs) passed quality control procedures. When HIV-1-exposed seronegative subjects (n = 486) were compared with newly seroconverted individuals (n = 313) and seroprevalent subjects (n = 478) who were already infected at enrollment, rs12407485 (G>A) in IL19 showed a robust association signal in adjusted logistic regression models (odds ratio = 0.64, P = 1.7 × 10⁻⁴, and q = 0.033). Sensitivity analyses demonstrated that (i) results from both cohorts and subgroups within each cohort were highly consistent; (ii) verification of HIV-1 infection status after enrollment was critical; and (iii) supporting evidence was readily obtained from Cox proportional hazards models. Data from public databases indicate that rs12407485 is part of an enhancer element for three transcription factors. Overall, these findings suggest that molecular features at the IL19 locus may modestly alter the establishment of HIV-1 infection.

Stress-related genes and heroin addiction: a role for a functional FKBP5 haplotype

BACKGROUND

Stress is a critical risk factor affecting both the development of and the relapse to drug addictions. Drug addictions are caused by genetic, environmental and drug-induced factors. The objective of this hypothesis-driven association study was to determine if genetic variants in stress-related genes are associated with heroin addiction.

METHODS

112 selected genetic variants in 26 stress-related genes were genotyped in 852 case subjects and 238 controls of predominantly European ancestry. The case subjects are former heroin addicts with a history of at least one year of daily multiple uses of heroin, treated at a methadone maintenance treatment program (MMTP). The two most promising SNPs were subsequently tested in an African-American sample comprising of 314 cases and 208 control individuals.

RESULTS

Nineteen single nucleotide polymorphisms (SNPs) in 9 genes (AVP, AVPR1A, CRHR1, CRHR2, FKBP5, GAL, GLRA1, NPY1R and NR3C2) showed nominally significant association with heroin addiction. The associations of two FKBP5 SNPs that are part of one haplotype block, rs1360780 (intron 2) and rs3800373 (the 3' untranslated region), remained significant after correction for multiple testing (Pcorrected=0.03; OR=2.35, Pcorrected=0.0018; OR=2.85, respectively). The two SNPs also showed nominally significant association (P<0.05) with heroin addiction in an independent African-American cohort. FKBP5 is a co-chaperone that regulates glucocorticoid sensitivity. These FKBP5 SNPs were previously associated with diverse affective disorders and showed functional differences in gene expression and stress response. This study also supports our and others' previous reports of association of the GAL SNP rs694066 and the AVPR1A SNPs rs11174811, rs1587097 and rs10784339 with heroin and general drug addiction, respectively.

CONCLUSIONS

This study suggests that variations in the FKBP5 gene contribute to the development of opiate addiction by modulating the stress response. These findings may enhance the understanding of the interaction between stress and heroin addiction.

Association of genetic variation in pharmacodynamic factors with methadone dose required for effective treatment of opioid addiction

AIM

The interindividual variability in the dose required for effective methadone maintenance treatment (MMT) for opioid addiction may be influenced in part by genetic variations in genes encoding pharmacodynamic factors of methadone. This study was conducted to identify some of these variants.

MATERIALS & METHODS

This study focused on 11 genes encoding components of the opioidergic (OPRM1, POMC and ARRB2), the dopaminergic (ANKK1 and DRD2) and the glutamatergic pathways (GRIN1 and GRIN2A), as well as the neurotrophin system (NGFB, BDNF, NTRK1 and NTRK2). The study includes 227 Israeli patients undergoing stable MMT.

RESULTS

Out of the 110 variants analyzed, 12 SNPs (in BDNF, NTRK2, OPRM1, DRD2 and ANKK1) were associated with methadone dose (nominal p < 0.05). Of these SNPs, ANKK1 rs7118900 and DRD2 rs2283265 are known to affect gene expression. Logistic regression of five representative SNPs discriminated between individuals requiring a methadone dose of >120 mg/day and <120 mg/day (p = 0.019), and showed moderate sensitivity and specificity (AUC of 0.63 in receiver operating characteristic analysis).

CONCLUSION

This data should stimulate further research on the potential influence and clinical significance of these variants on MMT.

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.

Inference of human continental origin and admixture proportions using a highly discriminative ancestry informative 41-SNP panel

Background

Accurate determination of genetic ancestry is of high interest for many areas such as biomedical research, personal genomics and forensics. It remains an important topic in genetic association studies, as it has been shown that population stratification, if not appropriately considered, can lead to false-positive and -negative results. While large association studies typically extract ancestry information from available genome-wide SNP genotypes, many important clinical data sets on rare phenotypes and historical collections assembled before the GWAS area are in need of a feasible method (i.e., ease of genotyping, small number of markers) to infer the geographic origin and potential admixture of the study subjects. Here we report on the development, application and limitations of a small, multiplexable ancestry informative marker (AIM) panel of SNPs (or AISNP) developed specifically for this purpose.

Results

Based on worldwide populations from the HGDP, a 41-AIM AISNP panel for multiplex application with the ABI SNPlex and a subset with 31 AIMs for the Sequenome iPLEX system were selected and found to be highly informative for inferring ancestry among the seven continental regions Africa, the Middle East, Europe, Central/South Asia, East Asia, the Americas and Oceania. The panel was found to be least informative for Eurasian populations, and additional AIMs for a higher resolution are suggested. A large reference set including over 4,000 subjects collected from 120 global populations was assembled to facilitate accurate ancestry determination. We show practical applications of this AIM panel, discuss its limitations for admixed individuals and suggest ways to incorporate ancestry information into genetic association studies.

Conclusion

We demonstrated the utility of a small AISNP panel specifically developed to discern global ancestry. We believe that it will find wide application because of its feasibility and potential for a wide range of applications.

Genetic ancestry inference using support vector machines, and the active emergence of a unique American population

We use genotype data from the Marshfield Clinical Research Foundation Personalized Medicine Research Project to investigate genetic similarity and divergence between Europeans and the sampled population of European Americans in Central Wisconsin, USA. To infer recent genetic ancestry of the sampled Wisconsinites, we train support vector machines (SVMs) on the positions of Europeans along top principal components (PCs). Our SVM models partition continent-wide European genetic variance into eight regional classes, which is an improvement over the geographically broader categories of recent ancestry reported by personal genomics companies. After correcting for misclassification error associated with the SVMs (<10%, in all cases), we observe a >14% discrepancy between insular ancestries reported by Wisconsinites and those inferred by SVM. Values of FST as well as Mantel tests for correlation between genetic and European geographic distances indicate minimal divergence between Europe and the local Wisconsin population. However, we find that individuals from the Wisconsin sample show greater dispersion along higher-order PCs than individuals from Europe. Hypothesizing that this pattern is characteristic of nascent divergence, we run computer simulations that mimic the recent peopling of Wisconsin. Simulations corroborate the pattern in higher-order PCs, demonstrate its transient nature, and show that admixture accelerates the rate of divergence between the admixed population and its parental sources relative to drift alone. Together, empirical and simulation results suggest that genetic divergence between European source populations and European Americans in Central Wisconsin is subtle but already under way.

CSK regulatory polymorphism is associated with systemic lupus erythematosus and influences B-cell signaling and activation

C-src tyrosine kinase, Csk, physically interacts with the intracellular phosphatase Lyp (PTPN22) and can modify the activation state of downstream Src kinases, such as Lyn, in lymphocytes. We identified an association of Csk with systemic lupus erythematosus (SLE) and refined its location to an intronic polymorphism rs34933034 (OR 1.32, p = 1.04 × 10⁻⁹). The risk allele is associated with increased CSK expression and augments inhibitory phosphorylation of Lyn. In carriers of the risk allele, B cell receptor (BCR)-mediated activation of mature B cells, as well as plasma IgM, are increased. Moreover, the fraction of transitional B cells is doubled in the cord blood of carriers of the risk allele compared to non-risk haplotypes due to an expansion of the late transitional cells, a stage targeted by selection mechanisms. This suggests that the Lyp-Csk complex increases susceptibility to lupus at multiple maturation and activation points of B cells.

An overview of the genetic structure within the Italian population from genome-wide data

In spite of the common belief of Europe as reasonably homogeneous at genetic level, advances in high-throughput genotyping technology have resolved several gradients which define different geographical areas with good precision. When Northern and Southern European groups were considered separately, there were clear genetic distinctions. Intra-country genetic differences were also evident, especially in Finland and, to a lesser extent, within other European populations. Here, we present the first analysis using the 125,799 genome-wide Single Nucleotide Polymorphisms (SNPs) data of 1,014 Italians with wide geographical coverage. We showed by using Principal Component analysis and model-based individual ancestry analysis, that the current population of Sardinia can be clearly differentiated genetically from mainland Italy and Sicily, and that a certain degree of genetic differentiation is detectable within the current Italian peninsula population. Pair-wise F(ST) statistics Northern and Southern Italy amounts approximately to 0.001 between, and around 0.002 between Northern Italy and Utah residents with Northern and Western European ancestry (CEU). The Italian population also revealed a fine genetic substructure underscoring by the genomic inflation (Sardinia vs. Northern Italy = 3.040 and Northern Italy vs. CEU = 1.427), warning against confounding effects of hidden relatedness and population substructure in association studies.

Two complementary perspectives on inter-individual genetic distance

This paper examines population structure through the prism of pairwise genetic distances. Two complementary perspectives, framed as two simple questions, are explored: Q1: What is the probability that a random pair of individuals from the same local population is more genetically dissimilar than a random pair from two distinct populations? Q2: On average, how genetically different are two individuals from the same local population, in comparison with two individuals chosen from any two distinct populations? Models are developed to provide quantitative answers for the two questions, given allele frequencies across any number of markers from two diploid populations. The probability from Q1 is shown to drop to zero with increasing number of genetic markers even for very closely-related populations and rare alleles. The average genetic dissimilarity of two individuals from distinct populations diverges from the average dissimilarity of two individuals from the same population by a percentage dependent on estimates of population differentiation. This perspective also suggests a measure of population distance based on the intuitive notion of pairwise genetic distance, along with a simple method of estimation. Results from recent empirical research on inter-individual genetic distance in human populations are analyzed in the context of the theoretical framework.

Estimating ancestral proportions in a multi-ethnic US sample: implications for studies of admixed populations

This study was designed to determine the ancestral composition of a multi-ethnic sample collected for studies of drug addictions in New York City and Las Vegas, and to examine the reliability of self-identified ethnicity and three-generation family history data. Ancestry biographical scores for seven clusters corresponding to world major geographical regions were obtained using STRUCTURE, based on genotypes of 168 ancestry informative markers (AIMs), for a sample of 1,291 African Americans (AA), European Americans (EA), and Hispanic Americans (HA) along with data from 1,051 HGDP-CEPH ‘diversity panel’ as a reference. Self-identified ethnicity and family history data, obtained in an interview, were accurate in identifying the individual major ancestry in the AA and the EA samples (approximately 99% and 95%, respectively) but were not useful for the HA sample and could not predict the extent of admixture in any group. The mean proportions of the combined clusters corresponding to European and Middle Eastern populations in the AA sample, revealed by AIMs analysis, were 0.13. The HA subjects, predominantly Puerto Ricans, showed a highly variable hybrid contribution pattern of clusters corresponding to Europe (0.27), Middle East (0.27), Africa (0.20), and Central Asia (0.14). The effect of admixture on allele frequencies is demonstrated for two single-nucleotide polymorphisms (118A > G, 17 C > T) of the mu opioid receptor gene (OPRM1). This study reiterates the importance of AIMs in defining ancestry, especially in admixed populations.

The MTAP-CDKN2A locus confers susceptibility to a naturally occurring canine cancer

BACKGROUND

Advantages offered by canine population substructure, combined with clinical presentations similar to human disorders, makes the dog an attractive system for studies of cancer genetics. Cancers that have been difficult to study in human families or populations are of particular interest. Histiocytic sarcoma is a rare and poorly understood neoplasm in humans that occurs in 15% to 25% of Bernese Mountain Dogs (BMD).

METHODS

Genomic DNA was collected from affected and unaffected BMD in North America and Europe. Both independent and combined genome-wide association studies (GWAS) were used to identify cancer-associated loci. Fine mapping and sequencing narrowed the primary locus to a single gene region.

RESULTS

Both populations shared the same primary locus, which features a single haplotype spanning MTAP and part of CDKN2A and is present in 96% of affected BMD. The haplotype is within the region homologous to human chromosome 9p21, which has been implicated in several types of cancer.

CONCLUSIONS

We present the first GWAS for histiocytic sarcoma in any species. The data identify an associated haplotype in the highly cited tumor suppressor locus near CDKN2A. These data show the power of studying distinctive malignancies in highly predisposed dog breeds.

IMPACT

Here, we establish a naturally occurring model of cancer susceptibility due to CDKN2 dysregulation, thus providing insight about this cancer-associated, complex, and poorly understood genomic region.

Molecular and population analysis of natural selection on the human haptoglobin duplication

Haptoglobin binds free haemoglobin that prevents oxidative damage produced by haemolysis. There is a copy number variant (CNV) in the haptoglobin gene (HP) consisting of two alleles, Hp1 (no duplication), and Hp2 (1.7kb duplication involving two exons). The spread of the Hp2 allele is believed to have taken place under selective pressures conferred by malaria resistance. However, molecular evidence is lacking and Hp did not emerge in genomewide SNPs surveys for evidence of selection. In Europe, there is geographical constancy of Hp2 frequency, indicative of absence of clinal pressures and that modern day European alleles represent a “snapshot” of their out-of-Africa migrations. In this work we test for signatures of natural selection acting on the Hp CNV in a sample from the UK population (Avon Longitudinal Study of Parents and Children, ALSPAC). We present here heterozygosity decay, pairwise F_ST values observed between ALSPAC and 301 populations from all five populated continents, extended haplotype homozygosity analyses involving the CNV and 80 SNPs surrounding the CNV ∼500kb in each direction, and linkage disequilibrium and pairwise haplotypic analyses involving 160 SNPs on chromosome 16q22.1. Taken together, our results represent the first molecular analysis of natural selection in the Hp CNV genetic region.

Bias in effect size of systemic lupus erythematosus susceptibility loci across Europe: a case-control study

Introduction

We aimed to investigate whether the effect size of the systemic lupus erythematosus (SLE) risk alleles varies across European subpopulations.

Methods

European SLE patients (n = 1,742) and ethnically matched healthy controls (n = 2,101) were recruited at 17 centres from 10 different countries. Only individuals with self-reported ancestry from the country of origin were included. In addition, participants were genotyped for top ancestry informative markers and for 25 SLE associated SNPs. The results were used to compare effect sizes between the Central Eureopan and Southern European subgroups.

Results

Twenty of the 25 SNPs showed independent association with SLE, These SNPs showed a significant bias to larger effect sizes in the Southern subgroup, with 15/20 showing this trend (P = 0.019) and a larger mean odds ratio of the 20 SNPs (1.46 vs. 1.34, P = 0.02) as well as a larger difference in the number of risk alleles (2.06 vs. 1.63, P = 0.027) between SLE patients and controls than for Central Europeans. This bias was reflected in a very significant difference in the cumulative genetic risk score (4.31 vs. 3.48, P = 1.8 × 10^-32). Effect size bias was accompanied by a lower number of SLE risk alleles in the Southern subjects, both patients and controls, the difference being more marked between the controls (P = 1.1 × 10^-8) than between the Southern and Central European patients (P = 0.016). Seven of these SNPs showed significant allele frequency clines.

Conclusion

Our findings showed a bias to larger effect sizes of SLE loci in the Southern Europeans relative to the Central Europeans together with clines of SLE risk allele frequencies. These results indicate the need to study risk allele clines and the implications of the polygenic model of inheritance in SLE.

Prediction of breed composition in an admixed cattle population

Swiss Fleckvieh was established in 1970 as a composite of Simmental (SI) and Red Holstein Friesian (RHF) cattle. Breed composition is currently reported based on pedigree information. Information on a large number of molecular markers potentially provides more accurate information. For the analysis, we used Illumina BovineSNP50 Genotyping Beadchip data for 90 pure SI, 100 pure RHF and 305 admixed bulls. The scope of the study was to compare the performance of hidden Markov models, as implemented in structure software, with methods conventionally used in genomic selection [BayesB, partial least squares regression (PLSR), least absolute shrinkage and selection operator (LASSO) variable selection)] for predicting breed composition. We checked the performance of algorithms for a set of 40 492 single nucleotide polymorphisms (SNPs), subsets of evenly distributed SNPs and subsets with different allele frequencies in the pure populations, using F(ST) as an indicator. Key results are correlations of admixture levels estimated with the various algorithms with admixture based on pedigree information. For the full set, PLSR, BayesB and structure performed in a very similar manner (correlations of 0.97), whereas the correlation of LASSO and pedigree admixture was lower (0.93). With decreasing number of SNPs, correlations decreased substantially only for 5% or 1% of all SNPs. With SNPs chosen according to F(ST) , results were similar to results obtained with the full set. Only when using 96 and 48 SNPs with the highest F(ST) , correlations dropped to 0.92 and 0.90 respectively. Reducing the number of pure animals in training sets to 50, 20 and 10 each did not cause a drop in the correlation with pedigree admixture.

Identification of CSK as a systemic sclerosis genetic risk factor through Genome Wide Association Study follow-up

Systemic sclerosis (SSc) is complex autoimmune disease affecting the connective tissue; influenced by genetic and environmental components. Recently, we performed the first successful genome-wide association study (GWAS) of SSc. Here, we perform a large replication study to better dissect the genetic component of SSc. We selected 768 polymorphisms from the previous GWAS and genotyped them in seven replication cohorts from Europe. Overall significance was calculated for replicated significant SNPs by meta-analysis of the replication cohorts and replication-GWAS cohorts (3237 cases and 6097 controls). Six SNPs in regions not previously associated with SSc were selected for validation in another five independent cohorts, up to a total of 5270 SSc patients and 8326 controls. We found evidence for replication and overall genome-wide significance for one novel SSc genetic risk locus: CSK [P-value = 5.04 × 10(-12), odds ratio (OR) = 1.20]. Additionally, we found suggestive association in the loci PSD3 (P-value = 3.18 × 10(-7), OR = 1.36) and NFKB1 (P-value = 1.03 × 10(-6), OR = 1.14). Additionally, we strengthened the evidence for previously confirmed associations. This study significantly increases the number of known putative genetic risk factors for SSc, including the genes CSK, PSD3 and NFKB1, and further confirms six previously described ones.

Implications for health and disease in the genetic signature of the Ashkenazi Jewish population

BACKGROUND

Relatively small, reproductively isolated populations with reduced genetic diversity may have advantages for genomewide association mapping in disease genetics. The Ashkenazi Jewish population represents a unique population for study based on its recent (< 1,000 year) history of a limited number of founders, population bottlenecks and tradition of marriage within the community. We genotyped more than 1,300 Ashkenazi Jewish healthy volunteers from the Hebrew University Genetic Resource with the Illumina HumanOmni1-Quad platform. Comparison of the genotyping data with that of neighboring European and Asian populations enabled the Ashkenazi Jewish-specific component of the variance to be characterized with respect to disease-relevant alleles and pathways.

RESULTS

Using clustering, principal components, and pairwise genetic distance as converging approaches, we identified an Ashkenazi Jewish-specific genetic signature that differentiated these subjects from both European and Middle Eastern samples. Most notably, gene ontology analysis of the Ashkenazi Jewish genetic signature revealed an enrichment of genes functioning in transepithelial chloride transport, such as CFTR, and in equilibrioception, potentially shedding light on cystic fibrosis, Usher syndrome and other diseases over-represented in the Ashkenazi Jewish population. Results also impact risk profiles for autoimmune and metabolic disorders in this population. Finally, residual intra-Ashkenazi population structure was minimal, primarily determined by class 1 MHC alleles, and not related to host country of origin.

CONCLUSIONS

The Ashkenazi Jewish population is of potential utility in disease-mapping studies due to its relative homogeneity and distinct genomic signature. Results suggest that Ashkenazi-associated disease genes may be components of population-specific genomic differences in key functional pathways.

Transancestral mapping of the MHC region in systemic lupus erythematosus identifies new independent and interacting loci at MSH5, HLA-DPB1 and HLA-G

Objectives

Systemic lupus erythematosus (SLE) is a chronic multisystem genetically complex autoimmune disease characterised by the production of autoantibodies to nuclear and cellular antigens, tissue inflammation and organ damage. Genome-wide association studies have shown that variants within the major histocompatibility complex (MHC) region on chromosome 6 confer the greatest genetic risk for SLE in European and Chinese populations. However, the causal variants remain elusive due to tight linkage disequilibrium across disease-associated MHC haplotypes, the highly polymorphic nature of many MHC genes and the heterogeneity of the SLE phenotype.

Methods

A high-density case-control single nucleotide polymorphism (SNP) study of the MHC region was undertaken in SLE cohorts of Spanish and Filipino ancestry using a custom Illumina chip in order to fine-map association signals in these haplotypically diverse populations. In addition, comparative analyses were performed between these two datasets and a northern European UK SLE cohort. A total of 1433 cases and 1458 matched controls were examined.

Results

Using this transancestral SNP mapping approach, novel independent loci were identified within the MHC region in UK, Spanish and Filipino patients with SLE with some evidence of interaction. These loci include HLA-DPB1, HLA-G and MSH5 which are independent of each other and HLA-DRB1 alleles. Furthermore, the established SLE-associated HLA-DRB1*15 signal was refined to an interval encompassing HLA-DRB1 and HLA-DQA1. Increased frequencies of MHC region risk alleles and haplotypes were found in the Filipino population compared with Europeans, suggesting that the greater disease burden in non-European SLE may be due in part to this phenomenon.

Conclusion

These data highlight the usefulness of mapping disease susceptibility loci using a transancestral approach, particularly in a region as complex as the MHC, and offer a springboard for further fine-mapping, resequencing and transcriptomic analysis.