A concerted neuron-astrocyte program declines in ageing and schizophrenia.

Human brains vary across people and over time; such variation is not yet understood in cellular terms. Here we describe a relationship between people's cortical neurons and cortical astrocytes. We used single-nucleus RNA sequencing to analyse the prefrontal cortex of 191 human donors aged 22-97 years, including healthy individuals and people with schizophrenia. Latent-factor analysis of these data revealed that, in people whose cortical neurons more strongly expressed genes encoding synaptic components, cortical astrocytes more strongly expressed distinct genes with synaptic functions and genes for synthesizing cholesterol, an astrocyte-supplied component of synaptic membranes. We call this relationship the synaptic neuron and astrocyte program (SNAP). In schizophrenia and ageing-two conditions that involve declines in cognitive flexibility and plasticity1,2-cells divested from SNAP: astrocytes, glutamatergic (excitatory) neurons and GABAergic (inhibitory) neurons all showed reduced SNAP expression to corresponding degrees. The distinct astrocytic and neuronal components of SNAP both involved genes in which genetic risk factors for schizophrenia were strongly concentrated. SNAP, which varies quantitatively even among healthy people of similar age, may underlie many aspects of normal human interindividual differences and may be an important point of convergence for multiple kinds of pathophysiology.

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.

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.

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.

A Neanderthal haplotype introgressed into the human genome confers protection against membranous nephropathy.

Class 2 HLA and PLA2R1 alleles are exceptionally strong genetic risk factors for membranous nephropathy (MN), leading, through an unknown mechanism, to a targeted autoimmune response. Introgressed archaic haplotypes (introduced from an archaic human genome into the modern human genome) might influence phenotypes through gene dysregulation. Here, we investigated the genomic region surrounding the PLA2R1 gene. We reconstructed the phylogeny of Neanderthal and modern haplotypes in this region and calculated the probability of the observed clustering being the result of introgression or common descent. We imputed variants for the participants in our previous genome-wide association study and we compared the distribution of Neanderthal variants between MN cases and controls. The region associated with the lead MN risk locus in the PLA2R1 gene was confirmed and showed that, within a 507 kb region enriched in introgressed sequence, a stringently defined 105 kb haplotype, intersecting the coding regions for PLA2R1 and ITGB6, is inherited from Neanderthals. Thus, introgressed Neanderthal haplotypes overlapping PLA2R1 are differentially represented in MN cases and controls, with enrichment In controls suggesting a protective effect.

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.

Concerted neuron-astrocyte gene expression declines in aging and schizophrenia.

Human brains vary across people and over time; such variation is not yet understood in cellular terms. Here we describe a striking relationship between people's cortical neurons and cortical astrocytes. We used single-nucleus RNA-seq to analyze the prefrontal cortex of 191 human donors ages 22-97 years, including healthy individuals and persons with schizophrenia. Latent-factor analysis of these data revealed that in persons whose cortical neurons more strongly expressed genes for synaptic components, cortical astrocytes more strongly expressed distinct genes with synaptic functions and genes for synthesizing cholesterol, an astrocyte-supplied component of synaptic membranes. We call this relationship the Synaptic Neuron-and-Astrocyte Program (SNAP). In schizophrenia and aging - two conditions that involve declines in cognitive flexibility and plasticity 1,2 - cells had divested from SNAP: astrocytes, glutamatergic (excitatory) neurons, and GABAergic (inhibitory) neurons all reduced SNAP expression to corresponding degrees. The distinct astrocytic and neuronal components of SNAP both involved genes in which genetic risk factors for schizophrenia were strongly concentrated. SNAP, which varies quantitatively even among healthy persons of similar age, may underlie many aspects of normal human interindividual differences and be an important point of convergence for multiple kinds of pathophysiology.

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.

Tumor Predisposing Post-Zygotic Chromosomal Alterations in Bladder Cancer-Insights from Histologically Normal Urothelium.

Bladder urothelial carcinoma (BLCA) is the 10th most common cancer with a low survival rate and strong male bias. We studied the field cancerization in BLCA using multi-sample- and multi-tissue-per-patient protocol for sensitive detection of autosomal post-zygotic chromosomal alterations and loss of chromosome Y (LOY). We analysed 277 samples of histologically normal urothelium, 145 tumors and 63 blood samples from 52 males and 15 females, using the in-house adapted Mosaic Chromosomal Alterations (MoChA) pipeline. This approach allows identification of the early aberrations in urothelium from BLCA patients. Overall, 45% of patients exhibited at least one alteration in at least one normal urothelium sample. Recurrence analysis resulted in 16 hotspots composed of either gains and copy number neutral loss of heterozygosity (CN-LOH) or deletions and CN-LOH, encompassing well-known and new BLCA cancer driver genes. Conservative assessment of LOY showed 29%, 27% and 18% of LOY-cells in tumors, blood and normal urothelium, respectively. We provide a proof of principle that our approach can characterize the earliest alterations preconditioning normal urothelium to BLCA development. Frequent LOY in blood and urothelium-derived tissues suggest its involvement in BLCA.

Evaluation of imputation performance of multiple reference panels in a Pakistani population.

Genotype imputation is crucial for GWAS, but reference panels and existing benchmarking studies prioritize European individuals. Consequently, it is unclear which publicly available reference panel should be used for Pakistani individuals, and whether ancestry composition or sample size of the panel matters more for imputation accuracy. Our study compared different reference panels to impute genotype data in 1814 Pakistani individuals, finding the best performance balancing accuracy and coverage with meta-imputation with TOPMed and the expanded 1000 Genomes (ex1KG) reference. Imputation accuracy of ex1KG outperformed TOPMed despite its 30-fold smaller sample size, supporting efforts to create future panels with diverse populations.