Genome-wide meta-analyses of restless legs syndrome yield insights into genetic architecture, disease biology and risk prediction.

Restless legs syndrome (RLS) affects up to 10% of older adults. Their healthcare is impeded by delayed diagnosis and insufficient treatment. To advance disease prediction and find new entry points for therapy, we performed meta-analyses of genome-wide association studies in 116,647 individuals with RLS (cases) and 1,546,466 controls of European ancestry. The pooled analysis increased the number of risk loci eightfold to 164, including three on chromosome X. Sex-specific meta-analyses revealed largely overlapping genetic predispositions of the sexes (rg = 0.96). Locus annotation prioritized druggable genes such as glutamate receptors 1 and 4, and Mendelian randomization indicated RLS as a causal risk factor for diabetes. Machine learning approaches combining genetic and nongenetic information performed best in risk prediction (area under the curve (AUC) = 0.82-0.91). In summary, we identified targets for drug development and repurposing, prioritized potential causal relationships between RLS and relevant comorbidities and risk factors for follow-up and provided evidence that nonlinear interactions are likely relevant to RLS risk prediction.

Platelet Receptor Glycoprotein VI-Dimer Is Overexpressed in Patients with Atrial Fibrillation at High Risk of Ischemic Stroke.

Introduction Atrial fibrillation (AF) increases the risk of ischemic stroke (IS). We hypothesized that the functional form of platelet receptor glycoprotein (GP) VI, GPVI-dimer, which binds to collagen and fibrin causing platelet activation, is overexpressed in patients with AF who have not had a stroke. Methods A total of 75 inpatients with AF were recruited. None were admitted with or had previously had thrombotic events, including IS or myocardial infarction. Platelet surface expression of total GPVI, GPVI-dimer, and the platelet activation marker P-selectin were quantitated by whole blood flow cytometry. Serum biomarkers were collected in AF patients. Results were compared against patients contemporaneously admitted to hospital with similar age and vascular risk-factor profiles without AF (noAF, n = 30). Results Patients with AF have similar total GPVI surface expression ( p = 0.58) and P-selectin exposure ( p = 0.73) on their platelets compared with noAF patients but demonstrate significantly higher GPVI-dimer expression ( p = 0.02 ). Patients with paroxysmal AF express similar GPVI-dimer levels compared with permanent AF and GPVI-dimer levels were not different between anticoagulated groups. Serum N-terminal pro b-type natriuretic peptide ( p < 0.0001 ) and high sensitivity C-reactive protein ( p < 0.0001 ) were significantly correlated with GPVI-dimer expression in AF platelets. AF was the only vascular risk factor that was independently associated with higher GPVI-dimer expression in the whole population ( p = 0.02 ) . Conclusion GPVI inhibition is being explored in clinical trials as a novel target for IS treatment. As GPVI-dimer is elevated in AF patients' platelets, the exploration of targeted GPVI-dimer inhibition for stroke prevention in patients at high risk of IS due to AF is supported.

Detection and annotation of transposable element insertions and deletions on the human genome using nanopore sequencing.

Repetitive sequences represent about 45% of the human genome. Some are transposable elements (TEs) with the ability to change their position in the genome, creating genetic variability both as insertions or deletions, with potential pathogenic consequences. We used long-read nanopore sequencing to identify TE variants in the genomes of 24 patients with antithrombin deficiency. We identified 7 344 TE insertions and 3 056 TE deletions, 2 926 were not previously described in publicly available databases. The insertions affected 3 955 genes, with 6 insertions located in exons, 3 929 in introns, and 147 in promoters. Potential functional impact was evaluated with gene annotation and enrichment analysis, which suggested a strong relationship with neuron-related functions and autism. We conclude that this study encourages the generation of a complete map of TEs in the human genome, which will be useful for identifying new TEs involved in genetic disorders.

A genome-wide association study of blood cell morphology identifies cellular proteins implicated in disease aetiology.

