Landscape of long non-coding RNAs in Trichophyton mentagrophytes-induced rabbit dermatophytosis lesional skin and normal skin

Emerging evidences suggest that long non-coding RNAs (lncRNAs) play important role in disease development. However, the role of rabbit lncRNAs in the pathogenesis of dermatophytosis remains elusive. The present study aimed to study and characterize lncRNA transcriptome in 8 T. mentagrophytes-induced female rabbit dermatophytosis lesional (TM) and 4 normal saline-infected (NS) skin biopsies using RNAseq. We identified 5883 lncRNAs in 12 strand-specific RNA-seq libraries and found 64 differentially expressed lncRNAs (q < 0.05) in TM relative to NS. As in other mammalian counterparts, rabbit lncRNAs were distributed in all chromosomes except the Y chromosome and were generally smaller in size and fewer in exon numbers compared to protein coding genes. Next, co-expression analysis revealed that 107 pairs between 32 DE lncRNAs and 96 protein coding genes showed a highly correlated expression (|r| > 0.8). Moreover, miRPara analysis of the lncRNAs revealed 173 lncRNAs with precursor sequences for 9561 probable novel miRNAs. Finally, q-PCR results validated the RNA-seq results with eight randomly selected lncRNAs. To the best of our knowledge, this is the first report on rabbit lncRNAs, and our results highlighted the potential role of lncRNAs in the pathogenesis of dermatophytosis.