Expression characteristics of pineal miRNAs at ovine different reproductive stages and the identification of miRNAs targeting the AANAT gene.

BACKGROUND: Many recent studies have shown that miRNAs play important roles in the regulation of animal reproduction, including seasonal reproduction. The pineal gland is a crucial hub in the regulation of seasonal reproduction. However, little is known about the expression characteristics of pineal miRNAs in different reproductive seasons (anestrus and breeding season). Therefore, the expression profiles and regulatory roles of ovine pineal miRNAs were investigated during different reproductive stages using Solexa sequencing technology and dual luciferase reporter assays. RESULTS: A total of 427 miRNAs were identified in the sheep pineal gland. Significant differences in miRNA expression were demonstrated between anestrus and the breeding season in terms of the frequency distributions of miRNA lengths, number of expressed miRNAs, and specifically and highly expressed miRNAs in each reproductive stage. KEGG analysis of the differentially expressed (DE) miRNAs between anestrus and the breeding season indicated that they are significantly enriched in pathways related to protein synthesis, secretion and uptake. Furthermore, transcriptome analysis revealed that many target genes of DE miRNAs in the ribosome pathway showed relatively low expression in the breeding season. On the other hand, analyses combining miRNA-gene expression data with target relationship validation in vitro implied that miR-89 may participate in the negative regulation of aralkylamine N-acetyltransferase (AANAT) mRNA expression by targeting its 3'UTR at a unique binding site. CONCLUSIONS: Our results provide new insights into the expression characteristics of sheep pineal miRNAs at different reproductive stages and into the negative regulatory effects of pineal miRNAs on AANAT mRNA expression.

Expression Profiling and Functional Characterization of miR-26a and miR-130a in Regulating Zhongwei Goat Hair Development via the TGF-ß/SMAD Pathway.

The Zhongwei goat is an important and unique goat breed indigenous to China. It has a natural hair curling phenotype at birth, but the degree of curling gradually decreases with growth. The molecular mechanism underlying the dynamic changes in the wool curvature in Zhongwei goats is poorly understood. MicroRNAs (miRNAs) play important roles in many biological processes, including hair growth and development. In this study, we selected skins from Zhongwei goats at different ages (45 and 108 days) that exhibited different levels of hair curvature and performed miRNA sequencing to explore the molecular mechanism of hair bending. In total, 28 significantly differentially expressed miRNAs (DE miRNAs) were identified in the three groups of samples between the two developmental stages. An analysis of the target genes of the above-mentioned DE miRNAs by the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses indicated that the DE miRNAs were involved in signal pathways which were previously associated with hair bending and hair follicle development, such as the TGF-ß/SMAD, PI3K-Akt, JAK-STAT, and MAPK pathways. A comprehensive analysis of the correlations between the miRNA-seq results and issued transcriptional findings indicated that SMAD1 was a target gene of miR-26a and SMAD5 was a target gene of miR-130a. Furthermore, goat dermal papilla cells were successfully isolated and purified to determine the role of miRNAs in follicle development in vitro. The study results demonstrated that miR-130a and miR-26a had significant effects on the proliferation of dermal papilla cells. In addition, the detection results of mRNA and protein levels indicate that the overexpression of miR-26a can promote the expression of related genes in the TGF-ß/SMAD pathway, while miR-130a has the opposite substitution effect. The dual luciferase report test showed that miR-26a targeted the SMAD1 gene and reduced the expression of the SMAD1 protein in hair papillary cells. Our results identified DE microRNAs which perhaps change at the time of hair straightening in Zhongwei goats and explore the role of miR-26a and miR-130a in dermal papilla cells proliferation. The present study provided a theoretical basis to explore the mechanisms underlying the Zhongwei hair growth and curly phenotype.

Epigenetic Regulation of miR-25 and Lnc107153 on Expression of Seasonal Estrus Key Gene CHGA in Sheep.

Pituitary pars tuberalis (PT) plays an important role as the transmission center in the seasonal reproduction of animals. It helps convert external photoperiod signals into intrinsic seasonal reproduction signals. In sheep PT, specific expression patterns of several genes (including short photoperiod-induced gene CHGA and long photoperiod genes EYA3 and TSHß) under different photoperiods are crucial characteristics during this signal transduction. Recent studies have revealed the role of epigenetics in regulating the expression of seasonal reproductive key genes. Therefore, we explored whether microRNAs and LncRNAs regulated the expressions of the above key genes. Firstly, the expression of miR-25 and CHGA showed a significant negative correlation in sheep PT. Results of the dual luciferase reporter assay and miR-25 overexpression indicated that miR-25 could inhibit the expression of CHGA by specifically binding to its 3'UTR region in pituitary cells. Then, expression negative correlation and dual luciferase reporter analyses were used to screen and identify the candidate LncRNA (Lnc107153) targeted by miR-25. Finally, the results of fluorescence in situ hybridization and Lnc107153 overexpression suggested that Lnc107153 and miR-25 were involved in the epigenetic regulation of CHGA expression. However, the expressions of EYA3 and TSHß were not regulated by miRNAs. These results will provide new insights into the epigenetic regulatory network of key genes in sheep seasonal reproduction.

