CrdR function in a curdlan-producing Agrobacterium sp. ATCC31749 strain

Background Agrobacterium sp. ATCC31749 is an efficient curdlan producer at low pH and under nitrogen starvation. The helix-turn-helix transcriptional regulatory protein (crdR) essential for curdlan production has been analyzed, but whether crdR directly acts to cause expression of the curdlan biosynthesis operon (crdASC) is uncertain. To elucidate the molecular function of crdR in curdlan biosynthesis, we constructed a crdR knockout mutant along with pBQcrdR and pBQNcrdR vectors with crdR expression driven by a T5 promoter and crdR native promoter, respectively. Also, we constructed a pAG with the green fluorescent protein (GFP) gene driven by a curdlan biosynthetic operon promoter (crdP) to measure the effects of crdR expression on curdlan biosynthesis. Results Compared with wild-type (WT) strain biomass production, the biomass of the crdR knockout mutant was not significantly different in either exponential or stationary phases of growth. Mutant cells were non-capsulated and planktonic and produced significantly less curdlan. WT cells were curdlan-capsulated and aggregated in the stationery phase. pBQcrdR transformed to the WT strain had a 38% greater curdlan yield and pBQcrdR and pBQNcrdR transformed to the crdR mutant strain recovered 18% and 105% curdlan titers of the WT ATCC31749 strain, respectively. Consistent with its function of promoting curdlan biosynthesis, curdlan biosynthetic operon promoter (crdP) controlled GFP expression caused the transgenic strain to have higher GFP relative fluorescence in the WT strain, and no color change was observed with low GFP relative fluorescence in the crdR mutant strain as evidenced by fluorescent microscopy and spectrometric assay. q-RT-PCR revealed that crdR expression in the stationary phase was greater than in the exponential phase, and crdR overexpression in the WT strain increased crdA, crdS, and crdC expression. We also confirmed that purified crdR protein can specifically bind to the crd operon promoter region, and we inferred that crdR directly acts to cause expression of the curdlan biosynthesis operon (crdASC). Conclusions CrdR is a positive transcriptional regulator of the crd operon for promoting curdlan biosynthesis in ATCC31749. The potential binding region of crdR is located within the −98 bp fragment upstream from the crdA start codon Electronic supplementary material The online version of this article (doi:10.1186/s12866-015-0356-1) contains supplementary material, which is available to authorized users.


Knockout of the CrdR (AGRO_0435 gene)
For AGRO_0435 gene knockout, a 2,533 bp fragment of the target gene (AGRO_0435) with up-and down-stream flanking sequences was PCR cloned using primers crdRup-1 and crdRdown-2. The amplified fragment, double digested with both SacI and HindIII, was inserted into the same sites of pUC19 to obtain the pUCcrdR flanking. The gentamicin (Gm) resistance gene expression cassette, obtained by PCR amplification with primers (Gm-1 and Gm-2) from pEX18Gm [37], was double digested with BamHI and SalI and then inserted into the same sites of the pUCcrdR flanking to obtain the suicide knockout vector pUC19T-crdR ( Figure 1A). Knockout plasmids were transformed with electroporation using an pairs of primers including crdR-1 and crdR-2, crdR-1 and crdRcheckGm-2, and crdRcheckGm-1 and crdR-2 ( Table 2).

Construction of pBQcrdR and pBQNcrdR for homogenous AGRO_0435 expression
A 437-bp full-length coding region of crdR was amplified with PCR using the primer pairs crdR-1 and crdR-2 ( Table 2) with genomic ATCC31749 DNA. The amplified crdR fragment, digested with BamHI and SacI, was ligated into the pBQ vector [38], creating the pBQcrdR vector (Table 1). To construct the vector for crdR expression driven by its native promoter of crdR, an AGRO_0435 fragment with up-and down-stream flanking sequences was PCR cloned using primers NcrdR-1 and NcrdR-2 (Table 2) The obtained 1,302 bp PCR fragment which was double digested with both SacI and XhoI was inserted into SalI and sacI sites of pBQ to create pBQNcrdR (Table 1, Figure 1).

