miR-29a/b1 Regulates the Luteinizing Hormone Secretion and Affects Mouse Ovulation

miR-29a/b1 was reportedly involved in the regulation of the reproductive function in female mice, but the underlying molecular mechanisms are not clear. In this study, female mice lacking miR-29a/b1 showed a delay in vaginal opening, irregular estrous cycles, ovulation disorder and subfertility. The level of luteinizing hormone (LH) was significantly lower in plasma but higher in pituitary of mutant mice. However, egg development was normal in mutant mice and the ovulation disorder could be rescued by the superovulation treatment. These results suggested that the LH secretion was impaired in mutant mice. Further studies showed that deficiency of miR-29a/b1 in mice resulted in an abnormal expression of a number of proteins involved in vesicular transport and exocytosis in the pituitary, indicating the mutant mice had insufficient LH secretion. However, the detailed mechanism needs more research.


INTRODUCTION
The miR-29 family consists of three related mature miRNAs, miR-29a, miR-29b and miR-29c, which are processed from two precursor sequences located at two distinct genomic clusters of miR-29a/b 1 and miR-29b 2 /c. Members of the miR-29 family are ubiquitously expressed, have considerable overall sequence homology with the same seed sequence. Although they have similar tissue expression patterns, miR-29a is the dominant member accounting for more than 50% of total miR-29 expressed in all tissues (1). miR-29 play important roles in regulating a number of physiological and pathological processes, including metabolism (1-3), inflammation (4,5), fibrosis (6), cancer (7) and neurodegeneration (8).
As a potential clinical marker or new form of nucleic acid drug, much attention has been paid to miR-29 research (9, 10). miR-29 deficiency causes a wide range of physiological defects in mice. Premature cardiac fibrosis and atherosclerotic plaque remodeling is considered as a result of abnormal expression of miR-29 target genes Col4a (11) and ECM (Col1a and Col5a) (12), and heart failure and metabolic disorders might be caused by up-regulating the target gene PCG1a (1). miR-29a responsible for repressing LPL in hepatocytes, contributes to physiological lipid distribution and protects hepatocytes from steatosis (13). Homozygous deletion of miR-29a/b 1 in mice led to decreased self-renewal and increased apoptosis in hematopoietic stem cells (HSCs) through up-regulating Dnmt3a (14). In addition, early puberty in hypothalamic miR-29 knockdown females is attributed to ectopic expression of Tbx21, a target gene of miR-29 (15). Reproduction in miR-29 brain-specifical knockdown mice was affected in a sex-dependent manner, with female mice exhibiting hyperfertility and males being subfertility (16); however, this result is inconsistent with the sterile phenotype reported in the miR-29a/b 1 knockout mice (1). Therefore, the relationship between miR-29a/b 1 and reproductive function is still not well understood.
In this work, we revealed that female miR-29a/b 1 knockout mice exhibit severe fertility problems. We proposed that the lack of miR-29a/b 1 in female mice may interfere with the secretion of luteinizing hormone in the pituitary, leading to ovulation failure and a subfertile phenotype.

Generation of miR-29a/b 1 Knockout Mice
A miR-29a/b 1 knockout mouse line was established using CRISPR/Cas9 gene editing technology and was supplied by Shanghai Center for Model Organisms (SMOC) (17). miR-29a/ b 1 −/− homozygous animals and their wild-type littermates were obtained by mating corresponding heterozygotes with each other. Genomic DNA was extracted from tail biopsies, using magnetic bead DNA isolation Kit (DE0596D, EmerTher, Shanghai). PCR was adopted for genotyping using 2 × Taq Plus Master Mix (P212-01, Vazyme) under the following conditions: denaturation at 98°C for 2 minutes, then 35 cycles of 98°C for 10 seconds, annealing at 63°C for 15 seconds, and extension at 68°C for 60 seconds. Primers used for genotyping are listed in Table S1.

