Search

The identification of genomic loci that are associated with osteoarthritis (OA) has provided a starting point for understanding how genetic variation activates catabolic processes in the joint. However, genetic variants can only alter gene expression and cellular function when the epigenetic environment is permissive to these effects. In this review, we provide examples of how epigenetic shifts at distinct life stages can alter the risk for OA, which we posit is critical for the proper interpretation of genome-wide association studies (GWAS). During development, intensive work on the growth and differentiation factor 5 ( GDF5) locus has revealed the importance of tissue-specific enhancer activity in controlling both joint development and the subsequent risk for OA. During homeostasis in adults, underlying genetic risk factors may help establish beneficial or catabolic “set points” that dictate tissue function, with a strong cumulative effect on OA risk. During aging, methylation changes and the reorganization of chromatin can “unmask” the effects of genetic variants. The destructive function of variants that alter aging would only mediate effects after reproductive competence and thus avoid any evolutionary selection pressure, as consistent with larger frameworks of biological aging and its relationship to disease. A similar “unmasking” may occur during OA progression, which is supported by the finding of distinct expression quantitative trait loci (eQTLs) in chondrocytes depending on the degree of tissue degradation. Finally, we propose that massively parallel reporter assays (MPRAs) will be a valuable tool to test the function of putative OA GWAS variants in chondrocytes from different life stages.

To address large gaps in our understanding of the molecular regulation of articular and growth plate cartilage development in humans, we used our directed differentiation approach to generate these distinct cartilage tissues from human embryonic stem cells. The resulting transcriptomic profiles of hESC-derived articular and growth plate chondrocytes were similar to fetal epiphyseal and growth plate chondrocytes, with respect to genes both known and previously unknown to cartilage biology. With the goal to characterize the regulatory landscapes accompanying these respective transcriptomes, we mapped chromatin accessibility in hESC-derived chondrocyte lineages, and mouse embryonic chondrocytes, using ATAC-sequencing. Integration of the expression dataset with the differentially accessible genomic regions revealed lineage-specific gene regulatory networks. We validated functional interactions of two transcription factors (TFs) (RUNX2 in growth plate chondrocytes and RELA in articular chondrocytes) with their predicted genomic targets. The maps we provide thus represent a framework for probing regulatory interactions governing chondrocyte differentiation. This work constitutes a substantial step towards comprehensive and comparative molecular characterizations of distinct chondrogenic lineages and sheds new light on human cartilage development and biology.

Individuals infected with the SARS-CoV-2 virus present with a wide variety of symptoms ranging from asymptomatic to severe and even lethal outcomes. Past research has revealed a genetic haplotype on chromosome 3 that entered the human population via introgression from Neanderthals as the strongest genetic risk factor for the severe response to COVID-19. However, the specific variants along this introgressed haplotype that contribute to this risk and the biological mechanisms that are involved remain unclear. Here, we assess the variants present on the risk haplotype for their likelihood of driving the genetic predisposition to severe COVID-19 outcomes. We do this by first exploring their impact on the regulation of genes involved in COVID-19 infection using a variety of population genetics and functional genomics tools. We then perform a locus-specific massively parallel reporter assay to individually assess the regulatory potential of each allele on the haplotype in a multipotent immune-related cell line. We ultimately reduce the set of over 600 linked genetic variants to identify four introgressed alleles that are strong functional candidates for driving the association between this locus and severe COVID-19. Using reporter assays in the presence/absence of SARS-CoV-2, we find evidence that these variants respond to viral infection. These variants likely drive the locus’ impact on severity by modulating the regulation of two critical chemokine receptor genes: CCR1 and CCR5. These alleles are ideal targets for future functional investigations into the interaction between host genomics and COVID-19 outcomes.

Epigenetic shifts are a hallmark of aging that impact transcriptional networks at regulatory level. These shifts may modify the effects of genetic regulatory variants during aging and contribute to disease pathomechanism. However, these shifts occur on the backdrop of epigenetic changes experienced throughout an individual’s development into adulthood; thus, the phenotypic, and ultimately fitness, effects of regulatory variants subject to developmental- versus aging-related epigenetic shifts may differ considerably. Natural selection therefore may act differently on variants depending on their changing epigenetic context, which we propose as a novel lens through which to consider regulatory sequence evolution and phenotypic effects. Here, we define genomic regions subjected to altered chromatin accessibility as tissues transition from their fetal to adult forms, and subsequently from early to late adulthood. Based on these epigenomic datasets, we examine patterns of evolutionary constraint and potential functional impacts of sequence variation (e.g., genetic disease risk associations). We find that while the signals observed with developmental epigenetic changes are consistent with stronger fitness consequences (i.e., negative selection pressures), they tend to have weaker effects on genetic risk associations for aging-related diseases. Conversely, we see stronger effects of variants with increased local accessibility in adult tissues, strongest in young adult when compared to old. We propose a model for how epigenetic status of a region may influence the effects of evolutionary relevant sequence variation, and suggest that such a perspective on gene regulatory networks may elucidate our understanding of aging biology.

To address large gaps in our understanding of the molecular regulation of articular and growth plate cartilage development in humans, we used our directed differentiation approach to generate these distinct cartilage tissues from human embryonic stem cells. The resulting transcriptomic profiles of hESC-derived articular and growth plate chondrocytes were similar to fetal epiphyseal and growth plate chondrocytes, with respect to genes both known and previously unknown to cartilage biology. With the goal to characterize the regulatory landscapes accompanying these respective transcriptomes, we mapped chromatin accessibility in hESC-derived chondrocyte lineages, and mouse embryonic chondrocytes, using ATAC-sequencing. Integration of the expression dataset with the differentially accessible genomic regions revealed lineage-specific gene regulatory networks. We validated functional interactions of two transcription factors (RUNX2 in growth plate chondrocytes and RELA in articular chondrocytes) with their predicted genomic targets. The maps we provide thus represent a framework for probing regulatory interactions governing chondrocyte differentiation. This work constitutes a substantial step towards comprehensive and comparative molecular characterizations of distinct chondrogenic lineages, and sheds new light on human cartilage development and biology.

