Role of selenoprotein S (SEPS1) -105G>A polymorphisms and PI3K/Akt signaling pathway in Kashin-Beck disease

Modulation of Ras signaling pathway by exosome miRNAs in T-2 toxin-induced chondrocyte injury

This study aims to investigate the impact of T-2 toxin on the regulation of downstream target genes and signaling pathways through exosome-released miRNA in the development of cartilage damage in Kashin-Beck disease (KBD). Serum samples from KBD patients and supernatant from C28/I2 cells treated with T-2 toxin were collected for the purpose of comparing the differential expression of exosomal miRNA using absolute quantitative miRNA-seq. Target genes of differential exosomal miRNAs were identified using Targetscan and Miranda databases, followed by GO and KEGG enrichment analyses. Validation of key indicators of chondrocyte injury in KBD was conducted using Real-time quantitative PCR (RT-qPCR) and Immunohistochemical staining (IHC). A total of 20 exosomal miRNAs related to KBD were identified in serum, and 13 in chondrocytes (C28/I2). The identified exosomal miRNAs targeted 48,459 and 60,612 genes, primarily enriched in cell organelles and membranes, cell differentiation, and cytoskeleton in the serum, and the cytoplasm and nucleus, metal ion binding in chondrocyte (C28/I2). The results of the KEGG enrichment analysis indicated that the Ras signaling pathway may play a crucial role in the pathogenesis of KBD. Specifically, the upregulation of hsa-miR-181a-5p and hsa-miR-21-3p, along with the downregulation of hsa-miR-152-3p and hsa-miR-186-5p, were observed. Additionally, T-2 toxin intervention led to a significant downregulation of RALA, REL, and MAPK10 expression. Furthermore, the protein levels of RALA, REL, and MAPK10 were notably decreased in the superficial and middle layers of cartilage tissues from KBD. The induction of differential expression of chondrocyte exosomal miRNAs by T-2 toxin results in the collective regulation of target genes RALA, REL, and MAPK10, ultimately mediating the Ras signaling pathway and causing a disruption in chondrocyte extracellular matrix metabolism, leading to chondrocyte injury.

"Alphabet" Selenoproteins: Their Characteristics and Physiological Roles

Selenium (Se) is a metalloid that is recognized as one of the vital trace elements in our body and plays multiple biological roles, largely mediated by proteins containing selenium-selenoproteins. Selenoproteins mainly have oxidoreductase functions but are also involved in many different molecular signaling pathways, physiological roles, and complex pathogenic processes (including, for example, teratogenesis, neurodegenerative, immuno-inflammatory, and obesity development). All of the selenoproteins contain one selenocysteine (Sec) residue, with only one notable exception, the selenoprotein P (SELENOP), which has 10 Sec residues. Although these mechanisms have been studied intensely and in detail, the characteristics and functions of many selenoproteins remain unknown. This review is dedicated to the recent data describing the identity and the functions of several selenoproteins that are less known than glutathione peroxidases (Gpxs), iodothyronine deiodinases (DIO), thioredoxin reductases (TRxRs), and methionine sulfoxide reductases (Msrs) and which are named after alphabetical letters (i.e., F, H, I, K, M, N, O, P, R, S, T, V, W). These "alphabet" selenoproteins are involved in a wide range of physiological and pathogenetic processes such as antioxidant defense, anti-inflammation, anti-apoptosis, regulation of immune response, regulation of oxidative stress, endoplasmic reticulum (ER) stress, immune and inflammatory response, and toxin antagonism. In selenium deficiency, the "alphabet" selenoproteins are affected hierarchically, both with respect to the particular selenoprotein and the tissue of expression, as the brain or endocrine glands are hardly affected by Se deficiency due to their equipment with LRP2 or LRP8.

Nano-elemental selenium particle developed via supramolecular self-assembly of chondroitin sulfate A and Na2SeO3 to repair cartilage lesions

Cartilage repair is a significant clinical issue due to its restricted ability to regenerate and self-heal after cartilage lesions or degenerative disease. Herein, a nano-elemental selenium particle (chondroitin sulfate A‑selenium nanoparticle, CSA-SeNP) is developed by the supramolecular self-assembly of Na₂SeO₃ and negatively charged chondroitin sulfate A (CSA) via electrostatic interactions or hydrogen bonds followed by in-situ reducing of l-ascorbic acid for cartilage lesions repair. The constructed micelle exhibits a hydrodynamic particle size of 171.50 ± 2.40 nm and an exceptionally high selenium loading capacity (9.05 ± 0.03 %) and can promote chondrocyte proliferation, increase cartilage thickness, and improve the ultrastructure of chondrocytes and organelles. It mainly enhances the sulfation modification of chondroitin sulfate by up-regulating the expression of chondroitin sulfate 4-O sulfotransferase-1, -2, -3, which in turn promotes the expression of aggrecan to repair articular and epiphyseal-plate cartilage lesions. The micelles combine the bio-activity of CSA with selenium nanoparticles (SeNPs), which are less toxic than Na₂SeO₃, and low doses of CSA-SeNP are even superior to inorganic selenium in repairing cartilage lesions in rats. Thus, the developed CSA-SeNP is anticipated to be a promising selenium supplementation preparation in clinical application to address the difficulty of healing cartilage lesions with outstanding repair effects.

