Biologic and clinical analysis of childhood gamma delta T-ALL identifies LMO2/STAG2 rearrangements as extremely high-risk.

Acute lymphoblastic leukemia expressing the gamma delta T cell receptor (yo T-ALL) is a poorly understood disease. We studied 200 children with yo T-ALL from 13 clinical study groups to understand the clinical and genetic features of this disease. We found age and genetic drivers were significantly associated with outcome. yo T-ALL diagnosed in children under three years of age was extremely high-risk and enriched for genetic alterations that result in both LMO2 activation and STAG2 inactivation. Mechanistically, using patient samples and isogenic cell lines, we show that inactivation of STAG2 profoundly perturbs chromatin organization by altering enhancer-promoter looping, resulting in deregulation of gene expression associated with T-cell differentiation. High throughput drug screening identified a vulnerability in DNA repair pathways arising from STAG2 inactivation, which can be targeted by Poly(ADP-ribose) polymerase (PARP) inhibition. These data provide a diagnostic framework for classification and risk stratification of pediatric yo T-ALL.

Regulatory variants of APOBEC3 genes potentially associate with COVID-19 severity in populations with African ancestry.

Since November 2019, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has caused the worldwide pandemic of the coronavirus disease 2019 (COVID-19), the impact of which is huge to the lives of world populations. Many studies suggested that such situation will continue due to the endless mutations in SARS-CoV-2 genome that result in complexity of the efforts for the control of SARS-CoV-2, since the special enrichment of nucleotide substitution C>U in SARS-CoV-2 sequences were discovered mainly due to the editing by human host factors APOBEC3 genes. The observation of SARS-CoV-2 variants Beta (B.1.351) and Omicron (B.1.1.529) firstly spreading in South Africa promoted us to hypothesize that genetic variants of APOBEC3 special in African populations may be attributed to the higher mutation rate of SARS-CoV-2 variants in Africa. Current study was conducted to search for functional variants of APOBEC3 genes associate with COVID-19 hospitalization in African population. By integrating data from the 1000 Genomes Project, Genotype-Tissue Expression (GTEx), and Host Genetics Initiative (HGI) of COVID-19, we identified potential functional SNPs close to APOBEC3 genes that are associated with COVID-19 hospitalization in African but not with other populations. Our study provides new insights on the potential contribution of APOBEC3 genes on the evolution of SARS-CoV-2 mutations in African population, but further replication is needed to confirm our results.

MAP3K19 regulatory variation in populations with African ancestry may increase COVID-19 severity.

To identify ancestry-linked genetic risk variants associated with COVID-19 hospitalization, we performed an integrative analysis of two genome-wide association studies and resolved four single nucleotide polymorphisms more frequent in COVID-19-hospitalized patients with non-European ancestry. Among them, the COVID-19 risk SNP rs16831827 shows the largest difference in minor allele frequency (MAF) between populations with African and European ancestry and also shows higher MAF in hospitalized COVID-19 patients among cohorts of mixed ancestry (odds ratio [OR] = 1.20, 95% CI: 1.10-1.30) and entirely African ancestry (OR = 1.30, 95% CI: 1.02-1.67). rs16831827 is an expression quantitative trait locus of MAP3K19. MAP3K19 expression is induced during ciliogenesis and most abundant in ciliated tissues including lungs. Single-cell RNA sequencing analyses revealed that MAP3K19 is highly expressed in multiple ciliated cell types. As rs16831827∗T is associated with reduced MAP3K19 expression, it may increase the risk of severe COVID-19 by reducing MAP3K19 expression.

Genomic landscape of Down syndrome-associated acute lymphoblastic leukemia.

