- 1Microbial Technology Department, CSIR-Central Institute of Medicinal and Aromatic Plants, Lucknow, India
- 2Biotechnology Division, CSIR-Central Institute of Medicinal and Aromatic Plants, Lucknow, India
- 3Analytical Chemistry Department, CSIR-Central Institute of Medicinal and Aromatic Plants, Lucknow, India
Secondary metabolite biosynthesis in medicinal plants is multi-step cascade known to be modulated by associated endophytes. While a single endophyte is not able to upregulate all biosynthetic steps, limiting maximum yield achievement. Therefore to compliment the deficient characteristics in an endophyte we tried consortium of endophytes to achieve maximum yield. Here, efforts were made to maximize the in planta morphine yield, using consortium of two endophytes; SM1B (Acinetobacter sp.) upregulating most of the genes of morphine biosynthesis except T6ODM and CODM, and SM3B (Marmoricola sp.) upregulating T6ODM and CODM in alkaloid-less Papaver somniferum cv. Sujata. Consortium-inoculation significantly increased morphine and thebaine content, and also increased the photosynthetic efficiency of poppy plants resulted in increased biomass, capsule weight, and seed yields compared to single-inoculation. The increment in morphine content was due to the modulation of metabolic-flow of key intermediates including reticuline and thebaine, via upregulating pertinent biosynthetic genes and enhanced expression of COR, key gene for morphine biosynthesis. This is the first report demonstrating the endophytic-consortium complimenting the functional deficiency of one endophyte by another for upregulating multiple genes of a metabolic pathway similar to transgenics (overexpressing multiple genes) for obtaining enhanced yield of pharmaceutically important metabolites.
Introduction
Increasing populations have threatened food security due to the availability of limited land for cultivation. Moreover, changing environmental conditions are also causing a considerable reduction in the crop productivity. It is estimated that by the year 2050, the world population may reach nine billion. Therefore, to satisfy the needs of the growing population diverse approaches have been tried to increase the food production at the same or even higher rates. Developing transgenics has so far received foremost attention by introducing desired traits to produce the desired product and improving the productivity in normal as well as changing environmental conditions including cultivation under biotic and abiotic stresses. Plant yield and tolerance to different environmental stresses are polygenic traits and regulated by multiple genes (Bohnert et al., 2006). Therefore, integration/expression/genetic manipulation of more than one genes in a plant became the preferred choice for improvement of plant yield and performance under stressful environmental conditions. However, these practices have limitations related to production cost, social acceptability, and sustainability (Halpin, 2005). Therefore, besides developing transgenic plants, alternate sustainable approaches need to be explored for achieving maximum plant yield.
Endophytes have emerged as a promising candidate for sustainable agriculture. Endophytes are plant-associated microbes residing within the plant without harming or causing any disease symptoms. They promote plant growth, protect plants from environmental stresses and are the promising source of therapeutic secondary metabolites (Kuklinsky-Sobral et al., 2004; Strobel et al., 2004; Waller et al., 2005; Rodriguez et al., 2008; Staniek et al., 2008; Kusari et al., 2012; Quecine et al., 2012; Otieno et al., 2015; Egamberdieva et al., 2017). Involvement of endophytes in the modulation of secondary metabolite biosynthesis of the host plant has also been established (Pandey et al., 2016a,b, 2018). Bacterial endophyte interaction with E. purpurea enhances its secondary metabolites which contributes to therapeutic properties of this medicinal plant (Maggini et al., 2017). An Artemisia annua endophytic actinobacterium Pseudonocardia strain YIM 63111 has been able to induce artemisinin production by upregulating the genes of artemisinin biosynthetic pathway, i.e., cytochrome P450 monooxygenase and cytochrome P450 oxidoreductase (Li et al., 2012). Also role of microbial community present in the cortical parenchymatous tissue of Vetiver in essential oil biosynthesis has been reported (Del Giudice et al., 2008). Inoculation with bacterial endophytes increases biomass and in planta content of terpenoid indole alkaloids like vindoline in Catharanthus roseus (Tiwari et al., 2013). Paenibacillus strain has been found to strongly influence the plant metabolites of in-vitro grown poplar plants (Scherling et al., 2009). Interestingly, endophytes are also able to produce secondary metabolites analogous to their host plants, e.g., camptothecin and analogs (Shweta et al., 2010), taxol (Soliman et al., 2011), vincristine, and vinblastine (Kumar et al., 2013). Fungal endophytes as symbiotically modified organism (SMO) have also been used to improve grass yield, provide resistance against pests, weeds and herbivores by accumulation of alkaloids and antioxidants (Gundel et al., 2013, 2018). It has been demonstrated in our laboratory that endophytes are involved in the improvement of the content of vindoline (intermediate for the biosynthesis of anticancerous vinblastine and vincristine) in Catharanthus roseus (Pandey et al., 2016a), withanolides in Withania somnifera (Pandey et al., 2018), and benzylisoquinoline alkaloids (BIAs) in Papaver somniferum (Pandey et al., 2016b). Previous studies have also demonstrated the diverse role of endophytes in a host plant (Di Fiore and Del Gallo, 1995; Saikkonen et al., 1998; Singh et al., 2011;Hardoim et al., 2015; Pandey et al., 2016a,b, 2018). Every endophyte has a specific characteristic, up/down-regulating certain specific genes responsible for plant growth, yield, and modulation of secondary metabolite biosynthesis. Although reports signifying the role of a particular endophyte modulating a single functional characteristic like growth promotion or disease resistance, stress tolerance or enhance secondary metabolite production are available, but the information on the possibility of strengthening multiple functional characteristics through the use of consortium is scarce. Therefore, the use of consortium of endophytes having different functional characteristics should be explored to maximize the yields and desired products.
Here, we demonstrated the use of consortium of selected endophytes associated with the capsule of opium poppy plant which is the site for BIA biosynthesis to maximize the yield of pharmaceutically important metabolite morphine. Opium poppy is an important medicinal plant and sole source of therapeutically important BIAs such as morphine, codeine, thebaine, papaverine, noscapine, and sanguinarine. Morphine belongs to the group “opioids” which are key drugs for the alleviation of moderate-to-severe acute pain associated with cancer and other diseases or after any surgery or physical trauma (Rosenblum et al., 2008; Whittle et al., 2011).
Previously we have established that the opium poppy has potential endophytes which can enhance in planta BIA biosynthesis. We demonstrated that different endophytes are associated with different parts of poppy plants performing an important role in a tissue-specific manner. In general, endophytes related with leaf were involved in modulating photosynthetic efficiency and endophytes associated with the capsule were drawn in modulating the BIA-biosynthesis of the poppy plant. It was observed that endophytes SM1B and SM3B isolated from capsule of alkaloid rich Sampada were unable to enhance the BIA production in P. somniferum cv. Sampada plants. However, these capsule-associated endophytes were able to increase the production of BIAs in alkaloid-less cv. Sujata. The cv. Sujata was developed from alkaloid-rich cv. Sampada with the aim to convert a narcotic “opium poppy” (i.e., cv. Sampada) into non-narcotic seed poppy (i.e., cv. Sujata) using gamma rays and EMS (ethyl methane sulphonate) (Sharma et al., 1999). The cv. Sujata is considered to be having the potential for commercial cultivation of seeds due to high nutritive value of 24% protein, >58% oil with high unsaturated fatty acids makes it a safe and potential food crop having protein rich seeds (Sharma et al., 2002). Previously we have established that the individual inoculation with endophytes SM1B and SM3B increased biomass upto 44% in cv. Sujata. SM1B inoculation in cv. Sujata plants increased content of Morphine by 1044%, Papaverine by 349% and Noscapine by 936%, and SM3B inoculation increased Morphine by 37%, Papaverine by 66%, Noscapine by 72%, and Thebiane by 154% (Pandey et al., 2016b). It was observed that endophyte SM1B could upregulate most of the genes of BIA biosynthesis except thebaine 6-O-demethylase (T6ODM) and codeine O-demethylase (CODM) (Pandey et al., 2016b). However, we could also identify another endophyte SM3B that was able to upregulate T6ODM and CODM expression. We, therefore, expected that the consortium of SM1B and SM3B will lead to the overexpression of all BIA biosynthetic pathway genes simultaneously and will produce more morphine than single inoculation. Here, efforts were made to maximize morphine yield using the consortium of two capsule-associated endophytes SM1B (Acinetobacter sp.) and SM3B (Marmoricola sp.) in the poppy plant by upregulating different genes of BIA biosynthesis.
Materials and Methods
Plant Materials and Growth Conditions
Seeds of Papaver somniferum cv. Sujata was obtained from the National Gene Bank for Medicinal and Aromatic Plants at CSIR-CIMAP, Lucknow. Poppy plants were grown in pots filled with 3.5 kg of autoclaved soil and vermicompost mixture (2:1, v/v) having dimension (17 cm height × 22 cm top diameter × 12 cm bottom diameter and 3.7 l volume) and irrigated with sterile water with an ordinary photoperiod in a greenhouse under natural light intensity at 20°C ± 2°C.
Antagonistic Activity of Endophytes
Antagonistic activity of both the selected endophytes was tested by in vitro antimicrobial assay. For this, the disc diffusion method of Bhunia et al. (1988) was followed with some modifications. SM1B and SM3B both were grown in nutrient broth for 24 h at 28°C. The final concentration of each bacterium was adjusted to 106 CFU/ml. 100 μl 106 CFU/ml cell suspension was spread on nutrient agar media and leave the plates to get dried on room temperature (Tajehmiri et al., 2014). The sterile paper disc (6 mm) were placed over nutrient agar media plates, seeded with indicator strains. Hundred microliter of culture was added to the sterile paper disc and incubated at 28°C for 48 h. After incubation, antagonistic activity was observed around the paper disc. As there was no zone of inhibition around the paper disc and both the bacteria were growing with each other which confers that there is no antagonistic activity between the two bacterial strains and they synergistically performed their function.
Endophyte Inoculation
For consortium development individual suspension (1 × 108 CFU ml−1) of the two endophytes (SM1B and SM3B) was prepared in PBS and were mixed at a CFU of 1 × 108ml−1 microbes each in 5 ml culture prepared in PBS. Endophyte-free seeds were treated by putting them in consortium suspension for 2 h and then used for sowing in pots. Pots were kept in the controlled condition in greenhouse. For each treatment the experiment was conducted in triplicates. Initially, there were 25 seeds sown in each pot, and when the seeds grown for 2 weeks. Only five healthy seedlings were kept in each pot. In order to boost up the soil with sufficient inoculum of endophyte, to ensure their presence in the soil. The pots were again inoculated with 10 ml pot−1 of the endophyte suspension (1 × 108 CFU ml−1). To make certain the presence of adequate numbers of endophyte in the soil, pots were inoculated for a second time after 2 weeks of germination of seeds. PBS was used in place of endophyte suspension in case of the control plants (uninoculated). To confirm the presence of specific endophytes, re-isolation of endophyte has been done from various parts of poppy plants at 60 days stage. No endophytes were found to be present at any stage from the uninoculated endophyte-free control plants. SM1B and SM3B inoculated plants (individual inoculation) were grown simultaneously with the consortium-inoculated (SM1B and SM3B) plants.
Colonization of Endophytes
Endophyte colonization was inspected in the capsule of poppy plants in triplicates after 90 days of growth. Tissues were surface sterilized and sterility check was performed as illustrated previously (Pandey et al., 2016b). Under sterile conditions, the surface sterilized tissues were thoroughly macerated, and serial dilutions were plated on nutrient agar with three replications each. The plates were then incubated at 28°C for 24–72 h and then the CFU were estimated. The obtained colonies were confirmed by the morphological characteristics followed by 16S rRNA sequencing.
Photosynthetic Pigment Content, Photosynthetic Rate, Stomatal Conductance, Transpiration Rate
Photosynthetic pigment (Chlorophyll) content and photosynthetic efficiency of fully expanded leaves (third from the top) of 60 days poppy plants (n = 6) were measured. Chlorophyll extraction was performed in chilled methanol, and the content was estimated (Lichtenthaler and Buschmann, 2001). Photosynthetic efficiency (net CO2 assimilation, transpiration rate and stomatal conductance) was measured in the attached leaves (third from the top) using the Portable Photosynthesis System (CIRAS-3, PP System, United States) attached with the Chlorophyll fluorescence module (CFM-3). For photosynthesis measurement, CO2 concentration in leaf cuvette was maintained at ambient CO2 (400 ppm), temperature (25°C) and 400 μmol photons m−2 s−1 light.
Biomass, Capsule Weight, and Seed Yield
Biomass of 60 days old poppy plant (n = 6) was measured by harvesting the entire shoots and dried at 70°C for 5 days, and their dry weights were measured. Capsule weight and seed yield were measured from 120 days mature poppy plants.
Benzylisoquinoline Alkaloids
BIAs analysis was performed from the capsules of 120 days old poppy plants. Capsules were completely dried and ground to fine powder. Powdered samples (1 g) were extracted three times in 100% methanol (30 ml) with 20 min sonication at 42°C. Extract was filtered through Whatman Number 1 filter paper and then the filtrates were concentrated in Rotavapour, dried, and re-dissolved in methanol. BIAs (morphine, codeine, thebaine, papaverine, noscapine and reticuline) were analyzed by isocratic reverse-phase HPLC using phosphate buffer (0.1 M NaH2PO4.2H2O, pH 3.5, 77%, v/v) and acetonitrile (23%, v/v), sample injection volume of 20 μl, at flow rate of 1.0 ml min−1 on Phenomenex, Luna C18 column®, Waters (4.6 × 250 mm, 5 μm). The analysis was performed with Empower Pro software (Waters, United States). Photodiode array data acquisition was in the range of 200–400 nm and the quantitation was done at 240 nm. Peaks of analyzed alkaloids were identified by comparing their retention times with commercially available authentic standards. Noscapine and reticuline were purchased from Sigma-Aldrich, St. Louis, United States. Morphine, codeine, thebaine and papaverine were gifted from CSIR-National Botanical Research Institute, Lucknow. The data was an average of two independent quantifications repeated in three biological replicates (n = 6).
qRT-PCR Analysis
Total RNA was isolated from the green capsule of un-inoculated endophyte-free control and endophytes-inoculated 90 days old poppy plants by using TRI-reagent (Sigma-Aldrich). RNA quantification was done using NanoDrop 1000 spectrophotometer (Thermo Fisher Scientific). Genomic-DNA contamination in the total RNA preparation was eliminated using RNase-free DNase I enzyme (Thermo Scientific). First-strand cDNA synthesis was performed with ∼5 μg of total RNA using the RevertAid First Strand cDNA Synthesis Kit with oligo(dT)18 primer (Thermo Scientific) following the manufacturer’s instructions. Previously described gene-specific primers were used for the analysis of relative quantification of 16 biosynthetic gene transcripts involved in BIA biosynthetic pathway (Pandey et al., 2016b). qRT-PCR was performed using SYBR-Green I detection on triplicate technical replicates of triplicate biological samples (n = 3) on an Applied Biosystems StepOnePlusTM Real-Time PCR System as described previously (Pandey et al., 2016b). PCR mixture contains 1 μl of 10 times diluted cDNA synthesis reaction, 300 nM each forward and reverse primers and 5 μl Power SYBR Green PCR Master Mix (Applied Biosystems), in a 10 μl reaction volume. PCR conditions were maintained at 95°C for 10 min, followed by 40 cycles of denaturation at 95°C for 15 s each and annealing at 60°C for 1 min each. The intensities of fluorescent signal were recorded on an Applied Biosystems StepOnePlusTM Real Time PCR System and analysis were carried out using StepOne software (Applied Biosystems). Specificity of RT-qPCR was evaluated by using the dissociation method (Applied Biosystems) subjecting all amplicons to a melt-curve analysis. The threshold cycle (Ct) for each gene was normalized against the Ct for β-actin (EB740770) of opium poppy, which was used as the endogenous reference transcript. Further, Mean Ct values were calculated from technical triplicates, and the relative levels of transcripts of different genes. BIAs biosynthetic genes in endophyte-inoculated plants were compared with the calibrator (non-inoculated control plant) using the relative quantification 2−ΔΔCt method and efficiency of PCR was also determined by the method of Livak and Schmittgen (2001) and Schmittgen and Livak (2008). For amplification serially diluted cDNA was used as template for different target genes and reference gene. To determine ΔCt (Ct target – Ct reference) the average Ct was calculated for both target gene and reference gene. A plot of the log cDNA dilution versus ΔCt was drawn, and the absolute value of the slope were found to be near zero for all the target genes that showed that the efficiencies of the target and reference genes were quite similar. Therefore, relative quantification 2−ΔΔCt method was applied to analyze the data (Livak and Schmittgen, 2001).
Statistical Analysis
Statistical analysis of data was carried out following ANOVA analysis of variance, suitable to completely randomized design (CRD), using SPSS (IBM, Chicago, IL, United States). Comparisons between means were carried out using Duncan’s multiple range tests at a significance level of P ≤ 0.05.
Results
Antagonistic Activity of Endophytes and Their Colonization in the Capsule of Poppy Plants
Antagonistic activity of SM1B and SM3B was tested. It was observed that both endophytes grew well together and did not inhibit the growth of each other (Supplementary Figure S1). Colonization of SM1B and SM3B endophytes was determined in the capsule of poppy plant inoculated with SM1B and SM3B individually as well as a consortium of SM1B and SM3B. It was observed that both endophytes could successfully colonize the capsule of P. somniferum plant (Figure 1). Colonization of SM1B was reduced in the poppy plants inoculated with a consortium of SM1B and SM3B compared to individual inoculation of SM1B. However, the colonization of SM3B was found to be similar in both individual (SM3B), and consortium inoculated plants (Supplementary Table S1).
Figure 1. Colonization of endophytes in the capsule of Papaver somniferum plants. CFU- colony forming unit. Values are the mean ± SE of three replicates (n = 3). Different letters (a–b) indicate statistically significant differences between treatments (Duncan’s multiple range test P < 0.05).
Effect of a Consortium of Endophytes on Photosynthetic Efficiency of the Plant
Photosynthetic efficiency of plants was studied by measuring the chlorophyll content, net CO2 assimilation rate, stomatal conductance, and transpiration rate. Application of endophytic consortium improved the measured parameters. Inoculation with the consortium of SM1B and SM3B in P. somniferum plants increased chlorophyll content, net CO2 assimilation rate, stomatal conductance, and transpiration rate by 27.1, 25.8, 39.3, and 46.5%, respectively, compared to endophyte-free control plants (Table 1). Individual inoculation of SM1B could increase the chlorophyll (21.43%), net CO2 assimilation rate (15.63%), stomatal conductance (30.36%), and transpiration rate (21.81%), however, SM3B inoculation could only increase the chlorophyll content (17.14%).
Table 1. Effect of inoculation of consortium of endophytes on different physiological parameters of Papaver somniferum plant.
Effect of a Consortium of Endophytes on Biomass, Capsule Weight, and Seed Yield
Application of consortium of endophytes significantly enhanced the biomass, capsule weight and seed yield of P. somniferum plants. Inoculation with a consortium of SM1B and SM3B increased the biomass of P. somniferum plant by 75% compared to endophyte-free control plants (Table 1). Consortium inoculated plants also had 66.4% higher capsule weight than that of endophyte-free control plants. More importantly, the consortium of endophytes could enhance the yield of seeds, another economic part of the plant, by 78.8% over control (Table 1). Individual inoculation (SM1B) could only increase the biomass by 31.7%, capsule weight by 30.3% and seed yield by 26.9% which was considerably lower than the consortium inoculated plants. However, SM3B inoculation could only increase the biomass (13.3%) (Table 1).
Effect of a Consortium of Endophytes on Alkaloid Content
To study the impact of a consortium of endophytes on BIA production, the content of morphine, codeine, thebaine, papaverine, and noscapine was estimated in the capsule of individual, consortium inoculated and non-inoculated endophyte-free control poppy plants (Figure 2). Consortium inoculation increased the content of morphine content (2250%) than that of endophyte-free control plants. Consortium inoculation also increased papaverine and noscapine content by 36.4 and 53.3%, respectively.
Figure 2. Effect of endophytes inoculation on the alkaloid content of Papaver somniferum. Morphine (A), codeine (B), thebaine (C), papaverine (D), and noscapine (E) content was measured from the dried capsule of 120 days old P. somniferum plants inoculated with endophytes SM1B and SM3B individually, and in the form of consortium (SM1B+SM3B). Non-inoculated endophyte free plants were used as a control. The data was an average of two independent quantifications repeated in three biological replicates (n = 6) and represented in % content/gram of dry weight (DW). Error bars represent standard errors. Different letters (a–d) indicate statistically significant differences between treatments (Duncan’s multiple range test P < 0.05).
Individual inoculation of SM1B and SM3B could increase the morphine (200–1150%) which was considerably lower than consortium inoculated plants; also an increase in the content of papaverine (27.3–263.6%) and noscapine (226.7–886.7%) was noticed over non-inoculated plants. Although consortium inoculated plants, had 1067% higher thebaine than that of endophyte-free control plants, inoculation of SM3B could enhance thebaine content only by 483% whereas inoculation SM1B did not affect the thebaine content. Individual inoculation of both the endophytes as well as consortium, on the other hand, decreased the content of codeine (17.6–98.9%).
Reticuline content was also found to be affected with the endophyte inoculations. Reticuline was not detected in SM1B-inoculated plants. However, SM3B-inoculated plants had a higher accumulation of reticuline compared to endophyte-free control plants (Figure 3). Consortium-inoculated plants had lower reticuline accumulation than that of SM3B, however, content was higher than that of endophyte-free control plants (Supplementary Table S2).
Figure 3. Effect of endophytes inoculation on the reticuline content of Papaver somniferum. Reticuline was measured from the dried capsule of 120 days old P. somniferum plants inoculated with endophytes SM1B and SM3B individually, and in the form of consortium (SM1B+SM3B). Non-inoculated endophyte free plants were used as a control. The data was an average of two independent quantifications repeated in three biological replicates (n = 6) and represented in % content/gram of dry weight (DW). Error bars represent standard errors. Different letters (a–c) indicate statistically significant differences between treatments (Duncan’s multiple range test P < 0.05). nd, not detected.
Expression of Genes of BIA Biosynthesis
To understand the mechanism of endophytes consortium-mediated increase in morphine content the expression of different genes of BIA biosynthesis was quantified using qRT-PCR. Transcript of a total of 16 genes was quantified. Results were normalized to endogenous gene of opium poppy reference transcript (β-actin) and are presented relative to the level in non-inoculated endophyte-free control plants (calibrator). RQ was calculated using the equation; RQ = 2−ΔΔCt. Expression of TyrAT was not affected significantly due to the individual as well as consortium inoculations. The expression of TYDC, NCS, and 6OMT in the consortium-inoculated plants was 4.4-, 5.3-, and 7.6-fold higher compared to non-inoculated endophyte-free control plants, respectively (Figure 4). Contrary, SM1B alone inoculated plants had only 2.5- and 3.2-fold higher expression of TYDC and NCS, respectively, while SM3B inoculation did not affect their expression significantly. 6OMT expression was lower in individual inoculated plants compared to consortium inoculated plants. Expression of CNMT and NMCH in consortium-inoculated plants was 8.9- and 9.2-fold higher than that of endophyte-free control plants, respectively (Figure 4). Individual inoculation of SM1B and SM3B could also increase the expression of CNMT (5.9–7.3-fold) and NMCH (4.2–6.3-fold), but the increment was lower compared to consortium inoculation.
Figure 4. Effect of endophytes inoculation on the expression of genes involved in reticuline biosynthesis. Total RNA was isolated from the green capsule of 90 days old P. somniferum plants inoculated with endophytes SM1B and SM3B individually, and in the form of consortium (SM1B+SM3B), reverse transcribed and used as a template for RT-qPCR with SYBR Green detection. The capsules of non-inoculated endophyte free plants were used as a control. Expression of TyrAT (A), TYDC (B), NCS (C), 6OMT (D), CNMT (E), NMCH (F) was analyzed. Results were normalized to actin (reference transcript) and are shown relative to the level in non-inoculated endophyte free control plants (calibrator). qRT-PCR was performed on triplicate technical replicates of triplicate biological samples (n = 3). Data are means ± SE (n = 3 biological replicates) and Y-axis represents relative quantity (RQ). RQ was calculated using the equation; RQ = 2−ΔΔCt. Different letters (a–d) indicate statistically significant differences between treatments (Duncan’s multiple range test P < 0.05).
Expression of genes involved in papaverine biosynthesis, i.e., N7OMT and 7OMT was increased by 7.2- and 2.6-fold, in consortium-inoculated plants, respectively (Figure 5). Individual inoculation could also upregulate N7OMT and 7OMT expression but were lower than that of consortium inoculated plants. In consortium-inoculated plants expression of genes involved in noscapine biosynthesis, i.e., BBE and TNMT was increased by 5.3- and 7.3-fold, respectively (Figure 5) whereas individual-inoculation could also increase BBE (5.5–11.8-fold) and TNMT (4.5–8-fold) expression.
Figure 5. Effect of endophytes inoculation on the expression of genes involved in papaverine and noscapine biosynthesis. Total RNA was isolated from the green capsule of 90 days old P. somniferum plants inoculated with endophytes SM1B and SM3B individually, and in the form of consortium (SM1B+SM3B), reverse transcribed and used as a template for RT-qPCR with SYBR Green detection. The capsules of non-inoculated endophyte free plants were used as a control. Expression of genes involved in papaverine biosynthesis N7OMT (A), 7OMT (B) and noscapine biosynthesis BBE (C), TNMT (D) was analyzed. Results were normalized to actin (reference transcript) and are shown relative to the level in non-inoculated endophyte free control plants (calibrator). qRT-PCR was performed on triplicate technical replicates of triplicate biological samples (n = 3). Data are means ± SE (n = 3 biological replicates) and Y-axis represents relative quantity (RQ). RQ was calculated using the equation; RQ = 2−ΔΔCt. Different letters (a–d) indicate statistically significant differences between treatments (Duncan’s multiple range test P < 0.05).
The increment in the expression of SalSyn, SalR, and SalAT mediated by consortium inoculation was higher than the individual endophyte inoculation (Figure 6). Consortium-inoculated plants had 7.0-, 5.9-, and 6.4-fold higher SalSyn, SalR, and SalAT expression compared to non-inoculated endophyte-free control plants. Individual inoculation of SM1B and SM3B could increase the expression of SalSyn (2.8–5.8-fold), SalR (2.7–2.8-fold) and SalAT (3.6-fold) to lower levels only. The consortium-inoculated plants and SM1B inoculated plants had reduced expression of T6ODM, however, SM3B inoculated plants had 2.4-fold higher T6ODM expression compare to endophyte-free control plants (Figure 6). Expression of CODM was not affected significantly due to the consortium- and SM1B-inoculation. However, SM3B inoculated plants had almost twofold higher CODM expression than the endophyte-free control plants, SM1B and consortium inoculated plants (Figure 6). Expression of COR, considered to be the most critical regulatory gene, was 2.3-fold higher in the consortium- inoculated plants whereas SM1B alone could increase it to the level of 0.9-folds only and SM3B inoculation reduced COR expression compare to endophyte-free control (Figure 6). Upregulation of multiple genes involved in morphine biosynthesis due to inoculation of a consortium of endophytes is represented in Figure 7 (Supplementary Table S3).
Figure 6. Effect of endophytes inoculation on the expression of genes involved in morphine biosynthesis. Total RNA was isolated from the green capsule of 90 days old P. somniferum plants inoculated with endophytes SM1B and SM3B individually, and in the form of consortium (SM1B+SM3B), reverse transcribed and used as a template for RT-qPCR with SYBR Green detection. The capsules of non-inoculated endophyte free plants were used as a control. Expression of SalSyn (A), SalR (B), SalAT (C), T6ODM (D), CODM (E), COR (F) was analyzed. Results were normalized to actin (reference transcript) and are shown relative to the level in non-inoculated endophyte free control plants (calibrator). qRT-PCR was performed on triplicate technical replicates of triplicate biological samples (n = 3). Data are means ± SE (n = 3 biological replicates) and Y-axis represents relative quantity (RQ). RQ was calculated using the equation; RQ = 2−ΔΔCt. Different letters (a–d) indicate statistically significant differences between treatments (Duncan’s multiple range test P < 0.05).
Figure 7. Upregulation of multiple genes involved in morphine biosynthesis due to inoculation of a consortium of endophytes. SM1B (Acinetobacter) endophytes upregulated the expression of most of the genes of BIAs pathway except T6ODM and CODM whereas SM3B (Marmoricola sp.) inoculation upregulated the T6ODM and CODM. Consortium-inoculation (SM1B+SM3B) could enhance the morphine production by upregulating, the expression of most of the BIA biosynthetic genes (TYDC, NCS, 6OMT, CNMT, NMCH, SalSyn, SalR, SalAT, COR) involved in morphine biosynthesis to a greater extent, than that of single-inoculations. Expression of different genes was presented in square boxes. Red color of boxes indicate upregulated expression (intensity of red color shows the level of expression, i.e., more red more expression and vice versa), yellow color shows the level of expression in non-inoculated endophyte free control plants and light yellow color shows downregulated expression. Enzyme abbreviations: TyrAT, tyrosine aminotransferase; TYDC, tyrosine/ dopa decarboxylase; NCS, norcoclaurine synthase; 6OMT, (S)-norcoclaurine 6-O-methyltransferase; CNMT, (S)- coclaurine N-methyltransferase; NMCH, (S)-N-methylcoclaurine 3-hydroxylase; 4-OMT, (S)-30-hydroxy-N-methylcoclaurine 4-O-methyltransferase; DRS, 1,2-dehydroreticuline synthase; SalSyn, salutaridine synthase; SalR, salutaridine reductase; SalAT, salutaridinol 7-O-acetyltransferase; T6ODM, thebaine 6-O-demethylase; COR, codeinone reductase; CODM, codeine O-demethylase.
Discussion
Due to limitations related to social acceptability and cost-effectiveness in genetic manipulation approaches such as the development of transgenics by overexpressing single or multiple genes simultaneously, there is a need of an alternate sustainable approach for improvement of crop production and a healthier environment. Microbes associated with plants such as plant growth promoting rhizobacteria (PGPRs) and endophytes have been used as promising candidates for sustainable agriculture. They are found to be involved in improving plant growth and protection from environmental stresses. We have demonstrated that PGPRs and endophytes can protect the plant from environmental stresses by modulating the phytohormone status of plants and by modulating the expression of stress-responsive genes (Bharti et al., 2016; Barnawal et al., 2016, 2017). Also, endophytes residing in most of the parts of the plant including root, stem, leaves, seeds, fruits, flower, and found to perform the tissue-specific function may also modulate primary and secondary metabolism of a plant (Pandey et al., 2016a). Endophytes are considered better candidates than PGPR because they function residing inside the plant and therefore, they do not face any competition with the resident soil microorganisms and also escape from the surrounding harsh environmental conditions. Diverse endophytes are available with specific characteristics related to growth promotion ability and stress tolerance such as an ability of nutrient acquisition, nitrogen fixation, phosphate solubilization, ACC-deaminase production, phytohormone production and suppression of plant pathogens that make them useful for improvement of plant growth and for providing tolerance to environmental stresses. A single endophyte may not be able to carry said multiple growth promotion activities. Therefore, to obtain maximum growth and crop yield consortium of endophytes having different characteristics will be a better approach than single-endophyte inoculation.
Previously, we have demonstrated that inoculation of an endophyte of opium poppy SM1B (Acinetobacter sp.) which was isolated from the capsule of an alkaloid-rich cultivar of P. somniferum (cv. Sampada) enhanced morphine content substantially in alkaloid-less cv. Sujata. However, the accumulation of intermediates of morphine synthesis; thebaine and codeine was very low in SM1B-inoculated plants (Pandey et al., 2016b). Also, another endophyte SM3B (Marmoricola sp.) inoculation could not increase the morphine content significantly, but substantially increased the thebaine (another important alkaloid) accumulation compared to endophyte-free control plants. Gene expression study of these endophytes (SM1B and SM3B) inoculated plants showed that most of the genes (except T6ODM and CODM) of BIA biosynthesis were upregulated in SM1B-inoculated plants. On the other hand, T6ODM and CODM could be upregulated by SM3B-inoculation. The present study therefore aimed at upregulation of all the genes involved in morphine biosynthesis by inoculating these two endophytes through complementarity. Single-inoculation of these endophytes enhanced the morphine content in a morphine-less cultivar of P. somniferum (cv. Sujata), however, this could not increase up to the level of an alkaloid-rich cultivar of P. somniferum (cv. Sampada) from where these endophytes were isolated. Therefore, the efficacy of combined inoculation of both the endophytes (SM1B and SM3B) in maximizing the production of morphine, was tested in P. somniferum cv. Sujata which is otherwise morphine-less variety.
Both the endophytes (SM1B and SM3B) were able to survive and perpetuate in combination and could colonize successfully in the capsule of the poppy plant when applied as a consortium. However, consortium inoculation slightly reduced the colonization of SM1B compared to individual SM1B inoculation. Consortium inoculation could enhance the growth of P. somniferum cv. Sujata plants compared to individual inoculations by increasing the content of photosynthetic pigments, photosynthetic rate, transpiration rate and stomatal conductance. The increased photosynthetic efficiency of the consortium-inoculated poppy plant could enhance the biomass (73.3%), capsule weight (66%), and seed yield (79%) which was substantially higher than single-inoculation.
Combined-inoculation of SM1B and SM3B could generate higher increments in morphine content compared to individual-inoculations by modulating the utilization of the key intermediates (reticuline and thebaine) of morphine biosynthesis. Enhanced morphine content was due to more thebaine production in consortium inoculated plants (0.07% content gDW−1) compared to individual-inoculation (0.005% content gDW−1 in case of SM1B and 0.035% content gDW−1 in case of SM3B). Higher thebaine production could result in lower accumulation of reticuline (upstream intermediate) in consortium inoculated plants compare to single inoculation of SM3B. However, very low accumulation of reticuline and thebaine indicates complete utilization of these intermediates for morphine biosynthesis in the case of SM1B-inoculated plants.
Consortium-inoculation could enhance the morphine production by upregulating, the expression of most of the BIA biosynthetic genes (TYDC, NCS, 6OMT, CNMT, NMCH, SalSyn, SalR, SalAT, COR) involved in morphine biosynthesis to a greater extent, than that of single-inoculations (Figure 7). Expression of TYDC, NCS, 6OMT, CNMT, NMCH which are the genes upstream to reticuline biosynthesis was higher in consortium inoculated plants, clearly indicating greater flux of BIA biosynthesis toward central intermediate, i.e., reticuline. However, due to the increased expression of SalSyn, SalR, and SalAT (which was higher in consortium-inoculated plants than that of single-inoculation); reticuline accumulation was found to be low than that of SM3B-inoculated plants. Due to higher expression of SalSyn, SalR, and SalAT than that of single inoculation; conversion of reticuline to thebaine was higher, and that could have resulted in more thebaine accumulation in consortium-inoculated plants compared to single-inoculated plants.
Here, we probably for the first time demonstrated that endophytes affecting the host metabolism at gene expression level having different modes/targets of action could be combined to compliment the inability of one endophyte to upregulate certain genes thus improving the biosynthesis of secondary metabolites. Also, by using the consortium of different endophytes with diverse functional growth attributes, we can compliment the functional deficiency of one by another. These endophytes would be a better alternative in place of transgenic plants and also could be explored for the development of designer plants. Here, in the present study, we also established that consortium of selected endophytes could be applied to overexpress multiple genes to increase the plant yield/desired product. The present study was conducted in morphine-less genotype Sujata to establish the role of combined inoculation of endophytes in enhancing morphine, and therefore the magnitude of enhancement was quite high considering very low concentrations inherently present in genotype Sujata. These enhancements, however, might not be possible in morphine/alkaloid rich variety(ies). It is, however, essential to understand the detailed mechanism of action of promising endophytes, and then a combined approach can be applied to attain maximum yields.
Author Contributions
AK, TR, and SP conceived and designed the experiments and wrote the manuscript. TR and MS performed the experiments. AK, KS, TR, SP, and AP analyzed the data. All authors read and approved the manuscript.
Funding
This work was supported by grant NWP BSC0117 (XII Five Year Plan Network Project) from the Council of Scientific and Industrial Research (CSIR), India.
Conflict of Interest Statement
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.
Acknowledgments
The authors are thankful to Director, CSIR-CIMAP, Lucknow, India for providing laboratory facilities and Director, CSIR-NBRI, Lucknow, India for providing BIA standards. SP greatly acknowledge CSIR for financial assistance in the form of Senior Research Associateship. TR is thankful to Department of Science and Technology, India for INSPIRE fellowship. AP is grateful to DST for INSPIRE Faculty award (IFA13-LSPA-20).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.00925/full#supplementary-material
References
Barnawal, D., Bharti, N., Tripathi, A., Pandey, S. S., Chanotiya, C. S., and Kalra, A. (2016). ACC-deaminase-producing endophyte Brachybacterium paraconglomeratum strain SMR20 ameliorates Chlorophytum salinity stress via altering phytohormone generation. J. Plant Growth Regul. 35, 553–564.
Barnawal, D., Pandey, S. S., Bharti, N., Pandey, A., Ray, T., Singh, S., et al. (2017). ACC deaminase-containing plant growth-promoting rhizobacteria protect Papaver somniferum from downy mildew. J. Appl. Microbiol. 122, 1286–1298. doi: 10.1111/jam.13417
Bharti, N., Pandey, S. S., Barnawal, D., Patel, V. K., and Kalra, A. (2016). Plant growth promoting rhizobacteria Dietzia natronolimnaea modulates the expression of stress responsive genes providing protection of wheat from salinity stress. Sci. Rep. 6:34768. doi: 10.1038/srep34768
Bhunia, A. K., Johnson, M. C., and Ray, B. (1988). Purification, characterization and antimicrobial spectrum of a bacteriocin produced by Pediococcus acidilactici. J. Appl. Bacteriol. 65, 261–268.
Bohnert, H. J., Gong, Q., Li, P., and Ma, S. (2006). Unraveling abiotic stress tolerance mechanisms–getting genomics going. Curr. Opin. Plant Biol. 9, 180–188. doi: 10.1016/j.pbi.2006.01.003
Del Giudice, L., Massardo, D. R., Pontieri, P., Bertea, C. M., Mombello, D., Carata, E., et al. (2008). The microbial community of Vetiver root and its involvement into essential oil biogenesis. Environ. Microbiol. 10, 2824–2841. doi: 10.1111/j.1462-2920.2008.01703.x
Di Fiore, S., and Del Gallo, M. (1995). “Endophytic bacteria: their possible role in the host plant,” in Azospirillum VI and Related Microorganisms, eds I. Fendrik, M. del Gallo, J. Vanderleyden, and M. de Zamaroczy (Heidelberg: Springer), 169–187.
Egamberdieva, D., Wirth, S. J., Shurigin, V. V., Hashem, A., and Abd Allah, E. F. (2017). Endophytic bacteria improve plant growth, symbiotic performance of chickpea (Cicer arietinum L.) and induce suppression of root rot caused by Fusarium solani under salt stress. Front. Microbiol. 8:1887. doi: 10.3389/fmicb.2017.01887
Gundel, P. E., Pérez, L. I., Helander, M., and Saikkonen, K. (2013). Symbiotically modified organisms: nontoxic fungal endophytes in grasses. Trends Plant Sci. 18, 420–427. doi: 10.1016/j.tplants.2013.03.003
Gundel, P. E., Seal, C., Biganzoli, F., Molina Montenegro, M., Vazquez-de-Aldana, B., Zabalgogeazcoa, I., et al. (2018). Occurrence of alkaloids in grass seeds symbiotic with vertically-transmitted epichloë fungal endophytes and its relationship with antioxidants. Front. Ecol. Evol. 6:211. doi: 10.3389/fevo.2018.00211
Halpin, C. (2005). Gene stacking in transgenic plants–the challenge for 21st century plant biotechnology. Plant Biotechnol. J. 3, 141–155. doi: 10.1111/j.1467-7652.2004.00113.x
Hardoim, P. R., Van Overbeek, L. S., Berg, G., Pirttilä, A. M., Compant, S., Campisano, A., et al. (2015). The hidden world within plants: ecological and evolutionary considerations for defining functioning of microbial endophytes. Microbiol. Mol. Biol. Rev. 79, 293–320. doi: 10.1128/MMBR.00050-14
Kuklinsky-Sobral, J., Araújo, W. L., Mendes, R., Geraldi, I. O., Pizzirani-Kleiner, A. A., and Azevedo, J. L. (2004). Isolation and characterization of soybean-associated bacteria and their potential for plant growth promotion. Environ. Microbiol. 6, 1244–1251.
Kumar, A., Patil, D., Rajamohanan, P. R., and Ahmad, A. (2013). Isolation, purification and characterization of vinblastine and vincristine from endophytic fungus Fusarium oxysporum isolated from Catharanthus roseus. PLoS One 8:e71805. doi: 10.1371/journal.pone.0071805
Kusari, S., Verma, V. C., Lamshoeft, M., and Spiteller, M. (2012). An endophytic fungus from Azadirachta indica A. Juss. that produces azadirachtin. World J. Microbiol. Biotechnol. 28, 1287–1294. doi: 10.1007/s11274-011-0876-2
Li, J., Zhao, G. Z., Varma, A., Qin, S., Xiong, Z., Huang, H. Y., et al. (2012). An endophytic pseudonocardia species induces the production of artemisinin in Artemisia annua. PLoS One 7:e51410. doi: 10.1371/journal.pone.0051410
Lichtenthaler, H. K., and Buschmann, C. (2001). “Chlorophylls and carotenoids: Measurement and characterization by UV-VIS spectroscopy,” in Current Protocols in Food Analytical Chemistry (CPFA), eds R. E. Wrolstad, T. E. Acree, H. An, E. A. Decker, M. H. Penner, D. S. Reid, et al. (New York, NY: Wiley).
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2_ΔΔCt method. Methods 25, 402–408.
Maggini, V., De Leo, M., Mengoni, A., Gallo, E. R., Miceli, E., Reidel, R. V. B., et al. (2017). Plant-endophytes interaction influences the secondary metabolism in Echinacea purpurea (L.) Moench: an in vitro model. Sci. Rep. 7:16924. doi: 10.1038/s41598-017-17110-w
Otieno, N., Lally, R. D., Kiwanuka, S., Lloyd, A., Ryan, D., Germaine, K. J., et al. (2015). Plant growth promotion induced by phosphate solubilizing endophytic Pseudomonas isolates. Front. Microbiol. 6:745. doi: 10.3389/fmicb.2015.00745
Pandey, S. S., Singh, S., Babu, C. V., Shanker, K., Srivastava, N. K., and Kalra, A. (2016a). Endophytes of opium poppy differentially modulate host plant productivity and genes for the biosynthetic pathway of benzylisoquinoline alkaloids. Planta 243, 1097–1114. doi: 10.1007/s00425-016-2467-9
Pandey, S. S., Singh, S., Babu, C. V., Shanker, K., Srivastava, N. K., Shukla, A. K., et al. (2016b). Fungal endophytes of Catharanthus roseus enhance vindoline content by modulating structural and regulatory genes related to terpenoid indole alkaloid biosynthesis. Sci. Rep. 6:26583. doi: 10.1038/srep26583
Pandey, S. S., Singh, S., Pandey, H., Srivastava, M., Ray, T., Soni, S., et al. (2018). Endophytes of Withania somnifera modulate in planta content and the site of withanolide biosynthesis. Sci. Rep. 8:5450. doi: 10.1038/s41598-018-23716-5
Quecine, M. C., Araújo, W. L., Rossetto, P. B., Ferreira, A., Tsui, S., Lacava, P. T., et al. (2012). Sugarcane growth promotion by the endophytic bacterium Pantoea agglomerans 33.1. Appl. Environ. Microbiol. 78, 7511–7518. doi: 10.1128/AEM.00836-12
Rodriguez, R. J., Henson, J., Van Volkenburgh, E., Hoy, M., Wright, L., Beckwith, F., et al. (2008). Stress tolerance in plants via habitat-adapted symbiosis. ISME J. 2, 404–416. doi: 10.1038/ismej.2007.106
Rosenblum, A., Marsch, L. A., Joseph, H., and Portenoy, R. K. (2008). Opioids and the treatment of chronic pain: controversies, current status, and future directions. Exp. Clin. Psychopharmacol. 16, 405–416. doi: 10.1037/a0013628
Saikkonen, K., Faeth, S. H., Helander, M., and Sullivan, T. J. (1998). Fungal endophytes: a continuum of interactions with host plants. Annu. Rev. Ecol. Syst. 29, 319–343.
Scherling, C., Ulrich, K., Ewald, D., and Weckwerth, W. (2009). A metabolic signature of the beneficial interaction of the endophyte Paenibacillus sp. Isolate and in vitro–grown poplar plants revealed by metabolomics. Mol. Plant Microbe Interact. 22, 1032–1037. doi: 10.1094/MPMI-22-8-1032
Schmittgen, T. D., and Livak, K. J. (2008). Analyzing real-time PCR data by the comparative CT method. Nat. Protoc. 3, 1101–1108. doi: 10.1038/nprot.2008.73
Sharma, J. R., Lal, R. K., Gupta, A. P., Misra, H. O., Pant, V., Singh, N. K., et al. (1999). Development of non-narcotic (opiumless and alkaloid-free) opium poppy, Papaver somniferum. Plant Breed. 118, 449–452.
Sharma, J. R., Lal, R. K., Gupta, M. M., Verma, R. K., and Misra, H. O. (2002). A novel non-narcotic seed variety sujata of opium poppy (Papaver somniferum). J. Med. Arom. Plant Sci. 24, 481–485.
Shweta, S., Zuehlke, S., Ramesha, B. T., Priti, V., Mohana Kumar, P., Ravikanth, G., et al. (2010). Endophytic fungal strains of Fusarium solani, from Apodytes dimidiata E. Mey. ex Arn (Icacinaceae) produce camptothecin, 10- hydroxycamptothecin and 9-methoxycamptothecin. Phytochemistry 71,117–122.
Singh, S. K., Strobel, G. A., Knighton, B., Geary, B., Sears, J., and Ezra, D. (2011). An endophytic phomopsis sp. Possessing bioactivity and fuel potential with its volatile organic compounds. Microb. Ecol. 61, 729–739. doi: 10.1007/s00248-011-9818-7
Soliman, S. S., Tsao, R., and Raizada, M. N. (2011). Chemical inhibitors suggest endophytic fungal paclitaxel is derived from both mevalonate and non-mevalonate-like pathways. J. Nat. Prod. 74, 2497–2504. doi: 10.1021/np200303v
Staniek, A., Woerdenbag, H. J., and Kayser, O. (2008). Endophytes: exploiting biodiversity for the improvement of natural product-based drug discovery. J. Plant Interact. 3, 75–93.
Strobel, G., Daisy, B., Castillo, U., and Harper, J. (2004). Natural products from endophytic microorganisms. J. Nat. Prod. 67, 257–268.
Tajehmiri, A., Darsanaki, R. K., Moslem, M. N., Lozoumi, Z., Kolavani, M. H., and Aliabadi, M. A. (2014). antimicrobial activity of lactobacillus spp. Isolated from commercial yoghurts against pathogenic bacteria. J. Pure Appl. Microbiol. 8, 2211–2215.
Tiwari, R., Awasthi, A., Mall, M., Shukla, A. K., Srinivas, K. S., Syamasundar, K. V., et al. (2013). Bacterial endophyte-mediated enhancement of in planta content of key terpenoid indole alkaloids and growth parameters of Catharanthus roseus. Ind. Crops Prod. 43, 306–310.
Waller, F., Achatz, B., Baltruschat, H., Fodor, J., Becker, K., Fischer, M., et al. (2005). The endophytic fungus Piriformospora indica reprograms barley to salt-stress tolerance, disease resistance, and higher yield. Proc. Natl. Acad. Sci. U.S.A. 102, 13386–13391.
Keywords: endophytes, consortium, benzylisoquinoline alkaloids, morphine, Papaver somniferum, opium
Citation: Ray T, Pandey SS, Pandey A, Srivastava M, Shanker K and Kalra A (2019) Endophytic Consortium With Diverse Gene-Regulating Capabilities of Benzylisoquinoline Alkaloids Biosynthetic Pathway Can Enhance Endogenous Morphine Biosynthesis in Papaver somniferum. Front. Microbiol. 10:925. doi: 10.3389/fmicb.2019.00925
Received: 08 January 2019; Accepted: 11 April 2019;
Published: 30 April 2019.
Edited by:
Sheridan Lois Woo, University of Naples Federico II, ItalyReviewed by:
Roberta Fulthorpe, University of Toronto Scarborough, CanadaPiippa Riitta Wäli, University of Oulu, Finland
Copyright © 2019 Ray, Pandey, Pandey, Srivastava, Shanker and Kalra. 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: Alok Kalra, YWxvay5rYWxyYUB5YWhvby5jb20=
†These authors have contributed equally to this work