Effects of Selenoprotein S Knockdown on Endoplasmic Reticulum Stress in ATDC5 Cells and Gene Expression Profiles in Hypertrophic Chondrocytes

Selenoprotein S (SelS), a member of the selenoprotein family, is mainly located on the endoplasmic reticulum (ER) membrane. SelS is involved in a variety of biological processes, including oxidative stress, inflammation, glucose metabolism regulation, and ER-associated protein degradation (ERAD). This study was designed to explore the role of SelS in chondrocytes. It was confirmed that SelS is a Se-sensitive selenoprotein in low-selenium rat and cell models. ER stress was not induced in SelS knockdown ATDC5 cells. However, treatment of ATDC5 cells with tunicamycin (Tm), an ER stress inducer, increased the expression of SelS, and knockdown of SelS aggravated ER stress induced by Tm, suggesting that SelS is a regulatory molecule involved in ER stress in chondrocytes. Both osteoarthritis and Kashin-Beck disease are osteochondral diseases associated with hypertrophic chondrocyte abnormalities. Therefore, ATDC5 cells were induced to hypertrophic chondrocytes. SelS was knocked down and RNA sequencing was performed. Bioinformatics analysis of the differentially expressed genes (DEGs) revealed that SelS knockdown affected a variety of biological processes, including cell adhesion, osteoclast differentiation, and extracellular matrix homeostasis. Collectively, this study verified that SelS is sensitive to selenium levels and is an ER stress-responsive molecule. Knocking down SelS can cause abnormal expression of adhesion molecules and matrix homeostasis disorder in hypertrophic chondrocytes.

Selenoprotein S: A versatile disordered protein

Selenoprotein S (selenos) is a small, intrinsically disordered membrane protein that is associated with various cellular functions, such as inflammatory processes, cellular stress response, protein quality control, and signaling pathways. It is primarily known for its contribution to the ER-associated degradation (ERAD) pathway, which governs the extraction of misfolded proteins or misassembled protein complexes from the ER to the cytosol for degradation by the proteasome. However, selenos's other cellular roles in signaling are equally vital, including the control of transcription factors and cytokine levels. Consequently, genetic polymorphisms of selenos are associated with increased risk for diabetes, dyslipidemia, and cardiovascular diseases, while high expression levels correlate with poor prognosis in several cancers. Its inhibitory role in cytokine secretion is also exploited by viruses. Since selenos binds multiple protein complexes, however, its specific contributions to various cellular pathways and diseases have been difficult to establish. Thus, the precise cellular functions of selenos and their interconnectivity have only recently begun to emerge. This review aims to summarize recent insights into the structure, interactome, and cellular roles of selenos.

Progress of Selenium Deficiency in the Pathogenesis of Arthropathies and Selenium Supplement for Their Treatment

Selenium, an essential trace element for human health, exerts an indispensable effect in maintaining physiological homeostasis and functions in the body. Selenium deficiency is associated with arthropathies, such as Kashin-Beck disease, rheumatoid arthritis, osteoarthritis, and osteoporosis. Selenium deficiency mainly affects the normal physiological state of bone and cartilage through oxidative stress reaction and immune reaction. This review aims to explore the role of selenium deficiency and its mechanisms existed in the pathogenesis of arthropathies. Meanwhile, this review also summarized various experiments to highlight the crucial functions of selenium in maintaining the homeostasis of bone and cartilage.

Serum proteomic analysis of differentially expressed proteins and pathways involved in the mechanism of endemic osteoarthritis

Kashin-Beck disease (KBD) is a chronic and endemic osteochondral disease and the etiology and pathogenic mechanism of KBD are still unknown. This study aimed to elucidate and screen KBD-associated proteins, which were differentially expressed between KBD patients and healthy controls. We combined protein fractionation and liquid chromatography-tandem mass spectrometry (LC-MS/MS) with a high-resolution mass spectrometer coupled with tandem mass tags (TMTs) to quantitatively analyze and screen KBD-associated proteins, which were differentially expressed between KBD patients and healthy controls. In addition, we used parallel reaction monitoring (PRM) to quantify proteins in serum from patients with KBD and healthy controls in order to verify the differentially expressed proteins in patients with KBD. We identified 224 differentially expressed proteins, including 11 up-regulated and 213 down-regulated proteins. Catalase (CAT) was observed to be significantly elevated in patients with KBD compared with control patients. Further, the fold difference of CAT is significantly elevated in PRM compared with label-free quantification. The results in this study suggest that CAT may be the reflection of the dynamic nature of KBD and could be considered as a novel pathogenic indicator for patients with KBD.

Dysregulation of Transcription Profile of Selenoprotein in Patients with Kashin-Beck Disease and Its Effect on Se Deficiency-Induced Chondrocyte Apoptosis

