• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Hydralazine represses Fpn ubiquitination to rescue injured neurons via competitive binding to UBA52

    2024-03-21 05:51:14ShengyouLiXueGaoYiZhengYujieYangJianboGaoDanGengLingliGuoTengMaYimingHaoBinWeiLiangliangHuangYitaoWeiBingXiaZhuojingLuoJinghuiHuang
    Journal of Pharmaceutical Analysis 2024年1期

    Shengyou Li , Xue Gao , Yi Zheng , Yujie Yang, Jianbo Gao, Dan Geng, Lingli Guo,Teng Ma, Yiming Hao, Bin Wei, Liangliang Huang, Yitao Wei, Bing Xia, Zhuojing Luo,Jinghui Huang

    Department of Orthopedics, Xijing Hospital, Fourth Military Medical University, Xi'an, 710032, China

    Keywords:Ferroptosis UBA52 Ferroportin Ubiquitination Hydralazine Peripheral nerve injury

    A B S T R A C T A major impedance to neuronal regeneration after peripheral nerve injury (PNI) is the activation of various programmed cell death mechanisms in the dorsal root ganglion.Ferroptosis is a form of programmed cell death distinguished by imbalance in iron and thiol metabolism, leading to lethal lipid peroxidation.However, the molecular mechanisms of ferroptosis in the context of PNI and nerve regeneration remain unclear.Ferroportin (Fpn), the only known mammalian nonheme iron export protein,plays a pivotal part in inhibiting ferroptosis by maintaining intracellular iron homeostasis.Here,we explored in vitro and in vivo the involvement of Fpn in neuronal ferroptosis.We first delineated that reactive oxygen species at the injury site induces neuronal ferroptosis by increasing intracellular iron via accelerated UBA52-driven ubiquitination and degradation of Fpn, and stimulation of lipid peroxidation.Early administration of the potent arterial vasodilator,hydralazine(HYD),decreases the ubiquitination of Fpn after PNI by binding to UBA52, leading to suppression of neuronal cell death and significant acceleration of axon regeneration and motor function recovery.HYD targeting of ferroptosis is a promising strategy for clinical management of PNI.

    1.Introduction

    Each year, millions of individuals worldwide experience peripheral nerve injury(PNI),which leads to long-term disability and high medical costs,resulting in negative socio-economic effects[1].PNI is usually caused by lacerations,contusions,and tractions that can arise from traumatic injuries, drug injection injuries, electrocution, infections, and disorders of metabolism [2].The peripheral nervous system possesses a certain ability to regenerate, enabling axonal regrowth after axon damage [3].Several studies have shown that treatment with natural medicinal compounds, drugs, and gene manipulation strategies can enable the recovery of nerve functions after mild-to-moderate injury [4].However, functional recovery is substantially lower after severe nerve injury.In cases of PNI,lesions at the interface between the central and peripheral nervous systems can result in neuronal loss of 35%-40% within two months after surgery, leading to a poor prognosis [5].Therefore, development of novel neuro-therapeutic approaches will require thorough comprehension of neurodegeneration/regeneration mechanisms,especially those undeciphered mechanisms potentially involved in PNI.

    Cell death is a basic biological process that is required for normal development, aging, and tissue homeostasis,and its dysregulation is observed in multiple disorders[6].Recent studies have identified ferroptosis in the pathophysiology of many diseases, which differs from other types of cell death at the levels of morphology,biochemistry, and genetics [7].In ferroptosis, iron-dependent signaling causes lipid peroxides to accumulate within cells,resulting in cell death.Ferroptosis was first observed after erastin treatment that selectively kill genetically engineered cells with an oncogenic mutation [8].Recent findings have demonstrated that ferroptosis is involved in multiple neuronal diseases [9] and its inhibition can prevent neurodegeneration [10].Nevertheless, the role of ferroptosis in the pathophysiology of PNI is yet to be elucidated.

    PNI is accompanied by oxidative damage and inflammation[11].Studies have recently shown that spinal nerve injury can generate reactive oxygen species(ROS)in microglial cells of the dorsal horn[12].In addition, the inhibition of oxidative stress improves functional recovery after PNI [13].Therefore, oxidative stress is considered as a primary contributor to neuronal injury and neuron loss[14].The number of residual sensory neurons following PNI is known to depend on the level of cell death within the dorsal root ganglia cells(DRGs).An in vivo study showed that up to 40%of DRGs undergo cell death within two months after PNI, compromising functional recovery [15].So far, previous studies have largely attributed neuronal loss after PNI to apoptosis.However, as our understanding of cell death has deepened over the past years,great attention has shifted to ferroptosis, a process that is triggered by ROS.Nevertheless, despite the well-defined role of ROS in PNI,whether ferroptosis is involved in the pathophysiology of nerve injury/repair has not been elucidated.In addition, the relationship between ferroptosis and functional recovery after PNI also remains unclear.The present study aimed to answer these questions.

    2.Material and methods

    2.1.Ethics approval

    The living conditions of all animals and all experimental procedures complied with the Guide for the Care and Use of Laboratory Animals (National Institutes of Health; Publication No.85-23,revised 1985)and were performed after approval from the Animal Research Committee of The Fourth Military Medical University,Xi'an,China(Approval No.:IACUC-20230062).All Sprague-Dawley(SD) rats used in the present study were sourced from the same center.

    2.2.Drugs

    H2O2and diphenyleneiodonium (DPI) were purchased from Sigma-Aldrich (St.Louis, MO, USA); and N-acetylcysteine (NAC),hydralazine(HYD),ferrostatin-1(Fer-1),and pyridoxal isonicotinoyl hydrozone(PIH)were purchased from MedChemExpress(Shanghai,China).Stock solutions of H2O2(13 mM) and NAC (30 mM) were dissolved in saline, while those of DPI (100 μM), Fer-1 (6 μM), and PIH(1 mM)were dissolved in dimethyl sulfoxide(DMSO).

    2.3.Development of sciatic nerve injury (SNI) model

    Adult male SD rats (200-250 g) were housed in groups of five per cage under controlled temperature(23±2°C),humidity(35%-60%), and 12-h light/12-h dark cycle.The rat model of SNI was developed using a modified version of a published protocol [16].Rats were anesthetized with intraperitoneal (i.p.) injection of 1%pentobarbital sodium (40 mg/kg; Amresco, Radnor, PA, USA).The left sciatic nerve was exposed with a 1 cm incision at the center of the thigh.The sciatic nerve was crushed at approximately 20 mm from the L4-L6 dorsal root ganglions in the orthogonal direction(2×20 s)to induce nerve injury.The incision was closed using nondegradable sutures.The same procedure was carried out in the sham group but without any induction of injury to the sciatic nerve.

    Rats were subsequently allocated to three groups at random:SNI with no treatment, immediate treatment, and delayed treatment, respectively.Treatment consisted of HYD (20 μg/kg body weight) via i.p.injections once daily for seven consecutive days,starting right after (immediate group) or one week after injury induction (delayed group).The sham and SNI with no treatment groups received equal volumes of saline via i.p.injections.The concentrations of HYD well below 0.3 mg/kg, which is lowest effective concentration to obtain antihypertensive effect.

    2.4.Detection of superoxide levels

    We examined superoxide levels by detecting 2-hydroxyethidine,a product of dihydroethidium(DHE;Invitrogen,Carlsbad,CA,USA)oxidation.Rats were administered with 3 mL of DHE(1 mg/mL,i.p.)2 h before and 8 h after operation,sacrificed 4 h later with sodium pentobarbital(40 mg/kg,i.p.)injection,and perfused with 0.9%saline and 200 mL of 4%paraformaldehyde(PFA;Servicebio,Wuhan,China).Tissue samples from each rat was post-fixed in 4%PFA for 6 h,dehydrated in 30%sucrose for 12 h and sectioned into 20-μm thick slices using a cryocut microtome (Leica CM3050S, Weztlar, Germany).Sections were mounted onto gelatin-coated slides with vector shield for imaging on a confocal microscope (A1+, Nikon,Tokyo, Japan).At least six fields from each section were analyzed using ImageJ(National Institutes of Health,Bethesda,MD,USA).

    2.5.Isolation and culture of DRGs

    DRGs were isolated from neonatal (postnatal day 0-1) SD rats using previously published protocols[17].DRGs were seeded onto rat tail collagen (Invitrogen) coated coverslips and cultured in neurobasal medium (Gibco,Grand Island,NY,USA)with 50 ng/mL nerve growth factor(R&D Systems,Minneapolis,MN,USA),2 mM Lglutamine (Gibco), 2% B27 (Gibco), and 1% penicillin-streptomycin(Procell, Wuhan, China).Fibroblasts in the cultures were eliminated with 10 mM cytosine arabinoside(Sigma-Aldrich)to obtain a pure DRGs culture.

    2.6.Assessment of cell viability

    DRGs were seeded in 24-well plates (3000 cells/well,n=4 per group) for 24 h, followed by H2O2treatment at the designated concentration for another 24 h.Erastin(20 μM;Sigma-Aldrich)was used as a positive control.Some cells were also treated with z-VADFMK (10 μM; MedChemExpress), Fer-1 (1 μM), or PIH (100 μM) in addition to H2O2.The LIVE/DEAD Cell Imaging Kit (Thermo Fisher Scientific Inc., Waltham, MA, USA) was used based on the manufacturer's instructions to estimate cell viability.Live (green) and dead (red)cells were quantified by ImageJ.

    2.7.Detection of intracellular ROS concentration

    ROS production was examined using 2,7-dichlorodihydrofluorescein diacetate (DCFH-DA; Beyotime; Shanghai, China).DRGs were cultivated in the DRGs culture medium for 24 h(8×104cells/well)and replaced with fresh medium containing the indicated concentrations of H2O2for a further 24 h.Subsequently,cells were incubated with 10 μM DCFH-DA for 15 min and imaged under a fluorescence microscope (Olympus, Tokyo, Japan).Fluorescence intensity was measured using ImageJ.

    2.8.Assessment of lactate dehydrogenase (LDH) and glutathione(GSH) levels

    We assessed the GSH and LDH concentrations using the S0055 and C0016 commercial kits from Beyotime, respectively, as described by the manufacturer.

    2.9.Apoptosis assay

    Apoptosis assay was conducted as published [18].Briefly, cells with and without H2O2treatment were stained with annexin-V(Beyotime) for 15 min, followed by propidium iodide (Beyotime)for 5 min at room temperature.Stained cells were visualized with flow cytometry.

    2.10.Heme measurement

    The QuantiChrom Heme Assay Kit (DIHM-250, BioAssay Systems, Bay Area, CA, USA) was used as described by the manufacturer to determine serum and tissue heme levels.Serum or tissue homogenates (10 μL each) were mixed with the appropriate reagents and the absorbance measured at 400 nm.Heme concentration of samples were derived from a standard curve of known heme concentrations.

    2.11.Lentivirus (Lv) and plasmid transfection

    Primary neurons cultured in 6-well plates with neurobasal media were transduced with Lv carrying ferroportin (Fpn) (pLVEf1a-Fpn-3Flag-Puro),vector only(pLV-Ef1a-3Flag-Puro),or short hairpin RNA (shRNA) targeting UBA52 (LV2-pGLVu6-shUBA52-puro) at an multiplicity of infection of 100.All Lv were purchased from Genechem(Shanghai,China).Stably transduced cells were selected with puromycin (50 μg/mL) after 72 h.Overexpression of Fpn was confirmed by quantitative real-time polymerase chain reaction(qPCR)and Western blot assays.All plasmid transfections were performed with Lipofectamine 2000 Transfection Reagent(Thermo Fisher Scientific Inc.)as described by the manufacturer.

    2.12.Assessment of iron levels

    Intracellular iron levels were estimated with FerroOrange probes (DojinDo, Kyushu Island, Japan).DRGs were cultivated in DRGs medium for 24 h (8 × 104cells/well), followed by replacement with medium containing erastin, H2O2, or HYD for another 24 h.Cells were stained with FerroOrange, imaged with a flurorescence microscope,and the color intensity assessed using ImageJ.

    2.13.Assessment of lipid peroxidation

    The levels of lipid peroxidation were examined using the C11-BODIPY581/591method (GlpBio, Montclair, CA, USA) and were further verified with detection of malondialdehyde (MDA, Beyotime) as instructed by the manufacturer.Cells were imaged with a fluorescence microscope and the number of green and red cells counted with ImageJ.

    2.14.qPCR

    Total RNA from the left DRGs (L4-L6) and sciatic nerve tissue was extracted and reverse transcribed with TRIzol (Thermo Fisher Scientific Inc.) and M-MLV Reverse Transcriptase (Thermo Fisher Scientific Inc.),respectively.qPCR was carried out on the 7500 qPCR system (Applied Biosystems) and the messenger RNA (mRNA)expression relative to β-actin was calculated with the 2-ΔΔCt method[19].Primer sequences for the evaluated genes are listed in Table S1.All qPCR analyses were conducted with no less than three independent samples.

    2.15.Western blot

    Tissue samples and cells were homogenized in ice-cold hypotonic lysis buffer containing proteinase inhibitors and resuspended in nuclear extract buffer containing proteinase inhibitors.Protein levels were quantified with the bicinchoninic acid protein assay kit(Beyotime).Proteins were denatured by boiling for 15 min, followed by sodium dodecyl-sulfate polyacrylamide gel electrophoresis(SDS-PAGE)separation and transfer to polyvinylidene fluoride membranes (Thermo Fisher Scientific Inc.).Membranes were blocked with 5% skim milk in tris buffered saline with tween-20(TBST; Servicebio) and incubated with the following primary antibodies at 4°C overnight: anti-SLC7A11 (ab175186; Abcam, Cambridge, UK), anti-glutathione peroxidase 4 (GPX4) (ab125066;Abcam),anti-transferrin receptor 1(TfR1,ab214039;Abcam),antiferritin H (FtH, ab183781; Abcam), anti-ferritin L (FtL, ab109373;Abcam), anti-Fpn (ab239583; Abcam), anti-GAP43 (ab16053,Abcam), and anti-β-actin (Cwbio, Taizhou, China).Subsequently,membranes were washed three times with TBST, incubated with the appropriate horseradish peroxidase-conjugated secondary antibodies at room temperature for 1 h, and the bound antibodies revealed with enhanced chemiluminescent substrates (Bio-rad).Each sample was assayed in triplicate.

    2.16.Mitochondrial morphology assays

    For ultrastructural analysis of mitochondrial morphology, cells were prepared for transmission electron microscopy (TEM) as follow:fixation in 2.5%glutaraldehyde in 0.1 M sodium cacodylate buffer (pH 7.4) and postfixed in 2% aqueous osmium tetroxide;gradual ethanol (30%-100%) and propylene oxide dehydration;Epon (Sigma-Aldrich) embedding at 60°C for 24 h; ultra-thin sectioning; section staining with 0.5% uranyl acetate (30 min) and 3%lead citrate(7 min)at 20°C;TEM analysis(Zeiss Libra 120 Plus,Carl Zeiss, Oberkochen, Germany); and imaging with a slow scan charge coupled device-camera and iTEM software (Olympus).

    2.17.Silver staining and qualitative protein analysis

    For silver staining,protein samples were separated on SDS-PAGE and processed with the Fast Silver Stain Kit (Beyotime) as instructed by the manufacturer.Qualitative protein analysis was performed by the Shanghai Bioprofile Institution (Shanghai, China).Protein samples were used for digestion, high performance liquid chromatography (HPLC), and mass spectrometry (MS) based on Bioprofile's protocol.The bioinformatics pipeline involved processing of raw data (filtering and quality control), alignment of experimental MS/MS data to published databases, and identification of proteins.

    2.18.Terminal deoxynucleotidyl transferase dUTP nick-end labeling(TUNEL) assay

    Percentage of cell death was determined by the TUNEL assay(Beyotime) as instructed by the manufacturer.Briefly, cells were fixed in 4%PFA,permeabilized with 0.3%Triton X-100(Servicebio)for 1 h, and incubated with the TUNEL detection solutions in the dark at 37°C for 1 h.Following nuclei staining with 4′,6-diamidino-2-phenylindole (DAPI; Servicebio), slides were coversliped and imaged with a fluorescence microscope (Olympus).

    2.19.Co-immunoprecipitation (Co-IP)

    Co-IP was performed as published[20].Cells were lysed with a weak radioimmunoprecipitation assay lysis buffer (Cell Signaling Technology,Danfoss,MA,USA),and the lysate pre-cleared with 50%protein A/G agarose for 1 h.2 μg of primary antibody was added to 500 mL of extracted proteins and incubated overnight at 4°C.Protein A/G agarose was added to pull down the immune complexes in a 4°C shaker for 4 h.Protein were eluted from the A/G beads through boiling in 1× loading buffer and used for immunoblotting analysis.

    2.20.Immunofluorescence staining

    Four weeks post-surgery, regenerated nerve segments were obtained from rat models of SNI,and sectioned longitudinally after paraffin embedding.The sections were incubated with the rabbit anti-GAP43 (1:400; ab16053, Abcam) primary antibody, followed by the goat anti-rabbit IgG cy3(1:200;Abcam)secondary antibody.Each incubation was preceded by a 5 min phosphate-buffered solution wash, and all experiments were repeated thrice.

    2.21.Functional assays

    Sciatic function index (SFI): to evaluate behavioral outcomes after nerve injury and repair,we tracked the gait of rats 1-4 weeks after surgery.Rats were allowed to walk inside a narrow passage(5.0 cm × 8.0 cm × 45.0 cm) with a white paper-lined floor and a dark goal box at one end.Rat paws were painted with a thin layer of acrylic paint to visualize and record their footprints as they walked down the track.The SFI was calculated as follows:

    In this formula, PL stands for print length and means the distance between the 3rd toe tip and the heel;TS stands for toe spread and means the distance between the 1st toe and 5th toe; and IT stands for intermediary toe spread and means the distance from the 2nd toe to the 4th toe.In addition, E means the distance measured from the operation-side of the paw, and N means the distance measured from the contralateral side of the paw.For example,EPLmeans the PL measured from the operation-side of the paw.SFI values could range from 0 (good recovery) to -100 (total dysfunction).

    2.21.1.Target muscle histology

    The operated hind limb was harvested four weeks post-surgery,and the gastrocnemius muscles extracted for weighing and normalization to the muscle on the unoperated side.Transverse sections (10 mm thick) of the isolated muscles were stained with the Masson trichrome stain.Six random fields were imaged across the muscle fiber sections under brightfield microscopy and the images analyzed with ImageJ.

    2.21.2.Electrophysiological assays

    Studies were performed at four weeks post-surgery as detailed in earlier studies[21].Rats were anesthetized with isoflurane(RWD Life Science Co., Ltd., Shenzhen, China) and the exposed sciatic nerve was clamped between two electrode placed 2 cm apart.Compound muscle action potentials (CMAP) from the gastrocnemius belly on the ipsilateral side were recorded using a multichannel electrophysiological recorder.The peak amplitude of the CMAP, nerve conduction velocity, and the latency of CMAP onset were quantified and compared across groups.

    2.22.Statistical analysis

    Fig.1.Sciatic nerve injury increases superoxide production.(A,B)Representative fluorescence images of dihydroethidium(DHE)staining of superoxide in sciatic nerves 2 h(A)and 12 h(B)following sciatic nerve injury(SNI),with and without N-acetylcysteine(NAC)treatment.(C)Representative fluorescence images of DHE staining in dorsal root ganglia cells(DRGs)12 h after H2O2 treatment.(D-F) Quantification of DHE fluorescence for experiments in Figs.1A-C, respectively (n = 6 per group).**P <0.01 and ***P <0.001 vs.sham or vehicle; ##P <0.01, and ###P <0.001 vs.SNI or H2O2 treatment.DPI: diphenyleneiodonium.

    All statistical analyses were performed using GraphPad Prism(Version 8.0,GraphPad)and the data expressed as mean±standard error of the mean.Data were analyzed and compared using oneway analysis of variance, with Tukey's post hoc tests.P <0.05 was designated as the threshold for statistical significance.

    3.Results

    3.1.ROS increase at the injury site and in DRG cell bodies after PNI

    Since ROS and redox stress are major contributors to cell death,we speculated that this may also occur in damaged DRG neurons and leads to neuronal cell death.Thus, we first evaluated ROS production and redox stress by detecting the expression of oxidized DHE in ex vivo DRG and sciatic nerve cultures.DHE fluorescence was significantly increased around the injury site as early as 2 h following SNI (Figs.1A and D) and could be detected in DRG 12 h after injury (Figs.1B and E), indicating that ROS levels in DRG neurons significantly increased after SNI.This injury-induced increase in ROS can be blocked by the oxidant scavenger,NAC,in both the sciatic nerve (Figs.1A and D) and the DRG (Figs.1B and E).In addition, injury-induced redox stress can also be demonstrated with H2O2treatment of uninjured sciatic nerve, which similarly increased DHE fluorescence and can be blocked by NAC (Figs.1C and F).Further supports for these findings can be seen in increased staining with the pan-ROS probe, hydrocyanine-Cy3, in SNI samples(Fig.S1),and the non-selective inhibitor of ROS-producing flavoenzymes, DPI, was as effective as NAC at reducing hydrocyanine-Cy3 staining.

    3.2.Ferroptosis is an important type of cell death in DRGs after H2O2 treatment

    Several methods have been reported to induce DRGs damage ranging from mechanical stimulation [22] to chemical stimulation[23].Our results showed that the level of ROS mounted up at the injury site, especially at DRG cell body after PNI.Compared with hypoxia and gamma rays,H2O2was applied in our study to trigger cellular oxidative stress in DRGs for its good stability, practical basis, and wide application in previous reports.In addition, H2O2treatment produced similar oxidative responses as SNI in vivo.Therefore, we established a model of oxidative stress using H2O2in vitro for further studies.First, DRGs were treated with various concentrations of H2O2to investigate the relationship between H2O2treatment and ROS generation.H2O2treatment increased ROS production in DRGs,as shown by increased DCFH-DA fluorescence intensity, in a dose-dependent manner (Figs.2A and C), and dramatically increased cell death (40.3% in DRGs) at the 1 mM concentration (Figs.2B and D), suggesting that oxidative injury occurred in DRGs after H2O2treatment.To determine the type of cell death that occur following H2O2treatment,we found that H2O2in the millimole range substantially increased apoptosis (Fig.S2A)and autophagy(Figs.S2B-D)in DRGs.Interestingly,H2O2-induced cell death was partially, but not completely, attenuated by z-VADFMK,an inhibitor of apoptosis(Figs.2B and D).This suggested that the mode of H2O2-induced cell death was not limited to apoptosis,and other forms of cell death, such as ferroptosis, may also be involved.

    Fig.2.H2O2-induced cell death is associated with ferroptosis.(A, C) Images (A) and quantification (C) of 2,7-dichlorodihydrofluorescein diacetate (DCFH-DA) staining of reactive oxygen species (ROS) in dorsal root ganglia cells (DRGs) exposed to different concentration of H2O2.(B, D) Images (B) and quantification (D) of the percentage of dead DRGs (red)with H2O2 (1 mM) or H2O2 + z-VAD-FMK (10 μM) treatment.(E, F) DRGs treated with H2O2 (1 mM) or erastin (20 μM) were subjected to imaging (E) and quantification (F) of FerroOrange staining of intracellular iron levels.(G) Transmission electron microscopy (TEM) analysis of mitochondrial ultrastructure.(H, I) Imaging (H) and quantification (I) of C11-BODIPY581/591 staining, with green fluorescence indicating lipid peroxidation.(J) Malondialdehyde (MDA) levels determination.(K) Glutathione (GSH) content determination.(L, M) Western blot (L) and quantification (M) of SLC7A11 and glutathione peroxidase 4 (GPX4) protein levels.*P <0.05, **P <0.01, and ***P <0.001 vs.control; #P <0.05 and##P <0.01 vs.H2O2 treatment.

    To test the possible involvement of ferroptosis, we first measured the changes in intracellular iron levels in DRG cultures after H2O2treatment and observed a 12.71-fold increase(Figs.2E and F).Morphologically, H2O2-treated cells displayed shrunken mitochondria,ruptured outer mitochondrial membrane,and light vacuoles that are likely related to the collapsed mitochondria(Fig.2G).These features are characteristic of ferroptosis [24], and comparable to those induced by erastin treatment (Fig.2G), a known ferroptosis inducer.In addition, H2O2treatment also elevated the production of membrane lipid peroxides(Figs.2H-J).It has been well established that ferroptosis is typically dependent on reduced expression of GPX4 and GSH, leading to excess lipid peroxidation [25].We observed that both H2O2and erastin treatment significantly reduced GSH levels (Fig.2K) and GPX4 expression (Figs.2L and M) in DRGs.We next examined the expression of SLC7A11, which transports the GSH precursor cysteine into the cytosol [26], and found that it was significantly inhibited by erastin but not by H2O2(Figs.2L and M).Together,these results indicate that ferroptosis is involved in H2O2-induced cell death in DRGs.

    3.3.Inhibition of ferroptosis increases the number of DRGs after H2O2 treatment

    To confirm that ferroptosis does indeed participate in H2O2-induced cell death,we performed studies with PIH,an iron chelator,and Fer-1, a strong selective inhibitor of ferroptosis.Both PIH and Fer-1 treatment significantly reduced H2O2-induced cell death and LDH levels (Figs.3A-C).H2O2-induced increase in lipid peroxidation and MDA levels were also reduced by PIH and Fer-1 (Figs.3D-F).Further studies showed that PIH and Fer-1 significantly increased GPX4 expression in H2O2-treated cell cultures (Figs.3G and H).The H2O2-induced morphological changes in mitochondria were also inhibited by PIH and Fer-1 (Fig.3I), indicating that ferroptosis does contribute to H2O2-induced cell death.

    3.4.Fpn is associated with ferroptosis in DRGs

    ROS may induce ferroptosis by degrading heme to release free iron or regulate the transportation of free iron.Evaluating possible effects on heme levels,we found that SNI did not alter the levels of heme in the DRG, sciatic nerve, liver,spleen,and serum compared with sham surgery (Figs.4A and B).Hemoglobin levels were also unaltered between the SNI and sham groups (Fig.4C).We then evaluated the expression of iron-related genes in DRG, including TfR1,FtH,FtL,and Fpn,and found no significant difference between the SNI and sham groups at the mRNA level(Fig.4D).Interestingly,the level of Fpn protein,the only known mammalian nonheme iron export protein[27],was reduced in the DRG after SNI(Figs.4E and F), suggesting that ROS-related ferroptosis may involve Fpn.To verify this,we found that H2O2administration caused aggregation of free iron in the DRGs, which can be suppressed by Fpn overexpression (Figs.4G and I).TUNEL assay indicated that Fpn overexpression can inhibit H2O2-induced DRG damage(Figs.4H and J),while Western blot showed that the H2O2-induced reduction in GPX4 level can be rescued by Fpn overexpression (Figs.4K and L).These data support the involvement of Fpn in H2O2-mediated ferroptosis and warrants further exploration.

    3.5.Ferroptosis is induced by Fpn ubiquitination and degradation in DRGs

    Although the mRNA expression of Fpn did not change in the sciatic nerve and DRG following SNI (Fig.4D), significant decrease in Fpn proteins was observed(Figs.4E and F),suggesting that redox regulation of Fpn protein may be required for neuronal ferroptosis.Thus, we next investigated possible mechanisms by which H2O2may regulate Fpn in vitro,such as via protein degradation pathways.Pretreatment of DRGs with the protein synthesis inhibitor, cycloheximide(CHX),reduced Fpn protein levels over time,and addition of H2O2enhanced Fpn degradation(Figs.5A-C).Next,we explored the two major protein degradation mechanisms, the autophagylysosomal pathway and ubiquitin-proteasomal pathway, with bafilomycin A1 (lysosome inhibitor) and bortezomib (proteasome inhibitor) treatments, respectively.The H2O2-induced Fpn degradation can be blocked by bortezomib alone or combined with bafilomycin A1, but not by bafilomycin A1 alone (Figs.5D and E),suggesting that H2O2primarily degrades Fpn via the ubiquitinproteasomal pathway.

    To confirm that Fpn can interact with ubiquitin, DRGs were cotransfected with plasmids expressing Flag-Fpn and HA-ubiquitin,and whole-cell lysates were used for Co-IP with anti-Flag or anti-HA antibody conjugated agarose beads.The Co-IP data showed that Fpn does indeed interact with ubiquitin with anti-HA, anti-Fpn, and anti-Flag antibodies, all able to detect ubiquitinated species of Fpn in the DRG lysates (Figs.5F-H).Next, using the proteosome inhibitor,Z-Leu-Leu-Leu-al(MG132),to enhance ubiquitin accumulation, we found increased ubiquitination of proteins in H2O2-treated DRGs (Fig.5I).Co-IP with anti-Fpn showed that ubiquitination of Fpn is also increased by MG132 treatment,but the effect of H2O2is greater, while combined treatment with MG132 and H2O2synergistically enhanced Fpn ubiquitination(Figs.5J and K), suggesting polyubiquitination of Fpn.This interaction is likely specific since Co-IP with anti-IgG antibody did not reveal specific bands.Altogether,these results indicate that H2O2may reduce Fpn by enhancing the ubiquitination of this protein and directing it toward proteasomal degradation.

    3.6.HYD inhibits Fpn degradation in DRGs by suppressing UBA52 binding to Fpn

    To identify the specific ubiquitin protein that regulate Fpn ubiquitination and degradation, proteins from neuronal lysates were pulled down with the Fpn antibody for silver staining and liquid chromatography(LC)-MS,which identified a single ubiquitin protein, UBA52 (Figs.6A and B).UBA52 co-localized with Fpn in DRGs (Fig.6C), pulled-down with Fpn, and is increased with H2O2treatment along with overall ubiquitination, while Fpn level decreases (Fig.6D).To further explore the function of UBA52, we knocked down UBA52 with shRNA and found that the morphological changes induced by H2O2in primary neurons, including thinning and fragmentation of the neurites, formation of large vacuoles and bright spots, as well as decreased cytoplasm, can be rescued by shUBA52 (Fig.S3A).Ultrastructural characteristics of ferroptosis, such as the altered mitochondrial structure in H2O2-treated cells, was also inhibited by shUBA52 administration (Fig.S3B).UBA52 can interact with different signaling proteins via different functional domains.To determine which domain of UBA52 interacts with Fpn, we generated three Myc-tagged UBA52 domains (amino acids (aa) 41-117, aa 120-167, and the whole aa 1-169;Fig.6E,top panel)to perform Co-IP with Fpn.Only domain(aa 41-117) showed specific interaction with Fpn (Fig.6E bottom panel).Interestingly, the UBA52 domain (aa 41-117) is conserved across species (Fig.S4).

    Fig.3.Iron chelation and ferrostatin-1(Fer-1)treatment inhibit H2O2-induced ferroptosis.(A,B)Dorsal root ganglia cells(DRGs)treated with H2O2(1 mM)alone or with pyridoxal isonicotinoyl hydrozone (PIH) (1 mM) or Fer-1 (6 μM) for 24 h were subjected to imaging (A) and quantification (B) of dead DRGs (red).(C) Lactate dehydrogenase (LDH) level determination.(D, E) Imaging (D) and quantification (E) of C11-BODIPY581/591 staining, with green fluorescence indicating lipid peroxidation.(F) Malondialdehyde (MDA) level determination.(G, H) Western blot (G) and quantification (H) of glutathione peroxidase 4 (GPX4) protein expression.(I) Transmission electron microscopy (TEM) analysis of mitochondrial morphology.**P <0.01 and ***P <0.001 vs.control; #P <0.05, ##P <0.01, and ###P <0.001 vs.H2O2 treatment.

    Since degradation of Fpn may contribute to neuronal death after injury, a potential therapeutic strategy is to reduce/prevent Fpn degradation, and disrupting UBA52-Fpn binding may achieve this.Thus, we interrogated the Drug Signatures Database (DSigDB),which contains detailed information on existing drugs,their targets and activity data,to find compounds that may disrupt the UBA52-Fpn binding.Based on the structure of UBA52, nine compounds were predicted,of which two,HYD and deltamethrin,were further analyzed.Utilizing the Protein Data Bank database, we evaluated the predicted protein structure of the UBA52 and calculated the energy required for interaction with UBA52.Both HYD and deltamethrin possessed sufficiently high binding energies to interact with UBA52(Figs.6G and H).However,deltamethrin is a common insecticide that has demonstrated neurological,cardiovascular,and reproductive toxicity in various aquatic organisms [28].Therefore,HYD was selected for subsequent biological testing.Conjugating HYD to biotin allowed us to perform an affinity pull-down assay(Fig.6I) and showed that HYD-biotin can bind to endogenous UBA52 in DRGs (Fig.6J).To examine whether HYD affect H2O2-induced Fpn degradation, we performed Co-IP assays after treatment with HYD.HYD treatment increased Fpn levels with or without H2O2stimulation and reduced the binding of UBA52 and ubiquitin to Fpn.However, HYD had little effect on UBA52 levels(Fig.6K).

    Fig.4.Sciatic nerve injury(SNI)promotes rapid and systemic iron overload independent of heme degradation.(A,B)Four days after SNI,heme concentrations were determined for various tissues (A) and serum (B).(C) Hemoglobin (HGB) levels were measured.(D-F) Expression of iron-related genes (D) and proteins (E, F) were measured.(G, I) Dorsal root ganglia cells(DRGs)transduced with lentivirus(Lv)overexpressing ferroportin(Fpn)and treated with H2O2 were subjected to imaging(G)and quantification(I)of intracellular iron levels.(H, J) imaging (H) and quantification (J) of terminal deoxynucleotidyl transferase dUTP nick-end labeling (TUNEL) staining.(K, L) Western blot(K) and quantification (L) of glutathione peroxidase 4(GPX4)protein levels.**P <0.01 and ***P <0.001 vs.sham or vector; ##P <0.01 and ###P <0.001 vs.H2O2+vector treatment.TfR1:transferrin receptor 1;FtH: ferritin H; FtL: ferritin L; DAPI: 4′,6-diamidino-2-phenylindole.

    As HYD can reverse the H2O2-induced Fpn degradation, it is likely to inhibit ferroptosis in H2O2-treated cells.To confirm this,DRGs were treated with H2O2or H2O2+HYD for 24 h and iron level,cell death, lipid peroxidation, and GPX4 levels were evaluated.H2O2-induced increases in iron production (Figs.7A and B), cell death (Figs.7C and D),and lipid peroxidation (Figs.7E and F), and decrease in GPX4 expression(Figs.7G and H)were all rescued with HYD treatment.Findings of the in vivo studies are also consistent with the above in vitro data.Significantly higher levels of iron were observed in the SNI group compared to the sham group, and HYD can reduce this iron accumulation in DRGs after SNI (Fig.S5A).In addition,Western blot and IHC results indicated that HYD can also rescue Fpn expression after SNI (Figs.S5B and C).These results revealed that HYD treatment significantly inhibits DRGs ferroptosis induced by ROS.

    3.7.HYD increases sciatic nerve regeneration by blocking ferroptosis

    The ability of HYD to rescue Fpn expression in vitro and in vivo suggest that it may have therapeutic potential.To verify this, HYD was administered to rats immediately or one week (delayed administration)after SNI(Fig.8A),followed by evaluation of axonal regeneration four weeks after injury.The length(from crush site)of axons positive for GAP43, a marker of axonal regeneration, was longest in the immediate HYD treatment group, followed by the delayed HYD treatment group and the untreated control group;differences among all groups were significant (Fig.8B), and Western blot quantifications confirmed these results (Figs.8C and D).These results demonstrate that inhibiting ferroptosis improves axonal regeneration following SNI.

    Gait is one of the primary indexes for evaluating sciatic nerve function and recovery.At two weeks after SNI, no significant difference was found between rats with and without HYD treatment,indicating a lack of functional recovery at this time-point across all groups.However, by four weeks after SNI, rats who had received HYD treatment showed significantly fuller footprints while walking(Fig.8E) and higher SFI values (Fig.8F), indicating that the inhibition of ferroptosis can promote functional neurological recovery following SNI.In addition, immediate HYD treatment produced better outcomes than delayed HYD treatment (Figs.8E and F),suggesting that inhibiting ferroptosis immediately after injury is more effective in promoting neurological recovery.

    Fig.5.H2O2 promotes ubiquitination-based degradation of ferroportin(Fpn).(A-C)Dorsal root ganglia cells(DRGs)pretreated with cycloheximide(CHX)(5 μg/mL)for 24 h were stimulated with H2O2 for different length of time and subject to Western blot analysis of Fpn expression (A), quantification of Fpn protein levels over time (B), and calculation of overall changes in Fpn expression (C).(D, E) Western blot (D) and quantification (E) of Fpn protein levels after treatment with CHX (5 μg/mL), bortezomib (Bort; 200 nM), bafilomycin A1(Baf A1;50 nM),or a combination of Bort and Baf A1 for 1 h before the introduction of 1 mM H2O2.(F-H)Co-immunoprecipitation(Co-IP)of Fpn-Flag and ubiquitin-HA with Western blot analysis of HA(F),Fpn(G),and Flag(H),depicting ubiquitinated species of Fpn.DRGs were exposed to Z-Leu-Leu-Leu-al(MG132;20 μM)for 6 h prior to exposure to H2O2(1 mM)for 24 h.(I) Total lysates were analyzed for ubiquitinated proteins via immunoblotting.(J) Co-IP with Fpn or IgG antibody and Western blot with anti-ubiquitin in MG132-treated and untreated cells.(K)Quantification of Western blot images from Fig.5J.*P <0.05, **P <0.01,and ***P <0.001 vs.control or MG132; ##P <0.01 vs.H2O2 treatment.GAPDH: glyceraldehyde-3-phosphate dehydrogenase; DMSO: dissolved in dimethyl sulfoxide; Ub-HA: HA labeled ubiquitin.

    The degree of gastrocnemius atrophy is also widely used to estimate the recovery of sciatic function.Masson trichrome staining of the gastrocnemius (Fig.8G) at four weeks after injury showed that HYD treatment significantly increased the average muscle fiber area(Fig.8H)and the relative wet weight ratio of the gastrocnemius (Fig.8I), with the effect of immediate treatment superior to delayed treatment.Hence, immediate inhibition of post-injury ferroptosis by HYD appeared to reduce target muscle atrophy.

    Fig.6.UBA52 is crucial for H2O2-induced ferroportin(Fpn)ubiquitination.(A)Dorsal root ganglia(DRG)cell lysates immunoprecipitated with Fpn antibody and treated with silver staining to reveal all proteins bound to Fpn antibody.(B)Mass spectrogram of UBA52.(C)Immunofluorescence staining showing co-localization of UBA52(green)and Fpn(red)in primary cortical neurons.(D) Immunoprecipitation (IP) with Fpn followed by Western blot with anti-ubiquitin, anti-Fpn and anti-UBA52 in dorsal root ganglia cells (DRGs).(E)Schematic representation of UBA52 fusion proteins (top panel).Co-IP showing interaction between UBA52 domain (amino acides (aa) 41-117) and Fpn (bottom panel).(F)Schematic of the experiment that predicts the docking with small molecule drugs.(G,H)Predication of binding between UBA52 and hydralazine(HYD)(G)or deltamethrin(H)as determined by molecular docking.(I) Schematic of the experiment that detects the binding of HYD-biotin to UBA52.(J) Pull-down assay to detect the binding of HYD-biotin to UBA52 and immunoblotting analysis of the input proteins.(K)IP with Fpn followed by Western blot with anti-ubiquitin shows that the addition of HYD reduces ubiquitination of Fpn in the presence or absence of H2O2.M:marker;DSigDB:Drug Signatures Database;GAPDH:glyceraldehyde-3-phosphate dehydrogenase;Myc:protein tag(protein sequence:EQKLISEEDL); IB: immunoblotting; SDS-PAGE: sodium dodecyl-sulfate polyacrylamide gel electrophoresis; Ub: ubiquitin.

    Fig.7.Hydralazine(HYD)inhibits H2O2-induced dorsal root ganglia cells(DRGs)death.(A,B)DRGs treated with H2O2(1 mM)with and without HYD for 24 h were subjected to imaging(A)and quantification(B)of intracellular iron levels.(C,D)imaging(C)and quantification (D) of dead cells (red).(E, F) Imaging (E) and quantification (F) of C11-BODIPY581/591 staining, with green fluorescence indicating lipid peroxidation.(G, H)Western blot (G) and quantification (H) of glutathione peroxidase 4 (GPX4) level.**P <0.01 and ***P <0.001 vs.control; #P <0.05 and ##P <0.01 vs.H2O2 treatment.

    We also evaluated functional recovery with electrophysiology.HYD treatment significantly increased the peak CMAP amplitude(Fig.8J) and nerve conduction velocity (Fig.8K), and reduced the CMAP latency (Fig.8L).In all cases, immediate HYD treatment resulted in better outcomes than delayed treatment.These data suggest that the inhibition of ferroptosis can promote functional motor recovery after SNI, with a more profound effect obtained with immediate HYD treatment.

    4.Discussion

    Excess levels of ROS can cause lipid, protein, and DNA damage,augmenting disease pathogenesis [29].Cell and organelle membranes are particularly vulnerable to ROS damage because they contain high levels of polyunsaturated fatty acids.In addition,peroxidation of membrane lipids not only damages cell components but also causes structural damage, finally resulting in cell rupture and death.In previous studies,Cox4i2 triggers an increase in reactive oxygen species, leading to ferroptosis and apoptosis in HHV7infected schwann cells, which is a glia cell in the peripheral nerves [30].Furthermore, knockdown of the lysosomal protein prosaposin strongly sensitizes neurons to oxidative stress by triggering the formation of lipofuscin, which traps iron, generating reactive oxygen species and triggering ferroptosis in central nervous system [31].In our study, ROS levels were found to be increased in both the sciatic nerve and DRG following SNI, in line with earlier reports[23].

    Previous studies have demonstrated that ROS production after spinal nerve transection mostly occur in the mitochondria of damaged neurons [32].Other studies have found that ROS is notably produced by macrophages and transported to the cell body via axons [33].Administration of H2O2onto the sciatic nerve significantly increases ROS in DRG.Hence,we speculated that ROS is produced by cells at the site of injury and delivered to the cell body by the axon.This conclusion is supported by the observation that administration of NAC at the injured site of the sciatic nerve can inhibit ROS production in DRG.Under normal conditions,GPX4 rapidly reduces lipid peroxides via GSH oxidation.However, GSH down-regulation owing to excess ROS accumulation leads to an increase in lipid peroxides, triggering cell death.In addition, we found that ROS increases iron accumulation, which can in turn damage the lipid and mitochondrial membranes.

    Iron is an essential cofactor for many cellular enzymes involved in myelin synthesis.Thus, imbalance in iron homeostasis will contribute to peripheral neuropathies.Recent studies showed that the accumulation of iron can increase ROS levels [34], and the expression of Fpn, which is involved in maintaining iron homeostasis,is decreased in the peripheral nerve system after SNI[35].In our study,SNI significantly increased the ROS level in DRGs,and the decreased Fpn expressions is associated with increased iron in the DRGs.However, the expression of Fpn was significantly decreased after H2O2treatment, indicating that ROS was unlikely to be induced by the accumulation of iron after SNI.

    We found that autophagy was significantly increased after H2O2treatment.As the autophagy-lysosomal pathway is one of the major mechanisms of protein degradation, we thus evaluated whether this pathway is involved in the degradation of Fpn after H2O2treatment.However, our data showed that ROS did not use the autophagy-lysosomal pathway for Fpn degradation.The other major protein degradation pathway involving proteosomes has been shown to degrade several plasma membrane metal transporters,including DMT1,SLC39A14,and SLC30A1[36].Fpn has also been shown to be degraded by this pathway following hepcidininduced endocytosis and ubiquitination [37].Similarly, we provide evidence in this study that ROS induced Fpn degradation principally via the ubiquitin-proteasomal pathway.

    The ubiquitin-proteasome system is essential in regulating mitochondrial physiology and controls the expression of different intramembrane space and inner mitochondria membrane proteins,such as succinate dehydrogenase subunit A, and other tricarboxylic acid-cycle metabolites[38].Recently,we reported that the ubiquitin protein, UBA52, facilitates the ubiquitination of chaperone HSP90 together with the E2 and E3 ligases, leading to inhibition of endoplasmic reticulum stress and Parkinson's diseasespecific pathological markers, such as tyrosine hydroxylase and αsynuclein[39].Our current study implicates UBA52 in SNI-induced ferroptosis as inhibition of UBA52 suppressed ferroptosis after SNI.We provided evidence of direct interaction between UBA52 and Fpn following H2O2treatment, indicating that UBA52 can effectively promote Fpn degradation after H2O2administration.This prompted us to explore potential targeting drugs to this pathway to mitigate SNI.

    Protein crystal structures not only highlights the potential mechanisms of action of specific proteins, but also offer insights into drugs and molecules that can modulate protein activity.Through such structural analysis of UBA52,we identified HYD and deltamethrin as potential molecules that can regulate the UBA52-Fpn interaction.Since deltamethrin is potentially toxic, we evaluated the function of HYD and found that it indeed has the potential to regulate Fpn ubiquitination by UBA52.Previous, HYD has been shown to support neuronal regeneration in a mouse model of spinal cord injury [40].However, it has not been demonstrated in protein degradation and peripheral nerve injury.Our results suggested that HYD can reduce Fpn degradation in DRG following injury, thus preventing accumulation of iron and subsequent ferroptosis.

    Recovery after PNI not only relies on enhanced axonal regeneration but also on reduced cell death, including reduced ferroptosis.Our study showed that HYD not only reduces sciatic nerve necrosis but also improves functional recovery following SNI.Previous studies found that primary sensory neurons started dying within 24 h of sciatic nerve transection[41],but neuronal loss only starts between 4 and 10 days following injury [42].This gap between the onset of neuronal death and neuronal loss reflects the time taken for neurons to move from the DNA fragmentation stage to complete involution.This gap may also be the key stage for preventing cell death in neurons.Here,we showed that immediate treatment with HYD following injury were far more effective in promoting axonal regeneration and functional recovery than delayed HYD treatment, suggesting that the early inhibition of ferroptosis is vital for axonal regeneration since ferroptosis begins 0-7 days following SNI.

    As a specific form of regulated cell death, ferroptosis plays crucial roles in multiple pathological conditions including ischemic organ injuries, cancer, and degenerative diseases [43];increasing evidence also implicates potential physiological functions of ferroptosis [44].Numerous organ injuries and degenerative pathologies are driven by ferroptosis [45].Multiple studies have established a role of ferroptosis in neuronal demise,including stroke and other brain injuries.An ex vivo study using rat hippocampal slice culture showed that glutamate-induced neuronal excitotoxic cell death can be blocked by radicaltrapping antioxidant ferrostatin 1 [45].As such, pharmacological modulation of ferroptosis, via both its induction and its inhibition, holds great potential for the treatment of drug-resistant organ injuries, cancers, and other degenerative diseases linked to extensive lipid peroxidation.Here, we demonstrated occurrence of ferroptosis in SNI arising from ROS-accelerated UBA52 ubiquitination and degradation of Fpn, leading to accumulation of intracellular iron and stimulation of lipid peroxidation.These findings suggest that UBA52 is a potential therapeutic target for peripheral nerve injury.Furthermore, we provided preliminarily evaluation of HYD, an U.S.Food and Drug Administrationapproved anti-hypertensive agent, as a potential treatment compound to inhibit ferroptosis in DRG by targeting UBA52.Clinical translation of HYD for restoration of peripheral nerve injury still requires further exploration and verification.Our study demonstrated a promising strategy of screening approved drug libraries for new indications.

    5.Conclusion

    In this study, ROS levels were found to be increased in both the sciatic nerve and DRG following SNI.Under normal conditions, GPX4 rapidly reduces lipid peroxides via GSH oxidation.However, GSH down-regulation owing to excess ROS accumulation leads to an increase in lipid peroxides,triggering ferroptosis.In addition, we found that ROS increases iron accumulation,which can in turn damage the lipid and mitochondrial membranes.In addition, the expression of Fpn was significantly decreased after H2O2treatment via the ubiquitin-proteasomal pathway.Furthermore, we found UBA52 can effectively promote Fpn degradation after H2O2administration.This prompted us to explore potential targeting drugs to this pathway to mitigate SNI.

    We found HYD can reduce Fpn degradation by targeting UBA52 in DRG following injury,thus preventing accumulation of iron and subsequent ferroptosis.Recovery after PNI not only relies on enhanced axonal regeneration but also on reduced cell death,including reduced ferroptosis.Our study showed that HYD not only reduces sciatic nerve necrosis but also improves functional recovery following SNI.In addition, we showed that immediate treatment with HYD following injury were far more effective in promoting axonal regeneration and functional recovery than delayed HYD treatment, suggesting that the early inhibition of ferroptosis is vital for axonal regeneration since ferroptosis begins 0-7 days following SNI.

    CRediT author statement

    Shengyou Li:Conceptualization,Data curation,Formal analysis,Investigation, Writing - Original draft preparation;Xue Gao:Data curation, Formal analysis, Writing - Original draft preparation;Yi Zheng:Data curation,Formal analysis,Validation;Yujie Yang:Data curation, Investigation;Jianbo Gao:Conceptualization, Supervision; Methodology;Dan GengandLingli Guo:Data curation,Formal analysis;Teng Ma:Data curation, Funding acquisition,Formal analysis;Yiming Hao:Conceptualization, Methodology;Bin Wei:Supervision, Formal analysis, Investigation;Liangliang Huang:Data curation, Funding acquisition;Yitao Wei:Formal analysis, Investigation;Bing Xia:Conceptualization, Funding acquisition, Writing - Reviewing and Editing;Zhuojing Luo:Project administration, Funding acquisition, Writing - Reviewing and Editing;Jinghui Huang:Conceptualization, Funding acquisition,Supervision, Project administration, Writing - Reviewing and Editing.

    Declaration of competing interest

    The authors declare that there are no conflicts of interest.

    Acknowledgments

    This work was supported by grants from the National Natural Science Foundation of China (Grant Nos.: 82122043, 81972052,81902213, 82201537, and 81730065), and the China Postdoctoral Science Foundation (Grant Nos.: 2021M693946 and 2019M653967).

    Appendix A.Supplementary data

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.jpha.2023.08.006.

    亚洲av国产av综合av卡| 国产成人免费无遮挡视频| 人人澡人人妻人| 免费一级毛片在线播放高清视频 | 99精国产麻豆久久婷婷| 精品国产一区二区久久| 国产免费视频播放在线视频| 久久久久久久大尺度免费视频| 一进一出好大好爽视频| 久久久国产成人免费| 天天添夜夜摸| 亚洲九九香蕉| 亚洲中文av在线| 久久免费观看电影| 欧美中文综合在线视频| 黑人巨大精品欧美一区二区mp4| 欧美性长视频在线观看| 在线十欧美十亚洲十日本专区| 亚洲欧美一区二区三区久久| 涩涩av久久男人的天堂| 肉色欧美久久久久久久蜜桃| 大片电影免费在线观看免费| 在线十欧美十亚洲十日本专区| 丁香六月欧美| 真人做人爱边吃奶动态| 王馨瑶露胸无遮挡在线观看| 国产日韩欧美亚洲二区| 人人妻人人澡人人爽人人夜夜| 18在线观看网站| 亚洲综合色网址| 亚洲专区字幕在线| 久久久久网色| 考比视频在线观看| a级毛片在线看网站| 精品午夜福利视频在线观看一区 | 国产在线视频一区二区| 亚洲人成伊人成综合网2020| 一本色道久久久久久精品综合| 丰满少妇做爰视频| 欧美在线黄色| 亚洲国产欧美一区二区综合| 久久久久国内视频| 国产欧美日韩一区二区精品| 久久久国产成人免费| 美女主播在线视频| 看免费av毛片| 成人国产一区最新在线观看| 精品人妻1区二区| 蜜桃国产av成人99| 亚洲成av片中文字幕在线观看| 国产91精品成人一区二区三区 | 国产一区二区三区视频了| 色老头精品视频在线观看| 亚洲av片天天在线观看| 在线观看舔阴道视频| 啦啦啦中文免费视频观看日本| 精品一区二区三区视频在线观看免费 | 午夜老司机福利片| 首页视频小说图片口味搜索| 免费久久久久久久精品成人欧美视频| 久久国产精品影院| 欧美亚洲日本最大视频资源| 精品少妇黑人巨大在线播放| 久久久精品国产亚洲av高清涩受| 亚洲av成人一区二区三| 岛国在线观看网站| 国产日韩欧美视频二区| 男女午夜视频在线观看| 亚洲熟女毛片儿| 国产精品麻豆人妻色哟哟久久| 久久久久久人人人人人| 美女视频免费永久观看网站| 少妇 在线观看| 在线观看人妻少妇| 丝袜喷水一区| 精品一品国产午夜福利视频| 国产日韩欧美在线精品| 国产免费av片在线观看野外av| 99精品在免费线老司机午夜| 久久精品亚洲av国产电影网| 欧美黄色淫秽网站| 久久人妻熟女aⅴ| 十八禁人妻一区二区| 久久久欧美国产精品| 国产免费现黄频在线看| 下体分泌物呈黄色| 蜜桃在线观看..| 国产av一区二区精品久久| 视频区欧美日本亚洲| 丰满饥渴人妻一区二区三| 欧美+亚洲+日韩+国产| 久久ye,这里只有精品| 99国产精品99久久久久| 91av网站免费观看| 黄色 视频免费看| 热re99久久精品国产66热6| 亚洲美女黄片视频| 我要看黄色一级片免费的| 婷婷成人精品国产| 精品一区二区三卡| 日韩熟女老妇一区二区性免费视频| 一级片免费观看大全| 国产精品秋霞免费鲁丝片| 久久影院123| 亚洲va日本ⅴa欧美va伊人久久| 日韩有码中文字幕| 国产免费现黄频在线看| 亚洲精品国产一区二区精华液| 一二三四在线观看免费中文在| 久久久欧美国产精品| 后天国语完整版免费观看| 亚洲色图 男人天堂 中文字幕| 亚洲 欧美一区二区三区| 好男人电影高清在线观看| 亚洲少妇的诱惑av| 亚洲七黄色美女视频| 国产精品 国内视频| 中文字幕最新亚洲高清| 精品亚洲乱码少妇综合久久| av网站免费在线观看视频| av天堂久久9| 亚洲精品中文字幕在线视频| 精品国产乱码久久久久久男人| 欧美另类亚洲清纯唯美| 成人国产av品久久久| 欧美一级毛片孕妇| 成年人免费黄色播放视频| 国产精品.久久久| 最新的欧美精品一区二区| 欧美黄色淫秽网站| 交换朋友夫妻互换小说| 极品人妻少妇av视频| 久久影院123| 老熟妇乱子伦视频在线观看| 黄色丝袜av网址大全| 91国产中文字幕| 老司机福利观看| 亚洲一区二区三区欧美精品| 久久热在线av| 国产精品影院久久| 国产精品国产av在线观看| 欧美乱妇无乱码| 精品国产一区二区久久| 国产视频一区二区在线看| 久久久久久免费高清国产稀缺| 国产精品.久久久| 精品久久蜜臀av无| 中文字幕制服av| 免费女性裸体啪啪无遮挡网站| 亚洲精品国产一区二区精华液| 黄网站色视频无遮挡免费观看| 日韩三级视频一区二区三区| 纵有疾风起免费观看全集完整版| 9色porny在线观看| 国产日韩欧美在线精品| 成年人黄色毛片网站| netflix在线观看网站| 男女下面进入的视频免费午夜| 18禁国产床啪视频网站| 99久久精品国产亚洲精品| 欧美激情久久久久久爽电影| 亚洲精品在线美女| 最近最新中文字幕大全免费视频| 真实男女啪啪啪动态图| 成人精品一区二区免费| 欧美午夜高清在线| 日本黄色片子视频| 精品久久久久久久久久免费视频| 亚洲自偷自拍图片 自拍| 国产伦精品一区二区三区视频9 | 成年女人毛片免费观看观看9| 国产精品一区二区精品视频观看| 欧美日韩中文字幕国产精品一区二区三区| 日本免费a在线| 色尼玛亚洲综合影院| 黄色 视频免费看| 免费观看人在逋| 中文在线观看免费www的网站| 天堂av国产一区二区熟女人妻| 久久久久国产精品人妻aⅴ院| 国产91精品成人一区二区三区| 午夜激情福利司机影院| 国产精品一区二区免费欧美| 村上凉子中文字幕在线| 麻豆成人午夜福利视频| 日韩欧美在线二视频| 黄色 视频免费看| 天天添夜夜摸| 美女 人体艺术 gogo| 国产三级在线视频| av欧美777| 99国产精品一区二区蜜桃av| 伦理电影免费视频| 亚洲专区国产一区二区| h日本视频在线播放| 久久精品国产99精品国产亚洲性色| 亚洲片人在线观看| 欧美性猛交黑人性爽| 夜夜爽天天搞| 熟妇人妻久久中文字幕3abv| 国产黄色小视频在线观看| 一区福利在线观看| 91在线精品国自产拍蜜月 | av女优亚洲男人天堂 | 99热只有精品国产| 18禁裸乳无遮挡免费网站照片| 19禁男女啪啪无遮挡网站| 后天国语完整版免费观看| 亚洲人成网站高清观看| 黄色片一级片一级黄色片| 日韩大尺度精品在线看网址| 国产激情偷乱视频一区二区| 久久精品91蜜桃| 听说在线观看完整版免费高清| 天天躁日日操中文字幕| 久久精品国产99精品国产亚洲性色| 久久午夜综合久久蜜桃| 亚洲自拍偷在线| 午夜免费激情av| 国产成+人综合+亚洲专区| 亚洲av美国av| 特大巨黑吊av在线直播| 免费一级毛片在线播放高清视频| 夜夜看夜夜爽夜夜摸| 午夜免费成人在线视频| 毛片女人毛片| 一卡2卡三卡四卡精品乱码亚洲| 最新美女视频免费是黄的| 国产av在哪里看| 禁无遮挡网站| 国产av一区在线观看免费| 97碰自拍视频| 精品欧美国产一区二区三| 狂野欧美激情性xxxx| 老汉色∧v一级毛片| 日日干狠狠操夜夜爽| 丁香六月欧美| 久久精品国产清高在天天线| 亚洲精品456在线播放app | 国产v大片淫在线免费观看| 国产高清videossex| 国内少妇人妻偷人精品xxx网站 | 噜噜噜噜噜久久久久久91| 免费大片18禁| 真人一进一出gif抽搐免费| 一个人看视频在线观看www免费 | 在线免费观看不下载黄p国产 | 成人三级黄色视频| 国产精品久久视频播放| 亚洲一区二区三区不卡视频| 女生性感内裤真人,穿戴方法视频| 欧美中文综合在线视频| 亚洲国产精品sss在线观看| 男女床上黄色一级片免费看| 中文字幕最新亚洲高清| 婷婷六月久久综合丁香| 成熟少妇高潮喷水视频| 亚洲自拍偷在线| 成人性生交大片免费视频hd| 亚洲在线自拍视频| 久久久成人免费电影| 亚洲av日韩精品久久久久久密| 国产高清videossex| 亚洲国产精品sss在线观看| 亚洲国产欧美人成| 母亲3免费完整高清在线观看| 久久久久久久久中文| 88av欧美| 国产人伦9x9x在线观看| 一a级毛片在线观看| xxxwww97欧美| 国产伦精品一区二区三区四那| 麻豆成人av在线观看| 在线观看美女被高潮喷水网站 | 俺也久久电影网| 国产亚洲av高清不卡| 舔av片在线| 久久久久九九精品影院| 老司机深夜福利视频在线观看| 国产亚洲欧美在线一区二区| 人妻丰满熟妇av一区二区三区| 久久婷婷人人爽人人干人人爱| 成在线人永久免费视频| 午夜日韩欧美国产| 综合色av麻豆| 操出白浆在线播放| 国产高清视频在线播放一区| 久久久国产成人免费| 99视频精品全部免费 在线 | 国内少妇人妻偷人精品xxx网站 | 巨乳人妻的诱惑在线观看| 亚洲成人久久爱视频| 黄色丝袜av网址大全| 欧美在线黄色| www国产在线视频色| 日韩欧美国产一区二区入口| 美女 人体艺术 gogo| 黄色 视频免费看| 国产成人影院久久av| 亚洲专区国产一区二区| 精品久久久久久久毛片微露脸| 欧美国产日韩亚洲一区| 久9热在线精品视频| 久久国产精品人妻蜜桃| 久久国产精品影院| 国产91精品成人一区二区三区| 真人一进一出gif抽搐免费| 免费在线观看日本一区| 国产精品久久久久久亚洲av鲁大| 网址你懂的国产日韩在线| 国产亚洲精品一区二区www| 久久国产乱子伦精品免费另类| 国产精品久久视频播放| 人妻丰满熟妇av一区二区三区| 亚洲 国产 在线| 欧美日韩乱码在线| 天堂网av新在线| svipshipincom国产片| 国产精品99久久99久久久不卡| 久久精品aⅴ一区二区三区四区| 亚洲人成电影免费在线| 午夜两性在线视频| 一个人免费在线观看的高清视频| 一级作爱视频免费观看| 国产探花在线观看一区二区| 国产探花在线观看一区二区| 国产欧美日韩精品一区二区| 欧美色欧美亚洲另类二区| 欧美日韩乱码在线| 久久久久精品国产欧美久久久| 欧美日韩综合久久久久久 | 一进一出好大好爽视频| 男人和女人高潮做爰伦理| 亚洲av电影在线进入| 亚洲av电影不卡..在线观看| 欧美黄色淫秽网站| 九九在线视频观看精品| 又紧又爽又黄一区二区| 国产不卡一卡二| 日本成人三级电影网站| 欧美大码av| 人妻夜夜爽99麻豆av| av中文乱码字幕在线| 极品教师在线免费播放| 久久久久久久久免费视频了| 欧美在线黄色| 亚洲人成电影免费在线| 国产精品亚洲美女久久久| 国产高清videossex| 男人舔奶头视频| 欧美不卡视频在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 国产私拍福利视频在线观看| 啦啦啦免费观看视频1| 亚洲精品久久国产高清桃花| 男女视频在线观看网站免费| 两个人视频免费观看高清| 亚洲狠狠婷婷综合久久图片| 亚洲欧洲精品一区二区精品久久久| av中文乱码字幕在线| 男插女下体视频免费在线播放| 日韩高清综合在线| 色综合婷婷激情| 亚洲国产欧美人成| 别揉我奶头~嗯~啊~动态视频| 人妻夜夜爽99麻豆av| 久久精品影院6| 中文字幕久久专区| 欧美日韩黄片免| 一本一本综合久久| 欧美大码av| 亚洲激情在线av| 精品午夜福利视频在线观看一区| 男人舔女人下体高潮全视频| 欧美中文综合在线视频| 国产精品av视频在线免费观看| 国产精品av久久久久免费| 欧美高清成人免费视频www| 亚洲成人免费电影在线观看| 18禁裸乳无遮挡免费网站照片| 一本精品99久久精品77| 夜夜看夜夜爽夜夜摸| 久久婷婷人人爽人人干人人爱| 91久久精品国产一区二区成人 | 无限看片的www在线观看| 国内毛片毛片毛片毛片毛片| 日本熟妇午夜| 亚洲美女视频黄频| 国产精华一区二区三区| 老熟妇仑乱视频hdxx| 操出白浆在线播放| 2021天堂中文幕一二区在线观| 91av网站免费观看| 激情在线观看视频在线高清| 精品久久久久久久末码| 亚洲欧美日韩高清在线视频| 亚洲精品美女久久av网站| 脱女人内裤的视频| 舔av片在线| 国产97色在线日韩免费| 一级a爱片免费观看的视频| 成人亚洲精品av一区二区| 黄色日韩在线| 亚洲18禁久久av| 国产真实乱freesex| 一级毛片女人18水好多| 露出奶头的视频| 亚洲专区字幕在线| 欧美在线一区亚洲| 日本黄色视频三级网站网址| 国产成人aa在线观看| 久久久久久久久中文| 听说在线观看完整版免费高清| 亚洲人成伊人成综合网2020| 日韩成人在线观看一区二区三区| 精品久久久久久成人av| 国产私拍福利视频在线观看| 日韩有码中文字幕| 村上凉子中文字幕在线| 国产精品 欧美亚洲| 国产精品免费一区二区三区在线| 国产精品野战在线观看| 亚洲国产精品999在线| 99久国产av精品| 老司机在亚洲福利影院| 青草久久国产| 18禁观看日本| 99久国产av精品| 亚洲专区国产一区二区| 欧美色视频一区免费| 1000部很黄的大片| 女人高潮潮喷娇喘18禁视频| 最近在线观看免费完整版| 亚洲色图 男人天堂 中文字幕| 亚洲最大成人中文| 亚洲一区二区三区色噜噜| 亚洲专区中文字幕在线| 亚洲av第一区精品v没综合| 观看免费一级毛片| 欧美高清成人免费视频www| 国语自产精品视频在线第100页| 国产激情欧美一区二区| xxx96com| 波多野结衣高清无吗| 桃红色精品国产亚洲av| 搞女人的毛片| 99国产精品一区二区三区| 国产精品一及| 久久婷婷人人爽人人干人人爱| 国产高清视频在线观看网站| 国产一区在线观看成人免费| 亚洲欧美激情综合另类| 日韩有码中文字幕| 白带黄色成豆腐渣| 欧美中文日本在线观看视频| 亚洲人成网站在线播放欧美日韩| 69av精品久久久久久| 天堂网av新在线| 欧美日韩国产亚洲二区| 国产熟女xx| 欧美日韩中文字幕国产精品一区二区三区| 午夜精品久久久久久毛片777| av在线蜜桃| 男女床上黄色一级片免费看| www.熟女人妻精品国产| 国产精品香港三级国产av潘金莲| 又爽又黄无遮挡网站| 日韩欧美在线二视频| av中文乱码字幕在线| 欧美日韩福利视频一区二区| 日韩欧美 国产精品| 久久精品国产亚洲av香蕉五月| 老司机午夜福利在线观看视频| 精品电影一区二区在线| 精品熟女少妇八av免费久了| 精品不卡国产一区二区三区| 国产精品亚洲一级av第二区| 两个人视频免费观看高清| 免费在线观看亚洲国产| ponron亚洲| 男人和女人高潮做爰伦理| 国产成人欧美在线观看| 国内少妇人妻偷人精品xxx网站 | 色尼玛亚洲综合影院| 久久精品aⅴ一区二区三区四区| 美女扒开内裤让男人捅视频| 亚洲成人精品中文字幕电影| 一个人免费在线观看的高清视频| 一级毛片高清免费大全| 久久久国产精品麻豆| 亚洲中文日韩欧美视频| 久久香蕉国产精品| 日韩大尺度精品在线看网址| h日本视频在线播放| 成人国产一区最新在线观看| 欧美色视频一区免费| 国产高潮美女av| 美女高潮喷水抽搐中文字幕| 欧美大码av| 宅男免费午夜| 九色国产91popny在线| 亚洲国产高清在线一区二区三| 一进一出抽搐gif免费好疼| 久久国产乱子伦精品免费另类| 日韩有码中文字幕| 中文亚洲av片在线观看爽| 精品乱码久久久久久99久播| 毛片女人毛片| 十八禁人妻一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲熟女毛片儿| 我要搜黄色片| 国产成人影院久久av| 国产亚洲精品av在线| 国产伦一二天堂av在线观看| 黄片小视频在线播放| 国产高清三级在线| 99热6这里只有精品| 久久久精品欧美日韩精品| 久久人妻av系列| 精品国产美女av久久久久小说| 夜夜躁狠狠躁天天躁| 欧美黑人欧美精品刺激| 熟女少妇亚洲综合色aaa.| 91字幕亚洲| 在线看三级毛片| 国产精品九九99| 欧美在线一区亚洲| 国产精品99久久久久久久久| 真人做人爱边吃奶动态| 国产高清视频在线观看网站| 男女之事视频高清在线观看| 熟女人妻精品中文字幕| 精品免费久久久久久久清纯| 亚洲一区二区三区不卡视频| 国产高清视频在线播放一区| 日日干狠狠操夜夜爽| 日韩国内少妇激情av| 18禁美女被吸乳视频| 在线十欧美十亚洲十日本专区| 精品国产乱码久久久久久男人| 熟妇人妻久久中文字幕3abv| 在线a可以看的网站| 老司机深夜福利视频在线观看| 女生性感内裤真人,穿戴方法视频| 国产aⅴ精品一区二区三区波| 国产一区在线观看成人免费| 成熟少妇高潮喷水视频| 色噜噜av男人的天堂激情| 久久亚洲精品不卡| 91麻豆av在线| 国产不卡一卡二| 久久久久精品国产欧美久久久| 久久久久久久久久黄片| 香蕉丝袜av| 激情在线观看视频在线高清| 99久久无色码亚洲精品果冻| 久99久视频精品免费| 久久欧美精品欧美久久欧美| 久久久久性生活片| 成人三级做爰电影| 久久精品夜夜夜夜夜久久蜜豆| 亚洲狠狠婷婷综合久久图片| 99国产精品99久久久久| 久久久水蜜桃国产精品网| 久久国产精品人妻蜜桃| 9191精品国产免费久久| 久久精品综合一区二区三区| 久久久久久大精品| 夜夜躁狠狠躁天天躁| 淫秽高清视频在线观看| 人人妻人人看人人澡| 国产高潮美女av| 琪琪午夜伦伦电影理论片6080| 亚洲一区二区三区不卡视频| 国产成人精品无人区| 女同久久另类99精品国产91| a级毛片在线看网站| 国产精品久久久久久久电影 | 国产亚洲欧美在线一区二区| 色综合婷婷激情| 搡老妇女老女人老熟妇| 国产成人av教育| 午夜福利高清视频| 国产精品免费一区二区三区在线| 亚洲中文日韩欧美视频| 免费看美女性在线毛片视频| 久久性视频一级片| 老熟妇仑乱视频hdxx| 国产97色在线日韩免费| 久久精品国产清高在天天线| 可以在线观看毛片的网站| 亚洲av第一区精品v没综合| 国产精品1区2区在线观看.| 给我免费播放毛片高清在线观看| 国产午夜精品论理片| 两个人的视频大全免费| 成年女人毛片免费观看观看9| 九九久久精品国产亚洲av麻豆 | 三级男女做爰猛烈吃奶摸视频| av天堂在线播放| 精品日产1卡2卡| 久久这里只有精品中国| 亚洲九九香蕉| 一二三四在线观看免费中文在| 国模一区二区三区四区视频 | 日本一二三区视频观看| 午夜久久久久精精品| 一本综合久久免费| 亚洲国产欧美一区二区综合| 国产一区二区三区视频了| 免费在线观看亚洲国产| 久久亚洲精品不卡| 国产综合懂色| 国产亚洲av高清不卡|