Next Article in Journal
Euphorbiasteroid Induces Apoptosis as Well as Autophagy through Modulating SHP-1/STAT3 Pathway in Hepatocellular Carcinoma Cells
Next Article in Special Issue
PSD95 as a New Potential Therapeutic Target of Osteoarthritis: A Study of the Identification of Hub Genes through Self-Contrast Model
Previous Article in Journal
FoxH1 Represses the Promoter Activity of cyp19a1a in the Ricefield Eel (Monopterus albus)
Previous Article in Special Issue
In Vitro and In Vivo Biocompatibility Assessment of a Thermosensitive Injectable Chitosan-Based Hydrogel for Musculoskeletal Tissue Engineering
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 18
10.3390/ijms241813703
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Expression Profiles of Long Non-Coding RNAs in the Articular Cartilage of Rats Exposed to T-2 Toxin

by
Fangfang Yu
,
Miao Wang
,
Kangting Luo
,
Lei Sun
,
Shuiyuan Yu
,
Juan Zuo
and
Yanjie Wang
*
School of Public Health, Zhengzhou University, Zhengzhou 450001, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(18), 13703; https://doi.org/10.3390/ijms241813703
Submission received: 15 July 2023 / Revised: 24 August 2023 / Accepted: 28 August 2023 / Published: 5 September 2023
(This article belongs to the Special Issue Advances in Molecular Research of Cartilage)
Browse Figures
Review Reports Versions Notes

Abstract

:
T-2 toxin could induce bone damage. But there is no specific mechanism about the long non-coding RNAs (lncRNAs) involved in T-2 toxin-induced articular cartilage injury. In this study, 24 SD rats were randomly divided into a control group and a T-2 group, which were administered 4% absolute ethanol and 100 ng/g · bw/day of T-2 toxin, respectively. After treatment for 4 weeks, safranin O/fast green staining identified the pathological changes in the articular cartilage of rats, and immunofluorescence verified the autophagy level increase in the T-2 group. Total RNA was isolated, and high-throughput sequencing was performed. A total of 620 differentially expressed lncRNAs (DE-lncRNAs) were identified, and 326 target genes were predicted. Enrichment analyses showed that the target genes of DE-lncRNAs were enriched in the autophagy-related biological processes and pathways. According to the autophagy database, a total of 23 autophagy-related genes were identified, and five hub genes (Foxo3, Foxo1, Stk11, Hdac4, and Rela) were screened using the Maximal Clique Centrality algorithm. The Human Protein Atlas database indicated that Rela and Hdac4 proteins were highly expressed in the bone marrow tissue, while Foxo3, Foxo1, and Stk11 proteins were reduced. According to Enrichr, etoposide and diatrizoic acid were identified as the key drugs. The real-time quantitative PCR results were consistent with the RNA sequencing (RNA-Seq) results. These results suggested that autophagy was involved in the rat articular cartilage lesions induced by T-2 toxin. The lncRNAs of NONRATG014223.2, NONRATG012484.2, NONRATG021591.2, NONRATG024691.2, and NONRATG002808.2, and their target genes of Foxo3, Foxo1, Stk11, Hdac4, and Rela, respectively, were the key regulator factors of autophagy.

1. Introduction

T-2 toxin is a kind of trichothecene mycotoxin produced by various Fusarium species and is known as the most toxic fungal secondary metabolite [1]. It mainly contaminates wheat, barley, maize, oats, and corn, which are closely related to human diet [2]. And it could induce oxidative stress [3], inflammatory response [4], apoptosis, and the autophagy of chondrocytes [5]. This autophagy is critical for the formation and development of cartilage [6]. It is a self-degradative process that is regulated by a series of autophagy-related genes (ATGs), and it could eliminate damaged organelles and maintain intracellular homeostasis [7,8]. Autophagy can be divided into four processes: (1) the activation of autophagy: the ATG1/ULK1-containing complex is involved in this process; (2) autophagosome formation: PI3K-ATG14, ATG18, and the ATG5-ATG12-ATG16L1 complex is involved in autophagosome expansion and maturation as well as in ATG3, ATG4, and ATG7-mediated lipidation of LC3; (3) the transport and fusion of autophagosomes with lysosomes; (4) degradation and recycling: hydrolases degrade the autophagosome, release degradation products into the cytosol, and then recycle them [9]. T-2 toxin is widely recognized as an important etiology of the Kashin–Beck disease (KBD) [10,11]. It is characterized by deep necrosis of the epiphysis plate and articular cartilage with a higher rate of disability, which greatly burdens the diseased person’s family and society [12,13]. Moreover, T-2 toxin also promotes proteoglycans degradation and activates catabolism via the ROS-NF-κB-HIF-2α pathway in the chondrocytes [14,15]. Therefore, it is important to identify the potential mechanism of T-2 toxin-induced articular cartilage damage.
Long non-coding RNAs (lncRNAs) are a class of non-coding RNAs longer than 200 nucleotides, which transcribe under the direction of RNA polymerase II/III [16]. LncRNAs account for a large proportion of the non-coding transcriptome [17]. Since lncRNAs cannot encode any proteins, they are generally considered as transcriptional “noise” by the scientific community [18]. However, increasing evidence suggests that lncRNAs play an important role in chromatin structure and function regulations, genetic transcription modulation, and RNA metabolism, as well as in cell autophagy, apoptosis, and embryonic development [19,20]. Previous studies demonstrated that lncRNAs are involved in many human diseases, including osteoarthritis and KBD [21]. For example, lncRNA HAGLR mediated the chondrocytes’ inflammatory response through the miR-130a-3p/JAK1 axis [22], and lncRNA HOTAIR accelerated mechanical-stimulation-induced apoptosis via regulating the miR-221/BBC3 axis in chondrocytes [23]. Furthermore, lncRNAs could also facilitate extracellular matrix synthesis by interacting with RNA-binding protein and activating the PI3K/AKT pathway [24,25]. However, less research has been focused on lncRNAs’ expression profiles’ change in rat articular cartilage induced by T-2 toxin.
In the current study, safranin O/fast green staining was used to measure the pathological changes in articular cartilage. The lncRNA profiles of SD rat articular cartilage in the T-2 group and the control group were obtained using high-throughput sequencing technology. Then, the target genes of differentially expressed lncRNAs (DE-lncRNAs) were predicted, and gene ontology (GO) and Kyoto Encyclopedia of Genes Genomes database (KEGG) enrichment analyses were performed. Immunofluorescence was performed to indicate the change in the autophagy level of the articular cartilage in the T-2 group. Consequently, a protein–protein interaction (PPI) network was constructed for screening the hub genes, the Human Protein Atlas (HPA) database was used to validate the change of hub proteins, and a real-time quantitative polymerase chain reaction (RT-qPCR) was used to validate the results of RNA sequencing (RNA-Seq).

2. Results

2.1. Safranin O/Fast Green Staining of Articular Cartilage

The safranin O/fast green staining of the articular cartilage in the T-2 group and the control group is shown in Figure 1. In the control group, the cartilage matrix was deep red, and the chondrocytes were closely arranged in a column. After being administered T-2 toxin for 4 weeks, the cartilage matrix was lightly stained, and the range of the empty lacunae increased, which suggested proteoglycan loss and chondrocytes death, respectively. Moreover, an irregular arrangement of chondrocytes can be observed in the T-2 group compared with the control group.

2.2. Quality Control of RNA Sequencing Data and Read Mapping