Kashin-Beck disease (KBD) is a chronic, degenerative osteoarthropathy related to selenium (Se) deficiency. Se participates in the synthesis of selenoprotein in the form of selenocysteine. In total, 25 selenoproteins, encoded by 25 genes, are currently found in humans; however, the effects of selenoprotein genes on chondrocyte apoptosis, particularly in apoptosis-related genes, remain poorly elucidated. Therefore, in the current study, the expression of selenoprotein genes and apoptosis-related genes were determined by RT-qPCR in patients and chondrocytes and the correlations between them were analyzed using Pearson and Spearman's rank correlation, and the chondrocyte apoptosis rate was detected by Annexin V-FITC/PI. The results showed that the mRNA levels of 17 selenoprotein genes were downregulated, whereas two genes were upregulated in patients with KBD. The BAX/BCL2 ratio and the mRNA levels of BAX and P53 were increased, but the mRNA levels of BCL2 and NF-κB p65 were decreased in patients with KBD. The mRNA levels of GPX2, GPX3, DIO1, TXNRD1, TXNRD3, and SPS2 were most closely associated with apoptosis-related genes in patients with KBD. Moreover, in the Se deficiency group, the mRNA levels of GPX3, DIO1, and TXNRD1 were downregulated and GPX activity was decreased, but the late apoptosis rate, the mRNA levels of BAX and P53, and the BAX/BCL2 ratio were increased; the opposite trend was observed in the Se supplement group. Collectively, these results indicate that selenoprotein transcription profile is dysregulated in patients with KBD. Furthermore, the expression of GPX3, DIO1, and TXNRD1 genes might be involved in the development of chondrocyte apoptosis by affecting antioxidant capacity.

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.

Patients with Osteoarthritis and Kashin-Beck Disease Display Distinct CpG Methylation Profiles in the DIO2, GPX3, and TXRND1 Promoter Regions

OBJECTIVE

We aimed to analyze deoxycytidine-deoxyguanosine dinucleotide (CpGs) methylation profiles in DIO2, GPX3, and TXNRD1 promoter regions in osteoarthritis (OA) and Kashin-Beck disease (KBD) patients.

METHODS

Blood samples were collected from 16 primary OA patients and corresponding 16 healthy individuals and analyzed for methylations in the CpGs of DIO2, GPX3, and TXNRD1 promoter regions using MALDI-TOF-MS. The methylation profiles of these regions were then compared between OA and KBD patients.

RESULTS

DIO2-1_CpG_2 and DIO2-1_CpG_3 methylations were significantly lower in OA than KBD patients (P < 0.05). A similar trend was observed for GPX3-1_CpG_4, GPX3-1_CpG_7, GPX3-1_CpG_8.9.10, GPX3-1_CpG_13.14.15 and GPX3-1_CpG_16 (P < 0.05) as well as TXNRD1-1_CpG_1 and TXNRD1-1_CpG_2 methylation between OA and KBD patients (P < 0.05). However, there was no difference in methylation levels of other CpGs between the 2 groups (P > 0.05).

CONCLUSION

OA and KBD patients display distinct methylation profiles in the CpG sites of DIO2, GPX3, and TXNRD1 promoter regions. These findings provide a strong background and new perspective for future studies on mechanisms underlying epigenetic regulation of selenoprotein genes associated with OA and KBD diseases.

Identification of HIF-1α/VEGFA signaling pathway and transcription factors in Kashin-Beck disease by integrated bioinformatics analysis

Kashin-Beck disease (KBD) is a chronic and endemic osteoarthropathy. The pathogenesis of KBD has yet to be fully elucidated, although previous studies have shown that its etiology may be associated with low selenium abundance and high exposure to mycotoxins, such as T-2 toxin. In the present study, the Comparative Toxicogenomics Database was used to identify key genes associated with KBD, T-2 toxin and selenium. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were used to identify the biological processes and pathways that key genes may be associated with. By searching the Search Tool for the Retrieval of Interacting Genes database and the Molecular Complex Detection plug-in with Cytoscape, it was possible to construct a KBD-associated protein-protein interaction (PPI) network, and screen the core modules and genes. Western blot analysis was subsequently used to verify the expression levels of hypoxia-inducible factor-1α (HIF-1α) and vascular endothelial growth factor A (VEGFA), two components that are associated with the HIF-1 signaling pathway in KBD disease. Via this approach, a total of 301 key genes were identified that were associated with KBD, T-2 toxin and selenium. The results of the GO and KEGG enrichment analyses demonstrated that these key genes were mainly involved in the process of apoptosis. Previous studies have demonstrated that excessive apoptosis of chondrocytes plays a crucial role in the pathophysiology of KBD, and that HIF-1α has an important role in chondrocyte apoptosis; therefore, the present study was focused on the expression level of HIF-1α in KBD. By analyzing the PPI network constructed from the key genes, a total of 10 core genes were obtained that may be associated with KBD. The results of western blotting experiments revealed that, after treating chondrocytes with different concentrations of T-2 toxin, the expression levels of HIF-1α and VEGFA were markedly downregulated. The iRegulon plug-in for Cytoscape was used to predict the transcription factors that may regulate HIF-1α and VEGFA in the HIF-1 signaling pathway. Using this approach, 10 core genes and 15 transcription factors were obtained. These results may help to clarify the pathogenesis of KBD, thereby providing further avenues for the therapeutic treatment of KBD.

Selenium in Endocrinology-Selenoprotein-Related Diseases, Population Studies, and Epidemiological Evidence