Animals
All animals were housed in a specific pathogen-free environment (12 h light/12 h dark with lights on at 7.00 h at 21 ± 2°C) with food and water ad libitum. This study was performed in strict accordance with institutional guidelines and approved by the Institutional Animal Care and Use Committee of Shanghai Model Organisms, and the IACUC permit number is 20090002.
Fertility Assessment 8-week-old miR-29a/b 1 KO and wild-type virgin female or male mice were bred with wild-type male or female mice with known fertility at a proportion of ♀2: ♂1, and vaginal plug formation was examined every morning for 20 consecutive days. Pregnant female mice were separated and pups were recorded, while nonpregnant mice continued to mate. Female or male mice that did not conceive within 1 month of mating were defined as infertile.

Sexual Maturity and Vaginal Smear
Female miR-29a/b 1 KO and wild-type mice from the age of 3week-old were examined twice daily with respect to on vaginal opening as a marker of rodent sexual maturity. The date of vaginal opening in each mouse was recorded. Female miR-29a/b 1 KO and wild-type (8-10 weeks old for each genotype) mice were caged individually for 3 weeks and at least two full estrous cycles were obtained in each mouse.
Vaginal smears were collected daily and the determination of estrous cycle was evaluated microscopically with the vaginal epithelium. The vaginal epithelium obtained from the vaginal opening by gently eluting 10 ml of physiological saline solution 2-4 times, then the vaginal epithelium transferred onto a microscopic slide and dried at room temperature and fixed with 100% methanol. The slides were stained with Wright's Giemsa (BASO) stain and examined with light microscopy. Proestrus cells are well-formed nucleated epithelial cells. Animals with 85% superficial epithelial cells were considered to be estrus. During metestrus, cornified squamous epithelial cells often in fragments, as well as leukocytes, may be observed. Otherwise, the predominant presence of leukocytes in the cytological smear was identified as diestrus.

Ovariectomy
Adult (8-10 weeks) miR-29a/b 1 KO and control females in diestrus morning were injected subcutaneously with pentobarbital (effective dose 320 mg/kg). Mice were deeply anesthetized and placed on a heating pad. The back skin was shaved and cleaned. About 1.0 cm long incision was made through the muscle layer above the ovaries on each side of the midline. Through the incision, the ovaries were gently pulled outside the body and removed by cauterization below the oviduct. The skin incision was closed with sutures. The mice were left to recover on a heating pad. Adult sham-operated mice were in diestrus on the day of recording as determined by vaginal cytology. Sham-treated animals were processed in the same way, except for the intact ovaries retained. Mice were killed 7 days post-surgery, and their serum were measured for LH and FSH levels.

Hormone Measurement
For hormone measurement, orbital blood was collected in the morning (10.00 h-11.00 h) and evening (18.00 h-19.00h) (18) from freely-moving conscious animals during randomly estrous cycle stages, and were kept at room temperature for 30 minutes. Serum was obtained by centrifuging for 15 minutes at 3000 g at 4°C and was stored at -80°C until analysis. Serum levels of hormone and pituitary proteins LH level were analyzed by Shanghai WESTANG BIO-TECH cooperation using enzymelinked immunosorbent assay (ELISA). The minimum detectable level of the LH assay was 0.1mIU/ml and the intra-and interassay coefficients of variation were 9.9% and 8.3%. The minimum detectable level of the FSH assay was 1mIU/ml and the intra-and inter-assay coefficients of variation were 9.8% and 8.6%, The minimum detectable level of the estrogen assay was 30 pg/ml and the intra-and inter-assay coefficients of variation were 9.3% and 8.5%, The minimum detectable level of the progesterone assay was 0.2ng/ml and the intra-and inter-assay coefficients of variation were 9.5% and 8.3%, The minimum detectable level of the testosterone assay was 0.1ng/ml and the intra-and interassay coefficients of variation were 9.4% and 8.2%, respectively.

GnRH Challenge
Animals received an intraperitoneally injection with 125ng/g (19) exogenous GnRH (L7134, Sigma-Aldrich, St Louis, MO, USA) or saline vehicle. Twenty minutes after GnRH or saline injection, orbital blood was collected, and the resultant serum samples were stored at −80°C for subsequent human-LH radioimmunoassay (RIA, performed by Beijing North Institute Biological Technology, Beijing, China) (20), with sensitivity and intra-and inter-assay coefficient of variation for LH of 0.5 mIU/ ml, 15% and 20%, respectively.

