Insight Into Pituitary lncRNA and mRNA at Two Estrous Stages in Small Tail Han Sheep With Different FecB Genotypes

Single-cell RNA-seq reveals the effects of the FecB mutation on the transcriptome profile in ovine cumulus cells

Genetic variations in the ovine ovulation rate, which are associated with the FecB mutation, provide useful models by which to explore the mechanisms regulating the development of mammalian antral follicles. In order to study the effects of the FecB mutation on cumulus cell differentiation, preovulatory follicles were aspirated and cumulus cells were isolated from three FecB genotypes (homozygous, heterozygous and wild type) of Small Tail Han (STH) sheep superstimulated with FSH. Transcriptome information from tens of thousands of cumulus cells was determined with the 10 × Genomics single-cell RNA-seq technology. Under the superovulation treatment, the observed number of preovulatory follicles in the ovaries of FecB carriers was still significantly higher than that in the wild-type (P < 0.05). The expression patterns of cumulus cells differed between FecB carriers and wild-type ewes. The screened cumulus cells could also be further divided into different cell clusters, and the differentiation states and fates of each group of cumulus cells also remained different, which supports the notion that heterogeneity in gene expression is prevalent in single cells. The oxidative phosphorylation pathway was significantly enriched in differentially expressed genes among the cell differentiation branch nodes of cumulus cells and among the differentially expressed genes of cumulus cells from the three genotypes. Combined with the important role of oxidative phosphorylation in the maturation of COCs, we suggest that the oxidative phosphorylation pathway of cumulus cells plays a crucial role in the differentiation process of cumulus cells and the mutation effect of the FecB gene.

Comparative whole-genome resequencing to uncover selection signatures linked to litter size in Hu Sheep and five other breeds

Hu sheep (HS), a breed of sheep carrying the FecB mutation gene, is known for its "year-round estrus and multiple births" and is an ideal model for studying the high fecundity mechanisms of livestock. Through analyzing and comparing the genomic selection features of Hu sheep and other sheep breeds, we identified a series of candidate genes that may play a role in Hu sheep's high fecundity mechanisms. In this study, we conducted whole-genome resequencing on six breeds and screened key mutations significantly correlated with high reproductive traits in sheep. Notably, the CC2D1B gene was selected by the fixation index (FST) and the cross-population composite likelihood ratio (XP-CLR) methods in HS and other five breeds. It was worth noting that the CC2D1B gene in HS was different from that in other sheep breeds, and seven missense mutations have been identified. Furthermore, the linkage disequilibrium (LD) analysis revealed a strong linkage disequilibrium in this specific gene region. Subsequently, by performing different grouping based on FecB genotypes in Hu sheep, genome-wide selective signal analysis screened several genes related to reproduction, such as BMPR1B and PPM1K. Besides, FST analysis identified functional genes related to reproductive traits, including RHEB, HSPA2, PPP1CC, HVCN1, and CCDC63. Additionally, a missense mutation was found in the CCDC63 gene and the haplotype was different between the high reproduction (HR) group and low reproduction (LR) group in HS. In summary, we discovered genetic differentiation among six distinct breeding sheep breeds at the whole genome level. Additionally, we identified a set of genes which were associated with reproductive performance in Hu sheep and visualized how these genes differed in different breeds. These findings laid a theoretical foundation for understanding genetic mechanisms behind high prolific traits in sheep.

Key mRNAs and lncRNAs of pituitary that affect the reproduction of FecB + + small tail han sheep

BACKGROUND

The pituitary directly regulates the reproductive process through follicle-stimulating hormone (FSH) and luteinizing hormone (LH). Transcriptomic research on the pituitaries of ewes with different FecB (fecundity Booroola) genotypes has shown that some key genes and lncRNAs play an important role in pituitary function and sheep fecundity. Our previous study found that ewes with FecB + + genotypes (without FecB mutation) still had individuals with more than one offspring per birth. It is hoped to analyze this phenomenon from the perspective of the pituitary transcriptome.

RESULTS

