Skip to main content

METHODS article

Front. Plant Sci., 21 January 2022
Sec. Technical Advances in Plant Science
This article is part of the Research Topic Plant Transformation View all 20 articles

Optimization of Tissue Culturing and Genetic Transformation Protocol for Casuarina equisetifolia

\r\nHuimin Ren&#x;Huimin Ren1†Yan Xu&#x;Yan Xu1†Xiaohong ZhaoXiaohong Zhao1Yan ZhangYan Zhang1Jamshaid HussainJamshaid Hussain2Fuqiang CuiFuqiang Cui1Guoning Qi*Guoning Qi1*Shenkui Liu*\r\nShenkui Liu1*
  • 1State Key Laboratory of Subtropical Silviculture, School of Forestry and Biotechnology, Zhejiang A & F University, Hangzhou, China
  • 2Department of Biotechnology, COMSATS University Islamabad, Abbottabad, Pakistan

Casuarina equisetifolia is widely used in agroforestry plantations for soil stabilization, ecosystem rehabilitation, reclamation, and coastal protection. Moreover, C. equisetifolia has remarkable resistance to typhoons, desert, low soil fertility, drought, and salinity, but not cold. Therefore, it is significant to breed high-quality Casuarina varieties to improve the tolerance and adaptability to cold weather by molecular techniques. The establishment of a rapid and efficient callus induction and regeneration system via tissue culture is pre-requisite for the genetic transformation of C. equisetifolia, which is so far lacking. In this study, we reported an efficient and rapid regeneration system using stem segment explants, in which callus induction was found to be optimal in a basal medium supplemented with 0.1 mg⋅L–1 TDZ and 0.1 mg⋅L–1 NAA, and proliferation in a basal medium containing 0.1 mg⋅L–1 TDZ and 0.5 mg⋅L–1 6-BA. For bud regeneration and rooting, the preferred plant growth regulator (PGR) in basal medium was 0.5 mg⋅L–1 6-BA, and a combination of 0.02 mg⋅L–1 IBA and 0.4 mg⋅L–1 IAA, respectively. We also optimized genetic a transformation protocol using Agrobacterium tumefaciens harboring the binary vector pCAMBIA1301 with β-glucuronidase (GUS) as a reporter gene. Consequently, 5 mg L–1 hygromycin, 20 mg L–1 acetosyringone (As), and 2 days of co-cultivation duration were optimized to improve the transformation efficiency. With these optimized parameters, transgenic plants were obtained in about 4 months. Besides that, Agrobacterium rhizogenes-mediated transformation involving adventitious root induction was also optimized. Our findings will not only increase the transformation efficiency but also shorten the time for developing transgenic C. equisetifolia plants. Taken together, this pioneer study on tissue culturing and genetic transformation of C. equisetifolia will pave the way for further genetic manipulation and functional genomics of C. equisetifolia.

Introduction

Casuarina equisetifolia, a symbiotic nitrogen-fixing tree, plays an important role in economic and environmental improvements in tropical and subtropical littoral zones of Asia, the Pacific, and Africa (Diem and Dommergues, 1990; Zhang et al., 2008; Zhong et al., 2010). It is widely used in agroforestry plantations for several purposes, such as soil stabilization, ecosystem rehabilitation, reclamation, and coastal protection (Pinyopusarerk and Williams, 2000). C. equisetifolia is also appreciated as a source of fuel wood, paper, timber, medicine, dye, pulp, tannin, and charcoal (Pan et al., 1996; Zhong et al., 2005, 2010; Diagne et al., 2013b; Potgieter et al., 2014).

Casuarina equisetifolia is one of the typical angiosperms with distinctive needle-like branchlets, due to which water loss by evapotranspiration is significantly low in C. equisetifolia. Therefore, C. equisetifolia is a valuable tree to adapt to survival in drought stressed areas (Le et al., 1996). Moreover, it is reported that C. equisetifolia establishes symbiotic relationships with soil microorganisms like Frankia for nitrogen fixation in Casuarina root nodules (Benson and Silvester, 1993; Diagne et al., 2013a) and mycorrhizas for increasing phosphorus uptake (Santi et al., 2003). These symbiotic relationships enhance its adaptability to harsh environments such as salinity, barrenness, and drought.

Since the 1950s, C. equisetifolia has been extensively planted in infertile and sandy arid areas along the South China coastline as a windbreak/shelterbelt. Nowadays, it is widely planted, and is one of the most important coastal forest species because of its remarkable resistance to typhoons, desert, low soil fertility, drought, and salinity (Tani and Sasakawa, 2003; Ye et al., 2019). However, low temperature is one of the main limitations for fast growth, cultivation expansion, and high yield of Casuarina species (Li et al., 2016). Therefore, it is significant to breed high-quality Casuarina varieties to improve its tolerance and adaptability to cold weather using molecular techniques (He et al., 2011). Although few studies have been conducted to understand how Casuarina species respond to low temperature at physiological level (He et al., 2011; Wu et al., 2012), little is known about the functional genes and molecular mechanism of cold response in Casuarina. Recently, whole genome sequencing of C. equis. ssp. Incana has been accomplished. Genome size is reported to be about 300 Mb (Zhong et al., 2005, 2010; Ye et al., 2019). This progress may greatly contribute in studies on molecular mechanisms underlying the phenotype and tolerance to diverse stresses in C. equisetifolia. Moreover, it will also pave the way for breeding new varieties by gene function analysis and genetic transformation.

A rapid and efficient callus induction and regeneration system via tissue culture is a pre-requisite and essential part for genetic transformation of C. equisetifolia, which has not been established so far. Previously, few studies have reported a regeneration and transformation system for Casuarina glauca Sieb. ex Spreng. and Casuarina cunninghamiana Miq. (Le et al., 1996; Franche et al., 1997; Santi et al., 2003; Obertello et al., 2005). However, owing to great morphological differences among Casuarina species (Ye et al., 2019), development of a reliable regeneration and transformation system is crucial to improve and produce new breeds. For plant genetic engineering, Agrobacterium-mediated plant transformation has been the most successful and widespread technique in recent decades (Gelvin, 2003b; Vain, 2007). Besides A. tumefaciens-mediated transformation, in recent years, A. rhizogenes-mediated transformation has been performed to induce adventitious roots, named “hairy roots” at the site of wounding and infection in woody plants (Ozyigit et al., 2013; Hwang et al., 2017).

