Analysis of pituitary transcriptomics indicates that lncRNAs are involved in the regulation of sheep estrus

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.

Comparative proteomics of ovaries elucidated the potential targets related to ovine prolificacy

Small Tail Han (STH) sheep, a unique Chinese breed, is recognized for its early maturity, year-round estrus, and prolificacy. However, the molecular mechanism of its high prolificacy has not been fully elucidated. The Proteomics approach is feasible and effective to reveal the proteins involved in the complex physiological processes of any organism. Given this, we performed the protein expression profiling of ovarian tissues during the luteal phase using polytocous STH sheep (litter size ≥2, three consecutive lambings) and monotocous STH sheep (litter size =1, three consecutive lambings) (PL vs. ML), and the follicular phase using polytocous STH sheep (litter size ≥2, three consecutive lambings) and monotocous STH sheep (litter size =1, three consecutive lambings) (PF vs. MF), respectively. Parallel Reaction Monitoring (PRM) was conducted to validate the differentially abundant proteins (DAPs). The tandem mass tag (TMT) quantitative proteomic results showed that a total of 5,237 proteins were identified, of which 49 and 44 showed differential abundance in the PL vs. ML and PF vs. MF groups, respectively. Enrichments analyses indicated that the DAPs including TIA1 cytotoxic granule-associated RNA-binding protein-like 1 (TIAL1), nicotinamide phosphoribosyltransferase (NAMPT), and cellular retinoic acid-binding protein 1 (CRABP1) were enriched at the luteal phase, while TIAL1, inhibin beta-a-subunit (A2ICA4), and W5PG55 were enriched at the follicular phase, potentially mediating reproductive processes in polytocous ewes. Furthermore, six DAPs were verified using PRM, confirming the accuracy of the TMT data acquired in this study. Together, our work expanded the database of indigenous sheep breeds and provided new ovarian candidate molecular targets, which will help in the study of the genetic mechanisms of ovine prolificacy.

Global Long Noncoding RNA Expression Profiling of MSTN and FGF5 Double-Knockout Sheep Reveals the Key Gatekeepers of Skeletal Muscle Development

Improving livestock and poultry growth rates and increasing meat production are urgently needed worldwide. Previously, we produced a myostatin (MSTN) and fibroblast growth factor 5 (FGF5) double-knockout (MF^-/-) sheep by CRISPR Cas9 system to improve meat production, and also wool production. Both MF^-/- sheep and the F1 generation (MF^+/-) sheep showed an obvious "double-muscle" phenotype. In this study, we identified the expression profiles of long noncoding RNAs (lncRNAs) in wild-type and MF^+/- sheep, then screened out the key candidate lncRNAs that can regulate myogenic differentiation and skeletal muscle development. These key candidate lncRNAs can serve as critical gatekeepers for muscle contraction, calcium ion transport and skeletal muscle cell differentiation, apoptosis, autophagy, and skeletal muscle inflammation, further revealing that lncRNAs play crucial roles in regulating muscle phenotype in MF^+/- sheep. In conclusion, our newly identified lncRNAs may emerge as novel molecules for muscle development or muscle disease and provide a new reference for MSTN-mediated regulation of skeletal muscle development.

Global gene expression analysis of the turkey hen hypothalamo-pituitary-gonadal axis during the preovulatory hormonal surge

The preovulatory hormonal surge (PS) consists of elevated circulating luteinizing hormone (LH) and progesterone levels and serves as the primary trigger for ovarian follicle ovulation. Increased LH and progesterone, produced by the pituitary and the granulosa layer of the largest ovarian follicle (F1), respectively, result from hypothalamic stimulation and steroid hormone feedback on the hypothalamo-pituitary-gonadal (HPG) axis. The hypothalamus, pituitary, F1 granulosa, and granulosa layer of the fifth largest follicle (F5) were isolated from converter turkey hens outside and during the PS and subjected to RNA sequencing (n = 6 per tissue). Differentially expressed genes were subjected to functional annotation using DAVID and IPA. A total of 12, 250, 1235, and 1938 DEGs were identified in the hypothalamus, pituitary, F1 granulosa, and F5 granulosa respectively (q<0.05, |fold change|>1.5, FPKM>1). Gene Ontology (GO) analysis revealed key roles for metabolic processes, steroid hormone feedback, and hypoxia induced gene expression changes. Upstream analysis identified a total of 4, 42, 126, and 393 potential regulators of downstream gene expression in the hypothalamus, pituitary, F1G, and F5G respectively, with a total of 63 potential regulators exhibiting differential expression between samples collected outside and during the PS (|z-score|>2). The results from this study serve to increase the current knowledge base surrounding the regulation of the PS in turkey hens. Through GO analysis, downstream processes and functions associated with the PS were linked to identified DEGs, and through upstream analysis, potential regulators of DEGs were identified for further analysis. Linking upstream regulators to the downstream PS and ovulation events could allow for genetic selection or manipulation of ovulation frequencies in turkey hens.

