Introduction

The world-wide soybean (Glycine max, L.) production has increased several fold starting from 1960s’ and currently soybean is one of the largest source of vegetable oil and animal feed (Pagano & Miransari, 2016). The demand for soybean will increase in future due to the rising population and the increased demand in both animal feed and products for human consumption. Growing of soybean meets a number of challenges such as biotic stress (micropathogens, pests and weeds). It is estimated that only in the US soybean diseases decrease the yield by up to 11% (Hartman et al., 2016).

Currently, there are no individual soybean genes related to the resistance to single micropathogenes, although the search for the quantitative trait loci associated with the resistance continues (Kandel et al., 2018; Zatybekov, Abugalieva, Didorenko, Rsaliyev, & Turuspekov, 2018). A large number of genes are associated with abiotic and biotic stress resistance in plants (Cabane, Afif, & Hawkins, 2012; Pérez-Clemente et al., 2013; Tran & Mochida, 2010). They are related to the various metabolic and signaling pathways including reactive oxygen species signaling (Baxter, Mittler, & Suzuki, 2014), plant hormones, inorganic ion fluxes (Kissoudis, Van de Wiel, Visser, & Van der Linden, 2014) as well as transcription factors (Birkenbihl, Liu, & Somssich, 2017). To date, genetic manipulation of a single or a few genes involved in signaling/regulatory pathways are used to improve environmental stress resistance of soybean plants (Homrich, Strohm, Weber, & Zanettini, 2012).

Lignin is a phenolic compound providing strength to the secondary cell wall in vascular plants, that is assigned to a broad range of physiological processes participating in plant growth, providing the rigidity to the cell walls, the natural mechanical barrier and defense against pathogen penetration. It acts as a barrier against pathogen penetration as well as a toxins produced by the pathogens, promoting plant immunity (Lee et al., 2019; Sattler & Harris, 2013). Lignin is produced by phenylpropanoid pathway that provides monolignol compounds such as p-coumaryl alcohol, coniferyl alcohol, and sinapyl alcohol which are precursors of p-hydroxyphenyl (H), guaiacyl (G), syringyl (S) monomeric units of lignin, respectively. Lignin composition was shown to change as a response to abiotic and biotic stresses (Liu, Luo, & Zheng, 2018; Moura, Bonine, Viana, Dornelas, & Mazzafera, 2010). The molecular mechanisms of lignin biosynthesis are well studied in model organisms such as arabidopsis (Arabidopsis thaliana) and poplar (Populus spp.) (Bonawitz & Chapple, 2010; Vanholme et al., 2012) as well as in some agricultural plants: rice (Oryza sativa) (Gui, Shen, & Li, 2011), grapevine (Vitis vinifera) (Cavallini et al., 2015) and others. The first step of lignin synthesis is deamination of phenylalanine by PAL enzymes to produce cinnamic acid – a step that is common between all metabolites of the phenylpropanoid pathway. Cinnamic acid is further converted in p-coumaric acid by C4H (cinnamate 4-hydroxylase) enzyme. Intermediate products are processed by other enzymes such as 4CL (4-coumarate-CoA ligase), COMT (caffeic acid methyltransferase), F5H (ferulate 5-hydroxylase) to synthesize aldehyde precursors of coniferyl and synapil alcohols. These are further reduced to the respective alcohols by CAD enzymes (cinnamyl alcohol dehydrogenases) (reviewed in Barros, Serk, Granlund, & Pesquet, 2015). Lignin synthesis is regulated by a number of transcription factors such as R2R3-MYB family (Bonawitz & Chapple, 2010). In Arabidopsis thaliana MYB46 and MYB83 were shown to regulate the key lignin synthesis genes (Kim, Kim, Ko, Kang, & Han, 2014). Similar functions are observed for Pinus taeda MYB46 homologs (Patzlaff et al., 2003).

Additionally, phenylpropanoid pathway is known to produce a large number of secondary metabolites that include flavanols, isoflavanoids, anthocyanins (Vogt, 2010). Some of these metabolites are demonstrated to be components of plant innate immunity, preventing the spread of the pathogens (Piasecka, Rey, & Bednarek, 2015). Genetic manipulation of the phenylpropanoid pathway may not only alter lignin synthesis but also enhance plant innate immunity. Therefore one of the strategies to increase the resistance of soybean plants is to increase the expression of genes involved in lignin biosynthesis (Lozovaya et al., 2007; Lozovaya, Lygin, Zernova, & Widholm, 2005; Lygin et al., 2013, 2010; Lygin, Rahman, Ulanov, Widholm, & Lozovaya, 2012). This approach aims to produce biotech products of second generation (Pérez-Clemente et al., 2013). These plants be more resistant to stress and could provide an increased yield, in comparison with biotech crops of first generation, such as Roundup Ready crops.

The most common techniques for soybean genetic transformation are particle bombardment (McCabe, Swain, Martinell, & Christou, 1988) and transformation mediated by Agrobacterium tumefaciens (Yin & Zhang, 2010) . Transformation of mature cotyledonary node or cotyledonary explants derived from mature soybean seed is currently regarded as the most effective methods (Li et al., 2017; Paz, Martinez, Kalvig, Fonger, & Wang, 2006). However, these methods have some critical disadvantages: they are genotype specific and include an expensive tissue culture step and plant regeneration in vitro, which takes about one year before transgenic first generation plants (T₁) are recovered and in most cases zero generation (T₀) plants are weak and require careful handling.

