Multi-ancestry polygenic risk scores for venous thromboembolism.

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality, with large disparities in incidence rates between Black and White Americans. Polygenic risk scores (PRSs) limited to variants discovered in genome-wide association studies in European-ancestry samples can identify European-ancestry individuals at high risk of VTE. However, there is limited evidence on whether high-dimensional PRS constructed using more sophisticated methods and more diverse training data can enhance the predictive ability and their utility across diverse populations. We developed PRSs for VTE using summary statistics from the International Network against Venous Thrombosis (INVENT) consortium genome-wide association studies meta-analyses of European- (71 771 cases and 1 059 740 controls) and African-ancestry samples (7482 cases and 129 975 controls). We used LDpred2 and PRS-CSx to construct ancestry-specific and multi-ancestry PRSs and evaluated their performance in an independent European- (6781 cases and 103 016 controls) and African-ancestry sample (1385 cases and 12 569 controls). Multi-ancestry PRSs with weights tuned in European-ancestry samples slightly outperformed ancestry-specific PRSs in European-ancestry test samples (e.g. the area under the receiver operating curve [AUC] was 0.609 for PRS-CSx_combinedEUR and 0.608 for PRS-CSxEUR [P = 0.00029]). Multi-ancestry PRSs with weights tuned in African-ancestry samples also outperformed ancestry-specific PRSs in African-ancestry test samples (PRS-CSxAFR: AUC = 0.58, PRS-CSx_combined AFR: AUC = 0.59), although this difference was not statistically significant (P = 0.34). The highest fifth percentile of the best-performing PRS was associated with 1.9-fold and 1.68-fold increased risk for VTE among European- and African-ancestry subjects, respectively, relative to those in the middle stratum. These findings suggest that the multi-ancestry PRS might be used to improve performance across diverse populations to identify individuals at highest risk for VTE.

A genetic association study of circulating coagulation factor VIII and von Willebrand factor levels.

ABSTRACT: Coagulation factor VIII (FVIII) and its carrier protein von Willebrand factor (VWF) are critical to coagulation and platelet aggregation. We leveraged whole-genome sequence data from the Trans-Omics for Precision Medicine (TOPMed) program along with TOPMed-based imputation of genotypes in additional samples to identify genetic associations with circulating FVIII and VWF levels in a single-variant meta-analysis, including up to 45 289 participants. Gene-based aggregate tests were implemented in TOPMed. We identified 3 candidate causal genes and tested their functional effect on FVIII release from human liver endothelial cells (HLECs) and VWF release from human umbilical vein endothelial cells. Mendelian randomization was also performed to provide evidence for causal associations of FVIII and VWF with thrombotic outcomes. We identified associations (P < 5 × 10-9) at 7 new loci for FVIII (ST3GAL4, CLEC4M, B3GNT2, ASGR1, F12, KNG1, and TREM1/NCR2) and 1 for VWF (B3GNT2). VWF, ABO, and STAB2 were associated with FVIII and VWF in gene-based analyses. Multiphenotype analysis of FVIII and VWF identified another 3 new loci, including PDIA3. Silencing of B3GNT2 and the previously reported CD36 gene decreased release of FVIII by HLECs, whereas silencing of B3GNT2, CD36, and PDIA3 decreased release of VWF by HVECs. Mendelian randomization supports causal association of higher FVIII and VWF with increased risk of thrombotic outcomes. Seven new loci were identified for FVIII and 1 for VWF, with evidence supporting causal associations of FVIII and VWF with thrombotic outcomes. B3GNT2, CD36, and PDIA3 modulate the release of FVIII and/or VWF in vitro.

CAVIN1-Mediated hERG Dynamics: A Novel Mechanism Underlying the Interindividual Variability in Drug-Induced Long QT.