Selenium (Se), apart from iodine, iron, and calcium, is one of the nutrient-derived key elements strongly affecting the endocrine system. However, no specific hormonal "feedback" regulation for Se status has yet been identified, in contrast to the fine-tuned hormone network regulating Ca2+ and phosphate balance or hepcidin-related iron status. Since its discovery as an essential trace element, the effects of Se excess or deficiency on the endocrine system or components of the hypothalamic-pituitary-periphery feedback circuits, the thyroid hormone axis, glucoregulatory and adrenal hormones, male and female gonads, the musculoskeletal apparatus, and skin have been identified. Analysis of the Se status in the blood or via validated biomarkers such as the hepatically derived selenoprotein P provides valuable diagnostic insight and a rational basis for decision making on required therapeutic or preventive supplementation of risk groups or patients. Endocrine-related epidemiological and interventional evidence linking Se status to beneficial or potentially adverse actions of selected selenoproteins mediating most of the (patho-) physiological effects are discussed in this mini-review. Autoimmune thyroid disease, diabetes and obesity, male fertility, as well as osteoporosis are examples for which observational or interventional studies have indicated Se effects. The currently prevailing concept relating Se and selenoproteins to "oxidative stress," reactive oxygen species, radical hypotheses, and related strategies of pharmacological approaches based on various selenium compounds will not be the focus. The crucial biological function of several selenoproteins in cellular redox-regulation and specific enzyme reactions in endocrine pathways will be addressed and put in clinical perspective.

The integrative analysis of DNA methylation and mRNA expression profiles confirmed the role of selenocompound metabolism pathway in Kashin-Beck disease

Kashin-Beck disease (KBD) is an endemic chronic osteochondropathy. The etiology of KBD remains unknown. In this study, we conducted an integrative analysis of genome-wide DNA methylation and mRNA expression profiles between KBD and normal controls to identify novel candidate genes and pathways for KBD. Articular cartilage samples from 17 grade III KBD patients and 17 healthy controls were used in this study. DNA methylation profiling of knee cartilage and mRNA expression profile data were obtained from our previous studies. InCroMAP was performed to integrative analysis of genome-wide DNA methylation profiles and mRNA expression profiles. Gene ontology (GO) enrichment analysis was conducted by online DAVID 6.7. The quantitative real-time polymerase chain reaction (qPCR), Western blot, immunohistochemistry (IHC), and lentiviral vector transfection were used to validate one of the identified pathways. We identified 298 common genes (such as COL4A1, HOXA13, TNFAIP6 and TGFBI), 36 GO terms (including collagen function, skeletal system development, growth factor), and 32 KEGG pathways associated with KBD (including Selenocompound metabolism pathway, PI3K-Akt signaling pathway, and TGF-beta signaling pathway). Our results suggest the dysfunction of many genes and pathways implicated in the pathogenesis of KBD, most importantly, both the integrative analysis and in vitro study in KBD cartilage highlighted the importance of selenocompound metabolism pathway in the pathogenesis of KBD for the first time.

The role of selenium metabolism and selenoproteins in cartilage homeostasis and arthropathies

As an essential nutrient and trace element, selenium is required for living organisms and its beneficial roles in human health have been well recognized. The role of selenium is mainly played through selenoproteins synthesized by the selenium metabolic system. Selenoproteins have a wide range of cellular functions including regulation of selenium transport, thyroid hormones, immunity, and redox homeostasis. Selenium deficiency contributes to various diseases, such as cardiovascular disease, cancer, liver disease, and arthropathy—Kashin–Beck disease (KBD) and osteoarthritis (OA). A skeletal developmental disorder, KBD has been reported in low-selenium areas of China, North Korea, and the Siberian region of Russia, and can be alleviated by selenium supplementation. OA, the most common form of arthritis, is a degenerative disease caused by an imbalance in matrix metabolism and is characterized by cartilage destruction. Oxidative stress serves as a major cause of the initiation of OA pathogenesis. Selenium deficiency and dysregulation of selenoproteins are associated with impairments to redox homeostasis in cartilage. We review the recently explored roles of selenium metabolism and selenoproteins in cartilage with an emphasis on two arthropathies, KBD and OA. Moreover, we discuss the potential of therapeutic strategies targeting the biological functions of selenium and selenoproteins for OA treatment.

Selenium, a micronutrient found in brazil nuts, shiitake mushrooms, and most meats, may aid in treating joint diseases, including the most common form of arthritis, osteoarthritis (OA). In addition to thyroid hormone metabolism and immunity, selenium is important in antioxidant defense. Oxidative damage can destroy cartilage and harm joints, and selenium deficiency is implicated in several joint diseases. Jin-Hong Kim at Seoul National University in South Korea and co-workers reviewed selenium metabolism, focusing on OA and and Kashin–Beck disease, a skeletal development disorder prevalent in selenium-deficient areas of northeast Asia. They report that selenium-containing proteins protect cells against oxidative damage and that selenium is crucial to cartilage production. Further investigation of selenium metabolism may point the way to new treatments for OA and other joint diseases.

Selenium deficiency causes immune damage by activating the DUSP1/NF-κB pathway and endoplasmic reticulum stress in chicken spleen