Several approaches are used to omit the complicated tissue culture step. One of them is floral dip method primarily used for A. thaliana transformation (Clough & Bent, 1998). In this method developing floral tissues are dipped into a solution containing A. tumefaciens, as well as surfactant and sucrose. However, this method is not suitable for legumes due to the flower structure. An alternative method is pollen-tube pathway transformation technique which is based on the removal of the recipient plant’s stigma shortly after pollination and introduction of an exogenous DNA solution into the severed style of the recipient plant (Ali, Bang, Chung, & Staub, 2015). Pollen-tube pathway soybean transformation was described by a number of authors however, it seems to be not efficient and irreducible (Liu, Su, An, & Yang, 2009; Li & Wu, 2007; Li, Nelson, Widholm, & Bent, 2002; Shou, Palmer, & Wang, 2002).

In this paper we propose a novel technique of germ-line genetic transformation via natural pollen tubes using A. tumefaciens mediated DNA transfer. We aim to use this method to target lignin content and composition in soybean plants which can significantly improve basic resistance to soybean pathogens. The ultimate goal of this work is to produce new strains of soybean for breeders and biodiversity with improved biotic, abiotic stress resistance and yield enhancement.

Material and methods

Plant material

Three cultivars of USA breeding: Jack and Jack x 4 (resistant to biotic stress), and Spencer (sensitive to micropathogens infection), 7 varieties of Kazakhstan breeding (Kazakhstan, Eureka, Vita, Tazhan, Zara, Perizat, Dannaya) were used for introduction of lignification genes. Ten biological replicates of each homozygous transgene and wild-type were grown simultaneously in a random block design, in the same environment in greenhouse. Plants were grown under long-day conditions (16 hours light, 28°C, and 65% humidity) to allow the development of health soybean plants.

Molecular vectors

Germ-line transformation was performed using A. tumefaciens, strain EHA 105 transformed with following vectors (Table 1).

Transformation in A. tumefaciens or E. coli (for glycerol plasmid stocks and plasmid propagation) was performed using electroporation.

Preparation of Agrobacterium tumefaciens suspension for flower stigma pipetting

Single colony of transformed A. tumefaciens was transferred from Petri dishes to 3 mL of LB medium with 100 mg L^-1 kanamycin (Sigma Aldrich), and incubated for 1 – 1.5 days at 28°C on a shaker at 150 rpm. Cells were harvested by centrifugation at 4500 g for 10 minutes and resuspended in LB medium without antibiotics. Pipetting of flower stigma was performed with the bacterial suspension 1010 cells mL^-1 (equal OD=1 when measured at ƛ=600 nm using Beckman Coulter DU 800 spectrophotometer (USA)). Surfactant Pluronic F68 (0.001%) (Thermo Fisher) was added to bacterial suspension to activate vir-zone of T-plasmid and enhance transformation efficiency.

Agrobacterium tumefaciens suspension pipetting in soybean flower stigma

Soybean flower petals were separated and flower was cut in the carpel remove 1/3 of stigma with small scissors. About 2 µL of A. tumefaciens suspension transformed with a vector containing a gene of interest was placed in the center of the open flower to the stigma by a microsyringe or a micropippete (Figure 1a). After 30 minutes flowers were inspected and if the suspension was not enough the second pipetting was performed. Treated flowers were marked by color threads correspondent to the genes of interest (Figure 1b). Pipetting time was optimized to introduce transformed A. tumefaciens at 3-4 stages of soybean flower development, when the anthers were fully formed and pollen tubes are growing into mature, but not yet dividing zygote (Figure 1c). To synchronize pipetting time all soybean flowers opened the day before were removed at 8:00. It was determined that the stigma allocation of substances that stimulate the germination of pollen starts at 9:00, at around 9:15 pollen is attached to the stigma and the pollen germination begins and continues till 9:30 (Figure 1d-f). Therefore, the best time for pipetting in a greenhouse conditions (at 28°C) was 9:00 am – 12:00 pm. In the field the time range can be 11:00 – 13:00 the same day, when the flowers have opened and the resulting pollen tubes are able to pass the desired genetic material (at the time of the pollen activity maximum). The timing of pipetting in the field may be dependent on temperature and solar radiation.

Molecular analysis

Plasmid DNA was extracted from an overnight culture by using QIAprep Spin Mini prep Kit QIAGEN (Qiagen, Valencia, CA, USA). Genomic DNA was extracted from fresh mature but young transgenic soybean leaves using the method (Weigel & Glazebrook, 2009)⁠ and was used in PCR. Following primers were used (Table 2).

Transgene presence in plant leaves tissues was confirmed by the PCR reaction using Taq DNA polymerase (New England Biolabs, Ipswich, MA, USA) with denaturation at 94 °C (1 min.), annealing at 48°C or 54°C (40 s) for PtMYB4., PAL5, F5H, CAD1 constructs and extension at 72°C (90 s), for 35 cycles.

Multiplex PCR was performed similarly, with the exception of the temperature of the annealing stage (60°C).

RNeasy Plant Mini Kit (Qiagen, Valencia, CA, USA) was used for extraction of total RNA from examples. DNA was digested with DNase I, using the TURBO DNA-free kit (Ambion/Life technologies, Carlsbad, CA, USA). cDNA was synthesized using oligo (dT) primer and reverse transcriptase enzyme, Promega GoScript kit, (Madison, WI, USA) for positive control or excluding reverse transcriptase enzyme as a negative control for genomic DNA contamination. To amplify PtMYB4, PAL5, F5H and CAD1 transcripts from cDNA same primers as for genomic DNA detection were used.

