Comparative Analysis of Differentially Expressed Genes in Chondrocytes from Rats Exposed to Low Selenium and T-2 Toxin

The aim of this study was to construct rat models of environmental risk factors for Kashin-Beck disease (KBD) with low selenium and T-2 toxin levels and to screen the differentially expressed genes (DEGs) between the rat models exposed to environmental risk factors. The Se-deficient (SD) group and T-2 toxin exposure (T-2) group were constructed. Knee joint samples were stained with hematoxylin-eosin, and cartilage tissue damage was observed. Illumina high-throughput sequencing technology was used to detect the gene expression profiles of the rat models in each group. Gene Ontology (GO) functional enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathway enrichment analysis were performed and five differential gene expression results were verified by quantitative real-time polymerase chain reaction (qRT-PCR). A total of 124 DEGs were identified from the SD group, including 56 upregulated genes and 68 downregulated genes. A total of 135 DEGs were identified in the T-2 group, including 68 upregulated genes and 67 downregulated genes. The DEGs were significantly enriched in 4 KEGG pathways in the SD group and 9 KEGG pathways in the T-2 group. The expression levels of Dbp, Pc, Selenow, Rpl30, and Mt2A were consistent with the results of transcriptome sequencing by qRT-PCR. The results of this study confirmed that there were some differences in DEGs between the SD group and the T-2 group and provided new evidence for further exploration of the etiology and pathogenesis of KBD.