Selenium (Se) is an essential trace element and its deficiency can lead to immune dysfunction. Many studies have investigated the immune damage caused by Se deficiency in chickens, but its mechanism still needs to be explored. In this study, we fed 1-day-old Hyline male chickens with Se deficient diets (the Se content was 0.008 mg kg^-1 of diet) and a basal diet (the Se content was 0.15 mg kg^-1 of diet). The spleen was collected at the sixth week and used for subsequent experiments. The pathological analysis showed that Se deficiency leads to the destruction of the normal nuclear structure of the spleen cell, and we can observe obvious chromatin condensation and nuclear debris. We constructed a transcriptome database and analyzed the abundance of various genes in the spleen by transcriptome sequence. The analysis of differentially expressed genes (DEGS) showed significant changes in 337 genes, including 210 up-regulations and 127 down-regulations after feeding Se deficient diets. Se deficiency can significantly change oxidative stress and inflammatory response genes in chicken spleen. This study confirmed that Se deficiency increased the IL-2 levels, whereas it down-regulated IL-17, IFN-γ and Foxp3, which indicates that the immune dysfunction of the spleen and Th1/Th2 is imbalanced. We also found that Se deficiency down-regulated some related genes for endoplasmic reticulum Ca²⁺ transport, leading to endoplasmic reticulum stress (ERS). Moreover, we determined that Se deficiency triggered the low expression of DUSP1/NF-κB. In summary, our results indicate that Se deficiency can inhibit the spleen immune function of chickens by regulating the DUSP1/NF-κB pathway and ERS, leading to spleen damage in chickens. Based on transcriptomics research, our results will help further study the harmful effects of Se deficiency.

Effects of methylation of deiodinase 3 gene on gene expression and severity of Kashin-Beck disease

Kashin-Beck disease (KBD) is a complex endemic osteoarthropathy, which mainly occurs in the northeast to southwest China. Iodothyronine deiodinases 3 (DIO3) is one of the selenoproteins, which is closely related to bone metabolism and unclear to KBD. This study aims to investigate the role and associated mechanisms of methylation and expression of DIO3 with disease severity in patients with KBD. We performed a bioinformatics analysis first to identify the biological mechanisms involved in selenoproteins. The methylation status of the DIO3 gene and DIO3 gene expression, as well as DIO3-related regulatory genes in patients with KBD, were analyzed. We found that 15 CpG sites of six selenoproteins were hypomethylated with 5-azacytidine treatment. DIO3 hypermethylation was associated with an increased risk of KBD and may lead to downregulation of DIO3 gene expression as well as be an indicator of the severity of KBD, which may provide a new insight for gene-environment correlations and interactions in etiology and pathogenesis of KBD.

Epitranscriptomic systems regulate the translation of reactive oxygen species detoxifying and disease linked selenoproteins

Here we highlight the role of epitranscriptomic systems in post-transcriptional regulation, with a specific focus on RNA modifying writers required for the incorporation of the 21st amino acid selenocysteine during translation, and the pathologies linked to epitranscriptomic and selenoprotein defects. Epitranscriptomic marks in the form of enzyme-catalyzed modifications to RNA have been shown to be important signals regulating translation, with defects linked to altered development, intellectual impairment, and cancer. Modifications to rRNA, mRNA and tRNA can affect their structure and function, while the levels of these dynamic tRNA-specific epitranscriptomic marks are stress-regulated to control translation. The tRNA for selenocysteine contains five distinct epitranscriptomic marks and the ALKBH8 writer for the wobble uridine (U) has been shown to be vital for the translation of the glutathione peroxidase (GPX) and thioredoxin reductase (TRXR) family of selenoproteins. The reactive oxygen species (ROS) detoxifying selenocysteine containing proteins are a prime examples of how specialized translation can be regulated by specific tRNA modifications working in conjunction with distinct codon usage patterns, RNA binding proteins and specific 3' untranslated region (UTR) signals. We highlight the important role of selenoproteins in detoxifying ROS and provide details on how epitranscriptomic marks and selenoproteins can play key roles in and maintaining mitochondrial function and preventing disease.

Systems Biology of Selenium and Complex Disease

Selenium is an essential trace element for maintenance of overall health, whose deficiency and dyshomeostasis have been linked to a variety of diseases and disorders. The majority of previous researches focused on characterization of genes encoding selenoproteins or proteins involved in selenium metabolism as well as their functions. Many studies in humans also investigated the relationship between selenium and complex diseases, but their results have been inconsistent. In recent years, systems biology and "-omics" approaches have been widely used to study complex and global variations of selenium metabolism and function in physiological and different pathological conditions. The present paper reviews recent progress in large-scale and systematic analyses of the relationship between selenium status or selenoproteins and several complex diseases, mainly including population-based cohort studies and meta-analyses, genetic association studies, and some other omics-based studies. Advances in ionomics and its application in studying the interaction between selenium and other trace elements in human health and diseases are also discussed.

Impact of the gene-gene interactions related to the HIF-1α signaling pathway with the knee osteoarthritis development

INTRODUCTION/OBJECTIVES

Articular cartilage is the target tissue of osteoarthritis (OA), and because it lacks capillary networks, the microenvironment is hypoxic. Hypoxia inducible factor-1 alpha (HIF-1α) regulates the homeostasis of this tissue. The aim of this study was to investigate whether genetic polymorphisms of the HIF-1α signaling pathway are involved in the development of knee OA.

METHOD

We performed a case-control association study and genotyped 134 knee OA patients and 267 healthy controls. All participants were genotyped in order to evaluate 42 SNPs from 22 genes involved in the HIF-1α signaling pathway using the OpenArray technology. Gene-gene interactions (epistasis) were analyzed using the multifactor dimensionality reduction (MDR) method.

RESULTS

