Identification and differential expression of microRNAs in ovaries of laying and Broody geese (Anser cygnoides) by Solexa sequencing

Comprehensive analysis of differentially expressed profiles of mRNA, lncRNA, and miRNA of Yili geese ovary at different egg-laying stages

Background

The development of the ovaries is an important factor that affects egg production performance in geese. Ovarian development is regulated by genes that are expressed dynamically and stage-specifically. The transcriptome profile analysis on ovarian tissues of goose at different egg laying stages could provide an important basis for screening and identifying key genes regulating ovarian development.

Results

In this study, 4 ovary tissues at each breeding period of pre-laying (PP), laying (LP), and ceased-laying period (CP), respectively, with significant morphology difference, were used for RNA extraction and mRNAs, lncRNAs, and miRNAs comparison in Yili geese. CeRNA regulatory network was constructed for key genes screening. A total of 337, 1136, and 525 differentially expressed DE mRNAs, 466, 925, and 742 DE lncRNAs and 258, 1131 and 909 DE miRNAs were identified between PP and LP, between CP and LP, and between CP and PP groups, respectively. Functional enrichment analysis showed that the differentially expressed mRNAs and non-coding RNA target genes were mainly involved in the cell process, cytokine-cytokine receptor interaction, phagosome, calcium signaling pathway, steroid biosynthesis and ECM-receptor interaction. Differential genes and non-coding RNAs, PDGFRB, ERBB4, LHCGR, MSTRG.129094.34, MSTRG.3524.1 and gga-miR-145–5p, related to reproduction and ovarian development were highly enriched. Furthermore, lncRNA-miRNA-mRNA regulatory networks related to ovary development were constructed.

Conclusions

Our study found dramatic transcriptomic differences in ovaries of Yili geese at different egg-laying stages, and a differential lncRNA-miRNA-mRNA regulatory network related to cell proliferation, differentiation and apoptosis and involved in stromal follicle development were established and preliminarily validated, which could be regarded as a key regulatory pathway of ovarian development in Yili geese.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08774-4.

Bone marrow-derived mesenchymal stem cells combined with gonadotropin therapy restore postnatal oogenesis of chemo-ablated ovaries in rats via enhancing very small embryonic-like stem cells

BACKGROUND

Very small embryonic-like stem cells (VSELs) are a rare population within the ovarian epithelial surface. They contribute to postnatal oogenesis as they have the ability to generate immature oocytes and resist the chemotherapy. These cells express markers of pluripotent embryonic and primordial germ cells.

OBJECTIVE

We aimed to explore the capability of VSELs in restoring the postnatal oogenesis of chemo-ablated rat ovaries treated with bone marrow-derived mesenchymal stem cells (BM-MSCs) combined with pregnant mare serum gonadotropin (PMSG).

METHODS

Female albino rats were randomly assigned across five groups: I (control), II (chemo-ablation), III (chemo-ablation + PMSG), IV (chemo-ablation + MSCs), and V (chemo-ablation + PMSG + MSCs). Postnatal oogenesis was assessed through measurement of OCT4, OCT4A, Scp3, Mvh, Nobox, Dazl4, Nanog, Sca-1, FSHr, STRA8, Bax, miR143, and miR376a transcript levels using qRT-PCR. Expression of selected key proteins were established as further confirmation of transcript expression changes. Histopathological examination and ovarian hormonal assessment were determined.

RESULTS

Group V displayed significant upregulation of all measured genes when compared with group II, III or IV. Protein expression confirmed the changes in transcript levels as group V displayed the highest average density in all targeted proteins. These results were confirmed histologically by the presence of cuboidal germinal epithelium, numerous primordial, unilaminar, and mature Graafian follicles in group V.

CONCLUSION

VSELs can restore the postnatal oogenesis in chemo-ablated ovaries treated by BM-MSCs combined with PMSG.

Novel miRNA identification and comparative profiling of miRNA regulations revealed important pathways in Jinding duck ovaries by small RNA sequencing