In this study, we report an efficient, rapid regeneration and two transformation systems for C. equisetifolia involving A. tumefaciens and A. rhizogenes.

Methods

Plant Materials and Growth Conditions

Young shoots were cut from 30- to 90-day-old seedlings of C. equis. ssp. Incana grown in a greenhouse at 25°C in a 16-h/8-h day/night photoperiod. The shoots were rinsed in running water for 2–3 h and disinfected using 75% ethanol (v/v) for 30 s, followed by rinsing with sterile water five times. Then, the shoots were further disinfected with 6–8% (v/v) sodium hypochlorite for 8 min and rinsed with sterile water five times. Finally, the sterile young shoots were cut into small pieces of about 1 cm in length and placed in a basal medium (1/2 MS, 30 g L–1 sucrose and 3.5 g L–1 phytagel) supplemented with different phytohormone combinations. The pH of media was adjusted to 5.7. For each hormone treatment, three replicates were performed, each with 14 explants.

For transformation using A. rhizogenes, seeds of C. equisetifolia were sown in pots filled with soil and grown in the greenhouse at 25°C in a 16-h/8-h day/night photoperiod. After 30 days, healthy seedlings of about 15 cm height were selected for agro-infiltration.

Optimization of Regeneration

For callus induction, five combinations of plant growth regulators (PGRs) were used as follows: C1: 0.5 mg⋅L–1 6-benzylaminopurine (6-BA) and 0.05 mg L–1 kinetin (KT); C2: 1.5 mg⋅L–1 6-BA and 0.03 mg⋅L–1 naphthylacetic acid (NAA); C3: 1.5 mg⋅L–1 6-BA and 0.02 IBA; C4: 0.1 mg⋅L–1 NAA and 0.1 mg⋅L–1 thidiazuron (TDZ); and C5: 0.5 mg⋅L–1 6-BA and 0.1 mg⋅L–1 TDZ. Pieces of sterile young shoots were placed in the medium at 28°C in the dark. After about 15 days, the effect of each hormone treatment on callus induction was observed.

For bud regeneration, good-quality calluses were shifted to the media with the following growth regulator combinations: B1: 0.5 mg⋅L–1 6-BA; B2: 0.5 mg⋅L–1 6-BA and 0.05 mg⋅L–1 KT; B3: 0.75 mg⋅L–1 6-BA and 0.03 NAA; B4: 1.5 mg⋅L–1 6-BA and 0.02 mg⋅L–1 IBA; and B5: 0.5 mg⋅L–1 6-BA and 0.1 mg⋅L–1 TDZ. To establish the preferred conditions, the rate of bud regeneration was recorded after growth for 30 days at 28°C under a 16-h light (60–80 μmol photons m–2 s–1)/8-h dark photoperiod. The buds (about 1 cm in length) were shifted to a rooting medium for root induction. The rate of root induction was calculated after 20-day cultivation in three different rooting media as follows: R1: 0.02 mg⋅L–1 IBA and 0.04 mg⋅L–1 IAA; R2: 0.4 mg⋅L–1 IAA; and R3: 1.5 mg⋅L–1 IAA.

Optimization of Hygromycin Concentration for Transformant Selection

Sterile young stems were placed in an optimized callus induction medium with 250 mg⋅L–1 carbenicillin (Car) and six different concentrations of hygromycin (Hyg: 0, 0.5, 1, 3, 5, and 10 mg⋅L–1). Each treatment consisted of 20 calluses, and 3 replicates were performed. The status of the calluses was observed, and the rate of callus induction was determined after 15 days of incubation in the dark.

Vector and Preparation of Agrobacterium tumefaciens for Transformation

Agrobacterium tumefaciens strain GV3101 was used for genetic transformation of C. equisetifolia. The binary expression vector pCAMBIA1301, with a β-glucuronidase (GUS) reporter gene under cauliflower mosaic virus (CaMV) 35S promoter, was mobilized into strain GV3101 using the freeze-thaw method (Holsters et al., 1978). A single colony of Agrobacteria carrying the pCAMBIA1301 binary vector was inoculated into an LB liquid containing 25 mg L–1 rifampicin and 50 mg L–1 kanamycin, and was grown overnight at 28°C with continuous shaking at 220 rpm. A 1-ml culture was transferred to 20-ml fresh liquid LB with appropriate antibiotics and 20 mg L–1 acetosyringone (As). The culture was incubated at 28°C until an OD600 of 0.6–0.8 was attained, and then centrifugated immediately at 5,000 rpm for 10 min. The pellet was resuspended in an approximately 50-mL ½ MS liquid containing 30 g L–1 sucrose and 20 mg L–1 As. The final concentration of the A. tumefaciens suspension was adjusted to an OD600 of 0.6 for transformation.

Evaluation of Factors Affecting Transformation Efficiency

For transformation, three co-cultivation durations (1, 2, 3 days) and four As concentrations (0, 5, 10, and 20 mg L–1) in an Agrobacterium infective suspension were evaluated. Stem segments of approximately 1 cm in length each were excised from the sterile young shoots and grown in an optimized callus induction medium. After 15 days, induced calluses were submerged into the Agrobacteria infective suspension and shaken at 28°C and 150 rpm for 2 h. Excess bacterial suspensions were removed from the calluses by drying on sterile filter papers. Subsequently, the calluses were placed in a basal medium with 0.1 mg L–1 NAA, 0.1 mg L–1 TDZ, and different concentration of As (0, 5, 10, and 20 mg L–1) for co-cultivation at 28°C. After co-cultivation, the calluses were washed 6–8 times with 250 mg L–1 carbenicillin water (each wash was 3 min) and submerged in 250 mg L–1 carbenicillin water for 8 min to decontaminate Agrobacteria. The induced calluses were then transferred into an optimized selection medium containing PGRs, hygromycin with different concentrations (0, 5, and 10 mg L–1), and carbenicillin (250 mg L–1). Each treatment consisted of 15 calluses in 3 replicates. All hygromycin-resistant shoots were further assayed by GUS histochemical staining.

β-Glucuronidase Staining Assay

