Gene Expression Profiling of Corpus luteum Reveals Important Insights about Early Pregnancy in Domestic Sheep

Role of chemokines in regulating luteal and uterine functions in pregnant cows

Pregnancy is intricately regulated by the interactions between various bioactive substances secreted by the conceptus, uterus, and corpus luteum (CL). Interferon-τ, synthesized and secreted by the conceptus, plays a central role in the interaction mechanism of maternal recognition in cows. Chemokines, chemotaxis mediators that are primarily secreted by immune cells, regulate various reproductive responses in various species. Although there are scattered reports on the potential roles of chemokines in the bovine CL and the uterus during the estrous cycle, there is little information on chemokines in these organs during pregnancy. Therefore, in this review, we discuss the possible physiological roles of chemokines in the CL and uterus of pregnant cows, focusing on our recent findings on chemokines and changes in their receptor expression in the CL and endometrium of cows at some stages of pregnancy.

New insights into how to induce and maintain embryonic diapause in the blastocyst

Embryonic diapause in mammals is a period of developmental pause of the embryo at the blastocyst stage. During diapause, the blastocyst has minimal cell proliferation, metabolic activity and gene expression. At reactivation, blastocyst development resumes, characterised by increases in cell number, biosynthesis and metabolism. Until recently, it has been unknown how diapause is maintained without any loss of blastocyst viability. This review focuses on recent progress in the identification of molecular pathways occurring in the blastocyst that can both cause and maintain the diapause state. A switch to lipid metabolism now appears essential to maintaining the diapause state and is induced by forkhead box protein O1. The forkhead box protein O transcription family is important for diapause in insects, nematodes and fish, but this is the first time a conclusive role has been established in mammals. Multiple epigenetic modifications are also essential to inducing and maintaining the diapause state, including both DNA and RNA methylation mechanisms. Finally, it now appears that diapause embryos, dormant stem cells and chemotherapeutic-resistant cancer cells may all share a universal system of quiescence.

Candidate genes associated with reproductive traits in rabbits

In the era of scientific advances and genetic progress, opportunities in the livestock sector are constantly growing. The application of molecular-based methods and approaches in farm animal breeding would accelerate and improve the expected results. The current work aims to comprehensively review the most important causative mutations in candidate genes that affect prolificacy traits in rabbits. Rabbits are a source of excellent-tasting meat that is high in protein and low in fat. Their early maturity and intensive growth are highly valued all over the world. However, improving reproductive traits and prolificacy in rabbits could be very tricky with traditional selection. Therefore, traditional breeding programs need new methods based on contemporary discoveries in molecular biology and genetics because of the complexity of the selection process. The study and implementation of genetic markers related to production in rabbits will help to create populations with specific productive traits that will produce the desired results in an extremely short time. Many studies worldwide showed an association between different genes and productive traits in rabbits. The study of these polymorphisms and their effects could be useful for molecular-oriented breeding, particularly marker-assisted selection programs in rabbit breeding.

Dynamic transcriptome analysis of Maiwa yak corpus luteum during the estrous cycle

Maiwa yak is a special breed of animal living on the Qinghai-Tibet Plateau, which has great economic value, but its fertility rate is low. The corpus luteum (CL) is a temporary tissue that plays a crucial role in maintaining the physiological cycle. However, little is known about the transcriptome profile in Maiwa yak CL. In the present study, the transcriptome of Maiwa yak CL at early (EYCL), middle (MYCL) and late-stages (LYCL) was studied employing high-throughput sequencing. A total of 25,922 transcripts were identified, including 22,277 known as well as 3,645 novel ones. Furthermore, 690 and 212 differentially expressed (DE) mRNAs were detected in the EYCL vs. MYCL and MYCL vs. LYCL groups, respectively. KEGG pathway enrichment analysis of DEGs illustrated that the most enriched pathway was PI3K-Akt pathway. Furthermore, twenty-six DEGs were totally found to be associated with different biological processes of CL development. One of these genes, PGRMC1, displayed a dynamical expression trend during the lifespan of yak CL. The knockdown of PGRMC1 in luteinized yak granulosa cells resulted in defective steroidogenesis. In conclusion, this study analyzed the transcriptome profiles in yak CL of different stages, and provided a novel database for analyzing the gene network in yak CL.HIGHLIGHTSThe manuscript analyzed the transcriptome profiles in yak CL during the estrous cycle.Twenty-six DEGs were found to be associated with the development or function of CL.One of the DEGs, PGRMC1, was found to be responsible for steroidogenesis in luteinized yak granulosa cells.

