lncRNA-miRNA-mRNA ceRNA Network Involved in Sheep Prolificacy: An Integrated Approach

Cell Cycle-Based Molecular Features via Synthetic Lethality and Non-Coding RNA Interactions in Cancer

BACKGROUND

The cell cycle, a critical and intricate biological process, comprises various phases, and its dysregulation plays a pivotal role in tumorigenesis and metastasis. The exploration of cell cycle-based molecular subtypes across pan-cancers, along with the application of synthetic lethality concepts, holds promise for advancing cancer therapies.

METHODS

A pan-cancer analysis was conducted to assess the cell cycle serves as a reliable signature for classifying molecular subtypes and to understand the potential clinical application of genes as potential drug targets based on synthetic lethality.

RESULTS

Molecular subtypes derived from cell cycle features in certain cancers, particularly kidney-related malignancies, exhibited distinct immune characteristics. Synthetic lethal interactions within the cell cycle pathway were common, with significant genetic interactions further identifying potential drug targets through the exploitation of genetic relationships with key driver genes. Additionally, miRNAs and lncRNAs may influence the cell cycle through miRNA:mRNA interactions and ceRNA networks, thereby enriching the genetic interaction landscape.

CONCLUSIONS

These findings suggest that the cell cycle pathway could serve as a promising molecular subtype signature to enhance cancer prognostication and offer potential targets for anticancer drug development through synthetic lethality.

Correction: Sadeghi et al. lncRNA-miRNA-mRNA ceRNA Network Involved in Sheep Prolificacy: An Integrated Approach. Genes 2022, 13, 1295

In the original publication [...].

Identification of ZNF652 as a Diagnostic and Therapeutic Target in Osteoarthritis Using Machine Learning

Purpose

Osteoarthritis (OA) is the most common degenerative joint disease. However, its etiology remains largely unknown. Zinc Finger Protein 652 (ZNF652) is a transcription factor implicated in various biological processes. Nevertheless, its role in OA has not been elucidated.

Methods

The search term "osteoarthritis" was utilized to procure transcriptome data relating to OA patients and healthy people from the Gene Expression Omnibus (GEO) database. Then a screening process was initiated to identify differentially expressed genes (DEGs). The DEGs were discerned using three distinct machine learning methods. The accuracy of these DEGs in diagnosing OA was evaluated using the Receiver Operating Characteristic (ROC) Curve. A competitive endogenous RNA (ceRNA) visualization network was established to delve into potential regulatory targets. The ZNF652 expression was confirmed in the cartilage of OA rats using quantitative reverse transcription polymerase chain reaction (qRT-PCR) and Western blotting (WB) and analyzed using an independent t-test.

Results

ZNF652 was identified as a DEG and exhibited the highest diagnostic value for OA according to the ROC analysis. The GO and KEGG enrichment analyses suggest that ZNF652 plays a vital role in OA development through processes including nitric oxide anabolism, macrophage proliferation, immune response, and the PI3K/Akt and the MAPK signaling pathways. The increased expression of ZNF652 in OA was validated in qRT-PCR (1.193 ± 0.005 vs 1.000 ± 0.005, p < 0.001) and WB (0.981 ± 0.055 vs 0.856 ± 0.026, p = 0.012) analysis.

Conclusion

ZNF652 was found to be related to OA pathogenesis and can potentially serve as a diagnostic and therapeutic target of OA. The underlying mechanism is that ZNF652 was related to nitric oxide anabolism, macrophage proliferation, various signaling pathways, and immune cells and their functions in OA. Nevertheless, the findings need to be confirmed in clinical trials and the molecular mechanism requires further study.

Molecular Investigation and Preliminary Validation of Candidate Genes Associated with Neurological Damage in Heat Stroke