Various tissues of transgenic young trees and control in different stages of transformation were stained for GUS staining analysis using the method described by Jefferson (Jefferson et al., 1987). After GUS staining, 70% (v/v) ethanol was used to remove the chlorophyll.

Efficient Hairy Root Transgenic System

A single colony of Agrobacterium rhizogenes strain K599 harboring vector pCAMBIA1300-35S:GFP was cultured in 4 ml of LB liquid medium with 50 mg L–1 streptomycin and 50 mg L–1 kanamycin, and with continuous shaking (180 rpm) at 28°C for about 8–12 h. The inoculum was then transferred to a 20-ml LB broth with the same antibiotics as 1:20. After the OD600 reached around 0.8, the cultures were centrifuged at 5,000 × g for 10 min at room temperature and then resuspended in 25 ml MES buffer containing 10 mM MES-KOH (pH 5.2), 10 mM MgCl2, and 100 μM As.

For transformation, seedlings of 5–8 cm height were gently picked from the pots and the residual soil on the roots was carefully removed. The stem between the cotyledons and roots was gently scratched using a clean blade, and then the entire plant was soaked into the infection solution for 2–4 h. After infection, the seedlings were replanted in the pots for growth. Scratched parts were exposed necessarily above the soil, so callus induction could easily be observed. Callus growth was observed about 15 days after infection in the wounded part of the stem. After about a month, when the hairy roots had developed well, the transgenic rate for the hairy roots was detected using LUYOR-3415RG (Hand-held Lamp). Transgenic efficiency was calculated as:

Transgenic root induction efficiency = (number of plants with induced hairy roots/total number of infected plants) × 100%;

Transformation efficiency = (positive transgenic hairy root lines/total number of induced hairy roots) × 100%.

Results

Improvement in Callus Induction and Regeneration of Casuarina equisetifolia

To establish a simple, rapid, and efficient regeneration system for C. equisetifolia, we tested 38 different combinations of plant growth regulators (PGR, cytokinin 6-BA, KT, TDZ, auxin IAA, IBA, an NAA) in a basal medium (1/2 MS medium, 30 g L–1 sucrose and 3.5 g L–1 phytagel) (Supplementary Table 1). Stem segments of approximately 1 cm in length each were excised from the sterile young shoots and grown in the abovementioned medium to induce callus (Supplementary Figure 1). After 15 days of growth in the dark in these media, five combinations exhibiting higher callus induction rate (over 80%) were selected for further experiments. The selected media were renamed as C1 (NO. 30: 0.5 mg⋅L–1 6-BA + 0.05 mg⋅L–1 KT), C2 (NO. 6: 1.5 mg⋅L–1 6-BA + 0.03 mg⋅L–1 NAA), C3 (NO. 26: 1.5 mg⋅L–1 6-BA + 0.02 mg⋅L–1 IBA), C4 (NO. 37: 0.1 mg⋅L–1 TDZ + 0.1 mg⋅L–1 NAA), and C5 (NO. 32: 0.1 mg⋅L–1 TDZ + 0.5 mg⋅L–1 6-BA). It was later observed that the calluses turned brown in the C2, C3, and C4 media (Figure 1A). The statistical results showed that the rate of callus induction was highest (around 98%) in the C4 medium, but that the frequency of callus browning was also higher in the C4 medium compared to that in the others (Figure 1B). On the contrary, the calluses in C5 had lower browning rate (Figure 1B). Thus, the C4 medium was used to induce calluses, while C5 was used for callus proliferation (Figures 1C,D).

FIGURE 1
www.frontiersin.org

Figure 1. Callus induction of Casuarina equisetifolia using the stem as explant. (A) Induced callus became brown. (B) Callus induction rate and browning callus rate in basal medium with different plant growth regulator (PGR) combinations. The results are presented as the mean and standard error from three independent experiments. (C) Inducement and (D) proliferation of calluses.

For budding induction, we found that budding can be successfully induced in medium numbers 7, 30, 5, 26, and 32, which were renamed as B1 (0.5 mg⋅L–1 6-BA), B2 (0.5 mg⋅L–1 6-BA + 0.05 mg⋅L–1 KT), B3 (0.75 mg⋅L–1 6-BA + 0.03 mg⋅L–1 NAA), B4 (1.5 mg⋅L–1 6-BA + 0.02 mg⋅L–1 IBA), and B5 (0.5 mg⋅L–1 6-BA + 0.1 mg⋅L–1 TDZ). Besides B1, which only contained 6-BA, the other media had 6-BA with KT, NAA, IBA, and TDZ. It took about 1–2 months for the calluses to develop adventitious buds (Figures 2A–E). With respect to the rate of budding induction, the B1 medium was found to be the best one (Figure 2F), indicating that only 6-BA was sufficient for budding induction while other hormones did not seem to play any obvious role. To sum it up, the basal medium supplemented with 0.5 mg⋅L–1 6-BA was found to be suitable for budding induction.

FIGURE 2
www.frontiersin.org

Figure 2. Adventitious bud induction of callus. (A–D) Different stages of induced adventitious bud growth. (E) Adventitious buds on a callus. (F) Bud induction rate in basal medium with different plant growth regulator (PGR) combinations. The results are presented as the mean and standard error from three independent experiments.

Next, we processed the buds for root regeneration. When the length of regenerated buds was more than 1 cm, these were divided and transferred to the rooting media. Among 38 tested combinations of growth hormones, 3 media (Nos. 20, 21, and 22) significantly promoted multiple root organogenesis. The statistical analysis of root induction rate showed that the rate in R2 (NO. 21: 0.02 mg⋅L–1 IBA + 0.4 mg⋅L–1 IAA) was 80%, while in R1 (NO. 20: 0.02 mg⋅L–1 IBA + 0.04 mg⋅L–1 IAA), and R3 (NO. 22: 0.02 mg⋅L–1 IBA + 1.5 mg⋅L–1 IAA) it was around 55 and 10%, respectively (Figure 3A). The number of induced roots in R2 was much higher than that in the R1 medium (Figure 3B). Roots emerged within a week after transfer in the rooting medium (Figure 3C). Finally, the whole seedlings, developed through tissue culture, were transferred into a basal medium supplemented with 0.5 g⋅L–1 activated carbon to accelerate and improve the quality of rooting (Zhao, 1999; Xu et al., 2006). Taken together, the regeneration of C. equisetifolia from explant to fully developed plantlet was accomplished in about 4 months, a greatly shortened duration for C. equisetifolia growth (Figure 4).