Superovulation and Oocyte Collection
Superovulation: To induce superovulation, 8-week-old mice were intraperitoneally injected with 5 IU pregnant mare serum gonadotropin (PMSG, Sigma) at afternoon (15:00h-16:00 h), followed by 5 IU human chorionic gonadotropin (hCG, Sigma) 48 hour later to trigger oocyte maturation and ovulation. Female mice were mated with 10-week-old fertile wild-type males 16 h after injection and checked for vaginal plug formation the next morning.
Oocyte collection: Super ovulated or natural mated mice with a visible plug were sacrificed by cervical dislocation, the ovaries were removed and the ampulla was collected. Oocytes were harvested in M2 media and quantified by microscopy (Nikon SMZ800) following brief digestion in hyaluronidase (800IU/ml, Sigma) to strip cumulus and pipetting for 30-60 s. Oocytes were washed 5 times with PBS. The washed oocytes were transferred to M16 media and cultured overnight, and two-cell stage embryos were counted in the next morning.

Histological Analysis and Follicles Count
Wild-type and miR-29a/b 1 KO were euthanized and transcardially perfused with cold saline, followed by 4% paraformaldehyde 0.1 M phosphate buffer (PFA). Brain, pituitary and ovary were collected and fixed overnight at 4°C. Paraffin-embedded ovary samples were serially sectioned at 4 mm-thick sections. Brain coronal (20 mm) slices were cut with a Leica CM1950 following by dehydration in 30% sucrose saline solution.
Pituitary and ovary stained with hematoxylin and eosin using standard histological techniques (Servicebio). Stained sections were scanned using LEICA CTR6000 with a 10X, 20X and 40X objective. Ovarian follicles at different developmental stages were classified and quantified in serial sections according to the Pedersen and Peters method (22). To avoid double counting of follicles across sections, only follicles containing oocyte with a clearly visible nucleus were scored (23), and follicles were counted in every fifth serial section. Any follicle also appearing in the adjacent lookup section was not counted. The entire section was analyzed without subsampling. Each ovary was coded with no information about genotype group for blind counters and prevent bias. The mean count per section was calculated. All follicle types were summed together to determine the total number of follicles.
For immunofluorescence, brain and pituitary sections permeabilized by incubation with 0.1% Triton X-100 in PBS for 10 minutes at room temperature. After permeabilization, the sections were washed three times in PBST, and blocked with 5% normal donkey serum in PBS for 1h at room temperature, then were incubated with mouse anti-Lutropin beta antibody (1:1000, SANTA CRUZ, sc-373941), rabbit anti-GnRH1 (1:500, Immunostar, PA1-121) or rabbit anti-GnRHR antibody (26) (1:100, Proteintech, 19950-1-AP) overnight at 4°C, followed by staining with Alexa Fluor 647-conjugated donkey anti-mouse antibody (Invitrogen Molecular Probes) or Alexa Fluor 594conjugated donkey anti-rabbit (Invitrogen Molecular Probes) antibody and DAPI dye to stain nuclei. The mouse liver and lung tissues were selected to negative control for GnRHR and Lutropin beta antibody respectively.
Stained sections were scanned using the 40X objective of a Zeiss Confocal microscope (LSM880). The area fractions of positive cells relative to entire area were determined using ImageJ (Fiji, NIH) software. Cell location was mapped to the atlas (27).