BackgroundDespite recent advances in the understanding of the genetic architecture of osteoarthritis (OA), only two genetic loci have been identified for OA of the hand, in part explained by the complexity of the different hand joints and heterogeneity of OA pathology.MethodsWe used data from the Rotterdam Study (RSI, RSII and RSIII) to create three hand OA phenotypes based on clustering patterns of radiographic OA severity to increase power in our modest discovery genome-wide association studies in the RS (n=8700), and sought replication in an independent cohort, the Framingham Heart Study (n=1203). We used multiple approaches that leverage different levels of information and functional data to further investigate the underlying biological mechanisms and candidate genes for replicated loci. We also attempted to replicate known OA loci at other joint sites, including the hips and knees.ResultsWe found two novel genome-wide significant loci for OA in the thumb joints. We identifiedWNT9Aas a possible novel causal gene involved in OA pathogenesis. Furthermore, several previously identified genetic loci for OA seem to confer risk for OA across multiple joints:TGFa,RUNX2,COL27A1,ASTN2,IL11andGDF5loci.ConclusionsWe identified a robust novel genetic locus for hand OA on chromosome 1, of whichWNT9Ais the most likely causal gene. In addition, multiple genetic loci were identified to be associated with OA across multiple joints. Our study confirms the potential for novel insight into the genetic architecture of OA by using biologically meaningful stratified phenotypes.

Dentine‐ and enamel‐forming cells secrete matrix in consistent rhythmic phases, resulting in the formation of successive microscopic growth lines inside tooth crowns and roots. Experimental studies of various mammals have proven that these lines are laid down in subdaily, daily (circadian), and multidaily rhythms, but it is less clear how these rhythms are initiated and maintained. In 2001, researchers reported that lesioning the so‐called master biological clock, the suprachiasmatic nucleus (SCN), halted daily line formation in rat dentine, whereas subdaily lines persisted. More recently, a key clock gene (Bmal1) expressed in the SCN in a circadian manner was also found to be active in dentine‐ and enamel‐ secretory cells. To probe these potential neurological and local mechanisms for the production of rhythmic lines in teeth, we reexamined the role of the SCN in growth line formation in Wistar rats and investigated the presence of daily lines in Bmal1 knockout mice (Bmal1 (−/−)). In contrast to the results of the 2001 study, we found that both daily and subdaily growth lines persisted in rat dentine after complete or partial SCN lesion in the majority of individuals. In mice, after transfer into constant darkness, daily rhythms continued to manifest as incremental lines in the dentine of each Bmal1 genotype (wild‐type, Bmal (+/–), and Bmal1(−/−)). These results affirm that the manifestation of biological rhythms in teeth is a robust phenomenon, imply a more autonomous role of local biological clocks in tooth growth than previously suggested, and underscore the need further to elucidate tissue‐specific circadian biology and its role in incremental line formation. Investigations of this nature will strengthen an invaluable system for determining growth rates and calendar ages from mammalian hard tissues, as well as documenting the early lives of fossil hominins and other primates.

