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.

How local reference panels improve imputation in French populations

Imputation servers offer the exclusive possibility to harness the largest public reference panels which have been shown to deliver very high precision in the imputation of European genomes. Many studies have nonetheless stressed the importance of 'study specific panels' (SSPs) as an alternative and have shown the benefits of combining public reference panels with SSPs. But such combined approaches are not attainable when using external imputation servers. To investigate how to confront this challenge, we imputed 550 French individuals using either the University of Michigan imputation server with the Haplotype Reference Consortium (HRC) panel or an in-house SSP of 850 whole-genome sequenced French individuals. With approximate geo-localization of both our target and SSP individuals we are able to pinpoint different scenarios where SSP-based imputation would be preferred over server-based imputation or vice-versa. This is achieved by showing to a high degree of resolution the importance of the proximity of the reference panel to target individuals; with a focus on the clear added value of SSPs for estimating haplotype phase and for the imputation of rare variants (minor allele-frequency below 0.01). Such benefits were most evident for individuals from the same geographical regions in France as the SSP individuals. Overall, only 42.3% of all 125,442 variants evaluated were better imputed with an SSP from France compared to an external reference panel, however this rises to 58.1% for individuals from geographic regions well covered by the SSP. By investigating haplotype sharing and population fine-structure in France, we show the importance of including SSP haplotypes for imputation but also that they should ideally be combined with large public panels. In the absence of the unattainable results from a combined panel of the HRC and our French SSP, we put forward a pragmatic solution where server-based and SSP-based imputation outcomes can be combined based on comparing posterior genotype probabilities. We show that such an approach can give a level of imputation accuracy in excess of what could be achieved with either strategy alone. The results presented provide detailed insights into the accuracy of imputation that should be expected from different strategies for European populations.

Left Ventricular Abnormal Substrate in Brugada Syndrome

BACKGROUND

Slow-conductive structural abnormalities located in the epicardium of the right ventricle (RV) underlie Brugada syndrome (BrS). The extent of such substrate in the left ventricle (LV) has not been investigated.

OBJECTIVES

This study sought to characterize the extent of epicardial substrate abnormalities in BrS.

METHODS

We evaluated 22 consecutive patients (mean age 46 ± 11 years, 21 male) referred for recurrent ventricular arrhythmias (mean 10 ± 13 episodes) in the setting of BrS. The patients underwent clinical investigations and wide genetic screening to identify SCN5A mutations and common risk variants. High-density biventricular epicardial mapping was performed to detect prolonged (>70 ms) fragmented electrograms, indicating abnormal substrate area.

RESULTS

All patients presented with abnormal substrate in the epicardial anterior RV (27 ± 11 cm2). Abnormal substrate was also identified on the LV epicardium in 10 patients (45%), 9 at baseline and 1 after ajmaline infusion, covering 15 ± 11 cm2. Of these, 4 had severe LV fascicular blocks. Patients with LV substrate had a longer history of arrhythmia (11.4 ± 6.7 years vs 4.3 ± 4.3 years; P = 0.003), longer PR (217 ± 24 ms vs 171 ± 14 ms; P < 0.001) and HV (60 ± 12 ms vs 46 ± 5 ms; P = 0.005) intervals, and abnormal substrate also extending into the inferior RV (100% vs 33%; P = 0.001). SCN5A mutation was present in 70% of patients with LV substrate (vs 25%; P = 0.035). SCN5A BrS patients with recurrent ventricular arrhythmias present a higher polygenic risk score compared with a nonselected BrS population (median of differences: -0.86; 95% CI: -1.48 to -0.27; P = 0.02).

CONCLUSIONS

A subset of patients with BrS present an abnormal substrate extending onto the LV epicardium and inferior RV that is associated with SCN5A mutations and multigenic variants.

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.

On the genes, genealogies, and geographies of Quebec

Population genetic models only provide coarse representations of real-world ancestry. We used a pedigree compiled from 4 million parish records and genotype data from 2276 French and 20,451 French Canadian individuals to finely model and trace French Canadian ancestry through space and time. The loss of ancestral French population structure and the appearance of spatial and regional structure highlights a wide range of population expansion models. Geographic features shaped migrations, and we find enrichments for migration, genetic, and genealogical relatedness patterns within river networks across regions of Quebec. Finally, we provide a freely accessible simulated whole-genome sequence dataset with spatiotemporal metadata for 1,426,749 individuals reflecting intricate French Canadian population structure. Such realistic population-scale simulations provide opportunities to investigate population genetics at an unprecedented resolution.

Heritable defects in telomere and mitotic function selectively predispose to sarcomas

Cancer genetics has to date focused on epithelial malignancies, identifying multiple histotype-specific pathways underlying cancer susceptibility. Sarcomas are rare malignancies predominantly derived from embryonic mesoderm. To identify pathways specific to mesenchymal cancers, we performed whole-genome germline sequencing on 1644 sporadic cases and 3205 matched healthy elderly controls. Using an extreme phenotype design, a combined rare-variant burden and ontologic analysis identified two sarcoma-specific pathways involved in mitotic and telomere functions. Variants in centrosome genes are linked to malignant peripheral nerve sheath and gastrointestinal stromal tumors, whereas heritable defects in the shelterin complex link susceptibility to sarcoma, melanoma, and thyroid cancers. These studies indicate a specific role for heritable defects in mitotic and telomere biology in risk of sarcomas.

Genetics and pathophysiology of mitral valve prolapse