Blood cells contain functionally important intracellular structures, such as granules, critical to immunity and thrombosis. Quantitative variation in these structures has not been subjected previously to large-scale genetic analysis. We perform genome-wide association studies of 63 flow-cytometry derived cellular phenotypes-including cell-type specific measures of granularity, nucleic acid content and reactivity-in 41,515 participants in the INTERVAL study. We identify 2172 distinct variant-trait associations, including associations near genes coding for proteins in organelles implicated in inflammatory and thrombotic diseases. By integrating with epigenetic data we show that many intracellular structures are likely to be determined in immature precursor cells. By integrating with proteomic data we identify the transcription factor FOG2 as an early regulator of platelet formation and α-granularity. Finally, we show that colocalisation of our associations with disease risk signals can suggest aetiological cell-types-variants in IL2RA and ITGA4 respectively mirror the known effects of daclizumab in multiple sclerosis and vedolizumab in inflammatory bowel disease.

The effects of pathogenic and likely pathogenic variants for inherited hemostasis disorders in 140 214 UK Biobank participants.

Rare genetic diseases affect millions, and identifying causal DNA variants is essential for patient care. Therefore, it is imperative to estimate the effect of each independent variant and improve their pathogenicity classification. Our study of 140 214 unrelated UK Biobank (UKB) participants found that each of them carries a median of 7 variants previously reported as pathogenic or likely pathogenic. We focused on 967 diagnostic-grade gene (DGG) variants for rare bleeding, thrombotic, and platelet disorders (BTPDs) observed in 12 367 UKB participants. By association analysis, for a subset of these variants, we estimated effect sizes for platelet count and volume, and odds ratios for bleeding and thrombosis. Variants causal of some autosomal recessive platelet disorders revealed phenotypic consequences in carriers. Loss-of-function variants in MPL, which cause chronic amegakaryocytic thrombocytopenia if biallelic, were unexpectedly associated with increased platelet counts in carriers. We also demonstrated that common variants identified by genome-wide association studies (GWAS) for platelet count or thrombosis risk may influence the penetrance of rare variants in BTPD DGGs on their associated hemostasis disorders. Network-propagation analysis applied to an interactome of 18 410 nodes and 571 917 edges showed that GWAS variants with large effect sizes are enriched in DGGs and their first-order interactors. Finally, we illustrate the modifying effect of polygenic scores for platelet count and thrombosis risk on disease severity in participants carrying rare variants in TUBB1 or PROC and PROS1, respectively. Our findings demonstrate the power of association analyses using large population datasets in improving pathogenicity classifications of rare variants.

A signature of platelet reactivity in CBC scattergrams reveals genetic predictors of thrombotic disease risk.

Genetic studies of platelet reactivity (PR) phenotypes may identify novel antiplatelet drug targets. However, such studies have been limited by small sample sizes (n < 5000) because of the complexity of measuring PR. We trained a model to predict PR from complete blood count (CBC) scattergrams. A genome-wide association study of this phenotype in 29 806 blood donors identified 21 distinct associations implicating 20 genes, of which 6 have been identified previously. The effect size estimates were significantly correlated with estimates from a study of flow cytometry-measured PR and a study of a phenotype of in vitro thrombus formation. A genetic score of PR built from the 21 variants was associated with the incidence rates of myocardial infarction and pulmonary embolism. Mendelian randomization analyses showed that PR was causally associated with the risks of coronary artery disease, stroke, and venous thromboembolism. Our approach provides a blueprint for using phenotype imputation to study the determinants of hard-to-measure but biologically important hematological traits.

Systematic Mendelian randomization using the human plasma proteome to discover potential therapeutic targets for stroke.

