Identification of Differentially Expressed Genes in Porcine Ovaries at Proestrus and Estrus Stages Using RNA-Seq Technique

Application of RNA-Seq Technology for Screening Reproduction-Related Differentially Expressed Genes in Tibetan and Yorkshire Pig Ovarian Tissue

The purpose of this study was to explore and verify genes that regulate the reproductive traits of Tibetan pigs at the mRNA level. The ovarian tissues of Tibetan pigs (TPs) and Yorkshire pigs (YPs) were selected as research objects, and cDNA libraries of the ovarian tissue transcripts of Tibetan pigs and Yorkshire pigs were successfully constructed by the RNA-Seq technique. A total of 651 differentially expressed genes (DEGs) were screened, including 414 up-regulated genes and 237 down-regulated genes. Through GO and KEGG enrichment analysis, it was found that these differentially expressed genes were significantly enriched in cell process, reproductive process, reproduction, cell proliferation, binding, and catalytic activity, as well as oxidative phosphorylation, endocrine resistance, thyroid hormone, Notch, and other signal transduction pathways. Genes significantly enriched in pathways closely related to reproductive regulation were analyzed and selected, and the AR, CYP11A1, CYP17A1, INHBA, ARRB2, EGFR, ETS1, HSD17B1, IGF1R, MIF, SCARB1, and SMAD4 genes were identified as important candidate genes. Twelve differentially expressed genes related to reproduction were verified by RT-qPCR. The results showed that the expression of the AR, CYP17A1, EGFR, ETS1, IGF1R, and SMAD4 genes was significantly higher in Tibetan pigs than in Yorkshire pigs, while the expression of the CYP11A1, INHBA, ARRB2, HSD17B, MIF, and SCARB1 genes in Tibetan pigs was significantly lower than in Yorkshire pigs. The purpose of this study is to provide a theoretical basis for exploring the molecular mechanism of reproductive trait effect genes and the application of molecular breeding in Tibetan pigs.

Deciphering Estrus Expression in Gilts: The Role of Alternative Polyadenylation and LincRNAs in Reproductive Transcriptomics

The fertility rate and litter size of female pigs are critically affected by the expression of estrus. The objective of this study was to elucidate the regulatory mechanisms of estrus expression by analyzing the differential expression of genes and long intergenic non-coding RNAs (lincRNA), as well as the utilization of alternative polyadenylation (APA) sites, in the vulva and vagina during the estrus and diestrus stages of Large White and indigenous Chinese Mi gilts. Our study revealed that the number of differentially expressed genes (DEG) in the vulva was less than that in the vagina, and the DEGs in the vulva were enriched in pathways such as "neural" pathways and steroid hormone responses, including the "Calcium signaling pathway" and "Oxytocin signaling pathway". The DEGs in the vagina were enriched in the "Metabolic pathways" and "VEGF signaling pathway". Furthermore, 27 and 21 differentially expressed lincRNAs (DEL), whose target genes were enriched in the "Endocrine resistance" pathway, were identified in the vulva and vagina, respectively. Additionally, we observed that 63 and 618 transcripts of the 3'-untranslated region (3'-UTR) were lengthened during estrus in the vulva and vagina, respectively. Interestingly, the genes undergoing APA events in the vulva exhibited species-specific enrichment in neural or steroid-related pathways, whereas those in the vagina were enriched in apoptosis or autophagy-related pathways. Further bioinformatic analysis of these lengthened 3'-UTRs revealed the presence of multiple miRNAs binding sites and cytoplasmic polyadenylation element (CPE) regulatory aspects. In particular, we identified more than 10 CPEs in the validated lengthened 3'-UTRs of the NFIX, PCNX4, CEP162 and ABHD2 genes using RT-qPCR. These findings demonstrated the involvement of APA and lincRNAs in the regulation of estrus expression in female pigs, providing new insights into the molecular mechanisms underlying estrus expression in pigs.

Transcriptomic and metabolomic analyses of the ovaries of Taihe black-bone silky fowls at the peak egg-laying and nesting period