FIGURE 3
www.frontiersin.org

Figure 3. Root regeneration of adventitious buds. (A) Root induction rate in basal medium with different plant growth regulator (PGR) combinations. The results are presented as the mean and standard error from three independent experiments. (B) Comparing roots in the R1 and R2 media. (C) Different stages of induced root.

FIGURE 4
www.frontiersin.org

Figure 4. Efficient regenerating system for Casuarina equisetifolia using the stem as explant. (A) Callus induction in the C4 medium. (B) Callus proliferation in the C5 medium. (C) Bud induction in the B14 medium. (D) Root induction in the R1 medium. (E) Root regeneration of adventitious bud. (F) Seedlings using suitable regenerating system. (G) Transplantation of regenerated seedling.

In addition to the above method, stem cuttings from aseptic seedling were also used for rapid regeneration of plants. The 8–10 cm long branchlets were cut and transferred to the rooting medium (also supplemented with activated carbon). It took about 7–10 days to generate calluses and 15 days to develop roots (Figure 5A). Then, these were transferred to the basal medium supplemented with activated carbon (Figures 5B,C). The whole process of regeneration from stem cuttings took only about 45 days (Figure 5A). The rate of rooting from cutting was up to 90% (data not shown). This method further shortened the growth cycle of C. equisetifolia to a large extent and was, therefore, an important achievement keeping in view plant transformation and subsequent regeneration.

FIGURE 5
www.frontiersin.org

Figure 5. Regenerating system for Casuarina equisetifolia by cutting. (A) Whole process of regenerating system by cutting. (B) Rooting of branchlet in the preferred root medium. (C) Branchlets with young roots were transferred to the basal medium supplemented with activated carbon.

Determination of Hygromycin and Carbenicillin Concentration for the Selection of Transformants

In order to determine the effective concentration for antibiotic selection, stem segment explants were grown in the callus induction medium supplemented with a series of hygromycin concentrations (0.5, 1, 3, 5, and 10 mg L–1). The results showed that the rate of callus induction declined with increasing in concentration of hygromycin, and reached the steady state at hygromycin a concentration of 5 mg L–1 (Figure 6A). Hence, 5 mg L–1 was determined to be optimum concentration of hygromycin for the selection of transgenic calluses.

FIGURE 6
www.frontiersin.org

Figure 6. Factors that affect transformation efficiency in Casuarina equisetifolia. (A) Hygromycin concentration for the selection of transformants. (B) Carbenicillin concentration for the inhibition of Agrobacterium. (C) Rate of callus and bud induction in the medium with different concentrations of acetosyringone (As) and hygromycin. The results are presented as the mean and standard error from three independent experiments. (D) Co-cultivation duration for transformation frequency. A0, A5, A10, and A20: concentration of As are 0, 5, 10, and 20 mg/L; H0, H5, and H10: concentration of hygromycin were 0, 5, and 10 mg/L.

In addition, we tested different concentrations of carbenicillin (Car), which was required to remove residual Agrobacterium GV3101 from the calluses. Various concentrations of carbenicillin (50, 100, 150, 200, 250, and 300 mg L–1) were supplemented into the Agrobacterium liquid medium for callus growth. We found that 100 mg L–1 carbenicillin was effective in inhibiting Agrobacterium growth (Figure 6B). In the case of stem explants, 250 mg L–1 carbenicillin was used to inhibit Agrobacterium without any effect on callus proliferation and bud regeneration.

Effect of Acetosyringone on Transformation

Acetosyringone is an inducer that activates Virgene on Ti plasmid and regulates T-DNA transfer during Agrobacterium-mediated plant transformation (Turk et al., 1994). It is useful to add As into a cocultivation medium for improving transformation efficiency (Rashid et al., 2011). In our study, we tested the callus induction rate and budding induction rate in a selected medium supplemented with different concentrations of As and hygromycin, including As 0 + Hyg 0 (without As and hygromycin, as control), As 0 + Hgy 5 mg L–1, As 0 + Hgy 10 mg L–1, As 5 mg L–1 + Hgy 5 mg L–1, As 5 mg L–1 + Hgy 10 mg L–1, As 10 mg L–1 + Hgy 5 mg L–1, As 10 mg L–1 + Hgy 10 mg L–1, As 20 mg L–1 + Hgy 5 mg L–1, and As 20 mg L–1 + Hgy 10 mg L–1. Compared with the control and other combinations, the combination with As 20 mg L–1 and Hgy 5 mg L–1 was found to be the best for transgenic callus and budding induction (Figure 6C).

Determination of Co-cultivation Duration

In addition to the above-mentioned factors, co-cultivation duration (1–3 days), was also optimized to improve the frequency of transformation. Agrobacterium GV3101 harboring the vector pCAMBIA1301 (with GUS as a reporter gene) was used for transformation. We found that up to 80% of the calluses were surrounded by Agrobacterium after 3 days of co-cultivation, and that most of the calluses grew well after 1 or 2 days (Figure 6D). Therefore, 2 days of co-cultivation was determined to be suitable for C. equisetifolia transformation.

Confirmation of Transgenic Plants by β-Glucuronidase Staining and Reverse Transcription-PCR

To confirm the positive transgenic plants, hygromycin resistant tissues in different stages, such as callus, bud, root, and seedling, were stained using a GUS buffer solution. A strong GUS staining signal was detected in different tissues of transgenic plants, but not in those of control plants (Figures 7A–D and Supplementary Figure 2). Furthermore, gene presence was also confirmed by reverse transcription-PCR (RT-PCR) amplification of the GUS reporter gene using GUS-specific primers and genomic DNA as template. A 625-bp band was detected only in transgenic lines but not in untransformed control plants (Figure 7E). These data confirmed the successful integration and expression of a foreign gene into the genome of C. equisetifolia.

FIGURE 7
www.frontiersin.org

Figure 7. β-Glucuronidase (GUS) expression in transgenic Casuarina equisetifolia plants. (A–D) GUS staining was observed in the (A) callus, (B) adventitious bud, (C) root, and (D) transgenic plantlet. (E) Analysis of GUS transcriptional levels in pCAMBIA1301 plasmid (P), non-transformed plants (N) and transgenic lines by RT-PCR analysis.

