Exome sequencing in undiagnosed congenital myopathy reveals new genes and refines genes-phenotypes correlations.

BACKGROUND: Congenital myopathies are severe genetic diseases with a strong impact on patient autonomy and often on survival. A large number of patients do not have a genetic diagnosis, precluding genetic counseling and appropriate clinical management. Our objective was to find novel pathogenic variants and genes associated with congenital myopathies and to decrease diagnostic odysseys and dead-end. METHODS: To identify pathogenic variants and genes implicated in congenital myopathies, we established and conducted the MYOCAPTURE project from 2009 to 2018 to perform exome sequencing in a large cohort of 310 families partially excluded for the main known genes. RESULTS: Pathogenic variants were identified in 156 families (50%), among which 123 families (40%) had a conclusive diagnosis. Only 44 (36%) of the resolved cases were linked to a known myopathy gene with the corresponding phenotype, while 55 (44%) were linked to pathogenic variants in a known myopathy gene with atypical signs, highlighting that most genetic diagnosis could not be anticipated based on clinical-histological assessments in this cohort. An important phenotypic and genetic heterogeneity was observed for the different genes and for the different congenital myopathy subtypes, respectively. In addition, we identified 14 new myopathy genes not previously associated with muscle diseases (20% of all diagnosed cases) that we previously reported in the literature, revealing novel pathomechanisms and potential therapeutic targets. CONCLUSIONS: Overall, this approach illustrates the importance of massive parallel gene sequencing as a comprehensive tool for establishing a molecular diagnosis for families with congenital myopathies. It also emphasizes the contribution of clinical data, histological findings on muscle biopsies, and the availability of DNA samples from additional family members to the diagnostic success rate. This study facilitated and accelerated the genetic diagnosis of congenital myopathies, improved health care for several patients, and opened novel perspectives for either repurposing of existing molecules or the development of novel treatments.

Constitutional DNA Polymorphisms Associated with the Plasma Imatinib Concentration in Chronic Myeloid Leukemia Patients.

The tyrosine kinase Inhibitor (TKI) imatinib is approved for the treatment of the chronic phase of chronic myeloid leukemia (CP-CML). Pharmacokinetic studies have highlighted the importance of inter-patient variability on imatinib plasma trough concentrations (ima[C]min). In the OPTIM-imatinib trial, we demonstrated that therapeutic drug monitoring (TDM) is able to improve the molecular response of CP-CML patients treated with imatinib. Here, we analyzed the constitutional exomes and RNAseq data of these patients. We performed an association analysis between the constitutional genetic variants of the patients and their ima[C]min, measured after 12 weeks of treatment with 400 mg once daily. Using linear regression, we identified 50 SNPs that showed excess heterozygosity depending on the ima[C]min. Ten SNPs were from non-coding sequences, and among the 40 remaining, 30 (from 25 genes) could be split into two categories. The first group of 16 SNPs concerns genes encoding extracellular matrix, cell junction, and membrane proteins. Coincidentally, cell adhesion proteins were also identified by RNA-seq as being overexpressed in patients with high ima[C]min. The other group of 14 SNPs were from genes encoding proteins involved in transcription/translation. Although most of the SNPs are intronic variants (28), we also identified missense (3), synonymous (4), 5'/3' (2), splicing (1), and upstream (4) variants. A haplotype analysis of four genes showed a significant association with high ima[C]min. None of the SNPs were significantly associated with the response. In conclusion, we identified a number of ima[C]min-associated SNPs, most of which correspond to genes encoding proteins that could play a role in the diffusion and transit of imatinib through membranes or epithelial barriers.

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.

Measuring the Efficiency of Purging by non-random Mating in Human Populations.

Human populations harbor a high concentration of deleterious genetic variants. Here, we tested the hypothesis that non-random mating practices affect the distribution of these variants, through exposure in the homozygous state, leading to their purging from the population gene pool. To do so, we produced whole-genome sequencing data for two pairs of Asian populations exhibiting different alliance rules and rates of inbreeding, but with similar effective population sizes. The results show that populations with higher rates of inbred matings do not purge deleterious variants more efficiently. Purging therefore has a low efficiency in human populations, and different mating practices lead to a similar mutational load.

