Genomic ancestry and social dynamics of the last hunter-gatherers of Atlantic France.

Since the early Holocene, western and central Europe was inhabited by a genetically distinct group of Western Hunter-Gatherers (WHGs). This group was eventually replaced and assimilated by the incoming Neolithic farmers. The western Atlantic façade was home to some of the last Mesolithic sites of mainland Europe, represented by the iconic open-air sites at Hoedic and Téviec in southern Brittany, France. These sites are known for the unusually well-preserved and rich burials. Genomic studies of Mesolithic European hunter-gatherers have been limited to single or a few individuals per site and our understanding of the social dynamics of the last Mesolithic hunter-gatherers of Europe and their interactions with incoming farmers is limited. We sequenced and analyzed the complete genomes of 10 individuals from the Late Mesolithic sites of Hoedic, Téviec, and Champigny, in France, four of which sequenced to between 23- and 8-times genome coverage. The analysis of genomic, chronological and dietary data revealed that the Late Mesolithic populations in Brittany maintained distinct social units within a network of exchanging mates. This resulted in low intra-group biological relatedness that prevented consanguineous mating, despite the small population size of the Late Mesolithic groups. We found no genetic ancestry from Neolithic farmers in the analyzed hunter-gatherers, even though some of them may have coexisted with the first farming groups in neighboring regions. Hence, contrary to previous conclusions based on stable isotope data from the same sites, the Late Mesolithic forager community was limited in mate-exchange to neighboring hunter-gatherer groups, to the exclusion of Neolithic farmers.

Brugada syndrome in Japan and Europe: a genome-wide association study reveals shared genetic architecture and new risk loci.

