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.

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.

Multi-ancestry genome-wide association study of kidney cancer identifies 63 susceptibility regions.

Here, in a multi-ancestry genome-wide association study meta-analysis of kidney cancer (29,020 cases and 835,670 controls), we identified 63 susceptibility regions (50 novel) containing 108 independent risk loci. In analyses stratified by subtype, 52 regions (78 loci) were associated with clear cell renal cell carcinoma (RCC) and 6 regions (7 loci) with papillary RCC. Notably, we report a variant common in African ancestry individuals ( rs7629500 ) in the 3' untranslated region of VHL, nearly tripling clear cell RCC risk (odds ratio 2.72, 95% confidence interval 2.23-3.30). In cis-expression quantitative trait locus analyses, 48 variants from 34 regions point toward 83 candidate genes. Enrichment of hypoxia-inducible factor-binding sites underscores the importance of hypoxia-related mechanisms in kidney cancer. Our results advance understanding of the genetic architecture of kidney cancer, provide clues for functional investigation and enable generation of a validated polygenic risk score with an estimated area under the curve of 0.65 (0.74 including risk factors) among European ancestry individuals.

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.

An extra X chromosome among adult women in the Million Veteran Program: A more benign perspective of trisomy X.

Despite affecting in 1 in every 1000 females, remarkably little is known about trisomy X syndrome (47,XXX), especially among older adults who are undiagnosed. In this study, we aimed to determine the prevalence of 47,XXX among females enrolled in the Million Veterans Program (MVP; mean age 50.2 ± 13.6 years), and compare broad health outcomes between females with 47,XXX and 46,XX matched controls. We identified 61 females with an additional X chromosome, corresponding to a prevalence of 103 per 100,000 females; 27.9% had been clinically diagnosed. Females with 47,XXX had taller stature (+6.1 cm, p < 0.001), greater rate of outpatient encounters (p = 0.026), higher odds of kidney disease (odds ratio [OR] = 12.3; 95% confidence interval [CI] 2.9-51.8), glaucoma (OR = 5.1; 95% CI 1.5-13.9), and congestive heart failure (OR = 5.6; 95% CI 1.4-24.2), and were more likely to be unemployed (p = 0.008) with lower annual income (p = 0.021) when compared with 46,XX controls of the same age and genetic ancestry. However, there were no differences in the rates of other encounter types, Charlson Comorbidity Index, all other medical and psychological diagnoses, military service history or quality of life metrics. In conclusion, in this aging and predominately undiagnosed sample, 47,XXX conferred few differences when compared with matched controls, offering a more reassuring perspective to the trisomy X literature.

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.

Prevalence, Morbidity, and Mortality of Men With Sex Chromosome Aneuploidy in the Million Veteran Program Cohort.

