Long noncoding RNAs profiling in ovary during laying and nesting in Muscovy ducks (Cairina moschata)

There are recent reports of the important functions of long noncoding RNAs (lncRNAs) in female reproductive and ovarian development. Studies in which there was characterization of lncRNAs in the ovaries of laying compared with nesting poultry, however, are limited. In this study, RNA libraries were constructed by obtaining sequencing data of ovarian tissues from laying and nesting Muscovy ducks. In the ovarian tissues of Muscovy ducks, a total of 334 differentially abundant mRNA transcripts (DEGs) and 36 differentially abundant lncRNA transcripts were identified in the nesting period, when compared with during the laying period. These results were subsequently validated by qRT-PCR using nine randomly-selected lncRNAs and six randomly-selected DAMTs. Furthermore, the cis- and trans-regulatory target genes of differentially abundant lncRNA transcripts were identified, and lncRNA-gene interaction networks of 34 differentially abundant lncRNAs and 263 DEGs were constructed. A total of 7601 lncRNAs neighboring 10,542 protein-coding genes were identified and found to be enriched in the Wnt signaling pathway and oocyte meiosis pathways associated with follicular development. Overall, only 11 cis-targets and 57 mRNA-mRNA except trans-targets were involved in the lncRNA-gene interaction networks. Based on the interaction networks, nine DEGs were trans-regulated by differentially abundant lncRNAs and 20 differentially abundant lncRNAs were hypothesized to have important functions in the regulation of broodiness in Muscovy ducks. In this study, a predicted interaction network of differentially abundant lncRNAs and DEGs in Muscovy ducks was constructed for the first time leading to an enhanced understanding of lncRNA and gene interactions regulating broodiness.