The MDR analysis showed epistasis between AKT2 (rs8100018) and IGF1 (rs2288377), AKT2 (rs8100018) and IGF1 (rs35767), IGF1 (rs35767) and COL2A1 (rs1793953), and between GSK3B (rs6438552) and IGF1 (rs35767) polymorphisms, with information gain values of 21.24%, 8.37%, 9.93%, and 5.73%, respectively. Additionally, our model allowed us to identify high- and low-risk genotypes among COL2A1 rs1793953, GSK3B rs6438552, AKT2 rs8100018, and IGF1 rs35767 polymorphisms.

CONCLUSIONS

Knowing the interactions of these polymorphisms involved in HIF-1α signaling pathway could provide a new diagnostic support tool to identify individuals at high risk of developing knee OA.

Genetic variants in the ITPR2 gene are associated with Kashin-Beck Disease in Tibetan

Background

Kashin‐Beck Disease (KBD) is a chronic, endemic osteoarthropathy. Inositol 1,4,5‐triphosphate receptor type 2 (ITPR2) gene is involved in chondrocytes. We speculated that single‐nucleotide polymorphisms (SNPs) in ITPR2 gene may have an association with KBD in Tibetan.

Methods

To prove this hypothesis, a total of eight SNPs (rs1049376, rs11048526, rs11048556, rs11048585, rs16931011, rs10842759, rs2230372, and rs7134213) were selected, and genotyped in 316 KBD patients and 320 controls. The association between ITPR2 variants and KBD risk was assessed by logistic regression analysis.

Results

The study identified significant association (p = 0.019) between KBD and a novel locus, ITPR2 SNP rs11048526 (OR = 1.49, 95% CI = 1.07–2.08). The variant A/G genotype frequency in rs11048526 had a significantly increased risk of KBD in co‐dominant model (OR = 3.70, 95% CI = 1.26–10.89, p = 0.018), dominant model (OR = 3.56, 95% CI = 1.22–10.40, p = 0.020), log‐additive model (OR = 3.00, 95% CI = 1.12–8.00, p = 0.029) after adjusted for age and gender.

Conclusion

The results indicate a potential association between ITPR2 and KBD risk in Tibetan. Further work is required to confirm these results and explore the mechanisms of these effects.

Pathway-based network analyses and candidate genes associated with Kashin-Beck disease

To perform a comprehensive analysis focusing on the biological functions and interactions of Kashin-Beck disease (KBD)-related genes to provide information towards understanding the pathogenesis of KBD.A retrospective, integrated bioinformatics analysis was designed and conducted. First, by reviewing the literature deposited in PubMed, we identified 922 genes genetically associated with KBD. Then, biological function and network analyses were conducted with Cytoscape software. Moreover, KBD specific molecular network analysis was conducted by Cytocluster using the Molecular Complex Detection Algorithm (MCODE).The biological function enrichment analysis suggested that collagen catabolic process, protein activation cascade, cellular response to growth factor stimulus, skeletal system development, and extrinsic apoptosis played important roles in KBD development. The apoptosis pathway, NF-kappa B signaling pathway, and the glutathione metabolism pathway were significantly enriched in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway network, suggesting that these pathways may play key roles in KBD occurrence and development. MCODE clusters showed that in top 3 clusters, 54 of KBD-related genes were included in the network and 110 candidate genes were discovered might be potentially related to KBD.The 110 candidate genes discovered in the current study may be related to the development of KBD. The expression changes of apoptosis and oxidative stress-related genes might serve as biomarkers for early diagnosis and treatment of KBD.

The study of GPX3 methylation in patients with Kashin-Beck Disease and its mechanism in chondrocyte apoptosis

OBJECTIVE

Selenium deficiency is a risk factor for Kashin-Beck Disease (KBD), an endemic osteoarthropathy. Although promoter hypermethylation of glutathione peroxidase 3 (GPX3) (a selenoprotein) has been identified in several cancers, little is known about promoter methylation and expression of GPX3 and their relation to selenium in KBD. The present study was thus conducted to investigate this research question.

METHODS

Methylation and expressions of GPX3 in whole blood drawn from 288 KBD patients and 362 healthy controls and in chondrocyte cell line were evaluated using methylation-specific PCR and qRT-PCR, respectively. The protein levels of PI3K/Akt/c-fos signaling in the whole blood and chondrocyte cell line were determined with Western blotting. Chondrocytes apoptosis were detected by Hoechst 33342 and Annexin V-FITC/PI staining.

RESULTS

GPX3 methylation was increased, GPX3 mRNA was decreased, and protein levels in the PI3K/Akt/c-fos signaling pathway were up-regulated in the whole blood collected from KBD patients as compared with healthy controls. Similar results were obtained for chondrocytes injured by oxidative stress. There was a significant, decreasing trend in GPX3 expression across groups of unmethylation, partial methylation, and complete methylation for GPX3, in sequence. Compared with unmethylation group, protein levels in PI3K/Akt/c-fos pathway were enhanced in partial and complete methylation groups. Treatment of chondrocytes with sodium selenite resulted in reduced methylation and increased expression of GPX3 as well as down-regulated level of PI3K/Akt/c-fos proteins.

CONCLUSIONS

The methylation and expression of GPX3 and expression of PI3K/Akt/c-fos pathway are altered in KBD and these changes are reversible by selenium supplementation.

Selenoprotein S: a therapeutic target for diabetes and macroangiopathy?