Heat stroke (HS) is a severe medical condition characterized by a systemic inflammatory response that may precipitate multi-organ dysfunction, with a particular predilection for inducing profound central nervous system impairments. We aim to employ bioinformatics techniques for the retrieval and analysis of genes associated with heat stroke-induced neurological damage. We performed a comprehensive analysis of the GSE64778 dataset from the Sequence Read Archive, resulting in the identification of 1178 significantly differentially expressed genes (DEGs). We retrieved 2914 genes associated with heat stroke from the GeneCards database and 2377 genes associated with heat stroke from the Comparative Toxicogenomics Database (CTD). The intersection of the top 300 DEGs in the GSE64778 dataset intersected with the search results of GeneCards and CTD, yielding 25 final candidates for DEGs associated with heat stroke. Gene Ontology functional annotation results indicated that the target genes were mainly involved in apoptosis, stress response, and negative regulation of cellular processes and function in processes such as protein dimerization and protein binding. The Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis revealed a predominant enrichment of candidate target genes within the PI3K-AKT signaling pathway. Subsequent protein-protein interaction network analysis highlighted HSP90aa1 as a central gene, indicating its pivotal role by possessing the highest number of edges among the genes enriched in the PI3K-AKT signaling pathway. Quantitative reverse transcription-polymerase chain reaction analysis performed on blood samples from patients validated the expression of Hsp90aa1 in individuals exhibiting early neurological damage in HS, consistent with the findings from the mRNA bioinformatics analysis. Additionally, the bioinformatics analysis of the upstream microRNAs (miRNAs) regulating HSP90aa1 and the target miRNAs associated with candidate long non-coding RNAs (lncRNAs) identified three lncRNAs, eight miRNAs, and one mRNA in the regulatory network. The DIANA Tools database and algorithms were employed for pathway enrichment and correlation analysis, revealing a significant association between LOC102547734 and MIR-206-3p, with the latter being identified as a target binding site Moreover, the analysis unveiled a correlation between MIR-206-3p and HSP90aa1, implicating the latter as a potential target binding site within the regulatory network.

N6-methyladenosine demethylase ALKBH5 homologous protein protects against cerebral I/R injury though suppressing SNHG3-mediated neural PANoptosis: Involvement of m6A-related macromolecules in the diseases of nervous system

In order to address this gap in knowledge, the present study utilized both in vivo and in vitro models to investigate the role of the m6A demethylase ALKBH5 in protecting against cerebral I/R injury by inhibiting PANoptosis (Pytoptosis, Ppoptosis, and Necroptosis) in an m6A-dependent manner. They observed that ALKBH5, the predominant m6A demethylase, was downregulated in these models, while SNHG3 and PANoptosis-related proteins (ZBP1, AIM2, Cappase-3, Caspase-8, cleaved Caspase-1, GSDMD-N, and p-MLKL) were elevated. Additionally, both ALKBH5 overexpression and SNHG3-deficiency were found to ameliorate PANoptosis and injury induced by OGD/reperfusion and OGD/RX in both mice tissues and astrocyte cells. Further experiments demonstrated that ALKBH5 induced m6A-demethylation in SNHG3, leading to its degradation. Low expression of SNHG3, on the other hand, prevented the formation of the SNHG3-ELAVL1-ZBP1/AIM2 complex, which in turn destabilized ZBP1 and AIM2 mRNA, resulting in the downregulation of these PANoptosis-related genes. Ultimately, the rescue experiments provided evidence that ALKBH5 protected against PANoptosis in cerebral I/R injury models through the inhibition of SNHG3.This study sheds light on the intricate molecular mechanisms involved in the pathogenesis of cerebral I/R injury and highlights the potential of m6A-related genes as therapeutic targets in this condition.

Integrated Analysis of Transcriptome Profiles and lncRNA-miRNA-mRNA Competing Endogenous RNA Regulatory Network to Identify Biological Functional Effects of Genes and Pathways Associated with Johne's Disease in Dairy Cattle