Plasma levels of complement components C5 and C9 are associated with thrombin generation.

BACKGROUND: The thrombin generation assay (TGA) evaluates the potential of plasma to generate thrombin over time, providing a global picture of an individual's hemostatic balance. OBJECTIVES: This study aimed to identify novel biological determinants of thrombin generation using a multiomics approach. METHODS: Associations between TGA parameters and plasma levels of 377 antibodies targeting 236 candidate proteins for cardiovascular risk were tested using multiple linear regression analysis in 770 individuals with venous thrombosis from the Marseille Thrombosis Association (MARTHA) study. Proteins associated with at least 3 TGA parameters were selected for validation in an independent population of 536 healthy individuals (Etablissement Français du Sang Alpes-Méditerranée [EFS-AM]). Proteins with strongest associations in both groups underwent additional genetic analyses and in vitro experiments. RESULTS: Eighteen proteins were associated (P < 1.33 × 10⁻4) with at least 3 TGA parameters in MARTHA, among which 13 demonstrated a similar pattern of associations in EFS-AM. Complement proteins C5 and C9 had the strongest associations in both groups. Ex vivo supplementation of platelet-poor plasma with purified C9 protein had a significant dose-dependent effect on TGA parameters. No effect was observed with purified C5. Several single nucleotide polymorphisms associated with C5 and C9 plasma levels were identified, with the strongest association for the C5 missense variant rs17611, which was associated with a decrease in C5 levels, endogenous thrombin potential, and peak in MARTHA. No association of this variant with TGA parameters was observed in EFS-AM. CONCLUSION: This study identified complement proteins C5 and C9 as potential determinants of thrombin generation. Further studies are warranted to establish causality and elucidate the underlying mechanisms.

Genome-wide investigation of exogenous female hormones, genetic variation, and venous thromboembolism risk.

BACKGROUND: Increased risk of venous thromboembolism (VTE) is a life-threatening side effect for users of oral contraceptives (OCs) or hormone therapy (HT). OBJECTIVES: To investigate the potential for genetic predisposition to VTE in OC or HT users, we conducted a gene-by-environment case-only meta-analysis of genome-wide association studies (GWAS). METHODS: Use or nonuse of OCs (7 studies) or HT (8 studies) at the time of the VTE event was determined by pharmacy records or self-report. A synergy index (SI) was modeled for each variant in each study and submultiplicative/supramultiplicative gene-by-environment interactions were estimated. The SI parameters were first meta-analyzed across OC and HT studies and subsequently meta-analyzed to obtain an overall estimate. The primary analysis was agnostic GWAS and interrogated all imputed genotypes using a P value threshold of <5.0 × 10-8; secondary analyses were candidate-based. RESULTS: The VTE case-only OC meta-analysis included 2895 OC users and 6607 nonusers; the case-only HT meta-analysis included 2434 HT users and 12 793 nonusers. In primary GWAS meta-analyses, no variant reached genome-wide significance, but the smallest P value approached statistical significance: rs9386463 (P = 5.03 × 10-8). We tested associations for 138 candidate variants and identified 2 that exceeded statistical significance (0.05/138 = 3.62 × 10-4): F5 rs6025 (P = 1.87 × 10-5; SI, 1.29; previously observed) and F11 rs2036914 (P = 2.0 × 10-4; SI, 0.91; new observation). CONCLUSION: The candidate variant approach to identify submultiplictive/supramultiplicative associations between genetic variation and OC and HT use identified a new association with common genetic variation in F11, while the agnostic interrogations did not yield new discoveries.

MYLK*FLNB and DOCK1*LAMA2 gene-gene interactions associated with rheumatoid arthritis in the focal adhesion pathway.