In summary, to get transgenic plants, explants were used to induce calluses in the C4 medium for 10–15 days at 28°C in the dark. Then, good-quality calluses were infected with Agrobacterium (GV3101) carrying the CaMV 35S:GUS binary vector at OD600 ≈0.8 for 2 h and co-cultivated in a co-cultivation medium for 2 days. Next, the infected calluses were transferred to a medium containing 0.5 mg L–1 6-BA, 0.1 mg L–1 TDZ, 250 mg L–1 Car, and 5 mg L–1 Hyg for resistant callus proliferation. After 2–4 weeks, a selection medium with 0.5 mg L–1 6-BA, 250 mg L–1 Car, and 5 mg L−1 Hyg was used for adventitious bud growth, which may last for 1–2 month. Afterward, buds with more than 1-cm length were shifted to the rooting medium supplemented with 0.02 mg L–1 IBA, 0.4 mg L–1 IAA, 250 mg L–1 Car, and 5 mg L–1 Hyg. Roots appeared in about 1–2 weeks after shifting the buds to the rooting medium. For better growth, the seedlings were planted into the basal medium supplemented with 0.5 g L–1 activated carbon, 5 mg L–1 Hyg and 250 mg L–1 Car to improve root growth (Figure 8 and Supplementary Table 2). When the height of transgenic plantlets reached about 10 cm, they were transplanted in the greenhouse at 28°C with an 8-h light/16-h dark photocycle. Besides, cuttings were also used for rapid regeneration of the transgenic plants.

FIGURE 8
www.frontiersin.org

Figure 8. Stepwise protocol for transforming Casuarina equisetifolia using stem segments as the explant.

Establishment of a Transgenic System for Casuarina equisetifolia Through Hairy Roots

An Agrobacterium-mediated transformation system has been the most successful and widespread method for plant genetic engineering in recent decades (Gelvin, 2003a; Vain, 2007; Matveeva and Lutova, 2014). An efficient root transgenic system has previously been reported for the major woody food crop pigeon pea (Meng et al., 2019). Owing to the superior abiotic stress tolerance of C. equisetifolia, we deemed it necessary to establish a fast, simple, and efficient root transgenic method for this plant. To do this, the vector pCAMBIA1-35S:GFP was introduced into A. rhizogenes strain K599 or C58C1. As shown in Table 1, A. rhizogenes strain K599 exhibited much higher efficiency in transgenic root induction (28.3%) and transformation (95.9%) than strain C58C1 (0%). The schematic diagram of the transgenic system for C. equisetifolia by A. rhizogenes-mediated transformation is shown in Figure 9A. C. equisetifolia seedlings of 5–8 cm height were scratched and infected using an Agrobacterium suspension (Figure 9Bi). After about 2 weeks of growth, obvious callus growth was observed around the wound site which then gradually expanded (Figure 9Bii). About 2 weeks later, small hairy roots grew from the calluses (Figures 9Biii–vi). Then, the original roots were cut off after transgenic hairy roots fully developed, and the seedlings were transplanted in a pot (Figures 9Bvii,viii). Furthermore, a strong GFP signal was detected in most of the transgenic roots (Figure 9C). It indicated that the GFP gene was successfully transformed into the hairy root of C. equisetifolia. Overall, it took about 2 months to generate the transgenic plants. The transformation protocol optimized in this study would be greatly helpful in studying the function of key stress-responsive genes in C. equisetifolia.

TABLE 1
www.frontiersin.org

Table 1. Transgenic hairy roots of Casuarina equisetifolia using dipping cut stem in bacterial suspension.

FIGURE 9
www.frontiersin.org

Figure 9. Hairy root transgenic system in Casuarina equisetifolia. (A) Schematic diagram of transgenic system for C. equisetifolia by Agrobacterium rhizogenes-mediated transformation. (B) Process of callus and hairy root regeneration. (i) The 5–8 cm height seedlings for infection; (ii–vii) different stages of callus and hairy root regeneration. (vi) The picture before cutting off the original root, and (vii) the picture of cutting off the original root after transgenic hairy roots fully developed; (viii) seedling with fully developed transgenic hairy roots. (C) GFP signal in transgenic hairy roots.

Discussion

Due to its potential for revegetation under arid conditions, C. equisetifolia has got considerable attention for research. Moreover, its whole genome has recently been sequenced (Ye et al., 2019). Therefore, it is important to develop a simple, rapid, and efficient genetic transformation system for C. equisetifolia, so the function of its key genes in vivo could be studied, especially under stress conditions such as drought, cold, and saline-alkali. Although some Casuarina species, such as C. glauca and Allocasuarina verticillata Lam. have been transformed successfully (Le et al., 1996; Franche et al., 1997; Santi et al., 2003), differences in these species limit the application of the protocols for the transformation of C. equisetifolia, which leads to impediment in the utilization of its genomic resources. Although it has been reported that GUS or GFP expression was observed in C. equisetifolia calluses and young shoots, rooted transgenic plants were not obtained successfully (Zhong et al., 2011). In this study, a novel regeneration and two efficient transformation systems were successfully developed for C. equisetifolia. From the explants to the development of transgenic plants using a A. tumefaciens-mediated (GV3101) method, it took about 4.5 months in total. On the other hand, plant transformation using the Agrobacterium rhizogenes-mediated (K599) hairy root development system took around 45 days to generate transgenic plants, which fairly shortened the time to produce transgenic plants and saved time for studying tree biology at the molecular level.

Plant Growth Regulators Exhibit Significant Effects on the Efficiency of Regeneration

Plant growth regulators not only perform a significant role in plant developmental processes but also play an important part in plant tissue culture. To establish an efficient regeneration system for C. equisetifolia, different kinds of auxin or cytokinin, such as NAA, IAA, IBA, 6-BA, TDZ, and KT, and various combinations of those were used in the basal medium. In this study, we found that cytokinin TDZ was more effective because of its higher callus induction efficiency (Figure 1B). Moreover, TDZ, in combination with 6-BA, could significantly decrease the rate of callus browning compared to its combination with NAA (Figure 1B). The best condition for adventitious bud induction in our study was the basal medium containing only 0.5 mg L–1 6-BA, in which bud induction rate exceeded 95% (Figure 2F). IBA, a main auxin, is widely used to induce adventitious roots in many plant species (Srivastava, 2002; Ludwig-Müller et al., 2005). In our research, the basal medium with 0.02 mg L–1 IBA and 0.4 mg L–1 IAA was found to be effective in C. equisetifolia root induction (Figures 3A,B).