Paratuberculosis or Johne's disease (JD), a chronic granulomatous gastroenteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP), causes huge economic losses and reduces animal welfare in dairy cattle herds worldwide. At present, molecular mechanisms and biological functions involved in immune responses to MAP infection of dairy cattle are not clearly understood. Our purpose was to integrate transcriptomic profiles and competing endogenous RNA (ceRNA) network analyses to identify key messenger RNAs (mRNAs) and regulatory RNAs involved in molecular regulation of peripheral blood mononuclear cells (PBMCs) for MAP infection in dairy cattle. In total, 28 lncRNAs, 42 miRNAs, and 370 mRNAs were identified by integrating gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. In this regard, we identified 21 hub genes (CCL20, CCL5, CD40, CSF2, CXCL8, EIF2AK2, FOS, IL10, IL17A, IL1A, IL1B, IRF1, MX2, NFKB1, NFKBIA, PTGS2, SOCS3, TLR4, TNF, TNFAIP3, and VCAM1) involved in MAP infection. Furthermore, eight candidate subnets with eight lncRNAs, 29 miRNAs, and 237 mRNAs were detected through clustering analyses, whereas GO enrichment analysis of identified RNAs revealed 510, 22, and 11 significantly enriched GO terms related to MAP infection in biological process, molecular function, and cellular component categories, respectively. The main metabolic-signaling pathways related to MAP infection that were enriched included the immune system process, defense response, response to cytokine, leukocyte migration, regulation of T cell activation, defense response to bacterium, NOD-like receptor, B cell receptor, TNF, NF-kappa B, IL-17, and T cell receptor signaling pathways. Contributions of transcriptome profiles from MAP-positive and MAP-negative sample groups plus a ceRNA regulatory network underlying phenotypic differences in the intensity of pathogenicity of JD provided novel insights into molecular mechanisms associated with immune system responses to MAP infection in dairy cattle.

The effect of LINC9137 targeting miR-140-3p-NKAIN3 signal axis on the development of goose testis sertoli cells

Sertoli cells (SC) are a type of important cells in the testes, which can provide transport proteins, regulatory proteins, growth factors, and other cytokines for the spermatogenic process. They participate in the regulation of the maturation and differentiation of spermatogenic cells and play an important supporting role in the migration, proliferation, and differentiation of germ cells at all levels in the testes. Previous studies found differential expression of LINC9137, miR-140-3p, and Sodium/Potassium Transporting ATPase Interacting 3 (NKAIN3) genesin high and low sperm motility goose testicular tissues. This study investigated the effects of the LINC9137-miR-140-3p-NKAIN3 signal axis on the proliferation and apoptosis of goose testicular sertoli cells at the cellular level, respectively. The results showed that through acridine orange staining, oil red O staining, Alkaline phosphatase (AKP) staining, and RT qPCR assay, it was comprehensively identified that the cultured testicular sertoli cells were purified in vitro. Through the dual luciferase activity detection test, it was found that LINC9137 has a targeted binding site with miR-140-3p and NKAIN3. In addition, this study found that overexpression of miR-140-3p significantly inhibited the expression of LINC9137 and NKAIN3 in sertoli cells, and their expression was significantly increased when miR-140-3p was interfered with. By measuring cell proliferation activity and apoptosis related gene expression, it was found that overexpression of LINC9137 decreased cell proliferation activity (P > 0.05), while the expression level of apoptosis factor Bcl2 Associated X Protein (Bax)/B-cell lymphoma-2 (Bcl2) increased (P > 0.05). On the contrary, when interfering with LINC9137, the cell proliferation activity of sertoli cells was significantly increased (P < 0.01), and the expression level of apoptosis factor Bax/Bcl2 was significantly reduced (P < 0.05); The effect of miR-140-3p on the proliferation and apoptosis of sertoli cells is opposite to that of LINC9137. Meanwhile, this study co transfected overexpressed LINC9137 and miR-140-3p plasmids into sertoli cells, and found that the effect of LINC9137 overexpression on supporting cell proliferation was weakened by miR-140-3p. This study elucidates the role and function of the LINC9137 miR-140-3p-NKAIN3 signaling axis in the development of goose testes and spermatogenesis, establishes a regulatory network related to spermatogenesis, and provides a theoretical basis for studying the genetic regulation of goose spermatogenesis.

