Cross platform analysis of transcriptomic data identifies ageing has distinct and opposite effects on tendon in males and females

Fast dual-echo estimation of apparent long T2 fraction using ultrashort echo time magnetic resonance imaging in tibialis tendons and its osteoporosis-related differences in women

Background

Tendon and bone comprise a critical interrelating unit. Bone loss, including that seen with osteopenia (OPe) or osteoporosis (OPo), may be associated with a reduction in tendon quality, though this remains incompletely investigated. Clinical magnetic resonance imaging (MRI) sequences cannot directly detect signals from tendons because of the very short T2. Clinical MRI may detect high-graded abnormalities by changes in the adjacent structures like bone. However, ultrashort echo time MRI (UTE-MRI) can capture high signals from all tendons. To determine if the long T2 fraction, as measured by a dual-echo UTE-MRI sequence, is a sensitive quantitative technique to the age- and bone-loss-related changes of the lower leg tendons.

Methods

This is a cross-sectional study conducted between January 2018 to February 2020 in the lower legs of 14 female patients with OPe [72±6 years old, body mass index (BMI) =25.8±6.2 kg/m2] and 31 female patients with OPo (73±6 years old, BMI=22.0±3.8 kg/m2), as well as 30 female subjects with normal bone (Normal, 35±18 years old, BMI =23.2±4.3 kg/m2), were imaged on a 3T clinical scanner using a dual-echo 3D Cones UTE sequence. We defined the apparent long T2 signal fraction (aFrac-LongT2) of tendons as the ratio between the signal at the second echo time (TE =2.2 ms) to the UTE signal. The average aFrac-LongT2 and the cross-sectional area were calculated for the anterior tibialis tendons (ATTs) and the posterior tibialis tendons (PTTs). The Kruskal-Wallis rank test was used to compare the differences in aFrac-LongT2 and the cross-sectional area of the tendons between the groups.

Results

The aFrac-LongT2 of the ATTs and PTTs were significantly higher in the OPo group compared with the Normal group (22.2% and 34.8% in the ATT and PTT, respectively, P<0.01). The cross-sectional area in the ATTs was significantly higher for the OPo group than in the Normal group (Normal/OPo difference was 28.7, P<0.01). Such a difference for PTTs did not reach the significance level. Mean aFrac-LongT2 and cross-sectional area in the OPe group were higher than the Normal group and lower than the OPo group. However, the differences did not show statistical significance, likely due to the higher BMI in the OPe group.

Conclusions

Dual-echo UTE-MRI is a rapid quantification technique, and aFrac-LongT2 values showed significant differences in tendons between Normal and OPo patients.

Small RNA signatures of the anterior cruciate ligament from patients with knee joint osteoarthritis

Introduction: The anterior cruciate ligament (ACL) is susceptible to degeneration, resulting in joint pain, reduced mobility, and osteoarthritis development. There is currently a paucity of knowledge on how anterior cruciate ligament degeneration and disease leads to osteoarthritis. Small non-coding RNAs (sncRNAs), such as microRNAs and small nucleolar RNA (snoRNA), have diverse roles, including regulation of gene expression. Methods: We profiled the sncRNAs of diseased osteoarthritic ACLs to provide novel insights into osteoarthritis development. Small RNA sequencing from the ACLs of non- or end-stage human osteoarthritic knee joints was performed. Significantly differentially expressed sncRNAs were defined, and bioinformatics analysis was undertaken. Results and Discussion: A total of 184 sncRNAs were differentially expressed: 68 small nucleolar RNAs, 26 small nuclear RNAs (snRNAs), and 90 microRNAs. We identified both novel and recognized (miR-206, -365, and -29b and -29c) osteoarthritis-related microRNAs and other sncRNAs (including SNORD72, SNORD113, and SNORD114). Significant pathway enrichment of differentially expressed miRNAs includes differentiation of the muscle, inflammation, proliferation of chondrocytes, and fibrosis. Putative mRNAs of the microRNA target genes were associated with the canonical pathways "hepatic fibrosis signaling" and "osteoarthritis." The establishing sncRNA signatures of ACL disease during osteoarthritis could serve as novel biomarkers and potential therapeutic targets in ACL degeneration and osteoarthritis development.

Mass Spectrometry-Based Atlas of Extracellular Matrix Proteins across 25 Mouse Organs