BACKGROUND: Drug-induced QT prolongation (diLQT) is a feared side effect that could expose susceptible individuals to fatal arrhythmias. The occurrence of diLQT is primarily attributed to unintended drug interactions with cardiac ion channels, notably the hERG (human ether-a-go-go-related gene) channels that generate the delayed-rectifier potassium current (IKr) and thereby regulate the late repolarization phase. There is an important interindividual susceptibility to develop diLQT, which is of unknown origin but can be reproduced in patient-specific induced pluripotent stem cell-derived cardiomyocytes (iPS-CMs). We aimed to investigate the dynamics of hERG channels in response to sotalol and to identify regulators of the susceptibility to developing diLQT. METHODS: We measured electrophysiological activity and cellular distribution of hERG channels after hERG blocker treatment in iPS-CMs derived from patients with highest sensitivity (HS) or lowest sensitivity (LS) to sotalol administration in vivo (ie, on the basis of the measure of the maximal change in QT interval 3 hours after administration). Specific small interfering RNAs and CAVIN1-T2A-GFP adenovirus were used to manipulate CAVIN1 expression. RESULTS: Whereas HS and LS iPS-CMs showed similar electrophysiological characteristics at baseline, the late repolarization phase was prolonged and IKr significantly decreased after exposure of HS iPS-CMs to low sotalol concentrations. IKr reduction was caused by a rapid translocation of hERG channel from the membrane to the cytoskeleton-associated fractions upon sotalol application. CAVIN1, essential for caveolae biogenesis, was 2× more highly expressed in HS iPS-CMs, and its knockdown by small interfering RNA reduced their sensitivity to sotalol. CAVIN1 overexpression in LS iPS-CMs using adenovirus showed reciprocal effects. We found that treatment with sotalol promoted translocation of the hERG channel from the plasma membrane to the cytoskeleton fractions in a process dependent on CAVIN1 (caveolae associated protein 1) expression. CAVIN1 silencing reduced the number of caveolae at the membrane and abrogated the translocation of hERG channel in sotalol-treated HS iPS-CMs. CAVIN1 also controlled cardiomyocyte responses to other hERG blockers, such as E4031, vandetanib, and clarithromycin. CONCLUSIONS: Our study identifies unbridled turnover of the potassium channel hERG as a mechanism supporting the interindividual susceptibility underlying diLQT development and demonstrates how this phenomenon is finely tuned by CAVIN1.

Testing for association with rare variants in the coding and non-coding genome: RAVA-FIRST, a new approach based on CADD deleteriousness score.

Rare variant association tests (RVAT) have been developed to study the contribution of rare variants widely accessible through high-throughput sequencing technologies. RVAT require to aggregate rare variants in testing units and to filter variants to retain only the most likely causal ones. In the exome, genes are natural testing units and variants are usually filtered based on their functional consequences. However, when dealing with whole-genome sequence (WGS) data, both steps are challenging. No natural biological unit is available for aggregating rare variants. Sliding windows procedures have been proposed to circumvent this difficulty, however they are blind to biological information and result in a large number of tests. We propose a new strategy to perform RVAT on WGS data: "RAVA-FIRST" (RAre Variant Association using Functionally-InfoRmed STeps) comprising three steps. (1) New testing units are defined genome-wide based on functionally-adjusted Combined Annotation Dependent Depletion (CADD) scores of variants observed in the gnomAD populations, which are referred to as "CADD regions". (2) A region-dependent filtering of rare variants is applied in each CADD region. (3) A functionally-informed burden test is performed with sub-scores computed for each genomic category within each CADD region. Both on simulations and real data, RAVA-FIRST was found to outperform other WGS-based RVAT. Applied to a WGS dataset of venous thromboembolism patients, we identified an intergenic region on chromosome 18 enriched for rare variants in early-onset patients. This region that was missed by standard sliding windows procedures is included in a TAD region that contains a strong candidate gene. RAVA-FIRST enables new investigations of rare non-coding variants in complex diseases, facilitated by its implementation in the R package Ravages.

Novel Cardiokine GDF3 Predicts Adverse Fibrotic Remodeling After Myocardial Infarction.

