Search

Pleiotropy is a widespread phenomenon that may increase insight into the etiology of biological and disease traits. Since genome-wide association studies frequently provide information on a single trait only, only univariate pleiotropy detection methods are applicable, with yet unknown comparative performance. Here, we compared five such methods with respect to their ability to detect pleiotropy, including meta-analysis, ASSET, conditional false discovery rate (cFDR), cross-phenotype Bayes (CPBayes), and pleiotropic analysis under the composite null hypothesis (PLACO), by performing extended computer simulations that varied the underlying etiological model for pleiotropy for a pair of traits, including the number of causal variants, degree of traits' overlap, effect sizes as well as trait prevalence, and varying sample sizes. Our results indicate that ASSET provides the best trade-off between power and protection against false positives. We then applied ASSET to a previously published International League Against Epilepsy (ILAE) consortium data set on complex epilepsies, comprising genetic generalized epilepsy and focal epilepsy cases and corresponding controls. We identified a novel candidate locus at 17q21.32 and confirmed locus 2q24.3, previously identified to act pleiotropically on both epilepsy subtypes by a mega-analysis. Functional annotation, tissue-specific expression, and regulatory function analysis as well as Bayesian colocalization analysis corroborated this result, rendering 17q21.32 a worthwhile candidate for follow-up studies on pleiotropy in epilepsies.

Region-based genome-wide scans are usually performed by use of a priori chosen analysis regions. Such an approach will likely miss the region comprising the strongest signal and, thus, may result in increased type II error rates and decreased power. Here, we propose a genomic exhaustive scan approach that analyzes all possible subsequences and does not rely on a prior definition of the analysis regions. As a prime instance, we present a computationally ultraefficient implementation using the rare-variant collapsing test for phenotypic association, the genomic exhaustive collapsing scan (GECS). Our implementation allows for the identification of regions comprising the strongest signals in large, genome-wide rare-variant association studies while controlling the family-wise error rate via permutation. Application of GECS to two genomic data sets revealed several novel significantly associated regions for age-related macular degeneration and for schizophrenia. Our approach also offers a high potential to improve genome-wide scans for selection, methylation, and other analyses.

To analyze family-based whole-genome sequence (WGS) data for complex traits, we developed a rare variant (RV) non-parametric linkage (NPL) analysis method, which has advantages over association methods. The RV-NPL differs from the NPL in that RVs are analyzed, and allele sharing among affected relative-pairs is estimated only for minor alleles. Analyzing families can increase power because causal variants with familial aggregation usually have larger effect sizes than those underlying sporadic diseases. Differing from association analysis, for NPL only affected individuals are analyzed, which can increase power, since unaffected family members can be susceptibility variant carriers. RV-NPL is robust to population substructure and admixture, inclusion of nonpathogenic variants, as well as allelic and locus heterogeneity and can readily be applied outside of coding regions. In contrast to analyzing common variants using NPL, where loci localize to large genomic regions (e.g., >50 Mb), mapped regions are well defined for RV-NPL. Using simulation studies, we demonstrate that RV-NPL is substantially more powerful than applying traditional NPL methods to analyze RVs. The RV-NPL was applied to analyze 107 late-onset Alzheimer disease (LOAD) pedigrees of Caribbean Hispanic and European ancestry with WGS data, and statistically significant linkage (LOD ≥ 3.8) was found with RVs in PSMF1 and PTPN21 which have been shown to be involved in LOAD etiology. Additionally, nominally significant linkage was observed with RVs in ABCA7, ACE, EPHA1, and SORL1, genes that were previously reported to be associated with LOAD. RV-NPL is an ideal method to elucidate the genetic etiology of complex familial diseases.

Sarcoidosis and chronic beryllium disease (CBD) are phenocopies, however the latter one has a clear trigger factor that is beryllium exposure. This study analyses single nucleotide polymorphisms (SNPs) in a large cohort for beryllium-exposed persons. SNPs were chosen for their relevance in sarcoidosis. Even though one of largest cohorts of beryllium-exposed persons was analysed, no statistically relevant association between any SNP and CBD could be verified. Notably, some SNPs exhibit inverse OR for beryllium sensitization and CBD with nominally statistical significance, which allows hypothesizing about pathophysiological role of genes for the disease triggering and development.

Hearing loss (HL) is one of the most common sensory defects affecting more than 466 million individuals worldwide. It is clinically and genetically heterogeneous with over 120 genes causing non-syndromic HL identified to date. Here, we performed exome sequencing (ES) on a cohort of Iranian families with no disease-causing variants in known deafness-associated genes after screening with a targeted gene panel. We identified likely causal variants in 20 out of 71 families screened. Fifteen families segregated variants in known deafness-associated genes. Eight families segregated variants in novel candidate genes for HL: DBH, TOP3A, COX18, USP31, TCF19, SCP2, TENM1, and CARMIL1. In the three of these families, intrafamilial locus heterogeneity was observed with variants in both known and novel candidate genes. In aggregate, we were able to identify the underlying genetic cause of HL in nearly 30% of our study cohort using ES. This study corroborates the observation that high-throughput DNA sequencing in populations with high rates of consanguineous marriages represents a more appropriate strategy to elucidate the genetic etiology of heterogeneous conditions such as HL.

The prediction of appearance traits by use of solely genetic information has become an established approach and a number of statistical prediction models have already been developed for this purpose. However, given limited knowledge on appearance genetics, currently available models are incomplete and do not include all causal genetic variants as predictors. Therefore such prediction models may benefit from the inclusion of additional information that acts as a proxy for this unknown genetic background. Use of priors, possibly informed by trait category prevalence values in biogeographic ancestry groups, in a Bayesian framework may thus improve the prediction accuracy of previously predicted externally visible characteristics, but has not been investigated as of yet. In this study, we assessed the impact of using trait prevalence-informed priors on the prediction performance in Bayesian models for eye, hair and skin color as well as hair structure and freckles in comparison to the respective prior-free models. Those prior-free models were either similarly defined either very close to the already established ones by using a reduced predictive marker set. However, these differences in the number of the predictive markers should not affect significantly our main outcomes. We observed that such priors often had a strong effect on the prediction performance, but to varying degrees between different traits and also different trait categories, with some categories barely showing an effect. While we found potential for improving the prediction accuracy of many of the appearance trait categories tested by using priors, our analyses also showed that misspecification of those prior values often severely diminished the accuracy compared to the respective prior-free approach. This emphasizes the importance of accurate specification of prevalence-informed priors in Bayesian prediction modeling of appearance traits. However, the existing literature knowledge on spatial prevalence is sparse for most appearance traits, including those investigated here. Due to the limitations in appearance trait prevalence knowledge, our results render the use of trait prevalence-informed priors in DNA-based appearance trait prediction currently infeasible.

EBioMedicine 72, 103588 (2021). doi:10.1016/j.ebiom.2021.103588, Published by Elsevier, Amsterdam [u.a.]

Ecuador is a multiethnic and pluricultural country with a complex history defined by migration and admixture processes. The present study aims to increase our knowledge on the Ecuadorian Native Amerindian groups and the unique South American Y-chromosome haplogroup C3-MPB373 through the analysis of up to 23 Y-chromosome STRs (Y-STRs) and several Y-SNPs in a sample of 527 Ecuadorians from 7 distinct populations and geographic areas, including Kichwa and non-Kichwa Native Amerindians, Mestizos and Afro-Ecuadorians. Our results reveal the presence of C3-MPB373 both in the Amazonian lowland Kichwa with frequencies up to 28 % and, for the first time, in notable proportions in Kichwa populations from the Ecuadorian highlands. The substantially higher frequencies of C3-MPB373 in the Amazonian lowlands found in Kichwa and Waorani individuals suggest a founder effect in that area. Notably, estimates for the time to the most recent common ancestor (TMRCA) in the range of 7.2-9.0 kya point to an ancient origin of the haplogroup and suggest an early Holocene expansion of C3-MPB373 into South America. Finally, the pairwise genetic distances (RST) separate the Kichwa Salasaka from all the other Native Amerindian and Ecuadorian groups, indicating a so far hidden diversity among the Kichwa-speaking populations and suggesting a more southern origin of this population. In sum, our study provides a more in-depth knowledge of the male genetic structure of the multiethnic Ecuadorian population, as well as a valuable reference dataset for forensic use.

