Genetic drivers and cellular selection of female mosaic X chromosome loss.

Mosaic loss of the X chromosome (mLOX) is the most common clonal somatic alteration in leukocytes of female individuals1,2, but little is known about its genetic determinants or phenotypic consequences. Here, to address this, we used data from 883,574 female participants across 8 biobanks; 12% of participants exhibited detectable mLOX in approximately 2% of leukocytes. Female participants with mLOX had an increased risk of myeloid and lymphoid leukaemias. Genetic analyses identified 56 common variants associated with mLOX, implicating genes with roles in chromosomal missegregation, cancer predisposition and autoimmune diseases. Exome-sequence analyses identified rare missense variants in FBXO10 that confer a twofold increased risk of mLOX. Only a small fraction of associations was shared with mosaic Y chromosome loss, suggesting that distinct biological processes drive formation and clonal expansion of sex chromosome missegregation. Allelic shift analyses identified X chromosome alleles that are preferentially retained in mLOX, demonstrating variation at many loci under cellular selection. A polygenic score including 44 allelic shift loci correctly inferred the retained X chromosomes in 80.7% of mLOX cases in the top decile. Our results support a model in which germline variants predispose female individuals to acquiring mLOX, with the allelic content of the X chromosome possibly shaping the magnitude of clonal expansion.

Autoimmune alleles at the major histocompatibility locus modify melanoma susceptibility.

Autoimmunity and cancer represent two different aspects of immune dysfunction. Autoimmunity is characterized by breakdowns in immune self-tolerance, while impaired immune surveillance can allow for tumorigenesis. The class I major histocompatibility complex (MHC-I), which displays derivatives of the cellular peptidome for immune surveillance by CD8+ T cells, serves as a common genetic link between these conditions. As melanoma-specific CD8+ T cells have been shown to target melanocyte-specific peptide antigens more often than melanoma-specific antigens, we investigated whether vitiligo- and psoriasis-predisposing MHC-I alleles conferred a melanoma-protective effect. In individuals with cutaneous melanoma from both The Cancer Genome Atlas (n = 451) and an independent validation set (n = 586), MHC-I autoimmune-allele carrier status was significantly associated with a later age of melanoma diagnosis. Furthermore, MHC-I autoimmune-allele carriers were significantly associated with decreased risk of developing melanoma in the Million Veteran Program (OR = 0.962, p = 0.024). Existing melanoma polygenic risk scores (PRSs) did not predict autoimmune-allele carrier status, suggesting these alleles provide orthogonal risk-relevant information. Mechanisms of autoimmune protection were neither associated with improved melanoma-driver mutation association nor improved gene-level conserved antigen presentation relative to common alleles. However, autoimmune alleles showed higher affinity relative to common alleles for particular windows of melanocyte-conserved antigens and loss of heterozygosity of autoimmune alleles caused the greatest reduction in presentation for several conserved antigens across individuals with loss of HLA alleles. Overall, this study presents evidence that MHC-I autoimmune-risk alleles modulate melanoma risk unaccounted for by current PRSs.

Inherited myeloproliferative neoplasm risk affects haematopoietic stem cells.

Myeloproliferative neoplasms (MPNs) are blood cancers that are characterized by the excessive production of mature myeloid cells and arise from the acquisition of somatic driver mutations in haematopoietic stem cells (HSCs). Epidemiological studies indicate a substantial heritable component of MPNs that is among the highest known for cancers1. However, only a limited number of genetic risk loci have been identified, and the underlying biological mechanisms that lead to the acquisition of MPNs remain unclear. Here, by conducting a large-scale genome-wide association study (3,797 cases and 1,152,977 controls), we identify 17 MPN risk loci (P < 5.0 × 10-8), 7 of which have not been previously reported. We find that there is a shared genetic architecture between MPN risk and several haematopoietic traits from distinct lineages; that there is an enrichment for MPN risk variants within accessible chromatin of HSCs; and that increased MPN risk is associated with longer telomere length in leukocytes and other clonal haematopoietic states-collectively suggesting that MPN risk is associated with the function and self-renewal of HSCs. We use gene mapping to identify modulators of HSC biology linked to MPN risk, and show through targeted variant-to-function assays that CHEK2 and GFI1B have roles in altering the function of HSCs to confer disease risk. Overall, our results reveal a previously unappreciated mechanism for inherited MPN risk through the modulation of HSC function.