Unravelling transmission ratio distortion across the bovine genome: identification of candidate regions for reproduction defects

BACKGROUND

Biological mechanisms affecting gametogenesis, embryo development and postnatal viability have the potential to alter Mendelian inheritance expectations resulting in observable transmission ratio distortion (TRD). Although the discovery of TRD cases have been around for a long time, the current widespread and growing use of DNA technologies in the livestock industry provides a valuable resource of large genomic data with parent-offspring genotyped trios, enabling the implementation of TRD approach. In this research, the objective is to investigate TRD using SNP-by-SNP and sliding windows approaches on 441,802 genotyped Holstein cattle and 132,991 (or 47,910 phased) autosomal SNPs.

RESULTS

The TRD was characterized using allelic and genotypic parameterizations. Across the whole genome a total of 604 chromosomal regions showed strong significant TRD. Most (85%) of the regions presented an allelic TRD pattern with an under-representation (reduced viability) of carrier (heterozygous) offspring or with the complete or quasi-complete absence (lethality) for homozygous individuals. On the other hand, the remaining regions with genotypic TRD patterns exhibited the classical recessive inheritance or either an excess or deficiency of heterozygote offspring. Among them, the number of most relevant novel regions with strong allelic and recessive TRD patterns were 10 and 5, respectively. In addition, functional analyses revealed candidate genes regulating key biological processes associated with embryonic development and survival, DNA repair and meiotic processes, among others, providing additional biological evidence of TRD findings.

CONCLUSIONS

Our results revealed the importance of implementing different TRD parameterizations to capture all types of distortions and to determine the corresponding inheritance pattern. Novel candidate genomic regions containing lethal alleles and genes with functional and biological consequences on fertility and pre- and post-natal viability were also identified, providing opportunities for improving breeding success in cattle.

Identification of hub genes associated with follicle development in multiple births sheep by WGCNA

Sheep exhibit a distinct estrous cycle that includes four different phases: proestrus, estrus, late estrus, and luteal phase. As the estrous cycle repeats, follicular development regularly alternates. We thus investigated ovarian transcriptome data from each of the four phases using weighted gene co-expression network analysis (WGCNA) to identify modules, pathways, and genes essential to follicle growth and development. We clustered mRNA and long non-coding RNA (lncRNA) into different modules by WGCNA, and calculated correlation coefficients between genes and Stages of the estrous cycle. Co-expression of the black module (cor = 0.81, P<0.001) and the yellow module (cor = 0.61, P<0.04) was found to be critical for follicle growth and development. A total of 2066 genes comprising the black and yellow modules was used for functional enrichment. The results reveal that these genes are mainly enriched in Cell cycle, PI3K-Akt signaling pathway, Oocyte meiosis, Apoptosis, and other important signaling pathways. We also identified seven hub genes (BUB1B, MAD2L1, ASPM, HSD3B1, WDHD1, CENPA, and MXI1) that may play a role in follicle development. Our study may provide several important new markers allowing in depth exploration of the genetic basis for multiparous reproduction in sheep.

Temperature Incubation Influences Gonadal Gene Expression during Leopard Gecko Development

During development, sexual differentiation results in physiological, anatomical and metabolic differences that implicate not only the gonads but also other body structures. Sex in Leopard geckos is determined by egg incubation temperature. Based on the premise that the developmental decision of gender does not depend on a single gene, we performed an analysis on E. macularius to gain insights into the genes that may be involved in gonads' sexual differentiation during the thermosensitive period. All the genes were identified as differentially expressed at stage 30 during the labile phase of sex differentiation. In this way, the expression of genes known to be involved in gonadal sexual differentiation, such as WNT4, SOX9, DMRT1, Erα, Erβ, GnRH, P450 aromatase, PRLand PRL-R, was investigated. Other genes putatively involved in sex differentiation were sought by differential display. Our findings indicate that embryo exposure to a sex-determining temperature induces differential expression of several genes that are involved not only in gonadal differentiation, but also in several biological pathways (ALDOC, FREM1, BBIP1, CA5A, NADH5, L1 non-LTR retrotransposons, PKM). Our data perfectly fit within the new studies conducted in developmental biology, which indicate that in the developing embryo, in addition to gonadal differentiation, sex-specific tissue and metabolic polarization take place in all organisms.

Chronic superphysiologic AMH promotes premature luteinization of antral follicles in human ovarian xenografts

