Comparative analysis of differentially expressed genes between the ovaries from pregnant and nonpregnant goats using RNA-Seq

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.

Construction of a circRNA- lincRNA-lncRNA-miRNA-mRNA ceRNA regulatory network identifies genes and pathways linked to goat fertility

Background: There is growing interest in the genetic improvement of fertility traits in female goats. With high-throughput genotyping, single-cell RNA sequencing (scRNA-seq) is a powerful tool for measuring gene expression profiles. The primary objective was to investigate comparative transcriptome profiling of granulosa cells (GCs) of high- and low-fertility goats, using scRNA-seq. Methods: Thirty samples from Ji'ning Gray goats (n = 15 for high fertility and n = 15 for low fertility) were retrieved from publicly available scRNA-seq data. Functional enrichment analysis and a literature mining approach were applied to explore modules and hub genes related to fertility. Then, interactions between types of RNAs identified were predicted, and the ceRNA regulatory network was constructed by integrating these interactions with other gene regulatory networks (GRNs). Results and discussion: Comparative transcriptomics-related analyses identified 150 differentially expressed genes (DEGs) between high- and low-fertility groups, based on the fold change (≥5 and ≤-5) and false discovery rate (FDR <0.05). Among these genes, 80 were upregulated and 70 were downregulated. In addition, 81 mRNAs, 58 circRNAs, 8 lincRNAs, 19 lncRNAs, and 55 miRNAs were identified by literature mining. Furthermore, we identified 18 hub genes (SMAD1, SMAD2, SMAD3, SMAD4, TIMP1, ERBB2, BMP15, TGFB1, MAPK3, CTNNB1, BMPR2, AMHR2, TGFBR2, BMP4, ESR1, BMPR1B, AR, and TGFB2) involved in goat fertility. Identified biological networks and modules were mainly associated with ovary signature pathways. In addition, KEGG enrichment analysis identified regulating pluripotency of stem cells, cytokine-cytokine receptor interactions, ovarian steroidogenesis, oocyte meiosis, progesterone-mediated oocyte maturation, parathyroid and growth hormone synthesis, cortisol synthesis and secretion, and signaling pathways for prolactin, TGF-beta, Hippo, MAPK, PI3K-Akt, and FoxO. Functional annotation of identified DEGs implicated important biological pathways. These findings provided insights into the genetic basis of fertility in female goats and are an impetus to elucidate molecular ceRNA regulatory networks and functions of DEGs underlying ovarian follicular development.

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.

The study of candidate genes in the improvement of egg production in ducks – a review

Duck is the second-largest poultry species aside from chicken. The rate of egg production is a major determinant of the economic income of poultry farmers. Among the reproductive organs, the ovary is a major part of the female reproductive system which is highly important for egg production. Based on the importance of this organ, several studies have been carried out to identify candidate genes at the transcriptome level, and also the expression level of these genes at different tissues or egg-laying conditions, and single nucleotide polymorphism (SNPs) of genes associated with egg production in duck. In this review, expression profile and association study analyses at SNPs level of different candidate genes with egg production traits of duck were highlighted. Furthermore, different studies on transcriptome analysis, Quantitative Trait Loci (QTL) mapping, and Genome Wide Association Study (GWAS) approach used to identify potential candidate genes for egg production in ducks were reported. This review would widen our knowledge on molecular markers that are associated or have a positive correlation to improving egg production in ducks, for the increasing world populace.

Granulosa cell transcriptomic study reveals the differential regulation of lncRNAs and mRNAs related to follicle development in goat

Mammalian follicle development is a complex biological process regulated by several factors. More than 99% of the follicles in goat ovaries will be atresia and only a few will eventually mature and ovulate. To investigate the potential long noncoding RNAs (lncRNAs) that regulate the expression of genes associated with follicular dominance or atresia, RNA-seq was performed on dominant follicles (DFs) and subordinate follicles (SFs) of granulosa cells from goats at the first follicular wave. A total of 92 differentially expressed lncRNAs and 676 differentially expressed mRNAs were detected in both types of follicles. The qRT-PCR results were consistent with the transcriptome sequencing data. Kyoto Encyclopedia of Genes and Genomes analysis of the differentially expressed mRNAs revealed that LHR and LDLR are associated with follicle dominance and are involved in the ovarian steroidogenesis pathway. The co-located mRNAs CALM2 and PPP1CA were significantly enriched during oocyte meiosis and in the cAMP and oxytocin signaling pathways. The co-expressed mRNAs were significantly enriched in the estrogen signaling pathway and in ovarian steroidogenesis and progesterone-mediated oocyte maturation. A co-expression network of lncRNAs, target genes, and differentially expressed genes was constructed. Follicle development-related genes, such as LDLR, NOTCH1, and FGF12, were included. These findings expand the lncRNA catalog and provide a basis for further studies on the mechanism of regulating follicular development in goats.