The poor reproductive performance of most local Chinese chickens limits the economic benefits and output of related enterprises. As an excellent local breed in China, Taihe black-bone silky fowl is in urgent need of our development and utilization. In this study, we performed transcriptomic and metabolomic analyses of the ovaries of Taihe black-bone silky fowls at the peak egg-laying period (PP) and nesting period (NP) to reveal the molecular mechanisms affecting reproductive performance. In the transcriptome, we identified five key differentially expressed genes (DEGs) that may affect the reproductive performance of Taihe black-bone silky fowl: BCHE, CCL5, SMOC1, CYTL1, and SCIN, as well as three important pathways: the extracellular region, Neuroactive ligand-receptor interaction and Cytokine-cytokine receptor interaction. In the metabolome, we predicted three important ovarian significantly differential metabolites (SDMs): LPC 20:4, Bisphenol A, and Cortisol. By integration analysis of transcriptome and metabolome, we identified three important metabolite-gene pairs: "LPC 20:4-BCHE", "Bisphenol A-SMOC1", and "Cortisol- SCIN". In summary, this study contributes to a deeper understanding of the regulatory mechanism of egg production in Taihe black-bone silky fowl and provides a scientific basis for improving the reproductive performance of Chinese local chickens.

Expression of genes regulating cell division in porcine follicular granulosa cells

BACKGROUND

Cell cycle regulation influences the proliferation of granulosa cells and affects many processes related to ovarian folliclular growth and ovulation. Abnormal regulation of the cell cycle can lead to many diseases within the ovary. The aim of this study was to describe the expression profile of genes within granulosa cells, which are related to the formation of the cytoskeleton, organization of cell organelles inside the cell, and regulation of cell division. Established in vitro primary cultures from porcine ovarian follicle granulosa cells were maintained for 48, 96, 144 h and evaluated via microarray expression analysis.

RESULTS

Analyzed genes were assigned to 12 gene ontology groups "actin cytoskeleton organization", "actin filament organization", "actin filament-based process", "cell-matrix adhesion", "cell-substrate adhesion", "chromosome segregation", "chromosome separation", "cytoskeleton organization", "DNA integrity checkpoint", "DNA replication initiation", "organelle fision", "organelle organization". Among the genes with significantly changed expression, those whose role in processes within the ovary are selected for consideration. Genes with increased expression include (ITGA11, CNN1, CCl2, TPM2, ACTN1, VCAM-1, COL3A1, GSN, FRMD6, PLK2). Genes with reduced expression inlcude (KIF14, TACC3, ESPL1, CDC45, TTK, CDC20, CDK1, FBXO5, NEK2-NIMA, CCNE2). For the results obtained by microarray expressions, quantitative validation by RT-qPCR was performed.

CONCLUSIONS

The results indicated expression profile of genes, which can be considered as new molecular markers of cellular processes involved in signaling, cell structure organization. The expression profile of selected genes brings new insight into regulation of physiological processes in porcine follicular granulosa cells during primary in vitro culture.

Analysis of serum reproductive hormones and ovarian genes in pubertal female goats

BACKGROUND

Age at puberty is an important factor affecting goat fertility, with endocrine and genetic factors playing a crucial role in the onset of puberty. To better understand the relationship between endocrine and genetic factors and mechanisms underlying puberty onset in goats, reproductive hormone levels were analyzed by ELISA and ultraperformance liquid chromatography-multiple reaction monitoring-multistage/mass spectrometry and RNA sequencing was performed to analyze ovarian genes.

RESULTS

Serum follicle stimulating hormone, luteinizing hormone, estradiol, 11-deoxycortisol, 11-deoxycorticosterone, corticosterone, cortisone, and cortisol levels were found to be higher but progesterone were lower in pubertal goats as compared to those in prepubertal goats (P < 0.05). A total of 18,139 genes were identified in cDNA libraries, and 75 differentially expressed genes (DEGs) were identified (|log₂ fold change|≥ 1, P ≤ 0.05), of which 32 were significantly up- and 43 were down-regulated in pubertal goats. Gene ontology enrichment analyses indicated that DEGs were mainly involved in "metabolic process," "signaling," "reproduction," and "growth." Further, DEGs were significantly enriched in 91 Kyoto Encyclopedia of Genes and Genomes pathways, including estrogen signaling pathway, steroid hormone biosynthesis, and cAMP signaling pathway. Bioinformatics analysis showed that PRLR and THBS1 were highly expressed in pubertal ovaries, and ZP3, ZP4, and ASTL showed low expression, suggesting their involvement in follicular development and lutealization.