The extracellular matrix (ECM) is a critical non-cellular component of multicellular organisms containing a variety of proteins, glycoproteins, and proteoglycans which have been implicated in a wide variety of essential biological processes, including development, wound healing, and aging. Due to low solubility, many ECM proteins have been underrepresented in previous proteomic datasets. Using an optimized three-step decellularization and ECM extraction method involving chaotrope extraction and digestion via hydroxylamine hydrochloride, we have generated coverage of the matrisome across 25 organs. We observe that the top 100 most abundant proteins from the ECM fractions of all tissues are generally present in all tissues, indicating that tissue matrices are principally composed of a shared set of ECM proteins. However, these proteins vary up to 4000-fold between tissues, resulting in highly unique matrix profiles even with the same primary set of proteins. A data reduction approach was used to reveal related networks of expressed ECM proteins across varying tissues, including basement membrane and collagen subtypes.

Research perspectives-Pipelines to human tendon transcriptomics

Tendon transcriptomics is a rapidly growing field in musculoskeletal biology. The ultimate aim of many current tendon transcriptomic studies is characterization of in vitro, ex vivo, or in vivo, healthy, and diseased tendon microenvironments to identify the underlying pathways driving human tendon pathology. The transcriptome interfaces between genomic, proteomic, and metabolomic signatures of the tendon cellular niche and the response of this niche to stimuli. Some of the greatest bottlenecks in tendon transcriptomics relate to the availability and quality of human tendon tissue, hence animal tissues are frequently used even though human tissue is most translationally relevant. Here, we review the variability associated with human donor and procurement factors, such as whether the tendon is cadaveric or a clinical remnant, and how these variables affect the quality and relevance of the transcriptomes obtained. Moreover, age, sex, and health demographic variables impact the human tendon transcriptome. Tendons present tissue-specific challenges for cell, nuclei, and RNA extraction that include a dense extracellular matrix, low cellularity, and therefore low RNA yield of variable quality. Consideration of these factors is particularly important for single-cell and single-nuclei resolution transcriptomics due to the necessity for unbiased and representative cell or nuclei populations. Different cell, nuclei, and RNA extraction methods, library preparation, and quality control methods are used by the tendon research community and attention should be paid to these when designing and reporting studies. We discuss the different components and challenges of human tendon transcriptomics, and propose pipelines, quality control, and reporting guidelines for future work in the field.

Lower Macromolecular Content in Tendons of Female Patients with Osteoporosis versus Patients with Osteopenia Detected by Ultrashort Echo Time (UTE) MRI

Tendons and bones comprise a special interacting unit where mechanical, biochemical, and metabolic interplays are continuously in effect. Bone loss in osteoporosis (OPo) and its earlier stage disease, osteopenia (OPe), may be coupled with a reduction in tendon quality. Noninvasive means for quantitatively evaluating tendon quality during disease progression may be critically important for the improvement of characterization and treatment optimization in patients with bone mineral density disorders. Though clinical magnetic resonance imaging (MRI) sequences are not typically capable of directly visualizing tendons, ultrashort echo time MRI (UTE-MRI) is able to acquire a high signal from tendons. Magnetization transfer (MT) modeling combined with UTE-MRI (i.e., UTE-MT-modeling) can indirectly assess macromolecular proton content in tendons. This study aimed to determine whether UTE-MT-modeling could detect differences in tendon quality across a spectrum of bone health. The lower legs of 14 OPe (72 ± 6 years) and 31 OPo (73 ± 6 years) female patients, as well as 30 female participants with normal bone (Normal-Bone, 36 ± 19 years), are imaged using UTE sequences on a 3T MRI scanner. Institutional review board approval is obtained for the study, and all recruited subjects provided written informed consent. A T1 measurement and UTE-MT-modeling are performed on the anterior tibialis tendon (ATT), posterior tibialis tendon (PTT), and the proximal Achilles tendon (PAT) of all subjects. The macromolecular fraction (MMF) is estimated as the main measure from UTE-MT-modeling. The mean MMF in all the investigated tendons was significantly lower in OPo patients compared with the Normal-Bone cohort (mean difference of 24.2%, p < 0.01), with the largest Normal-Bone vs. OPo difference observed in the ATT (mean difference of 32.1%, p < 0.01). Average MMF values of all the studied tendons are significantly lower in the OPo cohort compared with the OPe cohort (mean difference 16.8%, p = 0.02). Only the PPT shows significantly higher T1 values in OPo patients compared with the Normal-Bone cohort (mean difference 17.6%, p < 0.01). Considering the differences between OPo and OPe groups with similar age ranges, tendon deterioration associated with declining bone health was found to be larger than a priori detected differences caused purely by aging, highlighting UTE-MT MRI techniques as useful methods in assessing tendon quality over the course of progressive bone weakening.