Inflammatory response, oxidative stress, and endoplasmic reticulum (ER) stress are important pathophysiological bases of the occurrence and development of diabetes mellitus (DM) and macroangiopathy complications. Selenoprotein S (SELENOS) is involved in the regulation of these mechanisms; therefore, its association with DM and macroangiopathy has gradually received attention from scholars worldwide. SELENOS has different biological functions in different tissues and organs: it exerts antioxidant protection and has anti-ER stress effects in the pancreas and blood vessels, while it promotes the occurrence and development of insulin resistance in the liver, adipose tissue, and skeletal muscle. In addition, studies have confirmed that some SELENOS gene polymorphisms can influence the inflammatory response and are closely associated with the risk for developing DM and macroangiopathy. Therefore, comprehensive understanding of the association between SELENOS and inflammation, oxidative stress, and ER stress may better elucidate and supplement the pathogenic mechanisms of DM and macroangiopathy complications. Furthermore, in-depth investigation of the association of SELENOS function in different tissues and organs with DM and macroangiopathy may facilitate the development of new strategies for the prevention and treatment of DM and macrovascular complications. Here, we summarize the consensus and controversy regarding functions of SELENOS on currently available evidence.

Genetic association, mRNA and protein expression analysis identify ATG4C as a susceptibility gene for Kashin-Beck disease

OBJECTIVE

Recent study observed defective autophagy in chondrocytes with Kashin-Beck Disease (KBD). To clarify the potential role of autophagy-related ATG4C gene in the development of KBD, we conducted an integrative analysis of genetic association, messenger ribonucleic acid (mRNA) and protein expression of ATG4C in KBD patients.

METHODS

1026 subjects (559 KBD patients and 467 healthy cases) were enrolled in discovery association study. Four single nucleotide polymorphisms (SNPs) of ATG4C gene (rs11208030, rs4409690, rs12097658 and rs6587988) were genotyped by Sequenom MassARRAY platform. Association analysis was conducted by PLINK software. The significant SNPs of ATG4C were replicated using an independent sample of 899 subjects (including 90 KBD patients and 809 healthy controls). Ungenotyped SNPs in ATG4C gene were imputed by IMPUTE 2.0. Knee cartilage specimens were collected from five KBD patients and five healthy subjects. Quantitative real-time polymerase chain reaction (qRT-PCR) and western blot were performed to compare the mRNA and protein expression levels of ATG4C between KBD cartilage and control cartilage.

RESULTS

We observed significant association between KBD and rs11208030 (P value = 0.003), rs4409690 (P value = 0.004), rs12097658 (P value = 0.003) and rs6587988 (P value = 0.003) in both discovery and replication samples. The mRNA expression level of ATG4C (ratio = 0.168, P value = 0.007) in KBD chondrocytes was significantly lower than that in normal chondrocytes. Western blot (P value < 0.001) further confirmed the reduced expression of ATG4C protein in both KBD cartilage and chondrocytes.

CONCLUSION

Our results strongly suggest that ATG4C was a novel autophagy-related susceptibility gene of KBD.

Genome-wide association study identifies COL2A1 locus involved in the hand development failure of Kashin-Beck disease

Kashin-Beck disease (KBD) is a chronic osteochondropathy. The pathogenesis of growth and development failure of hand of KBD remains elusive now. In this study, we conducted a two-stage genome-wide association study (GWAS) of palmar length-width ratio (LWR) of KBD, totally including 493 study subjects. Affymetrix Genome Wide Human SNP Array 6.0 was applied for genome-wide SNP genotyping of 90 KBD patients. Association analysis was conducted by PLINK. Imputation analysis was performed by IMPUTE against the reference panel of the 1000 genome project. Two SNPs were selected for replication in an independent validation sample of 403 KBD patients. In the discovery GWAS, significant association was observed between palmar LWR and rs2071358 of COL2A1 gene (P value = 4.68 × 10⁻⁸). In addition, GWAS detected suggestive association signal at rs4760608 of COL2A1 gene (P value = 1.76 × 10⁻⁴). Imputation analysis of COL2A1 further identified 2 SNPs with association evidence for palmar LWR. Replication study observed significant association signals at both rs2071358 (P value = 0.017) and rs4760608 (P value = 0.002) of COL2A1 gene. Based on previous and our study results, we suggest that COL2A1 was a likely susceptibility gene involved in the hand development failure of KBD.

Identified molecular mechanism of interaction between environmental risk factors and differential expression genes in cartilage of Kashin-Beck disease

As environmental risk factors (ERFs) play an important role in the pathogenesis of Kashin-Beck disease (KBD), it is important to identify the interaction between ERFs and differentially expression genes (DEGs) in KBD. The environmental response genes (ERGs) were analyzed in cartilage of KBD in comparison to normal controls.We searched 5 English and 3 Chinese databases from inception to September 2015, to identify case-control studies that examined ERFs for KBD using integrative meta-analysis and systematic review. Total RNA was isolated from articular cartilage of KBD patients and healthy controls. Human whole genome microarray chip (Agilent) was used to analyze the amplified, labeled, and hybridized total RNA, and the validated microarray data were partially verified using real-time quantitative polymerase chain reaction (qRT-PCR). The ERGs were derived from the Comparative Toxicogenomics Database. The identified ERGs were subjected to KEGG pathway enrichment, biological process (BP), and interaction network analyses using the Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.7, and STRING.The trace elements (selenium and iodine), vitamin E, and polluted grains (T-2 toxin/HT-2 toxin, deoxynivalenol, and nivalenol) were identified as the ERFs for KBD using meta-analysis and review. We identified 21 upregulated ERGs and 7 downregulated ERGs in cartilage with KBD compared with healthy controls, which involved in apoptosis, metabolism, and growth and development. KEGG pathway enrichment analysis found that 2 significant pathways were involved with PI3K-Akt signaling pathway and P53 signaling pathway, and gene ontology function analysis found 3 BPs involved with apoptosis, death, and cell death in KBD cartilage.According to previous results and our own research, we suggest that the trace element selenium and vitamin E induce PI3K-Akt signaling pathway and the mycotoxins (T-2 toxin/HT-2 toxin and DON) induce P53 signaling pathway, contributing to the development of KBD, and chondrocyte apoptosis and cell death.