Single nucleotide polymorphisms in the 3′ UTR of follistatin-like 4 and scavenger receptor class B member 1 are associated with Dazu black goat litter size

The untranslated regions (UTR) of genes play crucial roles in regulating gene expression at the post-transcriptional level such as affecting mRNA stabilization. In this study, 26 single nucleotide polymorphisms (SNPs) and one deletion located in UTR were genotyped from 186 Dazu black goats via SNaPshot, and the correlation between genotype and litter size was analyzed. Results indicated that two SNP loci, SNP_chr17-20182525 and SNP_chr7-65652612, which are located at the 3′UTR of scavenger receptor class B member 1 and follistatin-like 4, are significantly (p<0.05) correlated with the litter size of first parity goats. SNP_chr7-65652612 is also significantly associated with the total litter size of first and second parity offspring (p<0.05). In conclusion, SNP_chr7-65652612 and SNP_chr17-20182525 have correlation with the litter size of Dazu black goat and they are potential genetic markers for litter size breeding.

Differential expression of FSHR and LHR genes and proteins during development of rabbit ovarian follicles

The development of an ovarian follicle is a complex process at the cellular and molecular level that is mainly regulated by follicle-stimulating hormone receptor (FSHR) and luteinizing hormone receptor (LHR). To elucidate the contribution of these receptors to ovarian follicle development, it is necessary to determine their expression profiles during this biological process. Therefore, this study aimed to investigate the relationship between ovarian development pattern and the differential ovarian expression pattern of FSHR and LHR genes as well as proteins at different developmental stages. Ovaries were collected from 30 New Zealand rabbits at day 0 (birth), week 2 (neonate), week 4 (cub), week 16 (maturity), and day 18 pregnancy. Ovarian histology, and gene as well as protein expression were determined using light microscopy, real-time PCR and western blotting, respectively. The results showed that the expression levels of FSHR mRNA and protein increased coincidently with age and the growth of ovarian follicles. The levels of LHR mRNA and protein remained low from the day of birth until week 4 and became significantly higher by week 16 coinciding with appearance of growing and antral follicles, which have a defined thecal layer. FSHR gene and protein expression decreased with pregnancy, whereas LHR increased, reaching a peak level during pregnancy. It can be concluded that changes in FSHR and LHR gene and protein expression could be related to the growth and development of follicles, indicating the regulatory role for these receptors in rabbit folliculogenesis.

Hypothalamic and ovarian transcriptome profiling reveals potential candidate genes in low and high egg production of white Muscovy ducks (Cairina moschata)

In China, the low egg production rate is a major challenge to Muscovy duck farmers. Hypothalamus and ovary play essential role in egg production of birds. However, there are little or no reports from these tissues to identify potential candidate genes responsible for egg production in White Muscovy ducks. A total of 1,537 laying ducks were raised; the egg production traits which include age at first egg (days), number of eggs at 300 d, and number of eggs at 59 wk were recorded. Moreover, 4 lowest (LP) and 4 highest producing (HP) were selected at 59 wk of age, respectively. To understand the mechanism of egg laying regulation, we sequenced the hypothalamus and ovary transcriptome profiles in LP and HP using RNA-Seq. The results showed that the number of eggs at 300 d and number of eggs at 59 wk in the HP were significantly more (P < 0.001) than the LP ducks. In total, 106.98G clean bases were generated from 16 libraries with an average of 6.68G clean bases for each library. Further analysis showed 569 and 2,259 differentially expressed genes (DEGs) were identified in the hypothalamus and ovary between LP and HP, respectively. The KEGG pathway enrichment analysis revealed 114 and 139 pathways in the hypothalamus and ovary, respectively which includes Calcium signaling pathway, ECM-receptor interaction, Focal adhesion, MAPK signaling pathway, Apoptosis and Apelin signaling pathways that are involved in egg production. Based on the GO terms and KEGG pathways results, 10 potential candidate genes (P2RX1, LPAR2, ADORA1, FN1, AKT3, ADCY5, ADCY8, MAP3K8, PXN, and PTTG1) were identified to be responsible for egg production. Further, protein-protein interaction was analyzed to show the relationship between these candidate genes. Therefore, this study provides useful information on transcriptome of hypothalamus and ovary of LP and HP Muscovy ducks.

Differential molecular and hormonal changes in oocytes, granulosa cells and follicular fluid of pregnant and non-pregnant camels