The 12 Small Tail Han Sheep were equally divided into polytocous sheep in the follicular phase (PF), polytocous sheep in the luteal phase (PL), monotocous sheep in the follicular phase (MF), and monotocous sheep in the luteal phase (ML). Pituitary tissues were collected after estrus synchronous treatment for transcriptomic analysis. A total of 384 differentially expressed genes (DEGs) (182 in PF vs. MF and 202 in PL vs. ML) and 844 differentially expressed lncRNAs (DELs) (427 in PF vs. MF and 417 in PL vs. ML) were obtained from the polytocous-monotocous comparison groups in the two phases. Functional enrichment analysis showed that the DEGs in the two phases were enriched in signaling pathways known to play an important role in sheep fecundity, such as calcium ion binding and cAMP signaling pathways. A total of 1322 target relationship pairs (551 pairs in PF vs. MF and 771 pairs in PL vs. ML) were obtained for the target genes prediction of DELs, of which 29 DEL-DEG target relationship pairs (nine pairs in PF vs. MF and twenty pairs in PL vs. ML). In addition, the competing endogenous RNA (ceRNA) networks were constructed to explore the regulatory relationships of DEGs, and some important regulatory relationship pairs were obtained.

CONCLUSION

According to the analysis results, we hypothesized that the pituitary first receives steroid hormone signals from the ovary and uterus and that VAV3 (Vav Guanine Nucleotide Exchange Factor 3), GABRG1 (Gamma-Aminobutyric Acid A Receptor, Gamma 1), and FNDC1 (Fibronectin Type III Domain Containing 1) played an important role in this process. Subsequently, the reproductive process was regulated by gonadotropins, and IGFBP1 (Insulin-like Growth Factor Binding Protein 1) was directly involved in this process, ultimately affecting litter size. In addition, TGIF1 (Transforming Growth Factor-Beta-Induced Factor 1) and TMEFF2 (Transmembrane Protein With EGF Like And Two Follistatin Like Domains 2) compensated for the effect of the FecB mutation and function by acting on TGF-β/SMAD signaling pathway, an important pathway for sheep reproduction. These results provided a reference for understanding the mechanism of multiple births in Small Tail Han Sheep without FecB mutation.

A novel lncRNA LOC105613571 binding with BDNF in pituitary promotes gonadotropin secretion by AKT/ERK-mTOR pathway in sheep associated with prolificacy

The pituitary is a vital endocrine organ for synthesis and secretion of gonadotropic hormones (FSH and LH), and the gonadotropin showed fluctuations in animals with different fecundity. Long non-coding RNAs (lncRNAs) have been identified as regulatory factors for the reproductive process. However, the profiles of lncRNAs and their roles involved in sheep fecundity remains unclear. In this study, we performed RNA-sequencing for the sheep pituitary gland associated with different fecundity, and identified a novel candidate lncRNA LOC105613571 targeting BDNF related to gonadotropin secretion. Our results showed that expression of lncRNA LOC105613571 and BDNF could be significantly upregulated by GnRH stimulation in sheep pituitary cells in vitro. Notably, either lncRNA LOC105613571 or BDNF silencing inhibited cell proliferation while promoted cell apoptosis. Moreover, lncRNA LOC105613571 knockdown could also downregulate gonadotropin secretion via inactivation AKT, ERK and mTOR pathway. In addition, co-treatment with GnRH stimulation and lncRNA LOC105613571 or BDNF knockdown showed the opposite effect on sheep pituitary cells in vitro. In summary, BDNF-binding lncRNA LOC105613571 in sheep regulates pituitary cell proliferation and gonadotropin secretion via the AKT/ERK-mTOR pathway, providing new ideas for the molecular mechanisms of pituitary functions.

CREB1 Is Involved in miR-134-5p-Mediated Endometrial Stromal Cell Proliferation, Apoptosis, and Autophagy

The successful establishment of endometrial receptivity is a key factor in ensuring the fertility of ewes and their economic benefits. Hu sheep have attracted attention due to their high fecundity and year-round estrus. In this study, we found that in the luteal phase, the uterine gland density, uterine coefficient, and number of uterine caruncles of high-fertility Hu sheep were higher than those of low-fertility Hu sheep. Thousands of differentially expressed genes were identified in the endometrium of Hu sheep with different fertility potential using RNA sequencing (RNA-Seq). Several genes involved in endometrial receptivity were screened using bioinformatics analysis. The qRT-PCR analysis further revealed the differential expression of cAMP reactive element binding protein-1 (CREB1) in the Hu sheep endometrium during the estrous cycle. Functionally, our results suggested that CREB1 significantly affected the expression level of endometrial receptivity marker genes, promoted cell proliferation by facilitating the transition from the G1 phase to the S phase, and inhibited cell apoptosis and autophagy. Moreover, we observed a negative linear correlation between miR-134-5p and CREB1 in the endometrium. In addition, CREB1 overexpression prevented the negative effect of miR-134-5p on endometrial stromal cell (ESC) growth. Taken together, these data indicated that CREB1 was regulated by miR-134-5p and may promote the establishment of uterine receptivity by regulating the function of ESCs. Moreover, this study provides new theoretical references for identifying candidate genes associated with fertility.