Trisomy 21, the genetic cause of Down syndrome (DS), is the most common congenital chromosomal anomaly. It is associated with a 20-fold increased risk of acute lymphoblastic leukemia (ALL) during childhood and results in distinctive leukemia biology. To comprehensively define the genomic landscape of DS-ALL, we performed whole-genome sequencing and whole-transcriptome sequencing (RNA-Seq) on 295 cases. Our integrated genomic analyses identified 15 molecular subtypes of DS-ALL, with marked enrichment of CRLF2-r, IGH::IGF2BP1, and C/EBP altered (C/EBPalt) subtypes compared with 2257 non-DS-ALL cases. We observed abnormal activation of the CEBPD, CEBPA, and CEBPE genes in 10.5% of DS-ALL cases via a variety of genomic mechanisms, including chromosomal rearrangements and noncoding mutations leading to enhancer hijacking. A total of 42.3% of C/EBP-activated DS-ALL also have concomitant FLT3 point mutations or insertions/deletions, compared with 4.1% in other subtypes. CEBPD overexpression enhanced the differentiation of mouse hematopoietic progenitor cells into pro-B cells in vitro, particularly in a DS genetic background. Notably, recombination-activating gene-mediated somatic genomic abnormalities were common in DS-ALL, accounting for a median of 27.5% of structural alterations, compared with 7.7% in non-DS-ALL. Unsupervised hierarchical clustering analyses of CRLF2-rearranged DS-ALL identified substantial heterogeneity within this group, with the BCR::ABL1-like subset linked to an inferior event-free survival, even after adjusting for known clinical risk factors. These results provide important insights into the biology of DS-ALL and point to opportunities for targeted therapy and treatment individualization.

Genome-wide screening of sex-biased genetic variants potentially associated with COVID-19 hospitalization.

Sex-biased difference in coronavirus disease 2019 (COVID-19) hospitalization has been observed as that male patients tend to be more likely to be hospitalized than female patients. However, due to the insufficient sample size and existed studies that more prioritized to sex-stratified COVID-19 genome-wide association study (GWAS), the searching for sex-biased genetic variants showing differential association signals between sexes with COVID-19 hospitalization was severely hindered. We hypothesized genetic variants would show potentially sex-biased genetic effects on COVID-19 hospitalization if they display significant differential association effect sizes between male and female COVID-19 patients. By integrating two COVID-19 GWASs, including hospitalized COVID-19 patients vs. general population separated into males (case = 1,917 and control = 221,174) and females (case = 1,343 and control = 262,886), we differentiated the association effect sizes of each common single nucleotide polymorphism (SNP) within the two GWASs. Twelve SNPs were suggested to show differential COVID-19 associations between sexes. Further investigation of genes (n = 58) close to these 12 SNPs resulted in the identification of 34 genes demonstrating sex-biased differential expression in at least one GTEx tissue. Finally, 5 SNPs are mapped to 8 genes, including rs1134004 (GADD45G), rs140657166 (TRIM29 and PVRL1), rs148143613 (KNDC1 and STK32C), rs2443615 (PGAP2 and TRIM21), and rs2924725 (CSMD1). The 8 genes display significantly differential gene expression in blood samples derived from COVID-19 patients compared to healthy controls. These genes are potential genetic factors contributing to sex differences in COVID-19 hospitalization and warranted for further functional studies.

Evaluation of a Functional Single Nucleotide Polymorphism of the SARS-CoV-2 Receptor ACE2 That Is Potentially Involved in Long COVID.

Since the occurrence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in December 2019, SARS-CoV-2 has led to a global coronavirus disease 2019 (COVID-19) pandemic. A better understanding of the SARS-CoV-2 receptor ACE2 at the genetic level would help combat COVID-19, particularly for long COVID. We performed a genetic analysis of ACE2 and searched for its common potential single nucleotide polymorphisms (SNPs) with minor allele frequency >0.05 in both European and Chinese populations that would contribute to ACE2 gene expression variation. We thought that the variation of the ACE2 expression would be an important biological feature that would strongly affect COVID-19 symptoms, such as "brain fog", which is highlighted by the fact that ACE2 acts as a major cellular receptor for SARS-CoV-2 attachment and is highly expressed in brain tissues. Based on the human GTEx gene expression database, we found rs2106809 exhibited a significant correlation with the ACE2 expression among multiple brain and artery tissues. This expression correlation was replicated in an independent European brain eQTL database, Braineac. rs2106809*G also displays significantly higher frequency in Asian populations than in Europeans and displays a protective effect (p = 0.047) against COVID-19 hospitalization when comparing hospitalized COVID-19 cases with non-hospitalized COVID-19 or SARS-CoV-2 test-negative samples with European ancestry from the UK Biobank. Furthermore, we experimentally demonstrated that rs2106809*G could upregulate the transcriptional activity of ACE2. Therefore, integrative analysis and functional experiment strongly support that ACE2 SNP rs2106809 is a functional brain eQTL and its potential involvement in long COVID, which warrants further investigation.