Stroke is the second leading cause of death with substantial unmet therapeutic needs. To identify potential stroke therapeutic targets, we estimate the causal effects of 308 plasma proteins on stroke outcomes in a two-sample Mendelian randomization framework and assess mediation effects by stroke risk factors. We find associations between genetically predicted plasma levels of six proteins and stroke (P ≤ 1.62 × 10-4). The genetic associations with stroke colocalize (Posterior Probability >0.7) with the genetic associations of four proteins (TFPI, TMPRSS5, CD6, CD40). Mendelian randomization supports atrial fibrillation, body mass index, smoking, blood pressure, white matter hyperintensities and type 2 diabetes as stroke risk factors (P ≤ 0.0071). Body mass index, white matter hyperintensity and atrial fibrillation appear to mediate the TFPI, IL6RA, TMPRSS5 associations with stroke. Furthermore, thirty-six proteins are associated with one or more of these risk factors using Mendelian randomization. Our results highlight causal pathways and potential therapeutic targets for stroke.

Long-Read Sequencing Identifies the First Retrotransposon Insertion and Resolves Structural Variants Causing Antithrombin Deficiency.

The identification of inherited antithrombin deficiency (ATD) is critical to prevent potentially life-threatening thrombotic events. Causal variants in SERPINC1 are identified for up to 70% of cases, the majority being single-nucleotide variants and indels. The detection and characterization of structural variants (SVs) in ATD remain challenging due to the high number of repetitive elements in SERPINC1. Here, we performed long-read whole-genome sequencing on 10 familial and 9 singleton cases with type I ATD proven by functional and antigen assays, who were selected from a cohort of 340 patients with this rare disorder because genetic analyses were either negative, ambiguous, or not fully characterized. We developed an analysis workflow to identify disease-associated SVs. This approach resolved, independently of its size or type, all eight SVs detected by multiple ligation-dependent probe amplification, and identified for the first time a complex rearrangement previously misclassified as a deletion. Remarkably, we identified the mechanism explaining ATD in 2 out of 11 cases with previous unknown defect: the insertion of a novel 2.4 kb SINE-VNTR-Alu retroelement, which was characterized by de novo assembly and verified by specific polymerase chain reaction amplification and sequencing in the probands and affected relatives. The nucleotide-level resolution achieved for all SVs allowed breakpoint analysis, which revealed repetitive elements and microhomologies supporting a common replication-based mechanism for all the SVs. Our study underscores the utility of long-read sequencing technology as a complementary method to identify, characterize, and unveil the molecular mechanism of disease-causing SVs involved in ATD, and enlarges the catalogue of genetic disorders caused by retrotransposon insertions.

Whole-exome sequencing identifies rare genetic variants associated with human plasma metabolites.

Metabolite levels measured in the human population are endophenotypes for biological processes. We combined sequencing data for 3,924 (whole-exome sequencing, WES, discovery) and 2,805 (whole-genome sequencing, WGS, replication) donors from a prospective cohort of blood donors in England. We used multiple approaches to select and aggregate rare genetic variants (minor allele frequency [MAF] < 0.1%) in protein-coding regions and tested their associations with 995 metabolites measured in plasma by using ultra-high-performance liquid chromatography-tandem mass spectrometry. We identified 40 novel associations implicating rare coding variants (27 genes and 38 metabolites), of which 28 (15 genes and 28 metabolites) were replicated. We developed algorithms to prioritize putative driver variants at each locus and used mediation and Mendelian randomization analyses to test directionality at associations of metabolite and protein levels at the ACY1 locus. Overall, 66% of reported associations implicate gene targets of approved drugs or bioactive drug-like compounds, contributing to drug targets' validating efforts.

Coagulation factor V is a T-cell inhibitor expressed by leukocytes in COVID-19.