CONCLUSIONS

To summarize, serum hormone changes and ovarian DEGs expression were investigated in our study. Further studies are warranted to comprehensively explore the functions of DEGs in goat puberty.

Comprehensive Transcriptome Analysis of Follicles from Two Stages of the Estrus Cycle of Two Breeds Reveals the Roles of Long Intergenic Non-Coding RNAs in Gilts

Simple Summary

This study provides new perspectives about the roles of lincRNAs in the estrus expression of gilts, which is correlated with ovarian steroid hormone and follicular development. Follicular tissues from two stages of the estrus cycle of Large White and Mi gilts were used for RNA-seq. Some genes and lincRNAs related to estrus expression in pigs were discovered. PPI and ceRNA networks related to the estrus expression were constructed. These results suggest that the estrus expression may be affected by lincRNAs and their target genes.

Abstract

Visible and long-lasting estrus expression of gilts and sows effectively sends a mating signal. To reveal the roles of Long Intergenic Non-coding RNAs (lincRNAs) in estrus expression, RNA-seq was used to investigate the lincRNAs expression of follicular tissues from Large White gilts at diestrus (LD) and estrus (LE), and Chinese Mi gilts at diestrus (MD) and estrus (ME). Seventy-three differentially expressed lincRNAs (DELs) were found in all comparisons (LE vs. ME, LD vs. LE, and MD vs. ME comparisons). Eleven lincRNAs were differentially expressed in both LD vs. LE and MD vs. ME comparisons. Fifteen DELs were mapped onto the pig corpus luteum number Quantitative Trait Loci (QTL) fragments. A protein–protein interaction (PPI) network that involved estrus expression using 20 DEGs was then constructed. Interestingly, three predicted target DEGs (PTGs) (CYP19A1 of MSTRG.10910, CDK1 of MSTRG.10910 and MSTRG.23984, SCARB1 of MSTRG.1559) were observed in the PPI network. A competitive endogenous RNA (ceRNA) network including three lincRNAs, five miRNAs, and five genes was constructed. Our study provides new insight into the lincRNAs associated with estrus expression and follicular development in gilts.

Time Series Ovarian Transcriptome Analyses of the Porcine Estrous Cycle Reveals Gene Expression Changes during Steroid Metabolism and Corpus Luteum Development

Simple Summary

The estrous cycle, which is divided into follicular and luteal phases based on ovulation, is influenced by reproductive hormones which affect reproduction and cause changes in the reproductive system of the pig. As the main reproductive organ, the ovary is involved in ovulation and changes in the corpus luteum. We aimed to identify dynamic changes in gene expression through differentially expressed gene profiling and to provide a comprehensive understanding of the molecular mechanisms that occur in the pig ovary during the estrous cycle. The transcriptome analysis revealed that the dynamic change in gene expression was more activated in the luteal phase than in the follicular phase. Functional analysis revealed that the metestrus and diestrus periods are important in preparation for pregnancy or the next estrous cycle after ovulation.

Abstract

The porcine estrous cycle is influenced by reproductive hormones, which affect porcine reproduction and result in physiological changes in the reproductive organs. The ovary is involved in ovulation, luteinization, corpus luteum development, and luteolysis. Here, we aimed to provide a comprehensive understanding of the gene expression patterns in porcine ovarian transcriptomes during the estrous cycle through differentially expressed genes profiling and description of molecular mechanisms. The transcriptomes of porcine ovary were obtained during the estrous cycle at three-day intervals from day 0 to day 18 using RNA-seq. At seven time points of the estrous cycle, 4414 DEG were identified; these were classified into three clusters according to their expression patterns. During the late metestrus and diestrus periods, the expression in cluster 1 increased rapidly, and steroid biosynthesis was significant in the pathway. Cluster 2 gene expression patterns represented the cytokine–cytokine receptor interaction in significant pathways. In cluster 3, the hedgehog signaling pathway was selected as the significant pathway. Our study exhibited dynamic gene expression changes with these three different patterns of cluster 1, 2, and 3. The results helped identify the functions and related significant genes especially during the late metestrus and diestrus periods in the estrous cycle.