Somatic and Germline Variant Calling from Next-Generation Sequencing Data.

Re-sequencing of the human genome by next-generation sequencing (NGS) has been widely applied to discover pathogenic genetic variants and/or causative genes accounting for various types of diseases including cancers. The advances in NGS have allowed the sequencing of the entire genome of patients and identification of disease-associated variants in a reasonable timeframe and cost. The core of the variant identification relies on accurate variant calling and annotation. Numerous algorithms have been developed to elucidate the repertoire of somatic and germline variants. Each algorithm has its own distinct strengths, weaknesses, and limitations due to the difference in the statistical modeling approach adopted and read information utilized. Accurate variant calling remains challenging due to the presence of sequencing artifacts and read misalignments. All of these can lead to the discordance of the variant calling results and even misinterpretation of the discovery. For somatic variant detection, multiple factors including chromosomal abnormalities, tumor heterogeneity, tumor-normal cross contaminations, unbalanced tumor/normal sample coverage, and variants with low allele frequencies add even more layers of complexity to accurate variant identification. Given the discordances and difficulties, ensemble approaches have emerged by harmonizing information from different algorithms to improve variant calling performance. In this chapter, we first introduce the general scheme of variant calling algorithms and potential challenges at distinct stages. We next review the existing workflows of variant calling and annotation, and finally explore the strategies deployed by different callers as well as their strengths and caveats. Overall, NGS-based variant identification with careful consideration allows reliable detection of pathogenic variant and candidate variant selection for precision medicine.

Epigenome-wide association study of posttraumatic stress disorder identifies novel loci in U.S. military veterans.

Posttraumatic stress disorder (PTSD) is a chronic and disabling psychiatric disorder prevalent in military veterans. Epigenetic mechanisms have been implicated in the etiology of PTSD, with DNA methylation being the most studied to identify novel molecular biomarkers associated with this disorder. We performed one of the largest single-sample epigenome-wide association studies (EWAS) of PTSD to date. Our sample included 1135 male European-American U.S. veterans who participated in the National Health and Resilience in Veterans Study (NHRVS). DNA was collected from saliva samples and the Illumina HumanMethylation EPIC BeadChip was used for the methylation analysis. PTSD was assessed using the PTSD Checklist. An EWAS was conducted using linear regression adjusted for age, cell-type proportions, first 10 principal components, and smoking status. After Bonferroni correction, we identified six genome-wide significant (GWS) CpG sites associated with past-month PTSD and three CpGs with lifetime PTSD (prange = 10-10-10-8). These CpG sites map to genes involved in immune function, transcription regulation, axonal guidance, cell signaling, and protein binding. Among these, SENP7, which is involved in transcription regulation and has been linked to risk-taking behavior and alcohol consumption in genome-wide association studies, replicated in an independent veteran cohort and was downregulated in medial orbitofrontal cortex of PTSD postmortem brain tissue. These findings suggest potential epigenetic biomarkers of PTSD that may help inform the pathophysiology of this disorder in veterans and other trauma-affected populations.

SPINK6 inhibits human airway serine proteases and restricts influenza virus activation.

SPINK6 was identified in human skin as a cellular inhibitor of serine proteases of the KLK family. Airway serine proteases are required to cleave hemagglutinin (HA) of influenza A viruses (IAVs) to initiate an infection in the human airway. We hypothesized that SPINK6 may inhibit common airway serine proteases and restrict IAV activation. We demonstrate that SPINK6 specifically suppresses the proteolytic activity of HAT and KLK5, HAT- and KLK5-mediated HA cleavage, and restricts virus maturation and replication. SPINK6 constrains the activation of progeny virions and impairs viral growth; and vice versa, blocking endogenous SPINK6 enhances HA cleavage and viral growth in physiological-relevant human airway organoids where SPINK6 is intrinsically expressed. In IAV-infected mice, SPINK6 significantly suppresses viral growth and improves mouse survival. Notably, individuals carrying the higher SPINK6 expression allele were protected from human H7N9 infection. Collectively, SPINK6 is a novel host inhibitor of serine proteases in the human airway and restricts IAV activation.

Genetic screening for hypertension and COVID-19 reveals functional variation of SPEG potentially associated with severe COVID-19 in women.