Following genetic transformation with the strain GV3101 harboring vector pCAMBIA1301, several factors were evaluated and optimized in our study. In previous studies, 200–250 mg/L cef has been used (Le et al., 1996; Franche et al., 1997; Santi et al., 2003). Here, 250 mg L–1 Car was added into the selection medium rather than Cef. This is because the green adventitious buds turned to brown after growth of about 1 month in the selection medium with Cef (250 mg L–1) (data not shown). Co-culture duration is also an important factor altering transformation efficiency by influencing Agrobacterium-plant cell interactions, and transformation efficiency for different species benefit from different co-cultivation (Le et al., 1996; Li et al., 2017; Song et al., 2020). In this study, 2 days of co-cultivation was found to be the best for infected callus growth to improve transformation efficiency (Figure 6D).

Acetosyringone, a phenolic compound, induces vir gene expression and promotes the efficiency of plant genetic transformation (Stachel et al., 1985, 1986; Nonaka et al., 2008). Our results demonstrated that the addition of As into the Agrobacterium suspension and co-cultivation medium obviously improved the rate of callus and bud induction (Figure 6C). The development of an efficient genetic transformation system will facilitate physiological and molecular biology studies on C. equisetifolia. Thus, it is important to shorten the time required for the development of C. equisetifolia transgenic plants. A. rhizogenes-mediated root transformation not only saves time for transgenic plant production, but also provides an efficient way to assess gene function. As each transgenic root originates from a single cell in A. rhizogenes-mediated root transformation, a great number of transformants can be obtained (Bercetche et al., 1987; Choi et al., 2004) and analyzed in a relatively short period of time (Kereszt et al., 2007). The Agrobacterium strain is one of the factors that influence the process of T-DNA delivery from Agrobacterium into plant cells (Crane et al., 2006). Again, transgenic root induction efficiency and transformation efficiency using K599 reached up to 28.3 and 95.9%,which were much higher than those using C58C1 (Table 1). Moreover, it took just about 45 days to obtain the transgenic plants. Thus, this system was highly efficient and convenient for abundant production of transgenic C. equisetifolia plants in a short time.

In summary, we established and optimized the protocol for regeneration and two different systems for genetic transformation of C. equisetifolia. Because of higher efficiency and simple propagation, these transformation systems can be used to produce C. equisetifolia transgenic lines with a variety of genes, which may be conductive to study its unique traits, especially tolerance to drought, salinity, and saline-alkali. Moreover, the establishment of these systems could make it possible to improve its cold stress tolerance by genetic engineering.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Author Contributions

HR: conceptualization, methodology, formal analysis, investigation, writing—original draft, review and editing, and funding acquisition. YX: methodology, formal analysis, investigation, and writing—review and editing. XZ: methodology, formal analysis, and investigation. YZ: methodology and formal analysis. JH: writing—review and editing. FC: investigation. GQ and SL: conceptualization, writing—review and editing, and funding acquisition. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by the State Key Laboratory of Subtropical Silviculture (Grant No: ZY20190201), the Natural Science Foundation of Zhejiang Province (Grant No: LY20C020001), the Open Project Program of State Key Laboratory of Crop Stress Biology for Arid Areas, NWFU, Yangling, Shanxi, China (Grant No: CSBAAKF2021008), and the Zhejiang Science and Technology Major Program on Agricultural New Variety Breeding (Grant No: 2016C02056-1).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.784566/full#supplementary-material

References

Benson, D., and Silvester, W. (1993). Biology of Frankia strains, actinomycete symbionts of actinorhizal plants. Microbiol. Rev. 57, 293–319. doi: 10.1128/mr.57.2.293-319.1993

PubMed Abstract | CrossRef Full Text | Google Scholar

Bercetche, J., Chriqui, D., Adam, S., and David, C. (1987). Morphogenetic and cellular reorientations induced by Agrobacterium rhizogenes (strains 1855, 2659 and 8196) on carrot, pea and tobacco. Plant Sci. 52, 195–210. doi: 10.1016/0168-9452(87)90053-7

CrossRef Full Text | Google Scholar

Choi, P. S., Kim, Y. D., Choi, K. M., Chung, H. J., Choi, D. W., and Liu, J. R. (2004). Plant regeneration from hairy-root cultures transformed by infection with Agrobacterium rhizogenes in Catharanthus roseus. Plant Cell Rep. 22, 828–831. doi: 10.1007/s00299-004-0765-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Crane, C., Wright, E., Dixon, R. A., and Wang, Z.-Y. (2006). Transgenic Medicago truncatula plants obtained from Agrobacterium tumefaciens-transformed roots and Agrobacterium rhizogenes-transformed hairy roots. Planta 223, 1344–1354. doi: 10.1007/s00425-006-0268-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Diagne, N., Arumugam, K., Ngom, M., Nambiar-Veetil, M., Franche, C., Narayanan, K. K., et al. (2013a). Use of Frankia and actinorhizal plants for degraded lands reclamation. Biomed. Res. Int. 2013:948258. doi: 10.1155/2013/948258

PubMed Abstract | CrossRef Full Text | Google Scholar

Diagne, N., Diouf, D., Svistoonoff, S., Kane, A., Noba, K., Franche, C., et al. (2013b). Casuarina in Africa: distribution, role and importance of arbuscular mycorrhizal, ectomycorrhizal fungi and Frankia on plant development. J. Environ. Manage. 128, 204–209. doi: 10.1016/j.jenvman.2013.05.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Diem, H. G., and Dommergues, Y. R. (1990). “16 - Current and potential uses and management of Casuarinaceae in the tropics and subtropics,” in The Biology of Frankia and Actinorhizal Plants, eds C. R. Schwintzer and J. D. Tjepkema (Cambridge: Academic Press), 317–342. doi: 10.1016/B978-0-12-633210-0.50021-6