Rheumatoid arthritis (RA) is a chronic, systemic autoimmune disease caused by a combination of genetic and environmental factors. Rare variants with low predicted effects in genes participating in the same biological function might be involved in developing complex diseases such as RA. From whole-exome sequencing (WES) data, we identified genes containing rare non-neutral variants with complete penetrance and no phenocopy in at least one of nine French multiplex families. Further enrichment analysis highlighted focal adhesion as the most significant pathway. We then tested if interactions between the genes participating in this function would increase or decrease the risk of developing RA disease. The model-based multifactor dimensionality reduction (MB-MDR) approach was used to detect epistasis in a discovery sample (19 RA cases and 11 healthy individuals from 9 families and 98 unrelated CEU controls from the International Genome Sample Resource). We identified 9 significant interactions involving 11 genes (MYLK, FLNB, DOCK1, LAMA2, RELN, PIP5K1C, TNC, PRKCA, VEGFB, ITGB5, and FLT1). One interaction (MYLK*FLNB) increasing RA risk and one interaction decreasing RA risk (DOCK1*LAMA2) were confirmed in a replication sample (200 unrelated RA cases and 91 GBR unrelated controls). Functional and genomic data in RA samples or relevant cell types argue the key role of these genes in RA.

Positive selection in the genomes of two Papua New Guinean populations at distinct altitude levels.

Highlanders and lowlanders of Papua New Guinea have faced distinct environmental stress, such as hypoxia and environment-specific pathogen exposure, respectively. In this study, we explored the top genomics regions and the candidate driver SNPs for selection in these two populations using newly sequenced whole-genomes of 54 highlanders and 74 lowlanders. We identified two candidate SNPs under selection - one in highlanders, associated with red blood cell traits and another in lowlanders, which is associated with white blood cell count - both potentially influencing the heart rate of Papua New Guineans in opposite directions. We also observed four candidate driver SNPs that exhibit linkage disequilibrium with an introgressed haplotype, highlighting the need to explore the possibility of adaptive introgression within these populations. This study reveals that the signatures of positive selection in highlanders and lowlanders of Papua New Guinea align closely with the challenges they face, which are specific to their environments.

A comparative genomics approach reveals a local genetic signature of Leishmania tropica in Morocco.

In Morocco, cutaneous leishmaniasis (CL) caused by Leishmania (L.) tropica is an important health problem. Despite the high incidence of CL in the country, the genomic heterogeneity of these parasites is still incompletely understood. In this study, we sequenced the genomes of 14 Moroccan isolates of L. tropica collected from confirmed cases of CL to investigate their genomic heterogeneity. Comparative genomics analyses were conducted by applying the recently established Genome Instability Pipeline (GIP), which allowed us to conduct phylogenomic and principal components analyses (PCA), and to assess genomic variations at the levels of the karyotype, gene copy number, single nucleotide polymorphisms (SNPs) and small insertions/deletions (INDELs) variants. Read-depth analyses revealed a mostly disomic karyotype, with the exception of the stable tetrasomy of chromosome 31. In contrast, we identified important gene copy number variations across all isolates, which affect known virulence genes and thus were probably selected in the field. SNP-based cluster analysis of the 14 isolates revealed a core group of 12 strains that formed a tight cluster and shared 45.1 % (87 751) of SNPs, as well as two strains (M3015, Ltr_16) that clustered separately from each other and the core group, suggesting the circulation of genetically highly diverse strains in Morocco. Phylogenetic analysis, which compared our 14 L. tropica isolates against 40 published genomes of L. tropica from a diverse array of locations, confirmed the genetic difference of our Moroccan isolates from all other isolates examined. In conclusion, our results indicate potential regional variations in SNP profiles that may differentiate Moroccan L. tropica from other L. tropica strains circulating in endemic countries in the Middle East. Our report paves the way for future research with a larger number of strains that will allow correlation of diverse phenotypes (resistance to treatments, virulence) and origins (geography, host species, year of isolation) to defined genomic signals such as gene copy number variations or SNP profiles that may represent interesting biomarker candidates.

TAD boundary deletion causes PITX2-related cardiac electrical and structural defects.