Construction of pAG vector with GFP expression driven by the crd operon promoter
The predicted crd promoter (crdP), which is a 607-bp fragment upstream from the start codon of crdA (ATG), was amplified from genomic ATCC31749 DNA with primers crdAPG-1 and crdAPG-2 ( Table 2). The GFP code sequence was amplified with primers GFP-1 and GFP-2 ( Table 2) from plasmid pEGFP (Clontech, Mountain View, CA) and the two fragments were fused via PCR amplification. The resultant fused fragment, digested with XhoI and BamHI, was inserted into the same sites of plasmid pQE81L to yield pQEAG. After digestion with SalI, the fragment containing the gram-negative broad host replicating origin and the kanamycin (Kan) resistant gene amplified from pBBR122 (Table 1) with primers pairs Rep-Kan-1 and Rep-Kan-2 (Table 2), was inserted into the XhoI site of pQEAG to yield pAG (Table 1).

Curdlan fermentation and yield analysis
A two-step fermentation protocol was used to measure curdlan yields. In brief, ATCC 31749 and modified strains were inoculated into test tubes containing 5 mL LB and grown overnight at 30°C with 200 revolutions per minute (rpm). About 2 mL each of the seed cultures (SC) were transferred into 500-mL flasks containing 100 mL LB with or without IPTG (final concentration 0.5 mM) at 25°C, 200 rpm for 4 h. Cells were collected by centrifugation (1000 × g for 10 min, 4°C) and cell pellets were added to 125 mL curdlan-producing medium in a 500-mL flask which was shaken at 200 rpm. Every 24 h for 5 days reaction, 15-mL samples were taken from the culture mixture and samples were centrifuged at 8,000 × g for 5 min to collect pellets. Pellets containing both cells and curdlan were resuspended in 15 mL NaOH solution (1 mol/L) for 2 h. Cells pellets were separated by centrifugation at 8,000 × g for 5 min and resulting curdlan was precipitated by the addition of 2.0 mol/L HCl and the pH was adjusted to 6.5. Curdlan was recovered by centrifugation, washed, and dried to a constant weight in an oven (80°C).
crdR, crdA, crdS, and crdC expression analysis using q-RT-PCR  (Table 2), were ligated into pMD18-T vectors respectively. Then, using those constructs as standard copies, q-RT-PCR quantification was performed using an Applied Biosystems 7500 fast realtime PCR system (Applied Biosystems, Grand Island, NY ) with SYBR Premix E*TaqII (TaKaRa). All samples were run in triplicate [39].