Cigarette smoking is the leading cause of preventable morbidity and mortality. Genetic variation contributes to initiation, regular smoking, nicotine dependence, and cessation. We present a Fagerström Test for Nicotine Dependence (FTND)-based genome-wide association study in 58,000 European or African ancestry smokers. We observe five genome-wide significant loci, including previously unreported loci MAGI2/GNAI1 (rs2714700) and TENM2 (rs1862416), and extend loci reported for other smoking traits to nicotine dependence. Using the heaviness of smoking index from UK Biobank (N = 33,791), rs2714700 is consistently associated; rs1862416 is not associated, likely reflecting nicotine dependence features not captured by the heaviness of smoking index. Both variants influence nearby gene expression (rs2714700/MAGI2-AS3 in hippocampus; rs1862416/TENM2 in lung), and expression of genes spanning nicotine dependence-associated variants is enriched in cerebellum. Nicotine dependence (SNP-based heritability = 8.6%) is genetically correlated with 18 other smoking traits (rg = 0.40–1.09) and co-morbidities. Our results highlight nicotine dependence-specific loci, emphasizing the FTND as a composite phenotype that expands genetic knowledge of smoking., There is strong genetic evidence for cigarette smoking behaviors, yet little is known on nicotine dependence (ND). Here, the authors perform a genome-wide association study on ND in 58,000 smokers, identifying five genome-wide significant loci.

Complex diseases are frequently modeled as following an additive model that excludes both intra- and inter-locus interaction, while at the same time reports on non-additive biological structures are ample, prominently featuring numerous metabolic and signaling pathways. Using extensive forward population simulations, we explored the impact of three basic pathway motifs on the relationship between epidemiological parameters, including disease prevalence, relative risk, sibling recurrence risk as well as causal variant number and allele frequency. We found that some but not all pathway motifs can shift the relationships between these parameters in comparison to the classical additive liability threshold model. The strongest deviations were observed with linear, cascade-like motifs that form an integral part of many reported pathways. We also modeled maturity-onset diabetes of the young (MODY) as a combination of different basic pathway motifs and observed a good concordance in epidemiological parameter values between our simulated data under this model and those reported in the literature. Given the widespread nature of pathways, including those in the etiology of human diseases, our results re-emphasize the need for non-additive interaction modeling of genetic variants to become an additional standard approach in analyzing human genetic data.

Hundreds of thousands of individuals have been genotyped in the past decades using genotyping arrays, representing both a valuable data resource for future biomedical research and a substantial investment in human genetic research. However, novel chip designs and their altered sets of single-nucleotide polymorphisms (SNPs) pose the question of how well established data resources, such as large samples of healthy controls genotyped on legacy arrays, can be combined with newer samples genotyped on those novel arrays using genotype imputation. We exemplarily investigated this question based on genotype data of 30 European and 30 African unrelated samples from the 1000 Genomes project and on markers present on two legacy SNP arrays, namely Affymetrix's Human SNP 6.0 and Illumina's 550k array, and three newer arrays, namely two Axiom arrays from Affymetrix and an OmniExpress array from Illumina. We cross-compared the imputation accuracy as well as efficacy and assessed genotype concordance among these arrays. Although the accuracy of genotype prediction was uniformly high across all arrays, the imputation efficacy, that is, the proportion of successfully imputed markers, differed considerably between array combinations in both sample sets, with legacy arrays showing a trend towards lower efficacy values compared with newer arrays when serving as imputation basis. We conclude that, given the substantial losses of markers covered by the legacy arrays, the re-genotyping of existing samples sets, in particular those of healthy population controls, would be a worthwhile endeavor to secure their continued use in the future.

The study of the genetic structure of different countries within Europe has provided significant insights into their demographic history and population structure. Although France occupies a particular location at the western part of Europe and at the crossroads of migration routes, few population genetic studies have been conducted so far with genome-wide data. In this study, we analyzed SNP-chip genetic data from 2184 individuals born in France who were enrolled in two independent population cohorts. Using FineSTRUCTURE, six different genetic clusters of individuals were found that were very consistent between the two cohorts. These clusters correspond closely to geographic, historical, and linguistic divisions of France, and contain different proportions of ancestry from Stone and Bronze Age populations. By modeling the relationship between genetics and geography using EEMS, we were able to detect gene flow barriers that are similar across the two cohorts and correspond to major rivers and mountain ranges. Estimations of effective population sizes also revealed very similar patterns in both cohorts with a rapid increase of effective population sizes over the last 150 generations similar to other European countries. A marked bottleneck is also consistently seen in the two datasets starting in the 14

Whole-genome and whole-exome sequencing of individual patients allow the study of rare and potentially causative genetic variation. In this study, we sequenced DNA of a trio comprising a boy with very-early-onset inflammatory bowel disease (veoIBD) and his unaffected parents. We identified a rare, X-linked missense variant in the NAPDH oxidase NOX1 gene (c.C721T, p.R241C) in heterozygous state in the mother and in hemizygous state in the patient. We discovered that, in addition, the patient was homozygous for a common missense variant in the CYBA gene (c.T214C, p.Y72H). CYBA encodes the p22phox protein, a cofactor for NOX1. Functional assays revealed reduced cellular ROS generation and antibacterial capacity of NOX1 and p22phox variants in intestinal epithelial cells. Moreover, the identified NADPH oxidase complex variants affected NOD2-mediated immune responses, and p22phox was identified as a novel NOD2 interactor. In conclusion, we detected missense variants in a veoIBD patient that disrupt the host response to bacterial challenges and reduce protective innate immune signaling via NOD2. We assume that the patient's individual genetic makeup favored disturbed intestinal mucosal barrier function.

DNA-based prediction of externally visible characteristics has become an established approach in forensic genetics, with the aim of tracing individuals who are potentially unknown to the investigating authorities but without using this prediction as evidence in court. While a number of prediction models have been proposed, use of prior probabilities in those models has largely been absent. Here, we aim at compiling information on the spatial distribution of eye and hair coloration in order to use this as prior knowledge to improve prediction accuracy. To this end, we conducted a detailed literature review and created maps showing the eye and hair pigmentation prevalence both by countries with available information and by interpolation in order to obtain prior estimates for populations without available data. Furthermore, we assessed the association between these two traits in a very large data set. A strong limitation was the quite low amount of available data, especially outside Europe. We hope that our results will facilitate the improvement of already existing and of novel prediction methods for pigmentation traits and induce further studies on the spatial distribution of these traits.

