skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on May 29, 2026

Title: Translational genomics of osteoarthritis in 1,962,069 individuals
Abstract Osteoarthritis is the third most rapidly growing health condition associated with disability, after dementia and diabetes1. By 2050, the total number of patients with osteoarthritis is estimated to reach 1 billion worldwide2. As no disease-modifying treatments exist for osteoarthritis, a better understanding of disease aetiopathology is urgently needed. Here we perform a genome-wide association study meta-analyses across up to 489,975 cases and 1,472,094 controls, establishing 962 independent associations, 513 of which have not been previously reported. Using single-cell multiomics data, we identify signal enrichment in embryonic skeletal development pathways. We integrate orthogonal lines of evidence, including transcriptome, proteome and epigenome profiles of primary joint tissues, and implicate 700 effector genes. Within these, we find rare coding-variant burden associations with effect sizes that are consistently higher than common frequency variant associations. We highlight eight biological processes in which we find convergent involvement of multiple effector genes, including the circadian clock, glial-cell-related processes and pathways with an established role in osteoarthritis (TGFβ, FGF, WNT, BMP and retinoic acid signalling, and extracellular matrix organization). We find that 10% of the effector genes express a protein that is the target of approved drugs, offering repurposing opportunities, which can accelerate translation.  more » « less
Award ID(s):
2205441 2054253
PAR ID:
10616920
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; « less
Corporate Creator(s):
; ; ; ; ; ; ; ;
Publisher / Repository:
Nature Portfolio, Springer Nature
Date Published:
Journal Name:
Nature
Volume:
641
Issue:
8065
ISSN:
0028-0836
Page Range / eLocation ID:
1217 to 1224
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract BackgroundUncovering the functional relevance underlying verbal declarative memory (VDM) genome-wide association study (GWAS) results may facilitate the development of interventions to reduce age-related memory decline and dementia. MethodsWe performed multi-omics and pathway enrichment analyses of paragraph (PAR-dr) and word list (WL-dr) delayed recall GWAS from 29,076 older non-demented individuals of European descent. We assessed the relationship between single-variant associations and expression quantitative trait loci (eQTLs) in 44 tissues and methylation quantitative trait loci (meQTLs) in the hippocampus. We determined the relationship between gene associations and transcript levels in 53 tissues, annotation as immune genes, and regulation by transcription factors (TFs) and microRNAs. To identify significant pathways, gene set enrichment was tested in each cohort and meta-analyzed across cohorts. Analyses of differential expression in brain tissues were conducted for pathway component genes. ResultsThe single-variant associations of VDM showed significant linkage disequilibrium (LD) with eQTLs across all tissues and meQTLs within the hippocampus. Stronger WL-dr gene associations correlated with reduced expression in four brain tissues, including the hippocampus. More robust PAR-dr and/or WL-dr gene associations were intricately linked with immunity and were influenced by 31 TFs and 2 microRNAs. Six pathways, including type I diabetes, exhibited significant associations with both PAR-dr and WL-dr. These pathways included fifteen MHC genes intricately linked to VDM performance, showing diverse expression patterns based on cognitive status in brain tissues. ConclusionsVDM genetic associations influence expression regulation via eQTLs and meQTLs. The involvement of TFs, microRNAs, MHC genes, and immune-related pathways contributes to VDM performance in older individuals. 
    more » « less
  2. Coxiella burnetii is an obligate intracellular bacterium that lives in a modified lysosome termed the Coxiella containing vacuole (CCV). C. burnetii is the causative agent of the zoonotic illness Q Fever, which primarily infects ruminant livestock and can spread to humans via inhalation. Acute Q fever is characterized by a flu like illness, while chronic disease is associated with more severe symptoms including endocarditis and chronic fatigue syndrome. C. burnetii has a biphasic life cycle consisting of a small cell variant (SCV) and large cell variant (LCV). The SCV is environmentally stable and can cause infection via inhalation of 1 - 10 infectious particles. Once taken up by the host cell, C. burnetii must manipulate the signaling pathways of the host cell to form the CCV where it switches to the LCV, the metabolic and replicative form of the bacterium. It primarily accomplishes this goal by utilizing a Type IV B Secretion System (T4BSS), which is unique to C. burnetii and Legionella pneumophila. The T4BSS, along with some other secretion mechanisms, secretes effector proteins into the host cell. These proteins then interfere with or modulate the host cell to recruit vacuoles, evade detection by the immune system, and prevent the host cell from initiating apoptosis. After about 6 days, the LCV will convert to SCV and then initiate host cell lysis to spread infection. This review looked at many eukaryotic cells signaling pathways and the interactions between C. burnetii and host proteins. These interactions are responsible for the modulation of host cell pathways necessary for CCV formation and C. burnetii survival. Understanding these interactions better will help with future treatments for C. burnetii infection. Further discoveries in these interactions are crucial for the future of C. burnetii research. 
    more » « less
  3. Schubert, Michael (Ed.)
    Histone variants are paralogs that replace canonical histones in nucleosomes, often imparting novel functions. However, how histone variants arise and evolve is poorly understood. Reconstruction of histone protein evolution is challenging due to large differences in evolutionary rates across gene lineages and sites. Here we used intron position data from 108 nematode genomes in combination with amino acid sequence data to find disparate evolutionary histories of the three H2A variants found inCaenorhabditis elegans: the ancient H2A.ZHTZ-1, the sperm-specific HTAS-1, and HIS-35, which differs from the canonical S-phase H2A by a single glycine-to-alanine C-terminal change. Although the H2A.ZHTZ-1protein sequence is highly conserved, its gene exhibits recurrent intron gain and loss. This pattern suggests that specific intron sequences or positions may not be important to H2A.Z functionality. For HTAS-1 and HIS-35, we find variant-specific intron positions that are conserved across species. Patterns of intron position conservation indicate that the sperm-specific variant HTAS-1 arose more recently in the ancestor of a subset ofCaenorhabditisspecies, while HIS-35 arose in the ancestor ofCaenorhabditisand its sister group, including the genusDiploscapter. HIS-35 exhibits gene retention in some descendent lineages but gene loss in others, suggesting that histone variant use or functionality can be highly flexible. Surprisingly, we find the single amino acid differentiating HIS-35 from core H2A is ancestral and common across canonicalCaenorhabditisH2A sequences. Thus, we speculate that the role of HIS-35 lies not in encoding a functionally distinct protein, but instead in enabling H2A expression across the cell cycle or in distinct tissues. This work illustrates how genes encoding such partially-redundant functions may be advantageous yet relatively replaceable over evolutionary timescales, consistent with the patchwork pattern of retention and loss of both genes. Our study shows the utility of intron positions for reconstructing evolutionary histories of gene families, particularly those undergoing idiosyncratic sequence evolution. 
    more » « less
  4. Rehfeldt, Florian (Ed.)
    Cellular response to stimulation governs tissue scale processes ranging from growth and development to maintaining tissue health and initiating disease. To determine how cells coordinate their response to such stimuli, it is necessary to simultaneously track and measure the spatiotemporal distribution of their behaviors throughout the tissue. Here, we report on a novel SpatioTemporal Response AnalysisIN Situ(STRAINS) tool that uses fluorescent micrographs, cell tracking, and machine learning to measure such behavioral distributions. STRAINS is broadly applicable to any tissue where fluorescence can be used to indicate changes in cell behavior. For illustration, we use STRAINS to simultaneously analyze the mechanotransduction response of 5000 chondrocytes—over 20 million data points—in cartilage during the 50 ms to 4 hours after the tissue was subjected to local mechanical injury, known to initiate osteoarthritis. We find that chondrocytes exhibit a range of mechanobiological responses indicating activation of distinct biochemical pathways with clear spatial patterns related to the induced local strains during impact. These results illustrate the power of this approach. 
    more » « less
  5. Genomics has grown exponentially over the last decade. Common variants are associated with physiological changes through statistical strategies such as Genome-Wide Association Studies (GWAS) and quantitative trail loci (QTL). Rare variants are associated with diseases through extensive filtering tools, including population genomics and trio-based sequencing (parents and probands). However, the genomic associations require follow-up analyses to narrow causal variants, identify genes that are influenced, and to determine the physiological changes. Large quantities of data exist that can be used to connect variants to gene changes, cell types, protein pathways, clinical phenotypes, and animal models that establish physiological genomics. This data combined with bioinformatics including evolutionary analysis, structural insights, and gene regulation can yield testable hypotheses for mechanisms of genomic variants. Molecular biology, biochemistry, cell culture, CRISPR editing, and animal models can test the hypotheses to give molecular variant mechanisms. Variant characterizations can be a significant component of educating future professionals at the undergraduate, graduate, or medical training programs through teaching the basic concepts and terminology of genetics while learning independent research hypothesis design. This article goes through the computational and experimental analysis strategies of variant characterization and provides examples of these tools applied in publications. © 2022 American Physiological Society. Compr Physiol 12:3303-3336, 2022. 
    more » « less