Real-Time Quantitative PCR
Total RNA was isolated using TRIzol (Tiangen Biotech, Beijing) according to the manufacturer's instructions and kept at -80°C subsequent for use. For microRNA measurement, 2 mg total RNA was transcribed into cDNA using the miRcute Plus miRNA First-Strand cDNA Synthesis Kit (Tiangen Biotech, Beijing). Expression level of mature miR-29a, miR-29b and miR-29c were measured using miRcute Plus miRNA qPCR Detection (Tiangen Biotech, Beijing). U6 snRNA was used for normalization.
For mRNA measurement, total RNA (2 mg) from each sample was transcribed by using EasyScript First-Strand cDNA Synthesis SuperMix (TransGen Biotech, Beijing), and mRNA levels of target genes were detected using TransStart Tip Green qPCR SuperMix (TransGen Biotech, Beijing) according to the manufacturer's instructions. Murine b-actin was used as a reference to normalize target gene expression levels. Real-time PCR amplification was performed using the Realplex system (Applied Biosystems QuantStudio3, ThermoFisher Scientific). The sequences of the specific primers used are listed in Supplementary Material, Table S1. RNA levels were calculated using the 2 −DCT method, where CT is the cycle threshold (28). Melting curve analysis for each primer set revealed only one peak for each product, and the sizes of PCR products were confirmed by comparing sizes with a commercial ladder after agarose gel electrophoresis. PCR products were further confirmed by sequencing.

Western Blot
Mice were euthanized and tissues were collected. Total tissue protein was extracted using RIPA buffer (ThermoFisher scientific) containing protease and phosphatase inhibitor cocktails (Selleck Chemicals). Protein concentration was quantified using the Enhanced BCA Protein Assay Kit (Beyotime). Protein (20 mg) from each sample was separated on 4%-20% SDS-PAGE (GenScript) and transferred onto nitrocellulose membranes (GE Healthcare). Membranes were blocked with Western BLoT Blocking Buffer (Protein Free) (Takara) for 1 h at room temperature and then incubated with primary antibodies, Lutropin beta (1:1000, SANTA CRUZ, sc-373941) or anti-b-actin (1:1000, Santa Cruz, sc-47778) diluted in Western BLoT Immuno Booster PF (Takara) at 4°C overnight. After washing with TBST three times, membranes were incubated with fluorescent-conjugated secondary antibody for 1 h (1:10000, LI-COR Biosciences). Quantitative detection of protein expression was then performed using the Odyssey Infrared Imaging system (LI-COR Biosciences) and analyzed with Image J software (National Institutes of Health, Bethesda, MD, USA).

Statistical Analysis
Data analysis was performed using GraphPad Prism 7 (GraphPad software Inc.). Data are expressed as the mean ± SEM. Difference in mean values between two groups were analyzed using the Student's t-test (continuous variables) or Mann-Whitney test (discrete variables). For comparisons involving more than two groups, ANOVA (continuous variables) or Kruskal-Wallis (discrete variables) with post hoc testing was used, and survival profiles were constructed by Kaplan-Meyer survival analysis. Statistically significant differences are shown with asterisks (*p < 0.05, **p < 0.01, ***p < 0.001, and **** p < 0.0001).

RESULTS
Genetic Ablation of miR-29a/b 1 Leads to Female Sterility A conventional miR-29a/b 1 knockout mouse line (miR-29a/b 1 KO) was previously established using CRISPR/Cas9 methods (17). The genotyping and the expression of miR-29 in different genotypes of mice were detected by PCR and real-time PCR, respectively ( Figure S1). To understand the role of miR-29a/b 1 in fertility, the reproductive ability of miR-29a/b 1 KO mice was evaluated. For data in Table 1, of the 25 females tested, 23 were sterile. The two pregnant miR-29a/b 1 KO female mice gave birth to two offspring each and were not subsequently pregnant again. Among males, 66.7% miR-29a/b 1 KO were still fertile. However, female miR-29a/b 1 KO mice exhibited serious reproductive problems. Vaginal plugs were checked to study the mating behavior. miR-29a/b 1 KO females had a significant lower mating frequency compared to wild-type females ( Table 1), suggesting abnormal sexual maturity and estrous cycle. Sexual maturity indicated by vaginal opening occurred 5 days later in miR-29a/b 1 KO female mice (postnatal day 28) compared to wild-type littermates (postnatal day 23) ( Figure 1A). At the time of puberty onset, the mutant mice are significantly lighter than wild-type mice ( Figure 1B). Meanwhile, abnormal estrous cycle with less time in estrus and metestrus and significantly more time in diestrus was observed in miR-29a/b 1 KO female mice. (Figures 1B, C and Figure S2). RT-PCR analysis revealed that expression of miR-29a periodically changed in pituitary and ovarian tissues ( Figure S3A), suggesting that miR-29a/b 1 may play a role in the estrous period in mammals. Taken together, these data illustrate that loss of miR-29a/b 1 induces growth retardation in mutant mice and subfertility in females.