Expression characteristics of piRNAs in ovine luteal phase and follicular phase ovaries.

PIWI-interacting RNAs (piRNAs), as a novel class of small non-coding RNAs that have been shown to be indispensable in germline integrity and stem cell development. However, the expressed characteristics and regulatory roles of piRNAs during different reproductive phases of animals remain unknown. In this study, we investigated the piRNAs expression profiles in ovaries of sheep during the luteal phase (LP) and follicular phase (FP) using the Solexa sequencing technique. A total of 85,219 and 1,27,156 piRNAs tags were identified in ovine ovaries across the two phases. Most expressed piRNAs start with uracil. piRNAs with a length of 24 nt or 27-29 nts accounted for the largest proportion. The obvious ping-pong signature appeared in the FP ovary. The piRNA clusters in the sheep ovary were unevenly distributed on the chromosomes, with high density on Chr 3 and 1. For genome distribution, piRNAs in sheep ovary were mainly derived from intron, CDS, and repeat sequence regions. Compared to the LP ovary, a greater number of expressed piRNA clusters were detected in the FP ovary. Simultaneously, we identified 271 differentially expressed (DE) piRNAs between LP and FP ovaries, with 96 piRNAs upregulated and 175 piRNAs downregulated, respectively. Functional enrichment analysis (GO and KEGG) indicated that their target genes were enriched in reproduction-related pathways including oocyte meiosis, PI3K-Akt, Wnt, and TGF-ß signaling pathways. Together, our results highlighted the sequence and expression characteristics of the piRNAs in the sheep ovary, which will help us understand the roles of piRNAs in the ovine estrus cycle.

Characteristics of piRNAs and their comparative profiling in testes of sheep with different fertility.

As a novel class of small RNAs, piRNAs are highly expressed in the animal gonads and their main known role is to inhibit transposon activity for ensuring the correctness and integrity of genome. In order to explore the characteristics of piRNAs in sheep testis and their possible regulatory roles on male reproduction, deep sequencing technology was used to sequence small RNAs and identify piRNAs in testes of sheep. The length of piRNAs in sheep testes showed a unimodal distribution between 26 and 31 nt, with a peak at 29 nt. These piRNAs exhibited obvious ping-pong signature and strand specificity. In the genome, they were mainly aligned to CDS, intron, repetitive sequence regions and unannotated regions. Furthermore, in transposon analysis, piRNAs were aligned predominantly to LINE, SINE, and LTR types of retrotransposon in sheep testes, and the piRNAs derived from each type showed obvious ping-pong signature. The piRNA clusters identified in sheep testes were mainly distributed on chromosomes 3, 7, 15, 17, 18 and 20. The results combining semen determination with pathway enrichment analysis implied that differentially expressed piRNAs between the testes of rams with different fertility might participate in spermatogenesis by regulating multiple pathways closely related to stabilization of blood-testis barrier and renewal and differentiation of spermatogonial stem cell. Taken together, the study provided new insights into the characteristics, origin and expression patterns of piRNAs in sheep testes tissue, which would help us better understand the role of piRNAs in sheep reproduction.

Detection of Novel Variations Related to Litter Size in BMP15 Gene of Luzhong Mutton Sheep (Ovis aries).

Litter size is an important economic trait in the mutton sheep industry. BMP15 is one of the key candidate genes for litter size in sheep. In this study, the entire ORF region of BMP15 was sequenced in 154 Luzhong mutton ewes, and the novel variations were determined. The association between polymorphism in BMP15 and litter size was analyzed using a general linear model. Six out of a total of thirteen variations were identified to be novel. Association analysis indicated that four (SNPs ENSOART00000010201.1:c.352+342C>A, c.352+1232T>C, c.352+1165A>G and c.353-2036T>A) were significantly associated with litter size. The joint analysis among three major genes (BMP15, BMPR1B and GDF9) exhibited significant interaction effects in three combinations (FecB and c.352+1232T>C of BMP15; FecB and c.352+1165A>G of BMP15; c.352+342C>A of BMP15 and ENSOART00000014382.1:c.994G>A of GDF9). For the SNPs c.352+1232T>C and c.352+342C>A, the global distribution of allele frequencies showed that the highest variation frequency occurs in Western Europe. In conclusion, the results demonstrated that BMP15 is a major gene for litter size in Luzhong mutton sheep and candidate SNPs associated with litter size were identified.