Analysis of circRNA-miRNA-mRNA regulatory network of embryonic gonadal development in Mulard duck

Abstract

The aim of the study was to explore the regulatory mechanism of differences in embryonic gonadal development between intergeneric distance hybrid offspring Mulard ducks and parent ducks. The morphological differences gonadal tissues of Muscovy ducks, Pekin ducks and Mulard ducks at 12.5-day embryonic age were observed by sectioning and hematoxylin-eosin (HE) staining. Then followed by transcriptome sequencing to screen for gonadal development-related differentially expressed circRNAs and mRNAs to construct a competitive endogenous RNA (ceRNA) regulatory network. Finally, qRT-PCR and luciferase reporter system were used to verify the sequencing data and targeting relationship of ceRNA pairs. The results showed that the seminiferous tubule lumen of Mulard ducks was not obvious, while there were obvious seminiferous tubules and tubular structures in testis of Pekin ducks and Muscovy ducks, with number and shape indicating maturity. There were 18 upregulated circRNAs and 16 downregulated circRNAs in Mulard ducks and Pekin ducks, respectively, and 39 upregulated circRNAs and 1 downregulated circRNA in Mulard ducks and Muscovy ducks, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis found that genes involves in dorso-ventral axis formation, for example, neurogenic locus notch homolog protein 1 (NOTCH1), were significantly enriched (P < 0.05). The novel_circ_0002265-gga-miR-122-5p-PAFAH1B2 regulatory network was constructed. The qRT-PCR results showed that the sequencing results were reliable. The dual-luciferase reporter assay showed that gga-miR-122-5p exists binding site of circ_0002265 and PAFAH1B2, indicating circ_0002265-gga-miR-122-5p-PAFAH1B2 targeting relationship. In summary, the embryonic gonadal development of intergeneric hybrid Mulard ducks may be regulated by differentially expressed circRNAs and genes, such as novel_circ_0000519, novel_circ_0003537, NOTCH1, FGFR2, PAFAH1B1, and PAFAH1B2, among which circ_0002265-gga-miR-122-5p-PAFAH1B2 may participate in the targeted regulation of gonadal development in Mulard ducks. The findings of this study are helpful for analyzing the mechanism of embryonic gonadal development differences in avians.

Keywords: Mulard duck; ceRNA; circRNA; embryonic stage; gonad.