Genome-wide analysis of circular RNAs involved in Marek's disease tumourigenesis in chickens.

Marek's disease (MD), induced by Marek's disease virus (MDV), is a lymphotropic neoplastic disease and causes huge economic losses to the poultry industry. Non-coding RNAs (ncRNAs) play important regulatory roles in disease pathogenesis. To investigate host circular RNA (circRNA) and microRNA (miRNA) expression profile, RNA sequencing was performed in tumourous spleens (TS), spleens from the survivors (SS) without any lesion after MDV infection, and non-infected chicken spleens (NS). A total of 2,169 circRNAs were identified and more than 80% of circRNAs were derived from exon. The flanking introns of 1,744 exonic circRNAs possessed 579 reverse complementary matches (RCMs), which mainly overlapped with chicken repeat 1 family (CR1F). It suggested that CR1F mediated the cyclization of exons by intron pairing. Out of 2,169 circRNAs, 113 were differentially expressed circRNAs (DECs). The Q-PCR and Rnase R digestion experiments showed circRNA possessed high stability compared with their linear RNAs. Integrated with previous transcriptome data, we profiled regulatory networks of circRNA/long non-coding RNA (lncRNA)-miRNA-mRNA. Extensive competing endogenous RNA (ceRNA) networks were predicted to be involved in MD tumourigenesis. Interestingly, circZMYM3, an intronic circRNA, interacted with seven miRNAs which targeted some immune genes, such as SWAP70 and CCL4. Gga-miR-155 not only interacted with circGTDC1 and circMYO1B, but also targeted immune-related genes, such as GATA4, which indicated the roles of non-coding RNAs played to mediate immune responsive genes. Collectively, this is the first study that integrated RNA expression profiles in MD model. Our results provided comprehensive interactions of ncRNAs and mRNA in MD tumourigenesis.