Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations.

Osteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100 independently associated risk variants across 11 osteoarthritis phenotypes, 52 of which have not been associated with the disease before. We report thumb and spine osteoarthritis risk variants and identify differences in genetic effects between weight-bearing and non-weight-bearing joints. We identify sex-specific and early age-at-onset osteoarthritis risk loci. We integrate functional genomics data from primary patient tissues (including articular cartilage, subchondral bone, and osteophytic cartilage) and identify high-confidence effector genes. We provide evidence for genetic correlation with phenotypes related to pain, the main disease symptom, and identify likely causal genes linked to neuronal processes. Our results provide insights into key molecular players in disease processes and highlight attractive drug targets to accelerate translation.

Identifying causal serum protein-cardiometabolic trait relationships using whole genome sequencing.

Cardiometabolic diseases, such as type 2 diabetes and cardiovascular disease, have a high public health burden. Understanding the genetically determined regulation of proteins that are dysregulated in disease can help to dissect the complex biology underpinning them. Here, we perform a protein quantitative trait locus (pQTL) analysis of 248 serum proteins relevant to cardiometabolic processes in 2893 individuals. Meta-analyzing whole-genome sequencing (WGS) data from two Greek cohorts, MANOLIS (n = 1356; 22.5× WGS) and Pomak (n = 1537; 18.4× WGS), we detect 301 independently associated pQTL variants for 170 proteins, including 12 rare variants (minor allele frequency < 1%). We additionally find 15 pQTL variants that are rare in non-Finnish European populations but have drifted up in the frequency in the discovery cohorts here. We identify proteins causally associated with cardiometabolic traits, including Mep1b for high-density lipoprotein (HDL) levels, and describe a knock-out (KO) Mep1b mouse model. Our findings furnish insights into the genetic architecture of the serum proteome, identify new protein-disease relationships and demonstrate the importance of isolated populations in pQTL analysis.

Transcriptomic and proteomic analysis of iris tissue and aqueous humor in juvenile idiopathic arthritis-associated uveitis.

Gene and protein expression profiles of iris biopsies, aqueous humor (AqH), and sera in patients with juvenile idiopathic arthritis-associated uveitis (JIAU) in comparison to control patients with primary open-angle glaucoma (POAG) and HLA-B27-positive acute anterior uveitis (AAU) were investigated. Via RNA Sequencing (RNA-Seq) and mass spectrometry-based protein expression analyses 136 genes and 56 proteins could be identified as being significantly differentially expressed (DE) between the JIAU and POAG group. Gene expression of different immunoglobulin (Ig) components as well as of the B cell-associated factors ID3, ID1, and EBF1 was significantly upregulated in the JIAU group as compared to POAG patients. qRT-PCR analysis showed a significantly higher gene expression of the B cell-related genes CD19, CD20, CD27, CD138, and MZB1 in the JIAU group. At the protein level, a significantly higher expression of Ig components in JIAU than in POAG was confirmed. The B cell-associated protein MZB1 showed a higher expression in JIAU patients than in POAG which was confirmed by western blot analysis. Using bead-based immunoassay analysis we were able to detect a significantly higher concentration of the B cell-activating and survival factors BAFF, APRIL, and IL-6 in the AqH of JIAU and AAU patients than in POAG patients. The intraocularly upregulated B cell-specific genes and proteins in iris tissue suggest that B cells participate in the immunopathology of JIAU. The intracameral environment in JIAU may facilitate local effector and survival functions of B cells, leading to disease course typical for anterior uveitis.

Rare genetic variants in SMAP1, B3GAT2, and RIMS1 contribute to pediatric venous thromboembolism.