Southern and northern blot

Total genomic DNA of transgenic plants was digested with BamHI endonuclease (Fermentas, ThermoFisher, USA). Fragment size was 1,5 kB for PAL5 transformants and 1,2 kb for CAD1 transformants. DNA was precipitated in isopropanol and sodium acetate solution, precipitate was washed with ethanol and dissolved in water. After DNA electrophoresis DNA was transferred to a nylon membrane (Amersham, UK) using denaturation alkaline buffer (1 N NaOH + 0.4 M NaOH + 0.6 M NaCl). Membrane was dried for 30 min., DNA was fixed in UV-light (265 nm) for 5 minutes. Membrane was used for subsequent hybridization with P³²-labeled probes prepared from PAL5 and GUS PCR products (as described in molecular analysis section) using NEBlot® Kit (New England Biolabs, Ipswich, MA, USA). After a prehybridization at 50°C hybridization was performed for 16h at 50°C in hybridization oven (Amersham, UK). Membrane was washed with 2 times in 1X SSC and 0.1% SDS at 50°С for 15 minutes and 3 times in 0.2X SSC and 0.1% SDS at 58°C. The intensity of radioactive label was measured using radioactivity counter. Exposure of the membrane with the X-ray film (Hyperfilm TMMP, Amersham, UK) was carried out with an amplifying screen. Same membrane was hybridized with the gus-probe.

For Northern blot total RNA was extracted from the same plants using guanidinium thiocyanate method. 5 ug RNA was subjected to agarose gel electrophoresis and was transferred to a nylon membrane (Amersham, UK). Membrane was prepared similarly to the membrane used for Southern blotting. Membrane was used for subsequent hybridization with P³²-labeled probes. Probes were prepared from PCR products (as described in the molecular analysis section) using NEBlot® Kit (New England Biolabs, Ipswich, MA, USA). Membrane was washed with 0.2X SSC and 1% SDS at 58°C for 1 hour.

Lignin analyses

Lignin content and composition was performed as described in (Van Acker et al., 2013). The lignin composition was investigated with thioacidolysis (Robinson & Mansfield, 2009). The monomers involved in β–O–4-ether bonds, released upon thioacidolysis, were detected with gas chromatography (GC) as their trimethylsilyl (TMS) ether derivatives on a Hewlett-Packard HP 6890 Series GC system (Agilent, Santa Clara, CA, USA) coupled with a HP-5973 mass-selective detector. The GC conditions were as described in the protocol. The quantitative evaluation was carried out based on the specific prominent ions for each compound. Response factors for H (hydroxyphenyl), G (guaiacyl), and S (syringyl) units were taken from. The sum of H, G, and S is a good estimate of the total thioacidolysis yield and, thus, the condensation degree of the lignin polymer.

Enzyme assays

PAL activity was determined by the conversion of L-phenylalanine to trans-cinnamic acid essentually as described previously (Edwards & Kessman, 1992)⁠ with a few modifications: 300 mg of fresh leaves were ground in liquid nitrogen and extracted with 1200 µL of the extraction buffer, containing 0.05 M TrisHCl (pH-8.8), 0.5% ascorbate, 10% glycerol and 10 mM 13-mercaptoethanol. 100 µL of crude enzyme extract was combined with 400 µL of 0.05 M Tris buffer (pH 8.8), containing 0.2 mM phenylalanine as substrate. The reaction mixture was incubated for 30, 60, 120 and 180 min. at 37°C. The reaction was stopped by adding of 100 µL 0.5M HCl. Cinnamic acid was extracted with 1 mL of toluene and the OD was measured at 290 nm, with toluene used as a blank. The calibration curve was obtained using trans-cinnamic acid from Sigma. Total protein was measured using Bio-Rad Protein Assay (Bradford method).

For CAD1 enzyme activity determination plant tissue (300 mg) was ground in liquid nitrogen and crude extract was prepared using the extraction buffer containing Tris HCL 0.05M, pH 8.8, 10% glycerol and 10 mM β-mercaptoethanol. The CAD assay was performed according to (Mansell, Gross, Stöckigt, Franke, & Zenk, 1974)⁠ with a few modifications. The enzyme activity was determined by the conversion of alcohol to aldehyde (reverse reaction). The formation of hydroxycinnamaldehyde at 37°C was monitored at 400 nm using the following molar extinction coefficients: 1.78 x 10⁴ M^-1 cm^-1 for coniferyl aldehyde and 14.7 x 103 M^-1 cm^-1 for sinapyl aldehyde (Hawkins & Boudet, 1994)⁠. The reaction mixture contained in a total volume of 500 µL : 25 mM Tris- HCl, pH 8.8, 200 µl NADP (2mM), 200 µL (2mM) coniferyl alcohol/sinapyl alcohol and 100 µL crude extract. The same mixture, except 100 µL of TrisHCl buffer instead of the crude extract, was used as a blank. Protein concentrations were determined by the Bradford dye-binding assay.

Results

Genetic transformation of soybean plants by Agrobacterium tumefaciens pipetting for introduction of genes of lignification