GFP expression
The constructed vector (pAG) was transformed into both wild-type ATCC31749 and a crdR mutant (ATCC3174 9ΔcrdR). Transformed bacterial cells were grown in LB for 12 h and curdlan-producing medium for 72 h at 30°C. GFP expression was observed under an optical microscope (Zeiss Observer Berlin, Germany), equipped with epi-fluorescence. Simultaneously, with excitation of 450-490 nm light, Green fluorescence of GFP was measured by a fluorospectrophotometer F97Pro (FProd, Shanghai, China) to collect the data of the emission spectrum and relative fluorescence of cells harvested from both bacterial cell -growing and curdlan-producing phases.
Expression and purification of 6 His-tagged crdR protein 6-His-tagged crdR was expressed in E. coli BL21 through pBQcrdR transformation. The resultant strain grew at 37°C in LB medium (OD 600 nm = 0.5-0.6), and crdR protein expression was induced by adding IPTG (final concentration = 0.5 mM). The culture was shaken at 30°C for 4 h at 220 rpm. Cells were harvested by centrifugation were immediately extracted or frozen at −80°C until they were used. 6-His-tagged protein was purified by affinity chromatography using One-Step His-Tagged Protein Miniprep Pack (TIANDZ, Shanghai, China) according to the manufacturer's instruction. Purified crdR dissolved in elution buffer was dialyzed with dialysis buffer (100 mM KAc, 1 mM CaCl 2, 1 mM MgAc 2 , 1 mM EDTA, and 1 mM dithiothreitol, 10% glycerol) overnight in a semipermeable membrane. Protein concentration was measured using an improved Coomassie assay with bovine serum albumin (BSA) as standard.
DNA binding analysis of CrdR by EMSA DNA fragments containing various lengths of the crd promoter (crdP, crdP142, crdP108, crdP98, crdP53, crdP13, and crdP1) and~450 bp upstream of the start codon (ATG) of crdR (relA(AGRO_1479) celA (AGRO_4469), and crdS (AGRO_1848) named relAP celAP and crdSP were obtained by PCR amplification with primers listed in Table 2 respectively, those fragments were purified with a DNA gel extraction kit (Sangon Biotech) respectively according to the manufacturer's protocol. A electrophoretic mobility shift (EMSA) binding assay was performed as previously described with slight modifications [40]. Briefly, 10 μL of 0.25-0.50 mg/mL purified His-tagged crdR in 4× EMSA buffer (15 mM HEPES, 100 mM NaCl, 1 mM EDTA, 5 mM MgCl 2 , 10% glycerol, 1 μg/mL poly dI-dC) was incubated with 10 μL of different purified target DNA fragments (0.5 μM) in ddH 2 O at room temperature for 30 min. DNA-protein complexes were loaded onto a 2% agarose gel and separated at 80 V for 1.5 h, and the gel was stained with SYBR Green I and visualized with a UV trans-illuminator (Upland, CA).

crdR knockout mutant construction and phenotypes
The crdR knockout mutant was constructed via homologous recombination by transformation of the suicide plasmid, pUC19T-crdR ( Figure 1A). After strains were selected on gentamicin (Gm) resistant LB plates, knockout mutants were confirmed with PCR amplification (Figure 2). Compared with wild-type (WT) ATCC31749, which is capsulated in the stationary phase, the crdR knockout strain (ATCC31749ΔcrdR) produced less curdlan (Figures 3  and 4) leading to motile and non-capsulated planktonic forms ( Figure 3) in both exponential and stationary phases. crdR expression driven by promoters of T5 and native crdR in both ATCC31749 and ATCC31749ΔcrdR strains, respectively, were obtained by transforming the constructs of pBQcrdR and pBQNcrdR ( Figure 1). The ATCC31749ΔcrdR/pBQNcrdR strain recovered its curdlan capsulated form of ATCC31749 ( Figure 3D and 3E). With a two-step flask-shaking process, curdlan production over 5 days [18,19] in 5 cultivated strains-ATCC31749, ATCC31749/pBQcrdR, ATCC31749ΔcrdR, ATCC31749Δ crdR/pBQcrdR and ATCC31749ΔcrdR/pBQNcrdR-was compared to assess biomass accumulation and curdlan yields. Data show that ( Figure 4B) the biomass of ATCC31749ΔcrdR, ATCC31749/pBQcrdR and ATCC3 1749ΔcrdR/pBQNcrdR were higher and that ATCC31 749ΔcrdR/pBQcrdR had less accumulation than WT ATCC31749 during cultivation days 2 to 4. By the 5 th day, however, cell biomasses of all strains were similar. Curdlan yields for ATCC31749, ATCC31749/pBQcrdR, ATCC31749ΔcrdR, ATCC31749ΔcrdR/pBQ-crdR and ATCC31749ΔcrdR/pBQNcrdR were 5.66 g/L, 7.80 g/L, 0.007 g/L, 1.13 g/L, and 5.91 g/L, respectively ( Figure 4A). Curdlan yield for the crdR overexpressed strain (ATC C31749/pBQcrdR) was 38% greater than that of WT. crdR controlled by the T5 promoter in ATCC31749ΔcrdR/ pBQcrdR could synthesize curdlan, but yields only reached 18% of WT yields; however crdR controlled by its native promoter in ATCC31749ΔcrdR/pBQcrdR recovered curdlan yields of crdR knockout. Low curdlan yields caused by expression of crdR controlled with T5 promoter suggests a more complex regulatory mechanism of crdR expression in the strains. Judged from these data, our observation of crdR is consistent with reports to suggest that crdR is an important regulator of curdlan biosynthesis [15,26].

Expression of curdlan biosynthesis genes responding to crdR overexpression
Because crdR is an important regulator of curdlan biosynthesis (Figures 3 and 4), we investigated whether crdR activates expression of crd operon genes. q-RT-PCR analysis was used to evaluate the effects of crdR on crdA, crdS, and crdC (genes of the curdlan biosynthetic operon) mRNA. Stationary phase cells favoring curdlan biosynthesis were compared with exponential phase cells favoring cell growth, and crdR native expression was found to be 29.2 copies/ng total RNA in the stationery phase and 14.0 copies/ng total RNA in the exponential phase in ATCC31749/pBQ. Correspondingly, expression of crdA, crdS, and crdC was at least 10 times greater in the stationary phase compared to the exponential phase ( Table 3). mRNA of crdR in ATCC31749/pBQcrdR, induced by 0.5 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) at 30°C for 2 h, was significantly increased compared to that in ATCC31749/pBQ. Corresponding mRNA of crdA, crdS and crdC in ATCC31749/pBQcrdR were more than twice greater when strains were cultivated in both growth and fermentation media. Data confirmed that crdR promotes curdlan production via activating expression of crd operon genes (Table 3).