Clotting Factor V (FV) is primarily synthesized in the liver and when cleaved by thrombin forms pro-coagulant Factor Va (FVa). Using whole blood RNAseq and scRNAseq of peripheral blood mononuclear cells, we find that FV mRNA is expressed in leukocytes, and identify neutrophils, monocytes, and T regulatory cells as sources of increased FV in hospitalized patients with COVID-19. Proteomic analysis confirms increased FV in circulating neutrophils in severe COVID-19, and immunofluorescence microscopy identifies FV in lung-infiltrating leukocytes in COVID-19 lung disease. Increased leukocyte FV expression in severe disease correlates with T-cell lymphopenia. Both plasma-derived and a cleavage resistant recombinant FV, but not thrombin cleaved FVa, suppress T-cell proliferation in vitro. Anticoagulants that reduce FV conversion to FVa, including heparin, may have the unintended consequence of suppressing the adaptive immune system.

Molecular Dissection of Structural Variations Involved in Antithrombin Deficiency.

Inherited antithrombin deficiency, the most severe form of thrombophilia, is predominantly caused by variants in SERPINC1. Few causal structural variants have been described, usually detected by multiplex ligation-dependent probe amplification or cytogenetic arrays, which only define the gain or loss and the approximate size and location. This study has done a complete dissection of the structural variants affecting SERPINC1 of 39 unrelated patients with antithrombin deficiency using multiplex ligation-dependent probe amplification, comparative genome hybridization array, long-range PCR, and whole genome nanopore sequencing. Structural variants, in all cases only affecting one allele, were deleterious and caused a severe type I deficiency. Most defects were deletions affecting exons of SERPINC1 (82.1%), but the whole cohort was heterogeneous, as tandem duplications, deletion of introns, or retrotransposon insertions were also detected. Their size was also variable, ranging from 193 bp to 8 Mb, and in 54% of the cases involved neighboring genes. All but two structural variants had repetitive elements and/or microhomologies in their breakpoints, suggesting a common mechanism of formation. This study also suggested regions recurrently involved in structural variants causing antithrombin deficiency and found three structural variants with a founder effect: the insertion of a retrotransposon, duplication of exon 6, and a 20-gene deletion. Finally, nanopore sequencing was determined to be the most appropriate method to identify and characterize all structural variants at nucleotide level, independently of their size or type.

Platelet surface receptor glycoprotein VI-dimer is overexpressed in stroke: The Glycoprotein VI in Stroke (GYPSIE) study results.

OBJECTIVES: Platelet activation underpins thrombus formation in ischemic stroke. The active, dimeric form of platelet receptor glycoprotein (GP) VI plays key roles by binding platelet ligands collagen and fibrin, leading to platelet activation. We investigated whether patients presenting with stroke expressed more GPVI on their platelet surface and had more active circulating platelets as measured by platelet P-selectin exposure. METHODS: 129 ischemic or hemorrhagic stroke patients were recruited within 8h of symptom onset. Whole blood was analyzed for platelet-surface expression of total GPVI, GPVI-dimer, and P-selectin by flow cytometry at admission and day-90 post-stroke. Results were compared against a healthy control population (n = 301). RESULTS: The platelets of stroke patients expressed significantly higher total GPVI and GPVI-dimer (P<0.0001) as well as demonstrating higher resting P-selectin exposure (P<0.0001), a measure of platelet activity, compared to the control group, suggesting increased circulating platelet activation. GPVI-dimer expression was strongly correlated circulating platelet activation [r2 = 0.88, P<0.0001] in stroke patients. Furthermore, higher platelet surface GPVI expression was associated with increased stroke severity at admission. At day-90 post-stroke, GPVI-dimer expression and was further raised compared to the level at admission (P<0.0001) despite anti-thrombotic therapy. All ischemic stroke subtypes and hemorrhagic strokes expressed significantly higher GPVI-dimer compared to controls (P<0.0001). CONCLUSIONS: Stroke patients express more GPVI-dimer on their platelet surface at presentation, lasting at least until day-90 post-stroke. Small molecule GPVI-dimer inhibitors are currently in development and the results of this study validate that GPVI-dimer as an anti-thrombotic target in ischemic stroke.

Machine learning optimized polygenic scores for blood cell traits identify sex-specific trajectories and genetic correlations with disease.

