Research | Open | Published:
HAS2-AS1 is a novel LH/hCG target gene regulating HAS2 expression and enhancing cumulus cells migration
Journal of Ovarian Researchvolume 12, Article number: 21 (2019)
The cumulus expansion process is one of the LH mediated ovulatory processes. Hyaluronan synthase 2 (HAS2) regulates the synthesis of hyaluronic acid, the main component of the cumulus expansion process. Recently, the lncRNA HAS2 antisense RNA 1 (HAS2-AS1) was identified in our global transcriptome RNA-sequencing of novel ovulation associated genes. The role of HAS2-AS1 in HAS2 regulation w.as studied previously with contradictive results in different models but not in the ovary. Taken together the induction of HAS2-AS1 and the important role of HAS2 in the cumulus expansion process, we hypothesize that HAS2-AS1 regulate HAS2 expression and function in the ovary. Therefore we undertook to study the expression, regulation, and possible functional role of HAS2-AS1 in the human ovary.
HAS2-AS1, located within the HAS2 gene that was highly regulated in our library. We found that HAS2-AS1 express mainly in cumulus cells (CCs). Furthermore, HAS2-AS1 showed low expression in immature CCs and a significant increase expression in mature CCs. Functional studies reveal that inhibition of HAS2-AS1 by siRNA caused decrease expression of HAS2. Furthermore, inhibition of HAS2-AS1 by siRNA results in decrease migration of granulosa cells.
Our results suggest that HAS2-AS1 is an LH/hCG target gene that plays a positive role in HAS2 expression and thus might play a role in regulating cumulus expansion and migration.
Critical to ovulation process is the cumulus-oocyte complex (COC) that undergos cumulus expansion. Genes that are induced following the LH surge and that are essential for proper expansion are hyaluronan synthase2 (HAS-2) and cyclooxygenase-2 (COX-2) that control the synthesis of hyaluronan and prostaglandins (such as PGE2) respectively [1, 2]. HAS2 is the main enzymes for synthesizing hyaluronic acid (HA). HA plays a critical role in vascular pathologies as it favors cell proliferation, migration, and the development of a pro-inflammatory state. In addition to hyaluronan structural backbone, there are several hyaluronan binding proteins, i.e. the proteoglycan versican , the serum-derived factor inter-αtrypsin inhibitor (IαI)  and the secreted protein tumor necrosis stimulated gene-6 (TSG-6) . Therefore, the production of hyaluronan is necessary for cumulus expansion . HAS2 was suggested as potential biomarkers in CCs for selection of oocytes and embryos in the IVF program .
Protein-coding genes account for about only 2% of the human genome, whereas the vast majority of transcripts are non-coding RNAs. Recent advances in RNA-sequencing technologies have led to the discovery of thousands of previously unannotated noncoding transcripts, including many long noncoding RNAs (lncRNAs). A growing volume of literature has proposed that lncRNAs are important factors in folliculogenesis and ovulation. However, the mechanism through which lncRNAs function remains largely unknown [8,9,10].
Recently, our lab applied global transcriptome RNA-sequencing approach to systematically identify novel ovulation associated genes and lncRNA . One of the lncRNAs which were differentially expressed is long non-protein coding HAS2 antisense RNA 1 (HAS2-AS1), located within HAS2 gene.
HAS2-AS1 has first discovered in human and mouse osteosarcoma cells . It was shown that HAS2-AS1 regulate the expression of HAS-2 in both directions. Initially, it was shown that the role of HAS2-AS1 is to suppress HAS2 mRNA levels and cell proliferation . Later, HAS2-AS1 was described in the renal proximal tubular epithelial cell , where its expression was found in correlation with HAS2 transcription, suggesting a positive effect of HAS2-AS1 on HAS2 mRNA expression. HAS2-AS1 was also described in human aortic smooth muscle cells where it was necessary for the induction of the HAS2 gene .
In view of the great importance of HAS2 in the ovulation process particularly in the stages of the cumulus expansion process and the unknown role of HAS2-AS1 in HAS2 regulation in the ovary, the aim of this study was to characterize HAS2-AS1 expression and regulation in vivo and in vitro in the human ovary, to elucidate its effect on HAS2 expression in human granulosa cells and to investigate its role in the ovulatory process.
As mention above, One of the lncRNA that show differential expression in our library of ovulation associated genes  was HAS2 long non-coding antisense 1 (HAS2-AS1). It showed low expression (below the threshold level) in immature CCs and a significant increase in expression in the mature CCs (36 reads, Fig. 1a).
To confirm the hCG dependent expression of HAS2-AS1 we compared its expression in CCs obtains from oocytes at different maturational stages. The results showed increased expression of HAS2-AS1 in CCs obtained from MII oocyte compared to CCs obtained from GV oocyte during stimulated IVF treatment (Fig. 1b).
The HAS2 role in cumulus expansion may suggest that its potential regulator, HAS2-AS1, will be expressed mainly in CCs. To characterize HAS2-AS1 expression in human ovarian follicle we compare HAS2-AS1 expression in CCs and mural granulosa cells (MGCs). The in vivo results indicate higher mRNA expression of HAS2-AS1 in CCs cells compared to MGCs (Fig. 1c).
The first hint for the effect of HAS2-AS1 on HAS2 expression was obtained in our library of ovulation associated genes  that show parallel induction of HAS2 and HAS2-AS1, suggest a positive effect. To confirm the positive effect of HAS2-AS1 on HAS2 expression, we transfected mural granulosa cell line (KGN) with HAS2-AS1 siRNA or control siRNA and 48 h later the cells were harvested. The results show that HAS2-AS1 levels (Fig. 2a) were reduced about 55% relative to control. Interestingly, HAS2 expression level (Fig. 2b) was also inhibited significantly. This result indicates and confirms that HAS2-AS1 has a positive effect on HAS2 expression.
To study the role of HAS2-AS1 on HAS2 function, we examine the functional role of HAS2-AS1 in cell migration that was previously described as a HAS2 function . Briefly, KGN cells, transfected with HAS2-AS1 siRNA or control siRNA, were grown to 95–100% confluency and starved overnight (12-14 h) in serum-free medium. A midline vertical scratch was performed and images were taken at 0 h and after 18 h, and 27 h of incubation. The results (Fig. 3) show that HAS2-AS1 siRNA transfected KGN cells migrate slower than the control siRNA transfected KGN cells. This result suggests a role for HAS2-AS1 in cell migration.
In summary, we showed in vivo (and ex vivo) that HAS2-AS1 is an hCG-induce gene and is positively effect HAS2 expression and function.
Small noncoding RNAs (ncRNAs) such as siRNAs, and miRNAs, are highly conserved and involved in transcriptional and post-transcriptional gene silencing through specific base pairing with their targets. In contrast, the lncRNAs mechanisms of gene expression regulation are numerous due to their diverse structural conformations, biochemical properties, and specific subcellular localization, resulting in both induction and suppression of their target gene [16,17,18].
Recently there has been a growing interest in the role of noncoding RNAs in the regulation of reproduction. Men and colleges  showed that StAR is a novel target of the microRNA let-7, which itself is regulated by the lncRNA H19. Overexpression of H19 stimulates StAR expression by antagonizing let-7, which inhibits StAR at the post-transcriptional level. Nakagawa  showed that the lncRNA Neat1 is required for corpus luteum formation and the establishment of pregnancy in mice. However, the information about the intriguing role of lncRNA in the ovulation process is still scarce and further studies are needed to learn about their role in the complex ovulatory process.
The fact that lncRNAs regulatory mechanisms are diverse may explain the opposite regulatory effects of HAS2-AS1 on HAS2 expression reported in different tissues. Our novel results show that HAS2-AS1 has a positive effect on the induction of HAS2 in human granulosa cells. The positive effect on the induction of HAS2 was shown in other tissues as well [21, 22] however, Chao et al.  demonstrated that HAS2-AS1 inhibits HAS2 expression in human osteosarcoma cells.
HAS2-AS1 is located on the opposite strand of the HAS2 gene, phenomena known as natural antisense (NAT) lncRNA. Although HAS2-AS1 share some sequence with HAS2, its action through direct binding is not likely because the HAS2-AS1 is expressed at much lower levels than HAS2 (data not shown). Our results showed that HAS2-AS1 positively regulates HAS2 expression. The molecular mechanism for this phenomenon might be similar to previous reports that transcription of the antisense RNA stabilizes or augments HAS2 mRNA expression via RNA/mRNA heteroduplex formation  and by altering the chromatin structure around the HAS2 proximal promoter . It is clear that additional studies are needed to understand the regulatory mechanisms of HAS2-AS1 on HAS2 expression and hence to understand the different effect in different tissues.
In this study, we suggest for the first time that HAS2-AS expression is induced by LH/hCG. It is known that HAS2 expression is induced by LH/hCG  however induced expression of HAS2-AS by LH/hCG is novel. This new finding raises intriguing questions regarding how LH/hCG effects HAS2 expression and the possibility of the synergistic influence of both LH/hCG and HAS2-AS on HAS2 expression and function.
During the hyaluronan synthesis, the cumulus cells lose contact with one another and migrate outward from the oocyte along the hyaluronan polymeric backbone, a process called expansion. Although the molecular mechanisms and factors that control COC expansion have been studied extensively, much less is known about what factors control cumulus cell movement. It was shown that calpain activity was associated with cumulus cell detachment and movement by mechanisms that did not alter HAS2 expression . Hyaluronic acid was shown to facilitate migration of several cells in other systems including human keratinocytes , human dental pulp stem cells , breast cancer cell lines , and human T-lymphocyte .
The role of cumulus expansion includes facilitates the extrusion of the oocyte through the ruptured follicle wall during ovulation and assists its capture by the oviductal fimbria and its entry into the oviduct. Cumulus expansion was also shown to be essential not only for ovulation but also for in vivo fertilization, and perhaps more specifically for sperm capacitation [30, 31].
We assume that inhibiting HAS2 activity using HAS2-AS1 as a target, will result in inhibition of the processes as mentioned above and therefore reducing ovulation and fertilization efficiency.
Ovulation includes several significant processes such as cumulus expansion, steroidogenesis and cytoplasmic and nuclear maturation of the oocyte. Interestingly, these processes are independent, and inhibition of one of these processes does not obligatory block the others. Accordingly, we speculated that by targeting cumulus expansion specifically, we might achieve ovulation inhibition without harming the other processes (such as steroidogenesis), and therefore less systemic side effects.
In summary, we show here for the first time that the human natural HAS-AS1 is an LH/hCG target gene that regulates HAS2 expression and function in human granulosa cells. It can be assumed that the use of this finding can lead to the development of targeted contraception with reduced systemic side effects.
Materials and methods
Normo-ovulatory young women (< 37 years of age) undergoing IVF due to male factor infertility or genetic indication for pre-implantation genetic diagnosis were selected for this study. Subjects carrier of Fragile X disorder, endometriosis and polycystic ovary syndrome (PCOS) were excluded. Ovarian stimulation was carried out as previously described . Briefly, short antagonist protocol was used, consisted of ovarian suppression with GnRH antagonists (0.25 mg/day, Cetrorelix, Cetrotide; Merck Serono International) and controlled ovarian hyperstimulation with human menopausal gonadotropin (HMG; Menopur) or recombinant FSH, either Gonal-F; Merck Serono or Puregon, MSD. The initial dose used was dependent upon age, body mass index, and previous IVF treatment history. When three or more follicles exceeded18 mm in diameter, 250 μg of hCG (Ovitrelle; Merck Serono) was administered to trigger ovulation. Transvaginal follicular aspiration was performed 36 h later with ultrasound guidance.
Cumulus granulosa cell collection and grouping
CCs were obtained during oocyte denudation for intracytoplasmic sperm injection (ICSI) procedures. After oocyte retrieval, CCs of each oocyte were removed using hyaluronidase (SAGE) and a glass denudation pipette (Swemed). The CCs were assigned according to the related oocyte maturation status, GV, MI, or MII. The percentage of GV oocytes was 5–10%. The CCs were washed in PBS and centrifuged at 5000×g for 5 min at room temperature. The CCs from 3 to 4 subjects were pooled to generate a single replicate and the resulting pellets were stored at − 80 °C until subjected to total RNA isolation.
Mural granulosa cell collection and grouping
MGCs were collected from aspirated follicular fluid during IVF procedures and re-suspended in a phosphate-buffered solution (PBS; Sigma-Aldrich). After allowing the cells to settle by gravity for a few minutes, the top portion of the medium was repeatedly aspirated until the medium proved clear. The cells were then centrifuged at 1000 rpm for 5 min at room temperature. Purified MGCs from large follicles (> 10 mm) of 3–4 subjects were pooled to generate a single replicate. All replicates were store at − 80 °C until subjected to total RNA purification.
RNA extraction and qPCR
Total RNA was extracted from MGCs using a Mini/Micro RNA Isolation I kit (Zymo Research) according to the manufacturer’s instructions. RNA purity and concentration were assessed using a NanoDrop spectrophotometer (NanoDrop 2000C, Thermo Scientific). Total RNA (25 ng) from each sample was used for cDNA synthesis with a high capacity reverse transcription kit (Applied Biosystems) according to the manufacturer’s instructions in a 10 μl total volume reaction. mRNA concentrations were analyzed by real-time PCR using the StepOnePlus real-time PCR system (Applied Biosystems). The real-time PCR mix contained 1 μl of cDNA, fast SYBR Green Master Mix (Applied Biosystems), and specific primers for HAS2 or HAS2-AS1 and β-actin (housekeeping gene) in a total volume of 10 μl. Cycling parameters were: 1 cycle at 95 °C for 20 s and 40 cycles each at 95 °C for 3 s and at 60 °C for 30 s. A melting curve analysis was performed at the end of each run to ensure a single amplicon. All samples were run in duplicates. qPCR results were analyzed with StepOne software. Relative gene expression was calculated using the delta-delta Ct method. Details of the primers are shown in Table 1.
KGN cells  were transfected, using Dharmafect1 (Dharmacon Inc.), with HAS2-AS1 siRNA (SMARTpool: Lincode Human HAS2-AS1 siRNA; NR_002835; R-187921-00-0005) or scrambled control (Dharmacon Lincode Non-targeting siRNA #1) according to the manufacture instruction. After 48 h cells were harvested and subjected to total RNA purification, cDNA and qRT-PCR with HAS2 and HAS2-AS1 primers.
In vitro scratch assay
KGN cells  were transfected, using Dharmafect1 (Dharmacon Inc.), with HAS2-AS1 siRNA (SMARTpool: Lincode Human HAS2-AS1 siRNA; NR_002835; R-187921-00-0005) or scrambled control (Dharmacon Lincode Non-targeting siRNA #1) according to the manufacture instruction. After 48 h cells were grown at 95–100% confluency, and starve overnight (12-14 h) in serum-free medium. During scratching assay, FBS in the medium was set to 0.2% to inhibit apoptosis and necrosis, and at the same time, inhibit proliferation. A midline vertical scratch was performed using a 10 μl tip with a blunt surface. Pictures were taken immediately (oh), after 18 h, and 27 h of incubation.
Hizaki H, Segi E, Sugimoto Y, Hirose M, Saji T, Ushikubi F, et al. Abortive expansion of the cumulus and impaired fertility in mice lacking the prostaglandin E receptor subtype EP(2). Proc Natl Acad Sci U S A. 1999;96(18):10501–6 PubMed PMID: 10468638. Epub 1999/09/01. eng.
Richards JS, Pangas SA. The ovary: basic biology and clinical implications. J Clin Invest. 2010;120(4):963–72 PubMed PMID: 20364094. Epub 2010/04/07. eng.
Russell DL, Robker RL. Molecular mechanisms of ovulation: co-ordination through the cumulus complex. Hum Reprod Update. 2007;13(3):289–312 PubMed PMID: 17242016. Epub 2007/01/24. eng.
Hess KA, Chen L, Larsen WJ. Inter-alpha-inhibitor binding to hyaluronan in the cumulus extracellular matrix is required for optimal ovulation and development of mouse oocytes. Biol Reprod. 1999;61(2):436–43. PubMed PMID: 10411524. Epub 1999/07/20. eng.
Fulop C, Szanto S, Mukhopadhyay D, Bardos T, Kamath RV, Rugg MS, et al. Impaired cumulus mucification and female sterility in tumor necrosis factor-induced protein-6 deficient mice. Development. 2003;130(10):2253–61 PubMed PMID: 12668637. Epub 2003/04/02. eng.
Richards JS. Ovulation: new factors that prepare the oocyte for fertilization. Mol Cell Endocrinol. 2005;234(1–2):75–9 PubMed PMID: 15836955. Epub 2005/04/20. eng.
Shepel EA, Voznesenskaya TY, Blashkiv TV, Yanchii RI. CUMULUS CELL GENES AS POTENTIAL BIOMARKERS OF OOCYTE AND EMBRYO DEVELOPMENTAL COMPETENCE. Fiziolohichnyi zhurnal (Kiev, Ukraine : 1994). 2016;62(1):107–13 PubMed PMID: 29537212. Epub 2016/01/01. eng.
Kurokawa R. Long noncoding RNA as a regulator for transcription. Prog Mol Subcell Biol. 2011;51:29–41 PubMed PMID: 21287132.
Yan B, Wang Z. Long noncoding RNA: its physiological and pathological roles. DNA Cell Biol. 2012;31(Suppl 1):S34–41 PubMed PMID: 22612272.
Kim ED, Sung S. Long noncoding RNA: unveiling hidden layer of gene regulatory networks. Trends Plant Sci. 2012;17(1):16–21 PubMed PMID: 22104407.
Yerushalmi GM, Salmon-Divon M, Yung Y, Maman E, Kedem A, Ophir L, et al. Characterization of the human cumulus cell transcriptome during final follicular maturation and ovulation. Mol Hum Reprod. 2014;20(8):719–35 PubMed PMID: 24770949. Epub 2014/04/29. eng.
Chao H, Spicer AP. Natural antisense mRNAs to hyaluronan synthase 2 inhibit hyaluronan biosynthesis and cell proliferation. J Biol Chem. 2005;280(30):27513–22 PubMed PMID: 15843373. Epub 2005/04/22. eng.
Michael DR, Phillips AO, Krupa A, Martin J, Redman JE, Altaher A, et al. The human hyaluronan synthase 2 (HAS2) gene and its natural antisense RNA exhibit coordinated expression in the renal proximal tubular epithelial cell. J Biol Chem. 2011;286(22):19523–32 PubMed PMID: 21357421. Pubmed Central PMCID: PMC3103331. Epub 2011/03/02. eng.
Vigetti D, Deleonibus S, Moretto P, Bowen T, Fischer JW, Grandoch M, et al. Natural antisense transcript for hyaluronan synthase 2 (HAS2-AS1) induces transcription of HAS2 via protein O-GlcNAcylation. J Biol Chem. 2014;289(42):28816–26 PubMed PMID: 25183006. Pubmed Central PMCID: PMC4200242. Epub 2014/09/04. eng.
Akison LK, Alvino ER, Dunning KR, Robker RL, Russell DL. Transient invasive migration in mouse cumulus oocyte complexes induced at ovulation by luteinizing hormone. Biol Reprod. 2012;86(4):125 PubMed PMID: 22238281. Epub 2012/01/13. eng.
Jandura A, Krause HM. The New RNA World: Growing Evidence for Long Noncoding RNA Functionality. Trends Genet. 2017;33(10). Epub 2017 Sep 1.
Liu SJ, et al. Science. 2016. https://doi.org/10.1126/science.aah7111. PMID:27980086.
Long Y, Wang X, Youmans DT, Cech TR. How do lncRNAs regulate transcription? Sci Adv. 2017;3(9) Epub 2017 Sep 27.
Men Y, Fan Y, Shen Y, Lu L, Kallen AN. The steroidogenic acute regulatory protein (StAR) is regulated by the H19/let-7 Axis. Endocrinology. 2017;158(2):402–9 PubMed PMID: 27813675. Pubmed Central PMCID: 5413078.
Nakagawa S, Shimada M, Yanaka K, Mito M, Arai T, Takahashi E, et al. The lncRNA Neat1 is required for corpus luteum formation and the establishment of pregnancy in a subpopulation of mice. Development. 2014;141(23):4618–27 PubMed PMID: 25359727. Pubmed Central PMCID: 4302932.
Kedem A, Yung Y, Yerushalmi GM, Haas J, Maman E, Hanochi M, et al. Anti Mullerian hormone (AMH) level and expression in mural and cumulus cells in relation to age. J Ovarian Res. 2014;7:113 PubMed PMID: 25500128. Pubmed Central PMCID: 4269874.
Yung Y, Aviel-Ronen S, Maman E, Rubinstein N, Avivi C, Orvieto R, et al. Localization of luteinizing hormone receptor protein in the human ovary. Mol Hum Reprod. 2014;20(9):844–9 PubMed PMID: 24874553. Epub 2014/05/31. eng.
Haas J, Ophir L, Barzilay E, Machtinger R, Yung Y, Orvieto R, et al. Standard human chorionic gonadotropin versus double trigger for final oocyte maturation results in different granulosa cells gene expressions: a pilot study. Fertil Steril. 2016;106(3):653–9 e1. PubMed PMID: 27341989. Epub 2016/06/28. eng.
Stock AE, Bouchard N, Brown K, Spicer AP, Underhill CB, Dore M, et al. Induction of hyaluronan synthase 2 by human chorionic gonadotropin in mural granulosa cells of equine preovulatory follicles. Endocrinology. 2002;143(11):4375–84 PubMed PMID: 12399434. Epub 2002/10/26. eng.
Kawashima I, Liu Z, Mullany LK, Mihara T, Richards JS, Shimada M. EGF-like factors induce expansion of the cumulus cell-oocyte complexes by activating calpain-mediated cell movement. Endocrinology. 2012;153(8):3949–59 PubMed PMID: 22673225. Pubmed Central PMCID: 3404342.
Shi H, Asher C, Yung Y, Kligman L, Reuveny E, Seger R, et al. Casein kinase 2 specifically binds to and phosphorylates the carboxy termini of ENaC subunits. Eur J Biochem. 2002;269(18):4551–8 PubMed PMID: 12230567. Epub 2002/09/17. eng.
Shi H, Asher C, Chigaev A, Yung Y, Reuveny E, Seger R, et al. Interactions of beta and gamma ENaC with Nedd4 can be facilitated by an ERK-mediated phosphorylation. J Biol Chem. 2002;277(16):13539–47 PubMed PMID: 11805112. Epub 2002/01/24. eng.
Yao Z, Flash I, Raviv Z, Yung Y, Asscher Y, Pleban S, et al. Non-regulated and stimulated mechanisms cooperate in the nuclear accumulation of MEK1. Oncogene. 2001;20(52):7588–96 PubMed PMID: 11753637. Epub 2001/12/26. eng.
Yung Y, Yao Z, Aebersold DM, Hanoch T, Seger R. Altered regulation of ERK1b by MEK1 and PTP-SL and modified Elk1 phosphorylation by ERK1b are caused by abrogation of the regulatory C-terminal sequence of ERKs. J Biol Chem. 2001;276(38):35280–9 PubMed PMID: 11463794. Epub 2001/07/21. eng.
Van Soom A, Tanghe S, De Pauw I, Maes D, de Kruif A. Function of the cumulus oophorus before and during mammalian fertilization. Reprod Domest Anim. 2002;37(3):144–51 PubMed PMID: 12071888. Epub 2002/06/20. eng.
Yokoo M, Sato E. Cumulus-oocyte complex interactions during oocyte maturation. Int Rev Cytol 2004;235:251–291. PubMed PMID: 15219785. Epub 2004/06/29. eng.
Hourvitz A, Lerner-Geva L, Elizur SE, Baum M, Levron J, David B, et al. Role of embryo quality in predicting early pregnancy loss following assisted reproductive technology. Reprod BioMed Online. 2006;13(4):504–9 PubMed PMID: 17007670. Epub 2006/09/30. eng.
Nishi Y, Yanase T, Mu Y, Oba K, Ichino I, Saito M, et al. Establishment and characterization of a steroidogenic human granulosa-like tumor cell line, KGN, that expresses functional follicle-stimulating hormone receptor. Endocrinology. 2001;142(1):437–45 PubMed PMID: 11145608.
This study was supported by grants from the Ministry of Health, Israel (3–00000-7410 to A.H.), and from the Chaim Sheba Medical Center (R&D Sheba grant 5/2014 to A.H.).
Availability of data and materials
The datasets during and/or analysed during the current study available from the corresponding author on reasonable request.
Ethics approval and consent to participate
The study was approved by the local Institutional Review Board (IRB) committee of Chaim Sheba Medical Center, Tel Hashomer (ethical approval number SMC-11-8707 and SMC-12-9342). Written informed consent was obtained from each patient who provided samples.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.