Comprehensive CircRNA Profiling and Selection of Key CircRNAs Reveal the Potential Regulatory Roles of CircRNAs throughout Ovarian Development and Maturation in Cynoglossus semilaevis

SIMPLE SUMMARY: CircRNAs: as molecules involved in gene regulation, have become a new research hotspot in the non-coding RNA field. CircRNAs show tissue- or developmental stage-specific patterns of expression and can influence the expression levels of their parental genes. Recent studies have documented the potential biological roles of circRNAs in the growth, development, reproduction and health of humans and animals. Tongue sole (Cynoglossus semilaevis) is a marine flatfish that is an economically important farmed species in China. The commercial aquaculture of tongue sole has developed in the last few years because wild resources have decreased. Reproduction is regulated by brain-pituitary-gonad-liver axis which limits the development of artificial tongue sole culture. However, the roles of circRNAs in the ovarian development and maturation of tongue sole has never been reported. The identification of the potential functions of circRNAs provides a foundation for understanding the genetic mechanisms that regulate oocyte growth and maturation, which will allow the efficiency of tongue sole reproduction to be improved. Moreover, our findings extend the knowledge about a new type of endogenous RNA involved in regulating the ovarian development and maturation of tongue sole. ABSTRACT: CircRNAs are novel endogenous non-coding small RNAs involved in the regulation of multiple biological processes. However, little is known regarding circRNAs in ovarian development and maturation of fish. Our study, for the first time, provides the genome-wide overview of the types and relative abundances of circRNAs in tongue sole tissues during three ovarian developmental stages. We detected 6790 circRNAs in the brain, 5712 in the pituitary gland, 4937 in the ovary and 4160 in the liver. Some circRNAs exhibit tissue-specific expression, and qRT-PCR largely confirmed 6 differentially expressed (DE) circRNAs. Gene Ontology and KEGG pathway analyses of DE mRNAs were performed. Some DE circRNA parental genes were closely associated with biological processes in key signalling pathways and may play essential roles in ovarian development and maturation. We found that the selected circRNAs were involved in 10 pathways. RNase R digestion experiment and Sanger sequencing verified that the circRNA had a ring structure and was RNase R resistant. qRT-PCR results largely confirmed differential circRNA expression patterns from the RNA-seq data. These findings indicate that circRNAs are widespread in terms of present in production-related tissues of tongue sole with potentially important regulatory roles in ovarian development and maturation.