A total of 58.10 million reads and 22,127 lncRNAs were obtained. After quality control, we obtained high-quality reads with Q30 > 95%, clean reads% > 99.0%, and the total bases of R1 plus R2 ranged from 12.6 to 14.1 Gb (Table 1). The data achieved the quality control standard with Q30 > 80%, clean reads% > 90.0%, and the total bases of R1 plus R2 were >12 Gb in 90% of the samples, which suggested the good quality of the RNA-Seq. After these reads were mapped to the reference database, the rate was from 89.7% to 90.8%. Taken together, these data indicated that the sequencing quality was good and sufficient for subsequent analysis.

2.3. Identification of Differentially Expressed lncRNAs

According to p-value < 0.05 and | log2 (fold change) | > 1, a total of 620 DE-lncRNAs were identified in the T-2 group compared with the control group, of which 360 were upregulated and 260 were downregulated. The volcano plot shows the distribution of DE-lncRNAs (Figure 2A). In addition, Figure 2B shows a hierarchical cluster analysis of DE-lncRNAs in the form of a heatmap, and the result showed that three samples (i.e., T-2_1, T-2_2, and T-2_3) from the T-2 group were clustered in the same cluster, while three samples (i.e., control_1, control_2, and control_3) from the control group were clustered in the other cluster. Significant differences were observed between the T-2 group and the control group, thereby validating the reliability of our sequencing data.

2.4. Functional Enrichment Analyses of Target Genes in the Cartilage Tissues Samples

According to minimum free energy, 326 target genes of DE-lncRNAs were predicted. Then, GO and KEGG analyses were performed on these genes. According to the gene ratio, the top 10 biological process (BP), cellular component (CC), and molecular function (MF) terms are presented in Figure 3A. BP terms were found to be significantly enriched in the regulation of actin filament polymerization, regulation of actin filament length, as well as bone development. CC terms were found to be enriched in the actin filament, autophagosome, cortical actin cytoskeleton, and mitochondrial outer membrane. Moreover, MF terms were related to protein kinase regulator activity, kinase regulator activity, and microtubule binding. Next, we performed KEGG analysis, and the top 20 KEGG pathways of the target genes are depicted in Figure 3B. The pathways were found to be significantly enriched in autophagy-animal, endocytosis, mitophagy-animal, the AMPK signaling pathway, and the FoxO signaling pathway.

2.5. T-2 Toxin Triggered Autophagy in Articular Cartilage

As shown in Figure 4, immunofluorescence assay results demonstrated that the LC3 protein level of articular cartilage increased in the T-2 group compared with the control group. These results suggested that the autophagy level was highly expressed in the articular cartilage administered with T-2 toxin.

2.6. PPI Network Analysis and the Validation of Hub Genes

A total of 23 ATGs were mapped according to the autophagy database. Then, they were imported in STRING, and the results were presented using Cytoscape. As shown in Figure 5A, the PPI network of the target genes consisted of 14 nodes and 16 edges. Additionally, the top five hub genes were ranked according to the Maximal Clique Centrality (MCC) algorithm (Figure 5B): Foxo3, Foxo1, Hdac4, Rela, and Stk11. The top five hub proteins were then screened using the HPA database and the results are shown in Figure 5C. The Rela and Hdac4 proteins were highly expressed in the bone marrow tissue, whereas the levels of Foxo3, Foxo1, and Stk11 reduced in the bone marrow tissue.

2.7. Prediction of Candidate Drugs

The potential intervention drugs of T-2 toxin-induced autophagy were predicted via the Enrichr platform. The top 10 candidate drugs are listed in Table 2. The results showed that etoposide (CTD 00005948) and diatrizoic acid (CTD 00005787) were the two key drugs.

2.8. Validation of DE-lncRNAs via RT-qPCR

To further validate the results of RNA-Seq, the lncRNAs of NORATG013414.2 and NONRATG026881.1 were randomly selected for RT-qPCR analysis. These results showed that the relative expression levels of NONRATG026881.1 and NONRATG013414.2 were all upregulated in the T-2 group compared with the control group (Figure 6A,B). According to log2 (fold change), the expression trends were consistent with the RNA-Seq results (Figure 6C).

3. Discussion