CrossRef Full Text | Google Scholar

Franche, C., Diouf, D., Le, Q. V., Bogusz, D., N’Diaye, A., Gherbi, H., et al. (1997). Genetic transformation of the actinorhizal tree Allocasuarina verticillata by Agrobacterium tumefaciens. Plant J. 11, 897–904. doi: 10.1046/j.1365-313X.1997.11040897.x

CrossRef Full Text | Google Scholar

Gelvin, S. B. (2003a). Agrobacterium-mediated plant transformation: the biology behind the ‘Gene-Jockeying’ tool. Microbiol. Mol. Biol. Rev. 67, 16–37. doi: 10.1128/MMBR.67.1.16-37.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Gelvin, S. B. (2003b). Improving plant genetic engineering by manipulating the host. Trends Biotechnol. 21, 95–98. doi: 10.1016/S0167-7799(03)00005-2

CrossRef Full Text | Google Scholar

He, G. P., Zhuo, R. Y., Chen, Y. C., Huang, Y. Q., and Peng, L. Q. (2011). Effect of low temperature on physiological index of cold-tolerant Casuarina glauca clones. For. Res. 24, 523–526.

Google Scholar

Holsters, M., de Waele, D., Depicker, A., Messens, E., van Montagu, M., and Schell, J. (1978). Transfection and transformation of Agrobacterium tumefaciens. Mol. Gen. Genet. 163, 181–187. doi: 10.1007/BF00267408

PubMed Abstract | CrossRef Full Text | Google Scholar

Hwang, H.-H., Yu, M., and Lai, E.-M. (2017). Agrobacterium-mediated plant transformation: biology and applications. Arabidopsis Book 15:e0186. doi: 10.1199/tab.0186

PubMed Abstract | CrossRef Full Text | Google Scholar

Jefferson, R. A., Kavanagh, T. A., and Bevan, M. W. (1987). GUS fusions: beta-glucuronidase as a sensitive and versatile gene fusion marker in higher plants. EMBO J. 6, 3901–3907. doi: 10.1002/j.1460-2075.1987.tb02730.x

CrossRef Full Text | Google Scholar

Kereszt, A., Li, D., Indrasumunar, A., Nguyen, C. D., Nontachaiyapoom, S., Kinkema, M., et al. (2007). Agrobacterium rhizogenes-mediated transformation of soybean to study root biology. Nat. Protoc. 2, 948–952. doi: 10.1038/nprot.2007.141

PubMed Abstract | CrossRef Full Text | Google Scholar

Le, Q. V., Bogusz, D., Gherbi, H., Lappartient, A., Duhoux, E., and Franche, C. (1996). Agrobacterium tumefaciens gene transfer to Casuarina glauca, a tropical nitrogen-fixing tree. Plant Sci. 118, 57–69. doi: 10.1016/0168-9452(96)04386-5

CrossRef Full Text | Google Scholar

Li, H. B., Li, N., Yang, S. Z., Peng, H. Z., Wang, L.-L., Wang, Y., et al. (2016). Transcriptomic analysis of Casuarina equisetifolia L. in responses to cold stress. Tree Genet. Genomes 13:7. doi: 10.1007/s11295-016-1090-z

CrossRef Full Text | Google Scholar

Li, S., Zhen, C., Xu, W., Wang, C., and Cheng, Y. (2017). Simple, rapid and efficient transformation of genotype Nisqually-1: a basic tool for the first sequenced model tree. Sci. rep. 7, 2638. doi: 10.1038/s41598-017-02651-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ludwig-Müller, J., Vertocnik, A., and Town, C. (2005). Analysis of indole-3-butyric acid-induced adventitious root formation on Arabidopsis stem segments. J. Exp. Bot. 56, 2095–2105. doi: 10.1093/jxb/eri208

PubMed Abstract | CrossRef Full Text | Google Scholar

Matveeva, T. V., and Lutova, L. A. (2014). Horizontal gene transfer from Agrobacterium to plants. Front. Plant Sci. 5:326. doi: 10.3389/fpls.2014.00326

PubMed Abstract | CrossRef Full Text | Google Scholar

Meng, D., Yang, Q., Dong, B., Song, Z., Niu, L., Wang, L., et al. (2019). Development of an efficient root transgenic system for pigeon pea and its application to other important economically plants. Plant Biotechnol. J. 17, 1804–1813. doi: 10.1111/pbi.13101

PubMed Abstract | CrossRef Full Text | Google Scholar

Nonaka, S., Yuhashi, K., Takada, K., Sugaware, M., Minamisawa, K., and Ezura, H. (2008). Ethylene production in plants during transformation suppresses vir gene expression in Agrobacterium tumefaciens. New Phytol. 178, 647–656. doi: 10.1111/j.1469-8137.2008.02400.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Obertello, M., Santi, C., Sy, M. O., Laplaze, L., Auguy, F., Bogusz, D., et al. (2005). Comparison of four constitutive promoters for the expression of transgenes in the tropical nitrogen-fixing tree Allocasuarina verticillata. Plant Cell Rep. 24, 540–548. doi: 10.1007/s00299-005-0963-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Ozyigit, I., Dogan, I., and Artam Tarhan, E. (2013). “Agrobacterium rhizogenes-mediated transformation and its biotechnological applications in crops,” in Crop Improvement, eds K. R. Hakeem, P. Ahmad, and M. Ozturk (Berlin: Springer Verlag), 1–48. doi: 10.1007/978-1-4614-7028-1_1

CrossRef Full Text | Google Scholar

Pan, Y., Li, Y., and Tan, T. (1996). Casuarina provenance test. For. Res. 9, 138–145.

Google Scholar

Pinyopusarerk, K., and Williams, E. R. (2000). Range-wide provenance variation in growth and morphological characteristics of Casuarina equisetifolia grown in Northern Australia. For. Ecol. Manage. 134, 219–232. doi: 10.1016/S0378-1127(99)00260-1

CrossRef Full Text | Google Scholar

Potgieter, L. J., Richardson, D. M., and Wilson, J. R. (2014). Casuarina cunninghamiana in the Western Cape, South Africa: determinants of naturalisation and invasion, and options for management. S. Afr. J. Bot. 92, 134–146. doi: 10.1016/j.sajb.2014.02.013