Discovery of alternatively spliced isoforms and long non-coding RNA in full length brain transcriptomes of anadromous Hilsa shad, Tenualosa ilisha (Hamilton, 1822)

BACKGROUND

Full length transcriptomes, achieved through long-read sequencing, along with the isoform analysis can reveal complexities in the gene expression profiles, as well as annotate the transcriptomes of non-model organisms.

METHODS AND RESULT

Full length transcripts of brain transcriptome of Tenualosa ilisha, Hilsa shad, were generated through PacBio single molecule real-time sequencing and were characterized. A total of 8.30 Gb clean reads were generated, with PacBio RSII, which resulted in 57,651 high quality consensus transcripts. After removing redundant reads, a total of 19,220 high-quality non-redundant transcripts and 17,341 full length ORF transcripts were classified to 7522 putative ortholog groups. Genes involved in various neural pathways were identified. In addition, isoform clusters and lncRNAs were discovered, along with Hilsa specific transcripts with coding frames and 29,147 SSRs in 944 transcripts (1141 annotated).

CONCLUSION

The present study provided, for the first time, a comprehensive view of the alternative isoforms of genes and transcriptome complexity in Hilsa shad brain and forms a rich resource for functional studies in brain of this anadromous fish.

Transcriptome analysis of messenger RNA and long noncoding RNA related to different developmental stages of tail adipose tissues of sunite sheep

The tail fat of sheep is the most typical deposited fat, and it can be widely used in human daily life, such as diet, cosmetics, and industrial raw materials. To understand the potential regulatory mechanism of different growth stages of tail fat in Sunite sheep, we performed high-throughput RNA sequencing to characterize the long noncoding RNA (lncRNA) and messenger RNA (mRNA) expression profiles of the sheep tail fat at the age of 6, 18, and 30 months. A total of 223 differentially expressed genes (DEGs) and 148 differentially expressed lncRNAs were found in the tail fat of 6-, 18-, and 30-month-old sheep. Based on functional analysis, we found that fat-related DEGs were mainly expressed at 6 months of age and gradually decreased at 18 and 30 months of age. The target gene prediction analysis shows that most of the lncRNAs target more than 20 mRNAs as their transregulators. Further, we obtained several fat-related differentially expressed target genes; these target genes interact with different differentially expressed lncRNAs at various ages and play an important role in the development of tail fat. Based on the DEGs and differentially expressed lncRNAs, we established three co-expression networks for each comparison group. Finally, we concluded that the development of the sheep tail fat is more active during the early stage of growth and gradually decreases with the increase in age. The mutual regulation of lncRNAs and mRNAs may play a key role in this complex biological process.

Enhancing Gonadotrope Gene Expression Through Regulatory lncRNAs

The world of long non-coding RNAs (lncRNAs) has opened up massive new prospects in understanding the regulation of gene expression. Not only are there seemingly almost infinite numbers of lncRNAs in the mammalian cell, but they have highly diverse mechanisms of action. In the nucleus, some are chromatin-associated, transcribed from transcriptional enhancers (eRNAs) and/or direct changes in the epigenetic landscape with profound effects on gene expression. The pituitary gonadotrope is responsible for activation of reproduction through production and secretion of appropriate levels of the gonadotropic hormones. As such, it exemplifies a cell whose function is defined through changes in developmental and temporal patterns of gene expression, including those that are hormonally induced. Roles for diverse distal regulatory elements and eRNAs in gonadotrope biology have only just begun to emerge. Here, we will present an overview of the different kinds of lncRNAs that alter gene expression, and what is known about their roles in regulating some of the key gonadotrope genes. We will also review various screens that have detected differentially expressed pituitary lncRNAs associated with changes in reproductive state and those whose expression is found to play a role in gonadotrope-derived nonfunctioning pituitary adenomas. We hope to shed light on this exciting new field, emphasize the open questions, and encourage research to illuminate the roles of lncRNAs in various endocrine systems.

Transcriptome Analysis Reveals Differentially Expressed Genes and Long Non-coding RNAs Associated With Fecundity in Sheep Hypothalamus With Different FecB Genotypes