GFP expression controlled by the crd operon promoter (crdP)
To confirm the effect of crdR on crd operon gene expression, a shuttle vector of pAG bearing GFP driven by the crdP was constructed. pAG was transformed into both ATCC31749 and its crdR mutant strain of ATCC317 49ΔcrdR. Data indicate that green fluorescence was undetectable by fluorescent microscopy in the crdR mutant strain. In contrast, strong green fluorescence was visible in  the WT ATCC31749 strain grown in fermentation media ( Figure 5A and 5B). Also, the ATCC31749/pAG strain was curdlan capsulated whereas ATCC31749ΔcrdR/ pAG was non-capsulated and motile ( Figure 5C and 5D). GFP expression detected by spectrophotometry was consistent with microscopic observations that the relative GFP florescence in ATCC31749/pAG exceeded that in ATCC31749ΔcrdR/pAG at both exponential and stationary phases ( Figure 5D). Because the GFP expression pattern represents expression profiles of crdA, crdS and crdC in engineered strains of both ATCC31749/pAG and ATCC31749ΔcrdR/pAG, the effects of crdR expression on GFP expression in those strains indicated that crdR might directly or indirectly interact with the crd operon promoter to regulate expression(s) of curdlan synthetic gene(s). Interestingly relative GFP florescence of exponential phase ATCC31749ΔcrdR/pAG was less than that measured in the stationary phase, suggesting that crdR may synergistically cooperate with other regulators to control crd operon gene expression.

crdR binding with' different crd operon promoter regions
Bioinformatic analysis of deduced amino acid sequences indicates that crdR has a conserved DNA-binding motif of a helix-turn-helix domain. To confirm that crdR protein can directly interact with the crd operon promoter region, with BSA as a negative protein control, DNAbinding analysis was performed with an electrophoretic mobility shift assay (EMSA) with amino terminal 6-Histagged crdR and DNA fragments containing its putative binding sites. 6 His-tagged crdR protein was expressed and purified into a single band (15 kDa) with SDS-PAGE from E. coli Bl21 which was transformed with pBQcrdR ( Figure 6B). The various crdR putative binding fragments, including serial regions of the crd promoter ranging from 607 bp upstream from ATG of crdA fused with or without of 395 bp downstream of ATG of crdA ( Figure 6A) and about 450 bp upstream from ATG of crdR, (p) ppGpp synthetase (relA), cellulose synthase catalytic subunit (celA) and curdlan synthase catalytic subunit (crdS), were amplified by PCR with the genomic DNA of ATCC31749 as a template ( Figure 6E). Data indicate that 6 His-tagged crdR protein cannot specifically bind 450 bp upstream of ATG of crdR, relA, celA, and crdS, but can bind to different promoter regions of the crd operon ( Figure 6C and 6D).
To locate crdR binding site(s) on the crd operon promoter, 607 bp upstream from ATG of the crdA sequence which can successfully drive GFP expression in the stationary phase ( Figure 5), was chosen and the fragment was shorted to 98 bp upstream from the ATG of crdA to measure binding abilities to crdR. Data indicate that those sequences did not reduce binding to crdR as evidenced by a band shift in EMSA by mixing with or without crdR ( Figure 6C). Then, the −98 bp fragments upstream from the crdA start codon was focused for further analysis by continual shortening to 13 bp upstream of ATG of crdA. Using the 395 bp coding sequence of crdA from 1 to 395 bp and BSA as a negative DNA and protein control, respectively, fragment mobility shifts containing −98, −53 *The data are means of three independent determinations, and the unite of value is copies per ng total RNA(copies/ng RNA) E: strain in exponential phase, S:strain in stationary phase.