It remains difficult to characterize the source of pain in knee joints either using radiographs or magnetic resonance imaging (MRI). We sought to determine if advanced machine learning methods such as deep neural networks could distinguish knees with pain from those without it and identify the structural features that are associated with knee pain. We constructed a convolutional Siamese network to associate MRI scans obtained on subjects from the Osteoarthritis Initiative (OAI) with frequent unilateral knee pain comparing the knee with frequent pain to the contralateral knee without pain. The Siamese network architecture enabled pairwise learning of information from two-dimensional (2D) sagittal intermediate-weighted turbo spin echo slices obtained from similar locations on both knees. Class activation mapping (CAM) was utilized to create saliency maps, which highlighted the regions most associated with knee pain. The MRI scans and the CAMs of each subject were reviewed by an expert radiologist to identify the presence of abnormalities within the model-predicted regions of high association. Using 10-fold cross-validation, our model achieved an area under curve (AUC) value of 0.808. When individuals whose knee WOMAC pain scores were not discordant were excluded, model performance increased to 0.853. The radiologist review revealed that about 86% of the cases that were predicted correctly had effusion-synovitis within the regions that were most associated with pain. This study demonstrates a proof of principle that deep learning can be applied to assess knee pain from MRI scans. • Our article is the first to leverage a deep learning framework to associate MR images of the knee with knee pain. • We developed a convolutional Siamese network that had the ability to fuse information from multiple two-dimensional (2D) MRI slices from the knee with pain and the contralateral knee of the same individual without pain to predict unilateral knee pain. • Our model achieved an area under curve (AUC) value of 0.808. When individuals who had WOMAC pain scores that were not discordant for knees (pain discordance

A tendon’s ordered extracellular matrix (ECM) is integral for transmitting force and highly prone to injury. How tendon cells, or tenocytes, embedded within this dense ECM mobilize and contribute to healing is unknown. Here, we identify a specialized Axin2+ population in mouse and human tendons that remains latent in homeostasis yet initiates the healing response and serves as a major source of tendon progenitors. Axin2+ tenocytes readily expand in vitro and express stem cell markers. In vivo, Axin2+-descendants are major functional contributors to repair: Axin2+ tenocytes proliferate, acquire injury responsive states, and re-adopt a tenocyte fate post-injury. Specific loss of Wnt secretion in Axin2+ cells alters their progenitor identity, disrupts their activation upon injury, and severely compromises any healing response. Our work highlights an unusual paradigm, wherein quiescent Axin2+ tenocytes self-regulate their identity and mobilization upon injury and provide the key initiating signal to organize a tendon-wide healing response.

Summary Despite the great diversity of vertebrate limb proportion and our deep understanding of the genetic mechanisms that drive skeletal elongation, little is known about how individual bones reach different lengths in any species. Here, we directly compare the transcriptomes of homologous growth cartilages of the mouse (Mus musculus) and bipedal jerboa (Jaculus jaculus), the latter of which has “mouse-like” arms but extremely long metatarsals of the feet. Intersecting gene-expression differences in metatarsals and forearms of the two species revealed that about 10% of orthologous genes are associated with the disproportionately rapid elongation of neonatal jerboa feet. These include genes and enriched pathways not previously associated with endochondral elongation as well as those that might diversify skeletal proportion in addition to their known requirements for bone growth throughout the skeleton. We also identified transcription regulators that might act as “nodes” for sweeping differences in genome expression between species. Among these, Shox2, which is necessary for proximal limb elongation, has gained expression in jerboa metatarsals where it has not been detected in other vertebrates. We show that Shox2 is sufficient to increase mouse distal limb length, and a nearby putative cis-regulatory region is preferentially accessible in jerboa metatarsals. In addition to mechanisms that might directly promote growth, we found evidence that jerboa foot elongation may occur in part by de-repressing latent growth potential. The genes and pathways that we identified here provide a framework to understand the modular genetic control of skeletal growth and the remarkable malleability of vertebrate limb proportion.

Although some variation introgressed from Neanderthals has undergone selective sweeps, little is known about its functional significance. We used a Massively Parallel Reporter Assay (MPRA) to assay 5,353 high-frequency introgressed variants for their ability to modulate the gene expression within 170 bp of endogenous sequence. We identified 2,548 variants in active putative cis-regulatory elements (CREs) and 292 expression-modulating variants (emVars). These emVars are predicted to alter the binding motifs of important immune transcription factors, are enriched for associations with neutrophil and white blood cell count, and are associated with the expression of genes that function in innate immune pathways including inflammatory response and antiviral defense. We combined the MPRA data with other data sets to identify strong candidates to be driver variants of positive selection including an emVar that may contribute to protection against severe COVID-19 response. We endogenously deleted two CREs containing expression-modulation variants linked to immune function, rs11624425 and rs80317430, identifying their primary genic targets as ELMSAN1, and PAN2 and STAT2, respectively, three genes differentially expressed during influenza infection. Overall, we present the first database of experimentally identified expression-modulating Neanderthal-introgressed alleles contributing to potential immune response in modern humans.

Evolutionary responses to selection for bipedalism and childbirth have shaped the human pelvis, a structure that differs substantially from that in apes. Morphology related to these factors is present by birth, yet the developmental-genetic mechanisms governing pelvic shape remain largely unknown. Here, we pinpoint and characterize a key gestational window when human-specific pelvic morphology becomes recognizable, as the ilium and the entire pelvis acquire traits essential for human walking and birth. We next use functional genomics to molecularly characterize chondrocytes from different pelvic subelements during this window to reveal their developmental-genetic architectures. We then find notable evidence of ancient selection and genetic constraint on regulatory sequences involved in ilium expansion and growth, findings complemented by our phenotypic analyses showing that variation in iliac traits is reduced in humans compared to African apes. Our datasets provide important resources for musculoskeletal biology and begin to elucidate developmental mechanisms that shape human-specific morphology.

Individuals infected with the SARS-CoV-2 virus present with a wide variety of phenotypes ranging from asymptomatic to severe and even lethal outcomes. Past research has revealed a genetic haplotype on chromosome 3 that entered the human population via introgression from Neanderthals as the strongest genetic risk factor for the severe COVID-19 phenotype. However, the specific variants along this introgressed haplotype that contribute to this risk and the biological mechanisms that are involved remain unclear. Here, we assess the variants present on the risk haplotype for their likelihood of driving the severe COVID-19 phenotype. We do this by first exploring their impact on the regulation of genes involved in COVID-19 infection using a variety of population genetics and functional genomics tools. We then perform an locus-specific massively parallel reporter assay to individually assess the regulatory potential of each allele on the haplotype in a multipotent immune-related cell line. We ultimately reduce the set of over 600 linked genetic variants to identify 4 introgressed alleles that are strong functional candidates for driving the association between this locus and severe COVID-19. These variants likely drive the locus’ impact on severity by putatively modulating the regulation of two critical chemokine receptor genes: CCR1 and CCR5. These alleles are ideal targets for future functional investigations into the interaction between host genomics and COVID-19 outcomes.

The angiotensin-converting enzyme 2 (ACE2) protein has been highly studied as a key catalytic regulator of the renin-angiotensin system (RAS), involved in fluid homeostasis and blood pressure modulation. In addition to its important physiological role as a broadly-expressed membrane-bound protein, ACE2 serves as a cell-surface receptor for some viruses - most notably, coronaviruses such as SARS-CoV and SARS-CoV-2. Differing levels of ACE2 expression may impact viral susceptibility and subsequent changes to expression may be a pathogenic mechanism of disease risk and manifestation. Therefore, an improved understanding of howACE2expression is regulated at the genomic and transcriptional level may help us understand not only how the effects of pre-existing conditions (e.g., chronic obstructive pulmonary disease) may manifest with increased COVID-19 incidence, but also the mechanisms that regulate ACE2 levels following viral infection. Here, we initially perform bioinformatic analyses of several datasets to generate hypotheses aboutACE2gene-regulatory mechanisms in the context of immune signaling and chronic oxidative stress. We then identify putative non-coding regulatory elements withinACE2intronic regions as potential determinants ofACE2expression activity. We perform functional validation of our computational predictions in vitro via targeted CRISPR-Cas9 deletions of the identifiedACE2 cis-regulatory elements in the context of both immunological stimulation and oxidative stress conditions. We demonstrate that intronicACE2regulatory elements are responsive to both immune signaling and oxidative-stress pathways, and this contributes to our understanding of how expression of this gene may be modulated at both baseline and during immune challenge. Our work supports the further pursuit of these putative mechanisms in our understanding, prevention, and treatment of infection and disease caused by ACE2-utilizing viruses such as SARS-CoV, SARS-CoV-2, and future emerging SARS-related viruses.Author SummaryThe recent emergence of the virus SARS-CoV-2 which has caused the COVID-19 pandemic has prompted scientists to intensively study how the virus enters human host cells. This work has revealed a key protein, ACE2, that acts as a receptor permitting the virus to infect cells. Much research has focused on how the virus physically interacts with ACE2, yet little is known on how ACE2 is turned on or off in human cells at the level of the DNA molecule. Understanding this level of regulation may offer additional ways to prevent or lower viral entry into human hosts. Here, we have examined the control of theACE2gene, the DNA sequence that instructs ACE2 protein receptor formation, and we have done so in the context of immune stimulation. We have indeed identified a number of DNA on/off switches forACE2that appear responsive to immuno-logical and oxidative stress. These switches may fine-tune howACE2is turned on or off before, during, and/or after infection by SARS-CoV-2 or other related coronaviruses. Our studies help pave the way for additional functional studies on these switches, and their potential therapeutic targeting in the future.

PURPOSE OF REVIEW: The goal of the review is to provide a comprehensive overview of the current understanding of mechanisms underlying variation in human stature. RECENT FINDINGS: Human height is an anthropometric trait that varies considerably within human populations as well as across the globe. Historically, much research focus was placed on understanding the biology of growth plate chondrocytes and how modifications to core chondrocyte proliferation and differentiation pathways potentially shaped height attainment in normal as well as pathological contexts. Recently, much progress has been made to improve our understanding regarding the mechanisms underlying the normal and pathological range of height variation within as well as between human populations, and today it is understood to reflect complex interactions among a myriad of genetic, environmental, and evolutionary factors. Indeed, recent improvements in genetics (e.g., GWAS) and breakthroughs in functional genomics (e.g., whole exome sequencing, DNA methylation analysis, ATAC-sequencing, and CRISPR) have shed light on previously unknown pathways/mechanisms governing pathological and common height variation. Additionally, use of an evolutionary perspective has also revealed important mechanisms that have shaped height across the planet. SUMMARY: This review provides an overview of the current knowledge of the biological mechanisms underlying height variation by highlighting new research findings on skeletal growth control with an emphasis on previously unknown pathways/mechanisms influencing pathological and common height variation. In this context, this review also discusses how evolutionary forces likely shaped the genomic architecture of height across the globe.

Hip geometry (HG) is an important predictor of fracture. We performed a meta-analysis of GWAS studies in adults to identify genetic variants that are associated with proximal femur geometry phenotypes. We analyzed four phenotypes: 1) Femoral neck length; 2) Neck-shaft angle; 3) Femoral neck width, and 4) Femoral neck section modulus, estimated from DXA scans using algorithms of hip structure analysis. In the Discovery stage, 10 cohort studies were included in the fixed-effect meta-analysis, with up to 18,719 men and women ages 16-93 years. Association analyses were performed with ~ 2.5 million polymorphismsunder an additive model adjusted for age and body mass index; an additional analysis also adjusted for height. Replication analyses of meta-GWAS significant loci (at genome-wide significance, GWS, threshold p≤5x10-8) were performed in 7 additional cohorts in-silico. In meta-analysis not adjusting for height (combined Discovery and Replication stages), GWS associations were found on chr. 4 (in HHIP), chr. 8 (C8orf34), chr. 13 (FAM10A4 and DLEU2), and chr. 20 (inDDX27). The height-adjusted meta-analysis showed significant associations at 5p15 (IRX1 and ADAMTS16); 5q35 near FGFR4; at 12p11 (in CCDC91); 11q13 (near LRP5 and PPP6R3 (rs7102273)). Several HG signals overlapped with bone mineral density (BMD), including JAG1 on chr. 20, near TNFRSF11B (chr. 8), SOX6 and LRP5 (chr. 11). Chr. 11 SNP rs7102273 was associated with any-type fracture (p = 7.5 x 10-5). We used bone transcriptome data and discovered several significant eQTLs, including rs7102273 and PPP6R3 expression (p=0.0007), and rs6556301 (intergenic, chr.5) and PDLIM7 expression(p=0.005). In conclusion, we found associations between HG measures and several genes being part of biological pathways relevant to BMD and fractures. The results provide a defined set of genes facilitating further experimental exploration and validation to understand biological mechanisms underlying human bone geometry and etiology of bone fragility.