BACKGROUND: Myocardial infarction (MI) induces a repair response that ultimately generates a stable fibrotic scar. Although the scar prevents cardiac rupture, an excessive profibrotic response impairs optimal recovery by promoting the development of noncontractile fibrotic areas. The mechanisms that lead to cardiac fibrosis are diverse and incompletely characterized. We explored whether the expansion of cardiac fibroblasts after MI can be regulated through a paracrine action of cardiac stromal cells. METHODS: We performed a bioinformatic secretome analysis of cardiac stromal PW1+ cells isolated from normal and post-MI mouse hearts to identify novel secreted proteins. Functional assays were used to screen secreted proteins that promote fibroblast proliferation. The expressions of candidates were subsequently analyzed in mouse and human hearts and plasmas. The relationship between levels of circulating protein candidates and adverse post-MI cardiac remodeling was examined in a cohort of 80 patients with a first ST-segment-elevation MI and serial cardiac magnetic resonance imaging evaluations. RESULTS: Cardiac stromal PW1+ cells undergo a change in paracrine behavior after MI, and the conditioned media from these cells induced a significant increase in the proliferation of fibroblasts. We identified a total of 12 candidates as secreted proteins overexpressed by cardiac PW1+ cells after MI. Among these factors, GDF3 (growth differentiation factor 3), a member of the TGF-ß (transforming growth factor-ß) family, was markedly upregulated in the ischemic hearts. Conditioned media specifically enriched with GDF3 induced fibroblast proliferation at a high level by stimulation of activin-receptor-like kinases. In line with the secretory nature of this protein, we next found that GDF3 can be detected in mice and human plasma samples, with a significant increase in the days after MI. In humans, higher GDF3 circulating levels (measured in the plasma at day 4 after MI) were significantly associated with an increased risk of adverse remodeling 6 months after MI (adjusted odds ratio, 1.76 [1.03-3.00]; P=0.037), including lower left ventricular ejection fraction and a higher proportion of akinetic segments. CONCLUSIONS: Our findings define a mechanism for the profibrotic action of cardiac stromal cells through secreted cardiokines, such as GDF3, a candidate marker of adverse fibrotic remodeling after MI. REGISTRATION: URL: https://www. CLINICALTRIALS: gov; Unique identifier: NCT01113268.

Impaired balance between neutrophil extracellular trap formation and degradation by DNases in COVID-19 disease.

BACKGROUND: Thrombo-inflammation and neutrophil extracellular traps (NETs) are exacerbated in severe cases of COVID-19, potentially contributing to disease exacerbation. However, the mechanisms underpinning this dysregulation remain elusive. We hypothesised that lower DNase activity may be associated with higher NETosis and clinical worsening in patients with COVID-19. METHODS: Biological samples were obtained from hospitalized patients (15 severe, 37 critical at sampling) and 93 non-severe ambulatory cases. Our aims were to compare NET biomarkers, functional DNase levels, and explore mechanisms driving any imbalance concerning disease severity. RESULTS: Functional DNase levels were diminished in the most severe patients, paralleling an imbalance between NET markers and DNase activity. DNase1 antigen levels were higher in ambulatory cases but lower in severe patients. DNase1L3 antigen levels remained consistent across subgroups, not rising alongside NET markers. DNASE1 polymorphisms correlated with reduced DNase1 antigen levels. Moreover, a quantitative deficiency in plasmacytoid dendritic cells (pDCs), which primarily express DNase1L3, was observed in critical patients. Analysis of public single-cell RNAseq data revealed reduced DNase1L3 expression in pDCs from severe COVID-19 patient. CONCLUSION: Severe and critical COVID-19 cases exhibited an imbalance between NET and DNase functional activity and quantity. Early identification of NETosis imbalance could guide targeted therapies against thrombo-inflammation in COVID-19-related sepsis, such as DNase administration, to avert clinical deterioration. TRIAL REGISTRATION: COVERAGE trial (NCT04356495) and COLCOV19-BX study (NCT04332016).

DNA-pools targeted-sequencing as a robust cost-effective method to detect rare variants: Application to dilated cardiomyopathy genetic diagnosis.

