Field synopsis and meta-analyses of genetic epidemiological evidence for Kashin–Beck disease, an endemic osteoarthropathy in China

Comprehensive comparative analysis of histopathology and gene expression in subchondral bone between kashin-beck disease and primary osteoarthritis

Kashin-Beck disease (KBD) is an endemic, degenerative osteoarthropathy that exhibits some similar characteristics to osteoarthritis (OA) but with different etiologies and pathogeneses. In addition to cartilage damage, microstructural changes of bone were observed in KBD. This study aimed to comparatively demonstrate the general histopathological changes, transcriptomics, and differentially expressed miRNAs of subchondral bone between KBD and OA. Tibial plateau subchondral bone samples were collected from eighteen patients with KBD and eighteen patients with OA. Histopathological changes were examined by hematoxylin-eosin (HE) staining, safranin O-fast green staining, and picrosirius red staining. RNA sequencing and miRNA array analysis were performed to screen the differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs), respectively. The subchondral bone samples of the tibial plateau of KBD and OA both showed increased thickness and sclerosis. A total of 179 DEGs and 124 DEMs were identified in subchondral bone between KBD and OA, which were involved in several vital GO terms and KEGG signaling pathways. Our results suggest that the pathological mechanisms of subchondral bone are different between KBD and OA, although they exhibit similar histopathological features. Integrated analysis revealed several genes such as ADAMTS14, SLC13A5, and CEACAM1, that may be crucial DEGs in subchondral bone between KBD and OA, suggesting that these genes could serve as potential differential diagnostic biomarkers for subchondral bone lesions in KBD and OA. These findings provide valuable information for further clarifying pathological changes in subchondral bone in KBD and OA.

Meta-analysis of Association Studies of Selenoprotein Gene Polymorphism and Kashin-Beck Disease: an Updated Systematic Review

To evaluate the association between selenoprotein gene polymorphisms and Kashin-Beck disease (KBD) susceptibility through a systematic review and updated meta-analysis. PubMed, Google Scholar, Cochrane library, and Chinese National Knowledge Infrastructure (CNKI) were electronically searched using the terms "selenoprotein" and "Kashin-Beck disease" or "KBD" with a search time from the establishment of the database to January 2021. The Newcastle-Ottawa Scale (NOS) was used for methodological quality evaluation of the included studies. Stata 14.0 software was used to pooled odds ratio (OR) and 95% confidence interval. There were a total of eight included case-control studies covering 2025 KBD patients and 1962 controls. Meta-analysis results show that the pooled odds ratios (OR) and 95% confidence intervals (CI) for DIO2 (rs225014) were 0.69 (0.52, 0.91), 0.69 (0.50, 0.96), and 0.72 (0.52, 0.99) in the allele, heterozygote, and dominant models, respectively. The OR and 95%CI for SEPS1 (-105G>A) were 2.47 (1.85, 3.29), 9.36 (4.58, 19.12), 2.17 (1.53, 3.08), and 8.60 (4.25, 17.38) in the allele, homozygote, dominant, and recessive models, respectively. In addition, the OR and 95%CI for Sep15 (rs5859) were 2.05 (1.06, 3.96) in the allele model. These results illustrate that there was a significant association between DIO2 (rs225014), SEPS1 (-105G>A), Sep15 (rs5859), and KBD. For GPX1 (rs1050450, rs1800668, rs3811699), DIO2 (rs225014, rs1352815, rs1388382), TrxR2 (rs1139793, rs5746841), GPX4 (rs713041, rs4807542), and SEPP1 (rs7579, 25191g/a), there was no significant statistical difference between the KBD and control groups (P>0.05). We conclude that the DIO2 (rs225014), SEPS1 (-105G>A), and Sep15 (rs5859) gene polymorphism are associated with susceptibility to KBD.

An integrative analysis of DNA methylation and transcriptome showed the dysfunction of MAPK pathway was involved in the damage of human chondrocyte induced by T-2 toxin

Background

T-2 toxin is thought to induce the growth plate and articular cartilage damage of Kashin-Beck disease (KBD), an endemic osteochondropathy in China. This study aims to explore the potential underlying mechanism of such toxic effects by integrating DNA methylation and gene expression profiles.

Methods