BACKGROUND AND AIMS: Brugada syndrome (BrS) is an inherited arrhythmia with a higher disease prevalence and more lethal arrhythmic events in Asians than in Europeans. Genome-wide association studies (GWAS) have revealed its polygenic architecture mainly in European populations. The aim of this study was to identify novel BrS-associated loci and to compare allelic effects across ancestries. METHODS: A GWAS was conducted in Japanese participants, involving 940 cases and 1634 controls, followed by a cross-ancestry meta-analysis of Japanese and European GWAS (total of 3760 cases and 11 635 controls). The novel loci were characterized by fine-mapping, gene expression, and splicing quantitative trait associations in the human heart. RESULTS: The Japanese-specific GWAS identified one novel locus near ZSCAN20 (P = 1.0 × 10-8), and the cross-ancestry meta-analysis identified 17 association signals, including six novel loci. The effect directions of the 17 lead variants were consistent (94.1%; P for sign test = 2.7 × 10-4), and their allelic effects were highly correlated across ancestries (Pearson's R = .91; P = 2.9 × 10-7). The genetic risk score derived from the BrS GWAS of European ancestry was significantly associated with the risk of BrS in the Japanese population [odds ratio 2.12 (95% confidence interval 1.94-2.31); P = 1.2 × 10-61], suggesting a shared genetic architecture across ancestries. Functional characterization revealed that a lead variant in CAMK2D promotes alternative splicing, resulting in an isoform switch of calmodulin kinase II-δ, favouring a pro-inflammatory/pro-death pathway. CONCLUSIONS: This study demonstrates novel susceptibility loci implicating potentially novel pathogenesis underlying BrS. Despite differences in clinical expressivity and epidemiology, the polygenic architecture of BrS was substantially shared across ancestries.

Dyslipidemia, inflammation, calcification, and adiposity in aortic stenosis: a genome-wide study.

AIMS: Although highly heritable, the genetic etiology of calcific aortic stenosis (AS) remains incompletely understood. The aim of this study was to discover novel genetic contributors to AS and to integrate functional, expression, and cross-phenotype data to identify mechanisms of AS. METHODS AND RESULTS: A genome-wide meta-analysis of 11.6 million variants in 10 cohorts involving 653 867 European ancestry participants (13 765 cases) was performed. Seventeen loci were associated with AS at P ≤ 5 × 10-8, of which 15 replicated in an independent cohort of 90 828 participants (7111 cases), including CELSR2-SORT1, NLRP6, and SMC2. A genetic risk score comprised of the index variants was associated with AS [odds ratio (OR) per standard deviation, 1.31; 95% confidence interval (CI), 1.26-1.35; P = 2.7 × 10-51] and aortic valve calcium (OR per standard deviation, 1.22; 95% CI, 1.08-1.37; P = 1.4 × 10-3), after adjustment for known risk factors. A phenome-wide association study indicated multiple associations with coronary artery disease, apolipoprotein B, and triglycerides. Mendelian randomization supported a causal role for apolipoprotein B-containing lipoprotein particles in AS (OR per g/L of apolipoprotein B, 3.85; 95% CI, 2.90-5.12; P = 2.1 × 10-20) and replicated previous findings of causality for lipoprotein(a) (OR per natural logarithm, 1.20; 95% CI, 1.17-1.23; P = 4.8 × 10-73) and body mass index (OR per kg/m2, 1.07; 95% CI, 1.05-1.9; P = 1.9 × 10-12). Colocalization analyses using the GTEx database identified a role for differential expression of the genes LPA, SORT1, ACTR2, NOTCH4, IL6R, and FADS. CONCLUSION: Dyslipidemia, inflammation, calcification, and adiposity play important roles in the etiology of AS, implicating novel treatments and prevention strategies.

Genome-wide association study reveals novel genetic loci: a new polygenic risk score for mitral valve prolapse.

AIMS: Mitral valve prolapse (MVP) is a common valvular heart disease with a prevalence of >2% in the general adult population. Despite this high incidence, there is a limited understanding of the molecular mechanism of this disease, and no medical therapy is available for this disease. We aimed to elucidate the genetic basis of MVP in order to better understand this complex disorder. METHODS AND RESULTS: We performed a meta-analysis of six genome-wide association studies that included 4884 cases and 434 649 controls. We identified 14 loci associated with MVP in our primary analysis and 2 additional loci associated with a subset of the samples that additionally underwent mitral valve surgery. Integration of epigenetic, transcriptional, and proteomic data identified candidate MVP genes including LMCD1, SPTBN1, LTBP2, TGFB2, NMB, and ALPK3. We created a polygenic risk score (PRS) for MVP and showed an improved MVP risk prediction beyond age, sex, and clinical risk factors. CONCLUSION: We identified 14 genetic loci that are associated with MVP. Multiple analyses identified candidate genes including two transforming growth factor-ß signalling molecules and spectrin ß. We present the first PRS for MVP that could eventually aid risk stratification of patients for MVP screening in a clinical setting. These findings advance our understanding of this common valvular heart disease and may reveal novel therapeutic targets for intervention.

Burden of rare variants in arrhythmogenic cardiomyopathy with right dominant form-associated genes provides new insights for molecular diagnosis and clinical management.

Arrhythmogenic cardiomyopathy with right dominant form (ACR) is a rare heritable cardiac cardiomyopathy disorder associated with sudden cardiac death. Pathogenic variants (PVs) in desmosomal genes have been causally related to ACR in 40% of cases. Other genes encoding nondesmosomal proteins have been described in ACR, but their contribution in this pathology is still debated. A panel of 71 genes associated with inherited cardiopathies was screened in an ACR population of 172 probands and 856 individuals from the general population. PVs and uncertain significance variants (VUS) have been identified in 36% and 18.6% of patients, respectively. Among the cardiopathy-associated genes, burden tests show a significant enrichment in PV and VUS only for desmosomal genes PKP2 (plakophilin-2), DSP (desmoplakin), DSC2 (desmocollin-2), and DSG2 (desmoglein-2). Importantly, VUS may account for 15% of ACR cases and should then be considered for molecular diagnosis. Among the other genes, no evidence of enrichment was detected, suggesting an extreme caution in the interpretation of these genetic variations without associated functional or segregation data. Genotype-phenotype correlation points to (1) a more severe and earlier onset of the disease in PV and VUS carriers, underlying the importance to carry out presymptomatic diagnosis in relatives and (2) to a more prevalent left ventricular dysfunction in DSP variant carriers.

Correction: AVPR1a and SLC6A4 Gene Polymorphisms Are Associated with Creative Dance Performance.

[This corrects the article DOI: 10.1371/journal.pgen.0010042.].

Association of Coding Variants in Hydroxysteroid 17-beta Dehydrogenase 14 (HSD17B14) with Reduced Progression to End Stage Kidney Disease in Type 1 Diabetes.

BACKGROUND: Rare variants in gene coding regions likely have a greater impact on disease-related phenotypes than common variants through disruption of their encoded protein. We searched for rare variants associated with onset of ESKD in individuals with type 1 diabetes at advanced kidney disease stage. METHODS: Gene-based exome array analyses of 15,449 genes in five large incidence cohorts of individuals with type 1 diabetes and proteinuria were analyzed for survival time to ESKD, testing the top gene in a sixth cohort (n=2372/1115 events all cohorts) and replicating in two retrospective case-control studies (n=1072 cases, 752 controls). Deep resequencing of the top associated gene in five cohorts confirmed the findings. We performed immunohistochemistry and gene expression experiments in human control and diseased cells, and in mouse ischemia reperfusion and aristolochic acid nephropathy models. RESULTS: Protein coding variants in the hydroxysteroid 17-ß dehydrogenase 14 gene (HSD17B14), predicted to affect protein structure, had a net protective effect against development of ESKD at exome-wide significance (n=4196; P value=3.3 × 10-7). The HSD17B14 gene and encoded enzyme were robustly expressed in healthy human kidney, maximally in proximal tubular cells. Paradoxically, gene and protein expression were attenuated in human diabetic proximal tubules and in mouse kidney injury models. Expressed HSD17B14 gene and protein levels remained low without recovery after 21 days in a murine ischemic reperfusion injury model. Decreased gene expression was found in other CKD-associated renal pathologies. CONCLUSIONS: HSD17B14 gene is mechanistically involved in diabetic kidney disease. The encoded sex steroid enzyme is a druggable target, potentially opening a new avenue for therapeutic development.

Functionally validated SCN5A variants allow interpretation of pathogenicity and prediction of lethal events in Brugada syndrome.

AIMS: The prognostic value of genetic variants for predicting lethal arrhythmic events (LAEs) in Brugada syndrome (BrS) remains controversial. We investigated whether the functional curation of SCN5A variations improves prognostic predictability. METHODS AND RESULTS: Using a heterologous expression system and whole-cell patch clamping, we functionally characterized 22 variants of unknown significance (VUSs) among 55 SCN5A mutations previously curated using in silico prediction algorithms in the Japanese BrS registry (n = 415). According to the loss-of-function (LOF) properties, SCN5A mutation carriers (n = 60) were divided into two groups: LOF-SCN5A mutations and non-LOF SCN5A variations. Functionally proven LOF-SCN5A mutation carriers (n = 45) showed significantly severer electrocardiographic conduction abnormalities and worse prognosis associated with earlier manifestations of LAEs (7.9%/year) than in silico algorithm-predicted SCN5A carriers (5.1%/year) or all BrS probands (2.5%/year). Notably, non-LOF SCN5A variation carriers (n = 15) exhibited no LAEs during the follow-up period. Multivariate analysis demonstrated that only LOF-SCN5A mutations and a history of aborted cardiac arrest were significant predictors of LAEs. Gene-based association studies using whole-exome sequencing data on another independent SCN5A mutation-negative BrS cohort (n = 288) showed no significant enrichment of rare variants in 16 985 genes including 22 non-SCN5A BrS-associated genes as compared with controls (n = 372). Furthermore, rare variations of non-SCN5A BrS-associated genes did not affect LAE-free survival curves. CONCLUSION: In vitro functional validation is key to classifying the pathogenicity of SCN5A VUSs and for risk stratification of genetic predictors of LAEs. Functionally proven LOF-SCN5A mutations are genetic burdens of sudden death in BrS, but evidence for other BrS-associated genes is elusive.

Rare Coding Variants in ANGPTL6 Are Associated with Familial Forms of Intracranial Aneurysm.

Intracranial aneurysms (IAs) are acquired cerebrovascular abnormalities characterized by localized dilation and wall thinning in intracranial arteries, possibly leading to subarachnoid hemorrhage and severe outcome in case of rupture. Here, we identified one rare nonsense variant (c.1378A>T) in the last exon of ANGPTL6 (Angiopoietin-Like 6)-which encodes a circulating pro-angiogenic factor mainly secreted from the liver-shared by the four tested affected members of a large pedigree with multiple IA-affected case subjects. We showed a 50% reduction of ANGPTL6 serum concentration in individuals heterozygous for the c.1378A>T allele (p.Lys460Ter) compared to relatives homozygous for the normal allele, probably due to the non-secretion of the truncated protein produced by the c.1378A>T transcripts. Sequencing ANGPTL6 in a series of 94 additional index case subjects with familial IA identified three other rare coding variants in five case subjects. Overall, we detected a significant enrichment (p = 0.023) in rare coding variants within this gene among the 95 index case subjects with familial IA, compared to a reference population of 404 individuals with French ancestry. Among the 6 recruited families, 12 out of 13 (92%) individuals carrying IA also carry such variants in ANGPTL6, versus 15 out of 41 (37%) unaffected ones. We observed a higher rate of individuals with a history of high blood pressure among affected versus healthy individuals carrying ANGPTL6 variants, suggesting that ANGPTL6 could trigger cerebrovascular lesions when combined with other risk factors such as hypertension. Altogether, our results indicate that rare coding variants in ANGPTL6 are causally related to familial forms of IA.

Mutations in DCHS1 cause mitral valve prolapse.

Mitral valve prolapse (MVP) is a common cardiac valve disease that affects nearly 1 in 40 individuals. It can manifest as mitral regurgitation and is the leading indication for mitral valve surgery. Despite a clear heritable component, the genetic aetiology leading to non-syndromic MVP has remained elusive. Four affected individuals from a large multigenerational family segregating non-syndromic MVP underwent capture sequencing of the linked interval on chromosome 11. We report a missense mutation in the DCHS1 gene, the human homologue of the Drosophila cell polarity gene dachsous (ds), that segregates with MVP in the family. Morpholino knockdown of the zebrafish homologue dachsous1b resulted in a cardiac atrioventricular canal defect that could be rescued by wild-type human DCHS1, but not by DCHS1 messenger RNA with the familial mutation. Further genetic studies identified two additional families in which a second deleterious DCHS1 mutation segregates with MVP. Both DCHS1 mutations reduce protein stability as demonstrated in zebrafish, cultured cells and, notably, in mitral valve interstitial cells (MVICs) obtained during mitral valve repair surgery of a proband. Dchs1(+/-) mice had prolapse of thickened mitral leaflets, which could be traced back to developmental errors in valve morphogenesis. DCHS1 deficiency in MVP patient MVICs, as well as in Dchs1(+/-) mouse MVICs, result in altered migration and cellular patterning, supporting these processes as aetiological underpinnings for the disease. Understanding the role of DCHS1 in mitral valve development and MVP pathogenesis holds potential for therapeutic insights for this very common disease.

Shared genetic risk between eating disorder- and substance-use-related phenotypes: Evidence from genome-wide association studies.

Eating disorders and substance use disorders frequently co-occur. Twin studies reveal shared genetic variance between liabilities to eating disorders and substance use, with the strongest associations between symptoms of bulimia nervosa and problem alcohol use (genetic correlation [rg ], twin-based = 0.23-0.53). We estimated the genetic correlation between eating disorder and substance use and disorder phenotypes using data from genome-wide association studies (GWAS). Four eating disorder phenotypes (anorexia nervosa [AN], AN with binge eating, AN without binge eating, and a bulimia nervosa factor score), and eight substance-use-related phenotypes (drinks per week, alcohol use disorder [AUD], smoking initiation, current smoking, cigarettes per day, nicotine dependence, cannabis initiation, and cannabis use disorder) from eight studies were included. Significant genetic correlations were adjusted for variants associated with major depressive disorder and schizophrenia. Total study sample sizes per phenotype ranged from ~2400 to ~537 000 individuals. We used linkage disequilibrium score regression to calculate single nucleotide polymorphism-based genetic correlations between eating disorder- and substance-use-related phenotypes. Significant positive genetic associations emerged between AUD and AN (rg = 0.18; false discovery rate q = 0.0006), cannabis initiation and AN (rg = 0.23; q < 0.0001), and cannabis initiation and AN with binge eating (rg = 0.27; q = 0.0016). Conversely, significant negative genetic correlations were observed between three nondiagnostic smoking phenotypes (smoking initiation, current smoking, and cigarettes per day) and AN without binge eating (rgs = -0.19 to -0.23; qs < 0.04). The genetic correlation between AUD and AN was no longer significant after co-varying for major depressive disorder loci. The patterns of association between eating disorder- and substance-use-related phenotypes highlights the potentially complex and substance-specific relationships among these behaviors.

RRAD mutation causes electrical and cytoskeletal defects in cardiomyocytes derived from a familial case of Brugada syndrome.

AIMS: The Brugada syndrome (BrS) is an inherited cardiac disorder predisposing to ventricular arrhythmias. Despite considerable efforts, its genetic basis and cellular mechanisms remain largely unknown. The objective of this study was to identify a new susceptibility gene for BrS through familial investigation. METHODS AND RESULTS: Whole-exome sequencing performed in a three-generation pedigree with five affected members allowed the identification of one rare non-synonymous substitution (p.R211H) in RRAD, the gene encoding the RAD GTPase, carried by all affected members of the family. Three additional rare missense variants were found in 3/186 unrelated index cases. We detected higher levels of RRAD transcripts in subepicardium than in subendocardium in human heart, and in the right ventricle outflow tract compared to the other cardiac compartments in mice. The p.R211H variant was then subjected to electrophysiological and structural investigations in human cardiomyocytes derived from induced pluripotent stem cells (iPSC-CMs). Cardiomyocytes derived from induced pluripotent stem cells from two affected family members exhibited reduced action potential upstroke velocity, prolonged action potentials and increased incidence of early afterdepolarizations, with decreased Na+ peak current amplitude and increased Na+ persistent current amplitude, as well as abnormal distribution of actin and less focal adhesions, compared with intra-familial control iPSC-CMs Insertion of p.R211H-RRAD variant in control iPSCs by genome editing confirmed these results. In addition, iPSC-CMs from affected patients exhibited a decreased L-type Ca2+ current amplitude. CONCLUSION: This study identified a potential new BrS-susceptibility gene, RRAD. Cardiomyocytes derived from induced pluripotent stem cells expressing RRAD variant recapitulated single-cell electrophysiological features of BrS, including altered Na+ current, as well as cytoskeleton disturbances.

SMMB: a stochastic Markov blanket framework strategy for epistasis detection in GWAS.

Motivation: Large scale genome-wide association studies (GWAS) are tools of choice for discovering associations between genotypes and phenotypes. To date, many studies rely on univariate statistical tests for association between the phenotype and each assayed single nucleotide polymorphism (SNP). However, interaction between SNPs, namely epistasis, must be considered when tackling the complexity of underlying biological mechanisms. Epistasis analysis at large scale entails a prohibitive computational burden when addressing the detection of more than two interacting SNPs. In this paper, we introduce a stochastic causal graph-based method, SMMB, to analyze epistatic patterns in GWAS data. Results: We present Stochastic Multiple Markov Blanket algorithm (SMMB), which combines both ensemble stochastic strategy inspired from random forests and Bayesian Markov blanket-based methods. We compared SMMB with three other recent algorithms using both simulated and real datasets. Our method outperforms the other compared methods for a majority of simulated cases of 2-way and 3-way epistasis patterns (especially in scenarii where minor allele frequencies of causal SNPs are low). Our approach performs similarly as two other compared methods for large real datasets, in terms of power, and runs faster. Availability and implementation: Parallel version available on https://ls2n.fr/listelogicielsequipe/DUKe/128/. Supplementary information: Supplementary data are available at Bioinformatics online.

New insights into mitral valve dystrophy: a Filamin-A genotype-phenotype and outcome study.

Aims: Filamin-A (FLNA) was identified as the first gene of non-syndromic mitral valve dystrophy (FLNA-MVD). We aimed to assess the phenotype of FLNA-MVD and its impact on prognosis. Methods and results: We investigated the disease in 246 subjects (72 mutated) from four FLNA-MVD families harbouring three different FLNA mutations. Phenotype was characterized by a comprehensive echocardiography focusing on mitral valve apparatus in comparison with control relatives. In this X-linked disease valves lesions were severe in men and moderate in women. Most men had classical features of mitral valve prolapse (MVP), but without chordal rupture. By contrast to regular MVP, mitral leaflet motion was clearly restricted in diastole and papillary muscles position was closer to mitral annulus. Valvular abnormalities were similar in the four families, in adults and young patients from early childhood suggestive of a developmental disease. In addition, mitral valve lesions worsened over time as encountered in degenerative conditions. Polyvalvular involvement was frequent in males and non-diagnostic forms frequent in females. Overall survival was moderately impaired in men (P = 0.011). Cardiac surgery rate (mainly valvular) was increased (33.3 ± 9.8 vs. 5.0 ± 4.9%, P < 0.0001; hazard ratio 10.5 [95% confidence interval: 2.9-37.9]) owing mainly to a lifetime increased risk in men (76.8 ± 14.1 vs. 9.1 ± 8.7%, P < 0.0001). Conclusion: FLNA-MVD is a developmental and degenerative disease with complex phenotypic expression which can influence patient management. FLNA-MVD has unique features with both MVP and paradoxical restricted motion in diastole, sub-valvular mitral apparatus impairment and polyvalvular lesions in males. FLNA-MVD conveys a substantial lifetime risk of valve surgery in men.

Seventy-five genetic loci influencing the human red blood cell.

Anaemia is a chief determinant of global ill health, contributing to cognitive impairment, growth retardation and impaired physical capacity. To understand further the genetic factors influencing red blood cells, we carried out a genome-wide association study of haemoglobin concentration and related parameters in up to 135,367 individuals. Here we identify 75 independent genetic loci associated with one or more red blood cell phenotypes at P < 10(-8), which together explain 4-9% of the phenotypic variance per trait. Using expression quantitative trait loci and bioinformatic strategies, we identify 121 candidate genes enriched in functions relevant to red blood cell biology. The candidate genes are expressed preferentially in red blood cell precursors, and 43 have haematopoietic phenotypes in Mus musculus or Drosophila melanogaster. Through open-chromatin and coding-variant analyses we identify potential causal genetic variants at 41 loci. Our findings provide extensive new insights into genetic mechanisms and biological pathways controlling red blood cell formation and function.

Testing the burden of rare variation in arrhythmia-susceptibility genes provides new insights into molecular diagnosis for Brugada syndrome.

The Brugada syndrome (BrS) is a rare heritable cardiac arrhythmia disorder associated with ventricular fibrillation and sudden cardiac death. Mutations in the SCN5A gene have been causally related to BrS in 20-30% of cases. Twenty other genes have been described as involved in BrS, but their overall contribution to disease prevalence is still unclear. This study aims to estimate the burden of rare coding variation in arrhythmia-susceptibility genes among a large group of patients with BrS. We have developed a custom kit to capture and sequence the coding regions of 45 previously reported arrhythmia-susceptibility genes and applied this kit to 167 index cases presenting with a Brugada pattern on the electrocardiogram as well as 167 individuals aged over 65-year old and showing no history of cardiac arrhythmia. By applying burden tests, a significant enrichment in rare coding variation (with a minor allele frequency below 0.1%) was observed only for SCN5A, with rare coding variants carried by 20.4% of cases with BrS versus 2.4% of control individuals (P = 1.4 × 10(-7)). No significant enrichment was observed for any other arrhythmia-susceptibility gene, including SCN10A and CACNA1C. These results indicate that, except for SCN5A, rare coding variation in previously reported arrhythmia-susceptibility genes do not contribute significantly to the occurrence of BrS in a population with European ancestry. Extreme caution should thus be taken when interpreting genetic variation in molecular diagnostic setting, since rare coding variants were observed in a similar extent among cases versus controls, for most previously reported BrS-susceptibility genes.

Utilizing the Jaccard index to reveal population stratification in sequencing data: a simulation study and an application to the 1000 Genomes Project.

MOTIVATION: Population stratification is one of the major sources of confounding in genetic association studies, potentially causing false-positive and false-negative results. Here, we present a novel approach for the identification of population substructure in high-density genotyping data/next generation sequencing data. The approach exploits the co-appearances of rare genetic variants in individuals. The method can be applied to all available genetic loci and is computationally fast. Using sequencing data from the 1000 Genomes Project, the features of the approach are illustrated and compared to existing methodology (i.e. EIGENSTRAT). We examine the effects of different cutoffs for the minor allele frequency on the performance of the approach. We find that our approach works particularly well for genetic loci with very small minor allele frequencies. The results suggest that the inclusion of rare-variant data/sequencing data in our approach provides a much higher resolution picture of population substructure than it can be obtained with existing methodology. Furthermore, in simulation studies, we find scenarios where our method was able to control the type 1 error more precisely and showed higher power. CONTACT: dmitry.prokopenko@uni-bonn.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

De Novo Truncating Mutations in the Kinetochore-Microtubules Attachment Gene CHAMP1 Cause Syndromic Intellectual Disability.

A rare syndromic form of intellectual disability with impaired speech was recently found associated with mutations in CHAMP1 (chromosome alignment-maintaining phosphoprotein 1), the protein product of which is directly involved in microtubule-kinetochore attachment. Through whole-exome sequencing in six unrelated nonconsanguineous families having a sporadic case of intellectual disability, we identified six novel de novo truncating mutations in CHAMP1: c.1880C>G p.(Ser627*), c.1489C>T; p.(Arg497*), c.1876_1877delAG; p.(Ser626Leufs*4), c.1043G>A; p.(Trp348*), c.1002G>A; p.(Trp334*), and c.958_959delCC; p.(Pro320*). Our clinical observations confirm the phenotypic homogeneity of the syndrome, which represents therefore a distinct clinical entity. Besides, our functional studies show that CHAMP1 protein variants are delocalized from chromatin and are unable to bind to two of its direct partners, POGZ and HP1. These data suggest a pathogenic mechanism of the CHAMP1-associated intellectual disability syndrome mediated by direct interacting partners of CHAMP1, several of which are involved in chromo/kinetochore-related disorders.

Variation in FTO contributes to childhood obesity and severe adult obesity.

We identified a set of SNPs in the first intron of the FTO (fat mass and obesity associated) gene on chromosome 16q12.2 that is consistently strongly associated with early-onset and severe obesity in both adults and children of European ancestry with an experiment-wise P value of 1.67 x 10(-26) in 2,900 affected individuals and 5,100 controls. The at-risk haplotype yields a proportion of attributable risk of 22% for common obesity. We conclude that FTO contributes to human obesity and hence may be a target for subsequent functional analyses.

Novel association approach for variable number tandem repeats (VNTRs) identifies DOCK5 as a susceptibility gene for severe obesity.

Variable number tandem repeats (VNTRs) constitute a relatively under-examined class of genomic variants in the context of complex disease because of their sequence complexity and the challenges in assaying them. Recent large-scale genome-wide copy number variant mapping and association efforts have highlighted the need for improved methodology for association studies using these complex polymorphisms. Here we describe the in-depth investigation of a complex region on chromosome 8p21.2 encompassing the dedicator of cytokinesis 5 (DOCK5) gene. The region includes two VNTRs of complex sequence composition which flank a common 3975 bp deletion, all three of which were genotyped by polymerase chain reaction and fragment analysis in a total of 2744 subjects. We have developed a novel VNTR association method named VNTRtest, suitable for association analysis of multi-allelic loci with binary and quantitative outcomes, and have used this approach to show significant association of the DOCK5 VNTRs with childhood and adult severe obesity (P(empirical)= 8.9 × 10(-8) and P= 3.1 × 10(-3), respectively) which we estimate explains ~0.8% of the phenotypic variance. We also identified an independent association between the 3975 base pair (bp) deletion and obesity, explaining a further 0.46% of the variance (P(combined)= 1.6 × 10(-3)). Evidence for association between DOCK5 transcript levels and the 3975 bp deletion (P= 0.027) and both VNTRs (P(empirical)= 0.015) was also identified in adipose tissue from a Swedish family sample, providing support for a functional effect of the DOCK5 deletion and VNTRs. These findings highlight the potential role of DOCK5 in human obesity and illustrate a novel approach for analysis of the contribution of VNTRs to disease susceptibility through association studies.