Evaluation of the relationship and genetic overlap between Kashin-Beck disease and body mass index

OBJECTIVES

Body mass index (BMI) is one of the major factors affecting the development of osteoarthritis (OA) but there is currently no information available regarding the relationship between BMI and Kashin-Beck disease (KBD). Our aim in this study was to investigate the relationship and genetic overlap between BMI and KBD.

METHOD

A total of 2050 Han Chinese subjects participated in this study. Using a cohort of 333 grade I KBD patients, logistic regression analysis was conducted to evaluate the correlation between BMI and KBD. Another independent sample of 1717 subjects was genotyped for a genome-wide association study (GWAS) using Affymetrix Human SNP 6.0 Arrays. Single nucleotide polymorphism (SNP) effect concordance analysis (SECA) was applied to the GWAS summaries of KBD and BMI for pleiotropy analysis. Genome-wide bivariate association analysis (GWBAA) of KBD and BMI was carried out to identify the genes with pleiotropic effects on KBD and BMI. The relevance of identified genes with KBD was validated by gene expression profiling and immunohistochemistry.

RESULTS

BMI correlated positively with knee movement disorder in KBD (coefficient β = 0.068, p = 0.045). SECA identified a significant pleiotropic effect (empirical p = 0.021) between KBD and BMI. In the GWBAA, the rs1893577 of the ADAMTS1 gene achieved the most significant association signal (p = 7.38 × 10^-9). ADAMTS1 was also up-regulated in KBD vs. normal (ratio = 2.64 ± 2.80) and KBD vs. OA (ratio = 2.31 ± 2.01). The rate of ADAMTS1-positive chondrocytes in KBD was significantly higher than that in OA (p < 0.05) and healthy controls (p < 0.05).

CONCLUSIONS

Our results suggest that ADAMTS1 is a novel susceptibility gene for KBD.

Membrane-bound selenoproteins

SIGNIFICANCE

Selenoproteins employ selenium to supplement the chemistry available through the common 20 amino acids. These powerful enzymes are affiliated with redox biology, often in connection with the detection, management, and signaling of oxidative stress. Among them, membrane-bound selenoproteins play prominent roles in signaling pathways, Ca(2+) regulation, membrane complexes integrity, and biosynthesis of lipophilic molecules.

RECENT ADVANCES

The number of selenoproteins whose physiological roles, protein partners, expression, evolution, and biosynthesis are characterized is steadily increasing, thus offering a more nuanced view of this specialized family. This review focuses on human membrane selenoproteins, particularly the five least characterized ones: selenoproteins I, K, N, S, and T.

CRITICAL ISSUES

Membrane-bound selenoproteins are the least understood, as it is challenging to provide the membrane-like environment required for their biochemical and biophysical characterization. Hence, their studies rely mostly on biological rather than structural and biochemical assays. Another aspect that has not received much attention is the particular role that their membrane association plays in their physiological function.

FUTURE DIRECTIONS

Findings cited in this review show that it is possible to infer the structure and the membrane-binding mode of these lesser-studied selenoproteins and design experiments to examine the role of the rare amino acid selenocysteine.

Exome sequencing identified FGF12 as a novel candidate gene for Kashin-Beck disease

The objective of this study was to identify novel causal genes involved in the pathogenesis of Kashin-Beck disease (KBD). A representative grade III KBD sib pair with serious skeletal growth and development failure was subjected to exome sequencing using the Illumina Hiseq2000 platform. The detected gene mutations were then filtered against the data of 1000 Genome Project, dbSNP database, and BGI inhouse database, and replicated by a genome-wide association study (GWAS) of KBD. Ninety grade II or III KBD patients with extreme KBD phenotypes and 1627 healthy controls were enrolled in the GWAS. Affymetrix Genome-Wide Human SNP Array 6.0 was applied for genotyping. PLINK software was used for association analysis. We identified a novel 106T>C at the 3'UTR of the FGF12 gene, which has not been reported by now. Sequence alignment observed high conversation at the mutated 3'UTR+106T>C locus across various vertebrates. In the GWAS of KBD, we detected nine SNPs of the FGF12 gene showing association evidence (P value < 0.05) with KBD. The most significant association signal was observed at rs1847340 (P value = 1.90 × 10(-5)). This study suggests that FGF12 was a susceptibility gene of KBD. Our results provide novel clues for revealing the pathogenesis of KBD and the biological function of FGF12.

Selenium, selenoproteins and neurodegenerative diseases

A review of selenium's essential role in normal brain function and its potential involvement in neurodegenerative diseases.