Discussion
Here, we report that the crdR, a homolog of PhrR of S. meliloti can activate curdlan synthetic gene expression in Agrobaterium sp. ATCC31749. To our knowledge, ours is the first report to depict molecular functions of the crdR gene. Our data indicate that curdlan yield in an over-expressing crdR strain increased 38% compared to the WT strain. Also, pBQNcrdR transformed to the crdR mutant strain recovered 105% curdlan synthesis of the WT strain ( Figure 4). Also, when pAG was transformed into both crdR mutant and WT strains GFP expression controlled by the crd promoter was undetectable by fluorescent microscopy with low relative fluorescence in the crdR mutant. In contrast, the WT strain had visible green color with high relative fluorescence. Finally, q-RT-PCR analysis indicated that crdR is highly expressed in the stationary phase and that overexpression of crdR in the WT strain significantly increased expression of crdA, crdS and crdC. These data agree with previous reports that crdR is key for regulating curdlan biosynthesis. Purified crdR from E. coli BL21 can also specifically bind to the promoter region of crd offering initial evidence that crdR is a positive transcriptional regulator of the crd operon in ATCC31749. The biomass accumulation in crdR mutant strains was not significantly different from the WT strain, suggesting that the crdR gene is not required for cell growth. Microscopic observation revealed that the crdR mutant was nearly curdlan deficient, resulting in mutant cells with non-capsulated planktonic forms. The WT ATCC31749 strain and the complementary strain of the crdR mutant, ATCC31749ΔcrdR/pBQNcrdR, accumulated curdlan in the stationary phase in culture media with low pH and limited nitrogen, leading to cells were capsulated and aggregated ( Figure 3). In addition, expression of crdR was higher in the stationary phase than in the exponential phase, and crdR expression further activated curdlan biosynthesis in the ATCC31749 strain to generate a biofilm. This suggests that curdlan may be critical for biofilm formation in ATCC31749 for improving stress tolerance to harsh conditions. Bioinformatic analysis indicated that crdR can be grouped into a conserved XRE-family of transcriptional factors that is comprised of HipB in E. coli, PhrR in S. meliloti, CH00371 in R. etli and PraR in A. caulinodans (Additional Figure 6 Binding ability of 6 His-tagged crdR to different DNA fragments. A: Different region of the crdP; B: purification of 6 His-tagged crdR; C, D: Binding ability 6 His-tagged crdR to different regions of crdP; E: Binding of 6 His-tagged crdR, to~450 bp upstream of ATG at different gene coding regions. B M : protein markers, B 1 : purified His-tagged crdR, B 2 : supernatant of pBQcrdR/E. coli Bl21; C 1 . C 4 , C 7 , and C 10 are 10 μL of 0.5 μM crdP 98, crdP102, crdP142 and crdP mixed with 10 L of 0.5 mg/mL purified 6 His-tagged crdR protein respectively, C 2 , C 5 , C 8 , and C 11 are same as C 1 . C 4 , C 7 , and C 10 , except 6 His-tagged crdR protein was reduced to 0.25 mg/mL; C 3 , C 8 , C 9 , and C 12 are same as C 1 . C 4 , C 7 , and C 10 , without 6 His-tagged crdR. D 1 . D 4 , D 7 , and D 10 are 10 μL of 0.5 μM crdP 98, crdP53, crdP13, and crdP1 only respectively, D 2 , D 5 , D 8 , and D 11 are same of .D 1 . D 4, D 7 , and D 10 mixed with 10 μL 0.5 mg/mL 6 His-tagged crdR protein; D 3 , D 8 , D 9 , and D 12 are same of .D 1 , D 4 , D 7 , and D 10 mixed with 10 μL 0.5 mg/mL BSA respectively, E 1 , E 4 , E 7 , and E 10 are crdRP, celAP, and crdSP, E 2 , E 5 , E 8 , and E 11 are same of E 1 , E 4 , E 7 , and E 10 mixed with 10 μL of 0.25 mg/mL 6 His-tagged crdR protein respectively, E 3 , E 6 , E 9 , and E 12 are same of E 1 , E 4 , E 7 , and E 10 mixed with 10 μL of 0.5 mg/mL BSA respectively. file 1) [27]. Apparently, diverse stress can induce expressions of XRE-family transcriptional factors. Combining with HipA, HipB, a crdR homologue of E. coli that mediates multidrug stress tolerance can bind to its cis elements with conserved sequences of TATCCN 8 GGATA (where N 8 indicates any 8 nucleotides). Genomic scanning indicated that there is no HipA counterpart in the ATCC31749 strain, and that there were no conservative TATCCN8G-GATA sequences in the promoter region of the crd operon. However, crdP does have three distinct hairpin structures located at the −10, −35, and −92 regions (Additional file 2), which are putative crdR binding sites. That purified crdR can bind to an amplified fragment containing the −92 region of crdP more than the −53, and −13 regions indicates that those region are likely the crdR binding site. HipAB is a heterodimer of a transcriptional repressive regulator [28], and crdR may play a different role as a transcriptional activator in the form of a homotetramer or homodimer, which must be confirmed with additional studies. PhrR in S. meliloti was affected by low pH, Cu 2 + , Zn 2 + and H 2 O 2 stresses [26]. Expression of CH00371 in R. etli was promoted by oxidation and osmotic shock [29,30]; whereas expression of PraR in A. caulinodans was increased by low nitrogen [31]. Primary data from transcriptome analysis obtained from RNAseq indicates that expression of crdR in the curdlan fermenting phase with both low pH and limited nitrogen was twice as great as that in the growing phase of bacterium in LB medium (data not shown) and q-RT-PCR analysis of crdR expression agreed with RNA-Seq data (Table 3). Thus, crdR expression should be triggered by stress factors as well.
Most organisms within the genera of Rhizobium, Azorhizobium, Bradyrhizobium, Mesorhizobium, and Sinorhizobium are symbiotic bacteria to various leguminous plants [32,33]. After bacteria enter plant tissues, their environment changes to low pH (normally 5.5), with limited nitrogen and sufficient carbohydrates [32,33], and this environmental shift affects strain morphology and physiology. To survive with limited nitrogen, genes related to bacterial nitrogen fixing and those for low pH tolerance are expressed. Transcriptional factors such as PhrR and its homologues are instrumental for expression of low pH tolerance genes. Also, with abundant carbohydrates from host photosynthesis, bacterial extracellular polysaccharides (EPS) are synthesized to produce a biofilm [34] to protect the bacteria from stress. Also, the biofilm provides a low-oxygen environment inside the bacteria to support nitrogen fixing reactions. This concept is in agreement with research with PraR from A. caulinodans ORS571, a homolog of crdR, which can mediate stem nodule formation, regulate expression of Reb genes, and increase nitrogen fixation of bacterial strains within stem nodules [31]. Therefore, biosynthesis of curdlan regulated by crdR may originate from ancestral characteristics for survival within a host plant. Currently, detailed regulatory mechanisms of crdR expression, controlled by disadvantageous conditions such as low pH and limited nitrogen remain unknown. However, reports regarding stringent response signal(p)ppGpp, which can induce EPS biosynthesis and bacterial biofilm formation [35] is worthy of study. We hypothesize that crdR expression may be activated by (p)ppGpp, which accumulates under stressful conditions.

