Rtal genesResultsIdentification of ovary-derived circRNAs through the onset of pubertyTo receive a dependable result, CIRI2 and find_circ software have been intersected to determine circRNAs. A total of 972 circRNAs candidates have been identified within the pubertal transition of pubertal ovaries (Further file 1) (Fig. 1a). Thereinto, 347, 293 and 827 circRNAs were identified in pre-, in- and post-puberty, respectively (Fig. 1b). These circRNAs have been much more distributed in Sus scrofa chromosome 1 (Fig. 1c). The expressions of circRNAs were the lowest inside the P/Q-type calcium channel custom synthesis post-puberty (Fig. 1d), compared to preand in-puberty. Meanwhile, the average genome distanceTo further investigate the biological functions of circRNAs involved in pubertal ovary, the parental genes of all 972 circRNAs in puberty were made use of for the analysis of KEGG pathway (Further file two). Clearly, the functional pathways have been considerably over-represented in pubertal ovaries, which includes steroid biosynthesis, autophagy-animal, MAPK signaling pathway, progesterone-mediated oocyte maturation and ras signaling pathway (Fig. 2a). In the steroid biosynthesis signaling pathway, “circ 14:149834024992789”, “14: 149845014992789” and “14:149892705001827” derived from FDFT1 gene were uniquely expressed in the post-puberty (Fig. 2b). In the autophagy-animal signaling pathway, “circ 1:19064825390654424” derived from HIF1A gene was uniquely expressed inside the in-puberty. In the MAPK signaling pathway, “circ 15:2512858125143159” and “circ 15:251403155143159” derived from MAP3K2 gene was uniquely expressed in the postpuberty. Within the progesterone-mediated oocyte maturation, “circ 16:504376150478728” derived from CPEB4 gene was uniquely expressed within the post-puberty. Within the ras signaling pathway, “circ 12:4367306943691261” derived from NF1 gene was expressed inside the in-puberty and post-puberty. Specifics of those circRNAs were shown in Additional file three. In addition, so as to further discover the circRNAs in pubertal genes, 10 puberty-related genes (ESR1, JAK2, NF1, ARNT, IGF1, KISS1, Gpr54, NKB, Mkrn3, GnRH) have been selected andPan et al. BMC Genomics(2021) 22:Page three ofFig. 1 Overview of the identified circRNAs by ULK2 Species RNA-seq analyses in ovaries of gilts. a CircRNAs have been identified by two algorithms. b The Venn diagram shows the number of distinctive and frequent circRNAs in pre-, in- and post-puberty. c Circos plot from the circRNAs distribution within the whole genome of gilts. The outermost circle represents the distribution of your number of circRNAs, the blue circle represents the distribution of expression amount of circRNAs, and also the red inner circle represents the length of circRNAs. d Expression amount of circRNAs in three stages, p 0.05. e Genomic dance of all detected circRNAs. f Transcript length of circRNAs. g Proportion of three kinds and exon number of the circRNAs in 3 stages. h The upset plot of three forms of circRNAs as well as the corresponding parental genesexplored by manual reviewing the literature and databases, and 5 circRNAs derived from 4 pubertal genes (ESR1, JAK2, NF1, ARNT) had been lastly discovered in the ovary of pubertal transition (Further file four). “circ 1: 143732324374308” (uniquely expressed within the prepuberty) and “circ 1:144163354457143” (expressed within the pre-, in- and post-puberty) were derived from ESR1; “circ 1:21691427516951002” (uniquely expressed inside the post-puberty) was derived from JAK2; “circ 12: 436730693691261” (expressed within the in- and postpuberty) was derived from NF1; “circ four:9836952098372553” (uniqu.