ObjectiveDevelop a bone shape measure that reflects the extent of cartilage loss and bone flattening in knee osteoarthritis (OA) and test it against estimates of disease severity.MethodsA fast region-based convolutional neural network was trained to crop the knee joints in sagittal dual-echo steady state MRI sequences obtained from the Osteoarthritis Initiative (OAI). Publicly available annotations of the cartilage and menisci were used as references to annotate the tibia and the femur in 61 knees. Another deep neural network (U-Net) was developed to learn these annotations. Model predictions were compared with radiologist-driven annotations on an independent test set (27 knees). The U-Net was applied to automatically extract the knee joint structures on the larger OAI dataset (9,434 knees). We defined subchondral bone length (SBL), a novel shape measure characterizing the extent of overlying cartilage and bone flattening, and examined its relationship with radiographic joint space narrowing (JSN), concurrent WOMAC pain and disability as well as subsequent partial or total knee replacement (KR). Odds ratios for each outcome were estimated using relative changes in SBL on the OAI dataset into quartiles.ResultMean SBL values for knees with JSN were consistently different from knees without JSN. Greater changes of SBL from baseline were associated with greater pain and disability. For knees with medial or lateral JSN, the odds ratios between lowest and highest quartiles corresponding to SBL changes for future KR were 5.68 (95% CI:[3.90,8.27]) and 7.19 (95% CI:[3.71,13.95]), respectively.ConclusionSBL quantified OA status based on JSN severity. It has promise as an imaging marker in predicting clinical and structural OA outcomes.