Exploring MicroRNA biogenesis, applications and bioinformatics analysis in livestock: A comprehensive review

Small non-coding RNAs called microRNAs (miRNAs) control the expression of genes post-transcriptionally. Their correlation with commercial economic traits including milk, meat and egg production, as well as their effective role in animal productivity, fertility, embryo survival and disease resistance, make them significant in livestock research. The miRNAs exhibit distinct spatial and temporal expression patterns, offering insights into their functional roles within cells and tissues. Aberrant miRNA production can disrupt vital cellular processes and genetic networks, contributing to conditions like metabolic disorders and viral diseases. These short RNA molecules are present in extracellular fluids, displaying remarkable stability against RNA degradation enzymes and extreme environmental conditions. miRNAs preservation is facilitated through packaging in lipid vesicles or complex formation with RNA-binding proteins. Numerous studies have illuminated the roles of miRNAs in diverse physiological processes, including embryonic stem cell differentiation, haematopoietic stem cell proliferation and differentiation and the coordinated development of organ systems. The integration of miRNA profiling, next-generation sequencing and bioinformatics analysis paves the way for transformative advancements in livestock research and industry. The present review underscores the applications of miRNAs in livestock, showcasing their potential to improve breeding strategies, diagnose diseases and enhance our understanding of fundamental biological processes.

LncRNA12097.1 contributes to endometrial cell growth by enhancing YES1 activating β-catenin via sponging miR-145-5p

Despite previous investigations elucidating the regulatory mechanisms of long non-coding RNAs (lncRNAs) in endometrial function and reproductive disorders, the precise pathways through which lncRNAs impact endometrial functions and fertility remain unclear. In this study, we performed an expression profile analysis of lncRNAs in the endometrial tissue of Hu sheep with different prolificacy, identifying 13,707 lncRNAs. We discovered a bidirectional lncRNA, designated lncRNA12097.1, exhibiting significant up-regulation exclusively in the endometrium of Hu sheep with high fecundity. Functional analyses revealed lncRNA12097.1 significantly enhanced proliferation and cell cycle progression in both endometrial epithelial cell (EEC) and stromal cells (ESC), while inhibiting apoptosis in these cell types. Mechanistically, we demonstrated a directly interaction between lncRNA12097.1 and miR-145-5p, with YES proto-oncogene 1 (YES1) being identified as a validated target of miR-145-5p. Interference with lncRNA12097.1 resulted in suppressed cell growth through down-regulation of YES1 expression, which could be rescued by miR-145-5p. Furthermore, lncRNA12097.1 functions as a competitive endogenous RNA (ceRNA) for miR-145-5p in ESCs, sequestering miR-145-5p and preventing its binding to the 3'UTR of YES1 mRNA. This interaction led to increased expression of YES1 and subsequent activation of downstream β-catenin signaling, thereby promoting ESC growth in Hu sheep. These findings provide novel molecular insights into the mechanisms underlying prolificacy in sheep.

Regulatory Roles of Long Non-Coding RNAs Relevant to Antioxidant Enzymes and Immune Responses of Apis cerana Larvae Following Ascosphaera apis Invasion

Long non-coding RNAs (lncRNAs) play an essential part in controlling gene expression and a variety of biological processes such as immune defense and stress-response. However, whether and how lncRNAs regulate responses of Apis cerana larvae to Ascosphaera apis invasion has remained unclear until now. Here, the identification and structural analysis of lncRNAs in the guts of A. cerana worker larvae were conducted, and the expression profile of larval lncRNAs during the A. apis infection process was then analyzed, followed by an investigation of the regulatory roles of differentially expressed lncRNAs (DElncRNAs) in the host response. In total, 76 sense lncRNAs, 836 antisense lncRNAs, 184 intron lncRNAs, 362 bidirectional lncRNAs, and 2181 intron lncRNAs were discovered in the larval guts. Additionally, 30 known and 9 novel lncRNAs were potential precursors for 36 and 11 miRNAs, respectively. In the three comparison groups, 386, 351, and 272 DElncRNAs were respectively identified, indicating the change in the overall expression pattern of host lncRNAs following the A. apis invasion. Analysis of cis-acting effect showed that DElncRNAs in the 4-, 5-, and 6-day-old comparison groups putatively regulated 55, 30, and 20 up- and down-stream genes, respectively, which were involved in a series of crucial functional terms and pathways, such as MAPK signaling pathway, and cell process. Analysis showed that 31, 8, and 11 DElncRNAs as potential antisense lncRNAs may interact with 26, 8, and 9 sense-strand mRNAs. Moreover, investigation of the competing endogenous RNA (ceRNA) network indicated that 148, 283, and 257 DElncRNAs were putatively regulated. The expression of target genes by targeting corresponding DEmiRNAs included those associated with antioxidant enzymes and immune responses. These results suggested that DElncRNAs played a potential part in the larval guts responding to the A. apis infection through a cis-acting manner and ceRNA mechanisms. Our findings deepen our understanding of interactions between A. cerana larvae and A. apis and offer a basis for clarifying the DElncRNA-mediated mechanisms underlying the host response to fungal invasion.

Identification of Differentially Expressed lncRNAs in Response to Blue Light and Expression Pattern Analysis of Populus tomentosa Hybrid Poplar 741

Poplar is an important shelterbelt, timber stand, and city tree species that has been the focus of forestry research. The regulatory role of the long non-coding RNA molecule (lncRNA; length > 200 nt) has been a research hotspot in plants. In this study, seedlings of 741 poplar were irradiated with LED blue and white light, and the Illumina HiSeq 2000 sequencing platform was used to identify lncRNAs. |logFC| > 1 and p < 0.05 were considered to indicate differentially expressed lncRNAs, and nine differentially expressed lncRNAs were screened, the target genes of which were predicted, and three functionally annotated target genes were obtained. The differentially expressed lncRNAs were identified as miRNA targets. Six lncRNAs were determined to be target sites for twelve mRNAs in six miRNA families. LncRNAs and their target genes, including lncRNA MSTRG.20413.1-ptc-miR396e-5p-GRF9, were verified using quantitative real-time polymerase chain reaction analysis, and the expression patterns were analyzed. The analysis showed that the ptc-miR396e-5p expression was downregulated, while lncRNA MSTRG.20413.1 and GRF9 expression was upregulated, after blue light exposure. These results indicate that lncRNAs interact with miRNAs to regulate gene expression and affect plant growth and development.

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.

Integrated Comparative Transcriptome and circRNA-lncRNA-miRNA-mRNA ceRNA Regulatory Network Analyses Identify Molecular Mechanisms Associated with Intramuscular Fat Content in Beef Cattle

Intramuscular fat content (IMF), one of the most important carcass traits in beef cattle, is controlled by complex regulatory factors. At present, molecular mechanisms involved in regulating IMF and fat metabolism in beef cattle are not well understood. Our objective was to integrate comparative transcriptomic and competing endogenous RNA (ceRNA) network analyses to identify candidate messenger RNAs (mRNAs) and regulatory RNAs involved in molecular regulation of longissimus dorsi muscle (LDM) tissue for IMF and fat metabolism of 5 beef cattle breeds (Angus, Chinese Simmental, Luxi, Nanyang, and Shandong Black). In total, 34 circRNAs, 57 lncRNAs, 15 miRNAs, and 374 mRNAs were identified by integrating gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Furthermore, 7 key subnets with 16 circRNAs, 43 lncRNAs, 7 miRNAs, and 237 mRNAs were detected through clustering analyses, whereas GO enrichment analysis of identified RNAs revealed 48, 13, and 28 significantly enriched GO terms related to IMF in biological process, molecular function, and cellular component categories, respectively. The main metabolic-signaling pathways associated with IMF and fat metabolism that were enriched included metabolic, calcium, cGMP-PKG, thyroid hormone, and oxytocin signaling pathways. Moreover, MCU, CYB5R1, and BAG3 genes were common among the 10 comparative groups defined as important candidate marker genes for fat metabolism in beef cattle. Contributions of transcriptome profiles from various beef breeds and a competing endogenous RNA (ceRNA) regulatory network underlying phenotypic differences in IMF provided novel insights into molecular mechanisms associated with meat quality.

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.