Dilated cardiomyopathy (DCM) is a heart disease characterized by left ventricular dilatation and systolic dysfunction. In 30% of cases, pathogenic variants, essentially private to each patient, are identified in at least one of almost 50 reported genes. Thus, while costly, exons capture-based Next Generation Sequencing (NGS) of a targeted gene panel appears as the best strategy to genetically diagnose DCM. Here, we report a NGS strategy applied to pools of 8 DNAs from DCM patients and validate its robustness for rare variants detection at 4-fold reduced cost. Our pipeline uses Freebayes to detect variants with the expected 1/16 allele frequency. From the whole set of detected rare variants in 96 pools we set the variants quality parameters optimizing true positives calling. When compared to simplex DNA sequencing in a shared subset of 50 DNAs, 96% of SNVs/InsDel were accurately identified in pools. Extended to the 384 DNAs included in the study, we detected 100 variants (ACMG class 4 and 5), mostly in well-known morbid gene causing DCM such as TTN, MYH7, FLNC, and TNNT2. To conclude, we report an original pool-sequencing NGS method accurately detecting rare variants. This innovative approach is cost-effective for genetic diagnostic in rare diseases.

Two SERPINC1 variants affecting N-glycosylation of Asn224 cause severe thrombophilia not detected by functional assays.

Antithrombin deficiency, the most severe congenital thrombophilia, might be underestimated, as some pathogenic variants are not detected by routine functional methods. We have identified 2 new SERPINC1 variants, p.Glu227Lys and p.Asn224His, in 4 unrelated thrombophilic patients with early and recurrent thrombosis that had normal antithrombin activity. In one case, the mutation was identified by whole genome sequencing, while in the 3 remaining cases, the mutation was identified by sequencing SERPINC1 based on a single functional positive finding supporting deficiency. The 2 variants shared a common functional defect, an impaired or null N-glycosylation of Asn224 according to a eukaryotic expression model. Carriers had normal anti-FXa or anti-FIIa activities but impaired anti-FVIIa activity and a detectable loss of inhibitory function when incubating the plasma for 1 hour at 41°C. Moreover, the ß glycoform of the variants, lacking 2 N-glycans, had reduced secretion, increased heparin affinity, no inhibitory activity, and a potential dominant-negative effect. These results explain the increased thrombin generation observed in carriers. Mutation experiments reflected the role that Lysine residues close to the N-glycosylation sequon have in impairing the efficacy of N-glycosylation. Our study shows new elements involved in the regulation of N-glycosylation, a key posttranslational modification that, according to our results, affects folding, secretion, and function, providing new evidence of the pathogenic consequence of an incorrect N-glycosylation of antithrombin. This study supports that antithrombin deficiency is underestimated and encourages the development of new functional and genetic tests to diagnose this severe thrombophilia.

Antithrombin, Protein C, and Protein S: Genome and Transcriptome-Wide Association Studies Identify 7 Novel Loci Regulating Plasma Levels.

BACKGROUND: Antithrombin, PC (protein C), and PS (protein S) are circulating natural anticoagulant proteins that regulate hemostasis and of which partial deficiencies are causes of venous thromboembolism. Previous genetic association studies involving antithrombin, PC, and PS were limited by modest sample sizes or by being restricted to candidate genes. In the setting of the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium, we meta-analyzed across ancestries the results from 10 genome-wide association studies of plasma levels of antithrombin, PC, PS free, and PS total. METHODS: Study participants were of European and African ancestries, and genotype data were imputed to TOPMed, a dense multiancestry reference panel. Each of the 10 studies conducted a genome-wide association studies for each phenotype and summary results were meta-analyzed, stratified by ancestry. Analysis of antithrombin included 25 243 European ancestry and 2688 African ancestry participants, PC analysis included 16 597 European ancestry and 2688 African ancestry participants, PSF and PST analysis included 4113 and 6409 European ancestry participants. We also conducted transcriptome-wide association analyses and multiphenotype analysis to discover additional associations. Novel genome-wide association studies and transcriptome-wide association analyses findings were validated by in vitro functional experiments. Mendelian randomization was performed to assess the causal relationship between these proteins and cardiovascular outcomes. RESULTS: Genome-wide association studies meta-analyses identified 4 newly associated loci: 3 with antithrombin levels (GCKR, BAZ1B, and HP-TXNL4B) and 1 with PS levels (ORM1-ORM2). transcriptome-wide association analyses identified 3 newly associated genes: 1 with antithrombin level (FCGRT), 1 with PC (GOLM2), and 1 with PS (MYL7). In addition, we replicated 7 independent loci reported in previous studies. Functional experiments provided evidence for the involvement of GCKR, SNX17, and HP genes in antithrombin regulation. CONCLUSIONS: The use of larger sample sizes, diverse populations, and a denser imputation reference panel allowed the detection of 7 novel genomic loci associated with plasma antithrombin, PC, and PS levels.