PURPOSE OF REVIEW: The ability to analyze the molecular events occurring within individual cells as opposed to populations of cells is revolutionizing our understanding of musculoskeletal tissue development and disease. Single cell studies have the great potential of identifying cellular subpopulations that work in a synchronized fashion to regenerate and repair damaged tissues during normal homeostasis. In addition, such studies can elucidate how these processes break down in disease as well as identify cellular subpopulations that drive the disease. This review highlights three emerging technologies: single cell RNA sequencing (scRNA-seq), Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq), and Cytometry by Time-Of-Flight (CyTOF) mass cytometry. RECENT FINDINGS: Technological and bioinformatic tools to analyze the transcriptome, epigenome, and proteome at the individual cell level have advanced rapidly making data collection relatively easy; however, understanding how to access and interpret the data remains a challenge for many scientists. It is, therefore, of paramount significance to educate the musculoskeletal community on how single cell technologies can be used to answer research questions and advance translation. SUMMARY: This article summarizes talks given during a workshop on “Single Cell Omics” at the 2020 annual meeting of the Orthopedic Research Society. Studies that applied scRNA-seq, ATAC-seq, and CyTOF mass cytometry to cartilage development and osteoarthritis are reviewed. This body of work shows how these cutting-edge tools can advance our understanding of the cellular heterogeneity and trajectories of lineage specification during development and disease.

The mechanical challenge of attaching elastic tendons to stiff bones is solved by the formation of a unique transitional tissue. Here, we show that murine tendon-to-bone attachment cells are bi-fated, activating a mixture of chondrocyte and tenocyte transcriptomes, under regulation of shared regulatory elements and Krüppel-like factors (KLFs) transcription factors. High-throughput bulk and single-cell RNA sequencing of humeral attachment cells revealed expression of hundreds of chondrogenic and tenogenic genes, which was validated by in situ hybridization and single-molecule ISH. ATAC sequencing showed that attachment cells share accessible intergenic chromatin areas with either tenocytes or chondrocytes. Epigenomic analysis revealed enhancer signatures for most of these regions. Transgenic mouse enhancer reporter assays verified the shared activity of some of these enhancers. Finally, integrative chromatin and motif analyses and transcriptomic data implicated KLFs as regulators of attachment cells. Indeed, blocking expression of bothKlf2andKlf4in developing limb mesenchyme impaired their differentiation.

Insulin-like growth factor 2 mRNA-binding protein 1 (IGF2BP1) plays essential roles in the proliferation of skeletal muscle satellite cells (MuSCs). Increasing evidence has shown that IGF2BP1 regulates the expression of noncoding RNAs and mRNAs. However, the related molecular network remains to be fully understood. Therefore, we performed RNA sequencing and analyzed the microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and mRNAs differentially expressed in goat MuSCs treated with IGF2BP1 overexpressing and empty vectors. A total of 36 miRNAs, 59 lncRNAs, and 44 mRNAs were differentially expressed caused by IGF2BP1. Expectedly, they were enriched in muscle development-related Rap1, PI3K-AKT, and FoxO signaling pathways. Finally, we constructed a lncRNA-miRNA-mRNA interaction network containing 30 lncRNAs, 15 miRNAs, and 34 mRNAs, in which several miRNAs, including miR-133a-3p, miR-204-5p, miR-125a-3p, miR-145-3p, and miR-423-5p, relate with cell growth and participate in muscle development. Overall, we constructed an IGF2BP1-related network, which provides new insight into the myogenic proliferation of goat.