Functional studies have revealed miRNAs play pivotal roles in ovulation and ovary development in mammalians, whereas little is known about the miRNA function in ducks. In this study, miRNA deep sequencing in the ovary tissues was carried out to obtain the miRNA profile from ovaries before oviposition (BO) and after oviposition (AO) in Jinding duck. Overall, an average of 23,128,075 and 26,020,523 reads were identified in the BO and AO samples, respectively, and 6739 miRNAs were identified from them through further mapping and analysis. Besides, 1570 miRNAs were identified as differentially expressed miRNAs compared with BO, including 493 miRNAs up-regulated and 1077 down-regulated in AO. Moreover, 2291 target genes were predicted from 443 significantly differentially expressed miRNAs. In addition, GO and KEGG pathway analysis indicated that target genes were enriched in some basic cell metabolism pathways as well as the productive pathways such as MAPK signaling pathway, gonadotropin-releasing hormone signaling pathway, TGF-beta signaling pathway which had been significantly changed. Our results helped to replenish the duck miRNA database and illustrate the potential mechanism of miRNA function in duck ovary development and reproduction process.

Selection of the Most Stable Endogenous Control Genes for Microrna Quantitation in Chicken Ovarian Follicles

MicroRNAs (miRNAs or miRs) belong to a class of small non-coding RNAs of 19 to 24 nucleotides long that act as negative gene regulators at the post-transcriptional level. Quantitative PCR (q-PCR) is a commonly used technique in the profiling of miRs, and identification of reliable endogenous controls is crucial for proper data normalisation. To date, no study has been performed on reference miRs for the normalisation of miR expression in chicken ovarian tissues. Therefore, the aim of the present study was to experimentally identify the most stable reference mirs for normalisation of miR q-PCR expression data in the chicken ovary. Relying on high-throughput sequencing, five putative reference miR (let-7a-3p, miR-140a-3p, miR-22-5p, miR-33-5p, miR-99a-3p) were identified and subsequently analysed in a total of 66 tissue samples. The stability of candidate endogenous controls validated by the most widely used algorithms, geNorm, NormFinder, and BestKeeper, showed that let-7a-3p, miR-140a-3p, and miR-22-5p are the most appropriate choice of reference genes. Application of different normalisation approaches to the relative quantitation of randomly chosen miR-1552-5p in chicken ovarian follicles indicated the impact of the selected reference genes on miR expression. Further, the results revealed a downregulation of miR-1552-5p. In summary, the three identified endogenous reference miRs are suitable for profiling the miR expression in ovarian tissues of laying hens. Our findings provide valuable information for future miR expression studies in the avian ovary.

Identification and comparative profiling of gonadal microRNAs in the adult pigeon (Columba livia)

1. MicroRNAs are small noncoding RNA molecules that play crucial roles in gene expression. However, the comparative profiling of testicular and ovarian microRNAs in birds are rarely reported, particularly in pigeon.2. In this study, Illumina next-generation sequencing technology was used to sequence miRNA libraries of the gonads from six healthy adult utility pigeons. A total of 344 conserved known miRNAs and 32 novel putative miRNAs candidates were detected. Compared with those of ovaries, 130 differentially expressed (DE) miRNAs were identified in the testes. Among them, 70 miRNAs showed down-regulation in the ovaries, while another 60 miRNAs were up-regulated.3. Combining the results of the expression of target gene measurements and pathway enrichment analyses, it was revealed that some DEmiRNAs from the gonad samples involved in sexual differentiation and development (such as cli-miR-210-3p and cli-miR-214-3p) could down-regulate AR (androgen receptor). Cli-miR-181b-5p, cli-miR-9622-3p and cli-miR-145-5p were highly expressed in both the ovaries and testes, which could co-target HOXC9, and were related to regulation of primary metabolic processes. KEGG enrichment analysis showed that DEmiRNAs may play biological and sex-related roles in pigeon gonads.4. The expression profiles of testicular and ovarian miRNA in adult pigeon gonads are presented for the first time, and the findings may contribute to a better understanding of gonadal expression in poultry.

Identification and characterization of a specific 13-miRNA expression signature during follicle activation in the zebrafish ovary