We performed A. tumefaciens-mediated transformation of 1070 flowers of 3 USA and 7 established and zoned in Kazakhstan soybean varieties. Overall we obtained 833 pods (78% efficiency) containing 2000 putative transgenic soybean seeds of zero generation (T₀), expressing the genes of interest (PtMYB4, PAL5, F5H, and CAD1), respectively. Seeds forming efficiency on average was 187%. We did not observe any significant differences by size, shape or mass of pods and seeds between T₀ and control plants (data is not shown). The average number of seeds in a pod at T₀ was 2.5 seeds. T₀ seeds were planted to obtain plants of the first generation (T₁) for the further molecular analysis.

Molecular analysis of T1 transgenic soybean plants

We analyzed the presence of the transgenes in the DNA and the presence of the transgene transcripts in T₁ soybean plants with PCR using genomic DNA and RT-PCR respectively. As a positive control for PCR of transgenic soybean plants, we used vector plasmid DNA. As a negative control for PCR, genomic DNA of untransformed wild parental plants was used.

The results of the standard PCR analysis of PtMYB4 gene insertion into soybean genome of T₁ plants shown in Figure 2a, b, and PAL5 gene insertion into genome of T1 plants – shown in Figure 2c. RT-PCR analysis of transgenic soybean lines of first generation demonstrated expression of PtMYB4 gene in RNA level (Figure 2 d, e).

Overall, the presence of PtMYB4 gene insertion into genome was confirmed in 45 transgenic soybean lines, presence of PAL5 – in 30 transgenic lines, F5H – in 17 lines, CAD1 – in 20 transgenic lines. From 2000 analyzed putative transgenic soybean seeds we generated 112 T₁ transgenic plants with average efficiency of transformation 5,6 %. Importantly we did not observe any significant difference between the varieties in term of transformation efficiency, which proves that the developed method of germ-line transformation is not genome-specific. Data for other lignification genes is not shown.

Seeds of T₁ transgenic soybean plants were planted to obtain of second-generation plants (T₂) for confirmation of stable transformation by molecular and biochemical analyses.

PCR molecular detection of T2 transgenic soybean plants

About 1000 T₂ transgenic soybean plants, obtained from confirmed by PCR T₁ transgenic plants, were selected randomly and screened by PCR and multiplex PCR (Figure 3). As shown in Figure 3 a-c, a large number of T₂ transgenic soybean lines demonstrates the presence of the gene of interest (PtMYB4 and PAL5) (Figure 3 a, b), and marker gene nptII (Figure 3c). PCR of soybean plants of second generation T₂ have confirmed stable insertion of valuable PtMYB4 and PAL5 genes and marker gene nptII (Figure 3 c) into soybean genome.

As shown in Figure 3d, all genetic elements of CsVMV/PtMYB4sens. construct: the valuable gene, reporter LicB - (lichenase), marker nptII, T-nos terminator confirmed to be inserted into soybean genome in second generation (T₂) of transgenic plants. Multiplex PCR shows the main genetic elements included in gene construction through alignment at annealing temperature 48°C for all genetic elements (valuable gene, reporter, marker, t-nos terminator) of the CsVMV/PtMYB4sens, 35SVTM/PAL5, C4H/F5H, CsVMV/CAD1 constructs. Multiplex PCR demonstrated the integration of the whole transgenic construction into genome of transgenic plant. This confirms the stable insertion of whole genetic construction containing the genes of interest into genome of soybean transgenic lines. Overall, we confirmed the stable integration of all elements of genetic constructions of CsVMV/PtMYB4sens., 35S/PAL5, С4H/F5H, Cs/CAD1 genes into soybean genome with the efficiency more than 50%, in second generation T₂ of transgenic soybean lines.

Southern and northern blot analyses of genomic DNA and total RNA

The same putative transgenic plant examples were subjected to the Southern and Northern blot analysis (Figure 4). Total genomic DNA digested with BamHI. Hybridization with the PAL-probe confirmed the presence of the transgenes in all tested lines except PBI 121 PAL5 #1 and #3, which is in agreement with the PCR data (Figure 4a). Further we hybridized the same membrane with the gus-probe allowed seeing the pattern of gene insertion and estimating the transgene copy number (Figure 4b). All tested lines have multiple copies of the transgenes. The fact that bands are much wider in PBI 101 PAL5 line suggests the presence of a tandem repeats. There were no bands in negative control lines.

The same transgenic lines were used for total genomic RNA isolation and a Northern blot analysis. Northern blot with total genomic RNA of PBI 101 PAL5 and PBI 121 PAL5 lines was performed with untreated samples and samples treated with water (Figure 4 c). Hybridization with the pal probe revealed the presence of the genomic PAL in all samples, including control and the presence of transgenic band in PBI 121 PAL5 line #4 (Figure 4 c). Hybridization of the same membrane with the gus probe resulted in a band at the same place, suggesting the presence of bisictronic mRNA due the expression of a transgene at the transcriptional level (Figure 4 d). This band appears much weaker in untreated samples due to the significantly lower amount of RNA. Hybridization with the nptII probe revealed the presence of the nptII band in all samples, except the control as expected (Figure 4 e). Thus, the soybean transgenic lines are capable of the expression of introduced genes at the transcriptional level.

Lignin content and composition

Further we analyzed the content of lignin in the cell wall in CAD1 and PAL5 transformants. Lignin content in control lines was 23.17% of the dry weight, 28.8% in Cs/CAD1 transformants and 29.99% in 35S/PAL5 transgenic lines. This confirms that lignin content is increased in the presence of the target lignification genes expression.

