Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 14, Article number: 7668 (2024 ) Cite this article Imida Pesticide
The phenylpyrazole insecticide fipronil has wide-ranging applications from agriculture to public health to control undesirable organisms. However, several studies have reported the residual environmental hazards of fipronil and demonstrated its harmful effects even in mammalian reproduction. Therefore, this study was conducted to demonstrate the mode of action of fipronil on mouse spermatozoa. We treated fipronil to spermatozoa and performed comprehensive function evaluations. Moreover, proteomic analyses were conducted to identify the alteration of protein expression levels in spermatozoa. Most of sperm motility and kinematic parameters and intracellular ATP levels were diminished, and the spontaneous acrosome reaction was promoted after treatment with fipronil. Proteomic analyses revealed altered expression levels of 14 proteins after treatment. These proteins have been reported to be associated with sperm-specific pathways, prominently the cytoskeleton of the sperm, “9 + 2” axoneme composition, metabolism, and fertility. Collectively, our results showed that fipronil alters sperm functional-related proteins and therefore influences male fertility. This study elucidates the possible reproductive toxic hazards associated with male infertility through aberrant suppression of sperm proteins.
Fipronil is a versatile phenylpyrazole insecticide with a wide range of applications such as public health, veterinary medicine, and agriculture. Fipronil disrupts the central nervous system of insects by targeting gamma-aminobutyric acid mediated chloride channel and glutamate-activated chloride channels, and it is more effective against insects than against mammals. Recently, considerable literature has accumulated regarding fipronil residues and its degradants detected in ecosystems such as urban surfaces, streams, groundwater, soil, and indoor dust1,2,3. According to existing research, fipronil residues and its degradants are prevalent in municipal wastewater treatment plants, river sediments, and urban estuarine sediments in USA4,5. Furthermore, in Vietnam, the threshold concentration of fipronil reviewed by the US EPA was detected in lake water6,7. A particular concern is that fipronil can be exposed and affect daily life such as through groundwater used for public supply8. In other words, although fipronil was reported to exhibit selective toxicity to undesirable organisms, the potential risks of human environmental exposure have become an increasing global concern. It is currently well established that fipronil exerts negative effects on various aspects of nontarget organisms, including mammals. A previous study reported that fipronil alters hormone levels and the estrous cycle length; fipronil-affected hormonal imbalance may result in a lower pregnancy rate in rats9. Moreover, aggravated ROS levels and DNA damage caused by fipronil increase apoptosis and cell cycle arrest during porcine oocyte maturation10. Even in embryos, it was found that fipronil depressed in vitro and in vivo developmental abilities in the mouse11. As well, the oxidant activity of fipronil in rats was found to suppress sperm production, reducing the epididymal sperm count12. Furthermore, fipronil increases DNA damage and apoptosis of rat spermatozoa, resulting in negative consequences on male fertility13. In particular, the previous study reported that fipronil may be involved in capacitation by binding with GABAA receptors, supplanting GABA, which may play a vital role as a putative modulator of sperm14. Despite these potential toxicological effects, studies demonstrating the molecular mechanisms of fipronil action on sperm function are still nascent. Therefore, in this study, we performed proteomic profiling of fipronil-treated spermatozoa and then applied bioinformatics to screen proteomic variations and the mechanisms of action of fipronil in mouse spermatozoa.
As shown in Table 1, all parameters of motility, i.e., total sperm motility (MOT) and progressive motility (PR) were significantly decreased with 1 µM fipronil treatment (p < 0.05). Among the kinematic parameters, average path velocity (VAP) was significantly decreased at lowest concentration of fipronil treatment (0.1 µM; p < 0.05). Curvilinear velocity (VCL) and straight-line velocity (VSL) showed a decreasing trend commencing at 1 µM (p < 0.05). In addition, beat cross frequency (BCF) and amplitude of lateral head displacement (ALH) decreased from 10 µM fipronil (p < 0.05). There was no significant difference in mean dance (DNM). In the measurements of capacitation status, the percentage of capacitated spermatozoa (B pattern) showed a sharp decline only at the highest concentration of fipronil (300 µM; p < 0.05). Meanwhile, the number of acrosome-reacted sperm (AR pattern) were drastically increased from 10 µM fipronil treatment, but there was no significant change in the non-capacitated pattern (F pattern) (p < 0.05; Table 1). The intracellular ATP level showed deficiencies at fipronil concentrations of 1, 10, 100, and 300 µM (p < 0.05; Fig. 1A). However, there was no statistical difference in cell viability between all treatment groups (Fig. 1B).
Comparison of intracellular ATP level and cell viability. (A) Differences in intracellular ATP levels after treatment with various concentration of fipronil. (B) Differences in cell viability between control and fipronil-treated groups. Data represent mean ± SEM, n = 3. Values with different superscripts (a and b) indicate significant difference between control and each treatment group (p < 0.05).
The measurement and identification of proteins and differential expression levels were accomplished through 2-DE, LC–MS/MS, and MASCOT database search. As shown in Table 2 and Supplementary Table S1, the proteomic analysis identified 14 differentially expressed proteins (DEPs) that were significantly downregulated (> threefold) due to fipronil treatment (p < 0.05). The normalized expression values of each DEP are shown in Fig. 2A.
Differentially expressed proteins (DEPs) and related pathways were measured and visualized. (A) Heat maps visualized the expression level of 14 DEPs after treatment with each concentration of fipronil. (B) GO and KEGG depict the discovered related pathways and protein counts in histogram.
The verification of DEPs was performed by western blotting and illustrated in Fig. 3. Among the 14 DEPs, we confirmed that the expression of cytochrome b5 domain-containing protein 1 (CYB5D1), axonemal dynein light intermediate polypeptide 1 (DNALI1), fatty acid-binding protein 9 (FABP9), 3-hydroxyisobutyrate dehydrogenase, mitochondrial (HIBADH), NADH dehydrogenase [ubiquinone] iron-sulfur protein 8, mitochondrial (NDUFS8), and radial spoke head protein 9 homolog (RSPH9) significantly decreased in all treatment groups. Moreover, the expression level of all six proteins were drastically decreased at the highest fipronil concentration (300 µM) relative to the control (p < 0.05; Fig. 3).
The histogram of fold change of the representative proteins relative to control. The expression ratios of representative 6 DEPs are depicted in the histogram. (A) NDUFS8, CYB5D1, HIBADH, DNALI1, RSPH9, and FABP9 were detected approximately at 26, 27, 35, 30, 57, and 15 kDa, respectively. (B) The expression level of NDUFS8. (C) The expression level of CYB5D1. (D) The expression level of HIBADH. (E) The expression level of DNALI1. (F) The expression level of RSPH9. (G) The expression level of FABP9. Data represent mean ± SEM, n = 3. Values with different superscripts (a, b, c, d, and e) indicate significant difference between control and each treatment group (p < 0.05).
The DEPs were analyzed by the DAVID database using GO and KEGG analyses. The seven biological processes, eight cellular components, and one molecular function related pathways were categorized (p < 0.05; Table 3 and Fig. 2B). In addition, KEGG enrichment analysis confirmed 12 related pathways, including metabolic pathways and oxidative phosphorylation process (p < 0.05; Table 3 and Fig. 2B).
The PPI network was constructed using Cytoscape with the STRING network. The results of the interaction analysis are summarized in Fig. 4. A total of 19 nodes and 40 edges were identified in Mus musculus. HIBADH, 10 kDa heat shock protein, mitochondrial (HSPE1), Isoform Nuclear of Phospholipid hydroperoxide glutathione peroxidase (GPX4), DNALI1, RSPH9, Isoaspartyl peptidase/L-asparaginase (ASRGL1), NDUFS8, Cytochrome c1, heme protein, mitochondrial (CYC1), and NADH dehydrogenase [ubiquinone] iron-sulfur protein 6, mitochondrial (NDUFS6) were found to interact with each other. The related signaling pathway of 10 proteins [RSPH9, DNALI1, Prelamin-A/C (LMNA), Microtubule-actin cross-linking factor 1 (MACF1), FABP9, HSPE1, ASRGL1, NDUFS8, CYC1, and GPX4] was constructed. All 10 DEPs interact with various proteins and are significantly related to cell processes such as sperm motility, reproductive cell development and maturation, and fertilization. Especially, it has been found that proteins LMNA and GPX4 are associated with diseases such as genetic disorder, embryo death, and infertility (p < 0.05; Fig. 5).
Construction of protein–protein interaction (PPI) network between DEPs. The schematic of PPI network among the 14 DEPs in Mus musculus is described using STRING. The 14 DEPs are circled in lemon and the 5 proteins in blue-green circles evidence new proteins that are functionally associated with 14 DEPs. A confidence cutoff for interaction was set to medium (0.400).
Construction of signaling pathway between DEPs. The signaling pathway of DEPs is illustrated using Pathway Studio. The DEPs are highlighted in blue background.
Fipronil is a widespread insecticidal ingredient that interrupts GABA and GluCl channels, impairing the central nervous system of the target organism. Although fipronil may be less toxic to mammals due to the absence of GluCl channels, recent studies have consistently shown that exposure to fipronil is associated with adverse effects in mammals9,10,11,12,13. Nevertheless, there is still scarce information regarding the potential impact of fipronil on sperm functions and the underlying molecular mechanisms. To explore the mechanism of fipronil action on the sperm proteome, we first evaluated the comprehensive sperm function.
Evaluation of sperm motility and kinematics revealed that fipronil treatment overall declined the sperm motility and kinematics parameters. Fipronil-treated spermatozoa exhibited decreased total motility, progressive motility, VCL, VSL, VAP, BCF, and ALH, although the total motility value was higher than the lower reference limit proposed by the WHO and thus estimated as normal motility values (Table 1). Similar to our results, a previous study demonstrated that perinatal exposure to fipronil reduced sperm motility and increased the number of nonmobile sperm in male rats15. Bae and Kwon14 also showed that fipronil suppressed the overall sperm motility and kinematics. Therefore, fipronil may affect sperm motility and kinematics. The highly motile specialized organelle flagella are components of the sperm tail and enable the sperm to produce locomotion. Frequently, axoneme deficiencies in sperm cause poor motility and results in infertility. The flagella of sperm are primarily organized with microtubules (MTs) composed of “9 + 2” axonemal complexes. The axonemal complexes consist of nine outer doublet microtubules, and two central pair interacts with and regulates ciliary movement. In the interaction, T-shaped radial spoke complexes with outer and inner dynein arms (ODA and IDA) assembling nine MTs with a central singlet. In this structure, sperm flagellar motility was coordinated by numerous dyneins in the axoneme using ATP16. In the present study, the expression levels of four DEPs, CFAP161, CYB5D1, DNALI1, and RSPH9 were significantly decreased, and CYB5D1, DNALI1, and RSPH9 positively correlated with sperm motility and kinematics, except DNM (Table 2, Supplementary Table S2, and Fig. 2A and 3). CYB5D1 is an RS component and coordinates the beating of flagella. A recent study showed that CYB5D1 mutation and knockdown zebrafish impaired coordinated ciliary beating in the olfactory epithelium and otic vesicle17. The component of dynein molecular motors, DNALI1, was previously reported to be expressed in testes, spermatids, and spermatocytes as well as the flagella of spermatozoa supporting axonemal dynein18,19. DNALI1 is localized along the IDA, which regulates flagellar beating, and the function of DNALI1 in the flagella of spermatozoa is still nascent, suggesting that it is involved in the generation of sperm beating20,21. In the DNALI1 knockdown zebrafish, ciliopathy phenotypes and depleted ODA has been observed, and revealed to be related to primary ciliary dyskinesia (PCD)22. PCD is a multisystemic cilia or flagella disorder of cells, including sperm, and closely related to male infertility. It is also known that RSPH9, which was differently expressed in our study, is linked to PCD, and a mutation of RSPH9 in patients resulted in RS abnormalities23. Previous studies have reported that in the “9 + 2” axonemal complex structure, the specific protein RSPH9 is a component of the axoneme and required for RS assembly and modulates the movement of axoneme in Chlamydomonas, mice, and humans24,25. Although the precise mechanism of Cilia- and flagella-associated protein 161 (CFAP161) remains nascent, CFAP161 was detected in testes26. A recent study reported that CFAP161 localized at the axoneme and was observed in the late stages of spermiogenesis and in spermatozoa. Hence, CFAP161-knockdown zebrafish lost ODAs, reducing the beating frequency of flagella27. Through this study, it was possible to confirm that the expression of CFAP161 was inhibited by fipronil. Synthetically, evidence shows that fipronil can alter sperm motility-related proteins. Consistent with this notion, the result of decreasing sperm motility and kinematics may be influenced by the diminished motility and the decreased expression of flagellum structure-related proteins induced by fipronil.
As a result of the capacitation status, during the capacitation process pattern was significantly decreased, and the non-capacitated spermatozoa tended to decrease, but no significant results were derived. Especially, the number of acrosome-reacted spermatozoa was anomalously increased, and according to the results of the present study, an approximately four-fold change was measured at the highest fipronil concentration (300 µM) relative to the control (Table 1). In previous research, the acrosome reaction was initiated by GABA in human spermatozoa28, and another earlier study has also demonstrated that fipronil binds with GABAA receptor instead of GABA1. Based on our study, we infer that spontaneous acrosome reaction was increased by fipronil. Furthermore, of the 14 identified DEPs, the expression of FABP9, HSPE1, and LMNA was decreased and related to sperm head. FABP9 is a testis- and germ cell-specific fatty acid-binding protein expressed during postnatal testicular development29,30 and identified as a perinuclear theca enriched in proteins in sperm31. Accordingly, Selvaraj et al.29 reported the presence of FABP9 in the post acrosomal region of the perinuclear theca, and Fabp9 knockout mice exhibited severe morphological abnormality primarily in the sperm head, even including fertile ones, and there was no significant change in litter size. Our western blot analysis revealed a significant decrease in the expression level of FABP9 in the 0.1 µM treatment group and a sharp decline in the 300 µM treatment group (Fig. 3G). HSPE1, also known as heat shock protein 10, has been expressed in the testes and epididymis of mouse spermatozoa32,33. Previously, HSPE1 was observed in the surface of sperm and predicted to play a role in capacitation, and after capacitation, the expression level of sperm surface HSPE1 was drastically increased33. However, as shown in Fig. 2A, the expression of HSPE1 significantly decreased despite an increase in the acrosome reaction. Prelamin-A/C is a component of the nuclear lamina; Shen et al.34 reported that LMNA plays an essential role in spermatogenesis for the structural and functional development of sperm head and acrosome. Moreover, spermatogenesis was found to be disrupted accompanied by increased apoptosis in Lmna–/– mice35. Based on these data, we infer that the increase in acrosome reaction indicates that fipronil-induced a spontaneously abnormal complex reaction of the acrosome. On this basis, it is speculated that the fipronil-induced silencing of sperm head-related DEPs affected the abnormally induced acrosome reaction (Table 1 and Fig. 3G). Therefore, we suggest a further study with more focus on the increased spontaneous acrosome reaction.
During the various phases of differentiation, the spermatozoa of individual taxa have a conformational pattern optimized for the penetration of the zona pellucida. During differentiation, a deficiency of the morphological perspective of sperm encompasses adverse outcomes on male fertility. From the structural point of view, numerous proteins compose the spermatozoa, and their complementary cooperation enables the sperm cell to maintain and generate the structure and movement biomechanically. Of the 14 identified DEPs, MACF1 has been known to be primarily involved in cytoskeleton (Table 2, Fig. 2A and 3). Previously, MACF1 was known as an important cytoskeletal crosslinking protein that binds actin filaments to microtubules36,37. Furthermore, another previous reported that deletion of Macf1 attenuates cytoskeletal arrangement38, and knockdown of Macf1 results in less stable F-actin and microtubule arrangement by the RhoA/ROCK1 signaling pathway in mice39. MACF1 mRNA was expressed in the mouse embryo37, and it is also involved in embryonic development based on the finding that Macf1−/− embryos died before birth40. Corroborating the previous findings, our results also elucidated that MACF1 binds with actin filament and is related to embryonic development (Fig. 5).
We discovered that, despite exhibiting a dose-dependently diminished tendency, the intracellular ATP level significantly decreased and reached a plateau of decline at 1 µM fipronil treatment (Fig. 1A). Of the 14 DEPs, the proteins CYC1, NDUFS6, NDUFS8 were related to sperm metabolism (Table 2, Fig. 2A, 3, and 5). ASRGL1 is an enzyme that catalyzes the hydrolysis of l -asparagine that was previously found to be expressed during oocyte development41 and in the sperm tail of rats and humans. Hashemitabar et al.42 also found that ASRGL1 was highly expressed in cases of asthenozoospermia. However, the gene expression in asthenozoospermia cases was downregulated compared with that in normozoospermia cases43. In our study, the protein expression level was downregulated after treatment with fipronil and positively correlated with sperm motility (Table 1, Supplementary Table S2, and Fig. 2A). CYC1 is a subunit of mitochondrial complex III of the mitochondrial electron transport chain, which is linked to oxidative phosphorylation and has been reported to be associated with sperm ROS production44,45. Downregulation of CYC1 reduces ATP production, whereas AMP production increases. Thus, knockdown of CYC1 sharply decreased the activation of mitochondrial complex III and activated AMPK in human breast cancer cells46. Moreover, an abnormal processing of CYC1 increased ROS production accompanied by affecting fertility in mice44. Therefore, it is plausible that diminished expression of CYC1 may related with deficiencies in ATP production (Fig. 1A). HIBADH plays a major role in catabolism, and a previous study showed that HIBADH, expressed during spermiogenesis in humans, is localized in the head and mid-piece of mature human spermatozoa, and diminished expression of HIBADH was related to low sperm motility47. Furthermore, HIBADH was expressed in the mid-piece of bull sperm and can be used as a motility-related marker48. Our results support the evidence from previous observations. In our study, the expression level of HIBADH positively correlated with sperm motility and intracellular ATP levels (Supplementary Table S2). NDUFS6 and NDUFS8 are subunits of the NADH dehydrogenase (ubiquinone) complex I that is involved in the respiratory electron transport chain. In Ndufs6−/− mice, bone marrow-derived mesenchymal stem cells exhibited increased intracellular and mitochondrial ROS levels49. NDUFS8 was upregulated in frozen-thawed sperm of Bos grunniens in contrast to that in fresh sperm50. Therefore, we speculate that decreased intracellular ATP levels are associated with the related DEPs. Although CYB5D1, DNALI1, FABP9, and HIBADH positively correlated with cell viability, no significant reduction was observed in cell viability between each of the treatment groups compared with the control (Fig. 1B and Supplementary Table S2).
The selenoprotein phospholipid hydroperoxide glutathione peroxidase GPX4 is an antioxidant enzyme that protects cells from ferrotopsis by preventing lipid peroxidation51,52. Previously, GPX4 was found to be primarily localized in the mitochondria of sperm and may play a role in antioxidant activity52. Moreover, spermatocyte specific Gpx4 knockout male mice showed abnormalities in spermatozoa such as mitochondrial dysfunction, decreased sperm cell count, and progressive motility, resulting in infertility53. Consistent with the literature, our study showed that GPX4 was related to flagellated sperm motility and fertilization (Fig. 5). These results reflect those of Ryu et al.54 who also reported that GPX4 was related to sperm abnormality, male and female infertility, and embryo death.
In the present study, we elucidated the alterations in sperm function and function-related protein expression in fipronil-affected spermatozoa. Briefly, fipronil diminished sperm motility, kinematics, sperm capacitation, and intracellular ATP level, and the most obvious finding of this study was that fipronil caused alteration in sperm proteins accompanied by an increase of acrosome reaction. Based on these results, we speculate that the interaction between diminished DEPs is involved in the functional and structural dysfunction of spermatozoa. These findings have significant implications for understanding how fipronil causes deficiencies in spermatozoa and how it is linked to male fertility. Therefore, these proteins may be used as novel predictive biomarkers in the estimation of sperm molecular deficiencies. Ultimately, our study demonstrated that fipronil negatively affects the expression of sperm function-related proteins as well as directly affects sperm physiological parameters such as motility and acrosome reaction. While our study provides valuable insights into the alterations in sperm function and protein expression, the exact mechanisms of fipronil on spermatozoa remain not fully understood. However, our findings confirmed that the identified DEPs are associated with specific mechanisms such as the flagellum and mitochondria. Therefore, further studies are suggested to reveal the exact mechanisms of fipronil focused on specific sperm pathways. Altogether, these findings shed new light on the molecular mechanisms of spermatozoa, contributing to an in-depth understanding of fipronil-induced sperm dysfunctions and emphasizing the potential toxicity of fipronil to male fertility.
All animal studies were conducted in compliance with the National Research Council's ‘Guide for the Care and Use of Laboratory Animals’ and followed the ARRIVE guidelines for ethical research. This study received ethical approval from the Institutional Animal Care and Use Committee of Kyungpook National University, with approval number KNU 2019–0089.
To induce sperm capacitation, 0.4% bovine serum albumin was added to modified Tyrode’s medium (97.84 mM NaCl, 1.42 mM KCl, 0.47 mM MgCl2·6H2O, 0.36 mM NaH2PO4·H2O, 5.56 mM d -glucose, 25 mM NaHCO3, 1.78 mM CaCl2.2H2O, 24.9 mM Na-lactate, 0.47 mM Na-pyruvate, 2 µg/mL gentamycin, and 0.005 mM phenol red). Unless otherwise stated, all chemicals and reagents were obtained from Sigma-Aldrich, St Louis, MO, USA.
ICR male mice (Nara Biotech, Seoul, Korea) were maintained individually providing ad libitum feeding (Cargill Agri purina, Inc., Seongnam, Korea) and water under the conditions of 21 °C ± 2 °C and 40%–60% humidity in 12-h light/dark cycles. Epididymal sperm were collected from ICR mice at the age of 8–12 weeks. Briefly, fat-removed cauda epididymis was placed in cell culture dishes with medium and punctured using a syringe to release spermatozoa. Then, the entire collected epididymal spermatozoa were incubated with various concentrations of fipronil (0.1, 1, 10, 100, and 300 μM), including a control, at 37 °C and 5% CO2 for 90 min to induce capacitation. The concentration of fipronil was used according to previous studies10,11.
The computer-assisted sperm analysis (CASA) system (FSA2016, Medical supply, Seoul, Korea), an OLYMPUS BX43 phase-contrast microscope (Olympus, Tokyo, Japan), and CMOS CAMERA, 2048 × 1536 (300 M pixel), 60 Frame (Medical supply) were used to evaluate sperm motility and kinematic parameters [MOT, sperm motility (%); Progressive, progressive motility (%); VCL, curvilinear velocity (μm/s); VSL, straight-line velocity (μm/s); VAP, average path velocity (μm/s); BCF, beat cross frequency (Hz); DNM, mean dance (µm); and ALH, amplitude of lateral head displacement (µm)]. Briefly, a 10 μL sample was transferred to a preheated Makler counting chamber (Sefi Medical Instruments, Haifa, Israel) and observed. Capacitation status was evaluated by combined Hoechst33258/chlortetracycline fluorescence staining. After the induction of capacitation, 135 µL of the sample was incubated with 15 μL of H33258 solution (10 μg H33258/mL PBS) for 10 min at RT. Then, 250 µL of 2% (w/v) polyvinylpyrrolidone was added and washed. The sample was then resuspended in 100 μL of PBS and 100 μL of CTC solution (750 mM CTC in 5 μL buffer: 20 mM Tris, 130 mM NaCl, and 5 mM cysteine, pH 7.4). After incubation for 20 min at 4 °C, an OLYMPUS BX43 epifluorescence illumination using ultraviolet BP 340–380/LP 425 and BP 450–490/LP 515 excitation/emission filters (Olympus) was used to evaluate the capacitation status. Subsequently, the spermatozoa were categorized into three groups (AR, acrosome-reacted; B, capacitated; and F, non-capacitated).
Intracellular ATP and cell viability levels were measured using ATP assay kit (Abcam, Cambridge, UK) and cell cytotoxicity assay kit (Abcam), respectively, according to the manufacturers’ instructions. The absorbance of intracellular ATP and cell viability was measured using GloMax Discover (Promega, Madison, WI, USA).
Two-dimensional electrophoresis was performed to elucidate differentially expressed proteins. To extract proteins, the samples were incubated in a rehydration buffer at 4 °C for 1 h (7 M urea, 2 M thiourea, 4% CHAPS (w/v), 0.05% Triton X-100, 24 μM PMSF, 1% octyl β-d -glucopyranoside, 20 mM DTT, 0.5% IPG buffer, and 0.005% bromophenol blue). Then, 250 µg of solubilized protein was placed on Immobiline DryStrip (pH 3–11 NL, 24 cm; Cytiva, Marlborough, MA, USA) for 12 h at 4 °C. First-dimension electrophoresis was performed using an IPGphor isoelectric focusing apparatus. The strips were isoelectric focused at 100 V for 1 h, 200 V for 1 h, 500 V for 1 h, 1000 V for 1 h, 5000 V for 1.5 h, 8000 V for 1.5 h, and 8000–90,000 V for 1 h. After isoelectric focusing, the strips were equilibrated with equilibration buffer A [6 M urea, 75 mM Tris–HCl (pH 8.8), 30% (v/v) glycerol, 2% (w/v) SDS, 0.002% (w/v) bromophenol blue, and 2% (w/v) DTT], and the second equilibration was performed using equilibration buffer B [6 M urea, 75 mM Tris–HCl (pH 8.8), 30% (v/v) glycerol, 2% (w/v) SDS, 0.002% (w/v) bromophenol blue, and 2.5% (w/v) iodoacetamide]. Second-dimension electrophoresis was performed using 12.5% (w/v) SDS-PAGE gels, with the strips being isoelectric focused at 100 V for 1 h and 500 V until the front of bromophenol blue began to move off the end of the gels. The gels were silver-stained according to the manufacturer’s instructions (Amersham Biosciences, Piscataway, NJ, USA), with the molecular weights ranging from 6.5 to 200 kDa and the pH levels ranging from 3 to 11. GS-800 calibrated Imaging Densitometer (Bio-Rad) was used to scan the gel for comparison and identification of the spots between each group. Then, the numerical value of protein expression level was computed using the PDQuest 8.0 software (Bio-Rad).
Proteins were subjected to in-gel trypsin digestion. Under shaking for 5 min, the excised gel spots were destained with 100 μL of destain solution (30 mM potassium ferricyanide, 100 mM sodium thiosulfate), and the solution was removed. Then, the gel spots were incubated with 200 mM ammonium bicarbonate for 20 min. The gel fragments were dried in a vacuum centrifuge after dehydration with 100 μL of acetonitrile. This procedure was repeated three times. The dried gel fragments were then rehydrated for 45 min on ice with 20 μL of 50 mM ammonium bicarbonate containing 0.2 μg of modified trypsin (Promega). After removing the solution, 70 μL of 50 mM ammonium bicarbonate was added. Protein digestion was performed overnight at 37 °C. Finally, a nano C18 column was used to desalinate the peptide solution.
Custom-made chromatographic columns were used for desalting and concentration of the peptide mixture before mass spectrometric analysis. A column consisting of 100–300 nL of Poros reverse-phase R2 material (20–30 μm bead size, PerSeptive Biosystems) was packed in a constricted GELoader tip (Eppendorf, Hamburg, Germany). A 10-mL syringe was used to force the liquid through the column by applying gentle air pressure. Next, 30 μL of the peptide mixture from the digested supernatant was diluted in 30 μL of 5% formic acid, loaded onto the column, and washed with 30 μL of 5% formic acid. For analyses by tandem mass spectrometry (MS/MS), the peptides were eluted with 1.5 μL of 50% methanol/49% H2O/1% formic acid.
The samples were resuspended in 0.1% formic acid in DW for using Ultimate 3000 (Thermo Fisher Scientific, Inc.). An autosampler was used to load 2-μL aliquots of the peptide solution into a C18 column (75 μm × 15 cm, particle size 2 μm) at a flow rate of 300 nL/min. The mobile phase A contained 0.1% formic acid in DW, and the mobile phase B contained 0.1% formic acid in 90% acetonitrile. The LC gradient elution started from 5 to 95% mobile phase B within 47.5 min and held at 95% for 5.0 min. Finally, back to 5% mobile phase B for 5.0 min. The mass spectrum (MS) scan range was 150 − 2000 m/z.
An MS/MS ion search was assigned as the ion search preference in the MASCOT software (version 2.4.1, Matrix Science, Boston, MA, USA). Peptide fragment files were obtained from the peptide peaks in ESI–MS by ESI–MS/MS. Trypsin was selected as an enzyme with a maximum of two missed cleavage sites. ESI-TRAP was selected as the instrument type. The peptide fragments were searched based on the database using the MASCOT software (version 2.4.1, Matrix Science) and FASTA search engine, and the search was limited to Mus musculus taxonomy in NCBInr and UniprotKB/swissprot databases. The mass tolerance was set at ± 10 ppm and ± 0.8 Da for the peptides and fragments, respectively. High scoring was defined as those greater than the default significance threshold in MASCOT (p < 0.05, peptide score, > 55).
The differentially expressed proteins were validated by western blotting analysis. Briefly, each group of samples were lysed using modified Laemmli sample buffer (315 mM Tris, 10% glycerol, 10% sodium dodecyl sulfate, 5% 2-mercaptoethanol, and 5% bromophenol blue). Then, the lysates were segregated by 12% SDS-PAGE (Mini PROTEIN Tetra Cell, Bio-Rad) and transferred to Immun-Blot polyvinylidene difluoride membranes (Bio-Rad). For blocking, the membranes were incubated with 3% ECL blocking agent (GE Healthcare, Chicago, IL, USA) in DPBS at RT for more than 2 h. Then, the membranes were washed with DPBS containing 0.01% Tween-20 (PBS-T) and incubated with primary antibodies diluted in 3% ECL blocking agent (GE Healthcare) [C6orf206 polyclonal antibody (1:300; Invitrogen, Thermo Fisher Scientific, Inc., Waltham, MA, USA), CYB5D1 polyclonal antibody (1:500; MyBioSource, San Diego, CA, USA), DNALI1 polyclonal antibody (1:4,000; MyBioSource), FABP9 monoclonal antibody (1:1000; Invitrogen), HIBADH antibody (1:500; MyBioSource), NDUFS8 polyclonal antibody (1:3000; Invitrogen)]. Anti-⍺-tubulin mouse antibody (1:5000; Abcam) was used as a loading control. After treatment, the membranes were washed with PBS-T. Anti-rabbit IgG, HRP-linked antibody (1:2000; Cell Signaling Technology, Danvers, MA, USA), or goat anti-mouse IgG H&L (HRP) (1:2000; Abcam) secondary antibody was diluted in 3% ECL blocking agent (GE Healthcare) and added. Protein bands were visualized using the iBright CL1500 imaging system (Invitrogen) with an ECL Substrate (Bio-Rad). Finally, the protein expression signals from each gel were measured using the Image Lab software, version 6.1.0 (Bio-Rad), and calculated relative to the control.
The biological process (BP), cellular component (CC), molecular function (MF), and pathway properties of DEPs were annotated using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases by Database for Annotation, Visualization, and Integrated Discovery 6.8 (DAVID) (https://david.ncifcrf.gov/home.jsp). The PPI network was constructed using the STRING applications (version 11.5) of Cytoscape (version 3.9.1)55. To visualize a signaling pathway, Pathway Studio (Version 12.4.0.5, Elsevier, Amsterdam, Netherlands) was implemented.
Data management and analysis were performed by one-way ANOVA and Pearson correlation using SPSS (Version 26.0, IBM, Armonk, NY, USA). Numerical data are represented as mean ± SEM. The level of significance was set at p < 0.05.
All the data generated or analysed during this study are included in this published article [and its supplementary information files].
Gan, J., Bondarenko, S., Oki, L., Haver, D. & Li, J. X. Occurrence of fipronil and its biologically active derivatives in urban residential runoff. Environ. Sci. Technol. 46, 1489–1495. https://doi.org/10.1021/es202904x (2012).
Article ADS CAS PubMed Google Scholar
Bonmatin, J. M. et al. Environmental fate and exposure; neonicotinoids and fipronil. Environ. Sci. Pollut. Res. Int. 22, 35–67. https://doi.org/10.1007/s11356-014-3332-7 (2015).
Article CAS PubMed Google Scholar
Testa, C. et al. Occurrence of Fipronil in residential house dust in the presence and absence of pets: a hint for a comprehensive toxicological assessment. J. Environ. Sci. Health Part B 54, 441–448. https://doi.org/10.1080/03601234.2019.1607133 (2019).
Maruya, K. A. et al. Bioanalytical and chemical-specific screening of contaminants of concern in three California (USA) watersheds. Heliyon 8, e09534. https://doi.org/10.1016/j.heliyon.2022.e09534 (2022).
Article PubMed PubMed Central Google Scholar
Lao, W., Tsukada, D., Greenstein, D. J., Bay, S. M. & Maruya, K. A. Analysis, occurrence, and toxic potential of pyrethroids, and fipronil in sediments from an urban estuary. Environ. Toxicol. Chem. 29, 843–851. https://doi.org/10.1002/etc.116 (2010).
Article CAS PubMed Google Scholar
EPA, U. Aquatic life benchmarks and ecological risk assessments for registered pesticides. (2020).
Wan, Y. et al. Neonicotinoids, fipronil, chlorpyrifos, carbendazim, chlorotriazines, chlorophenoxy herbicides, bentazon, and selected pesticide transformation products in surface water and drinking water from northern Vietnam. Sci. Total Environ. 750, 141507. https://doi.org/10.1016/j.scitotenv.2020.141507 (2021).
Article ADS CAS PubMed Google Scholar
Shi, L., Chen, L., Wan, Y., Zeng, H. & Xia, W. Spatial variation of fipronil and its derivatives in tap water and ground water from China and the fate of them during drinking water treatment in Wuhan, central China. Chemosphere 251, 126385. https://doi.org/10.1016/j.chemosphere.2020.126385 (2020).
Article CAS PubMed Google Scholar
Ohi, M., Dalsenter, P. R., Andrade, A. J. & Nascimento, A. J. Reproductive adverse effects of fipronil in Wistar rats. Toxicol. Lett. 146, 121–127. https://doi.org/10.1016/j.toxlet.2003.08.008 (2004).
Article CAS PubMed Google Scholar
Zhou, W. et al. Fipronil induces apoptosis and cell cycle arrest in porcine oocytes during in vitro maturation. Apoptosis 24, 718–729. https://doi.org/10.1007/s10495-019-01552-w (2019).
Article CAS PubMed Google Scholar
Sefcikova, Z. et al. Fipronil causes toxicity in mouse preimplantation embryos. Toxicology 410, 214–221. https://doi.org/10.1016/j.tox.2018.08.008 (2018).
Article CAS PubMed Google Scholar
Mazzo, M., Balieira, K. V. B., Bizerra, P. F. V. & Mingatto, F. E. Fipronil-induced decrease in the epididymal sperm count: oxidative effect and protection by vitamin E. Anim. Reprod. 15, 1223–1230. https://doi.org/10.21451/1984-3143-ar2017-0040 (2018).
Article PubMed PubMed Central Google Scholar
Khan, S., Jan, M. H., Kumar, D. & Telang, A. G. Firpronil induced spermotoxicity is associated with oxidative stress, DNA damage and apoptosis in male rats. Pestic. Biochem. Physiol. 124, 8–14. https://doi.org/10.1016/j.pestbp.2015.03.010 (2015).
Article ADS CAS PubMed Google Scholar
Bae, J. W. & Kwon, W. S. Investigating the effects of fipronil on male fertility: Insight into the mechanism of capacitation. Reprod. Toxicol. 94, 1–7. https://doi.org/10.1016/j.reprotox.2020.04.002 (2020).
Article CAS PubMed Google Scholar
de Barros, A. L. et al. Perinatal exposure to insecticide fipronil: effects on the reproductive system in male rats. Reprod. Fertil. Dev. 29, 1130–1143. https://doi.org/10.1071/rd15517 (2017).
Lin, J. & Nicastro, D. Asymmetric distribution and spatial switching of dynein activity generates ciliary motility. Science 360, eaar1968. https://doi.org/10.1126/science.aar1968 (2018).
Article CAS PubMed PubMed Central Google Scholar
Zhao, L. et al. Heme-binding protein CYB5D1 is a radial spoke component required for coordinated ciliary beating. Proc. Natl. Acad. Sci. USA https://doi.org/10.1073/pnas.2015689118 (2021).
Article PubMed PubMed Central Google Scholar
Dong, F. et al. Pih1d3 is required for cytoplasmic preassembly of axonemal dynein in mouse sperm. J. Cell Biol. 204, 203–213. https://doi.org/10.1083/jcb.201304076 (2014).
Article CAS PubMed PubMed Central Google Scholar
Rashid, S. et al. The murine Dnali1 gene encodes a flagellar protein that interacts with the cytoplasmic dynein heavy chain 1. Mol. Reprod. Dev. 73, 784–794. https://doi.org/10.1002/mrd.20475 (2006).
Article CAS PubMed Google Scholar
Inaba, K. & Mizuno, K. Sperm dysfunction and ciliopathy. Reprod. Med. Biol. 15, 77–94. https://doi.org/10.1007/s12522-015-0225-5 (2016).
Article CAS PubMed Google Scholar
Aprea, I. et al. Defects in the cytoplasmic assembly of axonemal dynein arms cause morphological abnormalities and dysmotility in sperm cells leading to male infertility. PLOS Genetics 17, e1009306. https://doi.org/10.1371/journal.pgen.1009306 (2021).
Article CAS PubMed PubMed Central Google Scholar
Austin-Tse, C. et al. Zebrafish ciliopathy screen plus human mutational analysis identifies C21orf59 and CCDC65 defects as causing primary ciliary dyskinesia. Am. J. Hum. Genetics 93, 672–686. https://doi.org/10.1016/j.ajhg.2013.08.015 (2013).
Castleman, V. H. et al. Mutations in radial spoke head protein genes RSPH9 and RSPH4A cause primary ciliary dyskinesia with central-microtubular-pair abnormalities. Am. J. Hum. Genetics 84, 197–209. https://doi.org/10.1016/j.ajhg.2009.01.011 (2009).
Zhu, L., Liu, H., Chen, Y., Yan, X. & Zhu, X. Rsph9 is critical for ciliary radial spoke assembly and central pair microtubule stability. Biol. Cell 111, 29–38. https://doi.org/10.1111/boc.201800060 (2019).
Article CAS PubMed Google Scholar
Abbasi, F. et al. RSPH6A is required for sperm flagellum formation and male fertility in mice. J. Cell Sci. https://doi.org/10.1242/jcs.221648 (2018).
Article PubMed PubMed Central Google Scholar
Shawki, H. H. et al. EFCAB2 is a novel calcium-binding protein in mouse testis and sperm. PLoS ONE 14, e0214687. https://doi.org/10.1371/journal.pone.0214687 (2019).
Article CAS PubMed PubMed Central Google Scholar
Beckers, A. et al. The highly conserved FOXJ1 target CFAP161 is dispensable for motile ciliary function in mouse and Xenopus. Sci. Rep. 11, 13333. https://doi.org/10.1038/s41598-021-92495-3 (2021).
Article ADS CAS PubMed PubMed Central Google Scholar
Shi, Q. X., Yuan, Y. Y. & Roldan, E. R. gamma-Aminobutyric acid (GABA) induces the acrosome reaction in human spermatozoa. Mol. Hum. Reprod. 3, 677–683. https://doi.org/10.1093/molehr/3.8.677 (1997).
Article CAS PubMed Google Scholar
Selvaraj, V. et al. Mice lacking FABP9/PERF15 develop sperm head abnormalities but are fertile. Dev. Biol. 348, 177–189. https://doi.org/10.1016/j.ydbio.2010.09.019 (2010).
Article CAS PubMed PubMed Central Google Scholar
Moradi , A. , Forootan , FS , Hosseini , M. & Pouresmaeili , F. Leydig cells expressing FABP9 in human testis .Hum.Antib.27 , 275–278 .https://doi.org/10.3233/HAB-190382 (2019).
Zhang, X. Z., Wei, L. L., Jin, H. J., Zhang, X. H. & Chen, S. R. The perinuclear theca protein Calicin helps shape the sperm head and maintain the nuclear structure in mice. Cell Rep. 40, 111049. https://doi.org/10.1016/j.celrep.2022.111049 (2022).
Article CAS PubMed Google Scholar
Hartman, D. J., Hoogenraad, N. J., Condron, R. & Høj, P. B. Identification of a mammalian 10-kDa heat shock protein, a mitochondrial chaperonin 10 homologue essential for assisted folding of trimeric ornithine transcarbamoylase in vitro. Proc. Natl. Acad. Sci. 89, 3394–3398. https://doi.org/10.1073/pnas.89.8.3394 (1992).
Article ADS CAS PubMed PubMed Central Google Scholar
Walsh, A. et al. Identification of the molecular chaperone, heat shock protein 1 (chaperonin 10), in the reproductive tract and in capacitating spermatozoa in the male mouse. Biol. Reprod. 78, 983–993. https://doi.org/10.1095/biolreprod.107.066860 (2008).
Article CAS PubMed Google Scholar
Shen, J. et al. Lamin A/C proteins in the spermatid acroplaxome are essential in mouse spermiogenesis. Reproduction 148, 479–487. https://doi.org/10.1530/REP-14-0012 (2014).
Article CAS PubMed Google Scholar
Alsheimer, M. et al. Disruption of spermatogenesis in mice lacking A-type lamins. J. Cell Sci. 117, 1173–1178. https://doi.org/10.1242/jcs.00975 (2004).
Article CAS PubMed Google Scholar
Gamper, I. et al. GAR22β regulates cell migration, sperm motility, and axoneme structure. Mol. Biol. Cell 27, 277–294. https://doi.org/10.1091/mbc.E15-06-0426 (2016).
Article CAS PubMed PubMed Central Google Scholar
Leung, C. L., Sun, D., Zheng, M., Knowles, D. R. & Liem, R. K. H. Microtubule actin cross-linking factor (Macf): A hybrid of dystonin and dystrophin that can interact with the actin and microtubule cytoskeletons. J. Cell Biol. 147, 1275–1286. https://doi.org/10.1083/jcb.147.6.1275 (1999).
Article CAS PubMed PubMed Central Google Scholar
Kodama, A., Karakesisoglou, I., Wong, E., Vaezi, A. & Fuchs, E. ACF7: an essential integrator of microtubule dynamics. Cell 115, 343–354. https://doi.org/10.1016/s0092-8674(03)00813-4 (2003).
Article CAS PubMed Google Scholar
Lin, X. et al. Knockdown of MACF1 inhibits the migration and cytoskeletal arrangement of pre-osteoclasts induced by simulated microgravity. Acta Astronaut. 190, 149–159. https://doi.org/10.1016/j.actaastro.2021.10.010 (2022).
Chen, H. J. et al. The role of microtubule actin cross-linking factor 1 (MACF1) in the Wnt signaling pathway. Genes Dev. 20, 1933–1945. https://doi.org/10.1101/gad.1411206 (2006).
Article CAS PubMed PubMed Central Google Scholar
Chavez, D. R., Lee, P. C. & Comizzoli, P. Oocyte meiotic competence in the domestic cat model: Novel roles for nuclear proteins BRD2 and NPM1. Front. Cell Dev. Biol. 9, 670021. https://doi.org/10.3389/fcell.2021.670021 (2021).
Article PubMed PubMed Central Google Scholar
Hashemitabar, M., Sabbagh, S., Orazizadeh, M., Ghadiri, A. & Bahmanzadeh, M. A proteomic analysis on human sperm tail: Comparison between normozoospermia and asthenozoospermia. J. Assist. Reprod. Genet. 32, 853–863. https://doi.org/10.1007/s10815-015-0465-7 (2015).
Article PubMed PubMed Central Google Scholar
Liu, X. X., Cai, L. & Liu, F. J. An in silico analysis of human sperm genes associated with asthenozoospermia and its implication in male infertility. Medicine 97, e13338. https://doi.org/10.1097/MD.0000000000013338 (2018).
Article CAS PubMed PubMed Central Google Scholar
Lu, B. et al. A mutation in the inner mitochondrial membrane peptidase 2-like gene (Immp2l) affects mitochondrial function and impairs fertility in mice1. Biol. Reprod. 78, 601–610. https://doi.org/10.1095/biolreprod.107.065987 (2008).
Article CAS PubMed Google Scholar
Dröse, S., Brandt, U. & Wittig, I. Mitochondrial respiratory chain complexes as sources and targets of thiol-based redox-regulation. Biochim. Biophys. Acta 1344–1354, 2014. https://doi.org/10.1016/j.bbapap.2014.02.006 (1844).
Han, Y. et al. CYC1 predicts poor prognosis in patients with breast cancer. Dis Mark. 2016, 3528064. https://doi.org/10.1155/2016/3528064 (2016).
Tasi, Y. C. et al. Characterization of 3-hydroxyisobutyrate dehydrogenase, HIBADH, as a sperm-motility marker. J. Assist. Reprod. Genet. 30, 505–512. https://doi.org/10.1007/s10815-013-9954-8 (2013).
Article PubMed PubMed Central Google Scholar
Zhang, S. et al. The g.-165 T>C rather than methylation is associated with semen motility in chinese holstein bulls by regulating the transcriptional activity of the HIBADH gene. PLoS ONE 10, 01276. https://doi.org/10.1371/journal.pone.0127670 (2015).
Zhang, Y. et al. Adult mesenchymal stem cell ageing interplays with depressed mitochondrial Ndufs6. Cell Death Dis. 11, 1075. https://doi.org/10.1038/s41419-020-03289-w (2020).
Article CAS PubMed PubMed Central Google Scholar
Fan, Y. et al. TMT-based quantitative proteomics analysis reveals the differential proteins between fresh and frozen-thawed sperm of yak (Bos grunniens). Theriogenology 200, 60–69. https://doi.org/10.1016/j.theriogenology.2023.01.024 (2023).
Article CAS PubMed Google Scholar
Maiorino, M., Conrad, M. & Ursini, F. GPx4, lipid peroxidation, and cell death: Discoveries, rediscoveries, and open issues. Antioxid. Redox. Signal 29, 61–74. https://doi.org/10.1089/ars.2017.7115 (2018).
Article CAS PubMed Google Scholar
Imai, H. et al. Failure of the expression of phospholipid hydroperoxide glutathione peroxidase in the spermatozoa of human infertile males. Biol. Reprod. 64, 674–683. https://doi.org/10.1095/biolreprod64.2.674 (2001).
Article CAS PubMed Google Scholar
Imai, H. et al. Depletion of selenoprotein GPx4 in spermatocytes causes male infertility in mice. J. Biol. Chem. 284, 32522–32532. https://doi.org/10.1074/jbc.M109.016139 (2009).
Ryu, D. Y., Rahman, M. S. & Pang, M. G. Determination of highly sensitive biological cell model systems to screen BPA-related health hazards using pathway studio. Int J. Mol. Sci. 18, 1909. https://doi.org/10.3390/ijms18091909 (2017).
Article CAS PubMed PubMed Central Google Scholar
Doncheva, N. T., Morris, J. H., Gorodkin, J. & Jensen, L. J. Cytoscape StringApp: Network analysis and visualization of proteomics data. J Proteome Res. 18, 623–632. https://doi.org/10.1021/acs.jproteome.8b00702 (2019).
Article CAS PubMed Google Scholar
This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (Ministry of Science and ICT) (NRF-2019R1F1A1049216).
Department of Animal Science and Biotechnology, Kyungpook National University, Sangju, Republic of Korea
Jeong-Won Bae & Woo-Sung Kwon
Research Institute for Innovative Animal Science, Kyungpook National University, Sangju, Republic of Korea
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
All authors contributed toward revisions that were critical to the intellectual content and approved the final version for publication. J.W.B.; Investigation, Validation, Methodology, Data curation, Writing—original draft. W.S.K; Investigation, Validation, Methodology, Data curation Writing—review & editing, Project administration, Funding acquisition.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Bae, JW., Kwon, WS. Proteomic analysis of fipronil-induced molecular defects in spermatozoa. Sci Rep 14, 7668 (2024). https://doi.org/10.1038/s41598-024-57876-4
DOI: https://doi.org/10.1038/s41598-024-57876-4
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
Monosultap 90% SP Sign up for the Nature Briefing: Translational Research newsletter — top stories in biotechnology, drug discovery and pharma.