Ovarian folliculogenesis is always of great interest in reproductive biology. However, the molecular mechanisms that control follicle development, particularly the early phase of follicle activation or recruitment, still remain poorly understood. In an attempt to decipher the gene networks and signaling pathways involved in such transition, we conducted a transcriptomic analysis (RNA-seq) on zebrafish primary growth (PG, stage I; inactive) and previtellogenic (PV, stage II; activated) follicles. A total of 118 unique microRNAs (miRNAs) (11 downregulated and 83 upregulated during PG/PV transition) and 56711 unique messenger RNAs (mRNAs) (1839 downregulated and 7243 upregulated during PG/PV transition) were identified. Real-time quantitative polymerase chain reaction analysis confirmed differential expression of 46 miRNAs from 66 candidates (66.67%). Among which, we chose to focus on 13 miRNAs (let-7a, -7b, -7c-5p, -7d-5p, -7h, -7i; miR-21, -23a-3p, -27c-3p, -107a-3p, -125b-5p, -145-3p, and -202-5p) that exhibited significant differential expression between PG and PV follicles (P ≤ 0.045*). With this 13-miRNA expression signature alone, PG follicles can be well differentiated from PV follicles by hierarchical clustering, suggesting their functional relevance during PG-to-PV transition. By overlaying predicted target genes and the differentially expressed mRNAs revealed by the RNA-seq analysis, especially those showing reciprocal miRNA-mRNA expression patterns, we shortlisted a panel of miRNA downstream targets for luciferase reporter validation. The reporter assay confirmed the interactions of let-7i:: atg4a (P = 0.01*), miR-202-5p::c23h20orf24 (P = 0.0004***), and miR-144-5p::ybx1 (P = 0.003**), implicating these potential miRNA-mRNA gene pairs in follicle activation during folliculogenesis. Our transcriptomic data analyses suggest that miRNA-mediated post-transcriptional control may represent an important mechanism underlying follicle activation.

Deep sequencing identified potential miRNAs involved in defence response, stress and plant growth characteristics of wild genotypes of cardamom

Cardamom has long been used as a food flavouring agent and in ayurvedic medicines for mouth ulcers, digestive problems and even depression. Extensive occurrence of pests and diseases adversely affect its cultivation and result in substantial reductions in total production and productivity. Numerous studies revealed the significant role of miRNAs in plant biotic stress responses. In the current study, miRNA profiling of cultivar and wild cardamom genotypes was performed using an Ion Proton sequencer. We identified 161 potential miRNAs representing 42 families, including monocot/tissue-specific and 14 novel miRNAs in both genotypes. Significant differences in miRNA family abundance between the libraries were observed in read frequencies. A total of 19 miRNAs (from known miRNAs) displayed a twofold difference in expression between wild and cultivar genotypes. We found 1168 unique potential targets for 40 known miRNA families in wild and 1025 potential targets for 42 known miRNA families in cultivar genotypes. The differential expression analysis revealed that most miRNAs identified were highly expressed in cultivars and, furthermore, lower expression of miR169 and higher expression of miR529 in wild cardamom proved evidence that wild genotypes have stronger drought stress tolerance and floral development than cultivars. Potential targets predicted for the newly identified miRNAs from the miRNA libraries of wild and cultivar cardamom genotypes involved in metabolic and developmental processes and in response to various stimuli. qRT-PCR confirmed miRNAs were differentially expressed between wild and cultivar genotypes. Furthermore, four target genes were validated experimentally to confirm miRNA-mRNA target pairing using RNA ligase-mediated 5' Rapid Amplification of cDNA Ends (5'RLM-RACE) PCR.

Deep sequencing identification of miRNAs in pigeon ovaries illuminated with monochromatic light

Background

The use of light of different wavelengths has grown popular in the poultry industry. An optimum wavelength is believed to improve pigeon egg production, but little is known about the role of microRNAs (miRNAs) in the effects of monochromatic light on ovarian pigeon function. Herein, we harvested ovaries from pigeons reared under monochromatic light of different wavelength and performed deep sequencing on various tissues using an Illumina Solexa high-throughput instrument.

Results

We obtained 66,148,548, 67,873,805, and 71,661,771 clean reads from ovaries of pigeons reared under red light (RL), blue light (BL), and white light (WL), respectively. We identified 1917 known miRNAs in nine libraries, of which 524 were novel. Three and five differentially expressed miRNAs were identified in BL vs. WL and RL vs. WL groups, respectively. Quantitative reverse transcription PCR was used to validate differentially expressed miRNAs (miR-200, miR-122, and miR-205b). In addition, 5824 target genes were annotated as differentially expressed miRNAs, most of which are involved in reproductive pathways including oestrogen signalling, cell cycle, and oocyte maturation. Notably, ovarian miR-205b expression was significantly negatively correlated with its target 11β-hydroxysteroid dehydrogenase type 1 (HSD11B1).

Conclusions

miRNA–mRNA network analysis suggests that miR-205b targeting of HSD11B1 plays a key role in the effects of monochromatic light on pigeon egg production. These findings indicate that monochromatic light shortens the oviposition interval of pigeons, which may be useful for egg production and pigeon breeding.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4831-6) contains supplementary material, which is available to authorized users.