Furthermore, we analyzed lignin composition in these transgenic lines. We observed an overall increase of lignin monomer content in 35S/PAL5 transgenic lines compared to the control lines. In the Cs/CAD1 transformants we observed an increase of S-monomers content compared to the control lines (Figure 5 a, b).

Additionally, we analyzed the activity of PAL and CAD enzymes expressed in the transgene plants to convert L-phenylalanine to trans-cinnamic acid and the conversion of coniferyl and synapil alcohols to the corresponding aldehydes, respectively. We observed a significant increase of PAL5 activity in PBI 101 PAL5 line #11 (Figure 5 c) as well as an increase of CAD1 activity in PBI 101 CAD1 line #5 (Figure 5 d). Correspondent with the transgene expression this proves the functional activity of the inserted genes of interest. All tested lines have multiple copies of the transgenes.

Discussion

We developed a novel germ-line technique for soybean genetic transformation using natural pollen tubes to introduce the key genes of lignification into mature but not yet divided zygote. This technique was used to perform a genetic transformation of 1070 flowers of 3 USA and 7 established and zoned in Kazakhstan soybean varieties. We obtained 833 pods (78% efficiency) with 2000 putative transgenic soybean seeds of zero generation (T₀). Seeds forming efficiency on average was 187%. We confirmed the expression of the key lignification genes as well as the expression of all elements in the constructs (CsVMV/PtMYB4sens., 35S/PAL5, С4H/F5H, Cs/CAD1) using PCR, RT-PCR, multiplex PCR in all 10 soybean varieties with the average efficiency of 5.6%. In the second generation of transgenic plants T₂ we observed a stable integration of all elements of genetic constructions of these genes into soybean genome with the efficiency more than 50%. All tested lines of PAL5 transformants have multiple copies of the transgenes as verified by Southern Blot and confirmed expression of transformed genes into soybean genome. Northern Blot analysis demonstrated expression of introduced genes (PAL, gus, nptII) at the transcriptional level. We observed the increase of lignin content from 23.17% of the dry weight in control to 28.8% in Cs/CAD1 transformants and 29.99% in PAL5 transformants. Additionally, PAL5 transformants demonstrated elevated levels of all lignin monomers (H-, G-, and S), and observed an increase of S-monomers in CAD1-transgenic lines compared to the control.

So, we have developed a novel soybean germ-line genetic transformation technique using natural pollen tubes as vectors for the Agrobacterium-mediated introduction of target transgenes into soybean genome. This method is based on previously described method of pollen-tube pathway transformation which is successfully used for crops like cotton or rice (Ali et al., 2015), however seem to be not efficient and irreducible in soybean (Shou et al., 2002). The major difference of the used method is using A. tumefaciens to introduce the vectors containing the genes of interest into the newly fertilized but not yet dividing zygote cell via pollen tubes.