The coronavirus disease 2019 (COVID-19) pandemic, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has led to more than 6.4 million deaths worldwide. The prevalent comorbidity between hypertension and severe COVID-19 suggests common genetic factors may affect the outcome of both diseases. As both hypertension and severe COVID-19 demonstrate sex-biased prevalence, common genetic factors between the two diseases may display sex-biased differential associations. By evaluating COVID-19 association signals of 172-candidate hypertension single nucleotide polymorphisms (SNPs) derived from more than 1 million European individuals in two sex-stratified severe COVID-19 genome-wide association studies from UK BioBank with European ancestry, we revealed one functional cis expression quantitative trait locus of SPEG (rs12474050) showing sex-biased association with severe COVID-19 in women. The risk allele rs12474050*T associates with higher blood pressure. In our study, we found it is significantly correlated with lower SPEG expression in muscle-skeletal but with higher expression in both brain cerebellum and cerebellar hemisphere. Additionally, nominal significances were detected for the association between rs12474050*T and lower SPEG expression in both heart left ventricle and atrial appendage; among these tissues, the SPEG expression is nominally significantly higher in females than in males. Further analysis revealed SPEG is mainly expressed in cardiomyocytes in heart and is upregulated upon SARS-CoV-2 infection, with significantly higher upregulation of SPEG only observed in female but not in male COVID-19 patients compared to both normal female and male individuals, suggesting upregulation of SPEG is a female-specific protective mechanism against COVID-19 induced heart damage. Taken together, our analyses suggest the involvement of SPEG in both hypertension and severe COVID-19 in women, which provides new insights for sex-biased effect of severe COVID-19 in women.

Identifying tumorigenic non-coding mutations through altered cis-regulation.

Identification of non-coding mutations driving tumorigenesis requires alternative approaches to coding mutations. Enriched associations between mutated regulatory elements and altered cis-regulation in tumors are a promising approach to stratify candidate non-coding driver mutations. Here we provide a bioinformatics pipeline to mine data from the Cancer Genomic Commons (GDC) for such associations. The pipeline integrates RNA and whole-genome sequencing with genotyping data to reveal putative non-coding driver mutations by cancer type. For complete information on the generation and use of this protocol, please refer to Cheng et al. (2021).

Enhancer Hijacking Drives Oncogenic BCL11B Expression in Lineage-Ambiguous Stem Cell Leukemia.

Lineage-ambiguous leukemias are high-risk malignancies of poorly understood genetic basis. Here, we describe a distinct subgroup of acute leukemia with expression of myeloid, T lymphoid, and stem cell markers driven by aberrant allele-specific deregulation of BCL11B, a master transcription factor responsible for thymic T-lineage commitment and specification. Mechanistically, this deregulation was driven by chromosomal rearrangements that juxtapose BCL11B to superenhancers active in hematopoietic progenitors, or focal amplifications that generate a superenhancer from a noncoding element distal to BCL11B. Chromatin conformation analyses demonstrated long-range interactions of rearranged enhancers with the expressed BCL11B allele and association of BCL11B with activated hematopoietic progenitor cell cis-regulatory elements, suggesting BCL11B is aberrantly co-opted into a gene regulatory network that drives transformation by maintaining a progenitor state. These data support a role for ectopic BCL11B expression in primitive hematopoietic cells mediated by enhancer hijacking as an oncogenic driver of human lineage-ambiguous leukemia. SIGNIFICANCE: Lineage-ambiguous leukemias pose significant diagnostic and therapeutic challenges due to a poorly understood molecular and cellular basis. We identify oncogenic deregulation of BCL11B driven by diverse structural alterations, including de novo superenhancer generation, as the driving feature of a subset of lineage-ambiguous leukemias that transcend current diagnostic boundaries.This article is highlighted in the In This Issue feature, p. 2659.

Molecular classification improves risk assessment in adult BCR-ABL1-negative B-ALL.