Identification of copy number variation in Tibetan sheep using whole genome resequencing reveals evidence of genomic selection

BACKGROUND

Copy number variation (CNV) is an important source of structural variation in the mammalian genome. CNV assays present a new method to explore the genomic diversity of environmental adaptations in animals and plants and genes associated with complex traits. In this study, the genome-wide CNV distribution characteristics of 20 Tibetan sheep from two breeds (10 Oula sheep and 10 Panou sheep) were analysed using whole-genome resequencing to investigate the variation in the genomic structure of Tibetan sheep during breeding.

RESULTS

CNVs were detected using CNVnator, and the overlapping regions of CNVs between individual sheep were combined. Among them, a total of 60,429 CNV events were detected between the indigenous sheep breed (Oula) and the synthetic sheep breed (Panou). After merging the overlapping CNVs, 4927 CNV regions (CNVRs) were finally obtained. Of these, 4559 CNVRs were shared by two breeds, and there were 368 differential CNVRs. Deletion events have a higher percentage of occurrences than duplication events. Functional enrichment analysis showed that the shared CNVRs were significantly enriched in 163 GO terms and 62 KEGG pathways, which were mainly associated with organ development, neural regulation, immune regulation, digestion and metabolism. In addition, 140 QTLs overlapped with some of the CNVRs at more than 1 kb, such as average daily gain QTL, body weight QTL, and total lambs born QTL. Many of the CNV-overlapping genes such as PPP3CA, SSTR1 and FASN, overlap with the average daily weight gain and carcass weight QTL regions. Moreover, VST analysis showed that XIRP2, ABCB1, CA1, ASPA and EEF2 differed significantly between the synthetic breed and local sheep breed. The duplication of the ABCB1 gene may be closely related to adaptation to the plateau environment in Panou sheep, which deserves further study. Additionally, cluster analysis, based on all individuals, showed that the CNV clustering could be divided into two origins, indicating that some Tibetan sheep CNVs are likely to arise independently in different populations and contribute to population differences.

CONCLUSIONS

Collectively, we demonstrated the genome-wide distribution characteristics of CNVs in Panou sheep by whole genome resequencing. The results provides a valuable genetic variation resource and help to understand the genetic characteristics of Tibetan sheep. This study also provides useful information for the improvement and breeding of Tibetan sheep in the future.

A meta-analysis of mRNA expression profiling studies in sheep with different FecB genotypes

The FecB mutation in the sheep BMPRIB is strongly correlated with high ovulation traits but its mechanism remains unclear. This study explored differentially expressed genes (DEGs) and their associated molecular mechanisms that may be involved in FecB mutation-induced high ovulation from the perspective of the hypothalamic-pituitary-gonadal (HPG) axis by conducting a systematic review and meta-analysis. The PubMed, EMBASE, CNKI, WanFang, and CBM databases were searched for eligible articles published before August 2022, focusing on mRNA sequencing of different tissues in the HPG axis in sheep with different FecB genotypes. A total of 6555 DEGs were identified from the analysis of six published articles and experimental results from our laboratory. The DEGs were screened by vote-counting rank and robust rank aggregation. Among these, in the follicular phase, FKBP5, CDCA7 and CRABP1 were upregulated in the hypothalamus. INSM2 was upregulated, while LDB3 was downregulated in the pituitary. CLU, SERPINA14, PENK, INHA and STAR were upregulated, while FERMT2 and NPY1R were downregulated in the ovary. On the HPG axis, TAC1 was upregulated and NPNT was downregulated. Many DEGs were found in sheep with different FecB genotypes. The genes FKBP5, CDCA7, CRABP1, INSM2, LDB3, CLU, SERPINA14, PENK, INHA, STAR, FERMT2, NPY1R, TAC1 and NPNT, may be associated with FecB mutation-induced high ovulation in different tissues. These candidate genes will further improve the mechanism of multiple fertility traits induced by the FecB mutation from the perspective of the HPG axis.

Two Retrotransposon Elements in Intron of Porcine BMPR1B Is Associated with Phenotypic Variation