It is well known that T-2 toxin is harmful to human health and causes cytotoxicity, immunotoxicity, and skeletal toxicity, among other harmful effects [26]. However, one of the most severe consequences of T-2 toxin exposure is articular cartilage damage [27]. However, the function of lncRNAs in the T-2 toxin-induced cartilage damage remains unclear. Therefore, this study investigated the change in lncRNAs’ expression profiles in rat articular cartilage using RNA-Seq. Subsequently, we performed functional enrichment analyses of DE-lncRNAs’ target genes. These results suggested that autophagy is a key mechanism for T-2 toxin-induced cartilage damage.
In this study, safranin O/fast green staining revealed the pathological changes in the rat articular cartilage samples that were exposed to T-2 toxin for 4 weeks. In total, we identified 620 DE-lncRNAs and predicted 326 target genes according to the articular cartilage profile. The 326 identified target genes were then subjected to enrichment analyses. According to the GO analysis, the target genes were mainly enriched in the negative regulation of cytoskeleton organization, regulation of actin filament polymerization, regulation of actin polymerization or depolymerization, and bone development process. As we all know, actin is the most abundant protein that participates in protein–protein interactions in numerous eukaryotic cells, and studies have reported that the balance between polymerization and depolymerization of actin is crucial for therapeutically controlling the chondrocyte phenotype [28]. Next, the KEGG pathway enrichment analysis result showed that autophagy-animal, endocytosis, mitophagy-animal, the AMPK signaling pathway, and the FoxO signaling pathway were significantly enriched. Autophagy is a process involving lysosome degradation under various stress conditions, such as DNA damage [29], toxin exposure [30], nutrient deficiency [31], and hypoxia [32]. Autophagy has been further divided into microautophagy, chaperone-mediated autophagy, and macroautophagy. LC3-associated endocytosis is regarded as a noncanonical autophagy pathway [33]. Moreover, previous studies have found that T-2 toxin exposure can induce an autophagic response that is protective in nature [34]. Thus, moderate autophagy could relieve T-2 toxin-induced toxicity by removing damaged mitochondria [35]. However, excessive activation of autophagy may lead to the degradation of basic cellular components, thereby leading to cell death, further inducing the formation of articular cartilage lesions [36]. In the process of mitophagy, PINK1 recruits and activates Parkin on the surface of dysfunctional mitochondria in response to oxidative stress. Then, Parkin ubiquitinates the proteins on the outer surface of the mitochondrial membrane. Fusion of autophagosome and lysosome can degrade the damaged mitochondrion, further relieving the damages induced by T-2 toxin exposure [37]. AMPK is a serine/threonine protein kinase that acts as an energy sensor in cells. Under pathological conditions, AMPK is phosphorylated and activated, thereby promoting autophagy via ATGs (i.e., ULK1, PIK3C3/VPS34, and mTORC1). Moreover, it can inhibit mitophagy via the regulation of transcription factors located downstream of Foxo3 [38]. Similarly, Foxo induces autophagy via transactivation of autophagy genes and regulation of autophagy activity. Post-translational modification and epigenetic modulation of Foxo can induce or inhibit autophagy [39,40]. Our KEGG enrichment result substantiated autophagy is involved in articular cartilage lesions exposed to T-2 toxin. Furthermore, immunofluorescence assay verified that the level of autophagy in articular cartilage was higher in the T-2 group relative to the control group.
PPI network identified five hub ATGs—Foxo3, Foxo1, Hdac4, Stk11, and Rela—which play important roles in autophagy. These ATGs have a variety of functions; for instance, they can deliver extracellular cargo to the lysosome, promote localization to the plasma membrane and trigger extracellular release of intracellular cargo, and coordinate with various other cell signaling pathways [8]. Foxo3 and Foxo1, as members of the Foxo family, induce autophagy by directly interacting with autophagy proteins [40]. In our study, Foxo3 and Foxo1 were regulated by the lncRNAs of NONRATG014223.2 and NONRATG012484.2, respectively. Foxo3 can activate the transcription of ULK1 and PIK3CA, which further activates the initiation of autophagy [41]. Under stress, Foxo1 is acetylated and binds to ATG7, thereby driving classical degradative autophagy via ATG8 lipidation or LC3-associated phagocytosis [42,43]. Studies have also demonstrated that target regulating Foxo1 in chondrocytes can prevent osteoarthritis via autophagy [44]. Hdac4 is a type of histone deacetylase located in the nucleus and cytoplasm [45] that we found to be regulated by lncRNA NONRATG024691.2. Moreover, Hdac4 can activate Foxo3 that subsequently interacts with ATG5 and LC3II, further inducing autophagy [46]. Stk11, also named LKB1, is a serine/threonine kinase that is widely expressed in tissues. In our study, Stk11 was found to be regulated by NONRATG021591.2. Stk11 activates the phosphorylation of the threonine 172 site on the AMPK α subunit. Next, AMPK-induced autophagy proceeds by inhibiting the mTOR signaling pathway and triggering ATG5 activity [47,48]. Rela (p65), an NF-κB subunit, always activates the NF-κB signaling pathway [49]. NF-κΒ triggers autophagy by directly inducing the expression of Beclin1, ATG5, and LC3 [50]. We found that Rela was regulated by NONRATG002808.2, and its specific regulatory mechanism is shown in Figure 7. Furthermore, HPA database indicated that the proteins of Hdac4 and Rela were highly expressed in the bone marrow tissue, but proteins of Foxo3, Foxo1, and Stk11 were reduced.
Next, we used the hub ATGs to predict potential drugs. The results suggested that etoposide and diatrizoic acid are the important pharmacological targets. Etoposide is a non-alkaloid lignan that is isolated from the dried roots and rhizomes of Podophyllum peltatum or Podophyllum emodi. It is used for the treatment of numerous cancer types [51,52]. Additionally, etoposide is also a derivative of podophyllotoxin that exhibits strong anti-inflammation and anti-oxidative activities [53]. Therefore, etoposide could be used to help scavenge the reactive oxygen and intervene in the autophagy response, thereby relieving T-2 toxin-induced cartilage tissue damage. Diatrizoic acid is an iodinated aromatic that is often used as a radiocontrast agent, and it can also induce mitophagy and oxidative stress via calcium dysregulation [54]. Studies have demonstrated that diatrizoic acid is a component of FocusClear™, which has been used for optical clearing of mouse brain, skin, bone, and lymph node cortex tissues [55,56]. Moreover, oral meglumine diatrizoate esophagography could screen esophageal fistulas during radiotherapy or chemoradiotherapy with minimal side effects [57]. Therefore, we speculated that diatrizoic acid could also be used for treating articular cartilage damage.
However, our study has certain limitations. First, an animal model was used to analyze DE-lncRNAs in the articular cartilage of rats exposed to T-2 toxin, but this was not tested in humans in a clinical setting. Second, the specific regulatory mechanisms between the lncRNAs and their target genes were not established during the experimental exploration. Further molecular biology studies are required to verify the regulatory mechanisms involved.

4. Materials and Methods

4.1. Establishment of Animal Models Treated with T-2 Toxin

Twenty-four specific-pathogen-free male Sprague–Dawley (SD) rats, weighing approximately 60–80 g at the age of 4 weeks, were purchased from the Henan provincial laboratory animal center. Rats were adaptively fed in a room with a humidity of 55% ± 5% and a temperature of 22 °C ± 2 °C under a 12 h light/dark cycle for 1 week [27,58]. Then, the rats were weighed and randomly divided into two groups: the control group rats were administered 4% absolute ethanol via gavage and the T-2 rats were administered 100 ng/g · bw/day T-2 toxin via gavage [27,59]. After 4 weeks of treatment, the articular cartilages of the twelve rats in each group were collected. The animal experiments were reviewed and approved by the Animal Ethics Research Committee of Zhengzhou University (ZZUIRB2021-69).

4.2. Safranin O/Fast Green Staining

Three articular cartilage samples from each group were fixed in 4% paraformaldehyde for 24 h and decalcified in EDTA decalcifying solution for 4 weeks. Next, each sample was embedded, sliced, and dewaxed to generate cartilage tissue sections. These sections were stained with Weigert solution and rinsed with water before placing in 0.2% fast green solution for 5 min, 1% ethylic acid solution for 15 s, and 0.1% safranin O solution for 5 min. Then, the samples were dehydrated, transparentized, and mounted with neutral balsam. The pathological changes in the articular cartilage between the T-2 group and the control group were scanned using a scanner (3DHISTECH, Budapest, Öv u. 3., Hungary).

4.3. Library Construction and RNA-Seq

For RNA-Seq, we first isolated the total RNA from the articular cartilage of three rats from each group using TRIzol (Thermo Fisher Scientific, Waltham, MA, USA). Qubit 2.0 fluorometer (Invitrogen, Carlsbad, CA, USA) was used to quantify the concentration of total RNA (sample concentration ≥ 100 ng/µL, total ≥ 2 μg), NanoDrop 2000 (Thermo Fisher Scientific, Waltham, MA, USA) was used to measure the RNA purity with OD260/280 in the range of 1.8–2.2 and OD260/230 ≥ 2.0. The RNA degradation was analyzed using agarose gel electrophoresis, and Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) was used to evaluate the RNA integrity (RNA integrity number ≥ 7).
To ensure the construction of high-quality RNA-Seq libraries, moderate amounts of RNA were sampled for ribosomal RNA (rRNA) depletion. Subsequently, RNAs were purified and digested into 100–300 bp fragments. The first complementary DNA (cDNA) chain was synthesized with short fragment primers and followed by the second cDNA chain synthesis. The 3′-ends of the resulting double-stranded cDNA were adenylated and ligated to Illumina sequencing adaptors using T4 DNA ligase. Subsequently, the initial library was purified using an Agencourt SPRIselect reagent kit (Beckman Coulter, Miami, FL, USA). The library with a peak fragment size of 300 bp was selected and then amplified using PCR. After amplification, a Qubit 2.0 (Invitrogen, Carlsbad, CA, USA) was used to quantify the DNA concentration of the constructed libraries. In addition, an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) was used to determine the size distribution of the library fragments. After quality inspection, an Illumina HiSeq2500 (Illumina, San Diego, CA, USA) was used to sequence the constructed libraries via a 2 × 150 bp paired sequencing strategy.

4.4. Identification of Differentially Expressed lncRNAs