Predicting appearance phenotypes from genotypes is relevant for various areas of human genetic research and applications such as genetic epidemiology, human history, anthropology, and particularly in forensics. Many appearance phenotypes, and thus their underlying genotypes, are highly correlated, with pigmentation traits serving as primary examples. However, all available genetic prediction models, including those for pigmentation traits currently used in forensic DNA phenotyping, ignore phenotype correlations. Here, we investigated the impact of appearance phenotype correlations on genetic appearance prediction in the exemplary case of three pigmentation traits. We used data for categorical eye, hair and skin colour as well as 41 DNA markers utilized in the recently established HIrisPlex-S system from 762 individuals with complete phenotype and genotype information. Based on these data, we performed genetic prediction modelling of eye, hair and skin colour via three different strategies, namely the established approach of predicting phenotypes solely based on genotypes while not considering phenotype correlations, and two novel approaches that considered phenotype correlations, either incorporating truly observed correlated phenotypes or DNA-predicted correlated phenotypes in addition to the DNA predictors. We found that using truly observed correlated pigmentation phenotypes as additional predictors increased the DNA-based prediction accuracies for almost all eye, hair and skin colour categories, with the largest increase for intermediate eye colour, brown hair colour, dark to black skin colour, and particularly for dark skin colour. Outcomes of dedicated computer simulations suggest that this prediction accuracy increase is due to the additional genetic information that is implicitly provided by the truly observed correlated pigmentation phenotypes used, yet not covered by the DNA predictors applied. In contrast, considering DNA-predicted correlated pigmentation phenotypes as additional predictors did not improve the performance of the genetic prediction of eye, hair and skin colour, which was in line with the results from our computer simulations. Hence, in practical applications of DNA-based appearance prediction where no phenotype knowledge is available, such as in forensic DNA phenotyping, it is not advised to use DNA-predicted correlated phenotypes as predictors in addition to the DNA predictors. In the very least, this is not recommended for the pigmentation traits and the established pigmentation DNA predictors tested here.

Leprosy, a chronic infectious disease caused by Mycobacterium leprae (M. leprae), was very common in Europe till the 16th century. Here, we perform an ancient DNA study on medieval skeletons from Denmark that show lesions specific for lepromatous leprosy (LL). First, we test the remains for M. leprae DNA to confirm the infection status of the individuals and to assess the bacterial diversity. We assemble 10 complete M. leprae genomes that all differ from each other. Second, we evaluate whether the human leukocyte antigen allele DRB1*15:01, a strong LL susceptibility factor in modern populations, also predisposed medieval Europeans to the disease. The comparison of genotype data from 69 M. leprae DNA-positive LL cases with those from contemporary and medieval controls reveals a statistically significant association in both instances. In addition, we observe that DRB1*15:01 co-occurs with DQB1*06:02 on a haplotype that is a strong risk factor for inflammatory diseases today., Leprosy, caused by infection with Mycobacterium leprae, was common in Europe in the Middle Ages. Here, Krause-Kyora et al. analyze ancient DNA from a medieval Danish leprosarium to assemble 10 complete bacterial genomes and perform association analysis of the DRB1*15:01 allele with risk of leprosy infection.

ObjectiveHomozygous alpha1-antitrypsin (AAT) deficiency increases the risk for developing cirrhosis, whereas the relevance of heterozygous carriage remains unclear. Hence, we evaluated the impact of the two most relevant AAT variants (‘Pi*Z’ and ‘Pi*S’), present in up to 10% of Caucasians, on subjects with non-alcoholic fatty liver disease (NAFLD) or alcohol misuse.DesignWe analysed multicentric case–control cohorts consisting of 1184 people with biopsy-proven NAFLD and of 2462 people with chronic alcohol misuse, both cohorts comprising cases with cirrhosis and controls without cirrhosis. Genotyping for the Pi*Z and Pi*S variants was performed.ResultsThe Pi*Z variant presented in 13.8% of patients with cirrhotic NAFLD but only in 2.4% of counterparts without liver fibrosis (pConclusionThe Pi*Z variant is the hitherto strongest single nucleotide polymorphism-based risk factor for cirrhosis in NAFLD and alcohol misuse, whereas the Pi*S variant confers only a weak risk in alcohol misusers. As 2%–4% of Caucasians are Pi*Z carriers, this finding should be considered in genetic counselling of affected individuals.

FOXO3 is consistently annotated as a human longevity gene. However, functional variants and underlying mechanisms for the association remain unknown. Here, we perform resequencing of the FOXO3 locus and single-nucleotide variant (SNV) genotyping in three European populations. We find two FOXO3 SNVs, rs12206094 and rs4946935, to be most significantly associated with longevity and further characterize them functionally. We experimentally validate the in silico predicted allele-dependent binding of transcription factors (CTCF, SRF) to the SNVs. Specifically, in luciferase reporter assays, the longevity alleles of both variants show considerable enhancer activities that are reversed by IGF-1 treatment. An eQTL database search reveals that the alleles are also associated with higher FOXO3 mRNA expression in various human tissues, which is in line with observations in long-lived model organisms. In summary, we present experimental evidence for a functional link between common intronic variants in FOXO3 and human longevity.

Acrocallosal syndrome (ACLS) is an autosomal recessive neurodevelopmental disorder caused by KIF7 defects and belongs to the heterogeneous group of ciliopathies related to Joubert syndrome (JBTS). While ACLS is characterized by macrocephaly, prominent forehead, depressed nasal bridge, and hypertelorism, facial dysmorphism has not been emphasized in JBTS cohorts with molecular diagnosis. To evaluate the specificity and etiology of ACLS craniofacial features, we performed whole exome or targeted Sanger sequencing in patients with the aforementioned overlapping craniofacial appearance but variable additional ciliopathy features followed by functional studies. We found (likely) pathogenic variants of KIF7 in 5 out of 9 families, including the original ACLS patients, and delineated 1000 to 4000-year-old Swiss founder alleles. Three of the remaining families had (likely) pathogenic variants in the JBTS gene C5orf42, and one patient had a novel de novo frameshift variant in SHH known to cause autosomal dominant holoprosencephaly. In accordance with the patients' craniofacial anomalies, we showed facial midline widening after silencing of C5orf42 in chicken embryos. We further supported the link between KIF7, SHH, and C5orf42 by demonstrating abnormal primary cilia and diminished response to a SHH agonist in fibroblasts of C5orf42-mutated patients, as well as axonal pathfinding errors in C5orf42-silenced chicken embryos similar to those observed after perturbation of Shh signaling. Our findings, therefore, suggest that beside the neurodevelopmental features, macrocephaly and facial widening are likely more general signs of disturbed SHH signaling. Nevertheless, long-term follow-up revealed that C5orf42-mutated patients showed catch-up development and fainting of facial features contrary to KIF7-mutated patients.

_Dear Editor_ We highly welcome the development of innovative approaches that help differentiate between monozygotic (MZ) twins based on DNA evidence in forensic casework, given that standard forensic STR profiling does not allow their individual identification. However, while Marqueta-Gracia et al. recently claimed in the Journal that almost half of the tested MZ twin pairs were discriminated based on the same three age-related CpG fragments, we present arguments and evidence demonstrating that this claim is not supported by their data. First, given the same chronological age of MZ twins per pair, applying age-related DNA methylation markers for discriminating MZ individuals within pairs appears counterintuitive. [...] Second, any data-driven claim on individual discrimination depends on the magnitude of technical error of the method used for data generation, which is especially crucial if the observed individual differences are small. [...] Third, Marqueta-Gracia et al. applied an approach to test for statistically significant inter-individual differences that appears inappropriate for the given purpose for several reasons. [...] Fourth, it appears peculiar that, although the authors have analyzed a total of six age-related CpG fragments, they based their final conclusions regarding MZ twin discrimination on only three of them (ITGA2B, ASPA, and ZIC5). [...] Fifth, cell type composition is a major challenge in epigenetic analysis. Differences in DNA methylation patterns between different cell types are well-known and typically require correction in epigenomewide association studies, an approach that should also be followed in targeted DNA methylation analysis. [...] In summary, convincing data evidence that i) age-related CpGs allow for the reliable discrimination of MZ twin individuals within pairs, ii) a small set of the same CpGs allows for reliable MZ twin discrimination across different pairs, and iii) a small set of the same CpGs allows for reliable MZ twin discrimination from forensic trace-type materials, all remain to be presented.