CrossRef Full Text | Google Scholar

Rashid, H., Chaudhry, Z., and Khan, M. H. (2011). Effect of explant plant source and acetosyringone concentration on transformation efficiency of wheat cultivars. Afr. J. Biotechnol. 10, 8737–8740. doi: 10.5897/ajb10.770

CrossRef Full Text | Google Scholar

Santi, C., Svistoonoff, S., Constans, L., Auguy, F., Duhoux, E., and Bogusz, D. (2003). Choosing a reporter for gene expression studies in transgenic actinorhizal plants of the Casuarinaceae family. Plant Soil 254, 229–237. doi: 10.1007/978-94-017-1601-7_24

CrossRef Full Text | Google Scholar

Song, G. Q., Han, X., Wiersma, A. T., Zong, X., Awale, H. E., and Kelly, J. D. (2020). Induction of competent cells for Agrobacterium tumefaciens-mediated stable transformation of common bean (Phaseolus vulgaris L.). PLoS One 15:e0229909. doi: 10.1371/journal.pone.0229909

PubMed Abstract | CrossRef Full Text | Google Scholar

Srivastava, L. M. (2002). Plant Growth And Development: Hormones And Environment. Cambridge: Academic press. doi: 10.1111/j.1755-0238.2010.00083_8.x

CrossRef Full Text | Google Scholar

Stachel, S. E., Messens, E., Van Montagu, M., and Zambryski, P. (1985). Identification of the signal molecules produced by wounded plant cells that activate T-DNA transfer in Agrobacterium tumefaciens. Nature 318, 624–629. doi: 10.1038/318624a0

CrossRef Full Text | Google Scholar

Stachel, S. E., Nester, E. W., and Zambryski, P. C. (1986). A plant cell factor induces Agrobacterium tumefaciens gene expression. Proc. Natl. Acad. Sci. U. S. A. 83, 379–383. doi: 10.1073/pnas.83.2.379

PubMed Abstract | CrossRef Full Text | Google Scholar

Tani, C., and Sasakawa, H. (2003). Salt tolerance of Casuarina equisetifolia and Frankia Ceq1 strain isolated from the root nodules of C. equisetifolia. Soil Sci. Plant Nutr. 49, 215–222. doi: 10.1080/00380768.2003.10410000

CrossRef Full Text | Google Scholar

Turk, S., Lange, R., Regensburg-Tuïnk, T. J., and Hooykaas, P. (1994). Localization of the VirA domain involved in acetosyringone-mediated vir gene induction in Agrobacterium tumefaciens. Plant Mol. Biol. 25, 899–907. doi: 10.1007/BF00028884

PubMed Abstract | CrossRef Full Text | Google Scholar

Vain, P. (2007). Thirty years of plant transformation technology development. Plant Biotechnol. J. 5, 221–229. doi: 10.1111/j.1467-7652.2006.00225.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, C., Zhang, Y., Ma, N., and Zhong, C. I. (2012). Response characteristics of Casuarina equisetifolia inoculated with mycorrhizal fungi under low temperature stress. Xibei Zhiwu Xuebao 32, 2068–2074.

Google Scholar

Xu, C., Bing, Z., Nan, G., Wang, X., and Wang, Y. (2006). Effect of the plant hormones and active carbon on rooting of in vitro shoot of saussurea involucrata kar. et kir. Chin. Agricult. Sci. Bull. 22, 41–43.

Google Scholar

Ye, G., Zhang, H., Chen, B., Nie, S., Liu, H., Gao, W., et al. (2019). De novo genome assembly of the stress tolerant forest species Casuarina equisetifolia provides insight into secondary growth. Plant J. 97, 779–794. doi: 10.1111/tpj.14159

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, L.-H., Lin, Y.-M., Ye, G.-F., Liu, X.-W., and Lin, G.-H. (2008). Changes in the N and P concentrations, N:P ratios, and tannin content in Casuarina equisetifolia branchlets during development and senescence. J. For. Res. 13, 302–311. doi: 10.1007/s10310-008-0081-9

CrossRef Full Text | Google Scholar

Zhao, Y. (1999). The use of the active carbon in the invitro tissue culture of pinellia lernatal. J. Changwei Teach. Coll. 18, 15–17.

Google Scholar

Zhong, C., Bai, J., and Zhang, Y. (2005). Introduction and conservation of casuarina trees in China. For. Res. 18, 345–350.

Google Scholar

Zhong, C., Jiang, Q., Zhang, Y., Nambiarveetil, M., Vaissayre, V., Bogusz, D., et al. (2011). “Genetic transformation of Casuarina equisetifolia by Agrobacterium tumefaciens,” in Improving Smallholder Livelihoods Through Improved Casuarina Productivity Proceedings of International Casuarina Workshop, eds C. Zhong, K. Pinyopusarerk, A. Kalinganire, and F. Claudine (Pékin: China Forestry Publishing House), 163–169.

Google Scholar

Zhong, C., Zhang, Y., Chen, Y., Jiang, Q., Chen, Z., Liang, J., et al. (2010). Casuarina research and applications in China. Symbiosis 50, 107–114. doi: 10.1007/s13199-009-0039-5

CrossRef Full Text | Google Scholar

Keywords: Casuarina equisetifolia, regeneration system, transformation system, Agrobacterium tumefaciens, Agrobacterium rhizogenes

Citation: Ren H, Xu Y, Zhao X, Zhang Y, Hussain J, Cui F, Qi G and Liu S (2022) Optimization of Tissue Culturing and Genetic Transformation Protocol for Casuarina equisetifolia. Front. Plant Sci. 12:784566. doi: 10.3389/fpls.2021.784566

Received: 28 September 2021; Accepted: 02 December 2021;
Published: 21 January 2022.

Edited by:

Horacio Esteban Hopp, University of Buenos Aires, Argentina

Reviewed by:

Sandra Sharry, National University of La Plata, Argentina
Gerardo Sánchez, Instituto Nacional de Tecnología Agropecuaria, Argentina

Copyright © 2022 Ren, Xu, Zhao, Zhang, Hussain, Cui, Qi and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Guoning Qi, gnqi@zafu.edu.cn; Shenkui Liu, shenkuiliu@nefu.edu.cn

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.