Small-tailed Han sheep, with different FecB genotypes, manifest distinct ovulation rates and fecundities, which are due to differences in reproductive hormones secreted by the hypothalamic-pituitary-ovarian axis. Nevertheless, the function of the hypothalamus against a FecB mutant background on increasing ovulation rate is rarely reported. Therefore, we determined the expression profiles of hypothalamus tissue collected from six wild-type (WW) and six FecB mutant homozygous (BB) ewes at the follicular and luteal phases by whole-transcriptome sequencing. We identified 53 differentially expressed mRNAs (DEGs) and 40 differentially expressed long non-coding RNAs (DELs) between the two estrus states. Functional annotation analysis revealed that one of the DEGs, PRL, was particularly enriched in the hypothalamic function, hormone-related, and reproductive pathways. The lncRNA-target gene interaction networks and KEGG analysis in combination suggest that the lncRNAs LINC-676 and WNT3-AS cis-acting on DRD2 and WNT9B in different phases may induce gonadotropin-releasing hormone (GnRH) secretion. Furthermore, there were differences of regulatory elements and WNT gene family members involved in the follicular-luteal transition in the reproductive process between wild-type (WNT7A) and FecB mutant sheep (WNT9B). We combined the DEG and DEL data sets screened from different estrus states and genotypes. The overlap of these two sets was identified to select the mRNAs and lncRNAs that have major effects on ovulation. Among the overlapping molecules, seven DEGs and four DELs were involved in the follicular-luteal transition regulated by FecB mutation. Functional annotation analysis showed that two DEGs (FKBP5 and KITLG) were enriched in melanogenesis, oxytocin, and GnRH secretion. LINC-219386 and IGF2-AS were highly expressed in the BB ewes compared with WW ewes, modulating their target genes (DMXL2 and IGF2) to produce more GnRH during follicular development, which explains why mutated ewes produced more mature follicles. These results from expression profiling of the hypothalamus with the FecB mutation at different estrus states provide new insights into how the hypothalamus regulates ovulation under the effect of the FecB mutation.

Identification of the conserved long non-coding RNAs in myogenesis

BACKGROUND

Our understanding of genome regulation is ever-evolving with the continuous discovery of new modes of gene regulation, and transcriptomic studies of mammalian genomes have revealed the presence of a considerable population of non-coding RNA molecules among the transcripts expressed. One such non-coding RNA molecule is long non-coding RNA (lncRNA). However, the function of lncRNAs in gene regulation is not well understood; moreover, finding conserved lncRNA across species is a challenging task. Therefore, we propose a novel approach to identify conserved lncRNAs and functionally annotate these molecules.

RESULTS

In this study, we exploited existing myogenic transcriptome data and identified conserved lncRNAs in mice and humans. We identified the lncRNAs expressing differentially between the early and later stages of muscle development. Differential expression of these lncRNAs was confirmed experimentally in cultured mouse muscle C2C12 cells. We utilized the three-dimensional architecture of the genome and identified topologically associated domains for these lncRNAs. Additionally, we correlated the expression of genes in domains for functional annotation of these trans-lncRNAs in myogenesis. Using this approach, we identified conserved lncRNAs in myogenesis and functionally annotated them.

CONCLUSIONS

With this novel approach, we identified the conserved lncRNAs in myogenesis in humans and mice and functionally annotated them. The method identified a large number of lncRNAs are involved in myogenesis. Further studies are required to investigate the reason for the conservation of the lncRNAs in human and mouse while their sequences are dissimilar. Our approach can be used to identify novel lncRNAs conserved in different species and functionally annotated them.

Use of gene expression profile to identify potentially relevant transcripts to myofibrillar fragmentation index trait

The myofibrillar fragmentation index (MFI) is an indicative trait for meat tenderness. Longissimus thoracis muscle samples from the 20 most extreme bulls (out of 80 bulls set) for MFI (high (n = 10) and low (n = 10) groups) trait were used to perform transcriptomic analysis, using RNA Sequencing (RNA-Seq). An average of 24.616 genes was expressed in the Nellore muscle transcriptome analysis. A total of 96 genes were differentially expressed (p value ≤ 0.001) between the two groups of divergent bulls for MFI. The HEBP2 and BDH1 genes were overexpressed in animals with high MFI. The MYBPH and MYL6, myosin encoders, were identified. The differentially expressed genes were related to increase mitochondria efficiency, especially in cells under oxidative stress conditions, and these also were related to zinc and calcium binding, membrane transport, and muscle constituent proteins, such as actin and myosin. Most of those genes were involved in metabolic pathways of oxidation-reduction, transport of lactate in the plasma membrane, and muscle contraction. This is the first study applying MFI phenotypes in transcriptomic studies to identify and understand differentially expressed genes for beef tenderness. These results suggest that differences detected in gene expression between high and low MFI animals are related to reactive mechanisms and structural components of oxidative fibers under the condition of cellular stress. Some genes may be selected as positional candidate genes to beef tenderness, MYL6, MYBPH, TRIM63, TRIM55, TRIOBP, and CHRNG genes. The use of MFI phenotypes could enhance results of meat tenderness studies.