Objective: Increasing availability of surgically resected brain tissue from patients with focal epilepsy and Focal Cortical Dysplasia (FCD) or low-grade glio-neuronal tumors has fostered large-scale genetic examination. However, assessment of pathogenicity of germline and somatic variants remains difficult. Here, we present a state of the art evaluation of reported genes and variants associated with epileptic brain lesions. Methods: We critically re-evaluated the pathogenicity for all neuropathology-associated variants reported to date in PubMed and ClinVar databases including 101 neuropathology-associated missense variants encompassing 11 disease-related genes. We assessed gene variant tolerance and classified all identified missense variants according to guidelines from the American College of Medical Genetics and Genomics (ACMG). We further extended the bioinformatic variant prediction by introducing a novel gene-specific deleteriousness ranking for prediction scores. Results: Application of ACMG guidelines and in silico gene variant tolerance analysis classified only seven out of 11 genes to be likely disease-associated according to the reported a disease mechanism, while 61 (60.4%) of 101 variants of those genes were classified as of uncertain significance (VUS), 37 (36.6%) as being likely pathogenic (LP) and 3 (3%) as being pathogenic (P). Significance: We concluded that the majority of neuropathology-associated variants reported to date do not have enough evidence to be classified as pathogenic. Interpretation of lesion-associated variants remains challenging and application of current ACMG guidelines is recommended for interpretation and prediction.

Genetic generalised epilepsy is the most common type of inherited epilepsy. Despite a high concordance rate of 80% in monozygotic twins, the genetic background is still poorly understood. We aimed to investigate the burden of rare genetic variants in genetic generalised epilepsy.For this exome-based case-control study, we used three different genetic generalised epilepsy case cohorts and three independent control cohorts, all of European descent. Cases included in the study were clinically evaluated for genetic generalised epilepsy. Whole-exome sequencing was done for the discovery case cohort, a validation case cohort, and two independent control cohorts. The replication case cohort underwent targeted next-generation sequencing of the 19 known genes encoding subunits of GABAStatistical comparison of 152 familial index cases with genetic generalised epilepsy in the discovery cohort to 549 ethnically matched controls suggested an enrichment of rare missense (Nonsyn) variants in the ensemble of 19 genes encoding GABAFunctionally relevant variants in genes encoding GABAEuroEPINOMICS (European Science Foundation through national funding organisations), Epicure and EpiPGX (Sixth Framework Programme and Seventh Framework Programme of the European Commission), Research Unit FOR2715 (German Research Foundation and Luxembourg National Research Fund).

Background The use of next-generation sequencing approaches in clinical diagnostics has led to a tremendous increase in data and a vast number of variants of uncertain significance that require interpretation. Therefore, prediction of the effects of missense mutations using in silico tools has become a frequently used approach. Aim of this study was to assess the reliability of in silico prediction as a basis for clinical decision making in the context of hereditary breast and/or ovarian cancer. Methods We tested the performance of four prediction tools (Align-GVGD, SIFT, PolyPhen-2, MutationTaster2) using a set of 236 BRCA1/2 missense variants that had previously been classified by expert committees. However, a major pitfall in the creation of a reliable evaluation set for our purpose is the generally accepted classification of BRCA1/2 missense variants using the multifactorial likelihood model, which is partially based on Align-GVGD results. To overcome this drawback we identified 161 variants whose classification is independent of any previous in silico prediction. In addition to the performance as stand-alone tools we examined the sensitivity, specificity, accuracy and Matthews correlation coefficient (MCC) of combined approaches. Results PolyPhen-2 achieved the lowest sensitivity (0.67), specificity (0.67), accuracy (0.67) and MCC (0.39). Align-GVGD achieved the highest values of specificity (0.92), accuracy (0.92) and MCC (0.73), but was outperformed regarding its sensitivity (0.90) by SIFT (1.00) and MutationTaster2 (1.00). All tools suffered from poor specificities, resulting in an unacceptable proportion of false positive results in a clinical setting. This shortcoming could not be bypassed by combination of these tools. In the best case scenario, 138 families would be affected by the misclassification of neutral variants within the cohort of patients of the German Consortium for Hereditary Breast and Ovarian Cancer. Conclusion We show that due to low specificities state-of-the-art in silico prediction tools are not suitable to predict pathogenicity of variants of uncertain significance in BRCA1/2. Thus, clinical consequences should never be based solely on in silico forecasts. However, our data suggests that SIFT and MutationTaster2 could be suitable to predict benignity, as both tools did not result in false negative predictions in our analysis. Electronic supplementary material The online version of this article (10.1186/s12920-018-0353-y) contains supplementary material, which is available to authorized users.

Individuals with Dupuytren disease (DD) are commonly seen by physicians and surgeons across multiple specialties. It is an increasingly common and disabling fibroproliferative disorder of the palmar fascia, which leads to flexion contractures of the digits, and is associated with other tissue-specific fibroses. DD affects between 5% and 25% of people of European descent and is the most common inherited disease of connective tissue. We undertook the largest GWAS to date in individuals with a surgically validated diagnosis of DD from the UK, with replication in British, Dutch, and German individuals. We validated association at all nine previously described signals and discovered 17 additional variants with p ≤ 5 × 10−8. As a proof of principle, we demonstrated correlation of the high-risk genotype at the statistically most strongly associated variant with decreased secretion of the soluble WNT-antagonist SFRP4, in surgical specimen-derived DD myofibroblasts. These results highlight important pathways involved in the pathogenesis of fibrosis, including WNT signaling, extracellular matrix modulation, and inflammation. In addition, many associated loci contain genes that were hitherto unrecognized as playing a role in fibrosis, opening up new avenues of research that may lead to novel treatments for DD and fibrosis more generally. DD represents an ideal human model disease for fibrosis research.