Whole-transcriptome analysis of atrophic ovaries in broody chickens reveals regulatory pathways associated with proliferation and apoptosis

Broodiness in laying hens results in atrophy of the ovary and consequently decreases productivity. However, the regulatory mechanisms that drive ovary development remain elusive. Thus, we collected atrophic ovaries (AO) from 380-day-old broody chickens (BC) and normal ovaries (NO) from even-aged egg-laying hens (EH) for RNA sequencing. We identified 3,480 protein-coding transcripts that were differentially expressed (DE), including 1,719 that were down-regulated and 1,761 that were up-regulated in AO. There were 959 lncRNA transcripts that were DE, including 56 that were down-regulated and 903 that were up-regulated. Among the116 miRNAs that were DE, 79 were down-regulated and 37 were up-regulated in AO. Numerous DE protein-coding transcripts and target genes for miRNAs/lncRNAs were significantly enriched in reproductive processes, cell proliferation, and apoptosis pathways. A miRNA-intersection gene-pathway network was constructed by considering target relationships and correlation of the expression levels between ovary development-related genes and miRNAs. We also constructed a competing endogenous RNA (ceRNA) network by integrating competing relationships between protein-coding genes and lncRNA transcripts, and identified several lncRNA transcripts predicted to regulate the CASP6, CYP1B1, GADD45, MMP2, and SMAS2 genes. In conclusion, we discovered protein-coding genes, miRNAs, and lncRNA transcripts that are candidate regulators of ovary development in broody chickens.

miRNAome expression profiles in the gonads of adult Melopsittacus undulatus

The budgerigar (Melopsittacus undulatus) is one of the most widely studied parrot species, serving as an excellent animal model for behavior and neuroscience research. Until recently, it was unknown how sexual differences in the behavior, physiology, and development of organisms are regulated by differential gene expression. MicroRNAs (miRNAs) are endogenous short non-coding RNA molecules that can post-transcriptionally regulate gene expression and play a critical role in gonadal differentiation as well as early development of animals. However, very little is known about the role gonadal miRNAs play in the early development of birds. Research on the sex-biased expression of miRNAs in avian gonads are limited, and little is known about M. undulatus. In the current study, we sequenced two small non-coding RNA libraries made from the gonads of adult male and female budgerigars using Illumina paired-end sequencing technology. We obtained 254 known and 141 novel miRNAs, and randomly validated five miRNAs. Of these, three miRNAs were differentially expressed miRNAs and 18 miRNAs involved in sexual differentiation as determined by functional analysis with GO annotation and KEGG pathway analysis. In conclusion, this work is the first report of sex-biased miRNAs expression in the budgerigar, and provides additional sequences to the avian miRNAome database which will foster further functional genomic research.

Identification and characterization of microRNA in the lung tissue of pigs with different susceptibilities to PCV2 infection

Porcine circovirus type 2 (PCV2) is the primary cause of post-weaning multisystemic wasting syndrome (PMWS) and other PCV-associated diseases. According to our previous RNA-sequencing analysis, the differences in the susceptibility to PCV2 infection depended on the genetic differences between the Laiwu (LW) and Yorkshire × Landrace crossbred (YL) pigs, but the cellular microRNA (miRNA) that are differentially expressed between the LW and YL pigs before and after PCV2 infection remain to be determined. In this study, high-throughput sequencing was performed to determine the abundance and differential expression of miRNA in lung tissues from PCV2-infected and PCV2-uninfected LW and YL pigs. In total, 295 known and 95 novel miRNA were identified, and 23 known and 25 novel miRNA were significantly differentially expressed in the PCV2-infected vs. PCV2-uninfected LW pigs and/or the PCV2-infected vs. PCV2-uninfected YL pigs. The expression levels of ssc-miR-122, ssc-miR-192, ssc-miR-451, ssc-miR-486, and ssc-miR-504 were confirmed by quantitative real-time PCR (qRT-PCR). Analysis of the potential targets of the four up-regulated miRNA (i.e., ssc-miR-122, ssc-miR-192, ssc-miR-451 and ssc-miR-486) identified pathways and genes that may be important for disease resistance. Among the up-regulated miRNA, ssc-miR-122 can repress the protein expression and viral DNA replication of PCV2 and down-regulate the expression of the nuclear factor of activated T-cells 5 (NFAT5) and aminopeptidase puromycin sensitive (NPEPPS) by binding to their 3′ untranslated region (3′UTR) in PK15 cells. Therefore, ssc-miR-122 may indirectly suppress PCV2 infection by targeting genes related to the host immune system, such as NFAT5 and NPEPPS.