Recent genome-wide association studies (GWAS) have confirmed known risk mutations for venous thromboembolism (VTE) and identified a number of novel susceptibility loci in adults. Here we present a GWAS in 212 nuclear families with pediatric VTE followed by targeted next-generation sequencing (NGS) to identify causative mutations contributing to the association. Three single nucleotide polymorphisms (SNPs) exceeded the threshold for genome-wide significance as determined by permutation testing using 100 000 bootstrap permutations (P < 10-5). These SNPs reside in a region on chromosome 6q13 comprising the genes small ARF GAP1 (SMAP1), an ARF6 guanosine triphosphatase-activating protein that functions in clathrin-dependent endocytosis, and ß-1,3-glucoronyltransferase 2 (B3GAT2), a member of the human natural killer 1 carbohydrate pathway. Rs1304029 and rs2748331 are associated with pediatric VTE with unpermuted/permuted values of P = 1.42 × 10-6/2.0 × 10-6 and P = 6.11 × 10-6/1.8 × 10-5, respectively. Rs2748331 was replicated (P = .00719) in an independent study sample coming from our GWAS on pediatric thromboembolic stroke (combined P = 7.88 × 10-7). Subsequent targeted NGS in 24 discordant sibling pairs identified 17 nonsynonymous coding variants, of which 1 located in SMAP1 and 3 in RIMS1, a member of the RIM family of active zone proteins, are predicted as damaging by Protein Variation Effect Analyzer and/or sorting intolerant from tolerant scores. Three SNPs curtly missed statistical significance in the transmission-disequilibrium test in the full cohort (rs112439957: P = .08326, SMAP1; rs767118962: P = .08326, RIMS1; and rs41265501: P = .05778, RIMS1). In conjunction, our data provide compelling evidence for SMAP1, B3GAT2, and RIMS1 as novel susceptibility loci for pediatric VTE and warrant future functional studies to unravel the underlying molecular mechanisms leading to VTE.

Standing genetic variation affects phenotypic heterogeneity in an SCN5A-mutation founder population with excess sudden cardiac death.

BACKGROUND: The Worm Study, ascertained from a multigeneration pedigree segregating a single amino acid deletion in SCN5A (c.4850_4852delTCT, p.(Phe1617del), rs749697698), is characterized by substantial phenotypic heterogeneity and overlap of sudden cardiac death, long-QT syndrome, cardiac conduction disease, Brugada syndrome, and isorhythmic atrioventricular dissociation. Linkage analysis for a synthetic trait derived from these phenotypes identified a single peak (logarithm of the odds [LOD] = 4.52) at the SCN5A/SCN10A/SCN11A locus on chromosome 3. OBJECTIVE: This study explored the role of additional genetic variation in the chromosome 3 locus as a source of phenotypic heterogeneity in the Worm Study population. METHODS: Genotypes underlying the linkage peak (n = 70) were characterized using microarrays. Haplotypes were determined using family-aware phasing and a population-specific reference panel. Variants with minor allele frequencies >0.10 were tested for association with cardiac conduction disease and isorhythmic dissociation using LAMP and logistic regression. RESULTS: Only 1 haplotype carried the p.Phe1617del/rs749697698 deletion, suggesting relatively recent development (∼18 generations); this haplotype contained 5 other missense variants spanning SCN5A/SCN10A/SCN11A. Noncarrier haplotypes (n = 74) ranged in frequency from 0.5% to 5%. Although no variants were associated with cardiac conduction disease, a homozygous missense variant in SCN10A was associated with isorhythmic dissociation after correction for multiple comparisons (odds ratio 11.23; 95% confidence interval 2.76-23.39; P = 1.2 × 10-4). This variant (rs12632942) was previously associated with PR interval. CONCLUSION: Our data suggest that other variants, alongside a pathogenic mutation, are associated with phenotypic heterogeneity. Single-mutation screening may be insufficient to predict electrical heart disease in patients and family members. In the Worm Study population, segregating a pathogenic SCN5A mutation, compound variation in the SCN5A/SCN10A/SCN11A locus determines arrhythmic outcome.

Gene-based whole genome sequencing meta-analysis of 250 circulating proteins in three isolated European populations.