It has been established that through binding to bone morphogenetic proteins (BMPs), bone morphogenetic protein receptor I B (BMPR1B) can mediate transforming growth factor β (TGF-β) signal transduction, and is involved in the regulation of several biological processes, such as bone and muscle formation and homeostasis, as well as folliculogenesis. Also known as FecB, BMPR1B has been reported as the major gene for sheep prolificacy. A number of previous studies have analyzed the relationship between single nucleotide polymorphisms (SNPs) in this gene and its related performance. In recent years, with the illustration of the effect of retrotransposon insertion on the expression of the proximal genes or phenotypic variation, retrotransposon insertion polymorphisms (RIPs) have been used as a novel type of molecular marker in the evaluation of evolution, population structure and breeding of plant and domestic animals. In this study, the RIPs in porcine BMPR1B gene were excavated, and thereafter verified using a comparative genome and polymerase chain reaction (PCR). The potential effects of phenotype, gene expression and functions related to RIPs were also explored. The results showed that 13 distinct RIPs were identified in introns of porcine BMPR1B. Among these, only BMPR1B-SINE-RIP9 and BMPR1B-LINE-RIP13 displayed a close relationship with the growth traits of Large White pigs. Moreover, the total number of BMPR1B-SINE+/+-RIP9 individuals born was found to be significantly higher than that of SINE−/− (p < 0.05). These two RIPs showed an obvious distribution pattern among Chinese indigenous breeds and Western commercial breeds. The expression of BMPR1B in ovaries of adult BMPR1B-SINE+/+-RIP9 Sushan pigs was found to be significantly higher in comparison to those of BMPR1B-SINE−/−-RIP9 (p < 0.05). SINE insertion of BMPR1B-SINE-RIP9 and LINE insertion of BMPR1B-LINE-RIP13 were observed to significantly increase the activity of Octamer binding transcription factor 4 (OCT4) minipromoter in CHO and C2C12 cells (p < 0.01). Therefore, these two RIPs could serve as useful molecular markers for modulating the growth or reproductive traits in assisted selection of pig breeding, while the mechanisms of the insertion function should be studied further.

Genetic Basis of Sexual Maturation Heterosis: Insights From Ovary lncRNA and mRNA Repertoire in Chicken

Sexual maturation is fundamental to the reproduction and production performance, heterosis of which has been widely used in animal crossbreeding. However, the underlying mechanism have long remained elusive, despite its profound biological and agricultural significance. In the current study, the reciprocal crossing between White Leghorns and Beijing You chickens were performed to measure the sexual maturation heterosis, and the ovary lncRNAs and mRNAs of purebreds and crossbreeds were profiled to illustrate molecular mechanism of heterosis. Heterosis larger than 20% was found for pubic space and oviduct length, whereas age at first egg showed negative heterosis in both crossbreeds. We identified 1170 known lncRNAs and 1994 putative lncRNAs in chicken ovary using a stringent pipeline. Gene expression pattern showed that nonadditivity was predominant, and the proportion of nonadditive lncRNAs and genes was similar between two crossbreeds, ranging from 44.24% to 49.15%. A total of 200 lncRNAs and 682 genes were shared by two crossbreeds, respectively. GO and KEGG analysis showed that the common genes were significantly enriched in the cell cycle, animal organ development, gonad development, ECM-receptor interaction, calcium signaling pathway and GnRH signaling pathway. Weighted gene co-expression network analysis (WGCNA) identified that 7 out of 20 co-expressed lncRNA-mRNA modules significantly correlated with oviduct length and pubic space. Interestingly, genes harbored in seven modules were also enriched in the similar biological process and pathways, in which nonadditive lncRNAs, such as MSTRG.17017.1 and MSTRG.6475.20, were strongly associated with nonadditive genes, such as CACNA1C and TGFB1 to affect gonad development and GnRH signaling pathway, respectively. Moreover, the results of real-time quantitative PCR (RT-qPCR) correlated well with the transcriptome data. Integrated with positive heterosis of serum GnRH and melatonin content detected in crossbreeds, we speculated that nonadditive genes involved in the GnRH signaling pathway elevated the gonad development, leading to the sexual maturation heterosis. We characterized a systematic landscape of ovary lncRNAs and mRNAs related to sexual maturation heterosis in chicken. The quantitative exploration of hybrid transcriptome changes lays foundation for genetic improvement of sexual maturation traits and provides insights into endocrine control of sexual maturation.