Genetic association studies for blood cell traits, which are key indicators of health and immune function, have identified several hundred associations and defined a complex polygenic architecture. Polygenic scores (PGSs) for blood cell traits have potential clinical utility in disease risk prediction and prevention, but designing PGS remains challenging and the optimal methods are unclear. To address this, we evaluated the relative performance of 6 methods to develop PGS for 26 blood cell traits, including a standard method of pruning and thresholding (P + T) and 5 learning methods: LDpred2, elastic net (EN), Bayesian ridge (BR), multilayer perceptron (MLP) and convolutional neural network (CNN). We evaluated these optimized PGSs on blood cell trait data from UK Biobank and INTERVAL. We find that PGSs designed using common machine learning methods EN and BR show improved prediction of blood cell traits and consistently outperform other methods. Our analyses suggest EN/BR as the top choices for PGS construction, showing improved performance for 25 blood cell traits in the external validation, with correlations with the directly measured traits increasing by 10%-23%. Ten PGSs showed significant statistical interaction with sex, and sex-specific PGS stratification showed that all of them had substantial variation in the trajectories of blood cell traits with age. Genetic correlations between the PGSs for blood cell traits and common human diseases identified well-known as well as new associations. We develop machine learning-optimized PGS for blood cell traits, demonstrate their relationships with sex, age, and disease, and make these publicly available as a resource.

G protein-coupled receptor kinase 5 regulates thrombin signaling in platelets via PAR-1.

The interindividual variation in the functional response of platelets to activation by agonists is heritable. Genome-wide association studies (GWASs) of quantitative measures of platelet function have identified fewer than 20 distinctly associated variants, some with unknown mechanisms. Here, we report GWASs of pathway-specific functional responses to agonism by adenosine 5'-diphosphate, a glycoprotein VI-specific collagen mimetic, and thrombin receptor-agonist peptides, each specific to 1 of the G protein-coupled receptors PAR-1 and PAR-4, in subsets of 1562 individuals. We identified an association (P = 2.75 × 10-40) between a common intronic variant, rs10886430, in the G protein-coupled receptor kinase 5 gene (GRK5) and the sensitivity of platelets to activate through PAR-1. The variant resides in a megakaryocyte-specific enhancer that is bound by the transcription factors GATA1 and MEIS1. The minor allele (G) is associated with fewer GRK5 transcripts in platelets and the greater sensitivity of platelets to activate through PAR-1. We show that thrombin-mediated activation of human platelets causes binding of GRK5 to PAR-1 and that deletion of the mouse homolog Grk5 enhances thrombin-induced platelet activation sensitivity and increases platelet accumulation at the site of vascular injury. This corroborates evidence that the human G allele of rs10886430 is associated with a greater risk for cardiovascular disease. In summary, by combining the results of pathway-specific GWASs and expression quantitative trait locus studies in humans with the results from platelet function studies in Grk5-/- mice, we obtain evidence that GRK5 regulates the human platelet response to thrombin via the PAR-1 pathway.

Circular RNAs exhibit limited evidence for translation, or translation regulation of the mRNA counterpart in terminal hematopoiesis.

Each day, about 1012 erythrocytes and platelets are released into the bloodstream. This substantial output from hematopoietic stem cells is tightly regulated by transcriptional and epigenetic factors. Whether and how circular RNAs (circRNAs) contribute to the differentiation and/or identity of hematopoietic cells is to date not known. We recently reported that erythrocytes and platelets contain the highest levels and numbers of circRNAs among hematopoietic cells. Here, we provide the first detailed analysis of circRNA expression during erythroid and megakaryoid differentiation. CircRNA expression not only significantly increased upon enucleation, but also had limited overlap between progenitor cells and mature cells, suggesting that circRNA expression stems from regulated processes rather than resulting from mere accumulation. To study circRNA function in hematopoiesis, we first compared the expression levels of circRNAs with the translation efficiency of their mRNA counterpart. We found that only one out of 2531 (0.04%) circRNAs associated with mRNA-translation regulation. Furthermore, irrespective of thousands of identified putative open reading frames, deep ribosome-footprinting sequencing, and mass spectrometry analysis provided little evidence for translation of endogenously expressed circRNAs. In conclusion, circRNAs alter their expression profile during terminal hematopoietic differentiation, yet their contribution to regulate cellular processes remains enigmatic.