Mitral valve prolapse (MVP) is a common condition affecting 2-3% of the general population, and the most complex form of valve pathology, with a complication rate up to 10-15% per year in advanced stages. Complications include mitral regurgitation which can lead to heart failure and atrial fibrillation, but also life-threatening ventricular arrhythmia and cardiovascular death. Sudden death has been recently brought to the forefront of MVP disease, increasing the complexity of management and suggesting that MVP condition is not properly understood. MVP can occur as part of syndromic conditions such as Marfan syndrome, but the most common form is non-syndromic, isolated or familial. Although a specific X-linked form of MVP was initially identified, autosomal dominant inheritance appears to be the primary mode of transmission. MVP can be stratified into myxomatous degeneration (Barlow), fibroelastic deficiency, and Filamin A-related MVP. While FED is still considered a degenerative disease associated with aging, myxomatous MVP and FlnA-MVP are recognized as familial pathologies. Deciphering genetic defects associated to MVP is still a work in progress; although FLNA, DCHS1, and DZIP1 have been identified as causative genes in myxomatous forms of MVP thanks to familial approaches, they explain only a small proportion of MVP. In addition, genome-wide association studies have revealed the important role of common variants in the development of MVP, in agreement with the high prevalence of this condition in the population. Furthermore, a potential genetic link between MVP and ventricular arrhythmia or a specific type of cardiomyopathy is considered. Animal models that allow to advance in the genetic and pathophysiological knowledge of MVP, and in particular those that can be easily manipulated to express a genetic defect identified in humans are detailed. Corroborated by genetic data and animal models, the main pathophysiological pathways of MVP are briefly addressed. Finally, genetic counseling is considered in the context of MVP.

Common Genetic Variation and Age of Onset of Anorexia Nervosa

Background

Genetics and biology may influence the age of onset of anorexia nervosa (AN). The aims of this study were to determine whether common genetic variation contributes to age of onset of AN and to investigate the genetic associations between age of onset of AN and age at menarche.

Methods

A secondary analysis of the Psychiatric Genomics Consortium genome-wide association study (GWAS) of AN was performed, which included 9335 cases and 31,981 screened controls, all from European ancestries. We conducted GWASs of age of onset, early-onset AN (<13 years), and typical-onset AN, and genetic correlation, genetic risk score, and Mendelian randomization analyses.

Results

Two loci were genome-wide significant in the typical-onset AN GWAS. Heritability estimates (single nucleotide polymorphism-h 2) were 0.01-0.04 for age of onset, 0.16-0.25 for early-onset AN, and 0.17-0.25 for typical-onset AN. Early- and typical-onset AN showed distinct genetic correlation patterns with putative risk factors for AN. Specifically, early-onset AN was significantly genetically correlated with younger age at menarche, and typical-onset AN was significantly negatively genetically correlated with anthropometric traits. Genetic risk scores for age of onset and early-onset AN estimated from independent GWASs significantly predicted age of onset. Mendelian randomization analysis suggested a causal link between younger age at menarche and early-onset AN.

Conclusions