While 3D chromatin organization in topologically associating domains (TADs) and loops mediating regulatory element-promoter interactions is crucial for tissue-specific gene regulation, the extent of their involvement in human Mendelian disease is largely unknown. Here, we identify 7 families presenting a new cardiac entity associated with a heterozygous deletion of 2 CTCF binding sites on 4q25, inducing TAD fusion and chromatin conformation remodeling. The CTCF binding sites are located in a gene desert at 1 Mb from the Paired-like homeodomain transcription factor 2 gene (PITX2). By introducing the ortholog of the human deletion in the mouse genome, we recapitulate the patient phenotype and characterize an opposite dysregulation of PITX2 expression in the sinoatrial node (ectopic activation) and ventricle (reduction), respectively. Chromatin conformation assay performed in human induced pluripotent stem cell-derived cardiomyocytes harboring the minimal deletion identified in family#1 reveals a conformation remodeling and fusion of TADs. We conclude that TAD remodeling mediated by deletion of CTCF binding sites causes a new autosomal dominant Mendelian cardiac disorder.

Estimating SARS-CoV-2 infection probabilities with serological data and a Bayesian mixture model.

The individual results of SARS-CoV-2 serological tests measured after the first pandemic wave of 2020 cannot be directly interpreted as a probability of having been infected. Plus, these results are usually returned as a binary or ternary variable, relying on predefined cut-offs. We propose a Bayesian mixture model to estimate individual infection probabilities, based on 81,797 continuous anti-spike IgG tests from Euroimmun collected in France after the first wave. This approach used serological results as a continuous variable, and was therefore not based on diagnostic cut-offs. Cumulative incidence, which is necessary to compute infection probabilities, was estimated according to age and administrative region. In France, we found that a "negative" or a "positive" test, as classified by the manufacturer, could correspond to a probability of infection as high as 61.8% or as low as 67.7%, respectively. "Indeterminate" tests encompassed probabilities of infection ranging from 10.8 to 96.6%. Our model estimated tailored individual probabilities of SARS-CoV-2 infection based on age, region, and serological result. It can be applied in other contexts, if estimates of cumulative incidence are available.

Biallelic NAA60 variants with impaired n-terminal acetylation capacity cause autosomal recessive primary familial brain calcifications.

Primary familial brain calcification (PFBC) is characterized by calcium deposition in the brain, causing progressive movement disorders, psychiatric symptoms, and cognitive decline. PFBC is a heterogeneous disorder currently linked to variants in six different genes, but most patients remain genetically undiagnosed. Here, we identify biallelic NAA60 variants in ten individuals from seven families with autosomal recessive PFBC. The NAA60 variants lead to loss-of-function with lack of protein N-terminal (Nt)-acetylation activity. We show that the phosphate importer SLC20A2 is a substrate of NAA60 in vitro. In cells, loss of NAA60 caused reduced surface levels of SLC20A2 and a reduction in extracellular phosphate uptake. This study establishes NAA60 as a causal gene for PFBC, provides a possible biochemical explanation of its disease-causing mechanisms and underscores NAA60-mediated Nt-acetylation of transmembrane proteins as a fundamental process for healthy neurobiological functioning.

Upstream open reading frame-introducing variants in patients with primary familial brain calcification.

More than 50% of patients with primary familial brain calcification (PFBC), a rare neurological disorder, remain genetically unexplained. While some causative genes are yet to be identified, variants in non-coding regions of known genes may represent a source of missed diagnoses. We hypothesized that 5'-Untranslated Region (UTR) variants introducing an AUG codon may initiate mRNA translation and result in a loss of function in some of the PFBC genes. After reannotation of exome sequencing data of 113 unrelated PFBC probands, we identified two upstream AUG-introducing variants in the 5'UTR of PDGFB. One, NM_002608.4:c.-373C>G, segregated with PFBC in the family. It was predicted to create an upstream open reading frame (ORF). The other one, NM_002608.4:c.-318C>T, was found in a simplex case. It was predicted to result in an ORF overlapping the natural ORF with a frameshift. In a GFP reporter assay, both variants were associated with a dramatic decrease in GFP levels, and, after restoring the reading frame with the GFP sequence, the c.-318C>T variant was associated with a strong initiation of translation as measured by western blotting. Overall, we found upstream AUG-introducing variants in the 5'UTR of PDGFB in 2/113 (1.7%) undiagnosed PFBC cases. Such variants thus represent a source of putative pathogenic variants.