miR-29a/b 1 Gene Knockout Leads to Decreased Plasma LH Level and Ovulation Disorder
Ovary and uteri weight in miR-29a/b 1 KO females were significantly reduced compared to wild-type females (Figures 2A, B), whereas, in males, testis and seminal pouch in miR-29a/b 1 KO mice and wild-type counterparts showed no difference ( Figure S4). Fertilized eggs were collected from the oviducts of wild-type and miR-29a/b 1 KO females with vaginal plug after mating with wild-type males. In 20 females miR-29a/b 1 KO mice, only 2 oocytes were found and with no two-cell embryos the next day, while among five wild-type mice, 34 oocytes and 19 two-cell embryos were collected ( Figure 3A). Histomorphometric analysis revealed that mutant ovaries contained normal primordial follicles, a similar number of secondary follicles with normal oocyte and a thick granulosa cell layer, indicating that the early follicles developed normally, but lacked corpora lutea formation ( Figures 3B-D). These results suggest that subfertility of the mutant female mice may be caused by an ovulation disorder. In females, hormonal control of the estrous cycle and ovulation is essential for the establishment of maturation and fertility in mammals (29). Thus, we examined hormone levels in the serum of wild-type and miR-29a/b 1 KO female mice. In the female miR-29a/b 1 KO mice, significant decreases in the serum LH and progesterone (P 4 ) ( Figure 4A) were observed, while there was no apparent difference in serum content of folliclestimulating hormone (FSH) or Testosterone (T) or Estradiol (E 2 ) compared to wild-type mice (Figures S5A-C). Cyp19a 1 and Cyp17a 1 , encoding enzymes involved in estradiol and testosterone synthesis, were expressed at identical levels in ovaries from the two groups of mice, while the Cyp11a, which essential to the level of sex hormones, was significantly decreased in ovaries from mutant mice ( Figure S5D). These results indicated that impaired corpora lutea formation in miR-29a/b 1 KO mice might be caused by a shortage of LH. This speculation was further confirmed by the superovulation experiment. Ovulation in the mutant mice was rescued by exogenous gonadotropin injection, indicating that responses to LH stimulation were not irreversibly lost in these mutant animals ( Figure 4B). Ovaries from superovulated adult miR-29a/b 1 KO mice showed normal morphology, and the corpora lutea were formed ( Figure 4C).
To determine whether the central regulated mechanisms mediating ovulation were altered in miR-29a/b 1 KO mice, females were subsequently treated with an intraperitoneally injection of 125ng/g GnRH or saline vehicle at 10.00 AM. Normal GnRH responsiveness was observed in miR-29a/b 1 KO pituitary, but serum LH level in miR-29a/b 1 KO females remained markedly below the levels observed in wild-type littermates ( Figure 4D). Furthermore, The GnRHRimmunoreactivity in the pituitary of miR-29a/b 1 KO mice was increased compared to wild-type mice ( Figures 5A, B), Again, to assess the impact of hyperstimulation with endogenous GnRH modulated by estrogen (30)(31)(32)(33), female control and miR-29a/b 1 KO animals were castrated or underwent a sham surgery.  Animals were euthanized after 7 days, and serum concentrations of LH and FSH were measured. Consistent with control females, castration resulted in an increase in both LH and FSH compared with sham-operated controls, however, the post-castration rise in LH secretion was blocked in miR-29a/b 1 KO females, while the FSH level was no significant differences in mutant mice serum from controls ( Figure 5C). LH levels overall were markable lower in miR-29a/b 1 KO females relative to controls. There was no apparent difference in Kiss1and Gnrh1, which stimulating secretion of gonadotropin releasing hormone from the hypothalamus (34-37) and luteinizing hormone from the pituitary (35), respectively ( Figures 5D, E). These results suggest that ovulation disorder in miR-29a/b 1 KO mice might be caused by dysregulation of related pituitary hormones, especially LH.

