In this study, circular RNAs (circRNAs) from Holstein cow mammary tissues were identified and compared between early lactation and non-lactation. After analysis, 10,684 circRNAs were identified, ranging from 48 to 99,406 bp, and the average size was 882 bp. The circRNAs were mainly distributed on chromosomes 1 to 11, and 89.89% of the circRNAs belonged to sense-overlapping circRNA. The exons contained with circRNAs ranged from 1 to 47 and were concentrated from 1 to 5. Compared with the non-lactating cows, 87 circRNAs were significantly differentially expressed in the peak lactation cows. There were 68 upregulated circRNAs and 19 downregulated circRNAs. Enrichment analysis of circRNAs showed that GO analysis mainly focused on immune response, triglyceride transport, T cell receptor signaling pathway, etc. Pathway analysis mainly focused on cytokine-cytokine receptor interaction, T helper 17 cell differentiation, fatty acid biosynthesis, the JAK-STAT signaling pathway, etc. Specific primers were designed for two proximal ends of the circRNA junction sites to allow for PCR validation of four randomly selected circRNAs and carry out circRNA-miRNA interaction research. This study revealed the expression profile and characteristics of circRNAs in mammary tissue from Holstein cows at early lactation and non-lactation, thus providing rich information for the study of circRNA functions and mechanisms, as well as potential candidate miRNA genes for studying lactation in Holstein cows.
Keywords: Holstein cows; RNA sequencing; circular RNA; mammary tissue.