An elevated level of interleukin-17A in a Senegalese malaria cohort is associated with rs8193038 IL-17A genetic variant.

Malaria infection is a multifactorial disease partly modulated by host immuno-genetic factors. Recent evidence has demonstrated the importance of Interleukin-17 family proinflammatory cytokines and their genetic variants in host immunity. However, limited knowledge exists about their role in parasitic infections such as malaria. We aimed to investigate IL-17A serum levels in patients with severe and uncomplicated malaria and gene polymorphism's influence on the IL-17A serum levels. In this research, 125 severe (SM) and uncomplicated (UM) malaria patients and 48 free malaria controls were enrolled. IL-17A serum levels were measured with ELISA. PCR and DNA sequencing were used to assess host genetic polymorphisms in IL-17A. We performed a multivariate regression to estimate the impact of human IL-17A variants on IL-17A serum levels and malaria outcomes. Elevated serum IL-17A levels accompanied by increased parasitemia were found in SM patients compared to UM and controls (P < 0.0001). Also, the IL-17A levels were lower in SM patients who were deceased than in those who survived. In addition, the minor allele frequencies (MAF) of two IL-17A polymorphisms (rs3819024 and rs3748067) were more prevalent in SM patients than UM patients, indicating an essential role in SM. Interestingly, the heterozygous rs8193038 AG genotype was significantly associated with higher levels of IL-17A than the homozygous wild type (AA). According to our results, it can be concluded that the IL-17A gene rs8193038 polymorphism significantly affects IL-17A gene expression. Our results fill a gap in the implication of IL-17A gene polymorphisms on the cytokine level in a malaria cohort. IL-17A gene polymorphisms also may influence cytokine production in response to Plasmodium infections and may contribute to the hyperinflammatory responses during severe malaria outcomes.

SARS-CoV-2 infection prevalence and associated factors among primary healthcare workers in France after the third COVID-19 wave.

Data on the SARS-CoV-2 infection among primary health care workers (PHCWs) are scarce but essential to reflect on policy regarding prevention and control measures. We assessed the prevalence of PHCWs who have been infected by SARS-CoV-2 in comparison with modeling from the general population in metropolitan France, and associated factors. A cross-sectional study was conducted among general practitioners (GPs), pediatricians, dental and pharmacy workers in primary care between May and August 2021. Participants volunteered to provide a dried-blood spot for SARS-CoV-2 antibody assessment and completed a questionnaire. The primary outcome was defined as the detection of infection-induced antibodies (anti-nucleocapsid IgG, and for non-vaccinees: anti-Spike IgG and neutralizing antibodies) or previous self-reported infection (positive RT-qPCR or antigenic test, or positive ELISA test before vaccination). Estimates were adjusted using weights for representativeness and compared with prediction from the general population. Poisson regressions were used to quantify associated factors. The analysis included 1612 PHCWs. Weighted prevalences were: 31.7% (95% CI 27.5-36.0) for GPs, 28.7% (95% CI 24.4-33.0) for pediatricians, 25.2% (95% CI 20.6-31.0) for dentists, and 25.5% (95% CI 18.2-34.0) for pharmacists. Estimates were compatible with model predictions for the general population. PHCWs more likely to be infected were: GPs compared to pharmacist assistants (adjusted prevalence ratio [aPR] = 2.26; CI 95% 1.01-5.07), those living in Île-de-France (aPR = 1.53; CI 95% 1.14-2.05), South-East (aPR = 1.57; CI 95% 1.19-2.08), North-East (aPR = 1.81; CI 95% 1.38-2.37), and those having an unprotected contact with a COVID-19 case within the household (aPR = 1.48; CI 95% 1.22-1.80). Occupational factors were not associated with infection. In conclusion, the risk of SARS-CoV-2 exposure for PHCWs was more likely to have occurred in the community rather than at their workplace.

DLK1/DIO3 locus upregulation by a ß-catenin-dependent enhancer drives cell proliferation and liver tumorigenesis.