A rare coding mutation in the MAST2 gene causes venous thrombosis in a French family with unexplained thrombophilia: The Breizh MAST2 Arg89Gln variant.

Rare variants outside the classical coagulation cascade might cause inherited thrombosis. We aimed to identify the variant(s) causing venous thromboembolism (VTE) in a family with multiple relatives affected with unprovoked VTE and no thrombophilia defects. We identified by whole exome sequencing an extremely rare Arg to Gln variant (Arg89Gln) in the Microtubule Associated Serine/Threonine Kinase 2 (MAST2) gene that segregates with VTE in the family. Free-tissue factor pathway inhibitor (f-TFPI) plasma levels were significantly decreased in affected family members compared to healthy relatives. Conversely, plasminogen activator inhibitor-1 (PAI-1) levels were significantly higher in affected members than in healthy relatives. RNA sequencing analysis of RNA interference experimental data conducted in endothelial cells revealed that, of the 13,387 detected expressed genes, 2,354 have their level of expression modified by MAST2 knockdown, including SERPINE1 coding for PAI-1 and TFPI. In HEK293 cells overexpressing the MAST2 Gln89 variant, TFPI and SERPINE1 promoter activities were respectively lower and higher than in cells overexpressing the MAST2 wild type. This study identifies a novel thrombophilia-causing Arg89Gln variant in the MAST2 gene that is here proposed as a new molecular player in the etiology of VTE by interfering with hemostatic balance of endothelial cells.

Addressing unmeasured confounders in cohort studies: Instrumental variable method for a time-fixed exposure on an outcome trajectory.

Instrumental variable methods, which handle unmeasured confounding by targeting the part of the exposure explained by an exogenous variable not subject to confounding, have gained much interest in observational studies. We consider the very frequent setting of estimating the unconfounded effect of an exposure measured at baseline on the subsequent trajectory of an outcome repeatedly measured over time. We didactically explain how to apply the instrumental variable method in such setting by adapting the two-stage classical methodology with (1) the prediction of the exposure according to the instrumental variable, (2) its inclusion into a mixed model to quantify the exposure association with the subsequent outcome trajectory, and (3) the computation of the estimated total variance. A simulation study illustrates the consequences of unmeasured confounding in classical analyses and the usefulness of the instrumental variable approach. The methodology is then applied to 6224 participants of the 3C cohort to estimate the association of type-2 diabetes with subsequent cognitive trajectory, using 42 genetic polymorphisms as instrumental variables. This contribution shows how to handle endogeneity when interested in repeated outcomes, along with a R implementation. However, it should still be used with caution as it relies on instrumental variable assumptions hardly testable in practice.

Association of LIfestyle for BRAin health risk score (LIBRA) and genetic susceptibility with incident dementia and cognitive decline.

INTRODUCTION: Evaluating whether genetic susceptibility modifies the impact of lifestyle-related factors on dementia is critical for prevention. METHODS: We studied 5170 participants from a French cohort of older persons free of dementia at baseline and followed for up to 17 years. The LIfestyle for BRAin health risk score (LIBRA) including 12 modifiable factors was constructed at baseline (higher score indicating greater risk) and was related to both subsequent cognitive decline and dementia incidence, according to genetic susceptibility to dementia (reflected by the apolipoprotein E [APOE] ε4 allele and a genetic risk score [GRS]). RESULTS: The LIBRA was associated with higher dementia incidence, with no significant effect modification by genetics (hazard ratio for one point score = 1.09 [95% confidence interval, 1.05; 1.13]) in APOE ε4 non-carriers and = 1.15 [1.08; 1.22] in carriers; P = 0.15 for interaction). Similar findings were obtained with the GRS and with cognitive decline. DISCUSSION: Lifestyle-based prevention may be effective whatever the genetic susceptibility to dementia.

A Platelet Function Modulator of Thrombin Activation Is Causally Linked to Cardiovascular Disease and Affects PAR4 Receptor Signaling.