OBJECTIVE To develop a bone shape measure that reflects the extent of cartilage loss and bone flattening in knee osteoarthritis (OA) and test it against estimates of disease severity. METHODS A fast region-based convolutional neural network was trained to crop the knee joints in sagittal dual-echo steady-state magnetic resonance imaging sequences obtained from the Osteoarthritis Initiative (OAI). Publicly available annotations of the cartilage and menisci were used as references to annotate the tibia and the femur in 61 knees. Another deep neural network (U-Net) was developed to learn these annotations. Model predictions were compared to radiologist-driven annotations on an independent test set (27 knees). The U-Net was applied to automatically extract the knee joint structures on the larger OAI data set (n = 9,434 knees). We defined subchondral bone length (SBL), a novel shape measure characterizing the extent of overlying cartilage and bone flattening, and examined its relationship with radiographic joint space narrowing (JSN), concurrent pain and disability (according to the Western Ontario and McMaster Universities Osteoarthritis Index), as well as subsequent partial or total knee replacement. Odds ratios (ORs) and 95% confidence intervals (95% CIs) for each outcome were estimated using relative changes in SBL from the OAI data set stratified into quartiles. RESULTS The mean SBL values for knees with JSN were consistently different from knees without JSN. Greater changes of SBL from baseline were associated with greater pain and disability. For knees with medial or lateral JSN, the ORs for future knee replacement between the lowest and highest quartiles corresponding to SBL changes were 5.68 (95% CI 3.90-8.27) and 7.19 (95% CI 3.71-13.95), respectively. CONCLUSION SBL quantified OA status based on JSN severity and shows promise as an imaging marker in predicting clinical and structural OA outcomes.

Type I interferons (IFNs) are our first line of defence against virus infection. Recent studies have suggested the ability of SARS-CoV-2 proteins to inhibit IFN responses. Emerging data also suggest that timing and extent of IFN production is associated with manifestation of COVID-19 severity. In spite of progress in understanding how SARS-CoV-2 activates antiviral responses, mechanistic studies into wildtype SARS-CoV-2-mediated induction and inhibition of human type I IFN responses are scarce. Here we demonstrate that SARS-CoV-2 infection induces a type I IFN response in vitro and in moderate cases of COVID-19. In vitro stimulation of type I IFN expression and signaling in human airway epithelial cells is associated with activation of canonical transcriptions factors, and SARS-CoV-2 is unable to inhibit exogenous induction of these responses. Furthermore, we show that physiological levels of IFNα detected in patients with moderate COVID-19 is sufficient to suppress SARS-CoV-2 replication in human airway cells., Graphical Abstract

The connection between different tissues is vital for the development and function of any organs and systems. In the musculoskeletal system, the attachment of elastic tendons to stiff bones poses a mechanical challenge that is solved by the formation of a transitional tissue, which allows the transfer of muscle forces to the skeleton without tearing. Here, we show that tendon-to-bone attachment cells are bi-fated, activating a mixture of chondrocyte and tenocyte transcriptomes, which is regulated by sharing regulatory elements with these cells and by Krüppel-like factors transcription factors (KLF).To uncover the molecular identity of attachment cells, we first applied high-throughput RNA sequencing to murine humeral attachment cells. The results, which were validated by in situ hybridization and single-molecule in situ hybridization, reveal that attachment cells express hundreds of chondrogenic and tenogenic genes. In search for the underlying mechanism allowing these cells to express these genes, we performed ATAC sequencing and found that attachment cells share a significant fraction of accessible intergenic chromatin areas with either tenocytes or chondrocytes. Epigenomic analysis further revealed transcriptional enhancer signatures for the majority of these regions. We then examined a subset of these regions using transgenic mouse enhancer reporter. Results verified the shared activity of some of these enhancers, supporting the possibility that the transcriptome of attachment cells is regulated by enhancers with shared activities in tenocytes or chondrocytes. Finally, integrative chromatin and motif analyses, as well as the transcriptome data, indicated that KLFs are regulators of attachment cells. Indeed, blocking the expression of Klf2 and Klf4 in the developing limb mesenchyme led to abnormal differentiation of attachment cells, establishing these factors as key regulators of the fate of these cells.In summary, our findings show how the molecular identity of bi-fated attachment cells enables the formation of the unique transitional tissue that connect tendon to bone. More broadly, we show how mixing the transcriptomes of two cell types through shared enhancers and a dedicated set of transcription factors can lead to the formation of a new cell fate that connects them.

Growth and Differentiation Factor 5 (GDF5) is a key risk locus for osteoarthritis (OA). However, little is known regarding regulation of Gdf5 expression following joint tissue damage. Here, we employed Gdf5-LacZ reporter mouse lines to assess the spatiotemporal activity of Gdf5 regulatory sequences in experimental OA following destabilisation of the medial meniscus (DMM) and after acute cartilage injury and repair. Gdf5 expression was upregulated in articular cartilage post-DMM, and was increased in human OA cartilage as determined by immunohistochemistry and microarray analysis. Gdf5 expression was also upregulated during cartilage repair in mice and was switched on in injured synovium in prospective areas of cartilage formation, where it inversely correlated with expression of the transcriptional co-factor Yes-associated protein (Yap). Indeed, overexpression of Yap suppressed Gdf5 expression in chondroprogenitors in vitro. Gdf5 expression in both mouse injury models required regulatory sequence downstream of Gdf5 coding exons. Our findings suggest that Gdf5 upregulation in articular cartilage and synovium is a generic response to knee injury that is dependent on downstream regulatory sequence and in progenitors is associated with chondrogenic specification. We propose a role for Gdf5 in tissue remodelling and repair after injury, which may partly underpin its association with OA risk.