Tendon and multiomics: advantages, advances, and opportunities

Tendons heal by fibrosis, which hinders function and increases re-injury risk. Yet the biology that leads to degeneration and regeneration of tendons is not completely understood. Improved understanding of the metabolic nuances that cause diverse outcomes in tendinopathies is required to solve these problems. 'Omics methods are increasingly used to characterize phenotypes in tissues. Multiomics integrates 'omic datasets to identify coherent relationships and provide insight into differences in molecular and metabolic pathways between anatomic locations, and disease stages. This work reviews the current literature pertaining to multiomics in tendon and the potential of these platforms to improve tendon regeneration. We assessed the literature and identified areas where 'omics platforms contribute to the field: (1) Tendon biology where their hierarchical complexity and demographic factors are studied. (2) Tendon degeneration and healing, where comparisons across tendon pathologies are analyzed. (3) The in vitro engineered tendon phenotype, where we compare the engineered phenotype to relevant native tissues. (4) Finally, we review regenerative and therapeutic approaches. We identified gaps in current knowledge and opportunities for future study: (1) The need to increase the diversity of human subjects and cell sources. (2) Opportunities to improve understanding of tendon heterogeneity. (3) The need to use these improvements to inform new engineered and regenerative therapeutic approaches. (4) The need to increase understanding of the development of tendon pathology. Together, the expanding use of various 'omics platforms and data analysis resulting from these platforms could substantially contribute to major advances in the tendon tissue engineering and regenerative medicine field.

Long noncoding RNA H19 accelerates tenogenic differentiation by modulating miR-140-5p/VEGFA signaling

Rotator cuff tear (RCT) is a common tendon injury, but the mechanisms of tendon healing remain incompletely understood. Elucidating the molecular mechanisms of tenogenic differentiation is essential to develop novel therapeutic strategies in clinical treatment of RCT. The long noncoding RNA H19 plays a regulatory role in tenogenic differentiation and tendon healing, but its detailed mechanism of action remains unknown. To elucidate the role of H19 in tenogenic differentiation and tendon healing, tendon-derived stem cells were harvested from the Achilles tendons of Sprague Dawley rats and a rat model of cuff tear was established for the exploration of the function of H19 in promoting tenogenic differentiation. The results showed that H19 overexpression promoted, while H19 silencing suppressed, tenogenic differentiation of tendon-derived stem cells (TDSCs). Furthermore, bioinformatic analyses and a luciferase reporter gene assay showed that H19 directly targeted and inhibited miR-140-5p to promote tenogenic differentiation. Further, inhibiting miR-140-5p directly increased VEGFA expression, revealing a novel regulatory axis between H19, miR-140-5p, and VEGFA in modulating tenogenic differentiation. In rats with RTC, implantation of H19-overexpressing TDSCs at the lesion promoted tendon healing and functional recovery. In general, the data suggest that H19 promotes tenogenic differentiation and tendon-bone healing by targeting miR-140-5p and increasing VEGFA levels. Modulation of the H19/miR-140-5p/VEGFA axis in TDSCs is a new potential strategy for clinical treatment of tendon injury.

Structure-function specialisation of the interfascicular matrix in the human achilles tendon

Tendon consists of highly aligned collagen-rich fascicles surrounded by interfascicular matrix (IFM). Some tendons act as energy stores to improve locomotion efficiency, but such tendons commonly obtain debilitating injuries. In equine tendons, energy storing is achieved primarily through specialisation of the IFM. However, no studies have investigated IFM structure-function specialisation in human tendons. Here, we compare the human positional anterior tibial tendon and energy storing Achilles tendons, testing the hypothesis that the Achilles tendon IFM has specialised composition and mechanical properties, which are lost with ageing. Data demonstrate IFM specialisation in the energy storing Achilles, with greater elasticity and fatigue resistance than in the positional anterior tibial tendon. With ageing, alterations occur predominantly to the proteome of the Achilles IFM, which are likely responsible for the observed trends towards decreased fatigue resistance. Knowledge of these key energy storing specialisations and their changes with ageing offers crucial insight towards developing treatments for tendinopathy.

Statement of significance