OBJECTIVE: Deep sequencing offers unparalleled access to rare variants in human populations. Understanding their role in disease is a priority, yet prohibitive sequencing costs mean that many cohorts lack the sample size to discover these effects on their own. Meta-analysis of individual variant scores allows the combination of rare variants across cohorts and study of their aggregated effect at the gene level, boosting discovery power. However, the methods involved have largely not been field-tested. In this study, we aim to perform the first meta-analysis of gene-based rare variant aggregation optimal tests, applied to the human cardiometabolic proteome. METHODS: Here, we carry out this analysis across MANOLIS, Pomak and ORCADES, three isolated European cohorts with whole-genome sequencing (total N = 4,422). We examine the genetic architecture of 250 proteomic traits of cardiometabolic relevance. We use a containerised pipeline to harmonise variant lists across cohorts and define four sets of qualifying variants. For every gene, we interrogate protein-damaging variants, exonic variants, exonic and regulatory variants, and regulatory only variants, using the CADD and Eigen scores to weigh variants according to their predicted functional consequence. We perform single-cohort rare variant analysis and meta-analyse variant scores using the SMMAT package. RESULTS: We describe 5 rare variant pQTLs (RV-pQTL) which pass our stringent significance threshold (7.45 × 10-11) and quality control procedure. These were split between four cis signals for MARCO, TEK, MMP2 and MPO, and one trans association for GDF2 in the SERPINA11 gene. We show that the cis-MPO association, which was not detectable using the single-point data alone, is driven by 5 missense and frameshift variants. These include rs140636390 and rs119468010, which are specific to MANOLIS and ORCADES, respectively. We show how this kind of signal could improve the predictive accuracy of genetic factors in common complex disease such as stroke and cardiovascular disease. CONCLUSIONS: Our proof-of-concept study demonstrates the power of gene-based meta-analyses for discovering disease-relevant associations complementing common-variant signals by incorporating population-specific rare variation.

Mapping the serum proteome to neurological diseases using whole genome sequencing.

Despite the increasing global burden of neurological disorders, there is a lack of effective diagnostic and therapeutic biomarkers. Proteins are often dysregulated in disease and have a strong genetic component. Here, we carry out a protein quantitative trait locus analysis of 184 neurologically-relevant proteins, using whole genome sequencing data from two isolated population-based cohorts (N = 2893). In doing so, we elucidate the genetic landscape of the circulating proteome and its connection to neurological disorders. We detect 214 independently-associated variants for 107 proteins, the majority of which (76%) are cis-acting, including 114 variants that have not been previously identified. Using two-sample Mendelian randomisation, we identify causal associations between serum CD33 and Alzheimer's disease, GPNMB and Parkinson's disease, and MSR1 and schizophrenia, describing their clinical potential and highlighting drug repurposing opportunities.

ADAMTS12, a new candidate gene for pediatric stroke.

We recently reported a family-based genome wide association study (GWAS) for pediatric stroke pointing our attention to two significantly associated genes of the ADAMTS (a disintegrin and metalloproteinase with thrombospondin motifs) gene family ADAMTS2 (rs469568, p = 8x10-6) and ADAMTS12 (rs1364044, p = 2.9x10-6). To further investigate these candidate genes, we applied a targeted resequencing approach on 48 discordant sib-pairs for pediatric stroke followed by genotyping of the detected non-synonymous variants in the full cohort of 270 offspring trios and subsequent fine mapping analysis. We identified eight non-synonymous SNPs in ADAMTS2 and six in ADAMTS12 potentially influencing the respective protein function. These variants were genotyped within a cohort of 270 affected offspring trios, association analysis revealed the ADAMTS12 variant rs77581578 to be significantly under-transmitted (p = 6.26x10-3) to pediatric stroke patients. The finding was validated in a pediatric venous thromboembolism (VTE) cohort of 189 affected trios. Subsequent haplotype analysis of ADAMTS12 detected a significantly associated haplotype comprising the originally identified GWAS variant. Several ADAMTS genes such as ADAMTS13 are involved in thromboembolic disease process. Here, we provide further evidence for ADAMTS12 to likely play a role in pediatric stroke. Further functional studies are warranted to assess the functional role of ADAMTS12 in the pathogenesis of stroke.

Whole-genome sequencing analysis of the cardiometabolic proteome.

