{'fraction_non_alphanumeric': 0.030540734251811738, 'fraction_numerical': 0.030301823684000955, 'mean_word_length': 5.009810959559704, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
MiR-181a Promotes Spermatogenesis by Targeting the S6K1 Pathway

Approximately 15% of couples suffer from infertility worldwide, and male factors contribute to about 30% of total sterility cases. However, there is little progress in treatments due to the obscured understanding of underlying mechanisms. Recently microRNAs have emerged as a key player in the process of spermatogenesis. Expression profiling of miR-181a was carried out in murine testes and spermatocyte culture system. In vitro cellular and biochemical assays were used to examine the effect of miR-181a and identify its target S6K1, as well as elucidate the function with chemical inhibitor of S6K1. Human testis samples analysis was employed to validate the findings. miR-181a level was upregulated during mouse spermatogenesis and knockdown of miR-181a attenuated the cell proliferation and G1/S arrest and increased the level of S6K1, which was identified as a downstream target of miR-181a. Overexpression of S6K1 also led to growth arrest of spermatocytes while inhibitor of S6K1 rescued the miR-181a knockdown-mediated cell proliferation defect. In human testis samples of azoospermia patients, low level of miR-181a was correlated with defects in the spermatogenic process. miR-181a is identified as a new regulator and high level of miR-181a contributes to spermatogenesis via targeting S6K1.


Introduction
Infertility besets about 15% of couples at childbearing age over the world, which equals about 48∼72 million couples according to the World Health Organization.
Female factors typically account for 50% of the cases while male factors are responsible for about 20∼30% of the cases. It is estimated that 2∼7% of all men exhibit different levels of sub-optimal sperm characteristics, such as alteration in sperm motility, morphology and sperm concentration (1,2). There is a deteriorating trend in the last few decades as continued decrease of sperm count has been reported from different regions of the world (3,4). The causes behind male infertility are complex, including the genetic mutation like Klinefelter syndrome with two X chromosomes and one Y chromosome, obstructive reproduction duct, immune disorder and inflammation, genital disorders, environmental exposure, and malfunction in the production of sperm (5,6). Azoospermia is one of the severe types of male infertility, it is a medical condition that no sperm is contained in semen due to various defects in the process of spermatogenesis or ductal system of testes (7,8). Of note, in the special case called ob-structive azoospermia, spermatogenesis is normal but no spermatozoa could be found in the ejaculate. Obstructive azoospermia is commonly seen as urologic disease for about 6.1∼13.6% of total male infertility patients (9,10). The treatment for male infertility is very limited and huge unmet needs exist (11)(12)(13)(14).
Germ cells undergo a fundamental and complicated process called spermatogenesis within the seminiferous tubule epithelium. Spermatogenesis includes one cycle of DNA replication and two consecutive cell divisions, during which dramatic changes occur in the cell cycle, growth, differentiation and morphology (15). Spermatogenesis is briefly divided into three successive phases: mitosis, meiosis and spermiogenesis. Spermatogonia first replicate the whole set of genes and divide into primary spermatocytes in mitosis like normal somatic cell division. Besides the self-renewal of spermatogonia, some offspring cells detach from the basement membrane and enter meiosis. Primary spermatocytes then undergo the second cell division without DNA replication and give rise to the haploid spermatids in the stage of meiosis. Finally, the spermatids further differentiate and undergo morphogenesis into spermatozoa in spermiogenesis. The spermatozoa are fully activated in the process called capacitation and are competent to fertilize the oocyte (16). During the process there is complex interplay between germ cells and supporting somatic cells, which requires cell and time-specific expression of a large set of genes and post-transcriptional regulation of genes and their products governing the diverse biological processes, including cell cycle and proliferation, DNA replication and chromatin remodeling, energy metabolism, mobility and morphology control (17). Many fundamental pathways are involved, such as MAPK (mitogen-activated protein kinase) (18), PI3K (phosphatidylinositol-4,5-bisphosphate 3-kinase) (19), TGF-β (transforming growth factor beta) (20) and mTOR (mammalian target of rapamycin) pathways (21). However, the detailed mechanism of gene regulation is still elusive.
Recently, microRNAs (miRNAs, miRs) have been implicated in the complex network of gene regulation in the process of spermatogenesis (22,23). MiRNAs are small non-coding RNAs of 19-23 nucleotides length which can bind to the target mRNAs through Watson-Crick base pairing and eventually cause the cleavage and degradation of mRNAs through the RISC (RNA-induced silencing complex) machinery (24). There have been multiple reports on miRNA's roles in spermatogenesis, such as miR-21, let-7, miR-17-92, miR-34c and miRNA-146 (25)(26)(27)(28). Some miRNAs are considered as potential biomarker for male infertility, for instance miR-141 is often upregulated in oligoasthenozoospermia (29). Genetic knockout studies suggest that loss of RISC-related components such as Dicer is linked with different defects of spermatogenesis (30), suggesting the importance of miRNA pathway in male germ cell development.
To shed more light on the role of miRNA in spermatogenesis, we studied the expression profile of miR-181a and established the connection between miR-181a and S6K1 in mouse spermatogenesis. miR-181a was highly expressed in testis tissues, especially the spermatogonia but not in spermatids. MiR-181a could downregulate S6K1 level and promote the cell proliferation of spermatogonia. Inhibition of miR-181a could suppress the cell proliferation which could be reversed by S6K1 inhibitor treatment. Interestingly, low level of miR-181a expression was observed in human male infertile patients with maturation arrest or hypospermatogenesis, compared to obstructive azoospermia patients which had normal spermatogenesis, suggesting the essential role of miR-181a to promote spermatogenesis.

Animals
C57BL/6 mice were from Shanghai SLAC Laboratory Animal Co. Ltd. Animals were housed in full ventilated cages in temperature controlled (20∼24℃) certified SPF grade facilities. All animals have free access to food and water unless indicated otherwise. All animal protocols were approved by the Institutional Animal Care and Use Committee of Qingdao Women and Children's Hospital, Qingdao University in accordance with the guidelines of the National Institutes of Health Guide for the Care and Use of Laboratory Animals. The testis tissues were col-lected from male mice at different stages after birth, various tissues from adult male mice were analysed for the expression of miR-181a, after the euthanasia of animals with CO2 inhalation and cervical dislocation.

RT-qPCR assay
Total RNA was extracted by RNeasy kit (Qiagen) and 25 ng total RNA was used for reverse transcription reaction. miRNA quantification was performed with TaqMan miRNA assay according to the manufacturer's instruction. Nucleolar small RNA RNU44 was used as internal control.

Western blotting
Total cell lysates were extracted with RIPA buffer and equal number of lysates were resolved on SDS-PAGE followed by standard immunoblotting procedures with indicated antibodies. The data were representative of three independent experiments.

MTT proliferation assay
GC-1 or GC-2 cells were plated at 2,000 cells/well density in 96 well plate and subject to different treatment/ transfection. After 48 h, cell medium was removed and after wash with PBS, fresh medium supplied with MTT reagent (0.3 mg/ml) was added to each well. The plate was scanned at plate reader (570 nm) for absorbance. Data were analysed with the Graphpad software.

Luciferase reporter assay
3'UTR of miR-181a (300∼500 nt) was inserted into pMirGlo vector, mutation at seed region was introduced as designed. 48h after transfection into GC-1 spg cells and various treatment, cell lysate was prepared and measured according to the instruction of manufacture (Promega dual luciferase kit). The results were analysed with Graphpad software.

Cell cycle profiling
Cells were treated with different conditions and then stained with propidium iodide (PI) solution after two brief washes with PBS and incubated at room temperature for 30 minutes with PI/RNase solution. Samples were analyzed on a FACScan flow cytometer (Becton Dickinson, San Jose, CA., USA). At least 20,000 events were collected on each sample. Cell cycle analysis of DNA histograms was performed with FlowJo software (Tree Star Inc., Ashland, OR., USA).

Human testicular samples analysis
Biopsy samples of 20 patients (aged 25∼40 years) diagnosed with Non-obstructive azoospermia (NOA) were obtained with approval of IRB committee in Qingdao Women and Children's Hospital, Qingdao University. Written informed consent were collected from patients before their enrolment into this study.

Statistics
All data were presented as means±S.D. Data was analysed by one-or two-way ANOVA analysis followed by a post hoc test. p＜0.05 was considered as statistically significant.

miR-181a was upregulated during spermatogenesis
To understand the role of miR-181a in mouse physiology, we profiled the expression level of miR-181a in different tissues of adult mice. As shown in Fig. 1A, the relative expression level of miR-181a of testis was the highest in adult male mice, followed by high levels in the kidney and thymus, compared to other tissues such as the brain, heart, spleen, muscle and skin. This indicated miR-181a was relevant to the function of testis. We next examined the time course of miR-181a expression in the different stages of postnatal mouse testes. From day 7 to 21, the miR-181a level showed a marginal increase from 100 to 130 with day 7 at relative level of 100. From day 28 to 42, the level increased dramatically to more than 300 (Fig.  1B). Maximum level was seen at more than 450 in the sample of day 42. This pattern matched the progress of spermatogenesis in mouse in which the mature spermatozoa are first produced in male mouse around 5 weeks after birth. The results demonstrated that miR-181a was enriched in male reproductive organ and upregulated in testis tissue during the spermatogenesis.

miR-181a promoted spermatogonia cell proliferation
To assess the role of miR-181a in spermatogenesis, we use the in vitro spermatogenic murine cell culture model  GC-1 and GC-2 cell lines. GC-1 spg cell resembles the transition between spermatogonia and primary spermatocytes and GC-2 spd (ts) cell is spermatocyte arrested at premeiotic stage. Firstly miR-181a expression level was determined in these cells. Compared to fibroblast NIH3T3 and myoblast C2C12 cell lines, GC-1 spg cells showed the highest level of miR-181a which was about 3-6 times higher than that in NIH3T3 and C2C12. MiR-181a level in GC-2 spd cells was about the 40% of that in GC-1 spg cells, as shown in Fig. 2A. Next, we knocked down the level of miR-181a by specific inhibitor in GC-1 spd cells to monitor the cell proliferation and cell cycle profile.
MTT cell proliferation assay result demonstrated the significantly reduced cell growth rate in miR-181a inhibitor transfected cells compared to scramble control (Fig. 2B). There was obvious G1-S arrest in miR-181a inhibitor transfected cells compared to control (Fig. 2C), as sig-nificantly less cells at S stage were detected in the cells transfected with miR-181a inhibitor. These results suggested that miR-181a was highly expressed in spermatogo- Luciferase reporter assay study of the interaction between miR-181a and 3'UTR of S6K1. GC1 spg cells were transfected with luciferase reporter genes conjugated with either wildtype (WT) or mutant (MUT) 3'UTR of S6K1. 24 hours later, the cells were transfected further with either miR-181a mimic or scramble control. Transfected cells were cultured for 48 hour and then harvested for the measurement of luciferase activities in a plate reader. (C) Western blotting analysis of mTOR1, S6K1, p-S6K1 in cells transfected with either miR-181a or scramble control. The representative image from at least three independent experiments was shown. The representative images from at least three independent experiments were shown. The data were presented as means±SD (*p＜0.05, **p＜0.01, ***p＜ 0.001 as compared with the control). nia stage and was required for cell proliferation.

S6K1 was downstream target of miR-181a in GC-1 cells
To find the downstream target of miR-181a, we use bioinformatic tool to identify potential match between the 3'UTR and seed region of miR-181a. One of the hits is S6K1 kinase which is an important protein kinase in mTOR signaling. mTOR is essential for cell survival, protein synthesis and energy metabolism. It has been reported to play a key role in spermatogenesis. To validate the in-silicon finding, we designed the mutant 3'UTR region of S6K1 and inserted into a luciferase reporter construct. Mutation at seed region on 3'UTR was introduced to disrupt the recognition of miR-181a as shown in Fig. 3A. Both of the luciferase assay reporters were transfected into GC-1 spg cells first, then the mimic of miR-181a or scramble control were transfected into cells. The luciferase assay result showed that miR-181a mimic could significantly suppress the reporter activity of wildtype S6K1 3'UTR. The suppression was abolished on mutant 3'UTR reporter construct, as shown in Fig. 3B. To further validate the result of luciferase assay, we examined the endogenous S6K1 protein level when miR-181a was transfected into GC-1 spg cells. Strikingly, S6K1 and phosphorylated S6K1 protein was completely undetectable in the miR-181a mimic transfected cells compared to scramble control (Fig. 3C), consistent with luciferase assay result. Collectively, these results confirmed that S6K1 was the downstream target of miR-181a.

S6K1 over-expression attenuated GC-1 spg cell proliferation
Next, we examine the function of S6K1 in GC-1 cell culture. In accordance with previous result, S6K1 protein level was much lower in GC-1 spg cells which had higher level of miR-181a, compared to GC-2 spd cells, as shown Fig. 4A. mTOR protein level was not changed between these two cell lines, suggesting the effect of miR-181a was very specific. To further assess the role of S6K1 in cell growth and differentiation, wildtype or kinase dead mutant form of S6K1 was over-expressed in GC-1 spg cells. As shown in MTT assay (Fig. 4B), cells overexpressed with wildtype S6K1 protein showed significantly reduced cell proliferation compared to cells with GFP control or S6K1 kinase dead mutant protein. Similarly, cells with wildtype S6K1 overexpression showed remarked G1/S arrest, compared to cells with GFP control or kinase dead S6K1 overexpression (Fig. 4C). This is in line with previous result in Fig. 2C, where inhibition of miR-181a showed G1/S arrest. These results demonstrated that S6K1 suppressed the spermatocyte proliferation in the manner dependent on its kinase activity.
Inhibiting S6K1 activity rescued the antagonizing effects of miR-181a inhibitor in GC-1 spg cell growth To further explore the relation between miR-181a and  inhibitor. In parallel, the miR-181a inhibitor transfected cells were also treated with S6K1 kinase inhibitor or DMSO. The growth rate of cells was monitored using MTT kit for 72 h. (C) Representative histogram data of the cell-cycle analysis of GC-1 spg cells in different treatment conditions as described above. The representative images from at least three independent experiments were shown. The data were presented as means±SD (*p＜0.05, **p＜0.01, ***p＜0.001 as compared with control).
S6K1, we introduced miR-181a inhibitor into GC-1 spg cells and observed the dramatic increase of S6K1 protein level compared to scramble control as shown in Fig. 5A, clearly suggesting that S6K1 was directly downregulated by miR-181a. We then examined whether S6K1 kinase inhibitor could affect miR-181a function. GC-1 spg cells Fig. 6. miR-181a was downregulated in the testis tissues of sterile men with non-obstructive azoospermia. miR-181a expression level in testis of patients with maturation arrest, hypo-spermatogenesis or non-obstructive azoospermia was quantified by real time PCR (n=10 for each group and **p＜0.01 as compared with non-obstructive azoospermia patient group). The expression levels of the mRNA were also quantified in normal testis tissues from post-mortal donors.
were transfected with miR-181a inhibitor or scramble control first, then S6K1 inhibitor or DMSO control vehicle was added to the cells with miR-181a inhibitor transfection. As shown in Fig. 5B, cells with miR-181a inhibitor showed decreased proliferation compared to cells with scramble control, this suppressive effect on cell growth was reverted by S6K1 inhibitor treatment. Likewise, miR-181a inhibitor could cause significant G1/S arrest as shown by cell cycle profiling in Fig. 5C, and S6K1 inhibitor treatment could recover the cell population in S phase to the level of scramble control. These results indicated that S6K1 was the downstream factor of miR-181a to regulate cell growth and inhibitor of S6K1 was effective to revert the suppressive role of low level of miR-181a.
miR-181a was downregulated in the testis tissues of sterile men with non-obstructive azoospermia To extend our finding in mouse cellular model of spermatogenesis, we assessed the level of miR-181a in human testis samples by RT-qPCR method. In non-obstructive azoospermia such as maturation arrest or hypo-spermatogenesis type of male infertility patients, the level of miR-181a was found at relatively low level, as shown in Fig. 6. However, miR-181a level was significantly higher in the patients of obstructive azoospermia in which the spermatogenesis is normal. In addition, we also included the expression levels of miRNA in normal testis tissues as comparison and observed a slightly lower expression than the sterile patients from maturation arrest or hypo-spermatogenesis types. This result indicated that high level of miR-181a was well correlated with normal spermatogenesis, in accordance with the previous results from in vitro cell model. Future study will be performed to establish the role of S6K1 in these patients and test the therapeutic potential of S6K1 kinase inhibitor in the non-obstructive azoospermia patients with low level of miR-181a.

Discussion
Male infertility is a common cause for couples who have difficulty to bear a child. Azoospermia is a severe type of male infertility with elusive mechanisms and limited available therapy. MiRNAs have been implicated to play a role in the spermatogenesis. Here we reported a new mechanism of miR-181a in this process. MiR-181a is highly enriched in male reproductive organ and increased expression level was observed in the maturation peak of spermatogenesis. Using a murine spermatocyte culture system, we demonstrate that miR-181a is required for cell proliferation and cell cycle progression. Mechanistically miR-181a could target S6K1 protein kinase of mTOR pathway by 3'UTR reporter assay and protein immunoblotting. Further we found S6K1 overexpression in spermatocyte could suppress the cell growth and kinase activity is essential for the suppression. On the other hand, S6K1 kinase inhibitor could revert the cell growth arrest by knocking down miR-181a level. In other words, S6K1 kinase inhibitor could compensate the effect caused by low level of miR-181a. Lastly human testis sample from male infertility patients revealed the correlation of high level of miR-181a with functional spermatogenesis. Our study established the role of miR-181a-S6K1 axis in spermatogenesis and pointed out the potential therapeutic application of miR-181a mimic and S6K1 protein kinase inhibitor in non-obstructive azoospermia.
Spermatogenesis involves the highly orchestrated waves of transcriptional and translational regulation of gene networks in the different cell types like gametes and somatic cells. Notably, only a small fraction of human genome is transcribed into coding mRNAs. Recently non-coding RNAs have attracted increased attention and play important functions in multiple biological processes. Among them, microRNAs are small non-coding RNAs which constitute the fine-tune control of huge number of coding genes as a supplement to conventional gene regulation. Many miRNAs have been reported to play important role in spermatogenesis, for example, miR-15a (31), miR-18 (32), among others as reviewed (22,33). A microarray profiling of miRNAs in mouse testes has mentioned the potential role of miR-181 family members (34). However, our study provides the thorough mechanistic insight of miR-181a in spermatogenesis and identifies the downstream target S6K1 which is the downstream key a serine/threonine kinase activated by mTOR. S6K1 is considered a master regulator of cellular metabolism and survival together with mTOR and its downstream effector ribosomal protein S6 (35)(36)(37). mTOR has been reported to be important for Sertoli cell polarity and spermatogenesis in human (38). A recent study of mTOR and S6K1 in rat spermatogenesis reported some seemingly contradictory findings with our results (21). Their conclusion that mTOR/ p70S6K1 promotes spermatogonia was based on the observation from 8-week, 38-week and 80-week old rats. The effect of decrease in mTOR signaling is marginal in their results and the significant change at 80-week rats is more about aging instead of spermatogenesis which should be active at much earlier stage of life. In addition, our result was supported by another miRNA study in mouse spermatogenesis that activation of mTOR pathway is detrimental to cell survival (39). The study using knockout miR-17-92 offered the solid proof on its role in spermatogenesis. What's more, gene disruption of S6K1 in mice does not cause male infertility (35,36), supporting our finding with S6K1 inhibitor (Fig. 5B and 5C). We think the discrepancy between our results and Xu et al. (21) report may be originated from different species used and different stages of sperm development under observation. The method difference is also noticeable. Xu et al. (21) observed that rapamycin treatment in 8-week male mice led to decreased sperm count which was consistent with human study with rapamycin. However, whole body administration of mTOR inhibitor rapamycin may have profound effects. It is reported that mTOR is essential for the polarity maintenance of the Sertoli nursing cells by controlling the tight gap junction in addition to regulation of cell metabolism (40). Our finding is based on the spermatocyte itself rather than focusing on the Sertoli cell and germ cell interaction. Loss of all members of miR-181 clusters leads to decreased survival, retarded growth and lymphocyte homeostasis in mice (41). This report is largely in line with our finding on cell proliferation, although there might be different functions for miR-181 in various type of tissues. Therefore, further tissue-specific knockout of miR-181a subtype in male reproductive organ is required to elucidate the mechanism of miR-181a in spermatogenesis in a more physiologically relevant condition.
Activation of mTOR/S6K1 protein kinase pathway is frequently observed in multiple cancers and thus targeting mTOR in cancer is under development (42). Inhibitor of S6K1 protein kinase has been evaluated in many different cancers (43)(44)(45). Inhibition of S6K1 caused the cell growth arrest and apoptosis in the setting of cancer cell lines. However, in our study we found inhibitor of S6K1 actually promoted the spermatocyte proliferation in the background of miR-181a inhibition. We speculate there might be cell type specific effect and miR-181a should have other targets besides S6K1 gene. The effect we seen was the compound consequence of inhibiting both miR-181a and S6K1. Further study is needed to dissect the detailed regulatory loop between miR-181a and S6K1, and other potential genes involved in the unique setting of spermatogenesis.
In summary, we reported the function of miR-181a in mouse germ cell differentiation and proliferation by the profiling of miR-181a expression in different tissues and development stages of male spermatogenesis and further unraveling the novel connections between miR-181a and S6K1 protein kinase as the mechanism. The arrest of cell growth caused by knockdown of miR-181a could be rescued by S6K1 kinase inhibition. Interestingly low level of miR-181a was found in sterile male testis samples compared to the ones with normal spermatogenesis. This finding opens a new avenue for the future study to investigate the inhibition of S6K1 in non-obstructive azoospermia patients with low level of miR-181a.