Dysregulated Pituitary LHb Release in miR-29a/b 1 KO Mice
LH is synthesized in and secreted by the pituitary. A lack of miR-29a/b 1 was confirmed in mutant pituitary tissues ( Figure S3B). The anterior pituitary undergoes rapid proliferation in neonatal mice, subsequently expanding the cells that produce factors required for growth and reproduction (38). Defective anterior pituitary development in animals contributes to many organismlevel developmental defects (39). However, there was no difference in pituitary structure, size or position of the anterior pituitary between wild-type and miR-29a/b 1 KO mice ( Figure 6A). No abnormalities were found upon pathological examination of mutant pituitary tissues ( Figure 6B). Notably, transcript levels of the Lhb gene in miR-29a/b 1 KO pituitary did not differ from control animals, but LH protein level and immunoreactivity were even higher in KO mice ( Figures 6C-F).
To further elucidate the effects of miR-29a/b 1 gene knockout on pituitary function, iTRAQ analysis was performed to compare proteomic changes in the pituitary between mutant and wildtype mice. Total pituitary protein from three biological replicates of each genotype were subjected to LC-MS/MS analysis. The hierarchical clustering profile of differential proteins is shown in the heat map ( Figure 7A). A total of 163 cellular proteins were statistically significant altered (p<0.05), including 75 upregulated proteins and 88 downregulated proteins ( Figure 7B, Table 2). LHb and FSHb were significantly increased in the pituitary of miR-29a/b 1 KO mice according to b/y ion signal intensity ( Figure 7C). Besides, TSHb and Cga were also markedly upregulated as a result of the miR-29a/b 1 deficiency ( Figures 7C, G). Through GO analysis, altered proteins identified in this study were found to be involved in a wide range of biological process, and most of the differential proteins were classified in the protein transport processes, which are essential for vesicle-mediated transport in the cytoplasm and exocytosis during plasma infusion (40-42) ( Figure 7D).
The intersected gene between upregulated expression and miR-29a targets through miRDB (http://mirdb.org) were analyzed, 11 potential direct target transcripts of miR-29a were discovered ( Figure 7E), and predicted target genes were expected to be upregulated in miRNA loss of-function models ( Figure 7G). Among them, collagen family Col1a1, Col4a2 and Col5a1 are target genes of miR-29a-3p, and promote cancer cells invasion and migration (43)(44)(45). In addition, miR-29a can promote the neurite outgrowth by targeting extracellular matrix-related genes like Fibrillin 1 (Fbn1) and hyaluronan and proteoglycan link protein 1 (Hapln1) (46,47), which dramatically increased in the pituitary of miR-29a/b 1 KO mice. Hdac4 (48), which is key epigenetic modified writer, may play important roles in the change of gene expression pattern in miR-29a/b 1 gene knockout mice, especially for down-regulated genes. For the 88 down-regulated proteins in the pituitary of miR-29a/b 1 KO mice, a considerable portion of them participate in vesicle-mediated transport and secretion (Ergic1, Fkbp2, Ssr3, Stat5a, Crhbp, Figure 7F). Notably, Ergic1, encodes a cycling membrane protein, and plays an important role in transport between endoplasmic reticulum and Golgi (49). Absence of trApg (SSr3) impairs protein translocation into the endoplasmic reticulum and affects transport (50). Myosins were reported as core players in the final stages of regulated secretory pathways (51). Treatment of pituitary cells with the myosin light chain (Myl2/3) kinase inhibitor, wortmannin, attenuated GnRH-induced LH release (52). Further validated by quantitative PCR (qPCR) that the mRNA transcripts of these genes, which were consistent with LC-MS/MS ( Figure 7G). These results indicated the deficiency of miR-29a/b 1 blocked proteins transportation, leading to impaired pituitary hormone secretion, especially LH released.

DISCUSSION
A lack of miR-29a/b 1 leads to female sterility in mice, which has been mentioned previously (1); however, the mechanisms underlying this result were not published or illustrated. In this work, we demonstrated that low serum LH level and ovulation disorder might be the direct cause of subfertility in female miR-29a/b 1 KO mice. This conclusion is further proved by the results that oocyte development is normal in the ovaries of mutant mice and normal eggs could be obtained through super-ovulated. Compared to wild-type mice, the pituitary gland in mutant mice stimulated with the same concentration of GnRH produced reduction LH secreted into the blood, indicating that   (53). Meanwhile the expression of LH protein was higher in mutant pituitaries than that in wild types. This suggests that knockout of miR-29a/b 1 results in deficits in LH secretion from the pituitary but not in LH synthesis stimulated by GnRH (54,55). Proteomic analysis of the pituitary showed that a large number of proteins related to cellular vesicle-mediated secretion and protein transport were significantly changed in miR-29a/b 1 KO mice. This effect seems to be omnidirectional, from the vesicle  transport between endoplasmic reticulum and Golgi apparatus, as well as the process of docking and priming of secretory vesicle on the cell membrane. As a result, many kinds of secretory proteins, including LHb, were accumulated in pituitary cells. These secreted proteins accounted for 44% of the upregulated proteins in the pituitary of mutant mice.
It is worth noting that FSH required for follicle growth and development and maturation of the ovum (56-58) was less affected by the knockout of miR-29a/b 1 . Different secretion modes between FSH and LH might be an important reason (59,60). LH is secreted via a regulated pathway, while FSH release is primarily constitutive and controlled by synthesis. Increased FSH protein level in the       Figure 7C), which may also explain the fertility of male mutant mice. There is much agreement that FSH influences the mitotic activity of the spermatogonia and promote cellular differentiation during the pubertal phase (61). Testosterone regulated by LH also plays a role for spermatogenesis, however, completely T-independent spermatogenesis is possible if high-dose FSH treatment (62). Of note, the use of intraventricular injection of miR-29 inhibitor or overexpression of an antisense sequence targeting miR-29 in the brain to knockdown expression of miR-29 leads to earlier puberty onset or hyperfertility (63). These findings are not consistent with our results. It is possible that lack of miR-29a/b 1 function throughout development could result in compensatory effects which may lead to differences between our results and the results of the above literature. The underlying reasons for the different effects between knockout and knockdown need to be further studied. In addition, it should be noted that KO mice also showed growth retardation (64). We found that the weight of KO mice remained light, even though they had reached sexual maturity. So, the causal relation between two events cannot be confirmed now. We speculated that growth retardation and delayed maturity may come from the same reason, which happened in pituitary or upstream signal of KO mice.
In conclusion, LH secretion was impaired by miR-29a/b 1 knockout which caused ovulation deficiency in the mutant mice. Further studies revealed the effect of miR-29a/b 1 on hormone secretion function in the pituitary. Our work provides novel mechanistic insights into the relationship of miR-29a/b 1 and reproduction, opening the possibility of clinical approaches to reproductive studies based on the regulatory circuitry of miR-29a/b 1 .

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

ETHICS STATEMENT
The animal study was reviewed and approved by Institutional Animal Care and Use Committee of Shanghai Engineering Research Center for Model Organisms, SMOC. Written informed consent was obtained from the owners for the participation of their animals in this study.

AUTHOR CONTRIBUTIONS
YG, RSu, RSh and JF designed research. YG and JF analyzed data. YG, YW, HS, HZ, LC, QH, ZhiW, and YT performed research. YG, LX and JF wrote the paper. HY, MZ and ZhuW contributed to discussion and the proof reading of the paper. All authors contributed to the article and approved the submitted version. Supplementary Figure 1 | Genotyping of miR-29a/b 1 KO knockout mice.