Integrative analysis of the plasma proteome and polygenic risk of cardiometabolic diseases.

Cardiometabolic diseases are frequently polygenic in architecture, comprising a large number of risk alleles with small effects spread across the genome1-3. Polygenic scores (PGS) aggregate these into a metric representing an individual's genetic predisposition to disease. PGS have shown promise for early risk prediction4-7 and there is an open question as to whether PGS can also be used to understand disease biology8. Here, we demonstrate that cardiometabolic disease PGS can be used to elucidate the proteins underlying disease pathogenesis. In 3,087 healthy individuals, we found that PGS for coronary artery disease, type 2 diabetes, chronic kidney disease and ischaemic stroke are associated with the levels of 49 plasma proteins. Associations were polygenic in architecture, largely independent of cis and trans protein quantitative trait loci and present for proteins without quantitative trait loci. Over a follow-up of 7.7 years, 28 of these proteins associated with future myocardial infarction or type 2 diabetes events, 16 of which were mediators between polygenic risk and incident disease. Twelve of these were druggable targets with therapeutic potential. Our results demonstrate the potential for PGS to uncover causal disease biology and targets with therapeutic potential, including those that may be missed by approaches utilizing information at a single locus.

Polygenic basis and biomedical consequences of telomere length variation.

Telomeres, the end fragments of chromosomes, play key roles in cellular proliferation and senescence. Here we characterize the genetic architecture of naturally occurring variation in leukocyte telomere length (LTL) and identify causal links between LTL and biomedical phenotypes in 472,174 well-characterized UK Biobank participants. We identified 197 independent sentinel variants associated with LTL at 138 genomic loci (108 new). Genetically determined differences in LTL were associated with multiple biological traits, ranging from height to bone marrow function, as well as several diseases spanning neoplastic, vascular and inflammatory pathologies. Finally, we estimated that, at the age of 40 years, people with an LTL >1 s.d. shorter than the population mean had a 2.5-year-lower life expectancy compared with the group with ≥1 s.d. longer LDL. Overall, we furnish new insights into the genetic regulation of LTL, reveal wide-ranging influences of LTL on physiological traits, diseases and longevity, and provide a powerful resource available to the global research community.

GoldVariants, a resource for sharing rare genetic variants detected in bleeding, thrombotic, and platelet disorders: Communication from the ISTH SSC Subcommittee on Genomics in Thrombosis and Hemostasis.

The implementation of high-throughput sequencing (HTS) technologies in research and diagnostic laboratories has linked many new genes to rare bleeding, thrombotic, and platelet disorders (BTPD), and revealed multiple genetic variants linked to those disorders, many of them being of uncertain pathogenicity when considering the accepted evidence (variant consequence, frequency in control datasets, number of reported patients, prediction models, and functional assays). The sequencing effort has also resulted in resources for gathering disease-causing variants associated with specific genes, but for BTPD, such well-curated databases exist only for a few genes. On the other hand, submissions by individuals or diagnostic laboratories to the variant database ClinVar are hampered by the lack of a submission process tailored to capture the specific features of hemostatic diseases. As we move toward the implementation of HTS in the diagnosis of BTPD, the Scientific and Standardization Committee for Genetics in Thrombosis and Haemostasis has developed and tested a REDCap-based interface, aimed at the community, to submit curated genetic variants for diagnostic-grade BTPD genes. Here, we describe the use of the interface and the initial submission of 821 variants from 30 different centers covering 14 countries. This open-access variant resource will be shared with the community to improve variant classification and regular bulk data transfer to ClinVar.