The human proteome is a crucial intermediate between complex diseases and their genetic and environmental components, and an important source of drug development targets and biomarkers. Here, we comprehensively assess the genetic architecture of 257 circulating protein biomarkers of cardiometabolic relevance through high-depth (22.5×) whole-genome sequencing (WGS) in 1328 individuals. We discover 131 independent sequence variant associations (P < 7.45 × 10-11) across the allele frequency spectrum, all of which replicate in an independent cohort (n = 1605, 18.4x WGS). We identify for the first time replicating evidence for rare-variant cis-acting protein quantitative trait loci for five genes, involving both coding and noncoding variation. We construct and validate polygenic scores that explain up to 45% of protein level variation. We find causal links between protein levels and disease risk, identifying high-value biomarkers and drug development targets.

Analysis of the genomic architecture of a complex trait locus in hypertensive rat models links Tmem63c to kidney damage.

Unraveling the genetic susceptibility of complex diseases such as chronic kidney disease remains challenging. Here, we used inbred rat models of kidney damage associated with elevated blood pressure for the comprehensive analysis of a major albuminuria susceptibility locus detected in these models. We characterized its genomic architecture by congenic substitution mapping, targeted next-generation sequencing, and compartment-specific RNA sequencing analysis in isolated glomeruli. This led to prioritization of transmembrane protein Tmem63c as a novel potential target. Tmem63c is differentially expressed in glomeruli of allele-specific rat models during onset of albuminuria. Patients with focal segmental glomerulosclerosis exhibited specific TMEM63C loss in podocytes. Functional analysis in zebrafish revealed a role for tmem63c in mediating the glomerular filtration barrier function. Our data demonstrate that integrative analysis of the genomic architecture of a complex trait locus is a powerful tool for identification of new targets such as Tmem63c for further translational investigation.

Targeted resequencing of a locus for heparin-induced thrombocytopenia on chromosome 5 identified in a genome-wide association study.

Immune-mediated heparin-induced thrombocytopenia (HIT) is the clinically most important adverse drug reaction (ADR) in response to heparin therapy characterized by a prothrombotic state despite a decrease in platelet count. We conducted a genome-wide association study in 96 suspected HIT cases and 96 controls to explore the genetic predisposition for HIT within a case-control pharmacovigilance study followed by replication in additional 86 cases and 86 controls from the same study. One single nucleotide polymorphism (SNP, rs1433265, P = 6.5 × 10-5, odds ratio (OR) 2.79) from 16 identified SNPs was successfully replicated (P = 1.5 × 10-4, OR 2.77; combined data set P = 2.7 × 10-8, OR 2.77) and remained the most strongly associated SNP after imputing locus genotypes. Fine mapping revealed a significantly associated risk-conferring haplotype (P = 4.9 × 10-6, OR 2.41). In order to find rare variants contributing to the association signals, we applied a targeted resequencing approach in a subgroup of 73 HIT patients and 23 controls for the regions with the 16 most strongly HIT-associated SNPs. C-alpha testing was applied to test for the impact of rare variants and we detected two candidate genes, the discoidin domain receptor tyrosine kinase 1 (DDR1, P = 3.6 × 10-2) and the multiple C2 and transmembrane domain containing 2 (MCTP2, P = 4.5 × 10-2). For the genes interactor of little elongation complex ELL subunit 1 (ICE1) and a disintegrin-like and metalloproteinase with thrombospondin type 1 motif, 16 (ADAMTS16) nearby rs1433265, we identified several missense variants. Although replication in an independent population is warranted, these findings provide a basis for future studies aiming to identify and characterize genetic susceptibility factors for HIT. KEY MESSAGES: We identified and validated a HIT-associated locus on chromosome 5. Targeted NGS analysis for rare variants identifies DDR1 and MCTP2 as novel candidates. In addition, missense variants for ADAMTS16 and ICE1 were identified in the locus.

Heritability in a SCN5A-mutation founder population with increased female susceptibility to non-nocturnal ventricular tachyarrhythmia and sudden cardiac death.