Correlates of suicidal behaviors and genetic risk among United States veterans with schizophrenia or bipolar I disorder.

Persons diagnosed with schizophrenia (SCZ) or bipolar I disorder (BPI) are at high risk for self-injurious behavior, suicidal ideation, and suicidal behaviors (SB). Characterizing associations between diagnosed health problems, prior pharmacological treatments, and polygenic scores (PGS) has potential to inform risk stratification. We examined self-reported SB and ideation using the Columbia Suicide Severity Rating Scale (C-SSRS) among 3,942 SCZ and 5,414 BPI patients receiving care within the Veterans Health Administration (VHA). These cross-sectional data were integrated with electronic health records (EHRs), and compared across lifetime diagnoses, treatment histories, follow-up screenings, and mortality data. PGS were constructed using available genomic data for related traits. Genome-wide association studies were performed to identify and prioritize specific loci. Only 20% of the veterans who reported SB had a corroborating ICD-9/10 EHR code. Among those without prior SB, more than 20% reported new-onset SB at follow-up. SB were associated with a range of additional clinical diagnoses, and with treatment with specific classes of psychotropic medications (e.g., antidepressants, antipsychotics, etc.). PGS for externalizing behaviors, smoking initiation, suicide attempt, and major depressive disorder were associated with SB. The GWAS for SB yielded no significant loci. Among individuals with a diagnosed mental illness, self-reported SB were strongly associated with clinical variables across several EHR domains. Analyses point to sequelae of substance-related and psychiatric comorbidities as strong correlates of prior and subsequent SB. Nonetheless, past SB was frequently not documented in health records, underscoring the value of regular screening with direct, in-person assessments, especially among high-risk individuals.

A multi-population phenome-wide association study of genetically-predicted height in the Million Veteran Program.

BACKGROUND: Height has been associated with many clinical traits but whether such associations are causal versus secondary to confounding remains unclear in many cases. To systematically examine this question, we performed a Mendelian Randomization-Phenome-wide association study (MR-PheWAS) using clinical and genetic data from a national healthcare system biobank. METHODS AND FINDINGS: Analyses were performed using data from the US Veterans Affairs (VA) Million Veteran Program in non-Hispanic White (EA, n = 222,300) and non-Hispanic Black (AA, n = 58,151) adults in the US. We estimated height genetic risk based on 3290 height-associated variants from a recent European-ancestry genome-wide meta-analysis. We compared associations of measured and genetically-predicted height with phenome-wide traits derived from the VA electronic health record, adjusting for age, sex, and genetic principal components. We found 345 clinical traits associated with measured height in EA and an additional 17 in AA. Of these, 127 were associated with genetically-predicted height at phenome-wide significance in EA and 2 in AA. These associations were largely independent from body mass index. We confirmed several previously described MR associations between height and cardiovascular disease traits such as hypertension, hyperlipidemia, coronary heart disease (CHD), and atrial fibrillation, and further uncovered MR associations with venous circulatory disorders and peripheral neuropathy in the presence and absence of diabetes. As a number of traits associated with genetically-predicted height frequently co-occur with CHD, we evaluated effect modification by CHD status of genetically-predicted height associations with risk factors for and complications of CHD. We found modification of effects of MR associations by CHD status for atrial fibrillation/flutter but not for hypertension, hyperlipidemia, or venous circulatory disorders. CONCLUSIONS: We conclude that height may be an unrecognized but biologically plausible risk factor for several common conditions in adults. However, more studies are needed to reliably exclude horizontal pleiotropy as a driving force behind at least some of the MR associations observed in this study.

A Phenome-Wide Association Study of genes associated with COVID-19 severity reveals shared genetics with complex diseases in the Million Veteran Program.