This study aimed to compare the expression of genes regulating follicles development, survival and steroid hormones secretion in oocytes and granulosa cells (GCs) and study the correlation between their expression and follicular fluid (FF) levels of progesterone (P4) in pregnant and non-pregnant camels. In total, 138 ovarian pairs from slaughtered camels were used. Gene expression and hormonal assay were determined using real-time polymerase chain reaction (PCR) and enzyme-linked immunosorbent assay (ELISA), respectively. The obtained results revealed that the number of follicles (3-8 mm) was significantly (P < 0.05) lower in pregnant, compared with non-pregnant, camels. P4 level in the FF was significantly (P < 0.05) higher in pregnant, compared with non-pregnant, camels. However, no significant (P > 0.05) difference was noticed in the oestradiol (E2) level. STAR, PTEN, IGF1 and BCL2 mRNA levels were significantly higher in GCs and significantly lower in oocytes of pregnant, compared with non-pregnant, camels. However, follicle-stimulating hormone receptor (FSHR) mRNA level was significantly lower in GCs and oocytes, and the BMP15 mRNA level was significantly lower in oocytes of pregnant, compared with non-pregnant, camels. P4 level in FF was positively correlated with STAR, PTEN, IGF1 and BCL2 mRNA levels in GCs and negatively correlated with BMP15 mRNA levels in oocytes and FSHR mRNA levels in GCs and oocytes of pregnant camels. It could be concluded that pregnancy-induced variations in oocytes and GC expression of BMP15, IGF1, FSHR, STAR, BCL2, and PTEN genes might be associated with a decrease in the number of follicles and an increase in the FF level of P4.

Analysis of long non-coding RNAs in epididymis of cattleyak associated with male infertility

Cattleyak (CY), is a cross breed between cattle and yak (YK), which display equal adaptability to the harsh environment as YK and much higher performances than YK. However, the CY is female fertile and male sterile. Previous studies were conducted on testes tissues to investigate the mechanism of male infertility in CY. There is no systematic research on genes, especially lncRNAs between CY and YK epididymis. In this study, Illumina Hiseq was performed to profile the epididymis transcriptome (lncRNA and mRNA) of CY and YK. In total 18859 lncRNAs were identified, from which lincRNAs 12458, antisense lncRNAs 2345, intronic lncRNAs 3101, and sense lncRNAs 955 respectively. We have identified 345 DE lncRNAs and 3008 DE mRNAs between YK and CY epididymis. Thirteen DEGs were validated by quantitative real-time PCR. Combing with DEG, 14 couples of lncRNAs and their target genes were both DE, and 6 of them including CCDC39, KCNJ16, NECTIN2, MRPL20, PSMC4, and DEFB112 show their potential infertility-related terms such as cellular motility, sperm maturation, sperm storage, cellular junction, folate metabolism, and capacitation. On the other hand, several down-regulated genes such as DEFB124, DEFB126, DEFB125, DEFB127, DEFB129, CES5A, TKDP1, CST3, RNASE9 and CD52 in CY compared to YK were involved in the immune response and sperm maturation. Therefore, comprehensive analysis for lncRNAs and their target genes may enhance our understanding of the molecular mechanisms underlying the process of sperm maturation in CY and may provide important resources for further research.

Evaluating the effect of TLR4-overexpressing on the transcriptome profile in ovine peripheral blood mononuclear cells

Background

Toll-like receptor 4 (TLR4) plays an important role in the elimination of Gram-negative bacteria infections and the initiation of antiinflammatory response. Using the technology of pronuclear microinjection, genetically modified (GM) sheep with TLR4 overexpression were generated. Previous studies have shown that these GM sheep exhibited a higher inflammatory response to Gram-negative bacteria infection than wild type (WT) sheep. In order to evaluate the gene expression of GM sheep and study the co-expressed and downstream genes for TLR4, peripheral blood mononuclear cells (PBMC) from TLR4-overexpressing (Tg) and wild type (WT) sheep were selected to discover the transcriptomic differences using RNA-Seq.

Result

An average of 18,754 and 19,530 known genes were identified in the Tg and WT libraries, respectively. A total of 338 known genes and 85 novel transcripts were found to be differentially expressed in the two libraries (p < 0.01). A differentially expressed genes (DEGs) enrichment analysis showed that the GO terms of inflammatory response, cell recognition, etc. were significantly (FDR < 0.05) enriched. Furthermore, the above DEGs were significantly (FDR < 0.05) enriched in the sole KEGG pathway of the Phagosome. Real-time PCR showed the OLR1, TLR4 and CD14 genes to be differentially expressed in the two groups, which validated the DEGs data.

Conclusions

The RNA-Seq results revealed that the overexpressed TLR4 in our experiment strengthened the ovine innate immune response by increasing the phagocytosis in PBMC.