Molecular assessment and transcriptome profiling of wild fish populations of Oryzias mekongensis and O. songkhramensis (Adrianichthyidae: Beloniformes) from Thailand

Among the fish of the genus Oryzias, two species are frequently used as model animals in biological research. In Thailand, Oryzias mekongensis is usually found in natural freshwater near the Mekong Basin in the northeast region, while O. songkhramensis inhabits the Songkhram Basin. For differential morphological identification, the coloured bands on the dorsal and ventral margins of the caudal fin are used to distinguish O. mekongensis from O. songkhramensis. However, these characteristics are insufficient to justify species differentiation, and little molecular evidence is available to supplement them. This study aimed to investigate the molecular population and transcriptome profiles of adult O. mekongensis and O. songkhramensis. In the molecular tree based on cytochrome b sequences, O. mekongensis exhibited four clades that were clearly distinguished from O. songkhramensis. Clade 1 of the O. mekongensis population was close to the Mekong River and lived in the eastern portion of the upper northeast region. Clade 2 was far from the Mekong River and inhabited the middle region of the Songkhram River. Clade 3 was positioned to the west of the Songkhram River, and clade 4 was to the south of the Songkhram River Basin. After RNA sequencing using an Illumina HiSeq 2500 platform, the gene category annotations hardly differentiated the species and were discussed in the text. Based on the present findings, population dispersal of these Oryzias species might be associated with geographic variations of the upper northeast region. Molecular genetics and transcriptome profiling might advance our understanding of the evolution of teleost fish.

Identification of Ovarian Circular RNAs and Differential Expression Analysis between MeiShan and Large White Pigs

MeiShan and Large White pigs differ in their female fecundity. However, the mechanisms behind the gene expression and regulation that cause these differences remain unclear. In this study, we profiled circRNAs and identified 5,879 circRNAs from the ovaries of MeiShan and Large White pigs. Eighty-five circRNAs were differentially expressed between the two pig breeds. Of these, 37 were up-regulated and 48 were down-regulated in MeiShan pigs. Gene ontology enrichment analysis suggested that the differentially expressed circRNA were involved in the hormone-mediated signaling pathway. We verified that circSCIN and its parent gene, scinderin (SCIN), were differentially expressed by reverse transcription and quantitative PCR (RT-qPCR). Luciferase assays demonstrated that circSCIN can target and sponge miR-133 and miR-148a/b. The identification of differentially expressed circRNAs (DECs) and their regulatory functions increased our understanding of the differences in reproductive efficiency between MeiShan and Large White pigs.

Serial gene co-expression network approach to mine biological meanings from integrated transcriptomes of the porcine endometrium during estrous cycle

The estrous cycle is a complex process regulated by several hormones. To understand the dynamic changes in gene expression that takes place in the swine endometrium during the estrous cycle relative to the day of estrus onset, we performed RNA-sequencing analysis on days 0, 3, 6, 9, 12, 15, and 18, resulting in the identification of 4495 differentially expressed genes (DEGs; Q ≤ 0.05 and |log₂FC| ≥ 1) at various phases in the estrous cycle. These DEGs were integrated into multiple gene co-expression networks based on different fold changes and correlation coefficient (R²) thresholds and a suitable network, which included 899 genes (|log2FC| ≥ 2 and R² ≥ 0.99), was identified for downstream analyses based on the biological relevance of the Gene Ontology (GO) terms enriched. The genes in this network were partitioned into 6 clusters based on the expression pattern. Several GO terms including cell cycle, apoptosis, hormone signaling, and lipid biosynthetic process were found to be enriched. Furthermore, we found 15 significant KEGG pathways, including cell adhesion molecules, cytokine-cytokine receptor signaling, steroid biosynthesis, and estrogen signaling pathways. We identified several genes and GO terms to be stage-specific. Moreover, the identified genes and pathways extend our understanding of porcine endometrial regulation during estrous cycle and will serve as a good resource for future studies.