© 2018 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals, Inc. We aimed to report the first genomewide association study (GWAS) meta-analysis of dual-energy X-ray absorptiometry (DXA)-derived hip shape, which is thought to be related to the risk of both hip osteoarthritis and hip fracture. Ten hip shape modes (HSMs) were derived by statistical shape modeling using SHAPE software, from hip DXA scans in the Avon Longitudinal Study of Parents and Children (ALSPAC; adult females), TwinsUK (mixed sex), Framingham Osteoporosis Study (FOS; mixed), Osteoporotic Fractures in Men study (MrOS), and Study of Osteoporotic Fractures (SOF; females) (total N = 15,934). Associations were adjusted for age, sex, and ancestry. Five genomewide significant (p < 5 × 10−9, adjusted for 10 independent outcomes) single-nucleotide polymorphisms (SNPs) were associated with HSM1, and three SNPs with HSM2. One SNP, in high linkage disequilibrium with rs2158915 associated with HSM1, was associated with HSM5 at genomewide significance. In a look-up of previous GWASs, three of the identified SNPs were associated with hip osteoarthritis, one with hip fracture, and five with height. Seven SNPs were within 200 kb of genes involved in endochondral bone formation, namely SOX9, PTHrP, RUNX1, NKX3-2, FGFR4, DICER1, and HHIP. The SNP adjacent to DICER1 also showed osteoblast cis-regulatory activity of GSC, in which mutations have previously been reported to cause hip dysplasia. For three of the lead SNPs, SNPs in high LD (r2 > 0.5) were identified, which intersected with open chromatin sites as detected by ATAC-seq performed on embryonic mouse proximal femora. In conclusion, we identified eight SNPs independently associated with hip shape, most of which were associated with height and/or mapped close to endochondral bone formation genes, consistent with a contribution of processes involved in limb growth to hip shape and pathological sequelae. These findings raise the possibility that genetic studies of hip shape might help in understanding potential pathways involved in hip osteoarthritis and hip fracture. © 2018 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals, Inc.

Recent studies have reported evidence suggesting that portions of contemporary human genomes introgressed from archaic hominin populations went to high frequencies due to positive selection. However, no study to date has specifically addressed the post-introgression population dynamics of these putative cases of adaptive introgression. Here, for the first time, we specifically define cases of immediate adaptive introgression (iAI), in which archaic haplotypes rose to high frequencies in humans as a result of a selective sweep that occurred shortly after the introgression event. We define these cases as distinct from instances of selection on standing introgressed variation (SI), in which an introgressed haplotype initially segregated neutrally and subsequently underwent positive selection. Using a geographically diverse dataset, we report novel cases of selection on introgressed variation in living humans and shortlisted among these cases those whose selective sweeps are more consistent with having been the product of iAI rather than SI. Many of these novel inferred iAI haplotypes have potential biological relevance, including three introgressed haplotypes that contain immune-related genes in West Siberians, South Asians, and West Eurasians. Overall, our results suggest that iAI may not represent the full picture of positive selection on archaically introgressed haplotypes in humans and that more work needs to be done to analyze the role of SI in the archaic introgression landscape of living humans.

Variants in GDF5 are associated with human arthritis and height reduction, but the causal mutations are still unknown. We surveyed GDF5 for regulatory regions in transgenic mice, and fine-mapped separate enhancers controlling expression in joints versus growing ends of long bones. A large downstream regulatory region harbors a novel growth enhancer (GROW1), which is required for normal GDF5 expression at ends of developing bones, and for normal bone lengths in vivo. Human GROW1 contains a common base pair change that reduces enhancer activity and colocalizes with peaks of positive selection in humans. The derived allele is rare in Africa, but common in Eurasia and found in Neandertals and Denisovans. Our studies suggest that an ancient regulatory variant in GROW1 has been repeatedly selected in northern environments, and that past selection on growth phenotypes explains the high frequency of a GDF5 haplotype that also increases arthritis susceptibility in many human populations., One Sentence Summary: A common regulatory variant at the GDF5 locus alters the function of a growth enhancer in Neandertals and modern humans.

Changes in cell proliferation define transitions from tissue growth to physiological homeostasis. In tendons, a highly organized extracellular matrix undergoes significant postnatal expansion to drive growth, but once formed, it appears to undergo little turnover. However, tendon cell activity during growth and homeostatic maintenance is less well defined. Using complementary methods of genetic H2B-GFP pulse-chase labeling and BrdU incorporation in mice, we show significant postnatal tendon cell proliferation, correlating with longitudinal Achilles tendon growth. Around day 21, there is a transition in cell turnover with a significant decline in proliferation. After this time, we find low amounts of homeostatic tendon cell proliferation from 3 to 20 months. These results demonstrate that tendons harbor significant postnatal mitotic activity, and limited, but detectable activity in adult and aged stages. It also points towards the possibility that the adult tendon harbors resident tendon progenitor populations, which would have important therapeutic implications., eLife digest Muscles are anchored to bones via fibrous structures called tendons, which are made from specialized cells and large proteins called collagens. Tendons in newborns have round cells that can easily divide to make new cells, allowing the tissue to grow. Compared to newborns, tendon cells in adults have a star-like shape and appear to have stopped dividing. Tendon cells in adults are less abundant than in newborns, with more of the tendon made up of collagen. Adult tendons also do not heal as well as young tendons following an injury. However, it was previously unknown when the characteristics of young tendons are lost. Additionally, it was unclear whether adult tendon cells completely stop dividing or simply do so more slowly, or if these changes in cell division are what causes adult tendons to heal less easily. To investigate this, Grinstein et al. used microscopy and cellular and molecular tools to examine tendon cell division in mice of different ages. The results indicated that tendon cells continue to divide after birth, but they do so more slowly as mice age. The researchers saw the most significant changes in the mice after they reached two weeks of age, which is also when the structure of the tendons begins to transition to that of an adult. Young mice, which are growing and learning to move, have tendons that produce new cells faster than adult tendons. This may be necessary for young tendons to grow and adapt to the strains of movement. As tendons age, most of their cells lose the ability to divide, which seems to prevent the tendons from healing fully after injury. Notably, Grinstein et al. found that some cells are still able to divide slowly in adult tendons, which may be important for healing. Further research is needed to examine the mechanisms involved in these changes and the factors that drive them, but a deeper understanding of tendon biology could lead to new therapies for treating tendon injuries.

Variation in pelvic morphology has a complex genetic basis and its patterning and specification is governed by conserved developmental pathways. Whether the mechanisms underlying the differentiation and specification of the pelvis also produce the morphological covariation on which natural selection may act, is still an open question in evolutionary developmental biology. We use high‐resolution quantitative trait locus (QTL) mapping in the F34 generation of an advanced intercross experiment (LG,SM‐G34) to characterize the genetic architecture of the mouse pelvis. We test the prediction that genomic features linked to developmental patterning and differentiation of the hind limb and pelvis and the regulation of chondrogenesis are overrepresented in QTL. We find 31 single QTL trait associations at the genome‐ or chromosome‐wise significance level coalescing to 27 pleiotropic loci. We recover further QTL at a more relaxed significance threshold replicating locations found in a previous experiment in an earlier generation of the same population. QTL were more likely than chance to harbor Pitx1 and Sox9 Class II chromatin immunoprecipitation‐seq features active during development of skeletal features. There was weak or no support for the enrichment of seven more categories of developmental features drawn from the literature. Our results suggest that genotypic variation is channeled through a subset of developmental processes involved in the generation of phenotypic variation in the pelvis. This finding indicates that the evolvability of complex traits may be subject to biases not evident from patterns of covariance among morphological features or developmental patterning when either is considered in isolation.

In tetrapods, the scapular and pelvic girdles perform the important function of anchoring the limbs to the trunk of the body and facilitating the movement of each appendage. This shared function, however, is one of relatively few similarities between the scapula and pelvis, which have significantly different morphologies, evolutionary histories, embryonic origins, and underlying genetic pathways. The scapula evolved in jawless fish prior to the pelvis, and its embryonic development is unique among bones in that it is derived from multiple progenitor cell populations, including the dermomyotome, somatopleure, and neural crest. Conversely, the pelvis evolved several million years later in jawed fish, and it develops from an embryonic somatopleuric cell population. The genetic networks controlling the formation of the pelvis and scapula also share similarities and differences, with a number of genes shaping only one or the other, while other gene products such as PBX transcription factors act as hierarchical developmental regulators of both girdle structures. Here, we provide a detailed review of the cellular processes and genetic networks underlying pelvis and scapula formation in tetrapods, while also highlighting unanswered questions about girdle evolution and development.

In bone, sclerostin is mainly osteocyte-derived and plays an important local role in adaptive responses to mechanical loading. Whether circulating levels of sclerostin also play a functional role is currently unclear, which we aimed to examine by two sample Mendelian Randomisation (MR). A genetic instrument for circulating sclerostin, derived from a genome wide association study (GWAS) meta-analysis of serum sclerostin in 10,584 European-descent individuals, was examined in relation to femoral neck bone mineral density (BMD; n= 32,744) in GEFOS, and estimated BMD by heel ultrasound (eBMD; n=426,824), and fracture risk (n=426,795), in UK Biobank. Our GWAS identified two novel serum sclerostin loci, B4GALNT3 (standard deviation (SD)) change in sclerostin per A allele (β=0.20, P=4.6×10−49), and GALNT1 (β=0.11 per G allele, P=4.4×10−11). B4GALNT3 is an N-acetyl-galactosaminyltransferase, adding a terminal LacdiNAc disaccharide to target glycocoproteins, found to be predominantly expressed in kidney, whereas GALNT1 is an enzyme causing mucin-type O-linked glycosylation. Using these two SNPs as genetic instruments, MR revealed an inverse causal relationship between serum sclerostin and femoral neck BMD (β= −0.12, 95%CI= −0.20 to −0.05) and eBMD (β= −0.12, 95%CI= −0.14 to −0.10), and a positive relationship with fracture risk (β= 0.11, 95%CI= 0.01 to 0.21). Colocalization analysis demonstrated common genetic signals within the B4GALNT3 locus for higher sclerostin, lower eBMD, and greater B4GALNT3 expression in arterial tissue (Probability>99%). Our findings suggest that higher sclerostin levels are causally related to lower BMD and greater fracture risk. Hence, strategies for reducing circulating sclerostin, for example by targeting glycosylation enzymes as suggested by our GWAS results, may prove valuable in treating osteoporosis.

Extensive variation in reproduction has arisen in Chinese Mongolian sheep during recent domestication. Hu and Small-tailed Han sheep, for example, have become non-seasonal breeders and exhibit higher fecundity than Tan and Ujumqin breeds. We therefore scanned reproduction-related single-nucleotide variations from methylated DNA-immunoprecipitation sequencing data generated from each of those four breeds to uncover potential mechanisms underlying this breed variation. We generated a high-quality map of single nucleotide variations (SNVs) in DNA methylation enriched regions, and found that the majority of variants are located within non-coding regions. We identified 359 SNVs within the Sheep Quantitative Trait Locus (QTL) database. Nineteen of these SNVs associated with the Aseasonal Reproduction QTL, and 10 out of the 19 reside close to genes with known reproduction functions. We also identified the well-known FecB mutation in high-fecundity sheep (Hu and Small-tailed Han sheep). When we applied these FecB finding to our breeding system, we improved lambing rate by 175%. In summary, this study provided strong candidate SNVs associated with sheep fecundity that can serve as targets for functional testing and to enhance selective breeding strategies. Mol. Reprod. Dev. 83: 958-967, 2016 © 2016 Wiley Periodicals, Inc.