FastQC software v0.12.1 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc, accessed on 13 August 2022) was used to filter out the adaptor and low-quality reads. The filter criteria as follows: (1) sequencing primers, (2) low quality sequences near the terminal (Q < 15), and (3) sequences with fragment length < 40 bp. Then, the high-quality clean reads were obtained and were mapped to the reference database using STAR software v2.7.11a (https://github.com/alexdobin/STAR, accessed on 13 August 2022. The results were subsequently analyzed using Picard v3.1.0 (http://broadinstitute.github.io/picard/, accessed on 13 August 2022). The expression levels of these lncRNAs were calculated based on fragments per kilobase of genes per million fragments mapped (FPKM) values. DE-lncRNAs were defined as p-value < 0.05 and | log2 (fold change) | > 1 via the Deseq2 software v1.40.2 (https://bioconductor.org/packages/release/bioc/html/DESeq2.html, accessed on 13 August 2022) [60]. Hierarchical cluster analysis was also used to determine the expression patterns of DE-lncRNAs.

4.5. Target Gene Prediction

The target genes of identified DE-lncRNAs were predicted using lncTar software v1.0 (http://www.cuilab.cn/lnctar, accessed on 13 August 2022). In this project, the minimum free energy of lnRNAs and mRNA binding sites were calculated to determine whether they could produce stable binding sites, and the default threshold of minimum free energy was −0.1.

4.6. GO and KEGG Pathway Enrichment Analyses

The target genes of identified DE-lncRNAs were subjected to GO and KEGG enrichment analyses using R software v4.3.1 (https://mirrors.tuna.tsinghua.edu.cn/CRAN/, accessed on 12 May 2023) cluster Profiler that annotates the functions of genes. In GO analysis, the vocabularies included the BP, CC, and MF. In KEGG analysis, the enrichment scores were used to rank the pathways and the pathways with p-value < 0.05 were considered significant.

4.7. Immunofluorescence Assay

The paraffin tissues of articular cartilage were sectioned, dewaxed in xylene, and dehydrated using gradient ethanol. After antigen retrieval, the sections were incubated with 3% BSA for 30 min. Then, they were incubated with primary antibodies against microtubule-associated protein 1 light chain 3 alpha (LC3, 1:400, 14600-1-AP, Proteintech, Wuhan, China) at 4 °C overnight, followed by rinsing the sections thrice in PBS solution and incubation with Cy3-conjugated secondary antibodies. Next, the DAPI staining solution was utilized to stain the nucleus of the cell. Immunofluorescence images were captured, and the fluorescence images were quantified using Image-Pro Plus 6.0 (Media Cybernetics, Rockville, MD, USA).

4.8. PPI Network Construction and Validation of Hub Genes

The autophagy database (http://www.tanpaku.org/autophagy/, accessed on 12 May 2023) is a manually compiled and managed database that provides up-to-date information regarding autophagy-related proteins and their homologs in 41 eukaryote species [61]. In the current study, it was used to map the autophagy-related target genes. The STRING database (https://cn.string-db.org/, accessed on 12 May 2023) was used to analyze the information regarding interactions among these mapped genes. These interactions were reflected by the comprehensive score. With a confidence score (0.4) set as the minimum interaction score, the autophagy-related target genes were used to construct a PPI network [62]. These results were presented and analyzed using Cytoscape v3.9.0 (https://cytoscape.org/, accessed on 12 May 2023) [63]. Next, we used cyto-Hubba, a Cytoscape plugin, to screen the top five hub genes via MCC [64].
The HPA database (https://www.proteinatlas.org/, accessed on 13 May 2023) consists of cell, tissues, and pathology atlas that provides the information regarding transcriptome and protein mapping for specific human tissues. In the current study, HPA database was used to obtain the data on hub protein expression levels in the bone marrow and skeletal muscle tissues.

4.9. Prediction of Candidate Pharmacological Targets

The identified five hub genes were applied to predict the candidate pharmacological targets using Enrichr v3.2.0 (https://maayanlab.cloud/Enrichr/, accessed on 13 May 2023). Enrichr is a comprehensive search engine that integrates biological knowledge with a comprehensive database of curated gene sets from previous studies [65]. Currently, Enrichr contains 211 gene set libraries, 431,749 terms, and 60,534,815 gene sets for analysis. Notably, it is important to evaluate protein–drug interactions to understand the structural features required for receptor sensitivity. The candidate pharmacological targets of the five hub genes were sorted by p-value from small to large, and the p-value of < 0.05 was considered statistically significant [66].

4.10. Verification Analysis

Total RNA from cartilage tissue was reverse-transcribed into cDNA using SuperScript IV reverse transcriptase (Thermo Fisher Scientific, Waltham, MA, USA). Two lncRNAs were then randomly selected for RT-qPCR using SYBR green (Gene Copoeia Inc., Rockville, MD, USA). GAPDH was used as an internal reference gene. The experimental conditions of the PCR amplification were set as follows: 94 °C for 5 min, followed by 40 cycles at 94 °C for 10 s, 55 °C for 20 s, and 72 °C for 20 s. The lncRNAs’ expression levels were calculated with the 2−ΔΔCt method. Three replicates were used for each group.

5. Conclusions

Autophagy was involved in the process of T-2 toxin-induced articular cartilage lesions. The lncRNAs of NONRATG014223.2, NONRATG012484.2, NONRATG021591.2, NONRATG024691.2, and NONRATG002808.2, and their target genes of Foxo3, Foxo1, Stk11, Hdac4, and Rela, respectively, were the key regulatory factors of autophagy.

Author Contributions

Investigation—F.Y., M.W., K.L., L.S., S.Y. and J.Z.; Methodology—F.Y., M.W., K.L., L.S., S.Y., J.Z. and Y.W.; Resources—F.Y. and Y.W.; Formal analysis, F.Y., M.W., K.L., L.S., S.Y., J.Z. and Y.W.; Supervision—F.Y., L.S., S.Y. and J.Z.; Validation—M.W. and K.L., Writing—original draft, F.Y. and M.W.; Writing—review and editing, M.W., F.Y. and Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Scientific Foundation of China (NSFC No. 82003400), Science and Technology Program of Henan Province (232102311079, 212102310622).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Animal Ethics Research Committee of Zhengzhou University (ZZUIRB2021-69, date: 9 July 2021).

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Jincai Luo from Sanmenxia Center for Disease Control and Prevention helped us review manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Janik, E.; Niemcewicz, M.; Podogrocki, M.; Ceremuga, M.; Stela, M.; Bijak, M. T-2 toxin-the most toxic trichothecene mycotoxin: Metabolism, toxicity, and decontamination strategies. Molecules 2021, 26, 6868. [Google Scholar] [CrossRef] [PubMed]
  2. Haidukowski, M.; Somma, S.; Ghionna, V.; Cimmarusti, M.T.; Masiello, M.; Logrieco, A.F.; Moretti, A. Deoxynivalenol and T-2 Toxin as Major Concerns in Durum Wheat from Italy. Toxins 2022, 14, 627. [Google Scholar] [CrossRef] [PubMed]
  3. Yin, H.; Han, S.; Chen, Y.; Wang, Y.; Li, D.; Zhu, Q. T-2 Toxin Induces Oxidative Stress, Apoptosis and Cytoprotective Autophagy in Chicken Hepatocytes. Toxins 2020, 12, 90. [Google Scholar] [CrossRef]
  4. Dai, C.; Das Gupta, S.; Wang, Z.; Jiang, H.; Velkov, T.; Shen, J. T-2 toxin and its cardiotoxicity: New insights on the molecular mechanism s and therapeutic implications. Food Chem. Toxicol. 2022, 167, 113262. [Google Scholar] [CrossRef] [PubMed]
  5. Yu, F.-F.; Lin, X.-L.; Wang, X.; Ping, Z.-G.; Guo, X. Comparison of Apoptosis and Autophagy in Human Chondrocytes Induced by the T-2 and HT-2 Toxins. Toxins 2019, 11, 260. [Google Scholar] [CrossRef]
  6. Kao, W.-C.; Chen, J.-C.; Liu, P.-C.; Lu, C.-C.; Lin, S.-Y.; Chuang, S.-C.; Wu, S.-C.; Chang, L.-H.; Lee, M.-J.; Yang, C.-D.; et al. The Role of Autophagy in Osteoarthritic Cartilage. Biomolecules 2022, 12, 1357. [Google Scholar] [CrossRef]
  7. Glick, D.; Barth, S.; MacLeod, K.F. Autophagy: Cellular and molecular mechanisms. J. Pathol. 2010, 221, 3–12. [Google Scholar] [CrossRef] [PubMed]
  8. Levine, B.; Kroemer, G. Biological Functions of Autophagy Genes: A Disease Perspective. Cell 2019, 176, 11–42. [Google Scholar] [CrossRef]
  9. Li, W.; He, P.; Huang, Y.; Li, Y.-F.; Lu, J.; Li, M.; Kurihara, H.; Luo, Z.; Meng, T.; Onishi, M.; et al. Selective autophagy of intracellular organelles: Recent research advances. Theranostics 2021, 11, 222–256. [Google Scholar] [CrossRef]
  10. Yu, F.-F.; Zuo, J.; Sun, L.; Yu, S.-Y.; Lei, X.-L.; Zhu, J.-H.; Zhou, G.-Y.; Guo, X.; Ba, Y. Animal models of Kashin-Beck disease exposed to environmental risk factors: Methods and comparisons. Ecotoxicol. Environ. Saf. 2022, 234, 113419. [Google Scholar] [CrossRef]
  11. Jia, Y.; Shi, S.; Cheng, B.; Cheng, S.; Liu, L.; Meng, P.; Yang, X.; Chu, X.; Wen, Y.; Zhang, F.; et al. Fluorine impairs carboxylesterase 1-mediated hydrolysis of T-2 toxin and increases its chondrocyte toxicity. Front. Nutr. 2022, 9, 935112. [Google Scholar] [CrossRef] [PubMed]
  12. Guo, X.; Ma, W.-J.; Zhang, F.; Ren, F.-L.; Qu, C.-J.; Lammi, M. Recent advances in the research of an endemic osteochondropathy in China: Kashin-Beck disease. Osteoarthr. Cartil. 2014, 22, 1774–1783. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, K.; Yu, J.; Liu, H.; Liu, Y.; Liu, N.; Cao, Y.; Zhang, X.; Sun, D. Endemic Kashin-Beck disease: A food-sourced osteoarthropathy. Semin. Arthritis Rheum. 2020, 50, 366–372. [Google Scholar] [CrossRef]
  14. Li, S.-Y.; Cao, J.-L.; Shi, Z.-L.; Chen, J.-H.; Zhang, Z.-T.; Hughes, C.E.; Caterson, B. Promotion of the articular cartilage proteoglycan degradation by T-2 toxin and selenium protective effect. J. Zhejiang Univ. B 2008, 9, 22–33. [Google Scholar] [CrossRef]
  15. Tian, J.; Yan, J.; Wang, W.; Zhong, N.; Tian, L.; Sun, J.; Min, Z.; Ma, J.; Lu, S. T-2 toxin enhances catabolic activity of hypertrophic chondrocytes through ROS-NF-κB-HIF-2α pathway. Toxicol. Vitr. 2012, 26, 1106–1113. [Google Scholar] [CrossRef]
  16. Bridges, M.C.; Daulagala, A.C.; Kourtidis, A. Lnccation: LncRNA localization and function. J. Cell Biol. 2021, 220, e202009045. [Google Scholar] [CrossRef] [PubMed]
  17. Yang, Z.; Jiang, S.; Shang, J.; Jiang, Y.; Dai, Y.; Xu, B.; Yu, Y.; Liang, Z.; Yang, Y. LncRNA: Shedding light on mechanisms and opportunities in fibrosis and aging. Ageing Res. Rev. 2019, 52, 17–31. [Google Scholar] [CrossRef] [PubMed]
  18. Murillo-Maldonado, J.M.; Riesgo-Escovar, J.R. The various and shared roles of lncRNAs during development. Dev. Dyn. 2019, 248, 1059–1069. [Google Scholar] [CrossRef]
  19. Statello, L.; Guo, C.J.; Chen, L.L.; Huarte, M. Gene regulation by long non-coding RNAs and its biological functions. Nat. Rev. Mol. Cell Biol. 2021, 22, 96–118. [Google Scholar] [CrossRef]
  20. Constanty, F.; Shkumatava, A. lncRNAs in development and differentiation: From sequence motifs to functional characterization. Development 2021, 148, 182741. [Google Scholar] [CrossRef]
  21. Dai, Y.; Jian, C.; Wang, X.; Dai, X. Comprehensive expression profiles of mRNAs, lncRNAs and miRNAs in Kashin-Beck disease identified by RNA-sequencing. Mol. Omics 2022, 18, 154–166. [Google Scholar] [CrossRef]
  22. Zuo, Y.; Xiong, C.; Gan, X.; Xie, W.; Yan, X.; Chen, Y.; Li, X. LncRNA HAGLR silencing inhibits IL-1β-induced chondrocytes inflammatory injury via miR-130a-3p/JAK1 axis. J. Orthop. Surg. Res. 2023, 18, 203. [Google Scholar] [CrossRef] [PubMed]
  23. Zheng, T.; Huang, J.; Lai, J.; Zhou, Q.; Liu, T.; Xu, Q.; Ji, G.; Ye, Y. Long non-coding RNA hotairincreased mechanical stimulation-induced apoptosis by regulating microRNA-221/BBC3 axis in C28/I2 cells. Bioengineered 2021, 12, 10734–10744. [Google Scholar] [CrossRef] [PubMed]
  24. Huang, H.; Yan, J.; Lan, X.; Guo, Y.; Sun, M.; Zhao, Y.; Zhang, F.; Sun, J.; Lu, S. LncRNA WDR11-AS1 Promotes Extracellular Matrix Synthesis in Osteoarthritis by Directly Interacting with RNA-Binding Protein PABPC1 to Stabilize SOX9 Expression. Int. J. Mol. Sci. 2023, 24, 817. [Google Scholar] [CrossRef]
  25. Xu, S.; Wang, J.; Jiang, J.; Song, J.; Zhu, W.; Zhang, F.; Shao, M.; Xu, H.; Ma, X.; Lyu, F. TLR4 promotes microglial pyroptosis via lncRNA-F630028o10Rik by activating PI3K/AKT pathway after spinal cord injury. Cell Death Dis. 2020, 11, 693. [Google Scholar] [CrossRef]
  26. Li, S.-J.; Zhang, G.; Xue, B.; Ding, Q.; Han, L.; Huang, J.-C.; Wu, F.; Li, C.; Yang, C. Toxicity and detoxification of T-2 toxin in poultry. Food Chem. Toxicol. 2022, 169, 113392. [Google Scholar] [CrossRef] [PubMed]
  27. Shi, L.; Liu, Q.; Yang, H.; Wang, Q.; Wang, J.; Fan, Y. Inflammation-related pathways involved in damaged articular cartilage of rats exposed to T-2 toxin based on RNA-sequencing analysis. Front. Genet. 2022, 13, 1079739. [Google Scholar] [CrossRef] [PubMed]
  28. Lauer, J.C.; Selig, M.; Hart, M.L.; Kurz, B.; Rolauffs, B. Articular Chondrocyte Phenotype Regulation through the Cytoskeleton and the Signaling Processes That Originate from or Converge on the Cytoskeleton: Towards a Novel Understanding of the Intersection between Actin Dynamics and Chondrogenic Function. Int. J. Mol. Sci. 2021, 22, 3279. [Google Scholar] [CrossRef]
  29. Juretschke, T.; Beli, P. Causes and consequences of DNA damage-induced autophagy. Matrix Biol. 2021, 100, 39–53. [Google Scholar] [CrossRef]
  30. Yang, Y.; Huang, X.; Yuan, W.; Xiang, Y.; Guo, X.; Wei, W.; Soberón, M.; Bravo, A.; Liu, K. Bacillus thuringiensis cry toxin triggers autophagy activity that may enhance cell death. Pestic. Biochem. Physiol. 2021, 171, 104728. [Google Scholar] [CrossRef]
  31. Kajikawa, M.; Fukuzawa, H. Algal Autophagy Is Necessary for the Regulation of Carbon Metabolism Under Nutrient Deficiency. Front. Plant Sci. 2020, 11, 36. [Google Scholar] [CrossRef] [PubMed]
  32. Filomeni, G.; De Zio, D.; Cecconi, F. Oxidative stress and autophagy: The clash between damage and metabolic needs. Cell Death Differ. 2015, 22, 377–388. [Google Scholar] [CrossRef] [PubMed]
  33. Fan, S.; Yue, L.; Wan, W.; Zhang, Y.; Zhang, B.; Otomo, C.; Li, Q.; Lin, T.; Hu, J.; Xu, P.; et al. Inhibition of autophagy by a small molecule through covalent modification of the LC3 protein. Angew. Chem. Int. Ed. Engl. 2021, 60, 26105–26114. [Google Scholar] [CrossRef] [PubMed]
  34. Sun, T.; Zhang, Q.; Li, M.; Tang, S.; Dai, C. T-2 Toxin Induces Apoptotic Cell Death and Protective Autophagy in Mouse Microglia BV2 Cells. J. Fungi 2022, 8, 761. [Google Scholar] [CrossRef] [PubMed]
  35. Li, W.; Jiang, W.S.; Su, Y.R.; Tu, K.W.; Zou, L.; Liao, C.R.; Wu, Q.; Wang, Z.H.; Zhong, Z.M.; Chen, J.T.; et al. PINK1/PARKIN-mediated mitophagy inhibits osteoblast apoptosis induced by advanced oxidation protein products. Cell Death Dis. 2023, 14, 88. [Google Scholar] [CrossRef]
  36. Dai, C.; Xiao, X.; Sun, F.; Zhang, Y.; Hoyer, D.; Shen, J.; Tang, S.; Velkov, T. T-2 toxin neurotoxicity: Role of oxidative stress and mitochondrial dysfunction. Arch. Toxicol. 2019, 93, 3041–3056. [Google Scholar] [CrossRef] [PubMed]
  37. Osman, S. PINK spots: Diseased mitochondria prepare for mitophagy. Nat. Struct. Mol. Biol. 2022, 29, 82. [Google Scholar] [CrossRef]
  38. Herzig, S.; Shaw, R.J. AMPK: Guardian of metabolism and mitochondrial homeostasis. Nat. Rev. Mol. Cell Biol. 2018, 19, 121–135. [Google Scholar] [CrossRef]
  39. Lee, S.; Dong, H.H. FoxO integration of insulin signaling with glucose and lipid metabolism. J. Endocrinol. 2017, 233, R67–R79. [Google Scholar] [CrossRef]
  40. Cheng, Z. The FoxO-autophagy axis in health and disease. Trends Endocrinol. Metab. 2019, 30, 658–671. [Google Scholar] [CrossRef]
  41. Guo, X.; Li, Z.; Zhu, X.; Zhan, M.; Wu, C.; Ding, X.; Peng, K.; Li, W.; Ma, X.; Lv, Z.; et al. A coherent FOXO3-SNAI2 feed-forward loop in autophagy. Proc. Natl. Acad. Sci. USA 2022, 119, e2118285119. [Google Scholar] [CrossRef] [PubMed]
  42. Zhao, Y.; Yang, J.; Liao, W.; Liu, X.; Zhang, H.; Wang, S.; Wang, D.; Feng, J.; Yu, L.; Zhu, W.-G. Cytosolic FoxO1 is essential for the induction of autophagy and tumour suppressor activity. Nat. Cell Biol. 2010, 12, 665–675. [Google Scholar] [CrossRef]
  43. Collier, J.J.; Suomi, F.; Oláhová, M.; McWilliams, T.G.; Taylor, R.W. Emerging roles of ATG7 in human health and disease. EMBO Mol. Med. 2021, 13, e14824. [Google Scholar] [CrossRef]
  44. Yue, J.; Aobulikasimu, A.; Sun, W.; Liu, S.; Xie, W.; Sun, W. Targeted regulation of FoxO1 in chondrocytes prevents age-related osteoarthritis via autophagy mechanism. J. Cell Mol. Med. 2022, 26, 3075–3082. [Google Scholar] [CrossRef] [PubMed]
  45. Huang, C.; Lin, Z.; Liu, X.; Ding, Q.; Cai, J.; Zhang, Z.; Rose, P.; Zhu, Y.Z. HDAC4 Inhibitors as Antivascular Senescence Therapeutics. Oxidative Med. Cell Longev. 2022, 2022, 3087916. [Google Scholar] [CrossRef] [PubMed]
  46. Hou, F.; Wei, W.; Qin, X.; Liang, J.; Han, S.; Han, A.; Kong, Q. The posttranslational modification of HDAC4 in cell biology: Mechanisms and potential targets. J. Cell Biochem. 2019, 121, 930–937. [Google Scholar] [CrossRef]
  47. Wani, A.; Al Rihani, S.B.; Sharma, A.; Weadick, B.; Govindarajan, R.; Khan, S.U.; Sharma, P.R.; Dogra, A.; Nandi, U.; Reddy, C.N.; et al. Crocetin promotes clearance of amyloid-β by inducing autophagy via the STK11/LKB1-mediated AMPK pathway. Autophagy 2021, 17, 3813–3832. [Google Scholar] [CrossRef]
  48. Liu, C.; Wang, X.; Wang, X.; Zhang, Y.; Min, W.; Yu, P.; Miao, J.; Shen, W.; Chen, S.; Zhou, S.; et al. A new LKB1 activator, piericidin analogue S14, retards renal fibrosis through promoting autophagy and mitochondrial homeostasis in renal tubular epithelial cells. Theranostics 2022, 12, 7158–7179. [Google Scholar] [CrossRef]
  49. Lu, J.; Zhang, J.; Jiang, H.; Hu, Z.; He, L.; Yang, J.; Xie, Y.; Wu, D.; Li, H.-Y.; Zeng, K.; et al. Vangl2 suppresses NF-κB signaling and ameliorates sepsis by targeting p65 for ndp52-mediated autophagic degradation. bioRxiv 2023. [Google Scholar] [CrossRef]
  50. Verzella, D.; Pescatore, A.; Capece, D.; Vecchiotti, D.; Ursini, M.V.; Franzoso, G.; Alesse, E.; Zazzeroni, F. Life, death, and autophagy in cancer: NF-κB turns up everywhere. Cell Death Dis. 2020, 11, 210. [Google Scholar] [CrossRef]
  51. Kluska, M.; Woźniak, K. Natural Polyphenols as Modulators of Etoposide Anti-Cancer Activity. Int. J. Mol. Sci. 2021, 22, 6602. [Google Scholar] [CrossRef]
  52. Ehl, S.; Astigarraga, I.; Greenwood, T.v.B.; Hines, M.; Horne, A.; Ishii, E.; Janka, G.; Jordan, M.B.; La Rosée, P.; Lehmberg, K.; et al. Recommendations for the Use of Etoposide-Based Therapy and Bone Marrow Transplantation for the Treatment of HLH: Consensus Statements by the HLH Steering Committee of the Histiocyte Society. J. Allergy Clin. Immunol. Pract. 2018, 6, 1508–1517. [Google Scholar] [CrossRef] [PubMed]
  53. Jin, L.; Song, Z.; Cai, F.; Ruan, L.; Jiang, R. Chemistry and Biological Activities of Naturally Occurring and Structurally Modified Podophyllotoxins. Molecules 2022, 28, 302. [Google Scholar] [CrossRef]
  54. Ward, D.B.; Brown, K.C.; Valentovic, M.A. Radiocontrast Agent Diatrizoic Acid Induces Mitophagy and Oxidative Stress via Calcium Dysregulation. Int. J. Mol. Sci. 2019, 20, 4074. [Google Scholar] [CrossRef] [PubMed]
  55. Yanina, I.Y.; Tanikawa, Y.; Genina, E.A.; Dyachenko, P.A.; Tuchina, D.K.; Bashkatov, A.N.; Dolotov, L.E.; Tarakanchikova, Y.V.; Terentuk, G.S.; Navolokin, N.A.; et al. Immersion optical clearing of adipose tissue in rats: Ex vivo and in vivo studies. J. Biophotonics 2022, 15, e202100393. [Google Scholar] [CrossRef] [PubMed]
  56. Tuchin, V.V.; Genina, E.A.; Tuchina, E.S.; Svetlakova, A.V.; Svenskaya, Y.I. Optical clearing of tissues: Issues of antimicrobial phototherapy and drug delivery. Adv. Drug Deliv. Rev. 2022, 180, 114037. [Google Scholar] [CrossRef]
  57. Wu, R.; Geng, L.; Zhao, Z.; Liao, D.; He, B.; Hu, H.; Lin, Y.; Li, M.; Xiang, M.; Zhang, Y.; et al. Clinical Application of Oral Meglumine Diatrizoate Esophagogram in Screening for Esophageal Fistula During Radiotherapy or Chemoradiotherapy for Esophageal Cancer. Front. Oncol. 2020, 10, 562147. [Google Scholar] [CrossRef]
  58. Guo, Z.; Chilufya, M.M.; Deng, H.; Qiao, L.; Liu, J.; Xiao, X.; Zhao, Y.; Lin, X.; Liu, H.; Xiang, R.; et al. Single and Combined Effects of Short-Term Selenium Deficiency and T-2 Toxin-Induced Kidney Pathological Injury through the MMPs/TIMPs System. Biol. Trace Element Res. 2023, 201, 4850–4860. [Google Scholar] [CrossRef]
  59. Chen, J.-H.; Xue, S.; Li, S.; Wang, Z.-L.; Yang, H.; Wang, W.; Song, D.; Zhou, X.; Chen, C. Oxidant damage in Kashin-Beck disease and a rat Kashin-Beck disease model by employing T-2 toxin treatment under selenium deficient conditions. J. Orthop. Res. 2012, 30, 1229–1237. [Google Scholar] [CrossRef]
  60. Shi, Y.; Tu, H.; Chen, X.; Zhang, Y.; Chen, L.; Liu, Z.; Sheng, J.; Han, S.; Yin, J.; Peng, B.; et al. The long non-coding RNA expression profile of Coxsackievirus A16 infected RD cells identified by RNA-seq. Virol. Sin. 2016, 31, 131–141. [Google Scholar] [CrossRef] [PubMed]
  61. Homma, K.; Suzuki, K.; Sugawara, H. The Autophagy Database: An all-inclusive information resource on autophagy that provides nourishment for research. Nucleic Acids Res. 2011, 39, D986–D990. [Google Scholar] [CrossRef] [PubMed]
  62. Liu, B.; Fu, X.; Du, Y.; Feng, Z.; Liu, X.; Li, Z.; Yu, F.; Zhou, G.; Ba, Y. In Silico Analysis of Ferroptosis-Related Genes and Its Implication in Drug Prediction against Fluorosis. Int. J. Mol. Sci. 2023, 24, 4221. [Google Scholar] [CrossRef] [PubMed]
  63. Doncheva, N.T.; Morris, J.H.; Gorodkin, J.; Jensen, L.J. Cytoscape StringApp: Network Analysis and Visualization of Proteomics Data. J. Proteome Res. 2019, 18, 623–632. [Google Scholar] [CrossRef] [PubMed]
  64. Wang, Y.; Zhang, X.; Duan, M.; Zhang, C.; Wang, K.; Feng, L.; Song, L.; Wu, S.; Chen, X. Identification of Potential Biomarkers Associated with Acute Myocardial Infarction by Weighted Gene Coexpression Network Analysis. Oxidative Med. Cell Longev. 2021, 2021, 5553811. [Google Scholar] [CrossRef]
  65. Kuleshov, M.V.; Jones, M.R.; Rouillard, A.D.; Fernandez, N.F.; Duan, Q.; Wang, Z.; Koplev, S.; Jenkins, S.L.; Jagodnik, K.M.; Lachmann, A.; et al. Enrichr: A comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016, 44, W90–W97. [Google Scholar] [CrossRef]
  66. Gu, Y.; Ma, X.; Li, J.; Ma, Y.; Zhang, Y. Identification of candidate targets for the diagnosis and treatment of atherosclerosis by bioinformatics analysis. Am. J. Transl. Res. 2021, 13, 4137–4151. [Google Scholar]
Ijms 24 13703 g001
Figure 1. The safranin O/fast green staining of articular cartilage in the T-2 group and the control group; the cartilage tissue was dyed red and the bone tissue was dyed green. The black dashed box represents partially enlarged range; the black solid box represents the empty lacunae increased range. The arrows represent the chondrocytes.
Figure 1. The safranin O/fast green staining of articular cartilage in the T-2 group and the control group; the cartilage tissue was dyed red and the bone tissue was dyed green. The black dashed box represents partially enlarged range; the black solid box represents the empty lacunae increased range. The arrows represent the chondrocytes.
Ijms 24 13703 g001
Ijms 24 13703 g002
Figure 2. Differentially expressed lncRNAs (DE-lncRNAs) and hierarchical cluster analysis. (A) Volcano plot of DE-lncRNAs in the T-2 group compared with the control group. (B) Heatmap of DE-lncRNAs in the T-2 group and the control group. Red represents upregulated lncRNAs, green represents downregulated lncRNAs, and blue represents not differentially expressed lncRNAs. The dotted line represents the threshold of DE-lncRNAs with p-value < 0.05 and | log2 (fold change) | > 1, n = 3. The color scale-bar from green to red indicates the increased expression levels of lncRNAs.
Figure 2. Differentially expressed lncRNAs (DE-lncRNAs) and hierarchical cluster analysis. (A) Volcano plot of DE-lncRNAs in the T-2 group compared with the control group. (B) Heatmap of DE-lncRNAs in the T-2 group and the control group. Red represents upregulated lncRNAs, green represents downregulated lncRNAs, and blue represents not differentially expressed lncRNAs. The dotted line represents the threshold of DE-lncRNAs with p-value < 0.05 and | log2 (fold change) | > 1, n = 3. The color scale-bar from green to red indicates the increased expression levels of lncRNAs.
Ijms 24 13703 g002
Ijms 24 13703 g003
Figure 3. GO and KEGG pathway enrichment analyses. (A) The top 10 biological process (BP), molecular function (CC), and molecular function (MF) terms of the target genes in the T-2 group vs. the control group. (B) The top 20 enriched pathways of the target genes in the T-2 group compared with the control group. Pathways with p-value < 0.05 were considered significant. Enrichment represents -log10 (p-value).
Figure 3. GO and KEGG pathway enrichment analyses. (A) The top 10 biological process (BP), molecular function (CC), and molecular function (MF) terms of the target genes in the T-2 group vs. the control group. (B) The top 20 enriched pathways of the target genes in the T-2 group compared with the control group. Pathways with p-value < 0.05 were considered significant. Enrichment represents -log10 (p-value).
Ijms 24 13703 g003
Ijms 24 13703 g004
Figure 4. Immunofluorescence staining of LC3 protein in the articular cartilage. Red represents the LC3 protein and blue represents DAPI. Scale bar is 20 μm. Data are represented as mean ± SEM (n = 3), * p-value < 0.05.
Figure 4. Immunofluorescence staining of LC3 protein in the articular cartilage. Red represents the LC3 protein and blue represents DAPI. Scale bar is 20 μm. Data are represented as mean ± SEM (n = 3), * p-value < 0.05.
Ijms 24 13703 g004
Ijms 24 13703 g005
Figure 5. Construction of PPI network and the validation of hub genes. (A) The PPI network among the autophagy-related target genes was constructed using STRING database and presented via Cytoscape. (B) The top five hub genes identified via the Maximal Clique Centrality algorithm. (C) Immunohistochemistry images of the five hub proteins in the bone marrow and skeletal muscle tissues from the Human Protein Atlas database. Scale bar is 200 μm.
Figure 5. Construction of PPI network and the validation of hub genes. (A) The PPI network among the autophagy-related target genes was constructed using STRING database and presented via Cytoscape. (B) The top five hub genes identified via the Maximal Clique Centrality algorithm. (C) Immunohistochemistry images of the five hub proteins in the bone marrow and skeletal muscle tissues from the Human Protein Atlas database. Scale bar is 200 μm.
Ijms 24 13703 g005
Ijms 24 13703 g006
Figure 6. Validation of RNA-Seq results using real-time quantitative PCR(RT-qPCR). The relative expression levels of NONRATG026881.1 (A) and NONRATG013414.2 (B) in the T-2 group vs. the control group. (C) The comparison of two DE-lncRNAs between RNA-Seq and RT-qPCR. Data are represented as mean ± SEM (n = 3).
Figure 6. Validation of RNA-Seq results using real-time quantitative PCR(RT-qPCR). The relative expression levels of NONRATG026881.1 (A) and NONRATG013414.2 (B) in the T-2 group vs. the control group. (C) The comparison of two DE-lncRNAs between RNA-Seq and RT-qPCR. Data are represented as mean ± SEM (n = 3).
Ijms 24 13703 g006
Ijms 24 13703 g007
Figure 7. The specific autophagy regulatory mechanism of articular cartilage induced by T-2 toxin. The lncRNA NONRATG014223.2 regulated Foxo3, activating ULK1, PI3CA-dependent autophagy. The lncRNAs of NONRATG024691.2, and NONRATG012484.2 severally regulated the target genes of Hdac4 and Foxo1, driving autophagy through acting on ATGs, respectively. The lncRNAs of NONRATG021591.2 and NONRATG002808.2 regulated target genes of Stk11 and Rela, and then acting on AMPK and NF-κB signaling to activate autophagy, respectively.
Figure 7. The specific autophagy regulatory mechanism of articular cartilage induced by T-2 toxin. The lncRNA NONRATG014223.2 regulated Foxo3, activating ULK1, PI3CA-dependent autophagy. The lncRNAs of NONRATG024691.2, and NONRATG012484.2 severally regulated the target genes of Hdac4 and Foxo1, driving autophagy through acting on ATGs, respectively. The lncRNAs of NONRATG021591.2 and NONRATG002808.2 regulated target genes of Stk11 and Rela, and then acting on AMPK and NF-κB signaling to activate autophagy, respectively.
Ijms 24 13703 g007
Table 1. The total quality and mapped results of RNA-Seq.
Table 1. The total quality and mapped results of RNA-Seq.
SamplesRaw DataClean DataClean Reads (%)Total Mapped Reads (%)
ReadsTotal Raw Bases (Gb)Q30 (%)ReadsTotal Clean Bases (Gb)Q30 (%)
T-2_147,053,79513.095.146,831,50212.995.499.590.4
T-2_250,495,16914.094.950,276,94214.095.199.690.8
T-2_351,064,32314.295.150,818,84014.195.499.590.7
control_149,245,93913.795.549,020,80613.695.799.590.2
control_247,073,83113.195.446,802,40513.095.799.489.7
control_345,589,58012.695.445,383,77312.695.599.590.0
Table 2. Prediction of the top 10 candidate drugs.
Table 2. Prediction of the top 10 candidate drugs.
Name of Drugsp-ValueAdjusted p-ValueGenes
etoposide CTD 000059481.37 × 10−67.65 × 10−4HDAC4, FOXO3, FOXO1, RELA
diatrizoic acid CTD 000057872.75 × 10−67.68 × 10−4FOXO3, RELA
mollugin CTD 000020587.64 × 10−67.99 × 10−4STK11, RELA
metformin BOSS8.19 × 10−67.99 × 10−4STK11, FOXO3, FOXO1
2-Butanone BOSS9.14 × 10−67.99 × 10−4STK11, FOXO3, FOXO1
rapamycin BOSS9.42 × 10−67.99 × 10−4STK11, FOXO3, FOXO1
dexamethasone BOSS1.00 × 10−57.99 × 10−4FOXO3, FOXO1, RELA
resveratrol BOSS1.38 × 10−58.81 × 10−4FOXO3, FOXO1, RELA
LY 294002 CTD 000030611.42 × 10−58.81 × 10−4FOXO3, FOXO1, RELA
rapamycin CTD 000073501.83 × 10−51.02 × 10−3FOXO3, FOXO1, RELA
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yu, F.; Wang, M.; Luo, K.; Sun, L.; Yu, S.; Zuo, J.; Wang, Y. Expression Profiles of Long Non-Coding RNAs in the Articular Cartilage of Rats Exposed to T-2 Toxin. Int. J. Mol. Sci. 2023, 24, 13703. https://doi.org/10.3390/ijms241813703

AMA Style

Yu F, Wang M, Luo K, Sun L, Yu S, Zuo J, Wang Y. Expression Profiles of Long Non-Coding RNAs in the Articular Cartilage of Rats Exposed to T-2 Toxin. International Journal of Molecular Sciences. 2023; 24(18):13703. https://doi.org/10.3390/ijms241813703

Chicago/Turabian Style

Yu, Fangfang, Miao Wang, Kangting Luo, Lei Sun, Shuiyuan Yu, Juan Zuo, and Yanjie Wang. 2023. "Expression Profiles of Long Non-Coding RNAs in the Articular Cartilage of Rats Exposed to T-2 Toxin" International Journal of Molecular Sciences 24, no. 18: 13703. https://doi.org/10.3390/ijms241813703

APA Style

Yu, F., Wang, M., Luo, K., Sun, L., Yu, S., Zuo, J., & Wang, Y. (2023). Expression Profiles of Long Non-Coding RNAs in the Articular Cartilage of Rats Exposed to T-2 Toxin. International Journal of Molecular Sciences, 24(18), 13703. https://doi.org/10.3390/ijms241813703

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop