Whole-Transcriptome Analysis of LncRNAs Mediated ceRNA Regulation in Granulosa Cells Isolated From Healthy and Atresia Follicles of Chinese Buffalo

MiR-23a promotes autophagy of yak cumulus cells to alleviate apoptosis via the apoptosis signal-regulating kinase 1/c-Jun N-terminal kinase pathway

The ultimate fate of Graafian follicles is ovulation or atresia which relies on the highly coordinated processes of apoptosis and autophagy in ovarian cells. Long non-coding RNA maternally expressed gene 3 (LncRNA MEG3), miR-23a, and apoptosis signal-regulating kinase 1 (ASK1) are factors associated with autophagy. However, whether these factors can regulate autophagy in cumulus cells (CCs) of yak is unclear. Here, miR-23a overexpression upregulated the LC3-II/LC3-I ratio and Beclin1 abundance while reducing p62 accumulation (p < 0.05). The monodansylcadaverine assay exhibited a marked increase in punctate green fluorescence, and the GFP-LC3B displayed increased yellow fluorescence (p < 0.05). The opposite effect was observed for miR-23a inhibitors. Furthermore, miR-23a overexpression downregulated the abundance of ASK1 mRNA and total ASK1 protein (t-ASK1), whereas miR-23a inhibitors up-regulated them (p < 0.05). The effects of miR-23a overexpression on ASK1 phosphorylated protein at serine 845 (P-845), total JNK (c-Jun N-terminal kinase) (t-JNK) and the JNK phosphorylated protein (p-JNK) were similar to those of t-ASK1 but elicited the opposite effect on ASK1 phosphorylated protein at serine 967 (P-967) (p < 0.05). We further demonstrated that ASK1 expression can be silenced by small-interfering RNA (siRNA), which had no significant effect on t-JNK abundance (p > 0.05) but significantly suppressed the p-JNK expression (p < 0.05). Silencing ASK1 significantly improved Beclin1 abundance and the LC3-II/LC3-I ratio, but decreased p62 abundance (p < 0.05). An increase in yellow GFP-LC3B puncta and green MDC staining puncta were observed (p < 0.05). Overexpression of LncRNA MEG3 significantly increased the expression of t-ASK1, P-845, and JNK and decreased the abundance of P-967 and miR-23a (p < 0.05). In addition, miR-23a upregulation reduced the number of the TUNEL-positive cells, and the addition of 8 mM 3-methyladenine (3-MA) reversed this downregulation (p < 0.05). Similar trends were observed for the Bax/Bcl2 ratio and cleaved-caspase3 abundance. In summary, miR-23a promotes autophagy by inhibiting ASK1 abundance, which reduces apoptosis of yak CCs. This effect can be inhibited by LncRNA MEG3, which has implications for decreasing abnormal Graafian follicular atresia and maintaining development.

LncRNA MEG3 regulates ASK1/JNK axis-mediated apoptosis and autophagy via sponging miR-23a in granulosa cells of yak tertiary follicles

Apoptosis and autophagy in granulosa cells (GCs) are highly related to follicular development and atresia. It has also been reported that they are related to LncRNA MEG3, miR-23a and apoptosis signal-regulating kinase 1 (ASK-1). However, their relationship to follicular development and the extent to which follicle stimulating hormone (FSH) or luteinizing hormone (LH) can regulate this process remain unknown. Here, we found that ASK1 and JNK were expressed in the GCs of gonadotropin-dependent follicles, and those levels were significantly higher (p < 0.05) in yak Tertiary follicles compared to that of Secondary follicles and Graafian follicles. Then, the effect of LncRNA MEG3 / miR-23a on apoptosis and autophagy via ASK1/JNK (c-Jun N-terminal kinase) in yak GCs was studied. Overexpressing LncRNA MEG3 reduced miR-23a levels and p-967 protein expression, but enhanced ASK1 and JNK mRNA levels as well as t-ASK1, p-845, t-JNK, and p-JNK proteins levels. And Up-regulation of LncRNA MEG3 promoted apoptosis while attenuating autophagy. The targeting relationship between miR-23a and the binding sites of LncRNA MEG3 and ASK1 was also confirmed with the dual luciferase reporter assay. And, the relationship between LncRNA MEG3 and miR-23a was observed as a negative feedback regulation, and changes in LncRNA MEG3 and miR-23a levels can alter the expression of ASK1/JNK axis in yaks GCs. In addition, FSH (10 μg/mL) or LH (100 μg/mL) ability to reverse the effects of LncRNA MEG3 on miR-23a levels and ASK1/JNK axis-mediated apoptosis and autophagy was verified in yak GCs. This is significantly beneficial for decreasing abnormal follicular atresia for yaks tertiary follicles.

MiR-29c-5p regulates the function of buffalo granulosa cells to induce follicular atresia by targeting INHBA