Dual antiplatelet therapy reduces ischemic events in cardiovascular disease, but it increases bleeding risk. Thrombin receptors PAR1 and PAR4 are drug targets, but the role of thrombin in platelet aggregation remains largely unexplored in large populations. We performed a genome-wide association study (GWAS) of platelet aggregation in response to full-length thrombin, followed by clinical association analyses, Mendelian randomization, and functional characterization including iPSC-derived megakaryocyte and platelet experiments. We identified a single sentinel variant in the GRK5 locus (rs10886430-G, p = 3.0 × 10-42) associated with increased thrombin-induced platelet aggregation (ß = 0.70, SE = 0.05). We show that disruption of platelet GRK5 expression by rs10886430-G is associated with enhanced platelet reactivity. The proposed mechanism of a GATA1-driven megakaryocyte enhancer is confirmed in allele-specific experiments. Utilizing further data, we demonstrate that the allelic effect is highly platelet- and thrombin-specific and not likely due to effects on thrombin levels. The variant is associated with increased risk of cardiovascular disease outcomes in UK BioBank, most strongly with pulmonary embolism. The variant associates with increased risk of stroke in the MEGASTROKE, UK BioBank, and FinnGen studies. Mendelian randomization analyses in independent samples support a causal role for rs10886430-G in increasing risk for stroke, pulmonary embolism, and venous thromboembolism through its effect on thrombin-induced platelet reactivity. We demonstrate that G protein-coupled receptor kinase 5 (GRK5) promotes platelet activation specifically via PAR4 receptor signaling. GRK5 inhibitors in development for the treatment of heart failure and cancer could have platelet off-target deleterious effects. Common variants in GRK5 may modify clinical outcomes with PAR4 inhibitors, and upregulation of GRK5 activity or signaling in platelets may have therapeutic benefits.

Association between ABO haplotypes and the risk of venous thrombosis: impact on disease risk estimation.

Genetic risk score (GRS) analysis is a popular approach to derive individual risk prediction models for complex diseases. In venous thrombosis (VT), such type of analysis shall integrate information at the ABO blood group locus, which is one of the major susceptibility loci. However, there is no consensus about which single nucleotide polymorphisms (SNPs) must be investigated when properly assessing association between ABO locus and VT risk. Using comprehensive haplotype analyses of ABO blood group tagging SNPs in 5425 cases and 8445 controls from 6 studies, we demonstrate that using only rs8176719 (tagging O1) to correctly assess the impact of ABO locus on VT risk is suboptimal, because 5% of rs8176719-delG carriers do not have an increased risk of developing VT. Instead, we recommend the use of 4 SNPs, rs2519093 (tagging A1), rs1053878 (A2), rs8176743 (B), and rs41302905 (O2), when assessing the impact of ABO locus on VT risk to avoid any risk misestimation. Compared with the O1 haplotype, the A2 haplotype is associated with a modest increase in VT risk (odds ratio, âˆ¼1.2), the A1 and B haplotypes are associated with an âˆ¼1.8-fold increased risk, whereas the O2 haplotype tends to be slightly protective (odds ratio, âˆ¼0.80). In addition, although the A1 and B blood groups are associated with increased von Willebrand factor and factor VIII plasma levels, only the A1 blood group is associated with ICAM levels, but in an opposite direction, leaving additional avenues to be explored to fully understand the spectrum of biological effects mediated by ABO locus on cardiovascular traits.

APOLD1 loss causes endothelial dysfunction involving cell junctions, cytoskeletal architecture, and Weibel-Palade bodies, while disrupting hemostasis.