In this study, C28/I2 chondrocytes were treated with T-2 toxin (5 ng/mL) for 24 h and 72 h. Global DNA methylation level of chondrocyte was tested by Enzyme-Linked Immuno Sorbent Assay. Genome-wide DNA methylation and expression profiles were detected using Illumina Infinium HumanMethylation850 BeadChip and RNA-seq technique, respectively. Differentially methylated genes (DMGs) and differentially expressed genes (DEGs) were identified mainly for two stages including 24 h group versus Control group and 72 h group versus 24 h group. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed by Metascape. DMGs and DEGs were further validated by Sequenom MassARRAY system and quantitative real-time polymerase chain reaction.

Results

The global DNA methylation levels of chondrocytes exposed to T-2 toxin were significantly increased (P < 0.05). For 24 h group versus Control group (24 VS C), 189 DEGs and 590 DMGs were identified, and 4 of them were overlapping. For 72 h group versus 24 h group (72 VS 24), 1671 DEGs and 637 DMGs were identified, and 45 of them were overlapping. The enrichment analysis results of DMGs and DEGs both showed that MAPK was the one of the mainly involved signaling pathways in the regulation of chondrocytes after T-2 toxin exposure (DEGs: P_24VSc = 1.62 × 10^− 7; P_72VS24 = 1.20 × 10^− 7; DMGs: P_24VSc = 0.0056; P_72VS24 = 3.80 × 10^− 5).

Conclusions

The findings depicted a landscape of genomic methylation and transcriptome changes of chondrocytes after T-2 toxin exposure and suggested that dysfunction of MAPK pathway may play important roles in the chondrocytes damage induced by T-2 toxin, which could provide new clues for understanding the potential biological mechanism of KBD cartilage damage induced by T-2 toxin.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-021-00404-3.

Comparison of the toxic mechanism of T-2 toxin and deoxynivalenol on human chondrocytes by microarray and bioinformatics analysis

T-2 toxin and deoxynivalenol (DON) are two representative mycotoxins that are commonly found in cereals and agricultural products. As T-2 toxin and DON are considered the cause of Kashin-Beck disease, a special osteoarticular disease, chondrocytes would be a vital target site for these toxins. To fully understand the toxicity effects of T-2 toxin and DON on chondrocytes, the present study investigated and compared the gene expression profiles and underlying mechanisms of T-2 toxin and DON on cultured human chondrocytes by microarray and bioinformatics analysis. Normal human chondrocytes were treated with T-2 toxin at 0.01 μg/ml and DON at 1.0 μg/ml for 72 h and analyzed by microarray using Affymetrix Human Gene Chip. Comprehensive analysis, including gene ontology, pathways and gene-gene networks was performed to identify the crucial gene functions, related signal pathways and key genes. A total of 175 and 237 differentially expressed genes were identified in human chondrocytes for T-2 toxin and DON treatment, respectively. Of these, 47 had the same expression tendencies in the two groups. The protein-protein interaction network analysis showed that the 10 hub genes were different between the two groups. Our results provide a comprehensive understanding of the toxic mechanism of T-2 toxin and DON on human chondrocytes and suggest that although T-2 toxin and DON showed some similar toxic mechanisms in human chondrocytes, they also had different toxic characteristics.

Endemic Kashin-Beck disease: A food-sourced osteoarthropathy

OBJECTIVE

Kashin-Beck disease (KBD) is an endemic osteoarthropathy, which causes disability and heavy socioeconomic burdens. The preventive measures have been taken in the past few decades. However, recent KBD-epidemiological trend and comprehensive effect of its preventive measures need to be evaluated.

METHODS

By employing typical survey, cross-sectional survey, case-control study, intervention trial, and national surveillance, the present study summarizes comprehensive role of KBD-preventive measures.

RESULTS

The endemic KBD is distributed in a long and narrow area of the world. The latest epidemic began in the late 1950s and lasted until the end of 1980s. Epidemiology of the KBD was characterized by early-onset, gender equality, agricultural area, regional discrepancy, family aggregation, annual fluctuation, etc. Multivariate regression analysis suggested that etiology of the KBD was food-related factors such as fungal contamination of grains, selenium deficiency, imbalance of protein intake, etc. A series of intervention measures for KBD control had been implemented since 1990s, and involved more than 300 million residents. National incidences were 22.1% in 1990, 16.0% in 1995, 12.3% in 2000, 5.5% in 2005, 0.38% in 2010, and 0.18 in 2015, respectively. Although new patients were annually decreased, it still affected 22,567,600 inhabitants and there were 574,925 patients in 2016.

CONCLUSIONS

Etiology of the KBD is food-sourced. Its decreased incidence may attribute to an effective implementation of preventive measures. It is possible to eradicate KBD from the earth in the near future.