Anti-Müllerian hormone (AMH) is produced by growing ovarian follicles and provides a diagnostic measure of reproductive reserve in women; however, the impact of AMH on folliculogenesis is poorly understood. We cotransplanted human ovarian cortex with control or AMH-expressing endothelial cells in immunocompromised mice and recovered antral follicles for purification and downstream single-cell RNA sequencing of granulosa and theca/stroma cell fractions. A total of 38 antral follicles were observed (19 control and 19 AMH) at long-term intervals (>10 weeks). In the context of exogenous AMH, follicles exhibited a decreased ratio of primordial to growing follicles and antral follicles of increased diameter. Transcriptomic analysis and immunolabeling revealed a marked increase in factors typically noted at more advanced stages of follicle maturation, with granulosa and theca/stroma cells also displaying molecular hallmarks of luteinization. These results suggest that superphysiologic AMH alone may contribute to ovulatory dysfunction by accelerating maturation and/or luteinization of antral-stage follicles.

TIMP1 may affect goat prolificacy by regulating biological function of granulosa cells

Tissue inhibitor of metalloproteinase 1 (TIMP1) is associated with animal reproductive processes, such as follicular growth, ovulation, luteinization, and embryo development in mammals. The purposes of this study were to explore the expression and localization of TIMP1 in the ovarian tissues and determine the effect of TIMP1 on the function of granulosa cells and the association of TIMP1 with lambing-related genes of the goats. Immunohistochemical analysis showed that TIMP1 protein was strongly expressed by granulosa cells. Enzyme-linked immunosorbent assay (ELISA) results showed that TIMP1 overexpression promoted the secretion of estradiol of granulosa cells after 12, 24, and 48 h of transfection. Moreover, in vitro experiments indicated that TIMP1 had the ability to promote the cell proliferation and elevate the transcriptional levels of four genes associated with goat prolificacy, including BMPR-1B, BMP15, GDF9, and FSHB, in granulosa cells. In conclusion, TIMP1 could be an important molecule in regulating reproductive performance of the goats by affecting estrogen secretion and cell proliferation, as well as the expression of lambing-related genes of granulosa cells in the goats.

Identification of differentially expressed genes associated with precocious puberty by suppression subtractive hybridization in goat pituitary tissues

The aim of this study was to identify genes related to precocious puberty expressed in the pituitary of goats at different growth stages by suppression subtractive hybridization (SSH). The pituitary glands from Jining Gray (JG) goats (early puberty) and Liaoning Cashmere (LC) goats (late puberty) at 30, 90, and 180 days were used in this study. To identify differentially expressed genes (DEGs) in the pituitary glands, mRNA was extracted from these tissues, and SSH libraries were constructed and divided into the following groups: juvenile group (30-JG vs. 30-LC, API), puberty group (90-JG vs. 180-LC, BPI), and control group (90-JG vs. 90-LC, EPI). A total of 60, 49, and 58 DEGs were annotated by 222 Gene Ontology (GO) terms and 75 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Most of the DEGs were significantly enriched in GO terms related to 'structural constituent of ribosome', 'translation' and 'GTP binding', and numerous DEGs were also significantly enriched in KEGG terms related to the Jak-STAT signaling and oocyte meiosis pathways. Candidate genes associated with precocious puberty and sexual development were screened from the SSH libraries. These genes were analyzed to determine if they were expressed in the pituitary tissues of the goats at different growth stages and to identify genes that may influence the hypothalamic-pituitary-gonadal (HPG) axis. In this study, we found precocious puberty-related genes (such as PRLP0, EIF5A, and YWHAH) that may be interesting from an evolutionary perspective and that could be investigated for use in future goat breeding programs. Our results provide a valuable dataset that will facilitate further research into the reproductive biology of goats.

Comparative Transcriptome Analysis of Milk Somatic Cells During Lactation Between Two Intensively Reared Dairy Sheep Breeds

In dairy sheep industry, milk production dictates the value of a ewe. Milk production is directly related to the morphology and physiology of the mammary gland; both being designated targets of breeding strategies. Although within a flock breeding parameters are mutual, large differences in milk production among individual ewes are usually observed. In this work, we tested two of the most productive dairy sheep breeds reared intensively in Greece, one local the Chios breed and one foreign the Lacaune breed. We used transcriptome sequencing to reveal molecular mechanisms that render the mammary gland highly productive or not. While highly expressed genes (caseins and major whey protein genes) were common among breeds, differences were observed in differentially expressed genes. ENSOARG00000008077, as a member of ribosomal protein 14 family, together with LPCAT2, CCR3, GPSM2, ZNF131, and ASIP were among the genes significantly differentiating mammary gland's productivity in high yielding ewes. Gene ontology terms were mainly linked to the inherent transcriptional activity of the mammary gland (GO:0005524, GO:0030552, GO:0016740, GO:0004842), lipid transfer activity (GO:0005319) and innate immunity (GO:0002376, GO:0075528, GO:0002520). In addition, clusters of genes affecting zinc and iron trafficking into mitochondria were highlighted for high yielding ewes (GO:0071294, GO:0010043). Our analyses provide insights into the molecular pathways involved in lactation between ewes of different performances. Results revealed management issues that should be addressed by breeders in order to move toward increased milk yields through selection of the desired phenotypes. Our results will also contribute toward the selection of the most resilient and productive ewes, thus, will strengthen the existing breeding systems against a spectrum of environmental threats.