Developing effective therapeutics or preventative measures for tendon injury necessitates the understanding of healthy tendon function and mechanics. By establishing structure-function relationships in human tendon and determining how these are affected by ageing, potential targets for therapeutics can be identified. In this study, we have used a combination of mechanical testing, immunolabelling and proteomics analysis to study structure-function specialisations in human tendon. We demonstrate that the interfascicular matrix is specialised for energy storing in the Achilles tendon, and that its proteome is altered with ageing, which is likely responsible for the observed trends towards decreased fatigue resistance. Knowledge of these key energy storing specialisations and their changes with ageing offers crucial insight towards developing treatments and preventative approaches for tendinopathy.

Gene expression analysis of combined RNA-seq experiments using a receiver operating characteristic calibrated procedure

Because of rapid advancements in sequencing technology, the experimental platforms of RNA-seq are updated frequently. It is quite common to combine data sets from several experimental platforms for analysis in order to increase the sample size and achieve more powerful tests for detecting the presence of differential gene expression. The data sets combined from different experimental platforms will have a complex data distribution, which causes a major problem in statistical modeling as well as in multiple testing. Although plenty of research have studied this problem by modeling the batch effects, there are no general and robust data-driven procedures for RNA-seq analysis. In this paper we propose a new robust procedure which combines the use of popular methods (packages) with a data-driven simulation (DDS). We construct the average receiver operating characteristic curve through the DDS to provide the calibrated levels of significance for multiple testing. Instead of further modifying the adjusted p-values, we calibrated the levels of significance for each specific method and mean effect model. The procedure was demonstrated with several popular RNA-seq analysis methods (edgeR, DEseq2, limma+voom). The proposed procedure relaxes the stringent assumptions of data distributions for RNA-seq analysis methods and is illustrated using colorectal cancer studies from The Cancer Genome Atlas database.

Epigenetic mechanisms in Tendon Ageing

INTRODUCTION

Tendon is a composite material with a well-ordered hierarchical structure exhibiting viscoelastic properties designed to transfer force. It is recognized that the incidence of tendon injury increases with age, suggesting a deterioration in homeostatic mechanisms or reparative processes. This review summarizes epigenetic mechanisms identified in ageing healthy tendon.

SOURCES OF DATA

We searched multiple databases to produce a systematic review on the role of epigenetic mechanisms in tendon ageing.

AREAS OF AGREEMENT

Epigenetic mechanisms are important in predisposing ageing tendon to injury.

AREAS OF CONTROVERSY

The relative importance of epigenetic mechanisms are unknown in terms of promoting healthy ageing. It is also unknown whether these changes represent protective mechanisms to function or predispose to pathology.

GROWING POINT

Epigenetic markers in ageing tendon, which are under-researched including genome-wide chromatin accessibility, should be investigated.

AREAS TIMELY FOR DEVELOPING RESEARCH

Metanalysis through integration of multiple datasets and platforms will enable a holistic understanding of the epigenome in ageing and its relevance to disease.

Spectrum of Tendon Pathologies: Triggers, Trails and End-State

The biggest compartment of the musculoskeletal system is the tendons and ligaments. In particular, tendons are dense tissues connecting muscle to bone that are critical for the integrity, function and locomotion of this system. Due to the increasing age of our society and the overall rise in engagement in extreme and overuse sports, there is a growing prevalence of tendinopathies. Despite the recent advances in tendon research and due to difficult early diagnosis, a multitude of risk factors and vague understanding of the underlying biological mechanisms involved in the progression of tendon injuries, the toolbox of treatment strategies remains limited and non-satisfactory. This review is designed to summarize the current knowledge of triggers, trails and end state of tendinopathies.

Age-associated changes in the response of tendon explants to stress deprivation is sex-dependent

Purpose of the Study: The incidence of tendon injuries increases dramatically with age, which presents a major clinical burden. While previous studies have sought to identify age-related changes in extracellular matrix structure and function, few have been able to explain fully why aged tissues are more prone to degeneration and injury. In addition, recent studies have also demonstrated that age-related processes in humans may be sex-dependent, which could be responsible for muddled conclusions in changes with age. In this study, we investigate short-term responses through an ex vivo explant culture model of stress deprivation that specifically questions how age and sex differentially affect the ability of tendons to respond to altered mechanical stimulus.Materials and Methods: We subjected murine flexor explants from young (4 months of age) and aged (22-24 months of age) male and female mice to stress-deprived culture conditions for up to 1 week and investigated changes in viability, cell metabolism and proliferation, matrix biosynthesis and composition, gene expression, and inflammatory responses throughout the culture period.Results and Conclusions: We found that aging did have a significant influence on the response to stress deprivation, demonstrating that aged explants have a less robust response overall with reduced metabolic activity, viability, proliferation, and biosynthesis. However, age-related changes appeared to be sex-dependent. Together, this work demonstrates that the aging process and the subsequent effect of age on the ability of tendons to respond to stress-deprivation are inherently different based on sex, where male explants favor increased activity, apoptosis, and matrix remodeling while female explants favor reduced activity and tissue preservation.