The CTNNB1 gene, encoding ß-catenin, is frequently mutated in hepatocellular carcinoma (HCC, ∼30%) and in hepatoblastoma (HB, >80%), in which DLK1/DIO3 locus induction is correlated with CTNNB1 mutations. Here, we aim to decipher how sustained ß-catenin activation regulates DLK1/DIO3 locus expression and the role this locus plays in HB and HCC development in mouse models deleted for Apc (ApcΔhep) or Ctnnb1-exon 3 (ß-cateninΔExon3) and in human CTNNB1-mutated hepatic cancer cells. We identified an enhancer site bound by TCF-4/ß-catenin complexes in an open conformation upon sustained ß-catenin activation (DLK1-Wnt responsive element [WRE]) and increasing DLK1/DIO3 locus transcription in ß-catenin-mutated human HB and mouse models. DLK1-WRE editing by CRISPR-Cas9 approach impaired DLK1/DIO3 locus expression and slowed tumor growth in subcutaneous CTNNB1-mutated tumor cell grafts, ApcΔhep HB and ß-cateninΔExon3 HCC. Tumor growth inhibition resulted either from increased FADD expression and subsequent caspase-3 cleavage in the first case or from decreased expression of cell cycle actors regulated by FoxM1 in the others. Therefore, the DLK1/DIO3 locus is an essential determinant of FoxM1-dependent cell proliferation during ß-catenin-driven liver tumorigenesis. Targeting the DLK1-WRE enhancer to silence the DLK1/DIO3 locus might thus represent an interesting therapeutic strategy to restrict tumor growth in primary liver cancers with CTNNB1 mutations.

De Novo Variants Found in Three Distinct Schizophrenia Populations Hit a Common Core Gene Network Related to Microtubule and Actin Cytoskeleton Gene Ontology Classes.

Schizophrenia (SZ) is a heterogeneous and debilitating psychiatric disorder with a strong genetic component. To elucidate functional networks perturbed in schizophrenia, we analysed a large dataset of whole-genome studies that identified SNVs, CNVs, and a multi-stage schizophrenia genome-wide association study. Our analysis identified three subclusters that are interrelated and with small overlaps: GO:0007017~Microtubule-Based Process, GO:00015629~Actin Cytoskeleton, and GO:0007268~SynapticTransmission. We next analysed three distinct trio cohorts of 75 SZ Algerian, 45 SZ French, and 61 SZ Japanese patients. We performed Illumina HiSeq whole-exome sequencing and identified de novo mutations using a Bayesian approach. We validated 88 de novo mutations by Sanger sequencing: 35 in French, 21 in Algerian, and 32 in Japanese SZ patients. These 88 de novo mutations exhibited an enrichment in genes encoding proteins related to GO:0051015~actin filament binding (p = 0.0011) using David, and enrichments in GO: 0003774~transport (p = 0.019) and GO:0003729~mRNA binding (p = 0.010) using Amigo. One of these de novo variant was found in CORO1C coding sequence. We studied Coro1c haploinsufficiency in a Coro1c+/- mouse and found defects in the corpus callosum. These results could motivate future studies of the mechanisms surrounding genes encoding proteins involved in transport and the cytoskeleton, with the goal of developing therapeutic intervention strategies for a subset of SZ cases.

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 GWAS meta-analyses of European- (71,771 cases and 1,059,740 controls) and African-ancestry samples (7,482 cases and 129,975 controls). We used LDpred2 and PRSCSx to construct ancestry-specific and multi-ancestry PRSs and evaluated their performance in an independent European- (6,261 cases and 88,238 controls) and African-ancestry sample (1,385 cases and 12,569 controls). Multi-ancestry PRSs with weights tuned in European- and African-ancestry samples, respectively, outperformed ancestry-specific PRSs in European- (PRSCSXEUR: AUC=0.61 (0.60, 0.61), PRSCSX_combinedEUR: AUC=0.61 (0.60, 0.62)) and African-ancestry test samples (PRSCSXAFR: AUC=0.58 (0.57, 0.6), PRSCSX_combined AFR: AUC=0.59 (0.57, 0.60)). 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 may be used to identify individuals at highest risk for VTE and provide guidance for the most effective treatment strategy across diverse populations.