Identification of Genomic Regions Influencing N-Metabolism and N-Excretion in Lactating Holstein- Friesians

Excreted nitrogen (N) of dairy cows contribute to environmental eutrophication. The main N-excretory metabolite of dairy cows is urea, which is synthesized as a result of N-metabolization in the liver and is excreted via milk and urine. Genetic variation in milk urea (MU) has been postulated but the complex physiology behind the trait as well as the tremendous diversity of processes regulating the N-metabolism impede the consistent determination of causal regions in the bovine genome. In order to map the genetic determinants affecting N-excretion, MU and eight other N-excretory metabolites in milk and urine were assessed in a genome-wide association study. Therefore phenotypes of 371 Holstein- Friesians were obtained in a trial on a dairy farm under near commercial conditions. Genotype data comprised SNP information of the Bovine 50K MD Genome chip (45,613 SNPs). Significantly associated genomic regions for MU concentration revealed GJA1 (BTA 9), RXFP1, and FRY1 (both BTA 12) as putative candidates. For milk urea yield (MUY) a promising QTL on BTA 17 including SH3D19 emerged, whereas RCAN2, CLIC5, ENPP4, and ENPP5 (BTA 23) are suggested to influence urinary urea concentration. Minor N-fractions in milk (MN) may be regulated by ELF2 and SLC7A11 (BTA 17), whilst ITPR2 and MYBPC1 (BTA 5), STIM2 (BTA 6), SGCD (BTA 7), SLC6A2 (BTA 18), TMCC2 and MFSD4A (BTA 16) are suggested to have an impact on various non-urea-N (NUN) fractions excreted via urine. Our results highlight genomic regions and candidate genes for N-excretory metabolites and provide a deeper insight into the predisposed component to regulate the N-metabolism in dairy cows.

Downregulated luteolytic pathways in the transcriptome of early pregnancy bovine corpus luteum are mimicked by interferon-tau in vitro

Background

Maintenance of the corpus luteum (CL) beyond the time of luteolysis is essential for establishing pregnancy. Identifying the distinct features of early pregnancy CL remains unresolved, hence we analyzed here the transcriptome of CL on day 18 pregnant (P) and non-pregnant (NP) cows using RNA-Seq. CL of P cows expressed ISGs, verifying exposure to the pregnancy recognition signal, interferon-tau (IFNT), whereas the CL of NP cows had elevated luteal progesterone levels, implying that luteolysis had not yet commenced.

Results

The DEGs, IPA, and metascape canonical pathways, along with GSEA analysis, differed markedly in the CL of P cows from those of NP cows, at the same day of the cycle. Both metascape and IPA identified similar significantly enriched pathways such as interferon alpha/beta, sonic hedgehog pathway, TNFA, EDN1, TGFB1, and PDGF. However, type-1 interferon and sonic hedgehog pathways were positively enriched whereas most of the enriched pathways were downregulated in the P compared to NP samples. Thirty-four % of these pathways are known to be elevated by PGF2A during luteolysis. Notably, selective DEGs in luteinized granulosa cells were modulated by IFNT in vitro in a similar manner to their regulation in the CL of P cows.

Conclusion

This study unraveled the unique transcriptomic signature of the IFNT-exposed, early pregnancy CL, highlighting the abundance of downregulated pathways known to be otherwise induced during luteolysis. These and IFNT-regulated in vitro pregnancy-specific DEGs suggest that IFNT contributes to the characteristics and maintenance of early pregnancy CL.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07747-3.

Identification and characterization of miRNAs during early pregnancy in domestic sheep

MicroRNA resources in sheep are limited compared with those in other domesticated mammalian species. By sequencing small RNAs of sheep corpus luteum and endometrium, we have generated the largest amount of miRNA‐seq data and compiled the most comprehensive list thus far of miRNAs (n = 599) in sheep. Additionally, we observed a highly conserved maternally imprinted cluster of miRNAs on chromosome 18 homologous to that found on chromosome 14 in human and several other eutherian mammals.