Transcriptome analysis revealed the possible regulatory pathways initiating female geese broodiness within the hypothalamic-pituitary-gonadal axis

Geese have the strongest tendency toward broodiness among all poultry. The mechanisms initiating broodiness within the goose hypothalamic-pituitary-gonadal axis (HPGA) are still unclear. Here, we reported the transcriptome differences between laying and initial nesting within the HPGA tissues of geese. We constructed a unigene database based on HPGA tissues and identified 128,148 unigenes, 100% of which have been annotated. By using Digital Gene Expression (DGE) sequencing, we screened 19, 110, 289, and 211 differentially expressed genes (DEGs) in the hypothalamus, pituitary gland, stroma ovarii, and follicles, respectively, between laying and nesting geese. Expression changes of hypocretin (HCRT) and pro-opiomelanocortin (POMC) in the hypothalamus of nesting geese may cause appetite reduction, which is possibly the first step and a prerequisite to initiate broodiness. In addition to prolactin (PRL), follicle-stimulating hormone (FSH) and luteinizing hormone (LH), genes including oxytocin-neurophysin (OXT), chordin-like protein 1 (CHRDL1) and growth hormone (GH), expressed in the pituitary gland, are new candidate molecules that may be involved in broodiness in geese. Heme oxygenase 1 (HMOX1) in the pituitary gland, the proto-oncogene c-Fos (FOS), heat shock protein 90-alpha (HSP90AA), and cyclin-dependent kinase 1 (CDK1) in the ovary that may consolidate and transduce signals regulating the HPGA during broodiness in geese.

miR-26a-5p Regulates TNRC6A Expression and Facilitates Theca Cell Proliferation in Chicken Ovarian Follicles

Ovarian theca cells play an indispensable role in ovarian follicular development and hormone secretion. miR-26a-5p was reported to be differentially expressed in mature and immature chicken ovaries in our previous study; however, the role of miR-26a-5p in regulating ovarian follicle function is still unclear. In this study, we demonstrated that the expression dynamics of TNRC6A mRNA in either chicken ovaries or follicles showed an opposite trend compared with that of chicken miR-26a-5p expression. miR-26a-5p inhibited TNRC6A mRNA expression by directly targeting its 3'-untranslated region in cultured chicken theca cells. Overexpression of miR-26a-5p promoted chicken follicular theca cell proliferation in vitro. Furthermore, overexpression of miR-26a-5p and knockdown of TNRC6A significantly upregulated the antiapoptotic BCL-2 gene. Taken together, this study revealed the expression dynamics of miR-26a-5p and TNRC6A in chicken ovaries and ovarian follicles and the relationship between the expression of miR-26a-5p and TNRC6A in chicken ovarian theca cells. These results suggest that miR-26a-5p facilitates chicken ovarian theca cell proliferation by targeting the TNRC6A gene.

Identification and profile of microRNAs in Xiang pig testes in four different ages detected by Solexa sequencing

To further understand the role of microRNA (miRNA) during testicular development, we constructed four small RNA libraries from the testes of the Chinese indigenous Xiang pig at four different ages, which were sequenced using high-throughput Solexa deep sequencing methods. It yielded over 23 million high-quality reads and 1,342,579 unique sequences. At two and three months of age, the proportion which represented miRNAs was the most abundant class of small RNAs, but it was gradually replaced by the category that represented piRNAs in adult testes. We identified 543 known and homologous conserved porcine miRNAs and 49 potential novel miRNAs. There were 306 known miRNAs which were co-expressed in four libraries. Six miRNAs and three potential novel miRNAs were validated in testes and sperms of Xiang pig by RT-qPCR method. Many clusters of mature miRNA variants were observed, in which let-7 family was the most abundant one. After comparison among libraries, 204 miRNAs were identified as being differentially expressed and likely involved in the development and spermatogenesis of pig testes. This work presented a general genome-wide expression profile of the testes-expressed small RNAs in different ages of pig testes. Our results suggested that miRNAs performed a role in the regulation of mRNAs in puberty pig testes while piRNAs likely functioned mainly in sexually mature pig testes.

Genome-wide identification of novel ovarian-predominant miRNAs: new insights from the medaka (Oryzias latipes)

