Multiple sclerosis (MS) is a chronic inflammatory and neurodegenerative disease. Previous genetic studies have identified susceptibility loci that primarily impact immune cells and microglia. Here we performed a multi-ancestry genome-wide association study of 20,831 MS cases and 729,220 controls and identified 236 susceptibility variants outside of the major histocompatibility complex, including four novel genomic loci. We also derived a polygenic score for MS; while optimized for European ancestry, it is informative for African American and Latino individuals. Integrating single-cell data from blood and brain tissue, we identified 76 candidate causal genes. Inhibitory neurons emerged as a key target cell type for MS-associated variants, with seven loci, including STAT3, displaying altered expression only in these cells. The STAT3 variant is also associated with cognition and white matter integrity in individuals with no MS and greater sNfL levels in individuals with MS, suggesting that MS susceptibility may reflect reduced central nervous system resilience to inflammatory challenges. Multi-ancestry genome-wide association analyses integrated with single-cell data from blood and brain tissue identify genetic risk variants influencing multiple sclerosis susceptibility and highlight inhibitory neurons as a key target cell type.
Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system characterized by demyelination disseminated in space and time. Here we performed a genome-wide association study (GWAS) using 688 MS cases and 205,199 controls from the Japanese population and identified significant associations in the major histocompatibility complex region and a population-specific risk variant in 11q24. Through cross-population GWAS meta-analyses using a total of 29,374 cases and 1,843,563 controls from 4 ancestral populations, we identified 22 novel susceptibility loci. Integration of GWAS and single-cell and single-nucleus RNA sequencing of peripheral blood mononuclear cells and subcortical lesions from patients with MS revealed enrichment of genetic risk factors for MS in CD4+ T helper cell lineage and regulatory T cells, as well as in endothelial cells. Furthermore, spatial transcriptomics of subcortical lesions demonstrated spatial and temporal heterogeneity in associations with MS genetic risk. Our study demonstrates the value of investigation of spatiocellular features of disease genetics across diverse populations and omics modalities. Multiancestry genome-wide association analyses of multiple sclerosis integrated with single-cell, single-nucleus and spatial transcriptomics data yield insights into the cellular, spatial and temporal heterogeneity underlying disease risk.
War threatens scientific continuity by severing the mentorship networks that reproduce expertise. The Ukrainian Biological Data Science Summer School, held in person in western Ukraine each year since 2023, shows how distributed training, recursive mentorship and non-extractive international collaboration can preserve biological data science capacity under systemic stress.
Tumor progression is driven by dynamic interactions between cancer cells and their surrounding microenvironment. Here we integrate high-resolution spatial transcriptomics and evolving lineage-tracing technologies to elucidate how tumor expansion, plasticity and metastasis co-evolve with microenvironmental remodeling in a Kras;Trp53-driven mouse model of lung adenocarcinoma. We find that subclonal expansion contributes to a hypoxic, immunosuppressive and fibrotic microenvironment that is associated with the emergence of prometastatic cancer cell states. We use tumor phylogeography to delineate intercellular interactions that are rewired in the expanding tumor niche and use co-culture systems to dissect how intercellular interactions and hypoxia influence cancer cell state. Furthermore, we find that metastases arise from spatially confined primary tumor subclones and remodel the distant metastatic niche into a fibrotic, collagen-rich microenvironment. Together, we present a comprehensive dataset integrating spatial assays and lineage tracing to elucidate how sequential changes in cancer cell state and microenvironmental structures cooperate to promote tumor progression.
Multiple myeloma (MM) is a plasma-cell malignancy with extensive genomic and transcriptional heterogeneity, limiting disease classification and precision therapy. Here we generated a clinically annotated, population-scale, single-cell atlas of MM from 341 individuals spanning the disease and treatment continuum. We identified five recurrent malignant transcriptional archetypes and an orthogonal proliferative program associated with genomic features, therapeutic resistance and clinical outcomes. Validation in the independent CoMMpass cohort demonstrated robustness, prognostic relevance and portability across platforms. We developed a single-cell, target-discovery pipeline prioritizing malignant enrichment, cell-type specificity and tissue restriction, identifying FCRL2 as a plasma-restricted or B cell-lineage-restricted surface target expressed by malignant plasma cells. FCRL2-targeted chimeric antigen receptor T cells demonstrated antigen-specific activity in vitro and survival benefit in vivo. Together, these data provide a clinically actionable blueprint for patient stratification and precision target nomination in plasma-cell malignancies.
Metazoan development relies on the coordinated establishment of diverse gene regulatory programs that drive the formation of specific cell types, tissues and organs. The temporal and spatial control of gene expression is achieved through the concerted activity of multiple classes of cis-regulatory elements encoded in the genome. Among these, enhancers enable the establishment of specific and precise gene expression patterns and control gene expression over long linear distances, a property often referred to as distance-independent regulatory activity. However, enhancer activity is, in fact, inversely correlated with linear genomic distance, and target gene expression and transcriptional precision decrease with increasing enhancer-promoter linear distances. Here, we highlight emerging insights into multiple mechanisms that enable enhancers to precisely and robustly activate gene expression across large genomic distances. Finally, we provide a more speculative perspective on the potential advantages that long-range regulation might confer during the establishment of developmental gene expression programs.
Price et al. forge the path for simultaneous mutagenesis of any lysine(s) across all histone H3 gene copies in mammalian cells.
Repeat-expansion disorders (REDs) are a mechanistically and clinically well-defined subgroup of rare diseases caused by the expansion of short tandem repeats (STRs). These expansions can exceed several kilobases and show complex features, such as noncanonical secondary structures, somatic instability, repeat interruptions and allele-specific methylation. These characteristics are highly relevant for understanding disease mechanisms, clinical variability, prognosis and potentially therapeutic decision-making, but cannot be fully resolved using traditional diagnostic methods or short-read sequencing technologies. By contrast, long-read sequencing (LRS) enables accurate investigation of STR complexity in a single assay, facilitates the discovery of new pathogenic repeat expansions and drives advances in diagnostics, clinical and basic research, which may allow for better patient stratification in future clinical trials. This Perspective discusses recent LRS-driven discoveries, methodological and bioinformatic advances, and emerging diagnostic applications to illustrate the potential of LRS in reshaping both research and clinical practice.
Head and neck squamous cell carcinoma (HNSCC) shows substantial intra- and inter-tumoral heterogeneity. We mapped tumor architecture across HPV-positive and HPV-negative HNSCC through spatial transcriptomics (n = 26). HPV-positive tumors display hypercellularity, higher lymphocyte presence, enriched hypoxia and reduced partial epithelial-to-mesenchymal transition (p-EMT) in malignant cells. We observed two distinct spatial architectures of p-EMT: p-EMT edge, where p-EMT is coupled to fibroblasts at the invasive front via TGFβ, and p-EMT core, in which tumor-infiltrating immunosuppressive macrophages and neutrophils induce p-EMT via oncostatin M in the core of tumor nests. These two p-EMT patterns were consistent across multiple samples from the same tumor, suggesting they are tumor-wide features. Together, these findings reveal that distinct interactions in the tumor microenvironment converge on a similar p-EMT cellular phenotype, but in a different spatial pattern that may have potential biological and clinical implications for our understanding of invasion, immune modulation and new targeted therapeutics for HNSCC.
We developed a high-content screening to investigate how Alzheimer's disease (AD) genetic risk factors may affect synaptic mechanisms in rat primary neuronal cultures. Of the target genes identified, we found that Plcg2 downregulation in mouse dentate gyrus neurons consistently disrupted dendritic morphology and synaptic function. In human neuronal cultures (hNCs), PLCG2 downregulation also impaired synaptic function and increased amyloid-β (Aβ) levels and Tau phosphorylation. Very rare PLCG2 loss-of-function (LoF) variants were associated with a tenfold increased AD risk. PLCG2 LoF carriers show low mRNA/protein PLCG2/PLCγ2 levels and the R953* LoF mutation compromised synaptic function and increased AD hallmarks in hNCs. Single-nucleus RNA sequencing analyses confirmed that the downregulation of PLCG2 impacted pathways related to synaptic and neuronal functions, potentially through neurexins in neurons. In conclusion, PLCγ2 downregulation could increase AD risk by impairing synaptic functions and by increasing Aβ levels and Tau phosphorylation in neurons.
Causal disease effect sizes of proximal single-nucleotide polymorphisms (SNPs) are widely assumed to be independent but could be correlated. Here we introduce a new method, linkage disequilibrium SNP-pair effect correlation regression (LDSPEC), to estimate the correlation of causal disease effect sizes of derived alleles between proximal SNPs; LDSPEC produced robust estimates in simulations. Analyzing 70 UK Biobank diseases and traits (average N = 305,646), we detected significantly non-zero SNP-pair effect correlations (for example, -0.37 ± 0.09 for low-frequency positive linkage disequilibrium 0-100-bp SNP pairs) that decayed with distance and varied with allele frequency and linkage disequilibrium between SNPs. SNP pairs with shared functions had stronger effect correlations that spanned longer genomic distances. Consequently, SNP heritability estimates were smaller than estimates of the sum of causal effect size variances across SNPs, particularly for certain functional annotations. We recapitulated our findings via forward simulations involving stabilizing selection, implicating the action of linkage masking, whereby haplotypes containing linked SNPs with opposite effects on disease have reduced effects on fitness and escape negative selection.
Enhancer RNAs interact with promoter-derived RNAs to dictate enhancer-promoter looping, but the RNA-binding protein that mediates this process has remained unidentified. Here we identify hnRNPK as a general structural regulator that preferentially binds to nascent RNAs transcribed from enhancer and promoter regions, promoting enhancer-promoter looping and transcriptional activation. We further show that hnRNPK forms phase-separated, cavity-containing condensates that encapsulate RNA polymerase II (Pol II) via its RPB3 subunit, facilitating chromatin looping and potentially enabling recruitment of Pol II from enhancers to promoters through protein dimerization. Notably, a mutation associated with Au-Kline syndrome in hnRNPK (c.953+1dupG) alters its condensates from a liquid-like to a gel-like state, leading to developmental defects in knock-in mice. Fibroblasts derived from these mutants display reduced enhancer-promoter looping and decreased Pol II recruitment at promoters of key developmental genes. These findings suggest that hnRNPK is a structural regulator of enhancer-promoter communication and highlight the importance of RNA-RNA interactions mediated by RNA-binding proteins in transcriptional regulation.
The incidence of colorectal cancer (CRC) has risen in recent decades, with a disproportionate increase observed among younger individuals in Japan and other countries. The etiological contribution of the gut microbiota to CRC pathogenesis is recognized, yet the mechanisms involved remain to be fully clarified. Here we integrated whole-genome sequencing (WGS) and transcriptome profiling of CRC with whole-genome metagenomic sequencing of fecal samples to interrogate host-microbiome interactions at high resolution. Application of interpretable artificial intelligence enabled the stratification of CRC into four distinct microbiome-informed subtypes. WGS analysis identified mutational signatures SBS88 and ID18, linked to colibactin exposure, as early clonal events detected in 44.8% of non-hypermutated patients. Notably, these signatures were significantly more frequent among patients born after the 1960s. Microbiome-based subclassification revealed subtype-specific clinical and molecular features. Collectively, our findings indicate that colibactin exposure constitutes a prevalent and potentially modifiable risk factor for CRC in the Japanese population.
Histone post-translational modifications (hPTMs) are key regulators of chromatin states, influencing gene expression, epigenetic memory and transposable element repression across eukaryotic genomes. While many hPTMs are evolutionarily conserved, the extent to which the chromatin states they define are similarly preserved remains unclear. Here we developed a combinatorial indexing chromatin immunoprecipitation followed by sequencing method to simultaneously profile specific hPTMs across diverse eukaryotic lineages, including amoebozoans, rhizarians, discobans and cryptomonads. Our analyses revealed highly conserved euchromatin states at active gene promoters and gene bodies. In contrast, we observed diverse configurations of repressive heterochromatin states associated with silenced genes and transposable elements, characterized by various combinations of hPTMs such as H3K9me3, H3K27me3 and/or different H3K79 methylations. These findings suggest that, while core hPTMs are ancient and broadly conserved, their functional readout has diversified throughout eukaryotic evolution, shaping lineage-specific chromatin landscapes.
