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.

Mass Spectrometry-Based Pipeline for Identifying RNA Modifications Involved in a Functional Process: Application to Cancer Cell Adaptation.

Cancer onset and progression are known to be regulated by genetic and epigenetic events, including RNA modifications (a.k.a. epitranscriptomics). So far, more than 150 chemical modifications have been described in all RNA subtypes, including messenger, ribosomal, and transfer RNAs. RNA modifications and their regulators are known to be implicated in all steps of post-transcriptional regulation. The dysregulation of this complex yet delicate balance can contribute to disease evolution, particularly in the context of carcinogenesis, where cells are subjected to various stresses. We sought to discover RNA modifications involved in cancer cell adaptation to inhospitable environments, a peculiar feature of cancer stem cells (CSCs). We were particularly interested in the RNA marks that help the adaptation of cancer cells to suspension culture, which is often used as a surrogate to evaluate the tumorigenic potential. For this purpose, we designed an experimental pipeline consisting of four steps: (1) cell culture in different growth conditions to favor CSC survival; (2) simultaneous RNA subtype (mRNA, rRNA, tRNA) enrichment and RNA hydrolysis; (3) the multiplex analysis of nucleosides by LC-MS/MS followed by statistical/bioinformatic analysis; and (4) the functional validation of identified RNA marks. This study demonstrates that the RNA modification landscape evolves along with the cancer cell phenotype under growth constraints. Remarkably, we discovered a short epitranscriptomic signature, conserved across colorectal cancer cell lines and associated with enrichment in CSCs. Functional tests confirmed the importance of selected marks in the process of adaptation to suspension culture, confirming the validity of our approach and opening up interesting prospects in the field.

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.

Physical modeling of ribosomes along messenger RNA: Estimating kinetic parameters from ribosome profiling experiments using a ballistic model.

Gene expression is the synthesis of proteins from the information encoded on DNA. One of the two main steps of gene expression is the translation of messenger RNA (mRNA) into polypeptide sequences of amino acids. Here, by taking into account mRNA degradation, we model the motion of ribosomes along mRNA with a ballistic model where particles advance along a filament without excluded volume interactions. Unidirectional models of transport have previously been used to fit the average density of ribosomes obtained by the experimental ribo-sequencing (Ribo-seq) technique in order to obtain the kinetic rates. The degradation rate is not, however, accounted for and experimental data from different experiments are needed to have enough parameters for the fit. Here, we propose an entirely novel experimental setup and theoretical framework consisting in splitting the mRNAs into categories depending on the number of ribosomes from one to four. We solve analytically the ballistic model for a fixed number of ribosomes per mRNA, study the different regimes of degradation, and propose a criterion for the quality of the inverse fit. The proposed method provides a high sensitivity to the mRNA degradation rate. The additional equations coming from using the monosome (single ribosome) and polysome (arbitrary number) ribo-seq profiles enable us to determine all the kinetic rates in terms of the experimentally accessible mRNA degradation rate.

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.

NICOTIANAMINE SYNTHASE activity affects nucleolar iron accumulation and impacts rDNA silencing and RNA methylation in Arabidopsis.

In plant cells, a large pool of iron (Fe) is contained in the nucleolus, as well as in chloroplasts and mitochondria. A central determinant for intracellular distribution of Fe is nicotianamine (NA) generated by NICOTIANAMINE SYNTHASE (NAS). Here, we used Arabidopsis thaliana plants with disrupted NAS genes to study the accumulation of nucleolar iron and understand its role in nucleolar functions and more specifically in rRNA gene expression. We found that nas124 triple mutant plants, which contained lower quantities of the iron ligand NA, also contained less iron in the nucleolus. This was concurrent with the expression of normally silenced rRNA genes from nucleolar organizer regions 2 (NOR2). Notably, in nas234 triple mutant plants, which also contained lower quantities of NA, nucleolar iron and rDNA expression were not affected. In contrast, in both nas124 and nas234, specific RNA modifications were differentially regulated in a genotype dependent manner. Taken together, our results highlight the impact of specific NAS activities in RNA gene expression. We discuss the interplay between NA and nucleolar iron with rDNA functional organization and RNA methylation.

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.

Multiplexed LC-MS/MS quantification of salivary RNA modifications in periodontitis.

OBJECTIVE: To analyse the salivary epitranscriptomic profiles as periodontitis biomarkers using multiplexed mass spectrometry (MS). BACKGROUND: The field of epitranscriptomics, which relates to RNA chemical modifications, opens new perspectives in the discovery of diagnostic biomarkers, especially in periodontitis. Recently, the modified ribonucleoside N6-methyladenosine (m6A) was revealed as a crucial player in the etiopathogenesis of periodontitis. However, no epitranscriptomic biomarker has been identified in saliva to date. MATERIALS AND METHODS: Twenty-four saliva samples were collected from periodontitis patients (n = 16) and from control subjects (n = 8). Periodontitis patients were stratified according to stage and grade. Salivary nucleosides were directly extracted and, in parallel, salivary RNA was digested into its constituent nucleosides. Nucleoside samples were then quantified by multiplexed MS. RESULTS: Twenty-seven free nucleosides were detected and an overlapping set of 12 nucleotides were detected in digested RNA. Among the free nucleosides, cytidine and three other modified nucleosides (inosine, queuosine and m6Am) were significantly altered in periodontitis patients. In digested RNA, only uridine was significantly higher in periodontitis patients. Importantly there was no correlation between free salivary nucleoside levels and the levels of those same nucleotides in digested salivary RNA, except for cytidine, m5C and uridine. This statement implies that the two detection methods are complementary. CONCLUSION: The high specificity and sensitivity of MS allowed the detection and quantification of multiple nucleosides from RNA and free nucleosides in saliva. Some ribonucleosides appear to be promising biomarkers of periodontitis. Our analytic pipeline opens new perspectives for diagnostic periodontitis biomarkers.