MicroRNAs (miRNAs) are small, highly conserved non-coding RNAs that play important roles in the regulation of many physiological processes. However, the role of miRNAs in vertebrate oocyte formation (i.e., oogenesis) remains poorly investigated. To gain new insights into the roles of miRNAs in oogenesis, we searched for ovarian-predominant miRNAs. Using a microarray displaying 3,800 distinct miRNAs originating from different vertebrate species, we identified 66 miRNAs that are expressed predominantly in the ovary. Of the miRNAs exhibiting the highest overabundance in the ovary, 20 were selected for further analysis. Using a combination of QPCR and in silico analyses, we identified 8 novel miRNAs that are predominantly expressed in the ovary, including 2 miRNAs (miR-4785 and miR-6352) that exhibit strict ovarian expression. Of these 8 miRNAs, 7 were previously uncharacterized in fish. The strict ovarian expression of miR-4785 and miR-6352 suggests an important role in oogenesis and/or early development, possibly involving a maternal effect. Together, these results indicate that, similar to protein-coding genes, a significant number of ovarian-predominant miRNA genes are found in fish.

Molecular cloning of the SMAD4 gene and its mRNA expression analysis in ovarian follicles of the Yangzhou goose (Anser cygnoides)

Mothers against decapentaplegic homolog 4 (SMAD4) is an important protein in animal reproduction. It plays pivotal roles in cellular pathways, including apoptosis. The expression profile of the SMAD4 gene in goose ovarian follicles has not been reported. In this study, the SMAD4 coding sequence was cloned from the Yangzhou goose. A phylogenetic analysis was performed and mRNA expression was examined in various tissues using quantitative real-time PCR. An alternative splice form of SMAD4, SMAD4-b having 1656 bp, was identified. SMAD4-a mRNA was widely expressed in various healthy tissues, whereas SMAD4-b was very weakly expressed. SMAD4 mRNA in the ovary and oviduct was significantly higher than that in the pituitary and hypothalamus. SMAD4 mRNA expression analysis in hierarchical follicles showed that the level of SMAD4 mRNA was higher in large white follicles and post-ovulatory follicles than in the other follicles. The results indicate that SMAD4 might be involved in the recruitment of hierarchical follicles.

High-throughput sequencing reveals differential expression of miRNAs in prehierarchal follicles of laying and brooding geese

Broodiness is the primary factor influencing egg production in geese, in which several genes and miRNAs participate. Detailed spatiotemporal profiles of miRNAs encompassing follicle development levels, however, are lacking. In this study, we collected preovulatory follicles (classified as small white follicles, large white follicles, and small yellow follicles) from brooding and laying geese and aimed to analyze microRNA (miRNA or miR) during folliculogenesis. High-throughput sequencing and bioinformatics analysis were used to identify the miRNAs involved in follicle development. The let7 family, miR-10 family, and miR-143 family were abundant in these libraries, and they have been suggested to play a housekeeping role during folliculogenesis. Joint comparisons revealed 23 upregulated and 21 downregulated miRNAs (in at least two comparisons of follicles during brooding and laying, P < 0.1) in the laying stage. Unlike reproduction pathways reported for ovaries, GO and KEGG analysis suggested pathways for cell apoptosis and proliferation, such as the regulation of actin cytoskeleton, endocytosis, axon guidance, pathways in cancer, tight junctions, focal adhesion, the MAPK signaling pathway, cytokine-cytokine receptor interactions, and the Wnt signaling pathway in folliculogenesis. This study revealed the miRNAs that were directly involved in follicular atresia, and our results added to the understanding of the functional involvement of miRNAs during specific stages of follicle development.

The dopamine β-hydroxylase gene in Chinese goose (Anas cygnoides): cloning, characterization, and expression during the reproductive cycle

BACKGROUND

Dopamine β-hydroxylase (DBH) is a critical enzyme in the biosynthesis of catecholamines. This enzyme's role in neuroendocrine regulation is well known, but there are some indications that it may also modulate reproduction and endocrine in mammals and birds. We selected goose (Anas cygnoides) as an ideal model species for investigating the role of DBH in avian reproduction.

RESULTS