Importance: The reported phenotypes of men with 47,XXY and 47,XYY syndromes include tall stature, multisystem comorbidities, and poor health-related quality of life (HRQOL). However, knowledge about these sex chromosome aneuploidy (SCA) conditions has been derived from studies in the less than 15% of patients who are clinically diagnosed and also lack diversity in age and genetic ancestry. Objectives: To determine the prevalence of clinically diagnosed and undiagnosed X or Y chromosome aneuploidy among men enrolled in the Million Veteran Program (MVP); to describe military service metrics of men with SCAs; and to compare morbidity and mortality outcomes between men with SCA with and without a clinical diagnosis vs matched controls. Design, Setting, and Participants: This cross-sectional study used a case-control recruitment design to select biological males enrolled in the MVP biobank in the US Veterans Administration health care system from 2011 to 2022. Cases were participants with 47,XXY syndrome or 47,XYY syndrome, matched 1:5 with controls based on sex, age, and genetic ancestry. Data were analyzed from January 2022 to December 2023. Exposure: Genomic identification of an additional X or Y chromosome. Main Outcomes and Measures: Outcomes of interest included prevalence of men with SCAs from genomic analysis; clinical SCA diagnosis; Charlson Comorbidity Index; rates of outpatient, inpatient, and emergency encounters per year; self-reported health outcomes; and standardized mortality ratio. Results: Of 595 612 genotyped males in the MVP, 862 had an additional X chromosome (47,XXY) and 747 had an extra Y chromosome (47,XYY), with the highest prevalence among men with East Asian (47,XXY: 10 of 7313 participants; 47,XYY: 14 of 7313 participants) and European (47,XXY: 725 of 427 143 participants; 47,XYY: 625 of 427 143 participants) ancestry. Mean (SD) age at assessment was 61 (12) years, at which point 636 veterans (74.X%) with 47,XXY and 745 veterans (99%) with 47,XYY remained undiagnosed. Individuals with 47,XXY and 47,XYY had similar military service history, all-cause standardized mortality ratio, and age of death compared with matched controls. Individuals with SCA, compared with controls, had higher Charlson Comorbidity Index scores (47,XXY: mean [SD], 4.30 [2.72] vs controls: mean [SD], 3.90 [2.47]; 47,XYY: mean [SD], 4.45 [2.90] vs controls: mean [SD], 3.82 [2.50]) and health care utilization (eg, median [IQR] outpatient encounters per year: 47,XXY, 22.6 [11.8-37.8] vs controls, 16.8 [9.4-28]; 47,XYY: 21.4 [12.4-33.8] vs controls: 17.0 [9.4-28.2]), while several measures of HRQOL were lower (eg, mean [SD] self-reported physical function: 47,XXY: 34.2 [12] vs control mean [SD] 37.8 [12.8]; 47,XYY: 36.3 [11.6] vs control 37.9 [12.8]). Men with a clinical diagnosis of 47,XXY, compared with individuals without a clinical diagnosis, had higher health care utilization (eg, median [IQR] encounters per year: 26.6 [14.9-43.2] vs 22.2 [11.3-36.0]) but lower Charlson Comorbidity Index scores (mean [SD]: 3.7 [2.7] vs 4.5 [4.1]). Conclusion and Relevance: In this case-control study of men with 47,XXY and 47,XYY syndromes, prevalence of SCA was comparable with estimates in the general population. While these men had successfully served in the military, they had higher morbidity and reported poorer HRQOL with aging. Longer longitudinal follow-up of this sample will be informative for clinical and patient-reported outcomes, the role of ancestry, and mortality statistics.

BCFtools/liftover: an accurate and comprehensive tool to convert genetic variants across genome assemblies.

MOTIVATION: Many genetics studies report results tied to genomic coordinates of a legacy genome assembly. However, as assemblies are updated and improved, researchers are faced with either realigning raw sequence data using the updated coordinate system or converting legacy datasets to the updated coordinate system to be able to combine results with newer datasets. Currently available tools to perform the conversion of genetic variants have numerous shortcomings, including poor support for indels and multi-allelic variants, that lead to a higher rate of variants being dropped or incorrectly converted. As a result, many researchers continue to work with and publish using legacy genomic coordinates. RESULTS: Here we present BCFtools/liftover, a tool to convert genomic coordinates across genome assemblies for variants encoded in the variant call format with improved support for indels represented by different reference alleles across genome assemblies and full support for multi-allelic variants. It further supports variant annotation fields updates whenever the reference allele changes across genome assemblies. The tool has the lowest rate of variants being dropped with an order of magnitude less indels dropped or incorrectly converted and is an order of magnitude faster than other tools typically used for the same task. It is particularly suited for converting variant callsets from large cohorts to novel telomere-to-telomere assemblies as well as summary statistics from genome-wide association studies tied to legacy genome assemblies. AVAILABILITY AND IMPLEMENTATION: The tool is written in C and freely available under the MIT open source license as a BCFtools plugin available at http://github.com/freeseek/score.