Responses of keratinocytes to Trichophyton mentagrophyte infection based on whole transcriptome analysis.

Abstract

Dermatophytosis is an intractable superficial mycosis in humans and animals mainly caused by Trichophyton mentagrophytes (T. mentagrophytes), with a global prevalence of about 20%. Keratinocytes are the most abundant participants in skin immunity, and they also play a role in the first-line defence against T. mentagrophytes. However, no studies of keratinocyte responses against T. mentagrophytes infection based on the whole transcriptome have been reported. Here, we systematically analysed changes in keratinocytes infected with T. mentagrophytes using whole transcriptome sequencing technology. The phenotypic changes in keratinocytes after infection with 1 × 105 conidia/mL T. mentagrophytes were observed by light microscopy, scanning electron microscopy, transmission electron microscopy and terminal deoxynucleotidyl transferase dUTP nick end labeling. RNA-sequencing (RNA-seq), small RNA-seq technology and related bioinformatics methods were used to systematically analyse the whole transcriptome changes in keratinocytes upon T. mentagrophytes stimulation. We found that T. mentagrophytes infection caused morphological changes, membrane damage, the formation of irregular organelles and keratinocyte apoptosis. A total of 204 differentially expressed (DE) circular RNAs (circRNAs), 868 DE long noncoding RNAs (lncRNAs), 2973 DE mRNAs and 209 DE micro RNAs (miRNAs) were identified between noninfected and T. mentagrophytes-infected keratinocytes. The expression level of selected RNAs was validated by quantitative real-time polymerase chain reaction (qRT-PCR). Functional enrichment analysis revealed that the parental genes of DE circRNAs were related to cell response, cell death and establishment of the skin barrier. Genes targeted by miRNA were involved in regulating the initiation of the immune response. Based on the expression level of circRNAs, lncRNAs, mRNAs and miRNAs, circRNA-miRNA-mRNA competing endogenous (ceRNA) networks comprised of 159 DE miRNAs, 141 DE circRNAs and 2307 DE mRNAs, and lncRNA-miRNA-mRNA ceRNA networks comprised of 790 DE lncRNAs, 190 DE miRNAs and 2663 DE mRNAs were constructed. The reliability of two selected ceRNA networks was verified using qRT-PCR. Further functional enrichment analysis revealed that the DE mRNAs interacting with circRNAs and lncRNAs in the ceRNA network mainly participated in fungal recognition, inflammation, the innate immune response and the death of keratinocytes. Our findings might provide new evidence on the pathogenesis of T. mentagrophytes-induced dermatophytosis, which is essential for identifying new therapeutic targets for dermatophytosis treatment.