Age-related decrease in collagen proton fraction in tibial tendons estimated by magnetization transfer modeling of ultrashort echo time magnetic resonance imaging (UTE-MRI)

Clinical magnetic resonance imaging (MRI) sequences are not often capable of directly visualizing tendons. Ultrashort echo time (UTE) MRI can acquire high signal from tendons thus enabling quantitative assessments. Magnetization transfer (MT) modeling combined with UTE-MRI-UTE-MT-modeling-can indirectly assess macromolecular protons in the tendon. This study aimed to determine if UTE-MT-modeling is a quantitative technique sensitive to the age-related changes of tendons. The legs of 26 young healthy (29 ± 6 years old) and 22 elderly (75 ± 8 years old) female subjects were imaged using UTE sequences on a 3T MRI scanner. Institutional review board approval was obtained, and all recruited subjects provided written informed consent. T1 and UTE-MT-modeling were performed on anterior tibialis tendons (ATT) and posterior tibialis tendons (PTT) as two representative human leg tendons. A series of MT pulse saturation powers (500-1500°) and frequency offsets (2-50 kHz) were used to measure the macromolecular fraction (MMF) and macromolecular T2 (T2MM). All measurements were repeated by three independent readers for a reproducibility study. MMF demonstrated significantly lower values on average in the elderly cohort compared with the younger cohort for both ATT (decreased by 16.8%, p = 0.03) and PTT (decreased by 23.0%, p < 0.01). T2MM and T1 did not show a significant nor a consistent difference between the young and elderly cohorts. For all MRI parameters, intraclass correlation coefficient (ICC) was higher than 0.98, indicating excellent consistency between measurements performed by independent readers. MMF serving as a surrogate measure for collagen content, showed a significant decrease in elderly leg tendons. This study highlighted UTE-MRI-MT techniques as a useful quantitative method to assess the impact of aging on human tendons.

Retinoic Acid Promotes Retinoic Acid Signaling by Suppression of Pitx1 In Tendon Cells: A Possible Mechanism of a Clubfoot-Like Phenotype Induced by Retinoic Acid

Background

The pathogenesis of idiopathic congenital clubfoot (CCF) is unknown. Although some familial patients have Pitx1 mutations, and the Pitx1+/− genotype causes a clubfoot-like phenotype in mice, the mechanism of Pitx1-induced CCF is unknown.

Material/Methods

We used tibialis anterior tendon samples to detect the expression of Pitx1 in idiopathic and neurogenic clubfoot patients. After obtaining Sprague-Dawley (SD) rat Achilles tendon cells, the expression of Pitx1 was knocked down by SiRNA. After 48 h of culture, mass spectrometry was used to quantitatively analyze proteins. Then, Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were used to assess the downstream pathway of PITX1. The relationship between Pitx1 and the promoter region of deacetylase 1 (Sirtuin-1 and Sirt1) was examined by luciferase and ChIP assays.

Results

We found that Pitx1 expression in the tendon samples of idiopathic CCF patients was downregulated. Mass spectrometry analysis revealed that the inhibition of Pitx1 induced the downregulation of Sirt1 expression in tendon cells. Luciferase and ChIP assays confirmed that Pitx1 binds to the promoter region of SIRT1 and promotes Sirt1 gene transcription. Further results showed that, after the inhibition of Pitx1 in tendon cells, CRABP2 acetylation increased, the nuclear import of CRABP2 was enhanced, and the expression of RARβ2 increased. After the inhibition of Pitx1, RARβ2 expression was further increased by RA treatment in tendon cells. In the presence of retinoic acid, the expression of Pitx1 was inhibited in tendon cells.

Conclusions

Pitx1 binds to the promoter region of SIRT1 and promotes the transcription of SIRT1. Positive feedback occurs between RA signaling and Pitx1.