Genomic classification has improved risk assignment of pediatric, but not adult B-lineage acute lymphoblastic leukemia (B-ALL). The international UKALLXII/ECOG-ACRIN E2993 (#NCT00002514) trial accrued 1229 adolescent/adult patients with BCR-ABL1- B-ALL (aged 14 to 65 years). Although 93% of patients achieved remission, 41% relapsed at a median of 13 months (range, 28 days to 12 years). Five-year overall survival (OS) was 42% (95% confidence interval, 39, 44). Transcriptome sequencing, gene expression profiling, cytogenetics, and fusion polymerase chain reaction enabled genomic subtyping of 282 patient samples, of which 264 were eligible for trial, accounting for 64.5% of E2993 patients. Among patients with outcome data, 29.5% with favorable outcomes (5-year OS 65% to 80%) were deemed standard risk (DUX4-rearranged [9.2%], ETV6-RUNX1/-like [2.3%], TCF3-PBX1 [6.9%], PAX5 P80R [4.1%], high-hyperdiploid [6.9%]); 50.2% had high-risk genotypes with 5-year OS of 0% to 27% (Ph-like [21.2%], KMT2A-AFF1 [12%], low-hypodiploid/near-haploid [14.3%], BCL2/MYC-rearranged [2.8%]); 20.3% had intermediate-risk genotypes with 5-year OS of 33% to 45% (PAX5alt [12.4%], ZNF384/-like [5.1%], MEF2D-rearranged [2.8%]). IKZF1 alterations occurred in 86% of Ph-like, and TP53 mutations in patients who were low-hypodiploid (54%) and BCL2/MYC-rearranged (33%) but were not independently associated with outcome. Of patients considered high risk based on presenting age and white blood cell count, 40% harbored subtype-defining genetic alterations associated with standard- or intermediate-risk outcomes. We identified distinct immunophenotypic features for DUX4-rearranged, PAX5 P80R, ZNF384-R/-like, and Ph-like genotypes. These data in a large adult B-ALL cohort treated with a non-risk-adapted approach on a single trial show the prognostic importance of genomic analyses, which may translate into future therapeutic benefits.

Cis-regulatory mutations with driver hallmarks in major cancers.

Despite the recent availability of complete genome sequences of tumors from thousands of patients, isolating disease-causing (driver) non-coding mutations from the plethora of somatic variants remains challenging, and only a handful of validated examples exist. By integrating whole-genome sequencing, genetic data, and allele-specific gene expression from TCGA, we identified 320 somatic non-coding mutations that affect gene expression in cis (FDR<0.25). These mutations cluster into 47 cis-regulatory elements that modulate expression of their subject genes through diverse molecular mechanisms. We further show that these mutations have hallmark features of non-coding drivers; namely, that they preferentially disrupt transcription factor binding motifs, are associated with a selective advantage, increased oncogene expression and decreased tumor suppressor expression.

Genome-wide association analyses of post-traumatic stress disorder and its symptom subdomains in the Million Veteran Program.

We conducted genome-wide association analyses of over 250,000 participants of European (EUR) and African (AFR) ancestry from the Million Veteran Program using electronic health record-validated post-traumatic stress disorder (PTSD) diagnosis and quantitative symptom phenotypes. Applying genome-wide multiple testing correction, we identified three significant loci in European case-control analyses and 15 loci in quantitative symptom analyses. Genomic structural equation modeling indicated tight coherence of a PTSD symptom factor that shares genetic variance with a distinct internalizing (mood-anxiety-neuroticism) factor. Partitioned heritability indicated enrichment in several cortical and subcortical regions, and imputed genetically regulated gene expression in these regions was used to identify potential drug repositioning candidates. These results validate the biological coherence of the PTSD syndrome, inform its relationship to comorbid anxiety and depressive disorders and provide new considerations for treatment.

Vagus nerve stimulation optimized cardiomyocyte phenotype, sarcomere organization and energy metabolism in infarcted heart through FoxO3A-VEGF signaling.

Vagus nerve stimulation (VNS) restores autonomic balance, suppresses inflammation action and minimizes cardiomyocyte injury. However, little knowledge is known about the VNS' role in cardiomyocyte phenotype, sarcomere organization, and energy metabolism of infarcted hearts. VNS in vivo and acetylcholine (ACh) in vitro optimized the levels of α/ß-MHC and α-Actinin positive sarcomere organization in cardiomyocytes while reducing F-actin assembly of cardiomyocytes. Consistently, ACh improved glucose uptake while decreasing lipid deposition in myocytes, correlating both with the increase of Glut4 and CPT1α and the decrease of PDK4 in infarcted hearts in vivo and myocytes in vitro, attributing to improvement in both glycolysis by VEGF-A and lipid uptake by VEGF-B in response to Ach. This led to increased ATP levels accompanied by the repaired mitochondrial function and the decreased oxygen consumption. Functionally, VNS improved the left ventricular performance. In contrast, ACh-m/nAChR inhibitor or knockdown of VEGF-A/B by shRNA powerfully abrogated these effects mediated by VNS. On mechanism, ACh decreased the levels of nuclear translocation of FoxO3A in myocytes due to phosphorylation of FoxO3A by activating AKT. FoxO3A overexpression or knockdown could reverse the specific effects of ACh on the expression of VEGF-A/B, α/ß-MHC, Glut4, and CPT1α, sarcomere organization, glucose uptake and ATP production. Taken together, VNS optimized cardiomyocytes sarcomere organization and energy metabolism to improve heart function of the infarcted heart during the process of delaying and/or blocking the switch from compensated hypertrophy to decompensated heart failure, which were associated with activation of both P13K/AKT-FoxO3A-VEGF-A/B signaling cascade.

Association of OPRM1 Functional Coding Variant With Opioid Use Disorder: A Genome-Wide Association Study.

Importance: With the current opioid crisis, it is important to improve understanding of the biological mechanisms of opioid use disorder (OUD). Objectives: To detect genetic risk variants for OUD and determine genetic correlations and causal association with OUD and other traits. Design, Setting, and Participants: A genome-wide association study of electronic health record-defined OUD in the Million Veteran Program sample was conducted, comprising 8529 affected European American individuals and 71 200 opioid-exposed European American controls (defined by electronic health record trajectory analysis) and 4032 affected African American individuals and 26 029 opioid-exposed African American controls. Participants were enrolled from January 10, 2011, to May 21, 2018, with electronic health record data for OUD diagnosis from October 1, 1999, to February 7, 2018. Million Veteran Program results and additional OUD case-control genome-wide association study results from the Yale-Penn and Study of Addiction: Genetics and Environment samples were meta-analyzed (total numbers: European American individuals, 10 544 OUD cases and 72 163 opioid-exposed controls; African American individuals, 5212 cases and 26 876 controls). Data on Yale-Penn participants were collected from February 14, 1999, to April 1, 2017, and data on Study of Addiction: Genetics and Environment participants were collected from 1990 to 2007. The key result was replicated in 2 independent cohorts: proxy-phenotype buprenorphine treatment in the UK Biobank and newly genotyped Yale-Penn participants. Genetic correlations between OUD and other traits were tested, and mendelian randomization analysis was conducted to identify potential causal associations. Main Outcomes and Measures: Main outcomes were International Classification of Diseases, Ninth Revision-diagnosed OUD or International Statistical Classification of Diseases and Related Health Problems, Tenth Revision-diagnosed OUD (Million Veteran Program), and DSM-IV-defined opioid dependence (Yale-Penn and Study of Addiction: Genetics and Environment). Results: A total of 114 759 individuals (101 016 men [88%]; mean [SD] age, 60.1 [12.8] years) were included. In 82 707 European American individuals, a functional coding variant (rs1799971, encoding Asn40Asp) in OPRM1 (µ-opioid receptor gene, the main biological target for opioid drugs; OMIM 600018) reached genome-wide significance (G allele: ß = -0.066 [SE = 0.012]; P = 1.51 × 10-8). The finding was replicated in 2 independent samples. Single-nucleotide polymorphism-based heritability of OUD was 11.3% (SE = 1.8%). Opioid use disorder was genetically correlated with 83 traits, including multiple substance use traits, psychiatric illnesses, cognitive performance, and others. Mendelian randomization analysis revealed the following associations with OUD: risk of tobacco smoking, depression, neuroticism, worry neuroticism subcluster, and cognitive performance. No genome-wide significant association was detected for African American individuals or in transpopulation meta-analysis. Conclusions and Relevance: This genome-wide meta-analysis identified a significant association of OUD with an OPRM1 variant, which was replicated in 2 independent samples. Post-genome-wide association study analysis revealed associated pleiotropic characteristics. Recruitment of additional individuals with OUD for future studies-especially those of non-European ancestry-is a crucial next step in identifying additional significant risk loci.