MicroRNAs (miRNAs) are involved in many physiological processes such as signal transduction, cell proliferation and apoptosis. Many studies have shown that miRNAs can regulate the process of follicular development. Our previous studies found that the expression of miR-29c-5p in buffalo atretic follicles was much higher than that in healthy follicles, suggesting that this miRNA may participate in the process of buffalo follicular atresia. In this study, we aim to explore to the role and molecular mechanisms of miR-29c-5p on the functions of buffalo granulosa cells (GCs). GCs cultured in vitro were transfected with miR-29c-5p mimics and its inhibitor, respectively, and it was found that the mimics significantly increased the apoptotic rate of GCs. They also inhibited the proliferation of GCs and the secretion of steroid hormones. The effect of the inhibitor was opposite to that of the mimics. MiR-29c-5p was subsequently shown to target the inhibin subunit beta A, (INHBA). Overexpression of INHBA could promote the production of activin A and inhibin A, and then reverse the effect of miR-29c-5p on buffalo GCs. In conclusion, these results suggest that miR-29c-5p promotes apoptosis and inhibits proliferation and steroidogenesis by targeting INHBA in buffalo GCs. This may ultimately promote atresia in buffalo follicles.

Targeted metabolomics analysis of bile acids and cell biology studies reveal the critical role of glycodeoxycholic acid in buffalo follicular atresia

The follicular fluid of mammals has a high abundance of bile acids and these have proven to be closely related to the follicular atresia. However, the origin and content of bile acids in follicular fluid and its mechanisms on follicular atresia remain largely unknown. In this work, we analyzed the origin of bile acids in buffalo follicles by using cell biology studies, and quantified the subspecies of bile acids in follicular fluid from healthy follicles (HF) and atretic follicles (AF) by targeted metabolomics. The function of differential bile acids on follicular granulosa cells was also studied. The results showed that the bile acids transporters were abundantly expressed in ovarian tissues, but the rate-limiting enzymes were not, which was consistent with the inability of cultured follicular cells to convert cholesterol into bile acids. Targeted metabolomics analysis revealed thirteen differential subspecies of bile acids between HF and AF. The free bile acids were significant down-regulated and their conjugated forms were significantly up-regulated in AF as compared to HF. Finally, cell biological validation found a specific differentially conjugated bile acid, glycodeoxycholic acid (GDCA), which could promote follicular granulosa cell apoptosis and reduce steroid hormone secretion. In summary, our studies suggest that bile acids in buffalo follicles are transported from the blood rather than being synthesized within the follicles. The conjugated bile acids such as GDCA, accumulate in buffalo follicles, and may accelerate atresia by promoting apoptosis of granulosa cells and inhibiting steroid hormone production. These results will provide new clues for studying the physiological role and mechanism of bile acids involved in buffalo follicular atresia.

Whole-Genome-Based Web Genomic Resource for Water Buffalo (Bubalus bubalis)

Water buffalo (Bubalus bubalis), belonging to the Bovidae family, is an economically important animal as it is the major source of milk, meat, and drought in numerous countries. It is mainly distributed in tropical and subtropical regions with a global population of approximately 202 million. The advent of low cost and rapid sequencing technologies has opened a new vista for global buffalo researchers. In this study, we utilized the genomic data of five commercially important buffalo breeds, distributed globally, namely, Mediterranean, Egyptian, Bangladesh, Jaffrarabadi, and Murrah. Since there is no whole-genome sequence analysis of these five distinct buffalo breeds, which represent a highly diverse ecosystem, we made an attempt for the same. We report the first comprehensive, holistic, and user-friendly web genomic resource of buffalo (BuffGR) accessible at http://backlin.cabgrid.res.in/buffgr/, that catalogues 6028881 SNPs and 613403 InDels extracted from a set of 31 buffalo tissues. We found a total of 7727122 SNPs and 634124 InDels distributed in four breeds of buffalo (Murrah, Bangladesh, Jaffarabadi, and Egyptian) with reference to the Mediterranean breed. It also houses 4504691 SSR markers from all the breeds along with 1458 unique circRNAs, 37712 lncRNAs, and 938 miRNAs. This comprehensive web resource can be widely used by buffalo researchers across the globe for use of markers in marker trait association, genetic diversity among the different breeds of buffalo, use of ncRNAs as regulatory molecules, post-transcriptional regulations, and role in various diseases/stresses. These SNPs and InDelscan also be used as biomarkers to address adulteration and traceability. This resource can also be useful in buffalo improvement programs and disease/breed management.

circ01592 regulates unsaturated fatty acid metabolism through adsorbing miR-218 in bovine mammary epithelial cells