The study aims to determine the shared genetic architecture between COVID-19 severity with existing medical conditions using electronic health record (EHR) data. We conducted a Phenome-Wide Association Study (PheWAS) of genetic variants associated with critical illness (n = 35) or hospitalization (n = 42) due to severe COVID-19 using genome-wide association summary data from the Host Genetics Initiative. PheWAS analysis was performed using genotype-phenotype data from the Veterans Affairs Million Veteran Program (MVP). Phenotypes were defined by International Classification of Diseases (ICD) codes mapped to clinically relevant groups using published PheWAS methods. Among 658,582 Veterans, variants associated with severe COVID-19 were tested for association across 1,559 phenotypes. Variants at the ABO locus (rs495828, rs505922) associated with the largest number of phenotypes (nrs495828 = 53 and nrs505922 = 59); strongest association with venous embolism, odds ratio (ORrs495828 1.33 (p = 1.32 x 10-199), and thrombosis ORrs505922 1.33, p = 2.2 x10-265. Among 67 respiratory conditions tested, 11 had significant associations including MUC5B locus (rs35705950) with increased risk of idiopathic fibrosing alveolitis OR 2.83, p = 4.12 × 10-191; CRHR1 (rs61667602) associated with reduced risk of pulmonary fibrosis, OR 0.84, p = 2.26× 10-12. The TYK2 locus (rs11085727) associated with reduced risk for autoimmune conditions, e.g., psoriasis OR 0.88, p = 6.48 x10-23, lupus OR 0.84, p = 3.97 x 10-06. PheWAS stratified by ancestry demonstrated differences in genotype-phenotype associations. LMNA (rs581342) associated with neutropenia OR 1.29 p = 4.1 x 10-13 among Veterans of African and Hispanic ancestry but not European. Overall, we observed a shared genetic architecture between COVID-19 severity and conditions related to underlying risk factors for severe and poor COVID-19 outcomes. Differing associations between genotype-phenotype across ancestries may inform heterogenous outcomes observed with COVID-19. Divergent associations between risk for severe COVID-19 with autoimmune inflammatory conditions both respiratory and non-respiratory highlights the shared pathways and fine balance of immune host response and autoimmunity and caution required when considering treatment targets.

Regulatory variants in TCF7L2 are associated with thoracic aortic aneurysm.

Thoracic aortic aneurysm (TAA) is characterized by dilation of the aortic root or ascending/descending aorta. TAA is a heritable disease that can be potentially life threatening. While 10%-20% of TAA cases are caused by rare, pathogenic variants in single genes, the origin of the majority of TAA cases remains unknown. A previous study implicated common variants in FBN1 with TAA disease risk. Here, we report a genome-wide scan of 1,351 TAA-affected individuals and 18,295 control individuals from the Cardiovascular Health Improvement Project and Michigan Genomics Initiative at the University of Michigan. We identified a genome-wide significant association with TAA for variants within the third intron of TCF7L2 following replication with meta-analysis of four additional independent cohorts. Common variants in this locus are the strongest known genetic risk factor for type 2 diabetes. Although evidence indicates the presence of different causal variants for TAA and type 2 diabetes at this locus, we observed an opposite direction of effect. The genetic association for TAA colocalizes with an aortic eQTL of TCF7L2, suggesting a functional relationship. These analyses predict an association of higher expression of TCF7L2 with TAA disease risk. In vitro, we show that upregulation of TCF7L2 is associated with BCL2 repression promoting vascular smooth muscle cell apoptosis, a key driver of TAA disease.

Design and Pilot Results from Million Veteran Program Return Of Actionable Genetic Results (MVP-ROAR) Study.

BACKGROUND: As a mega-biobank linked to a national healthcare system, the Million Veteran Program (MVP) can directly improve the health care of participants. To determine the feasibility and outcomes of returning medically actionable genetic results to MVP participants, the program launched the MVP Return Of Actionable Results (MVP-ROAR) Study, with familial hypercholesterolemia (FH) as an exemplar actionable condition. METHODS: The MVP-ROAR Study consists of a completed single-arm pilot phase and an ongoing randomized clinical trial (RCT), in which MVP participants are recontacted and invited to receive clinical confirmatory gene sequencing testing and a telegenetic counseling intervention. The primary outcome of the RCT is 6-month change in low-density lipoprotein cholesterol (LDL-C) between participants receiving results at baseline and those receiving results after 6 months. RESULTS: The pilot developed processes to identify and recontact participants nationally with probable pathogenic variants in low-density lipoprotein receptor (LDLR) on the MVP genotype array, invite them to clinical confirmatory gene sequencing, and deliver a telegenetic counseling intervention. Among participants in the pilot phase, 8 (100%) had active statin prescriptions after 6 months. Results were shared with 16 first-degree family members. Six-month ΔLDL-C (low-density lipoprotein cholesterol) after the genetic counseling intervention was -37 mg/dL (95% CI: -12 to -61; p=0.03). The ongoing RCT will determine between-arm differences in this primary outcome. CONCLUSION: While underscoring the importance of clinical confirmation of research results, the pilot phase of the MVP-ROAR Study marks a turning point in MVP and demonstrates the feasibility of returning genetic results to participants and their providers. The ongoing RCT will contribute to understanding how such a program might improve patient health care and outcomes.

Modeling the longitudinal changes of ancestry diversity in the Million Veteran Program.

BACKGROUND: The Million Veteran Program (MVP) participants represent 100 years of US history, including significant social and demographic changes over time. Our study assessed two aspects of the MVP: (i) longitudinal changes in population diversity and (ii) how these changes can be accounted for in genome-wide association studies (GWAS). To investigate these aspects, we divided MVP participants into five birth cohorts (N-range = 123,888 [born from 1943 to 1947] to 136,699 [born from 1948 to 1953]). RESULTS: Ancestry groups were defined by (i) HARE (harmonized ancestry and race/ethnicity) and (ii) a random-forest clustering approach using the 1000 Genomes Project and the Human Genome Diversity Project (1kGP + HGDP) reference panels (77 world populations representing six continental groups). In these groups, we performed GWASs of height, a trait potentially affected by population stratification. Birth cohorts demonstrate important trends in ancestry diversity over time. More recent HARE-assigned Europeans, Africans, and Hispanics had lower European ancestry proportions than older birth cohorts (0.010 < Cohen's d < 0.259, p < 7.80 × 10-4). Conversely, HARE-assigned East Asians showed an increase in European ancestry proportion over time. In GWAS of height using HARE assignments, genomic inflation due to population stratification was prevalent across all birth cohorts (linkage disequilibrium score regression intercept = 1.08 ± 0.042). The 1kGP + HGDP-based ancestry assignment significantly reduced the population stratification (mean intercept reduction = 0.045 ± 0.007, p < 0.05) confounding in the GWAS statistics. CONCLUSIONS: This study provides a characterization of ancestry diversity of the MVP cohort over time and compares two strategies to infer genetically defined ancestry groups by assessing differences in controlling population stratification in genome-wide association studies.

Genotyping Array Design and Data Quality Control in the Million Veteran Program.

The Million Veteran Program (MVP), initiated by the Department of Veterans Affairs (VA), aims to collect biosamples with consent from at least one million veterans. Presently, blood samples have been collected from over 800,000 enrolled participants. The size and diversity of the MVP cohort, as well as the availability of extensive VA electronic health records, make it a promising resource for precision medicine. MVP is conducting array-based genotyping to provide a genome-wide scan of the entire cohort, in parallel with whole-genome sequencing, methylation, and other 'omics assays. Here, we present the design and performance of the MVP 1.0 custom Axiom array, which was designed and developed as a single assay to be used across the multi-ethnic MVP cohort. A unified genetic quality-control analysis was developed and conducted on an initial tranche of 485,856 individuals, leading to a high-quality dataset of 459,777 unique individuals. 668,418 genetic markers passed quality control and showed high-quality genotypes not only on common variants but also on rare variants. We confirmed that, with non-European individuals making up nearly 30%, MVP's substantial ancestral diversity surpasses that of other large biobanks. We also demonstrated the quality of the MVP dataset by replicating established genetic associations with height in European Americans and African Americans ancestries. This current dataset has been made available to approved MVP researchers for genome-wide association studies and other downstream analyses. Further data releases will be available for analysis as recruitment at the VA continues and the cohort expands both in size and diversity.

A MUC5B Gene Polymorphism, rs35705950-T, Confers Protective Effects Against COVID-19 Hospitalization but Not Severe Disease or Mortality.

Rationale: A common MUC5B gene polymorphism, rs35705950-T, is associated with idiopathic pulmonary fibrosis (IPF), but its role in severe acute respiratory syndrome coronavirus 2 infection and disease severity is unclear. Objectives: To assess whether rs35705950-T confers differential risk for clinical outcomes associated with coronavirus disease (COVID-19) infection among participants in the Million Veteran Program (MVP). Methods: The MUC5B rs35705950-T allele was directly genotyped among MVP participants; clinical events and comorbidities were extracted from the electronic health records. Associations between the incidence or severity of COVID-19 and rs35705950-T were analyzed within each ancestry group in the MVP followed by transancestry meta-analysis. Replication and joint meta-analysis were conducted using summary statistics from the COVID-19 Host Genetics Initiative (HGI). Sensitivity analyses with adjustment for additional covariates (body mass index, Charlson comorbidity index, smoking, asbestosis, rheumatoid arthritis with interstitial lung disease, and IPF) and associations with post-COVID-19 pneumonia were performed in MVP subjects. Measurements and Main Results: The rs35705950-T allele was associated with fewer COVID-19 hospitalizations in transancestry meta-analyses within the MVP (Ncases = 4,325; Ncontrols = 507,640; OR = 0.89 [0.82-0.97]; P = 6.86 × 10-3) and joint meta-analyses with the HGI (Ncases = 13,320; Ncontrols = 1,508,841; OR, 0.90 [0.86-0.95]; P = 8.99 × 10-5). The rs35705950-T allele was not associated with reduced COVID-19 positivity in transancestry meta-analysis within the MVP (Ncases = 19,168/Ncontrols = 492,854; OR, 0.98 [0.95-1.01]; P = 0.06) but was nominally significant (P < 0.05) in the joint meta-analysis with the HGI (Ncases = 44,820; Ncontrols = 1,775,827; OR, 0.97 [0.95-1.00]; P = 0.03). Associations were not observed with severe outcomes or mortality. Among individuals of European ancestry in the MVP, rs35705950-T was associated with fewer post-COVID-19 pneumonia events (OR, 0.82 [0.72-0.93]; P = 0.001). Conclusions: The MUC5B variant rs35705950-T may confer protection in COVID-19 hospitalizations.

Harmonizing Genetic Ancestry and Self-identified Race/Ethnicity in Genome-wide Association Studies.

Large-scale multi-ethnic cohorts offer unprecedented opportunities to elucidate the genetic factors influencing complex traits related to health and disease among minority populations. At the same time, the genetic diversity in these cohorts presents new challenges for analysis and interpretation. We consider the utility of race and/or ethnicity categories in genome-wide association studies (GWASs) of multi-ethnic cohorts. We demonstrate that race/ethnicity information enhances the ability to understand population-specific genetic architecture. To address the practical issue that self-identified racial/ethnic information may be incomplete, we propose a machine learning algorithm that produces a surrogate variable, termed HARE. We use height as a model trait to demonstrate the utility of HARE and ethnicity-specific GWASs.

Genome-wide association study of cognitive performance in U.S. veterans with schizophrenia or bipolar disorder.

Cognitive impairment is a frequent and serious problem in patients with various forms of severe mental illnesses (SMI), including schizophrenia (SZ) and bipolar disorder (BP). Recent research suggests genetic links to several cognitive phenotypes in both SMI and in the general population. Our goal in this study was to identify potential genomic signatures of cognitive functioning in veterans with severe mental illness and compare them to previous findings for cognition across different populations. Veterans Affairs (VA) Cooperative Studies Program (CSP) Study #572 evaluated cognitive and functional capacity measures among SZ and BP patients. In conjunction with the VA Million Veteran Program, 3,959 European American (1,095 SZ, 2,864 BP) and 2,601 African American (1,095 SZ, 2,864 BP) patients were genotyped using a custom Affymetrix Axiom Biobank array. We performed a genome-wide association study of global cognitive functioning, constructed polygenic scores for SZ and cognition in the general population, and examined genetic correlations with 2,626 UK Biobank traits. Although no single locus attained genome-wide significance, observed allelic effects were strongly consistent with previous studies. We observed robust associations between global cognitive functioning and polygenic scores for cognitive performance, intelligence, and SZ risk. We also identified significant genetic correlations with several cognition-related traits in UK Biobank. In a diverse cohort of U.S. veterans with SZ or BP, we demonstrate broad overlap of common genetic effects on cognition in the general population, and find that greater polygenic loading for SZ risk is associated with poorer cognitive performance.

Impact of genotyping errors on statistical power of association tests in genomic analyses: A case study.

A key step in genomic studies is to assess high throughput measurements across millions of markers for each participant's DNA, either using microarrays or sequencing techniques. Accurate genotype calling is essential for downstream statistical analysis of genotype-phenotype associations, and next generation sequencing (NGS) has recently become a more common approach in genomic studies. How the accuracy of variant calling in NGS-based studies affects downstream association analysis has not, however, been studied using empirical data in which both microarrays and NGS were available. In this article, we investigate the impact of variant calling errors on the statistical power to identify associations between single nucleotides and disease, and on associations between multiple rare variants and disease. Both differential and nondifferential genotyping errors are considered. Our results show that the power of burden tests for rare variants is strongly influenced by the specificity in variant calling, but is rather robust with regard to sensitivity. By using the variant calling accuracies estimated from a substudy of a Cooperative Studies Program project conducted by the Department of Veterans Affairs, we show that the power of association tests is mostly retained with commonly adopted variant calling pipelines. An R package, GWAS.PC, is provided to accommodate power analysis that takes account of genotyping errors (http://zhaocenter.org/software/).

Development and utility of a clinical research informatics application for participant recruitment and workflow management for a return of results pilot trial in familial hypercholesterolemia in the Million Veteran Program.

Objective: The development of clinical research informatics tools and workflow processes associated with re-engaging biobank participants has become necessary as genomic repositories increasingly consider the return of actionable research results. Materials and Methods: Here we describe the development and utility of an informatics application for participant recruitment and enrollment management for the Veterans Affairs Million Veteran Program Return Of Actionable Results Study, a randomized controlled pilot trial returning individual genetic results associated with familial hypercholesterolemia. Results: The application is developed in Python-Flask and was placed into production in November 2021. The application includes modules for chart review, medication reconciliation, participant contact and biospecimen logging, survey recording, randomization, and documentation of genetic counseling and result disclosure. Three primary users, a genetic counselor and two research coordinators, and 326 Veteran participants have been integrated into the system as of February 23, 2023. The application has successfully handled 3367 task requests involving greater than 95 000 structured data points. Specifically, application users have recorded 326 chart reviews, 867 recruitment telephone calls, 158 telephone-based surveys, and 61 return of results genetic counseling sessions, among other available study tasks. Conclusion: The development of usable, customizable, and secure informatics tools will become increasingly important as large genomic repositories begin to return research results at scale. Our work provides a proof-of-concept for developing and using such tools to aid in managing the return of results process within a national biobank.

A multi-ancestry GWAS of Fuchs corneal dystrophy highlights the contributions of laminins, collagen, and endothelial cell regulation.

Fuchs endothelial corneal dystrophy (FECD) is a leading indication for corneal transplantation, but its molecular etiology remains poorly understood. We performed genome-wide association studies (GWAS) of FECD in the Million Veteran Program followed by multi-ancestry meta-analysis with the previous largest FECD GWAS, for a total of 3970 cases and 333,794 controls. We confirm the previous four loci, and identify eight novel loci: SSBP3, THSD7A, LAMB1, PIDD1, RORA, HS3ST3B1, LAMA5, and COL18A1. We further confirm the TCF4 locus in GWAS for admixed African and Hispanic/Latino ancestries and show an enrichment of European-ancestry haplotypes at TCF4 in FECD cases. Among the novel associations are low frequency missense variants in laminin genes LAMA5 and LAMB1 which, together with previously reported LAMC1, form laminin-511 (LM511). AlphaFold 2 protein modeling, validated through homology, suggests that mutations at LAMA5 and LAMB1 may destabilize LM511 by altering inter-domain interactions or extracellular matrix binding. Finally, phenome-wide association scans and colocalization analyses suggest that the TCF4 CTG18.1 trinucleotide repeat expansion leads to dysregulation of ion transport in the corneal endothelium and has pleiotropic effects on renal function.

Neighborhood Deprivation, Race and Ethnicity, and Prostate Cancer Outcomes Across California Health Care Systems.

Importance: Non-Hispanic Black (hereafter, Black) individuals experience worse prostate cancer outcomes due to socioeconomic and racial inequities of access to care. Few studies have empirically evaluated these disparities across different health care systems. Objective: To describe the racial and ethnic and neighborhood socioeconomic status (nSES) disparities among residents of the same communities who receive prostate cancer care in the US Department of Veterans Affairs (VA) health care system vs other settings. Design, Setting, and Participants: This cohort study obtained data from the VA Central Cancer Registry for veterans with prostate cancer who received care within the VA Greater Los Angeles Healthcare System (VA cohort) and from the California Cancer Registry (CCR) for nonveterans who received care outside the VA setting (CCR cohort). The cohorts consisted of all males with incident prostate cancer who were living within the same US Census tracts. These individuals received care between 2000 and 2018 and were followed up until death from any cause or censoring on December 31, 2018. Data analyses were conducted between September 2022 and December 2023. Exposures: Health care setting, self-identified race and ethnicity (SIRE), and nSES. Main Outcomes and Measures: The primary outcome was all-cause mortality (ACM). Cox proportional hazards regression models were used to estimate hazard ratios for associations of SIRE and nSES with prostate cancer outcomes in the VA and CCR cohorts. Results: Included in the analysis were 49 461 males with prostate cancer. Of these, 1881 males were in the VA cohort (mean [SD] age, 65.3 [7.7] years; 833 Black individuals [44.3%], 694 non-Hispanic White [hereafter, White] individuals [36.9%], and 354 individuals [18.8%] of other or unknown race). A total of 47 580 individuals were in the CCR cohort (mean [SD] age, 67.0 [9.6] years; 8183 Black individuals [17.2%], 26 206 White individuals [55.1%], and 13 191 individuals [27.8%] of other or unknown race). In the VA cohort, there were no racial disparities observed for metastasis, ACM, or prostate cancer-specific mortality (PCSM). However, in the CCR cohort, the racial disparities were observed for metastasis (adjusted odds ratio [AOR], 1.36; 95% CI, 1.22-1.52), ACM (adjusted hazard ratio [AHR], 1.13; 95% CI, 1.04-1.24), and PCSM (AHR, 1.15; 95% CI, 1.05-1.25). Heterogeneity was observed for the racial disparity in ACM in the VA vs CCR cohorts (AHR, 0.90 [95% CI, 0.76-1.06] vs 1.13 [95% CI, 1.04-1.24]; P = .01). No evidence of nSES disparities was observed for any prostate cancer outcomes in the VA cohort. However, in the CCR cohort, heterogeneity was observed for nSES disparities with ACM (AHR, 0.82; 95% CI, 0.80-0.84; P = .002) and PCSM (AHR, 0.86; 95% CI, 0.82-0.89; P = .007). Conclusions and Relevance: Results of this study suggest that racial and nSES disparities were wider among patients seeking care outside of the VA health care system. Health systems-related interventions that address access barriers may mitigate racial and socioeconomic disparities in prostate cancer.

Rare genetic variation in VE-PTP is associated with central serous chorioretinopathy, venous dysfunction and glaucoma.

Central serous chorioretinopathy (CSC) is a fluid maculopathy whose etiology is not well understood. Abnormal choroidal veins in CSC patients have been shown to have similarities with varicose veins. To identify potential mechanisms, we analyzed genotype data from 1,477 CSC patients and 455,449 controls in FinnGen. We identified an association for a low-frequency (AF=0.5%) missense variant (rs113791087) in the gene encoding vascular endothelial protein tyrosine phosphatase (VE-PTP) (OR=2.85, P=4.5×10-9). This was confirmed in a meta-analysis of 2,452 CSC patients and 865,767 controls from 4 studies (OR=3.06, P=7.4×10-15). Rs113791087 was associated with a 56% higher prevalence of retinal abnormalities (35.3% vs 22.6%, P=8.0×10-4) in 708 UK Biobank participants and, surprisingly, with varicose veins (OR=1.31, P=2.3×10-11) and glaucoma (OR=0.82, P=6.9×10-9). Predicted loss-of-function variants in VEPTP, though rare in number, were associated with CSC in All of Us (OR=17.10, P=0.018). These findings highlight the significance of VE-PTP in diverse ocular and systemic vascular diseases.