Male-pattern baldness (MPB) is a common and highly heritable trait characterized by androgen-dependent, progressive hair loss from the scalp. Here, we carry out the largest GWAS meta-analysis of MPB to date, comprising 10,846 early-onset cases and 11,672 controls from eight independent cohorts. We identify 63 MPB-associated loci (P, Male-pattern baldness is a common condition in which hair is progressively lost from the scalp. Here, the authors find 23 new genetic variants associated with this condition and suggest that it is not an isolated trait but may share an underlying biological basis with various diseases.

Genome-wide association studies (GWAS) have revealed genetic determinants of iron metabolism, but correlation of these with clinical phenotypes is pending. Homozygosity for HFE C282Y is the predominant genetic risk factor for hereditary hemochromatosis (HH) and may cause liver cirrhosis. However, this genotype has a low penetrance. Thus, detection of yet unknown genetic markers that identify patients at risk of developing severe liver disease is necessary for better prevention. Genetic loci associated with iron metabolism (TF, TMPRSS6, PCSK7, TFR2 and Chr2p14) in recent GWAS and liver fibrosis (PNPLA3) in recent meta-analysis were analyzed for association with either liver cirrhosis or advanced fibrosis in 148 German HFE C282Y homozygotes. Replication of associations was sought in additional 499 Austrian/Swiss and 112 HFE C282Y homozygotes from Sweden. Only variant rs236918 in the PCSK7 gene (proprotein convertase subtilisin/kexin type 7) was associated with cirrhosis or advanced fibrosis (P = 1.02 × 10(-5)) in the German cohort with genotypic odds ratios of 3.56 (95% CI 1.29-9.77) for CG heterozygotes and 5.38 (95% CI 2.39-12.10) for C allele carriers. Association between rs236918 and cirrhosis was confirmed in Austrian/Swiss HFE C282Y homozygotes (P = 0.014; ORallelic = 1.82 (95% CI 1.12-2.95) but not in Swedish patients. Post hoc combined analyses of German/Swiss/Austrian patients with available liver histology (N = 244, P = 0.00014, ORallelic = 2.84) and of males only (N = 431, P = 2.17 × 10(-5), ORallelic = 2.54) were consistent with the premier finding. Association between rs236918 and cirrhosis was not confirmed in alcoholic cirrhotics, suggesting specificity of this genetic risk factor for HH. PCSK7 variant rs236918 is a risk factor for cirrhosis in HH patients homozygous for the HFE C282Y mutation.

China has repeatedly been the subject of genetic studies to elucidate its prehistoric and historic demography. While some studies reported a genetic distinction between Northern and Southern Han Chinese, others showed a more clinal picture of small differences within China. Here, we investigated the distribution of Y chromosome variation along administrative as well as ethnic divisions in the mainland territory of the People's Republic of China, including 28 administrative regions and 19 recognized Chinese nationalities, to assess the impact of recent demographic processes. To this end, we analyzed 37,994 Y chromosomal 17-marker haplotype profiles from the YHRD database with respect to forensic diversity measures and genetic distance between groups defined by administrative boundaries and ethnic origin. We observed high diversity throughout all Chinese provinces and ethnicities. Some ethnicities, including most prominently Kazakhs and Tibetans, showed significant genetic differentiation from the Han and other groups. However, differences between provinces were, except for those located on the Tibetan plateau, less pronounced. This discrepancy is explicable by the sizeable presence of Han speakers, who showed high genetic homogeneity all across China, in nearly all studied provinces. Furthermore, we observed a continuous genetic North-South gradient in the Han, confirming previous reports of a clinal distribution of Y chromosome variation and being in notable concordance with the previously observed spatial distribution of autosomal variation. Our findings shed light on the demographic changes in China accrued by a fast-growing and increasingly mobile population.

Iran, despite its size, geographic location and past cultural influence, has largely been a blind spot for human population genetic studies. With only sparse genetic information on the Iranian population available, we pursued its genome-wide and geographic characterization based on 1021 samples from eleven ethnic groups. We show that Iranians, while close to neighboring populations, present distinct genetic variation consistent with long-standing genetic continuity, harbor high heterogeneity and different levels of consanguinity, fall apart into a cluster of similar groups and several admixed ones and have experienced numerous language adoption events in the past. Our findings render Iran an important source for human genetic variation in Western and Central Asia, will guide adequate study sampling and assist the interpretation of putative disease-implicated genetic variation. Given Iran’s internal genetic heterogeneity, future studies will have to consider ethnic affiliations and possible admixture., Author summary Based on genome-wide genotype data on over 1000 samples from eleven ethnic groups present in Iran and by comparison to reference data sets of both extant populations and ancient DNA samples, we show that the Iranian population comprises distinct genetic variation with respect to populations in close geographic proximity, a cluster of genetically largely overlapping ethnic groups as well as a number of strongly admixed groups. These observations, also corroborated by f3 migration statistics and other approaches, indicate genetic continuity of and limited influx into the cluster groups over several millennia, despite Iran’s geographic position at a crossroads in West Asia. They also suggest, correspondingly, several instances of language adoption instead of demic replacement in the past. Future human genetic studies, both with a focus on population and medical genetics, will have to consider differences in heterogeneity, consanguinity and degree of admixture between the ethnic groups for an adequate design and interpretation.

Background & Aims Genome-wide association studies (GWAS) have identified 140 Crohn's disease (CD) susceptibility loci. For most loci, the variants that cause disease are not known and the genes affected by these variants have not been identified. We aimed to identify variants that cause CD through detailed sequencing, genetic association, expression, and functional studies. Methods We sequenced whole exomes of 42 unrelated subjects with CD and 5 healthy subjects (controls) and then filtered single nucleotide variants by incorporating association results from meta-analyses of CD GWAS and in silico mutation effect prediction algorithms. We then genotyped 9348 subjects with CD, 2868 subjects with ulcerative colitis, and 14,567 control subjects and associated variants analyzed in functional studies using materials from subjects and controls and in vitro model systems. Results We identified rare missense mutations in PR domain-containing 1 ( PRDM1 ) and associated these with CD. These mutations increased proliferation of T cells and secretion of cytokines on activation and increased expression of the adhesion molecule L-selectin. A common CD risk allele, identified in GWAS, correlated with reduced expression of PRDM1 in ileal biopsy specimens and peripheral blood mononuclear cells (combined P = 1.6 × 10 −8 ). We identified an association between CD and a common missense variant, Val248Ala, in nuclear domain 10 protein 52 ( NDP52 ) ( P = 4.83 × 10 −9 ). We found that this variant impairs the regulatory functions of NDP52 to inhibit nuclear factor κB activation of genes that regulate inflammation and affect the stability of proteins in Toll-like receptor pathways. Conclusions We have extended the results of GWAS and provide evidence that variants in PRDM1 and NDP52 determine susceptibility to CD. PRDM1 maps adjacent to a CD interval identified in GWAS and encodes a transcription factor expressed by T and B cells. NDP52 is an adaptor protein that functions in selective autophagy of intracellular bacteria and signaling molecules, supporting the role of autophagy in the pathogenesis of CD.

Antigen presentation in chronic skin disorders is mediated through the interleukin (IL)-12/IL-23 pathway and, hence, through the IL-12 receptor. Recent evidence suggesting dysregulated antigen presentation in skin lesions of hidradenitis suppurativa (HS) led to investigate the role of single nucleotide polymorphisms (SNPs) of the gene IL-12RB1 coding for the IL12-Rβ1 receptor subunit. Genomic DNA was isolated from 139 patients and 113 healthy controls; nine SNPs in the transcribed region of IL12RB1 were genotyped. No significant differences of genotype and allele frequencies were found between the two groups. Two common haplotypes were recognized, namely h1 and h2. Carriage of h2 related with minor frequency alleles was associated with a greater risk for the acquisition of Hurley III disease stage and with the involvement of a greater number of skin areas. Patients with the h1 haplotype presented disease at an older age. This is the genetic study enrolling the largest number of patients with HS to date. Although SNPs of IL12RB1 do not seem to convey genetic predisposition, they are associated directly with the phenotype of the disease.

BACKGROUND/OBJECTIVES: Fatty liver disease (FLD) is an important intermediate trait along the cardiometabolic disease spectrum and strongly associates with type 2 diabetes. Knowledge of biological pathways implicated in FLD is limited. An untargeted metabolomic approach might unravel novel pathways related to FLD. SUBJECTS/METHODS: In a population-based sample (n=555) from Northern Germany, liver fat content was quantified as liver signal intensity using magnetic resonance imaging. Serum metabolites were determined using a non-targeted approach. Partial least squares regression was applied to derive a metabolomic score, explaining variation in serum metabolites and liver signal intensity. Associations of the metabolomic score with liver signal intensity and FLD were investigated in multivariable-adjusted robust linear and logistic regression models, respectively. Metabolites with a variable importance in the projection >1 were entered in in silico overrepresentation and pathway analyses. RESULTS: In univariate analysis, the metabolomics score explained 23.9% variation in liver signal intensity. A 1-unit increment in the metabolomic score was positively associated with FLD (n=219; odds ratio: 1.36; 95% confidence interval: 1.27-1.45) adjusting for age, sex, education, smoking and physical activity. A simplified score based on the 15 metabolites with highest variable importance in the projection statistic showed similar associations. Overrepresentation and pathway analyses highlighted branched-chain amino acids and derived gamma-glutamyl dipeptides as significant correlates of FLD. CONCLUSIONS: A serum metabolomic profile was associated with FLD and liver fat content. We identified a simplified metabolomics score, which should be evaluated in prospective studies.European Journal of Clinical Nutrition advance online publication, 5 April 2017; doi:10.1038/ejcn.2017.43.

Dupuytren's disease, a fibromatosis of the connective tissue in the palm, is a common complex disease with a strong genetic component. Up to date nine genetic loci have been found to be associated with the disease. Six of these loci contain genes that code for Wnt signalling proteins. In spite of this striking first insight into the genetic factors in Dupuytren's disease, much of the inherited risk in Dupuytren's disease still needs to be discovered. The already identified loci jointly explain similar to 1% of the heritability in this disease. To further elucidate the genetic basis of Dupuytren's disease, we performed a genome-wide meta-analysis combining three genome-wide association study (GWAS) data sets, comprising 1,580 cases and 4,480 controls. We corroborated all nine previously identified loci, six of these with genome-wide significance (p-value

Typing of X chromosomal short tandem repeat (X STR) markers has become a standard element of human forensic genetic analysis. Joint consideration of many X STR markers at a time increases their discriminatory power but, owing to physical linkage, requires inter-marker recombination rates to be accurately known. We estimated the recombination rates between 15 well established X STR markers using genotype data from 158 families (1041 individuals) and following a previously proposed likelihood-based approach that allows for single-step mutations. To meet the computational requirements of this family-based type of analysis, we modified a previous implementation so as to allow multi-core parallelization on a high-performance computing system. While we obtained recombination rate estimates larger than zero for all but one pair of adjacent markers within the four previously proposed linkage groups, none of the three X STR pairs defining the junctions of these groups yielded a recombination rate estimate of 0.50. Corroborating previous studies, our results therefore argue against a simple model of independent X chromosomal linkage groups. Moreover, the refined recombination fraction estimates obtained in our study will facilitate the appropriate joint consideration of all 15 investigated markers in forensic analysis.

The vast majority of patients with Nijmegen Breakage Syndrome (NBS) are of Slavic origin and carry a deleterious deletion (c.657del5; rs587776650) in the NBN gene on chromosome 8q21. This mutation is essentially confined to Slavic populations and may thus be considered a Slavic founder mutation. Notably, not a single parenthood of a homozygous c.657del5 carrier has been reported to date, while heterozygous carriers do reproduce but have an increased cancer risk. These observations seem to conflict with the considerable carrier frequency of c.657del5 of 0.5% to 1% as observed in different Slavic populations because deleterious mutations would be eliminated quite rapidly by purifying selection. Therefore, we propose that heterozygous c.657del5 carriers have increased reproductive success, i.e., that the mutation confers heterozygote advantage. In fact, in our cohort study of the reproductive history of 24 NBS pedigrees from the Czech Republic, we observed that female carriers gave birth to more children on average than female non-carriers, while no such reproductive differences were observed for males. We also estimate that c.657del5 likely occurred less than 300 generations ago, thus supporting the view that the original mutation predated the historic split and subsequent spread of the 'Slavic people'. We surmise that the higher fertility of female c.657del5 carriers reflects a lower miscarriage rate in these women, thereby reflecting the role of the NBN gene product, nibrin, in the repair of DNA double strand breaks and their processing in immune gene rearrangements, telomere maintenance, and meiotic recombination, akin to the previously described role of the DNA repair genes BRCA1 and BRCA2.

Rationale: Sarcoidosis is a complex inflammatory disease with a heterogeneous clinical picture. Among others, an acute and chronic clinical course can be distinguished, for which specific genetic risk factors are known. Objectives: To identify additional risk loci for sarcoidosis and its acute and chronic subforms, we analyzed imputed data from a genomewide association scan for these phenotypes. Methods: After quality control, the genome-wide association scan comprised nearly 1.3 million imputed single-nucleotide polymorphisms based on an Affymetrix 6.0 Gene Chip dataset of 564 German sarcoidosis cases, including 176 acute and 354 chronic cases and 1,575 control subjects. Measurements and Main Results: We identified chromosome 11q13.1 (rs479777) as a novel locus influencing susceptibility to sarcoidosis with genome-wide significance. The marker was significantly associated in three distinct German case-control populations and in an additional German family sample with odds ratios ranging from 0.67 to 0.77. This finding was further replicated in two independent European case-control populations from the Czech Republic (odds ratio, 0.75) and from Sweden (odds ratio, 0.79). In a meta-analysis of the included European case-control samples the marker yielded a P value of 2.68 x 10(-18). The locus was previously reported to be associated with Crohn disease, psoriasis, alopecia areata, and leprosy. For sarcoidosis, fine-mapping and expression analysis suggest KCNK4, PRDX5, PCLB3, and most promising CCDC88B as candidates for the underlying risk gene in the associated region. Conclusions: This study provides striking evidence for association of chromosome 11q13.1 with sarcoidosis in Europeans, and thus identified a further genetic risk locus shared by sarcoidosis, Crohn disease and psoriasis.

So far, the annotation of translation initiation sites (TISs) has been based mostly upon bioinformatics rather than experimental evidence. We adapted ribosomal footprinting to puromycin-treated cells to generate a transcriptome-wide map of TISs in a human monocytic cell line. A neural network was trained on the ribosomal footprints observed at previously annotated AUG translation initiation codons (TICs), and used for the ab initio prediction of TISs in 5062 transcripts with sufficient sequence coverage. Functional interpretation suggested 2994 novel upstream open reading frames (uORFs) in the 5′ UTR, 1406 uORFs overlapping with the coding sequence, and 546 N-terminal protein extensions. The TIS detection method was validated on the basis of previously published alternative TISs and uORFs. Among primates, TICs in newly annotated TISs were significantly more conserved than control codons, both for AUGs and near-cognate codons. The transcriptome-wide map of novel candidate TISs derived as part of the study will shed further light on the way in which human proteome diversity is influenced by alternative translation initiation and regulation.

Objective: Sporadic, genetically complex essential tremor (ET) is one of the most common movement disorders and may lead to severe impairment of the quality of life. Despite high heritability, the genetic determinants of ET are largely unknown. We performed the second genome-wide association study (GWAS) for ET to elucidate genetic risk factors of ET. Methods: Using the Affymetrix Genome-Wide SNP Array 6.0 (1000K) we conducted a two-stage GWAS in a total of 990 subjects and 1,537 control subjects from Europe to identify genetic variants associated with ET. Results: We discovered association of an intronic variant of the main glial glutamate transporter ( SLC1A2 ) gene with ET in the first-stage sample (rs3794087, p = 6.95 × 10 −5 , odds ratio [OR] = 1.46). We verified the association of rs3794087 with ET in a second-stage sample ( p = 1.25 × 10 −3 , OR = 1.38). In the subgroup analysis of patients classified as definite ET, rs3794087 obtained genome-wide significance ( p = 3.44 × 10 −10 , OR = 1.59) in the combined first- and second-stage sample. Genetic fine mapping using nonsynonymous single nucleotide polymorphisms (SNPs) and SNPs in high linkage disequilibrium with rs3794087 did not reveal any SNP with a stronger association with ET than rs3794087. Conclusions: We identified SLC1A2 encoding the major glial high-affinity glutamate reuptake transporter in the brain as a potential ET susceptibility gene. Acute and chronic glutamatergic overexcitation is implied in the pathogenesis of ET. SLC1A2 is therefore a good functional candidate gene for ET.

Psoriasis (PS) and Crohn disease (CD) have been shown to be epidemiologically, pathologically, and therapeutically connected, but little is known about their shared genetic causes. We performed meta-analyses of five published genome-wide association studies on PS (2,529 cases and 4,955 controls) and CD (2,142 cases and 5,505 controls), followed up 20 loci that showed strongest evidence for shared disease association and, furthermore, tested cross-disease associations for previously reported PS and CD risk alleles in additional6,115 PS cases, 4,073 CD cases, and 10,100 controls. We identified seven susceptibility loci outside the human leukocyte antigen region (9p24 near JAK2, 10q22 at ZMIZ1, 11q13 near PRDX5, 16p13 near SOCS1, 17q21 at STAT3, 19p13 near FUT2, and 22q11 at YDJC) shared between PS and CD with genome-wide significance (p < 5× 10(-8)) and confirmed four already established PS and CD risk loci (IL23R, IL12B, REL, and TYK2). Three of the shared loci are also genome-wide significantly associated with PS alone (10q22 at ZMIZ1, p(rs1250544) = 3.53×10(-8), 11q13 near PRDX5, p(rs694739) = 3.71× 10(-09), 22q11 at YDJC, p(rs181359) = 8.02× 10(-10)). In addition, we identified one susceptibility locus for CD (16p13 near SOCS1, p(rs4780355) = 4.99× 10(-8)). Refinement of association signals identified shared genome-wide significant associations for exonic SNPs at 10q22 (ZMIZ1) and insilico expression quantitative trait locus analyses revealed that the associations at ZMIZ1 and near SOCS1 have a potential functional effect on gene expression. Our results show the usefulness of joint analyses of clinically distinct immune-mediated diseases and enlarge the map of shared genetic risk loci.

Several polymorphisms of the transcription factor 4 ( TCF4 ) have been shown to increase the risk for schizophrenia, particularly TCF4 rs9960767. This polymorphism is associated with impaired sensorimotor gating measured by prepulse inhibition—an established endophenotype of schizophrenia. We therefore investigated whether TCF4 polymorphisms also affect another proposed endophenotype of schizophrenia, namely sensory gating assessed by P50 suppression of the auditory evoked potential. Although sensorimotor gating and sensory gating are not identical, recent data suggest that they share genetic fundamentals. In a multicenter study at six academic institutions throughout Germany, we applied an auditory P50 suppression paradigm to 1,821 subjects (1,023 never-smokers, 798 smokers) randomly selected from the general population. Samples were genotyped for 21 TCF4 polymorphisms. Given that smoking is highly prevalent in schizophrenia and affects sensory gating, we also assessed smoking behavior, cotinine plasma concentrations, exhaled carbon monoxide, and the Fagerström Test (FTND). P50 suppression was significantly decreased in carriers of schizophrenia risk alleles of the TCF4 polymorphisms rs9960767, rs10401120rs, rs17597926, and 17512836 ( P < 0.0002–0.00005). These gene effects were modulated by smoking behavior as indicated by significant interactions of TCF4 genotype and smoking status; heavy smokers (FTND score ≥4) showed stronger gene effects on P50 suppression than light smokers and never-smokers. Our finding suggests that sensory gating is modulated by an interaction of TCF4 genotype with smoking, and both factors may play a role in early information processing deficits also in schizophrenia. Consequently, considering smoking behavior may facilitate the search for genetic risk factors for schizophrenia.

AimInvolvement of TLR2 in the pathophysiology of periodontitis has widely been discussed, but hitherto, no validated genetic associations were reported. Previous association studies lacked sufficient statistical power and adequate haplotype information to draw unambiguous conclusions. The aim of this study was to comprehensively investigate TLR2 linkage disequilibrium (LD) regions for their potential associations with periodontitis in two large analysis populations of aggressive (AgP) and chronic periodontitis (CP) of North West European descent.Materials and methodsThe study population comprised 598 AgP patients, 914 CP patients and 1804 healthy controls. Analysis of TLR2 LD regions was performed with haplotype tagging SNPs (tagSNPs) using SNPlex and TaqMan genotyping assays. Genotypic, dominant, multiplicative, and recessive genetic models were tested. The genotypes were adjusted for the covariates smoking, diabetes, and gender. Resequencing was performed by Sanger technology.ResultsUpon covariate adjustment and correction for multiple testing, no tagSNPs showed significant associations with AgP or CP. Targeted resequencing of exon 3 in 47 AgP cases identified carriership of two common and three rare variants.ConclusionCommon LD regions of TLR2 do not show genetic associations with periodontitis in the North West European population. Resequencing of exon 3 could not identify disease-associated rare variants in TLR2.

Alcohol misuse is the leading cause of cirrhosis and the second most common indication for liver transplantation in the Western world. We performed a genome-wide association study for alcohol-related cirrhosis in individuals of European descent (712 cases and 1,426 controls) with subsequent validation in two independent European cohorts (1,148 cases and 922 controls). We identified variants in the MBOAT7 (P = 1.03 × 10(-9)) and TM6SF2 (P = 7.89 × 10(-10)) genes as new risk loci and confirmed rs738409 in PNPLA3 as an important risk locus for alcohol-related cirrhosis (P = 1.54 × 10(-48)) at a genome-wide level of significance. These three loci have a role in lipid processing, suggesting that lipid turnover is important in the pathogenesis of alcohol-related cirrhosis.

We conducted a case-control genome-wide association study (GWAS) of human longevity, comparing 664,472 autosomal SNPs in 763 long-lived individuals (LLI; mean age: 99.7 years) and 1085 controls (mean age: 60.2 years) from Germany. Only one association, namely that of SNP rs4420638 near the APOC1 gene, achieved genome-wide significance (allele-based P = 1.8 x 10(-10)). However, logistic regression analysis revealed that this association, which was replicated in an independent German sample, is fully explicable by linkage disequilibrium with the APOE allele epsilon 4, the only variant hitherto established as a major genetic determinant of survival into old age. Our GWAS failed to identify any additional autosomal susceptibility genes. One explanation for this lack of success in our study would be that GWAS provide only limited statistical power for a polygenic phenotype with loci of small effect such as human longevity. A recent GWAS in Dutch LLI independently confirmed the APOE-longevity association, thus strengthening the conclusion that this locus is a very, if not the most, important genetic factor influencing longevity. (C) 2011 Elsevier Ireland Ltd. All rights reserved.

Epidemiological studies have indicated a relationship between coronary heart disease (CHD) and periodontitis. Recently, CDKN2BAS was reported as a shared genetic risk factor of CHD and aggressive periodontitis (AgP), but the causative variant has remained unknown. To identify and validate risk variants in different European populations, we first explored 150 kb of the genetic region of CDKN2BAS including the adjacent genes CDKN2A and CDKN2B, covering 51 tagging single nucleotide polymorphisms (tagSNPs) in AgP and chronic periodontitis (CP) in individuals of Dutch origin (n=313). In a second step, we tested the significant SNP associations in an independent AgP and CP population of German origin (n=1264). For the tagSNPs rs1360590, rs3217992, and rs518394, we could validate the associations with AgP before and after adjustment for the covariates smoking, gender and diabetes, with SNP rs3217992 being the most significant (OR 1.48, 95% CI 1.19 to 1.85; p=0.0004). We further showed in vivo gene expression of CDKN2BAS, CDKN2A, CDKN2B, and CDK4 in healthy and inflamed gingival epithelium (GE) and connective tissue (CT), and detected a significantly higher expression of CDKN2BAS in healthy CT compared to GE (p=0.004). After 24 h of stimulation with Porphyromonas gingivalis in Streptococcus gordonii pre-treated gingival fibroblast (HGF) and cultured gingival epithelial cells (GECs), we observed a 25-fold and fourfold increase of CDKN2BAS gene expression in HGFs (p=0.003) and GECs (p=0.004), respectively. Considering the global importance of CDKN2BAS in the disease risk of CHD, this observation supports the theory of inflammatory components in the disease physiology of CHD.

Sarcoidosis is a complex granulomatous inflammatory disorder that shares several clinical and pathogenic features with inflammatory bowel disease (IBD). Postulating a common genetic basis of inflammatory diseases, we tested 106 single-nucleotide polymorphisms (SNPs) that are known or have been suggested to be associated with IBD for a potential association with sarcoidosis and its acute and chronic subphenotypes. We genotyped 1,996 German sarcoidosis patients, comprising 648 acutely and 1,161 chronically affected individuals, 2,622 control subjects, and 342 German trios with affected offspring using SNPlex™ technology. The nonsynonymous SNP rs11209026 (Arg381Gln) in the interleukin (IL)-23 receptor ( IL23R ) gene was associated with chronic sarcoidosis (OR 0.63; p = 5.58×10 −5 ), which was supported by the result of a transmission disequilibrium test analysis in the independent family sample (OR 0.50; p = 0.031). Marker rs12035082 located at chromosome 1q24.3 was found to be associated with the acute subphenotype (OR 1.36; p = 6.80×10 −7 ) and rs916977 ( HERC2 locus; OR 1.30; p = 4.49×10 −5 ) was associated with sarcoidosis. Our results highlight the potential importance of the IL-23 signalling pathway for the development of chronic sarcoidosis. The finding links sarcoidosis pathogenesis to other inflammatory conditions and may contribute to new hypotheses on disease mechanisms.

The availability of high-density panels of genetic polymorphisms has led to the discovery of extended regions of apparent autozygosity in the human genome. At the genotype level, these regions present as sizeable stretches, or 'runs', of homozygosity (ROH). Here, we investigated both the genomic and the geographic distribution of ROHs in a large European sample of individuals originating from 23 subpopulations. The genomic ROH distribution was found to be characterized by a pattern of highly significant non-uniformity that was virtually identical in all subpopulations studied. Some 77 chromosomal regions contained ROHs at considerable frequency, thereby forming 'ROH islands' that were not explicable by high linkage disequilibrium alone. At the geographic level, the number and cumulative length of ROHs followed a prominent South to North gradient in agreement with expectations from European population history. The individual ROH length, in contrast, showed only minor and unsystematic geographic variation. While our findings are thus consistent with a larger effective population size in Southern than in Northern Europe, combined with a higher historic population density and mobility, they also indicate that the patterns of meiotic recombination in humans must have been very similar throughout the continent. Extending previous reports of a strong correlation between geography and identity-by-state, our data show that the genomic identity-by-descent patterns of Europeans are also clinal. As a consequence, the planning, design and interpretation of ROH-based genetic studies must take sample origin into account in order for such studies to be sensible and valid.

Essential tremor (ET) is one of the most common movement disorders. Former association studies focussing on candidate genes in ET found a number of risk variants but most of them were not replicated. Recently, a genome-wide association study revealed two intronic sequence variants in the LINGO1 gene associated with ET. Here, we have confirmed association between sequence variants in the LINGO1 gene and the ET phenotype in independent German and French ET samples. The odds ratios for the identified intronic markers rs8030859 (P = 1.0x10(-4)), rs9652490 (P = 9.1x10(-4)), and rs11856808 (P = 3.6x10(-2)) were 1.72 (CI 1.31-2.26), 1.61 (CI 1.21-2.14), and 1.30 (CI 1.02-1.66), respectively, in our German sample. LINGO1 is an interesting candidate gene because it plays a key role in central nervous system biology, is selectively expressed in the nervous system, and is an inhibitor of oligodendrocyte differentiation and neuronal myelination. Our study gives further evidence that LINGO1 acts as a susceptibility gene for ET.

Background & Aims We aimed to characterize the genetic susceptibility to primary sclerosing cholangitis (PSC) by means of a genome-wide association analysis of single nucleotide polymorphism (SNP) markers. Methods A total of 443,816 SNPs on the Affymetrix SNP Array 5.0 (Affymetrix, Santa Clara, CA) were genotyped in 285 Norwegian PSC patients and 298 healthy controls. Associations detected in this discovery panel were re-examined in independent case-control panels from Scandinavia (137 PSC cases and 368 controls), Belgium/The Netherlands (229 PSC cases and 735 controls), and Germany (400 cases and 1832 controls). Results The strongest associations were detected near HLA-B at chromosome 6p21 (rs3099844: odds ratio [OR], 4.8; 95% confidence interval [CI], 3.6–6.5; P = 2.6 × 10 −26 ; and rs2844559: OR, 4.7; 95% CI, 3.5–6.4; P = 4.2 × 10 −26 in the discovery panel). Outside the HLA complex, rs9524260 at chromosome 13q31 showed significant associations in 3 of 4 study panels. Lentiviral silencing of glypican 6, encoded at this locus, led to the up-regulation of proinflammatory markers in a cholangiocyte cell line. Of 15 established ulcerative colitis susceptibility loci, significant replication was obtained at chromosomes 2q35 and 3p21 (rs12612347: OR, 1.26; 95% CI, 1.06–1.50; and rs3197999: OR, 1.22; 95% CI, 1.02–1.47, respectively), with circumstantial evidence supporting the G-protein–coupled bile acid receptor 1 and macrophage-stimulating 1 , respectively, as the likely disease genes. Conclusions Strong HLA associations and a subset of genes involved in bile homeostasis and other inflammatory conditions constitute key components of the genetic architecture of PSC.

Recent studies indicate a mutual epidemiological relationship between coronary heart disease (CHD) and periodontitis. Both diseases are associated with similar risk factors and are characterized by a chronic inflammatory process. In a candidate-gene association study, we identify an association of a genetic susceptibility locus shared by both diseases. We confirm the known association of two neighboring linkage disequilibrium regions on human chromosome 9p21.3 with CHD and show the additional strong association of these loci with the risk of aggressive periodontitis. For the lead SNP of the main associated linkage disequilibrium region, rs1333048, the odds ratio of the autosomal-recessive mode of inheritance is 1.99 (95% confidence interval 1.33–2.94; P = 6.9×10−4) for generalized aggressive periodontitis, and 1.72 (1.06–2.76; P = 2.6×10−2) for localized aggressive periodontitis. The two associated linkage disequilibrium regions map to the sequence of the large antisense noncoding RNA ANRIL, which partly overlaps regulatory and coding sequences of CDKN2A/CDKN2B. A closely located diabetes-associated variant was independent of the CHD and periodontitis risk haplotypes. Our study demonstrates that CHD and periodontitis are genetically related by at least one susceptibility locus, which is possibly involved in ANRIL activity and independent of diabetes associated risk variants within this region. Elucidation of the interplay of ANRIL transcript variants and their involvement in increased susceptibility to the interactive diseases CHD and periodontitis promises new insight into the underlying shared pathogenic mechanisms of these complex common diseases., Author Summary Coronary heart disease (CHD) and periodontitis are the most widespread diseases in the Western industrialized world and pose a substantial health threat to populations worldwide. CHD is a leading cause for premature death, and periodontitis is the major cause for tooth loss in adults over 40 years. Both diseases are associated with similar risk factors such as smoking, diabetes, and gender, and both diseases are further characterized by a chronic inflammatory process. In the last year, several genome studies have identified a region of the human genome near the CDKN2A and CDKN2B genes as having an influence on CHD. We show that this genetic region, being the most important susceptibility locus for CHD to date, is also associated with a substantial risk increase of aggressive periodontitis. The associated genetic region maps to a genomic region that codes for an “antisense RNA,” which partly overlaps regulatory and coding sequences of genes CDKN2A/CDKN2B. The interplay between these common inflammatory complex diseases could be partially due to the shared genetic risk variants of this antisense RNA.

Genome-wide association studies have contributed significantly to the genetic dissection of complex diseases. In order to increase the power of existing marker sets even further, methods have been proposed to predict individual genotypes at un-typed loci from other marker sets by imputation, usually employing HapMap data as a reference. Although various imputation algorithms have been used in practice already, a comprehensive evaluation and comparison of these approaches, using genome-wide SNP data from one and the same population is still lacking. We therefore investigated four publicly available programs for genotype imputation (BEAGLE, IMPUTE, MACH, and PLINK) using data from 449 German individuals genotyped in our laboratory for three genome-wide SNP sets [Affymetrix 5.0 (500 k), Affymetrix 6.0 (1,000 k), and Illumina 550 k]. We observed that HapMap-based imputation in a northern European population is powerful and reliable, even in highly variable genomic regions such as the extended MHC on chromosome 6p21. However, while genotype predictions were found to be highly accurate with all four programs, the number of SNPs for which imputation was actually carried out ('imputation efficacy') varied substantially. BEAGLE, IMPUTE, and MACH yielded nearly identical trade-offs between imputation accuracy and efficacy whereas PLINK performed consistently poorer. We nevertheless recommend either MACH or BEAGLE for practical use because these two programs are more user-friendly and generally require less memory than IMPUTE.