Vascular homeostasis is impaired in various diseases thereby contributing to the progression of their underlying pathologies. The endothelial immediate early gene Apolipoprotein L domain-containing 1 (APOLD1) helps to regulate endothelial function. However, its precise role in endothelial cell biology remains unclear. We have localized APOLD1 to endothelial cell contacts and to Weibel-Palade bodies (WPB) where it associates with von Willebrand factor (VWF) tubules. Silencing of APOLD1 in primary human endothelial cells disrupted the cell junction-cytoskeletal interface, thereby altering endothelial permeability accompanied by spontaneous release of WPB contents. This resulted in an increased presence of WPB cargoes, notably VWF and angiopoietin-2 in the extracellular medium. Autophagy flux, previously recognized as an essential mechanism for the regulated release of WPB, was impaired in the absence of APOLD1. In addition, we report APOLD1 as a candidate gene for a novel inherited bleeding disorder across three generations of a large family in which an atypical bleeding diathesis was associated with episodic impaired microcirculation. A dominant heterozygous nonsense APOLD1:p.R49* variant segregated to affected family members. Compromised vascular integrity resulting from an excess of plasma angiopoietin-2, and locally impaired availability of VWF may explain the unusual clinical profile of APOLD1:p.R49* patients. In summary, our findings identify APOLD1 as an important regulator of vascular homeostasis and raise the need to consider testing of endothelial cell function in patients with inherited bleeding disorders without apparent platelet or coagulation defects.

Association of ABO blood groups with venous thrombosis recurrence in middle-aged patients: insights from a weighted Cox analysis dedicated to ambispective design.

BACKGROUND: In studies of time-to-events, it is common to collect information about events that occurred before the inclusion in a prospective cohort. When the studied risk factors are independent of time, including both pre- and post-inclusion events in the analyses, generally referred to as relying on an ambispective design, increases the statistical power but may lead to a selection bias. In the field of venous thromboembolism (VT), ABO blood groups have been the subject of extensive research due to their substantial effect on VT risk. However, few studies have investigated their effect on the risk of VT recurrence. Motivated by the study of the association of genetically determined ABO blood groups with VT recurrence, we propose a methodology to include pre-inclusion events in the analysis of ambispective studies while avoiding the selection bias due to mortality. METHODS: This work relies on two independent cohorts of VT patients, the French MARTHA study built on an ambispective design and the Dutch MEGA study built on a standard prospective design. For the analysis of the MARTHA study, a weighted Cox model was developed where weights were defined by the inverse of the survival probability at the time of data collection about the events. Thanks to the collection of information on the vital status of patients, we could estimate the survival probabilities using a delayed-entry Cox model on the death risk. Finally, results obtained in both studies were then meta-analysed. RESULTS: In the combined sample totalling 2,752 patients including 993 recurrences, the A1 blood group has an increased risk (Hazard Ratio (HR) of 1.18, p = 4.2 × 10-3) compared with the O1 group, homogeneously in MARTHA and in MEGA. The same trend (HR = 1.19, p = 0.06) was observed for the less frequent A2 group. CONCLUSION: The proposed methodology increases the power of studies relying on an ambispective design which is frequent in epidemiologic studies about recurrent events. This approach allowed to clarify the association of ABO blood groups with the risk of VT recurrence. Besides, this methodology has an immediate field of application in the context of genome wide association studies.

NPTX1 mutations trigger endoplasmic reticulum stress and cause autosomal dominant cerebellar ataxia.

With more than 40 causative genes identified so far, autosomal dominant cerebellar ataxias exhibit a remarkable genetic heterogeneity. Yet, half the patients are lacking a molecular diagnosis. In a large family with nine sampled affected members, we performed exome sequencing combined with whole-genome linkage analysis. We identified a missense variant in NPTX1, NM_002522.3:c.1165G>A: p.G389R, segregating with the phenotype. Further investigations with whole-exome sequencing and an amplicon-based panel identified four additional unrelated families segregating the same variant, for whom a common founder effect could be excluded. A second missense variant, NM_002522.3:c.980A>G: p.E327G, was identified in a fifth familial case. The NPTX1-associated phenotype consists of a late-onset, slowly progressive, cerebellar ataxia, with downbeat nystagmus, cognitive impairment reminiscent of cerebellar cognitive affective syndrome, myoclonic tremor and mild cerebellar vermian atrophy on brain imaging. NPTX1 encodes the neuronal pentraxin 1, a secreted protein with various cellular and synaptic functions. Both variants affect conserved amino acid residues and are extremely rare or absent from public databases. In COS7 cells, overexpression of both neuronal pentraxin 1 variants altered endoplasmic reticulum morphology and induced ATF6-mediated endoplasmic reticulum stress, associated with cytotoxicity. In addition, the p.E327G variant abolished neuronal pentraxin 1 secretion, as well as its capacity to form a high molecular weight complex with the wild-type protein. Co-immunoprecipitation experiments coupled with mass spectrometry analysis demonstrated abnormal interactions of this variant with the cytoskeleton. In agreement with these observations, in silico modelling of the neuronal pentraxin 1 complex evidenced a destabilizing effect for the p.E327G substitution, located at the interface between monomers. On the contrary, the p.G389 residue, located at the protein surface, had no predictable effect on the complex stability. Our results establish NPTX1 as a new causative gene in autosomal dominant cerebellar ataxias. We suggest that variants in NPTX1 can lead to cerebellar ataxia due to endoplasmic reticulum stress, mediated by ATF6, and associated to a destabilization of NP1 polymers in a dominant-negative manner for one of the variants.