The method has a high efficiency of T₁ plants generation containing the genes of interest (5.6% on average, compared to 2.8% dependent on genotype in transformation in cotyledonary nodes (Jia et al., 2015), or 3.8 (1.4–8.7%) (Paz et al., 2006). We confirmed the presence of the transgenes using standard molecular biology techniques such as PCR, RT-PCR and Southern and Northern blotting. One of the important findings is that the transformation efficiency seems to be independent on the genotype which is a big advantage of the method. However, further investigations on the efficiency of A. tumefaciens-mediated transduction using natural pollen tubes should be performed, including the selection of the most transductable genotypes and more rigorous analysis of the selected plants. We found several important advantages of the suggested method. Firstly, this method is similar with natural pollen hybridization and transgenic plants are generated within a short period of time (90–120 days). No tissue culture and plant regeneration steps are required which makes this technique much less expensive and simpler compared to the conventional approaches. The plants derived in this study are able to produce normal health seeds, comparable with control plants. We observed that the transgenic plants are healthy and are very similar to the plants derived using conventional breeding. Also, due to the transformation of the germ-line these plants are less likely chimeric therefore the genes of interest are transmitted to the next generation.

Importantly, as the method uses only single and simple transformation step it can be performed in greenhouses or in the field. This can greatly reduce the soybean transformation expenses on tissue culture and recovering steps thus the overall cost of transgenic plant production.

The described method was used to produce soybean plants with altered lignin concentrations by introduction of the genes related to the lignin synthesis. We created a number of transgenic cell lines stably expressing PtMYB4, PAL5, F5H, and CAD1. Importantly, we observed the an increased amount of lignin in the plant lines harboring PAL5 and CAD1 genes which are required for proper lignin formation. This demonstrates that the used technique is capable of producing transgenic lines with the desired phenotype. The next step is to determine if the changes in lignin content and lignin composition will provide the resistance to various pathogenes. Earlier it was shown that soybean plants with low stem lignin concentrations are more resistant to sclerotinia stem rot – one of the most common soybean fungal diseases (Peltier, Hatfield, & Grau, 2009). Therefore more studies should be performed to characterized the resistance profile of the derived transgenic plant lines. Other important line of the further studies is to identify the transgenic constructs that will increase soybean resistance and will not compromise plant qualities such as yield and nutritional value.

Conclusion

In summary, we developed a novel simple and effective germ-line transformation method for soybean genetic transformation and were able to introduce genes related to lignification which may be important for soybean biotic stress and disease resistance improvement. This simple and effective genetic transformation technique is expected to be beneficial for the Kazakhstan and US soybean growers and also world-wild and will help to generate great variation in phenotype, increase of biodiversity and establishment of different valuable soybean genotypes with enhanced stress tolerance and productivity for breeders.

Acknowledgements

Funding for this research was provided by project "Improving the sustainability of soy to biotic stress by genetic engineering of phenylpropanoid cycle" has made in the framework of sub-program 101 "Grant funding for research." Priority "intellectual potential of the country", MES, Kazakhstan, 2013-2015. As well as by Fulbright grant for senior visitor researcher Olga Kershanskaya, the US State Department for Scientific Research at the University of Illinois at Urbana-Champaign in the laboratory of JM Widholm, 2013-2014. This research also has been funded by the Science Committee of the Ministry of Education and Science of the Republic of Kazakhstan (Grant No. AP05132774. 2018-2020).

We express our deep gratitude to the staff of Lab. Prof. J.M. Widholm, Crop Science Dep., University of Illinois at Urbana-Champaign, USA, for providing genetic constructs of lignification genes, participation and assistance in the experiments.

References

Ali, A., Bang, S. W., Chung, S.-M., & Staub, J. E. (2015). Plant transformation via pollen tube-mediated gene transfer. Plant Molecular Biology Reporter, 33, 742–747. doi: 10.1007/s11105-014-0839-5

Barros, J., Serk, H., Granlund, I., & Pesquet, E. (2015). The cell biology of lignification in higher plants. Annals of Botany, 115(7), 1053-1074. doi: 10.1093/aob/mcv046

Baxter, A., Mittler, R., & Suzuki, N. (2014). ROS as key players in plant stress signalling. Journal of Experimental Botany, 65(5), 1229–1240. doi: 10.1093/jxb/ert375

Birkenbihl, R. P., Liu, S., & Somssich, I. E. (2017). Transcriptional events defining plant immune responses. Current Opinion in Plant Biology, 38, 1–9. doi: 10.1016/j.pbi.2017.04.004

Bonawitz, N. D., & Chapple, C. (2010). The genetics of lignin biosynthesis: connecting genotype to phenotype. Annual Review of Genetics, 44, 337–363. doi: 10.1146/annurev-genet-102209-163508

Cabane, M., Afif, D., & Hawkins, S. (2012). Lignins and abiotic stresses. In L. Jouanin & C. Lapierre (Eds.), Advances in Botanical Research Lignins (Vol. 61, p. 219–262). Burlington: Academic Press.

Cavallini, E., Matus, J. T., Finezzo, L., Zenoni, S., Loyola, R., Guzzo, F., … Tornielli, G. B. (2015). The phenylpropanoid pathway is controlled at different branches by a set of R2R3-MYB C2 repressors in grapevine. Plant Physiology, 167(4), 1448–1470. doi: 10.1104/pp.114.256172

Clough, S. J., & Bent, A. F. (1998). Floral dip: a simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana. The Plant Journal, 16(6), 735–743. doi: 10.1046/j.1365-313X.1998.00343.x

Edwards, R., & Kessman, H. (1992). Isoflavonoid phytoalexins and their biosynthetic enzymes. In S. J. Gurr, M. J. McPherson, & D. J. Bowles (Eds.), Molecular plant pathology: a practical approach (2nd ed., p. 45–62). Oxon: Oxford University Press.

Gui, J., Shen, J., & Li, L. (2011). Functional characterization of evolutionarily divergent 4-coumarate: Coenzyme a ligases in rice. Plant Physiology, 157(2), 574–586. doi: 10.1104/pp.111.178301

Hartman, G. L., Rupe, J. C., Sikora, E. J., Domier, L. L., Davis, J. A., & Steffey, K. L. (2016). Part 4: Soybean disease and pest management strategies. In G. L. Hartman, J. C. Rupe, E. J. Sikora, L. L. Domier, J. A. Davis & K. L. Steffey (Eds.), Compendium of soybean diseases and pests (5th ed., p. 167–173). Saint Paul: The American Phytopathological Society.

Hawkins, S. W, & Boudet, A. M. (1994). Purification and characterization of cinnamyl alcohol dehydrogenase isoforms from the periderm of Eucalyptus gunnii Hook. Plant Physiology, 104(1), 75–84. doi: 10.1104/pp.104.1.75

Homrich, M. S., Strohm, B. W., Weber, R. L. M., & Zanettini, M. H. B. (2012). Soybean genetic transformation: a valuable tool for the functional study of genes and the production of agronomically improved plants. Genetics and Molecular Biology, 35(4), 998-1010. doi: 10.1590/S1415-47572012000600015

Jia, Y., Yao, X., Zhao, M., Zhao, Q., Du, Y., Yu, C., & Xie, F. (2015). Comparison of soybean transformation efficiency and plant factors affecting transformation during the Agrobacterium infection process. International Journal of Molecular Sciences, 16(8), 18522–18543. doi: 10.3390/ijms160818522

Kandel, R., Chen, C. Y., Grau, C. R., Dorrance, A. E., Liu, J. Q., Wang, Y., & Wang, D. (2018). Soybean resistance to white mold: evaluation of soybean germplasm under different conditions and validation of QTL. Frontiers in Plant Science, 9, 505. doi: 10.3389/fpls.2018.00505

Kim, W.-C., Kim, J. Y., Ko, J.-H., Kang, H., & Han, K.-H. (2014). Identification of direct targets of transcription factor MYB46 provides insights into the transcriptional regulation of secondary wall biosynthesis. Plant Molecular Biology, 85(6), 589–599. doi: 10.1007/s11103-014-0205-x

Kissoudis, C., Van de Wiel, C., Visser, R. G. F., & Van der Linden, G. (2014). Enhancing crop resilience to combined abiotic and biotic stress through the dissection of physiological and molecular crosstalk. Frontiers in Plant Science, 5, 207. doi: 10.3389/fpls.2014.00207

Lee, M.-H., Jeon, H. S., Kim, S. H., Chung, J. H., Roppolo, D., Lee, H.-J., … Park, O. K. (2019). Lignin based barrier restricts pathogens to the infection site and confers resistance in plants. The EMBO Journal, 38(23), e101948. doi: 10.15252/embj.2019101948

Li, H., & Wu, T.-L. (2007). Transfoming Agrobacterium into soybean by means of pollen tube pathway induced by CaCl2. Soybean Science, 26(5), 55–59. doi: 10.11861/J.ISSN.1000-9841.2015.01.0036

Li, S., Cong, Y., Liu, Y., Wang, T., Shuai, Q., Chen, N., … Li, Y. (2017). Optimization of Agrobacterium-mediated transformation in soybean. Frontiers in Plant Science, 8, 246. doi: 10.3389/fpls.2017.00246

Li, Z., Nelson, R. L., Widholm, J. M., & Bent, A. (2002). Soybean transformation via the pollen tube pathway. Soybean Genetics Newsletter, 29, 1–11.

Liu, J., Su, Q., An, L., & Yang, A. (2009). Transfer of a minimal linear marker-free and vector-free smGFP cassette into soybean via ovary-drip transformation. Biotechnology Letters, 31, 295–303. doi: 10.1007/s10529-008-9851-x

Liu, Q., Luo, L., & Zheng, L. (2018). Lignins: biosynthesis and biological functions in plants. International Journal of Molecular Sciences, 19(2), 335. doi: 10.3390/ijms19020335

Lozovaya, V. V., Lygin, A. V., Zernova, O. V., Ulanov, A. V., Li, S., Hartman, G. L., & Widholm, J. M. (2007). Modification of phenolic metabolism in soybean hairy roots through down regulation of chalcone synthase or isoflavone synthase. Planta, 225(3), 665–679. doi: 10.1007/s00425-006-0368-z

Lozovaya, V. V., Lygin, A. V., Zernova, O. V., & Widholm, J. M. (2005). Genetic engineering of plant root disease resistance by modification of the phenylpropanoid pathway. Plant Biosystems, 139(1), 20–23. doi: 10.1080/11263500500055247

Lygin, A. V., Hill, C. B., Zernova, O. V., Crull, L., Widholm, J. M., Hartman, G. L., & Lozovaya, V. V. (2010). Response of soybean pathogens to glyceollin. Phytopathology, 100(9), 897–903. doi: 10.1094/PHYTO-100-9-0897

Lygin, A. V., Rahman, M. M. A., Ulanov, A. V., Widholm, J. M., & Lozovaya, V. V. (2012). Polyethylene glycol treatment promotes metabolic events associated with maize callus morphogenic competence. Phytochemistry, 82, 46–55. doi: 10.1016/j.phytochem.2012.07.002

Lygin, A. V., Zernova, O. V., Hill, C. B., Kholina, N. A., Widholm, J. M., Hartman, G. L., & Lozovaya, V. V. (2013). Glyceollin is an important component of soybean plant defense against Phytophthora sojae and Macrophomina phaseolina. Phytopathology, 103(10), 984–994. doi: 10.1094/PHYTO-12-12-0328-R

Mansell, R. L., Gross, G. G., Stöckigt, J., Franke, H., & Zenk, M. H. (1974). Purification and properties of cinnamyl alcohol dehydrogenase from higher plants involved in lignin biosynthesis. Phytochemistry, 13(11), 2427–2435. doi: 10.1016/S0031-9422(00)86917-4

McCabe, D. E., Swain, W. F., Martinell, B. J., & Christou, P. (1988). Stable transformation of soybean (Glycine max) by particle acceleration. Nature Biotechnology, 6(8), 923–926. doi: 10.1038/nbt0888-923

Meyer, K., Cusumano, J. C., Somerville, C., & Chapple, C. C. (1996). Ferulate-5-hydroxylase from Arabidopsis thaliana defines a new family of cytochrome P450-dependent monooxygenases. PNAS - Proceedings of the National Academy of Sciences of the United States of America, 93(14), 6869–6874. doi: 10.1073/pnas.93.14.6869

Meyer, K., Shirley, A. M., Cusumano, J. C., Lelong, D. A. B., & Chapple, C. (1998). Lignin monomer composition is determined by the expression of a cytochrome P450-dependent monooxygenase in Arabidopsis. PNAS - Proceedings of the National Academy of Sciences of the United States of America, 95(12), 6619–6623. doi: 10.1073/pnas.95.12.6619

Moura, J. C. M. S., Bonine, C. A. V., Viana, J. d. O. F., Dornelas, M. C., & Mazzafera, P. (2010). Abiotic and biotic stresses and changes in the lignin content and composition in plants. JIPB - Journal of Integrative Plant Biology, 52(4), 360–376. doi: 10.1111/j.1744-7909.2010.00892.x

Pagano, M. C., & Miransari, M. (2016). The importance of soybean production worldwide. In M. Miransari (Ed.), Abiotic and biotic stresses in soybean production (p. 1–26). London: Academic Press. doi: 10.1016/B978-0-12-801536-0.00001-3

Patzlaff, A., McInnis, S., Courtenay, A., Surman, C., Newman, L. J., Smith, C., … Campbell, M. M. (2003). Characterisation of a pine MYB that regulates lignification. The Plant Journal, 36(6), 743–754. doi: 10.1046/j.1365-313X.2003.01916.x

Paz, M. M., Martinez, J. C., Kalvig, A. B., Fonger, T. M., & Wang, K. (2006). Improved cotyledonary node method using an alternative explant derived from mature seed for efficient Agrobacterium-mediated soybean transformation. Plant Cell Reports, 25, 206–213. doi: 10.1007/s00299-005-0048-7

Peltier, A. J., Hatfield, R. D., & Grau, C. R. (2009). Soybean stem lignin concentration relates to resistance to Sclerotinia sclerotiorum. Plant Disease, 93(2), 149–154. doi: 10.1094/PDIS-93-2-0149

Pérez-Clemente, R. M., Vives, V., Zandalinas, S. I., López-Climent, M. F., Muñoz, V., & Gómez-Cadenas, A. (2013). Biotechnological approaches to study plant responses to stress. Biomed Research International, 1–10. doi: 10.1155/2013/654120

Piasecka, A., Rey, N. J., & Bednarek, P. (2015). Secondary metabolites in plant innate immunity: conserved function of divergent chemicals. New Phytologist, 206(3), 948-964. doi: 10.1111/nph.13325

Raes, J., Rohde, A., Christensen, J. H., Van De Peer, Y., & Boerjan, W. (2003). Genome-wide characterization of the lignification toolbox in Arabidopsis. Plant Physiology, 133(3), 1051–1071. doi: 10.1104/pp.103.026484

Ritter, H., & Schulz, G. E. (2004). Structural basis for the entrance into the phenylpropanoid metabolism catalyzed by phenylalanine ammonia-lyase. Plant Cell, 16(12), 3426–3436. doi: 10.1105/tpc.104.025288

Robinson, A. R., & Mansfield, S. D. (2009). Rapid analysis of poplar lignin monomer composition by a streamlined thioacidolysis procedure and near-infrared reflectance-based prediction modeling. The Plant Journal, 58(4), 706–714. doi: 10.1111/j.1365-313X.2009.03808.x

Sattler, S. E., & Harris, D. L. F. (2013). Modifying lignin to improve bioenergy feedstocks: strengthening the barrier against pathogens? Frontiers in Plant Science, 4, 70. doi: 10.3389/fpls.2013.00070

Shou, H., Palmer, R. G., & Wang, K. (2002). Irreproducibility of the soybean pollen-tube pathway transformation procedure. Plant Molecular Biology Reporter, 20(4), 325-334. doi: 10.1007/BF02772120

Tran, L.-S. P., & Mochida, K. (2010). Functional genomics of soybean for improvement of productivity in adverse conditions. Functional & Integrative Genomics, 10(4), 447-462. doi: 10.1007/s10142-010-0178-z

Van Acker, R., Vanholme, R., Storme, V., Mortimer, J. C., Dupree, P., & Boerjan, W. (2013). Lignin biosynthesis perturbations affect secondary cell wall composition and saccharification yield in Arabidopsis thaliana. Biotechnology for Biofuels, 6(46). doi: 10.1186/1754-6834-6-46

Vanholme, R., Storme, V., Vanholme, B., Sundin, L., Christensen, J. H., Goeminne, G., … Boerjan, W. (2012). A systems biology view of responses to lignin biosynthesis perturbations in Arabidopsis. The Plant Cell, 24(9), 3506–3529. doi: 10.1105/tpc.112.102574

Vogt, T. (2010). Phenylpropanoid Biosynthesis. Molecular Plant, 3(1), 2–20. doi: 10.1111/J.1365–313X.2009.03996.X

Wanner, L. A., Li, G., Ware, D., Somssich, I. E., & Davis, K. R. (1995). The phenylalanine ammonia-lyase gene family in Arabidopsis thaliana. Plant Molecular Biology, 27(2), 327–338. doi: 10.1007/BF00020187

Weigel, D., & Glazebrook, J. (2009). Dellaporta miniprep for plant DNA isolation. Cold Spring Harbor Protocols, 4(3). doi: 10.1101/pdb.prot5178

Weissgerber, T. L., Milic, N. M., Winham, S. J., Garovic, V. D. (2015). Beyond bar and line graphs: time for a new data presentation paradigm. PLOS Biology, 13(4), e1002128. doi: 10.371/ journal.bbio. 1002128

Yin, X., & Zhang, Z. J. (2010). Recent patents on plant transgenic technology. Recent Patents on Biotechnology, 4(2), 98–111. doi: 10.2174/187220810791110688

Zatybekov, A., Abugalieva, S., Didorenko, S., Rsaliyev, A., & Turuspekov, Y. (2018). GWAS of a soybean breeding collection from South East and South Kazakhstan for resistance to fungal diseases. Vavilov Journal of Genetics and Breeding, 22(5), 536–543. doi: 10.18699/VJ18.392

Notas de autor