BackgroundParatuberculosis (Johne’s disease, JD) is a chronic and enteric disease in a range of ruminants, often caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection, leading to substantial economic losses worldwide. Yet, the molecular underpinning of paratuberculosis remains elusive. Here we performed RNA sequencing (RNA-seq) and small RNA sequencing (sRNA-seq) of the jejunum tissues from the Holstein cows with three distinct statuses of paratuberculosis, i.e., healthy, subclinical, and clinical to screen potential genes, lncRNAs, and miRNAs associated with the resistance or susceptibility to MAP infection and build ceRNA regulatory networks via miRNAs.ResultsWe applied whole transcriptome sequencing analysis to examine the jejunum tissue in nine Holstein cows. Starting with 19,994 expressed genes, 13,529 lncRNAs, and 735 miRNAs, we screened out differentially expressed genes (DEGs), lncRNAs, and miRNAs of three comparison groups, i.e., clinical vs. healthy, subclinical vs. healthy, and clinical vs. subclinical, subsequently identifying ceRNA pairs. Ultimately, we detected 76, 74, and 24 DEGs, 19, 39, and 10 lncRNAs, as well as 28, 61, and 20 miRNAs across the three comparison groups, respectively. Through integrating these DEGs with functional annotation, previously reported QTLs, and GWAS results, we proposed eight genes (LYZ, LYZ1, BOLA-DQB, BOLA-DQA1, TAP, CATD, VNN1, and PPARG), six lncRNAs, 48 miRNAs, and 107 ceRNA pairs implying their potential associations with susceptibility to MAP infection.ConclusionThe present study provided a global view of the dynamics in transcriptomes of the bovine jejunum tissues in terms of JD status. Our results demonstrated that not only mRNAs but also lncRNAs and miRNAs played important roles in regulating MAP infection in dairy cattle. This study provided a valuable resource for understanding the molecular basis of JD, potentially contributing to the genetic improvement of JD resistance in dairy cattle.
Read full abstract