Our results provide evidence consistent with a common variant genetic basis for age of onset and implicate biological pathways regulating menarche and reproduction.

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 in desmosomal genes have been causally related to ACR in 40% of cases. Other genes encoding non desmosomal 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. Pathogenic variants (PV) and variants of uncertain significance (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, DSP, DSC2 and DSG2. 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. This article is protected by copyright. All rights reserved.

Familial Recurrence Patterns in Congenitally Corrected Transposition of the Great Arteries: An International Study

BACKGROUND

Congenitally corrected transposition of the great arteries (ccTGA) is a rare disease of unknown cause. We aimed to better understand familial recurrence patterns.

METHODS

An international, multicentre, retrospective cohort study was conducted in 29 tertiary hospitals in 6 countries between 1990 and 2018, entailing investigation of 1043 unrelated ccTGA probands.

RESULTS

Laterality defects and atrioventricular block at diagnosis were observed in 29.9% and 9.3%, respectively. ccTGA was associated with primary ciliary dyskinesia in 11 patients. Parental consanguinity was noted in 3.4% cases. A congenital heart defect was diagnosed in 81 relatives from 69 families, 58% of them being first-degree relatives, including 28 siblings. The most prevalent defects in relatives were dextro-transposition of the great arteries (28.4%), laterality defects (13.6%), and ccTGA (11.1%); 36 new familial clusters were described, including 8 pedigrees with concordant familial aggregation of ccTGA, 19 pedigrees with familial co-segregation of ccTGA and dextro-transposition of the great arteries, and 9 familial co-segregation of ccTGA and laterality defects. In one family co-segregation of ccTGA, dextro-transposition of the great arteries and heterotaxy syndrome in 3 distinct relatives was found. In another family, twins both displayed ccTGA and primary ciliary dyskinesia.

CONCLUSIONS

ccTGA is not always a sporadic congenital heart defect. Familial clusters as well as evidence of an association between ccTGA, dextro-transposition of the great arteries, laterality defects and in some cases primary ciliary dyskinesia, strongly suggest a common pathogenetic pathway involving laterality genes in the pathophysiology of ccTGA.

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.

Genome-wide association analyses identify new Brugada syndrome risk loci and highlight a new mechanism of sodium channel regulation in disease susceptibility

Brugada syndrome (BrS) is a cardiac arrhythmia disorder associated with sudden death in young adults. With the exception of SCN5A, encoding the cardiac sodium channel Na_V1.5, susceptibility genes remain largely unknown. Here we performed a genome-wide association meta-analysis comprising 2,820 unrelated cases with BrS and 10,001 controls, and identified 21 association signals at 12 loci (10 new). Single nucleotide polymorphism (SNP)-heritability estimates indicate a strong polygenic influence. Polygenic risk score analyses based on the 21 susceptibility variants demonstrate varying cumulative contribution of common risk alleles among different patient subgroups, as well as genetic associations with cardiac electrical traits and disorders in the general population. The predominance of cardiac transcription factor loci indicates that transcriptional regulation is a key feature of BrS pathogenesis. Furthermore, functional studies conducted on MAPRE2, encoding the microtubule plus-end binding protein EB2, point to microtubule-related trafficking effects on Na_V1.5 expression as a new underlying molecular mechanism. Taken together, these findings broaden our understanding of the genetic architecture of BrS and provide new insights into its molecular underpinnings.

Computational estimates of annular diameter reveal genetic determinants of mitral valve function and disease

The fibrous annulus of the mitral valve plays an important role in valvular function and cardiac physiology, while normal variation in the size of cardiovascular anatomy may share a genetic link with common and rare disease. We derived automated estimates of mitral valve annular diameter in the 4-chamber view from 32,220 MRI images from the UK Biobank at ventricular systole and diastole as the basis for GWAS. Mitral annular dimensions corresponded to previously described anatomical norms, and GWAS inclusive of 4 population strata identified 10 loci, including possibly novel loci (GOSR2, ERBB4, MCTP2, MCPH1) and genes related to cardiac contractility (BAG3, TTN, RBFOX1). ATAC-Seq of primary mitral valve tissue localized multiple variants to regions of open chromatin in biologically relevant cell types and rs17608766 to an algorithmically predicted enhancer element in GOSR2. We observed strong genetic correlation with measures of contractility and mitral valve disease and clinical correlations with heart failure, cerebrovascular disease, and ventricular arrhythmias. Polygenic scoring of mitral valve annular diameter in systole was predictive of risk mitral valve prolapse across 4 cohorts. In summary, genetic and clinical studies of mitral valve annular diameter revealed genetic determinants of mitral valve biology, while highlighting clinical associations. Polygenic determinants of mitral valve annular diameter may represent an independent risk factor for mitral prolapse. Overall, computationally estimated phenotypes derived at scale from medical imaging represent an important substrate for genetic discovery and clinical risk prediction.

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.

Genome-Wide Association Meta-Analysis Supports Genes Involved in Valve and Cardiac Development to Associate With Mitral Valve Prolapse

BACKGROUND

Mitral valve prolapse (MVP) is a common cardiac valve disease, which affects 1 in 40 in the general population. Previous genome-wide association study has identified 6 risk loci for MVP. But these loci explained only partially the genetic risk for MVP. We aim to identify additional risk loci for MVP by adding data set from the UK Biobank.

METHODS

We also incorporated 434 MVP cases and 4527 controls from the UK Biobank for discovery analyses. Genetic association was conducted using SNPTEST and meta-analyses using METAL. We used Functional Mapping and Annotation of Genome-Wide Association Studies for post-genome-wide association study annotations and Multi-marker Analysis of GenoMic Annotation for gene-based and gene-set analyses.

RESULTS

We found Trans-Omics for Precision Medicine imputation to perform better in terms of accuracy in the lower ranges of minor allele frequency below 0.1. Our updated meta-analysis included UK Biobank study for ≈8 million common single-nucleotide polymorphisms (minor allele frequency >0.01) and replicated the association on Chr2 as the top association signal near TNS1. We identified an additional risk locus on Chr1 (SYT2) and 2 suggestive risk loci on chr8 (MSRA) and chr19 (FBXO46), all driven by common variants. Gene-based association using Multi-marker Analysis of GenoMic Annotation revealed 6 risk genes for MVP with pronounced expression levels in cardiovascular tissues, especially the heart and globally part of enriched GO terms related to cardiac development.

CONCLUSIONS

We report an updated meta-analysis genome-wide association study for MVP using dense imputation coverage and an improved case-control sample. We describe several loci and genes with MVP spanning biological mechanisms highly relevant to MVP, especially during valve and heart 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.

Genome-wide association study of intracranial aneurysms identifies 17 risk loci and genetic overlap with clinical risk factors

Rupture of an intracranial aneurysm leads to subarachnoid hemorrhage, a severe type of stroke. To discover new risk loci and the genetic architecture of intracranial aneurysms, we performed a cross-ancestry, genome-wide association study in 10,754 cases and 306,882 controls of European and East Asian ancestry. We discovered 17 risk loci, 11 of which are new. We reveal a polygenic architecture and explain over half of the disease heritability. We show a high genetic correlation between ruptured and unruptured intracranial aneurysms. We also find a suggestive role for endothelial cells by using gene mapping and heritability enrichment. Drug-target enrichment shows pleiotropy between intracranial aneurysms and antiepileptic and sex hormone drugs, providing insights into intracranial aneurysm pathophysiology. Finally, genetic risks for smoking and high blood pressure, the two main clinical risk factors, play important roles in intracranial aneurysm risk, and drive most of the genetic correlation between intracranial aneurysms and other cerebrovascular traits.

Genome-wide association study identifies 18 new susceptibility variants loci associated with Brugada Syndrome

 

The Brugada Syndrome (BrS) is characterized by ST-segment elevation in the right precordial leads and is associated with an increased risk of sudden cardiac death. The disorder was initially described as a monogenic primary cardiac electrical disease. However, mutations in SCN5A, encoding the cardiac sodium channel (NaV1.5), which is the major gene associated with the disorder are found in only around 20% of cases and are associated with low penetrance. Furthermore many cases did not display familial aggregation. Based on a previous GWAS conducted on 312 BrS patients and the discovery of the unexpected strong effect of 3 common variants, we proposed that the BrS may comprise a more complex inheritance model. We conducted a genome-wide association study on 2820 individuals with BrS and 10001 ancestry-matched controls to uncover additional genetic loci that modulate susceptibility to BrS, to characterize further the BrS genetic architecture and to uncover new molecular mechanisms.

We identified 21 susceptibility variants that passed the genome-wide statistical significance threshold (P&lt;5.10–8), of which 18 were novel. Eight were located at the SCN5A-SCN10A locus, illustrating the central role of NaV1.5 in the disease. Interestingly, 9 occur in the vicinity of genes known to play a crucial role in cardiac development (HEY2, TBX20, GATA4, ZFPM2, WT1, TBX5, IRX3, IRX5) and / or control cardiac ion channel expression. Of note, 2 others signals occurred in the vicinity of microtubule / cytoskeleton associated proteins (MAPRE2 and MYO18B). Through studies in zebrafish and in human iPSC-derived cardiomyocytes, we demonstrate a role of MAPRE2 on NaV1.5 function.

We identified 18 new susceptibility variants associated with BrS and uncovered a new pathophysiological molecular mechanism underlying BrS susceptibility. We provided further support for a complex genetic architecture underlying susceptibility for the disorder.

Funding Acknowledgement

Type of funding source: Public Institution(s). Main funding source(s): H2020 - Marie Sklodowska Curie IF grant, Rising star grant from the Pays de la Loire regional council

GENECHOC study: genetic markers of arrhythmic risk in heart failure

Background

Ventricular arrhythmic events are responsible for 50% of death in heart failure but no reliable predictive marker is known to discriminate patients at risk of fatal arrhythmia. Interestingly, familial predisposition has been reported suggesting a role of genetic factors.

Purpose

Identify genetic markers increasing the arrhythmic risk in heart failure population.

Method

We prospectively included heart failure patients with left ventricular ejection fraction (LVEF) under 35% and a cardioverter defibrillator in primary prevention in 22 French centres between 2009 and 2017. Patients were followed for 72 months and divided into two groups: cases with an arrhythmic event during follow-up and controls. A Genome Wide Association Study (GWAS) was done. Single Nucleotide Polymorphisms (SNPs) genotyping was performed on Affymetrix Axiom Precision Medicine Research Array plates. To complement the directly genotyped SNPs we performed large-scale imputation based on the Haplotype Reference Consortium European ancestry panel leading to a dataset of 7,5 million of SNPs.

Results

332 cases and 567 controls were included (86% men, mean age at implantation 52±11 years). 78% of patients had ischaemic cardiopathy, 20% had dilated cardiomyopathy. Mean LVEF was 27±5%. No statistical difference was found between cases and controls on clinical parameters, biological results, electrocardiographic measures. No locus shows genome-wide significant association (p&lt;5.10–8) on the GWAS analysis. However, 16 signals with a p-value between 5.10–8 and 5.10–5 were investigated. eQTL and chromatin conformation point to 35 genes with cardiac expression previously associated with heart failure, cardiomyopathies, cardiogenesis, arrhythmias and inflammation. Variants identified point to regulatory regions of the genome and may then propose a molecular mechanism predisposing patients to arrhythmias.

Conclusion

No locus raises genome-wide significance, but several signals with a nominal p-value point to relevant genes and pathways. Replication of the GWAS is ongoing on a cohort of 156 new patients with a less severe cardiopathy implanted with a cardioverter defibrillator in secondary prevention.

Funding Acknowledgement

Type of funding source: Public hospital(s). Main funding source(s): Nantes University Hospital

Heritability of aortic valve stenosis and evidence of bicuspid enrichment in families with tricuspid aortic valve stenosis

Background

Calcific aortic valve stenosis (CAVS) is the most common valvular heart disease. Although a familial component has been described, the heritability remains unknown. This study aimed to determine the heritability of CAVS and to compare the prevalence of bicuspid aortic valve (BAV) relatives among families with tricuspid and bicuspid aortic valve phenotype and to a control cohort.

Methods

Probands were recruited following aortic valve replacement in three French university hospitals. All relatives who agreed to participate were enrolled in the study. Families were defined as tricuspid aortic valve (TAV) or BAV families based on the aortic valve phenotype observed in the proband. As control group, we used probands diagnosed with mitral valve prolapse and their families. All recruited probands and relatives underwent comprehensive transthoracic echocardiography to determine aortic valve morphology and CAVS status. Additive genetic heritability was performed with SOLAR genetic statistics package. Prevalence of BAV relatives in each sub-group and Prevalence Ratio were calculated.

Results

Within the 2371 relatives coming from 138 CAVS families, the heritability of CAVS was significant (h2=0.49, p&lt;0.001). Furthermore, TAV stenosis heritability was only significant in TAV families (h2=0.20, p=0.013). The prevalence of BAV was significantly higher in TAV families and BAV families than control families (5.96 [95% CI: 4.2–8.3%], [n=32 of 537] vs. 2.49 [95% CI: 1.3–4.7%], [n=9 of 361]; PR 2.64 [95% CI: 1.41–5.91]; p=0.005 and 10.20 [95% CI: 7.0–14.5%], [n=26 of 255] vs. 2.49 [95% CI: 1.3–4.7%], [n=9 of 361]; PR 4.63 [95% CI: 2.35–13.44]; p&lt;0.001; respectively).

Conclusions

The heritability of aortic valve stenosis in this population confirms the genetic effect in determining aortic stenosis for all aortic valve phenotype and for tricuspid aortic valve. The bicuspid enrichment in TAV families suggests common genetic predisposition of aortic valve stenosis in some families. Taking all together, these findings encourage echocardiography screening of relatives of affected individuals.

Funding Acknowledgement

Type of funding source: None

Genome-wide association meta-analysis in 652,134 participants identifies 9 novel susceptibility loci for aortic stenosis

Background

Aortic stenosis (AS) is the most common form of incident valvular heart disease. While valve replacement is effective, the absence of an approved medical therapy provides no alternatives to patients with contraindications or mild disease. An improved understanding of the genetics of AS could identify targets for pharmacological intervention.

Methods

An inverse variance-weighted, fixed effects meta-analysis of the association of 11,591,806 variants with AS was undertaken using data from 10 European cohorts totalling 652,134 participants (13,758 cases of AS). We queried publicly available datasets to characterize the functional consequences of genome-wide significant variants, conducted a phenome-wide association study to assess their association with other outcomes, and constructed polygenic risk scores to examine their association with AS. We also performed gene- and gene-set enrichment analyses, estimated genetic correlation with cardiovascular traits, and assessed whether five lipid or immunological biomarkers were causally associated with AS using Mendelian randomization.

Results

Eighteen independent variants at 16 loci attained genome-wide significance in the meta-analysis, including variants at all seven previously reported loci. Many of the significant variants were intronic or intergenic, and the phenome-wide association study revealed extensive pleiotropy with apolipoprotein B, C-reactive protein, and other cardiovascular and immunological traits. A weighted polygenic risk score composed of the 18 variants was strongly associated with AS (adjusted OR per SD, 1.38; 95% CI, 1.33 to 1.44; p=4.6×10–57), and improved the discriminatory ability for AS when added to a model that contained clinical risk factors (difference in the area under the curve p=2.0×10–11). Gene-based approaches indicated higher IL6R expression in the blood among AS cases compared to controls (p=3.1×10–6), and the association of LDLR with AS (p=2.3×10–10). Gene set analyses revealed that genes bound by the transcription factor TCF7 or micro-RNAs miR-21, miR-219, miR-491, and miR-19 were differentially expressed in the liver depending on AS status (p≤5.7×10–4), suggesting disease development may be mediated by tissue-specific transcriptional and post-transcriptional regulation. Mendelian randomization supported a causal association of five lipid and immunological biomarkers with AS, including low-density lipoprotein cholesterol (OR per mmol/L, 1.61; 95% CI, 1.48 to 1.75; p=1.3×10–30).

Conclusions

Evidence from large-scale genetic analyses indicate that lipid metabolism, inflammation, and calcification are key contributors to AS.

Funding Acknowledgement

Type of funding source: Public Institution(s). Main funding source(s): Canadian Institutes of Health Research, National Institutes of Health

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.

Lipoprotein-associated phospholipase A2 activity, genetics and calcific aortic valve stenosis in humans

BACKGROUND

Lipoprotein-associated phospholipase A2 (Lp-PLA2) activity has been shown to predict calcific aortic valve stenosis (CAVS) outcomes. Our objective was to test the association between plasma Lp-PLA2 activity and genetically elevated Lp-PLA2 mass/activity with CAVS in humans.

METHODS AND RESULTS

Lp-PLA2 activity was measured in 890 patients undergoing cardiac surgery, including 476 patients undergoing aortic valve replacement for CAVS and 414 control patients undergoing coronary artery bypass grafting. After multivariable adjustment, Lp-PLA2 activity was positively associated with the presence of CAVS (OR=1.21 (95% CI 1.04 to 1.41) per SD increment). We selected four single nucleotide polymorphisms (SNPs) at the PLA2G7 locus associated with either Lp-PLA2 mass or activity (rs7756935, rs1421368, rs1805017 and rs4498351). Genetic association studies were performed in eight cohorts: Quebec-CAVS (1009 cases/1017 controls), UK Biobank (1350 cases/349 043 controls), European Prospective Investigation into Cancer and Nutrition-Norfolk (504 cases/20 307 controls), Genetic Epidemiology Research on Aging (3469 cases/51 723 controls), Malmö Diet and Cancer Study (682 cases/5963 controls) and three French cohorts (3123 cases/6532 controls), totalling 10 137 CAVS cases and 434 585 controls. A fixed-effect meta-analysis using the inverse-variance weighted method revealed that none of the four SNPs was associated with CAVS (OR=0.99 (95% CI 0.96 to 1.02, p=0.55) for rs7756935, 0.97 (95% CI 0.93 to 1.01, p=0.11) for rs1421368, 1.00 (95% CI 1.00 to 1.01, p=0.29) for rs1805017, and 1.00 (95% CI 0.97 to 1.04, p=0.87) for rs4498351).

CONCLUSIONS

Higher Lp-PLA2 activity is significantly associated with the presence of CAVS and might represent a biomarker of CAVS in patients with heart disease. Results of our genetic association study suggest that Lp-PLA2 is however unlikely to represent a causal risk factor or therapeutic target for CAVS.

Genetic and In Vitro Inhibition of PCSK9 and Calcific Aortic Valve Stenosis

The authors investigated whether PCSK9 inhibition could represent a therapeutic strategy in calcific aortic valve stenosis (CAVS). A meta-analysis of 10 studies was performed to determine the impact of the PCSK9 R46L variant on CAVS, and the authors found that CAVS was less prevalent in carriers of this variant (odds ratio: 0.80 [95% confidence interval: 0.70 to 0.91]; p = 0.0011) compared with noncarriers. PCSK9 expression was higher in the aortic valves of patients CAVS compared with control patients. In human valve interstitials cells submitted to a pro-osteogenic medium, PCSK9 levels increased and a PCSK9 neutralizing antibody significantly reduced calcium accumulation.

Primary cilia defects causing mitral valve prolapse

Mitral valve prolapse (MVP) affects 1 in 40 people and is the most common indication for mitral valve surgery. MVP can cause arrhythmias, heart failure, and sudden cardiac death, and to date, the causes of this disease are poorly understood. We now demonstrate that defects in primary cilia genes and their regulated pathways can cause MVP in familial and sporadic nonsyndromic MVP cases. Our expression studies and genetic ablation experiments confirmed a role for primary cilia in regulating ECM deposition during cardiac development. Loss of primary cilia during development resulted in progressive myxomatous degeneration and profound mitral valve pathology in the adult setting. Analysis of a large family with inherited, autosomal dominant nonsyndromic MVP identified a deleterious missense mutation in a cilia gene, DZIP1 A mouse model harboring this variant confirmed the pathogenicity of this mutation and revealed impaired ciliogenesis during development, which progressed to adult myxomatous valve disease and functional MVP. Relevance of primary cilia in common forms of MVP was tested using pathway enrichment in a large population of patients with MVP and controls from previously generated genome-wide association studies (GWAS), which confirmed the involvement of primary cilia genes in MVP. Together, our studies establish a developmental basis for MVP through altered cilia-dependent regulation of ECM and suggest that defects in primary cilia genes can be causative to disease phenotype in some patients with MVP.

Genetic Association Analyses Highlight IL6, ALPL, and NAV1 As 3 New Susceptibility Genes Underlying Calcific Aortic Valve Stenosis

BACKGROUND

Calcific aortic valve stenosis (CAVS) is a frequent and life-threatening cardiovascular disease for which there is currently no medical treatment available. To date, only 2 genes, LPA and PALMD, have been identified as causal for CAVS. We aimed to identify additional susceptibility genes for CAVS.

METHODS

A GWAS (genome-wide association study) meta-analysis of 4 cohorts, totaling 5115 cases and 354 072 controls of European descent, was performed. A TWAS (transcriptome-wide association study) was completed to integrate transcriptomic data from 233 human aortic valves. A series of post-GWAS analyses were performed, including fine-mapping, colocalization, phenome-wide association studies, pathway, and tissue enrichment as well as genetic correlation with cardiovascular traits.

RESULTS

In the GWAS meta-analysis, 4 loci achieved genome-wide significance, including 2 new loci: IL6 (interleukin 6) on 7p15.3 and ALPL (alkaline phosphatase) on 1p36.12. A TWAS integrating gene expression from 233 human aortic valves identified NAV1 (neuron navigator 1) on 1q32.1 as a new candidate causal gene. The CAVS risk alleles were associated with higher mRNA expression of NAV1 in valve tissues. Fine-mapping identified rs1800795 as the most likely causal variant in the IL6 locus. The signal identified colocalizes with the expression of the IL6 RNA antisense in various tissues. Phenome-wide association analyses in the UK Biobank showed colocalized associations between the risk allele at the IL6 lead variant and higher eosinophil count, pulse pressure, systolic blood pressure, and carotid artery procedures, implicating modulation of the IL6 pathways. The risk allele at the NAV1 lead variant colocalized with higher pulse pressure and higher prevalence of carotid artery stenosis. Association results at the genome-wide scale indicated genetic correlation between CAVS, coronary artery disease, and cardiovascular risk factors.

CONCLUSIONS

Our study implicates 3 new genetic loci in CAVS pathogenesis, which constitute novel targets for the development of therapeutic agents.

Population structure of modern-day Italians reveals patterns of ancient and archaic ancestries in Southern Europe

European populations display low genetic differentiation as the result of long-term blending of their ancient founding ancestries. However, it is unclear how the combination of ancient ancestries related to early foragers, Neolithic farmers, and Bronze Age nomadic pastoralists can explain the distribution of genetic variation across Europe. Populations in natural crossroads like the Italian peninsula are expected to recapitulate the continental diversity, but have been systematically understudied. Here, we characterize the ancestry profiles of Italian populations using a genome-wide dataset representative of modern and ancient samples from across Italy, Europe, and the rest of the world. Italian genomes capture several ancient signatures, including a non-steppe contribution derived ultimately from the Caucasus. Differences in ancestry composition, as the result of migration and admixture, have generated in Italy the largest degree of population structure detected so far in the continent, as well as shaping the amount of Neanderthal DNA in modern-day populations.

Genome-wide association study identifies eight risk loci and implicates metabo-psychiatric origins for anorexia nervosa

Characterized primarily by a low body-mass index, anorexia nervosa is a complex and serious illness1, affecting 0.9-4% of women and 0.3% of men2-4, with twin-based heritability estimates of 50-60%5. Mortality rates are higher than those in other psychiatric disorders6, and outcomes are unacceptably poor7. Here we combine data from the Anorexia Nervosa Genetics Initiative (ANGI)8,9 and the Eating Disorders Working Group of the Psychiatric Genomics Consortium (PGC-ED) and conduct a genome-wide association study of 16,992 cases of anorexia nervosa and 55,525 controls, identifying eight significant loci. The genetic architecture of anorexia nervosa mirrors its clinical presentation, showing significant genetic correlations with psychiatric disorders, physical activity, and metabolic (including glycemic), lipid and anthropometric traits, independent of the effects of common variants associated with body-mass index. These results further encourage a reconceptualization of anorexia nervosa as a metabo-psychiatric disorder. Elucidating the metabolic component is a critical direction for future research, and paying attention to both psychiatric and metabolic components may be key to improving outcomes.

Genetic Variation in LPA, Calcific Aortic Valve Stenosis in Patients Undergoing Cardiac Surgery, and Familial Risk of Aortic Valve Microcalcification

Importance

Genetic variants at the LPA locus are associated with both calcific aortic valve stenosis (CAVS) and coronary artery disease (CAD). Whether these variants are associated with CAVS in patients with CAD vs those without CAD is unknown.

Objective

To study the associations of LPA variants with CAVS in a cohort of patients undergoing heart surgery and LPA with CAVS in patients with CAD vs those without CAD and to determine whether first-degree relatives of patients with CAVS and high lipoprotein(a) (Lp[a]) levels showed evidence of aortic valve microcalcification.

Design, Setting, and Participants

This genetic association study included patients undergoing cardiac surgery from the Genome-Wide Association Study on Calcific Aortic Valve Stenosis in Quebec (QUEBEC-CAVS) study and patients with CAD, patients without CAD, and control participants from 6 genetic association studies: the UK Biobank, the European Prospective Investigation of Cancer (EPIC)-Norfolk, and Genetic Epidemiology Research on Aging (GERA) studies and 3 French cohorts. In addition, a family study included first-degree relatives of patients with CAVS. Data were collected from January 1993 to September 2018, and analysis was completed from September 2017 to September 2018.

Exposures

Case-control studies.

Main Outcomes and Measures

Presence of CAVS according to a weighted genetic risk score based on 3 common Lp(a)-raising variants and aortic valve microcalcification, defined as the mean tissue to background ratio of 1.25 or more, measured by fluorine 18-labeled sodium fluoride positron emission tomography/computed tomography.

Results

This study included 1009 individuals undergoing cardiac surgery and 1017 control participants in the QUEBEC-CAVS cohort; 3258 individuals with CAVS and CAD, 41 100 controls with CAD, 2069 individuals with CAVS without CAD, and 380 075 control participants without CAD in the UK Biobank, EPIC-Norfolk, and GERA studies and 3 French cohorts combined; and 33 first-degree relatives of 17 patients with CAVS and high Lp(a) levels (≥60 mg/dL) and 23 control participants with normal Lp(a) levels (<60 mg/dL). In the QUEBEC-CAVS study, each SD increase of the genetic risk score was associated with a higher risk of CAVS (odds ratio [OR], 1.35 [95% CI, 1.10-1.66]; P = .003). Each SD increase of the genetic risk score was associated with a higher risk of CAVS in patients with CAD (OR, 1.30 [95% CI, 1.20-1.42]; P < .001) and without CAD (OR, 1.33 [95% CI, 1.14-1.55]; P < .001). The percentage of individuals with a tissue to background ratio of 1.25 or more or CAVS was higher in first-degree relatives of patients with CAVS and high Lp(a) (16 of 33 [49%]) than control participants (3 of 23 [13%]; P = .006).

Conclusions and Relevance

In this study, a genetically elevated Lp(a) level was associated with CAVS independently of the presence of CAD. These findings support further research on the potential usefulness of Lp(a) cascade screening in CAVS.

Genome-Wide Association Study-Driven Gene-Set Analyses, Genetic, and Functional Follow-Up Suggest GLIS1 as a Susceptibility Gene for Mitral Valve Prolapse

Background Mitral valve prolapse (MVP) is a common heart valve disease, the most frequent indication for valve repair or replacement. MVP is characterized by excess extracellular matrix secretion and cellular disorganization, which leads to bulky valves that are unable to coapt correctly during ventricular systole resulting in mitral regurgitation, and it is associated with sudden cardiac death. Here we aim to characterize globally the biological mechanisms underlying genetic susceptibility to MVP to better characterize its triggering mechanisms. Methods We applied i-GSEA4GWAS and DEPICT, two pathway enrichment tools to MVP genome-wide association studies. We followed-up the association with MVP in an independent dataset of cases and controls. This research was conducted using the UK Biobank Resource. Immunohistochemistry staining for Glis1 (GLIS family zinc finger 1) was conducted in developing heart of mice. Knockdown of Glis1 using morpholinos was performed in zebrafish animals 72 hours postfertilization. Results We show that genes at risk loci are involved in biological functions relevant to actin filament organization, cytoskeleton biology, and cardiac development. The enrichment for positive regulation of transcription, cell proliferation, and migration motivated the follow-up of GLIS1, a transcription factor from the Krüppel-like zinc finger family. In combination with previously available data, we now report a genome-wide significant association with MVP (odds ratio, 1.20; P=4.36×10-10), indicating that Glis1 is expressed during embryonic development predominantly in nuclei of endothelial and interstitial cells of mitral valves in mouse. We also show that Glis1 knockdown causes atrioventricular regurgitation in developing hearts in zebrafish. Conclusions Our findings define globally molecular and cellular mechanisms underlying common genetic susceptibility to MVP and implicate established and unprecedented mechanisms. Through the GLIS1 association and function, we point at regulatory functions during cardiac development as common mechanisms to mitral valve degeneration.

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.

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.

An Expanded Genome-Wide Association Study of Type 2 Diabetes in Europeans

To characterize type 2 diabetes (T2D)-associated variation across the allele frequency spectrum, we conducted a meta-analysis of genome-wide association data from 26,676 T2D case and 132,532 control subjects of European ancestry after imputation using the 1000 Genomes multiethnic reference panel. Promising association signals were followed up in additional data sets (of 14,545 or 7,397 T2D case and 38,994 or 71,604 control subjects). We identified 13 novel T2D-associated loci (P < 5 × 10^-8), including variants near the GLP2R, GIP, and HLA-DQA1 genes. Our analysis brought the total number of independent T2D associations to 128 distinct signals at 113 loci. Despite substantially increased sample size and more complete coverage of low-frequency variation, all novel associations were driven by common single nucleotide variants. Credible sets of potentially causal variants were generally larger than those based on imputation with earlier reference panels, consistent with resolution of causal signals to common risk haplotypes. Stratification of T2D-associated loci based on T2D-related quantitative trait associations revealed tissue-specific enrichment of regulatory annotations in pancreatic islet enhancers for loci influencing insulin secretion and in adipocytes, monocytes, and hepatocytes for insulin action-associated loci. These findings highlight the predominant role played by common variants of modest effect and the diversity of biological mechanisms influencing T2D pathophysiology.

Identification of novel risk loci for restless legs syndrome in genome-wide association studies in individuals of European ancestry: a meta-analysis

BACKGROUND

Restless legs syndrome is a prevalent chronic neurological disorder with potentially severe mental and physical health consequences. Clearer understanding of the underlying pathophysiology is needed to improve treatment options. We did a meta-analysis of genome-wide association studies (GWASs) to identify potential molecular targets.

METHODS

In the discovery stage, we combined three GWAS datasets (EU-RLS GENE, INTERVAL, and 23andMe) with diagnosis data collected from 2003 to 2017, in face-to-face interviews or via questionnaires, and involving 15 126 cases and 95 725 controls of European ancestry. We identified common variants by fixed-effect inverse-variance meta-analysis. Significant genome-wide signals (p≤5 × 10-8) were tested for replication in an independent GWAS of 30 770 cases and 286 913 controls, followed by a joint analysis of the discovery and replication stages. We did gene annotation, pathway, and gene-set-enrichment analyses and studied the genetic correlations between restless legs syndrome and traits of interest.

FINDINGS

We identified and replicated 13 new risk loci for restless legs syndrome and confirmed the previously identified six risk loci. MEIS1 was confirmed as the strongest genetic risk factor for restless legs syndrome (odds ratio 1·92, 95% CI 1·85-1·99). Gene prioritisation, enrichment, and genetic correlation analyses showed that identified pathways were related to neurodevelopment and highlighted genes linked to axon guidance (associated with SEMA6D), synapse formation (NTNG1), and neuronal specification (HOXB cluster family and MYT1).

INTERPRETATION

Identification of new candidate genes and associated pathways will inform future functional research. Advances in understanding of the molecular mechanisms that underlie restless legs syndrome could lead to new treatment options. We focused on common variants; thus, additional studies are needed to dissect the roles of rare and structural variations.

FUNDING

Deutsche Forschungsgemeinschaft, Helmholtz Zentrum München-Deutsches Forschungszentrum für Gesundheit und Umwelt, National Research Institutions, NHS Blood and Transplant, National Institute for Health Research, British Heart Foundation, European Commission, European Research Council, National Institutes of Health, National Institute of Neurological Disorders and Stroke, NIH Research Cambridge Biomedical Research Centre, and UK Medical Research Council.

Significant Locus and Metabolic Genetic Correlations Revealed in Genome-Wide Association Study of Anorexia Nervosa

OBJECTIVE

The authors conducted a genome-wide association study of anorexia nervosa and calculated genetic correlations with a series of psychiatric, educational, and metabolic phenotypes.

METHOD

Following uniform quality control and imputation procedures using the 1000 Genomes Project (phase 3) in 12 case-control cohorts comprising 3,495 anorexia nervosa cases and 10,982 controls, the authors performed standard association analysis followed by a meta-analysis across cohorts. Linkage disequilibrium score regression was used to calculate genome-wide common variant heritability (single-nucleotide polymorphism [SNP]-based heritability [h²_SNP]), partitioned heritability, and genetic correlations (r_g) between anorexia nervosa and 159 other phenotypes.

RESULTS

Results were obtained for 10,641,224 SNPs and insertion-deletion variants with minor allele frequencies >1% and imputation quality scores >0.6. The h²_SNP of anorexia nervosa was 0.20 (SE=0.02), suggesting that a substantial fraction of the twin-based heritability arises from common genetic variation. The authors identified one genome-wide significant locus on chromosome 12 (rs4622308) in a region harboring a previously reported type 1 diabetes and autoimmune disorder locus. Significant positive genetic correlations were observed between anorexia nervosa and schizophrenia, neuroticism, educational attainment, and high-density lipoprotein cholesterol, and significant negative genetic correlations were observed between anorexia nervosa and body mass index, insulin, glucose, and lipid phenotypes.

CONCLUSIONS

Anorexia nervosa is a complex heritable phenotype for which this study has uncovered the first genome-wide significant locus. Anorexia nervosa also has large and significant genetic correlations with both psychiatric phenotypes and metabolic traits. The study results encourage a reconceptualization of this frequently lethal disorder as one with both psychiatric and metabolic etiology.

DoEstRare: A statistical test to identify local enrichments in rare genomic variants associated with disease

Next-generation sequencing technologies made it possible to assay the effect of rare variants on complex diseases. As an extension of the "common disease-common variant" paradigm, rare variant studies are necessary to get a more complete insight into the genetic architecture of human traits. Association studies of these rare variations show new challenges in terms of statistical analysis. Due to their low frequency, rare variants must be tested by groups. This approach is then hindered by the fact that an unknown proportion of the variants could be neutral. The risk level of a rare variation may be determined by its impact but also by its position in the protein sequence. More generally, the molecular mechanisms underlying the disease architecture may involve specific protein domains or inter-genic regulatory regions. While a large variety of methods are optimizing functionality weights for each single marker, few evaluate variant position differences between cases and controls. Here, we propose a test called DoEstRare, which aims to simultaneously detect clusters of disease risk variants and global allele frequency differences in genomic regions. This test estimates, for cases and controls, variant position densities in the genetic region by a kernel method, weighted by a function of allele frequencies. We compared DoEstRare with previously published strategies through simulation studies as well as re-analysis of real datasets. Based on simulation under various scenarios, DoEstRare was the sole to consistently show highest performance, in terms of type I error and power both when variants were clustered or not. DoEstRare was also applied to Brugada syndrome and early-onset Alzheimer's disease data and provided complementary results to other existing tests. DoEstRare, by integrating variant position information, gives new opportunities to explain disease susceptibility. DoEstRare is implemented in a user-friendly R package.

Dysfunction of the Voltage-Gated K+ Channel β2 Subunit in a Familial Case of Brugada Syndrome

BACKGROUND

The Brugada syndrome is an inherited cardiac arrhythmia associated with high risk of sudden death. Although 20% of patients with Brugada syndrome carry mutations in SCN5A, the molecular mechanisms underlying this condition are still largely unknown.

METHODS AND RESULTS

We combined whole-exome sequencing and linkage analysis to identify the genetic variant likely causing Brugada syndrome in a pedigree for which SCN5A mutations had been excluded. This approach identified 6 genetic variants cosegregating with the Brugada electrocardiographic pattern within the pedigree. In silico gene prioritization pointed to 1 variant residing in KCNAB2, which encodes the voltage-gated K(+) channel β2-subunit (Kvβ2-R12Q). Kvβ2 is widely expressed in the human heart and has been shown to interact with the fast transient outward K(+) channel subunit Kv4.3, increasing its current density. By targeted sequencing of the KCNAB2 gene in 167 unrelated patients with Brugada syndrome, we found 2 additional rare missense variants (L13F and V114I). We then investigated the physiological effects of the 3 KCNAB2 variants by using cellular electrophysiology and biochemistry. Patch-clamp experiments performed in COS-7 cells expressing both Kv4.3 and Kvβ2 revealed a significant increase in the current density in presence of the R12Q and L13F Kvβ2 mutants. Although biotinylation assays showed no differences in the expression of Kv4.3, the total and submembrane expression of Kvβ2-R12Q were significantly increased in comparison with wild-type Kvβ2.

CONCLUSIONS

Altogether, our results indicate that Kvβ2 dysfunction can contribute to the Brugada electrocardiographic pattern.

Endo and exo aortic morphometry (endoluminal) in the aortic arch and its branches

Our study was conducted by the examination of angioCT’s, performing external measurements of aortic arch and the exo and endoaortic measurements of its three. The diameter of the thoracic aorta prior to the origin of the brachiocephalic arterial trunk was found with an average of 31.65 mm and below the left subclavian artery origin we found an average diameter of 24.3 mm. The brachiocephalic arterial trunk had an average diameter of 11.575 mm, 6.05 at carotid artery level and 9.05 mm at the level of the left subclavian artery. The endoaortic average diameter of the brachiocephalic arterial trunk: horizontally, 13.0 mm and vertically 11.7 mm; left common carotid artery horizontal diameter was 10.5 mm and 9.7 mm vertically and the left subclavian artery have 14.1 mm horizontally and 10.8 mm vertically.