Whole-exome sequencing identifies rare variants in STAB2 associated with venous thromboembolic disease.

Deep vein thrombosis and pulmonary embolism, collectively defined as venous thromboembolism (VTE), are the third leading cause of cardiovascular death in the United States. Common genetic variants conferring increased varying degrees of VTE risk have been identified by genome-wide association studies (GWAS). Rare mutations in the anticoagulant genes PROC, PROS1 and SERPINC1 result in perinatal lethal thrombosis in homozygotes and markedly increased VTE risk in heterozygotes. However, currently described VTE variants account for an insufficient portion of risk to be routinely used for clinical decision making. To identify new rare VTE risk variants, we performed a whole-exome study of 393 individuals with unprovoked VTE and 6114 controls. This study identified 4 genes harboring an excess number of rare damaging variants in patients with VTE: PROS1, STAB2, PROC, and SERPINC1. At STAB2, 7.8% of VTE cases and 2.4% of controls had a qualifying rare variant. In cell culture, VTE-associated variants of STAB2 had a reduced surface expression compared with reference STAB2. Common variants in STAB2 have been previously associated with plasma von Willebrand factor and coagulation factor VIII levels in GWAS, suggesting that haploinsufficiency of stabilin-2 may increase VTE risk through elevated levels of these procoagulants. In an independent cohort, we found higher von Willebrand factor levels and equivalent propeptide levels in individuals with rare STAB2 variants compared with controls. Taken together, this study demonstrates the utility of gene-based collapsing analyses to identify loci harboring an excess of rare variants with functional connections to a complex thrombotic disease.

Single-Cell Study of Two Rat Models of Pulmonary Arterial Hypertension Reveals Connections to Human Pathobiology and Drug Repositioning.

Rationale: The cellular and molecular landscape and translational value of commonly used models of pulmonary arterial hypertension (PAH) are poorly understood. Single-cell transcriptomics can enhance molecular understanding of preclinical models and facilitate their rational use and interpretation.Objectives: To determine and prioritize dysregulated genes, pathways, and cell types in lungs of PAH rat models to assess relevance to human PAH and identify drug repositioning candidates.Methods: Single-cell RNA sequencing was performed on the lungs of monocrotaline (MCT), Sugen-hypoxia (SuHx), and control rats to identify altered genes and cell types, followed by validation using flow-sorted cells, RNA in situ hybridization, and immunofluorescence. Relevance to human PAH was assessed by histology of lungs from patients and via integration with human PAH genetic loci and known disease genes. Candidate drugs were predicted using Connectivity Map.Measurements and Main Results: Distinct changes in genes and pathways in numerous cell types were identified in SuHx and MCT lungs. Widespread upregulation of NF-κB signaling and downregulation of IFN signaling was observed across cell types. SuHx nonclassical monocytes and MCT conventional dendritic cells showed particularly strong NF-κB pathway activation. Genes altered in SuHx nonclassical monocytes were significantly enriched for PAH-associated genes and genetic variants, and candidate drugs predicted to reverse the changes were identified. An open-access online platform was developed to share single-cell data and drug candidates (http://mergeomics.research.idre.ucla.edu/PVDSingleCell/).Conclusions: Our study revealed the distinct and shared dysregulation of genes and pathways in two commonly used PAH models for the first time at single-cell resolution and demonstrated their relevance to human PAH and utility for drug repositioning.

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.