Full-length cDNA encoding DBH was cloned from Zhedong goose using reverse transcription PCR and rapid amplification of cDNA ends. The cDNA consisted of a 126-base pair (bp) 5'-untranslated region (UTR), a 379-bp 3'-UTR, and an 1896-bp open reading frame encoding a polypeptide of 631 amino acids. The deduced amino acid sequence of gDBH shared high homology with an analogue from other birds and contained three conserved domains from a mono-oxygenase family including a DOMON domain and two Cu2_mono-oxygen domains. Real-time quantitative PCR analysis showed that gDBH mRNA was expressed in both reproductive and endocrine tissues of Zhedong goose, specifically in the hypothalamus, pituitary, ovary, and oviduct. More DBH mRNA of reproductive and endocrine tissues was detected at ovulation than at oviposition in Zhedong goose. Evidence of opposite trend of gDBH expression was found between the hypothalamus-pituitary and oviduct during the ovulation phase and the broody phase. In addition, we assessed DBH mRNA expression during ovulation in two breeds of geese that differ in egg production. The reproductive and endocrine tissues of Yangzhou geese with higher egg production had more gDBH expression than Zhedong geese. Finally, the five non-synonymous SNP(c.1739 C > T, c.1760G > T, c.1765A > G, c.1792 T > C and c.1861G > C) were identified in the coding region of DBH gene between Zhedong goose and Yangzhou goose.

CONCLUSIONS

We conclude that goose DBH mRNA show obvious periodically variation in reproductive and endocrine tissues during the reproductive cycle in geese.

Goose broodiness is involved in granulosa cell autophagy and homeostatic imbalance of follicular hormones

Broodiness is observed in most domestic fowls and influences egg production. The goose is one of the most important waterfowls, having strong broody behavior. However, whether autophagy and follicular internal environment play a role in the broodiness behavior of goose is unknown. In this report, we analyzed the follicular internal environment and granulosa cell autophagy of goose follicles. The results show that the contents of hormones, including prolactin (PRL), progesterone (P4), and estradiol (E2), increased in broody goose follicles. Most importantly, the level of granulosa cell autophagy in broody goose follicles was elevated, detected by electron microscopy and western blotting. Also, the expressions of positive regulators of autophagy, including miR-7, miR-29, miR-100, miR-181, PRLR, LC3, p53,Beclin1, Atg9, and Atg12, were up-regulated and the expressions of negative regulators of autophagy, including miR-34b and miR-34c, were down-regulated in broody goose follicles. Our results suggest that goose broodiness is involved in increased granulosa cell autophagy and homeostasis imbalance of internal environment in the follicles. This work contributes to our knowledge of goose broodiness and may influence egg production.

Molecular cloning and expression patterns of the cholesterol side chain cleavage enzyme (CYP11A1) gene during the reproductive cycle in goose (Anas cygnoides)

Background

CYP11A1, a gene belonging to the family 11 of cytochrome P450, encodes a crucial steroidogenic enzyme that catalyzes the initial step in the production of all classes of steroids. Many studies show that CYP11A1 plays a role in ovary function. However, the role of CYP11A1 in goose reproductive cycle remains largely unknown.

Results

In this study, full-length CYP11A1 cDNA of Zhedong goose was obtained using reverse transcription polymerase chain reaction (RT-PCR) and rapid amplification of cDNA ends (RACE). The cDNA consisted of a 96-base pair (bp) 5′untranslated region (UTR), a 179-bp 3′UTR and a 1509-bp open reading frame. The open reading frame encodes a putative 503 amino acid protein that shares high homology with CYP11A1 of other birds. The amino acid sequence possesses conserved domains of the P450 superfamily, which include the steroid-binding domain and the heme-binding region. Real-time quantitative polymerase chain reaction (qPCR) analysis revealed CYP11A1 mRNA was expressed ubiquitously in every Zhedong goose tissue analyzed, including the heart, liver, glandular stomach, lung, spleen, kidney, intestinum tenue, intestinum crassum, cerebrum, cerebellum, muscle, oviduct, pituitary, hypothalamus and ovary.. The relatively low levels of CYP11A1 mRNA were detected in pituitary, ovary and oviduct tissues at ovulation when compared with levels at oviposition. Interestingly, higher expression was observed in ovary and oviduct tissues during brooding. Lastly, higher mRNA expression of Yangzhou geese was detected during the ovulation period than that of Zhedong geese.

Conclusions

Our findings reveal the sequence characterization and expression patterns of the CYP11A1 gene during the goose reproductive cycle, which may provides correlative evidence that CYP11A1 expression is important in reproduction activity.

Coordinated microRNA and messenger RNA expression profiles for understanding sexual dimorphism of gonads and the potential roles of microRNA in the steroidogenesis pathway in Nile tilapia (Oreochromis niloticus)