Long non-coding RNA MSTRG.5970.28 regulates proliferation and apoptosis of goose follicle granulosa cells via the miR-133a-3p/ANOS1 pathway

The development of follicles in the ovaries is a critical determinant of poultry egg production. There are existing studies on the follicular development patterns in poultry, but the specific regulatory mechanisms still need further study. In a previous study, we identified long non-coding RNA (lncRNA) MSTRG.5970.28, anosmin 1 (ANOS1), and its predicted target miR-133a-3p that may be associated with goose ovary development. However, the function of MSTRG.5970.28 in goose granulosa cells and its regulatory mechanisms affecting granulosa cell proliferation and apoptosis have not been reported. In the present study, MSTRG.5970.28 and miR-133a-3p overexpression and interference vectors were constructed. Combined with reverse-transcription real-time quantitative PCR (RT-qPCR), a dual luciferase activity assay, Cell Counting Kit-8 (CCK-8), and flow cytometric analysis, we investigated the role of the MSTRG.5970.28-miR-133a-3p-ANOS1 axis in goose follicular granulosa cells and the associated regulatory mechanisms. MSTRG.5970.28 was found to be localized in the cytoplasm and its expression was influenced by reproductive hormones. The targeting relationship among MSTRG.5970.28, ANOS1, and miR-133a-3p were verified by a dual luciferase activity assay. CCK-8 and apoptosis assays showed that MSTRG.5970.28 inhibited the proliferation and promoted apoptosis of goose granulosa cells. The regulatory role of miR-133a-3p on granulosa cell proliferation and apoptosis was opposite to MSTRG.5970.28. We found that the proliferative and apoptotic effects of granulosa cells caused by MSTRG.5970.28 overexpression were attenuated by miR-133a-3p. MSTRG.5970.28 functions as a competitive endogenous RNA that regulates ANOS1 expression by sponging miR-133a-3p and thus exerts regulatory functions in granulosa cells. In sum, the present study identified lncRNA MSTRG.5970.28 as associated with goose ovary development, which affects the expression of ANOS1 by targeting miR-133a-3p, thereby influencing the proliferation and apoptosis of goose granulosa cells.

Epigenetic Regulation of miR-25 and Lnc107153 on Expression of Seasonal Estrus Key Gene CHGA in Sheep

Pituitary pars tuberalis (PT) plays an important role as the transmission center in the seasonal reproduction of animals. It helps convert external photoperiod signals into intrinsic seasonal reproduction signals. In sheep PT, specific expression patterns of several genes (including short photoperiod-induced gene CHGA and long photoperiod genes EYA3 and TSHβ) under different photoperiods are crucial characteristics during this signal transduction. Recent studies have revealed the role of epigenetics in regulating the expression of seasonal reproductive key genes. Therefore, we explored whether microRNAs and LncRNAs regulated the expressions of the above key genes. Firstly, the expression of miR-25 and CHGA showed a significant negative correlation in sheep PT. Results of the dual luciferase reporter assay and miR-25 overexpression indicated that miR-25 could inhibit the expression of CHGA by specifically binding to its 3'UTR region in pituitary cells. Then, expression negative correlation and dual luciferase reporter analyses were used to screen and identify the candidate LncRNA (Lnc107153) targeted by miR-25. Finally, the results of fluorescence in situ hybridization and Lnc107153 overexpression suggested that Lnc107153 and miR-25 were involved in the epigenetic regulation of CHGA expression. However, the expressions of EYA3 and TSHβ were not regulated by miRNAs. These results will provide new insights into the epigenetic regulatory network of key genes in sheep seasonal reproduction.