The composition of fatty acids plays a key role in regulating milk flavor and quality. Therefore, to improve the quality of milk, it is particularly important to study the regulatory mechanism of fatty acid metabolism in dairy cows. In this study, the expression profiles at non-lactation, peak-lactation, mid-lactation and late-lactation were constructed by high-throughput sequencing. Considering non-lactation as the control group and the other points as the experimental groups, the differentially expressed genes were screened. ELOVL5 was significantly upregulated and was selected for subsequent analyses. Bioinformatics prediction, a dual-luciferase assay, qPCR analysis and western blot analysis were used for verification. The results showed that ELOVL5 was a downstream target gene of miR-218 that regulated milk fat metabolism. A dual-luciferase assay and expression level analysis showed that circ01592 can directly bind to miR-218 and that overexpression of circ01592 (pcDNA-circ01592) significantly reduced the expression of miR-218 and enhanced the expression of ELOVL5, the target gene of miR-218 in BMECs. A functional study of BMECs showed that circ01592 promoted the synthesis of TAG and increased the content of UFA. The function of miR-218 was opposite to that of circ01592. Overall, the data showed that circ01592 promoted TAG synthesis and fatty acid composition by binding miR-218, alleviating the inhibitory effect of miR-218 on ELOVL5 expression. These mechanisms provide a new research approach and theoretical basis for improving milk quality.

Follicle-stimulating hormone regulates glycolysis of water buffalo follicular granulosa cells through AMPK/SIRT1 signalling pathway

Glycolysis in follicular granulosa cells (GCs) is the primary source of energy metabolism substrate of oocytes and is closely related to follicular development in mammals. Many physiological functions of GCs are regulated by follicle-stimulating hormone (FSH). In contrast, whether FSH regulates the glycolysis of GCs and its mechanism remains unclear. This study explored the correlation between FSH concentration and glycolysis level of GCs from different diameters of water buffalo follicles, and further explored the mechanism of FSH regulation in glycolysis in vitro cultured GCs. Results showed the variation trend of lactic acid concentration in follicular fluid and the expression level of glycolysis-related genes in GCs were consistent with the variation trend of FSH concentration in follicular fluid from follicles with different diameters. When GCs were treated with FSH in vitro, the expression level of glycolysis-related genes, lactate production and glucose uptake increased correspondingly (p < .05). Furthermore, we found that expression trend of AMPK/Sirtuin1 (SIRT1) pathway-related genes in GCs was consistent with the expression trend of glycolysis-related genes and was positively correlated with FSH concentrations in vivo or cultured in vitro. Activation of SIRT1 increased the expression level of glycolytic key proteins and lactic acid production in GCs, while inhibition of SIRT1 showed the opposite effect. In general, glycolysis in water buffalo GCs in vivo or cultured in vitro was positively correlated with FSH concentration. AMPK/SIRT1 pathway plays an important role in the regulation of FSH on glycolysis in GCs. Our findings will enrich the understanding of FSH regulating the development of water buffalo follicles.

Ovarian Aging: Role of Pituitary-Ovarian Axis Hormones and ncRNAs in Regulating Ovarian Mitochondrial Activity

The number of mitochondria in the oocyte along with their functions (e.g., energy production, scavenger activity) decline with age progression. Such multifaceted functions support several processes during oocyte maturation, ranging from energy supply to synthesis of the steroid hormones. Hence, it is hardly surprising that their impairment has been reported in both physiological and premature ovarian aging, wherein they are crucial players in the apoptotic processes that arise in aged ovaries. In any form, ovarian aging implies the progressive damage of the mitochondrial structure and activities as regards to ovarian germ and somatic cells. The imbalance in the circulating hormones and peptides (e.g., gonadotropins, estrogens, AMH, activins, and inhibins), active along the pituitary-ovarian axis, represents the biochemical sign of ovarian aging. Despite the progress accomplished in determining the key role of the mitochondria in preserving ovarian follicular number and health, their modulation by the hormonal signalling pathways involved in ovarian aging has been poorly and randomly explored. Yet characterizing this mechanism is pivotal to molecularly define the implication of mitochondrial dysfunction in physiological and premature ovarian aging, respectively. However, it is fairly difficult considering that the pathways associated with ovarian aging might affect mitochondria directly or by altering the activity, stability and localization of proteins controlling mitochondrial dynamics and functions, either unbalancing other cellular mediators, released by the mitochondria, such as non-coding RNAs (ncRNAs). We will focus on the mitochondrial ncRNAs (i.e., mitomiRs and mtlncRNAs), that retranslocate from the mitochondria to the nucleus, as active players in aging and describe their role in the nuclear-mitochondrial crosstalk and its modulation by the pituitary-ovarian hormone dependent pathways. In this review, we will illustrate mitochondria as targets of the signaling pathways dependent on hormones and peptides active along the pituitary/ovarian axis and as transducers, with a particular focus on the molecules retrieved in the mitochondria, mainly ncRNAs. Given their regulatory function in cellular activities we propose them as potential diagnostic markers and/or therapeutic targets.