Sexual dimorphism is a widespread phenomenon in animals. However, the potential role of microRNAs (miRNAs) in regulating this dimorphism is not fully understood. In our study, we used an integrated approach to identify functional targets of miRNA by combining the paired expression profiles of miRNAs and messenger RNAs (mRNAs) in ovaries and testes of young Nile tilapia, Oreochromis niloticus. The results revealed that 67 upregulated and nine downregulated miRNAs and 2299 upregulated and 3260 downregulated genes were identified in the ovary compared with those in the testis (P < 0.01). The target genes of differentially expressed miRNAs were predicted and overlapped with the differentially expressed mRNAs. Furthermore, Kyoto Encyclopedia of Genes and Genomes pathway analyses were conducted in these coincident genes. By correlating miRNA-mRNA and predicting computational target, two types of negatively regulatory miRNA-mRNA correlations (upregulated or downregulated miRNA and downregulated or upregulated mRNA) were obtained. Seven functional miRNA-target gene pairs, miR-17-5p/DMRT1, miR-20a/DMRT1, miR-138/CYP17A2, miR-338/CYP17A2, miR-200a/CYP17A2, miR-456/AMH, and miR-138/AMH, were predicted at the sequence level and further detected by real-time polymerase chain reaction on the basis of the significantly negative relationships. Our results suggest that the integrated analysis of miRNA and mRNA expression profiling can provide novel insights into the molecular mechanism of sexual dimorphism.

Identification of MicroRNAs in Zebrafish Spermatozoa

MicroRNAs (miRNAs) participate in almost all biological processes. Plenty of evidences show that some testis- or spermatozoa-specific miRNAs play crucial roles in the process of gonad and germ cell development. In this study, the spermatozoa miRNA profiles were investigated through a combination of illumina deep sequencing and bioinformatics analysis in zebrafish. Deep sequencing of small RNAs yielded 11,820,680 clean reads. By mapping to the zebrafish genome, we identified 400 novel and 204 known miRNAs that could be grouped into 104 families. Furthermore, we selected the six highest expressions of known miRNAs to detect their expression patterns in different tissues by stem-loop quantitative real-time polymerase chain reaction. We found that among the six miRNAs, dre-miR-202-5p displayed specific and high expression in zebrafish spermatozoa and testis. Fluorescence in situ hybridization analysis indicated that dre-miR-202-5p was predominantly expressed in all kind of germ cells at different spermatogenetic stages, including spermatogonia and spermatozoa, but barely expressed in the germ cells in the ovary. This sex-biased expression pattern suggests that dre-miR-202-5p might be related to spermatogenesis and the functioning of spermatozoa. The identification of miRNAs in zebrafish spermatozoa and germ cells offers new insights into the spermatogenesis and spermatozoa in the teleost and other vertebrates.

MicroRNA Expression Profile in Penile Cancer Revealed by Next-Generation Small RNA Sequencing

Penile cancer (PeCa) is a relatively rare tumor entity but possesses higher morbidity and mortality rates especially in developing countries. To date, the concrete pathogenic signaling pathways and core machineries involved in tumorigenesis and progression of PeCa remain to be elucidated. Several studies suggested miRNAs, which modulate gene expression at posttranscriptional level, were frequently mis-regulated and aberrantly expressed in human cancers. However, the miRNA profile in human PeCa has not been reported before. In this present study, the miRNA profile was obtained from 10 fresh penile cancerous tissues and matched adjacent non-cancerous tissues via next-generation sequencing. As a result, a total of 751 and 806 annotated miRNAs were identified in normal and cancerous penile tissues, respectively. Among which, 56 miRNAs with significantly different expression levels between paired tissues were identified. Subsequently, several annotated miRNAs were selected randomly and validated using quantitative real-time PCR. Compared with the previous publications regarding to the altered miRNAs expression in various cancers and especially genitourinary (prostate, bladder, kidney, testis) cancers, the most majority of deregulated miRNAs showed the similar expression pattern in penile cancer. Moreover, the bioinformatics analyses suggested that the putative target genes of differentially expressed miRNAs between cancerous and matched normal penile tissues were tightly associated with cell junction, proliferation, growth as well as genomic instability and so on, by modulating Wnt, MAPK, p53, PI3K-Akt, Notch and TGF-β signaling pathways, which were all well-established to participate in cancer initiation and progression. Our work presents a global view of the differentially expressed miRNAs and potentially regulatory networks of their target genes for clarifying the pathogenic transformation of normal penis to PeCa, which research resource also provides new insights into future investigations aimed to explore the in-depth mechanisms of miRNAs and other small RNAs including piRNAs in penile carcinogenesis regulation and effective target-specific theragnosis.