Conclusions
In this study, we confirmed that crdR regulates curdlan synthesis by activating expressions of its biosynthetic genes. Ours is the first work to identify XER family transcriptional factor which can activate EPS biosynthesis. crdR may be a multiple-effect regulator controlling expression(s) of the curdlan synthesis gene(s) in ATCC31749 under oxidative stress/low pH and/or limited nitrogen with abundant sugar. This function of curdlan regulation indicates that curdlan biosynthesis of ATCC31749 under harsh conditions may have evolutionary origin.
Additional file 2: A98 pb promoter sequence of the crd operon (crdP) upstream from starting codon of crdA (ATG) with its secondary structure predicted by DNAMAN 98 bp fragment upstream ATG) of the crdA have 3 putative crdR binding sites, which has 3 hairpin structures located in −92, −35, and −10 regions.
Authors' contributions XY constructed vectors pAG and pBQ-crdR, analyzed curdlan yield and measured GFP expression, performed q-RT-PCR analysis, crdR purification, and EMSA and drafted the manuscript. CZ constructed vectors of pUC19T-crdR and pBQNcrdR, assisted with crdR knockout and edited the paper. LY assisted with crdR knockout. LZ assisted with statistics and bioinformatic analysis. CL assisted with curdlan yield analysis. ZL assisted with experimental conditions and manuscript editing. ZM conceptualized the study and coordinated the study design as well as assisted with the paper. All authors read and approved the final manuscript.