BACKGROUND: Heritable cardiac-sodium channel dysfunction is associated with various arrhythmia syndromes, some predisposing to ventricular fibrillation. Phenotypic diversity among carriers of identical-by-descent mutations is often remarkable, suggesting influences of genetic modifiers. OBJECTIVE: The purpose of this study was to identify a unique SCN5A-mutation founder population with mixed clinical phenotypes and sudden cardiac death, and to investigate the heritability of electromechanical traits besides the SCN5A-mutation effect. METHODS: The 16-generation founder population segregating SCN5A c.4850_4852delTCT, p.(Phe1617del), was comprehensively phenotyped. Variance component analysis was used to evaluate the mutation's effects and assess heritability. RESULTS: In 45 p.(Phe1617del) positives, the mutation associated strongly with QTc prolongation (472 ± 60 ms vs 423 ± 35 ms in 26 mutation negatives; P <.001; odds ratio for long-QT syndrome 22.4; 95% confidence interval 4.5-224.2; P <.001) and electromechanical window (EMW) negativity (-29 ± 47 ms vs 34 ± 26 ms; P <.001). Overlapping phenotypes including conduction delay and Brugada syndrome were noted in 19. Polymorphic ventricular tachyarrhythmias occurred mostly in the daytime, after arousal-evoked heart-rate acceleration and repolarization prolongation. Cox proportional hazards regression analysis revealed female gender as an independent risk factor for cardiac events (hazard ratio 5.1; 95% confidence interval 1.6-16.3; P = .006). p.(Phe1617del) was an important determinant of QTcbaseline, QTcmax, and EMW, explaining 18%, 28%, and 37%, respectively, of the trait's variance. Significant heritability was observed for PQ interval (P = .003) after accounting for the p.(Phe1617del) effect. CONCLUSION: This SCN5A-p.(Phe1617del) founder population with phenotypic divergence and overlap reveals long-QT syndrome-related and arousal-evoked ventricular tachyarrhythmias with a female preponderance. Variance component analysis indicates additional genetic variance for PQ interval hidden in the genome, besides a dominant p.(Phe1617del) effect on QTc and EMW.

Rare Variants in the ADAMTS13 Von Willebrand Factor-Binding Domain Contribute to Pediatric Stroke.

BACKGROUND: Recently, we reported a gene network of ADAMTS (A Disintegrin-like and Metalloprotease with Thrombospondin motifs) genes as central component of the genetic risk contributing to pediatric stroke. ADAMTS13 is a prime example for such a key component as it cleaves von Willebrand factor multimers, reduces platelet adhesion and aggregation, and downregulates thrombus formation and inflammation. METHODS AND RESULTS: We characterized the genetic architecture of ADAMTS13 through targeted next-generation sequencing of 48 affected children and their unaffected siblings and identified in total 241 variants (single nucleotide polymorphisms or insertions/deletions) in the ADAMTS13 gene. From these, based on significance in the sibship disequilibrium test (P<0.05) or protein-altering properties, we selected 21 common variants covering the complete ADAMTS13 gene for genotyping in 270 trios and subsequent association analyses. Transmission disequilibrium testing was performed for affection status and ADAMTS13 activity levels using PLINK and FBAT, respectively. Ten single nucleotide polymorphisms were significantly associated with pediatric stroke (P<0.05 to P<0.001), 2 of which (rs2285489 and rs28793911) were also significantly associated with ADAMTS13 levels (P=0.0004 and P=0.0092). The resulting protective haplotype H1.1. (T:U 95.5: 144.4; P=0.0016) is associated with increased ADAMTS13 levels (age-adjusted P=0.0108). Haplotype association using a sliding window approach assigns this association to the ADAMTS13 von Willebrand factor-binding domain (P=1.2×10(-4)). CONCLUSIONS: Our data provide a link between the genetic architecture of ADAMTS13, ADAMTS13 levels, and stroke susceptibility. Altogether, these studies render ADAMTS13 an attractive candidate for functional studies and may contribute to personalized diagnosis and treatment options in future.

An information theoretic thresholding method for detecting protein colocalizations in stacks of fluorescence images.

In this note, we present a new method that allows us to determine threshold values for separating presence and absence of proteins in a stack of fluorescence images describing a spatial distribution of proteins across a biological object (like a slice of nervous tissue, a sample of blood cells, etc.). We apply this method to stacks of fluorescence images and find that the resulting threshold values are almost identical with threshold values found using completely independent methods based on technological and biological aspects of the images in question.