
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39169119
70272
10.1038/s41598-024-70272-2
Article
Overexpression of Ginkbilobin-2 homologous domain gene improves tolerance to Phytophthora cinnamomi in somatic embryos of Quercus suber
Serrazina Susana 1
Martínez MªTeresa 2
Soudani Serine 3
Candeias Gonçalo 1
Berrocal-Lobo Marta 3
Piñeiro Pablo 2
Malhó Rui 1
Costa Rita Lourenço 4
Corredoira Elena elenac@mbg.csic.es

2
1 https://ror.org/01c27hj86 grid.9983.b 0000 0001 2181 4263 Faculdade de Ciências, BioISI-Biosystems & Integrative Sciences Institute, Universidade de Lisboa, Lisbon, Portugal
2 grid.502190.f 0000 0001 2292 6080 Misión Biológica de Galicia, Sede Santiago de Compostela, Consejo Superior de Investigaciones Científicas (MBG-CSIC), Avda Vigo S/N, 15705 Santiago de Compostela, La Coruña Spain
3 https://ror.org/03n6nwv02 grid.5690.a 0000 0001 2151 2978 Centro Para La Biodiversidad y Desarrollo Sostenible (CBDS), ETSIMFMN, Universidad Politécnica de Madrid, Ciudad Universitaria S/N, 28040 Madrid, Spain
4 https://ror.org/01fqrjt38 grid.420943.8 0000 0001 0190 2100 Instituto Nacional de Investigação Agrária E Veterinária I.P., Oeiras, Portugal
21 8 2024
21 8 2024
2024
14 1935718 5 2024
14 8 2024
© The Author(s) 2024, corrected publication 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
In recent decades an extensive mortality and decline of Quercus suber populations mainly caused by Phytophthora cinnamomi has been observed. In the current study, a chestnut gene homologous to ginkbilobin-2 (Cast_Gnk2-like), which in Ginkgo biloba codifies an antifungal protein, was transferred into cork oak somatic embryos of three different embryogenic lines by Agrobacterium mediated transformation. The transformation efficiency varied on the genotype from 2.5 to 9.2%, and a total of 22 independent transformed lines were obtained. The presence of Cast_Gnk2-like gene in transgenic embryos was verified in all lines by PCR. The number of transgene copies was estimated by qPCR in embryogenic lines with high proliferation ability and it varied between 1 and 5. In addition, the expression levels of Cast_Gnk2-like gene were determined in the embryogenic lines, with higher levels in lines derived from the genotype ALM6-WT. Transgenic plants were obtained from all transgenic lines and evaluated after cold storage of the somatic embryos for 2 months and subsequent transfer to germination medium. In vitro tolerance tests made under controlled conditions and following zoospore treatment showed that plants overexpressing Cast_Gnk2-like gene improved tolerance against Pc when compared to wild type ones.

Keywords

Agrobacterium tumefaciens
Cork oak
Cysteine-rich repeat secretory protein
Genetic transformation
In vitro tolerance assays
Oak decline disease
Subject terms

Biological techniques
Biotechnology
Plant sciences
FCT PORTUGAL2022.06990.PTDC Serrazina Susana http://dx.doi.org/10.13039/501100001871 Fundação para a Ciência e a Tecnologia UIDB/04046/2020 and UIDP/04046/2020 Ministerio de Ciencia, Innovación y Universidades (MICIN, Spain)AEI/10.13039/501100011033 UPM-Algerian Government from MHESRDecree No. 14–196 Soudani Serine Consejo Superior de Investigaciones Cientificas (CSIC)Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature.

issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Cork oak (Quercus suber L.) is a very familiar evergreen oak from the western Mediterranean, where it occupies an area estimated to be over 2.1 million hectares1. Portugal, Spain, Morocco and Algeria together host more than 90% of the species' distribution area2. Cork oak trees can appear forming forests called Sobreirais in Portuguese and Alcornocales in Spanish or forming part of the agroforestry ecosystems called Montados in Portuguese and Dehesas in Spanish. Cork oak forests and dehesas constitute ecologically and economically sustainable systems as they serve as an important tool in desertification prevention and generate high levels of biodiversity3. From an economic point of view, cork oak is the only species used for industrial cork production4. Cork has been used by humans for more than 2,000 years5 and is the sixth most economically important non-woody forest material at the global level. More than 80% of the world cork is produced in the Iberian Peninsula2,6. Cork is lightweight, water-resistant, flexible, chemically stable and fire-resistant7. It is also renewable, recyclable and biodegradable8. The greatest application of cork is the production of stoppers for the wine industry and champagne, but cork is also used as a construction material, for floors and thermal insulation2,6. Its applications are so broad that even NASA uses it for the insulation of spaceships, and surfboards and skateboards are also made from cork9.

However, since the early 90’s cork oaks have been seriously affected by the oak decline dieback which seriously threatened their populations with a resulting decrease in cork production (losses ranged from 40 to 89%)10,11. Oak decline is a multifactorial disease in which many abiotic and biotic factors such as drought, floods, frost, over-grazing, low regeneration ability, insect pests, and pathogens interact in a scenario of global climate change12. Among them, the oomycete Phytophthora cinnamomi Rands (Pc) is considered the main cause13. It seems that environmental changes provoke physiological decline of the trees, and this decline significantly increases the tree vulnerability to Pc, insects and other pathogens14. To date, traditional agronomic practices have not been effective in controlling infection caused by this oomycete and the generation of tolerant genotypes seems the most suitable option. However, conventional improvement programs for forest species face an enormous limitation: the long generation time of forest species, which in the case of oak species can exceed 30 years, and the high level of heterozygosis15,16. In the short future, the European Parliamentary Research Service (EPRS) recommends the relaxation of the rules for the use of transgenic plants obtained by targeted mutagenesis, cisgenesis or intragenesis17. Consequently, genetic transformation with antifungal genes isolated in the same species or a closely related species could be a good alternative to conventional breeding efforts.

The efficient introduction of new characters in plants involves the development of an effective transformation protocol and an available regeneration method that permits the generation of a plant from the transformed cells18. In hardwood species, somatic embryogenesis (SE) is considered the best regeneration method for genetic engineering and cryopreservation procedures19. In the last two decades, considerable advances have been reported on the development of SE in cork oak, even from explants derived from adult trees20,21. In contrast, genetic transformation has been less studied, especially for the overexpression of genes of interest15. Until now, genetic transformation of cork oak has been achieved using somatic embryos from both mature and juvenile materials introducing marker genes22,23 as well as a gene to confer herbicide resistance24. However, resistance has only been tested in somatic embryos at laboratory level, as plantlet conversion was not reported. As far as we are aware, the only report of genetic transformation in cork oak with genes to confer resistance to Pc is the overexpression of CsTL1, which codifies a thaumatin-like protein, a pathogenesis-related protein25.

Since the specific genes involved in oomycete defense have not yet been identified in cork oak, it is necessary to look for alternatives to induce some type of resistance against Pc. There are many indications that the overexpression of genes that codify antifungal agents may confer tolerance to pathogen attacks26–28. Under the term antifungal agents are pathogenesis-related proteins, antimicrobial proteins and peptides, DUF26-containing proteins, and other proteins such as polygalacturonase-inhibiting proteins or albumins29. DUF26-containing proteins are characterized as containing cysteine-rich domains with a conserved C-X8-C-X2-C motif in which lies its ability to fight plant pathogens and especially fungi29,30. The most well-known DUF26-containing protein is ginkbilobin-2 (Gnk2) from Ginkgo biloba seeds, which is predictly secreted and inhibits the growth in vitro of Fusarium oxysporum, Tricoderma reesei and Candida albicans31. Gnk2 showed similarity to the extracellular domain of plant cysteine-rich receptor-like kinases, but not to pathogenesis-related proteins32. Its antifungal activity is attributed to the binding of the DUF26 domain to sugar residues in the fungal cell wall33. A candidate gene for Pc resistance, homologous to Gnk2, was identified in the root transcriptome of the resistant chestnut Castanea crenata Sieb. and Zucc.34. In C. crenata, Cast_Gnk2-like had a high expression in non-inoculated roots and was distinctly up-regulated 48 h after Pc inoculation, contrasting with the feeble expression in C. sativa35. Furthermore, the recombinant protein encoded by Cast_Gnk2-like delays Pc mycelial growth in vitro36. Quercus ilex L. plants overexpressing Cast_Gnk2-like gene were able to survive more days in the presence of Pc when compared to non-transformed controls37. The studies carried out with Cast_Gnk2-like are revealing an important role in the defense of susceptible Fagaceae to Pc. Here, we describe the development of transgenic cork oak somatic embryos and plantlets overexpressing the Cast_Gnk2-like gene through Agrobacterium-mediated transformation. We also evaluate the resistance of transgenic plants to inoculation with P. cinnamomi.

Results

Somatic embryo genetic transformation

Somatic embryos from three embryogenic lines were precultured for one day, cocultured for 5 days with the bacteria and cultivated on kanamycin-containing medium during 10 weeks. The first kan-resistant somatic embryos were observed after 4 weeks of culture on kanamycin-containing medium but only in line TGR3 (Fig. 1a). The first kan-resistant somatic embryos in lines ALM6 and ALM80 were detected later, after 8 weeks on selection medium. The best percentages of kan-resistant explants (14.2%) were also obtained with TGR3 line, although without significant differences with the values achieved with lines ALM80 (11.7%) and ALM6 (10%) (Table 1). The transformation efficiency (i.e. the explants GFP +) was also higher in line TGR3 (9.2%) than in lines ALM80 (6.7%) and ALM6 (2.5%) (Table1; Fig. 1b,c), but similar to the percentage of kan-resistant explants, no significant differences were observed between the three lines (Table 1). The negative controls (non-transformed explants but cultured on kanamycin-containing medium) died without forming callus or new somatic embryos, which indicates that the selection pressure is adequate.Figure 1 Genetic transformation of somatic embryos of cork oak. (a) Kanamycin resistant explants observed after ten weeks on selection medium. Diameter of the Petri dish, 90 mm. (b) Group of transgenic somatic embryos observed under white light. (c) The same group of transgenic somatic embryos observed under blue light showing green fluorescence. (d) GFP expression in leaves isolated from a non-transgenic plant (left) and a transgenic plant (right). (e) Expression of GFP in roots from a non-transgenic plant (left) and from a transgenic plant (right). (f) Transgenic plantlets (left and right) and non-transgenic plant (center) obtained from somatic embryos of genotype TGR3 following to two months of cold storage and eight weeks on germination medium. Bar, 1 cm.

Table 1 Effect of genotype on the percentage of kanamycin-resistant explants, on the transformation efficiency and the selection efficiency following culture in selection medium of cork oak somatic embryos transformed with the strain EHA105pKWG2D-GIN.

Embryogenic line	Kanamycin resistant explants (%)1	Transformation efficiency (%)2	
TGR3	14.2 ± 2.8	9.2 ± 2.5	
ALM80	11.7 ± 3.1	6.7 ± 2.7	
ALM6	10.0 ± 3.5	2.5 ± 1.3	
ANOVA I	ns	ns	
Values are means ± SE of 12 Petri dishes with 10 explants each. 1Percentage of initial explants that are kanamycin-resistant after 10 weeks in selection medium; 2Percentage of initial explants that are GFP + after 14 weeks in selection medium.

ns: not significant.

After GFP evaluation, one somatic embryo at cotyledonary stage was isolated from each kan-resistant and GFP + explant. The selected embryos were proliferated independently on kanamycin-containing medium to establish the different transgenic lines. After, 2–3 subculture periods a total of 22 transgenic lines were obtained: 11 transgenic lines in TGR3, 8 transgenic lines in ALM80 and 3 transgenic lines in ALM6, that were maintained through secondary embryogenesis.

Once the different lines were established, the ability to form cotyledonary embryos larger than 5 mm (i.e. somatic embryos used for plant regeneration) was evaluated. In general, lines from ALM6 genotype showed more somatic embryos at cotyledonary stage than lines derived from genotypes TGR3 and ALM80 (Fig. 2). Among the three transgenic lines generated from genotype ALM6, the embryo production was significantly higher (p ≤ 0.05) in lines ALM6-GIN1 (8.5 se) and ALM6-GIN3 (8.3 se) (Fig. 2a). In genotype ALM80, the best results were obtained with lines ALM80-GIN2 (4.6 se) and ALM80-GIN9 (3.8 se) but without significant differences between the different lines generated (Fig. 2b). By contrast, significant differences (p ≤ 0.05) were observed between transgenic lines generated from TGR3 genotype and the embryo production was higher in lines TGR3-GIN1 and TGR3-GIN2 (Fig. 2c). Hence, these transformed lines from the 3 genotypes were selected for molecular analysis.Figure 2 Production of somatic embryos cotyledonary stage from non-transgenic and transgenic lines generated from ALM6 (a), ALM80 (b) and TGR3 (c) genotypes. In each genotype, columns represent means ± error standard of 3 Petri dishes with 8 explants each. The statistical analysis of these data by ANOVA I is show in the table.

Molecular analysis of transgenic lines

The presence of transgenes in the kan-resistant lines was confirmed by PCR analysis. The 472 bp of NPTII fragment, 740 bp of GFP fragment and 890 bp and 1227 bp of Cast_Gnk2-like fragments (i.e. the gene was amplified in both transcriptional senses; see Supplementary Information 1) were successfully amplified from all the transgenic lines and in the plasmid (positive control) confirming presence of the genes in the genome of somatic embryos (Fig. 3). No amplification product was observed in non-transgenic control somatic embryos of the three cork oak embryogenic lines evaluated (Fig. 3).Figure 3 Confirmation by PCR amplification of the presence of Cast_Gnk2-like gene. Cast_Gnk2-like-R (a), Cast_Gnk2-like-F (b), GFP (c) and NPTII (d) fragments in the transformed cork oak lines. M: 600 bp molecular weight marker; P: plasmid; wt: non-transformed line; lanes 1–2: amplification of the transformed lines of ALM6-GIN1 and -GIN3; lane 3–4: amplification of the transformed line of ALM80-GIN2 and -GIN9; lanes 5–6: amplification of the transformed line of TGR3-GIN1 and -GIN2; arrows indicate 600 pb.

Cast_Gnk2-like copy number in the transformed lines was estimated in somatic embryos with the amplification of a CaMV35S promoter fragment (Table 2). The correlation coefficient of the standard curve that plots the CT with the copy number was 0.89, despite the high efficiency obtained in the qPCR reactions (see Supplementary Information 1). The value of R2 may reflect the high variation in the detection of fewer copies, close to the limit of equipment detection. The highest number of transgene copies (5) was obtained for ALM6-GIN1, followed by TGR3 lines (3), ALM6-GIN3 (2) and ALM80 lines (1) (Table 2).Table 2 Estimated copy number in cork oak somatic embryogenic lines transformed with Cast_Gnk2-like gene.

Line	CT Mean	CT Standard deviation	Estimated copy number 1	
ALM6-WT	25.47	0.83	0	
ALM6-GIN1	18.60	0.83	5	
ALM6-GIN3	20.76	0.55	2	
ALM80-WT	26.30	0.77	0	
ALM80-GIN2	21.63	0.20	1	
ALM80-GIN9	21.27	1.38	1	
TGR3-WT	25.81	1.51	0	
TGR3-GIN1	20.02	0.79	3	
TGR3-GIN2	20.25	0.89	3	
1Values from the correlation between the quantity of plasmid and the number of copies in cork oak, as described in Material and Methods.

Transgene expression in each transformed line was compared to corresponding non-transformed genotypes in somatic embryos (Fig. 4). Cast_Gnk2-like oligos used in qPCR were based on the coding sequence of the original species, Castanea crenata. Q.suber and C. crenata are of the Fagaceae family, and when a BLASTn is executed in GenBank with Cast_Gnk2-like ORF as Query, the orthologous gene Q. suber cysteine-rich repeat secretory protein 38-like (Supplementary Information 2). Therefore, endogenous and transgenes transcripts may be amplified in the same reaction. But Cast_Gnk2-like had CaMV35S as the promoter, which has a constitutive action, so an overexpression of the transgene in transformed lines and an expression level higher than in the non-transformed lines was expected. This was true and significant for ALM6, TGR3-GIN1, and ALM80-GIN9 transformed lines (Fig. 4). ALM80-WT and TGR3-WT showed evident basal expression levels (4.3 and 1.8 respectively), and the transformed lines ALM80-GIN2 and TGR3-GIN2 didn’t show an overexpression of the transgene.Figure 4 Expression analysis of Cast_Gnk2-like gene in somatic embryos of cork oak. Expression level in transformed lines and non-transformed lines (WT) was normalized to Q. suber Elongation factor 1-α and Tubulin reference genes. Bars are mean ± standard error (n = 3). Asterisks indicate significant differences in the expression when compared to non-transformed lines with p < 0.05 (Tukey HSD).

Plant regeneration of transgenic lines

Plant regeneration ability was assessed in transgenic lines exhibiting elevated production of cotyledonary embryos compared to their respective non-transformed counterparts. Cotyledonary embryos were stored for two months at 4 °C under darkness conditions (see Supplementary Information 3) and subsequently they were cultured on germination medium. Regeneration was evaluated as simultaneous development of shoot and root (i.e. plant conversion), or as the development of only root (see Supplementary Information 3). After 8 weeks in the germination medium, only root development was observed in all tested cork oak lines, indicating uniformity in this aspect across the lines (Table 3). However, differences in the values of all evaluated parameters were observed among the lines (Table 3). The higher conversion frequencies were observed in line ALM6-WT and their corresponding transgenic lines with values ranging 83 to 92% (Table 3). Line TGR3-WT and their corresponding transgenic lines were the ones that showed the lowest conversion percentages (27–33%), whereas line ALM80-WT and their corresponding transgenic lines presented intermediate conversion values (23–53%). The quality of the plants evaluated in terms of root length, shoot length and number of leaves was higher in almost all transgenic lines compared to their wild-type counterparts.Table 3 Plant regeneration performance with only root development and conversion into plantlets (root and shoot development) of different cork oak transformed and untransformed lines after 8 weeks on germination medium.

Line	Root only	Conversion (shoot + root)	
(%)	RL (mm)	(%)	RL (mm)	SL (mm)	LN	
LINES ALM6	
ALM6-WT	12.5 ± 6.2	45.0 ± 2.2	83.3 ± 7.4	71.4 ± 9.3	17.5 ± 0.9	3.7 ± 0.4	
ALM6-GIN1	12.5 ± 9.7	80.0 ± 0.0	83.3 ± 9.1	119.4 ± 9.0	15.8 ± 1.5	7.9 ± 1.0	
ALM6-GIN3	4.2 ± 3.2	68.0 ± 0.0	91.7 ± 3.7	107.3 ± 13.4	20.8 ± 1.7	6.6 ± 0.6	
ANOVA I	ns	ns	ns	p ≤ 0.05	ns	ns	
LINES ALM80	
ALM80-WT	30.0 ± 8.7	31.2 ± 5.0	53.3 ± 8.7	44.3 ± 5.3	14.2 ± 2.4	6.3 ± 0.6	
ALM80-GIN2	53.3 ± 11.9	61.0 ± 16.0	23.3 ± 7.6	74.0 ± 19.5	16.9 ± 3.2	6.6 ± 0.3	
ALM80-GIN9	56.7 ± 12.1	38.6 ± 4.1	30.0 ± 9.9	34.9 ± 6.6	16.5 ± 1.7	8.1 ± 0.7	
ANOVA I	ns	ns	ns	ns	ns	ns	
LINES TGR3	
TGR3-WT	70.0 ± 10.7	40.9 ± 10.7	30.0 ± 12.8	50.0 ± 15.8	10.5 ± 1.7	5.2 ± 0.6	
TGR3-GIN1	53.3 ± 5.6	63.3 ± 11.9	33.3 ± 8.2	57.0 ± 16.8	15.3 ± 2.4	6.7 ± 1.3	
TGR3-GIN2	66.7 ± 10.5	49.4 ± 8.0	26.7 ± 12.1	114.9 ± 15.1	15.6 ± 0.3	5.3 ± 0.9	
ANOVA I	ns	ns	ns	ns	ns	ns	
Each value represents the mean ± standard error of 5 replications with 6 explants in each replicate. RL Root length, SL Shoot length, LN Leaf number, ns not significant.

GFP expression was also detected in roots and leaves isolated from transgenic plants (Fig. 1d,e). As expected, no green fluorescence was observed in roots and leaves isolated from non-transformed regenerated plants (Fig. 1d,e). No phenotypic differences were found relativeto untransformed control plants, indicating that the GFP gene does not provoke deleterious or toxic effects on cork oak cells (Fig. 1f).

Tolerance assay to Pc

To determine under controlled axenic conditions, if the expression of the Cast-Gnk2-like gene might improve tolerance to Pc, plants of ALM6-WT, ALM6-GIN1 and ALM6-GIN3 were inoculated with zoospores (see methods section below). Root necrosis and disease symptoms were evident on wild type plants relative to transformed plants in all assays performed (Fig. 5a,b). At four days post-inoculation, wild type plants decayed while transgenic plants had lighter root necrosis (Fig. 5c). After 4 days, all non-transgenic plants die and the tolerance experiments were concluded to obtain enough non-necrotized tissue to achieve RNA from the roots. At this time, callose deposition was higher in transgenic plants when compared to non-transgenic ones (Fig. 5d; Supplementary Information 4). Root necrosis was in accordance to the loss of fresh weight observed, being higher in wild-type plants in comparison to transformed plants (Fig. 5e,f). The transcriptional levels of Cast-Gnk2-like gene were quantified using qPCR, revealing a correlation between the expression of the Cast_Gnk2-like gene and increased plant tolerance to Pc in ALM6 plants compared to wild types, as illustrated in Fig. 5g.Figure 5 Disease symptoms of Quercus suber responding to Phytophthora cinnamomi Rands (Pc). ALM6-WT plants and in transgenic ALM6-GIN1 plants were inoculated with zoospores (107 zs/ml). a Symptoms on explants at time “0” (T0), and after 4 days post-inoculation (4 dpi). Bar: (1 cm). b Detail of cell death symptoms on fresh leaves and stained with Trypan Blue (TB). Bars: (1 mm). c Kinetic of disease symptoms expressed in percent related to the total number of explants per treatment (n = 6), at 1, 2, 3 and 4 dpi, where: 0. No symptoms, 1. Leaf chlorosis and light necrosis on roots, 2. Apparent necrosis on leaves and light necrosis on roots, 3. High necrosis on leaves and roots, 4. Decayed explant. d Callose deposition on roots (1, differentiation, 2, elongation, and 3, meristematic area), detected by Aniline Blue staining on WT (left) and transgenic ALM6-GIN 1 (middle), and ALM6-GIN 3 (right), explants, inoculated (left panels) with zoospores or non-inoculated (right panels). (100 × magnification). Bar: (1 mm). e Root necrosis percentage on explants at 4dpi. f Fresh weight lost expressed in percentage related to controls at 4dpi. g Quantification of the expression of Cast_GNK2-like gene in WT and transgenic plants infected (Pc) and non-infected (C), data were normalized with β-Tubulin reference gene. The * in e denotes a statistically significant difference using variance check (P-value ≤ 0.05) and Tukey HSD test.

Discussion

Oak decline syndrome, mainly caused by P. cinnamomi, has provoked serious damage in cork oak populations and therefore considerable economic losses. To date, cultural practices and chemical methods to control Pc were scarcely effective38. Additionally, fungicides are environmentally detrimental and expensive39. In this context, genetic transformation offers an alternative not only generating transgenic plants with enhanced disease tolerance, but also to provide us a method to explore the function of putative genes involved on the response mechanisms of plants to pathogens40,41. Although transgenic strategies have an evident potential to enhance cork oak disease tolerance, to date, only a gene with antifungal activity has been evaluated in this species25. In the present study, the first challenge was the overexpression of Cast_Gnk2-like gene, a novel type of antifungal protein, isolated from C. crenata transcriptome by using somatic embryos as target material. An efficient in vitro regeneration system is imperative to produce transgenic plants42. In this regard, somatic embryogenesis is considered the ideal system to address genetic transformation in woody plants, especially in highly recalcitrant species43,44. In the present report, genetic transformation was achieved in the three embryogenic lines evaluated although results obtained here showed that transformation frequencies with Cast_Gnk2-like gene varied widely with the genotype of embryogenic lines (2.5–9.2%). The scarce protocols previously described in the literature for overexpression of interesting genes in cork oak were also highly dependent on the genotype of the embryogenic lines. In line with our current investigation, Alvarez et al. 45 documented significant variances among the six embryogenic lines during transformation with marker genes. Similarly, Cano et al. 25 identified the influence of genotype on the transformation process using the same embryogenic lines featured in this report, particularly when overexpressing the CsTL1 gene, encoding a thaumatin-like protein. As in this previous report, genotype TGR3 produced the best transformation frequencies, which reinforces the role of genotype in transformation. The effect of genotype on transformation rates was also mentioned for other members of Fagaceae family such as chestnut46–50, pedunculated oak51 and holm oak52. Although transformation rates are affected by the genotype, the values obtained in our report were higher than those described for other species using the same gene. Previously, we achieved transformation values of 2.5% using nodular embryogenic structures of holm oak transformed with Cast_Gnk2-like gene37. McGuigan et al.50 obtained 14 putative embryogenic lines expressing the Cast_Gnk2-like gene. However, only 3 were confirmed as containing the T-DNA insertion using RITA® Temporary Immersion Bioreactors to transform American chestnut somatic embryos.

As for holm oak, in the present study, GFP and NTPII genes were also successfully introduced and expressed in somatic embryos of the three embryogenic lines evaluated. This double selection with kan and GFP provided, respectively, a good selection and excellent visual screening of transgenic tissues in cork oak, as high levels of GFP expression were observed in both somatic embryos and in transgenic plants. Also, in both tissue types, non-toxic or deleterious effects were detected. This result points out that transgenic cork oak embryos could be selected using GFP expression alone even though the addition of kanamycin could improve the efficiency. This efficiency explains why GFP has replaced β‐glucuronidase (GUS) method as reporter procedure in transgenic plants due to its rapid detection, its high sensitivity, which does not harm the plant in which it is expressed and allows live visualization of tissues, using a UV light and without tissue destruction53,54.

The estimated Cast_Gnk2-like gene copies in transgenic lines varied from 1 to 5. Similarly, Alvarez et al.24 found variations from 1 to 4 copies of the inserted phosphinothricin acetyl transferase among transformed cork oak lines. A relationship between the number of transgenes and the expression level was not found. Other factors besides the copy number influence the expression level of the transgenes. One of them, the promoter CaMV35S, was used to confer high and constitutive expression of Cast_Gnk2-like. In effect, transgenic lines with different levels of Cast_Gnk2-like overexpression were obtained for all the somatic embryo genotypes, which was an important achievement for the Cast_Gnk2-like functional analysis approach. CaMV35S-driven transgene expression can be highly variable in the different transformants with the same construct55. Other factors that influence transgene expression are the position and the integration of the T-DNA in the genome. For instance, if multiple T-DNAs are integrated in one locus, they can be organized as direct or inverted repeats, resulting in different levels of expression55.

Remarkably, gene expression levels observed in non-transformed ALM80 and TGR3 indicate a high basal level of a endogenous Q. suber ginkbilobin2-like. The alignment of the aminoacid sequences of such endogenous protein with Cast_Gnk2-like show a high similarity between them (95%, Supplementary Information 2). In common there is a predicted signal peptide for secretion and 2 conserved cysteine-rich domains (C-X8-C-X2-C). The second domain is equal, although in the first domain there is one dissimilar aminoacid (lysine for the endogenous protein, threonine for C. crenata). The 2 aminoacids have different R-groups, lysine has an electrically charged side chain and threonine has a polar uncharged side chain. Could this difference contribute to confer different resistance levels to Q. suber and C. crenata? Our results indicate that the low basal expression level of the gene in non-transformed ALM6 is not sufficient to confer tolerance to P. cinnamomi. Additional assays will be necessary to determine if Q. suber ginkbilobin2-like has a relevant role against the pathogen.

In addition to Cast_Gnk2-like gene overexpression in Q. suber somatic embryos, a second challenge was to generate viable plants from transgenic somatic embryos to use for tolerance experiments. The ability to regenerate complete plants from transformed cotyledonary-stage embryos is one of the main bottlenecks in woody plants as it has been mentioned in previous studies, as in the genetic transformation of somatic embryos of European chestnut48 or pedunculate oak51. Nevertheless, in the present report suitable conversion percentages have been obtained, especially in transgenic lines derived from ALM6 genotype. Moreover, the plants obtained showed a good quality which allowed us to have adequate material for the evaluation of their tolerance to P. cinnamomi.

The selection criteria for conducting tolerance analysis in transgenic lines often hinge on the expression levels of the transgene. However, determining the optimal criterion remains a subject of debate. Some studies advocate for evaluating transgenic lines with higher protein expression levels, arguing that elevated expression increases the likelihood of achieving tolerance against pathogens25,56–58. Conversely, other researchers have failed to observe a clear correlation between tolerance and gene expression levels59. In our study, we opted to select embryogenic lines for tolerance analysis based on their high conversion rates and root length.

Overexpression of the Cast_Gnk2-like gene potentially confers improved tolerance to P. cinnamomi in cork oak plants using an in vitro test. Evaluating biotic and abiotic tolerance in plants poses significant challenges, particularly in woody species60. In vitro testing has emerged as a valuable method for assessing tolerance, offering the advantages of handling a larger number of specimens in a confined space and precise control over environmental conditions61. While in vitro tests have predominantly been employed for assessing tolerance to abiotic stress, there are instances of their application in evaluating tolerance to biotic stress, albeit primarily in herbaceous species such as carnation62, Arabidopsis63, and rice64. In our study, we conducted an in vitro tolerance test by inoculating cork oak plant roots with a standardized stock of zoospores under controlled conditions to minimize additional variability that may arise from using mycelium fragments 65, building upon prior methodologies66–68. Zoospore inoculation offers a faster and precise mimicry of P. cinnamomi infection in natural conditions. Our in vitro assays revealed that plants overexpressing the Cast_Gnk2-like gene exhibited reduced disease symptoms compared to untransformed plants, with approximately half the root necrosis observed in wild-type plants, along with increased callose deposition. Callose deposition is recognized as a key stress response in plants69, with its role in defense responses being contingent upon the pathogen's mode of infection. For instance, callose deposition in avocado roots' parenchyma and cortex is crucial for tolerance by impeding root zoospore colonization70, a phenomenon observed in maize as well 71. Callose deposition at plasmodesmata serves to impede cell-to-cell spreading induced by viruses and blast fungus72, restricting the spread of Phytophthora brassicae infection73 and potentially limiting cell-to-cell trafficking, thus safeguarding the plant host. However, mutants of the Powdery Mildew Resistance (PMR4) callose synthase gene exhibited enhanced tolerance to oidium74, suggesting a context-dependent role of callose. Cahill and Weste75 demonstrated in over 30 species, including both forest and non-forest species, that P. cinnamomi infection induced by zoospore inoculation led to callose deposition and papillae production around germ tube entry, exclusively in tolerant hosts, while being absent in susceptible ones. Our results, obtained under controlled conditions, indicate that callose deposition is induced in transgenic lines post-infection, whereas it was scarcely detectable in wild-type roots, the latter of which exhibited near-total root death. We found a positive correlation between callose deposition induction and the overexpression of the Cast_Gnk2-like gene, suggesting that this gene alters callose deposition patterns, thereby enhancing plant tolerance. Contrasting with inoculated plants, callose deposition in non-inoculated transgenic plants remained lower than in wild-type plants, indicating that callose deposition is modulated by gene overexpression. While overexpression of stress-related genes may contribute to root callose deposition, the correlation with increased tolerance is not always straightforward76, likely dependent on the pathogen's infection mechanism. On other hand, besides the role in stress response, callose also plays a crucial role in plant development, e.g. in cell plate formation during cell division (Chen and Kim77). In wild type non-inoculated plants, callose deposition can be interpreted as an adaptation to the in vitro inoculation, focusing on growth and cell division. In contrast, wild type inoculated plants do not exhibit callose signals because the cells have decayed due to infection. Our findings suggest that Cast_Gnk2-like gene overexpression in Q. suber contributes to callose deposition in response to P. cinnamomi, though the molecular underpinnings necessitate further investigation.

To confirm the high levels of Cast_Gnk2-like gene expression in ALM6-GIN 1 and GIN 3 lines during inoculation assays, qPCR was performed. As anticipated, both lines exhibited consistently high gene expression levels throughout infection, with a slight non-significant increase observed in infected plants. No significant differences were observed between the two lines, likely due to the already high levels of overexpression in both.

To date, only two studies have investigated the overexpression of a GNK-2 gene, and in line with our results, both reported enhanced tolerance. For instance, overexpression of ginkbilobin2-1 (GNK2-1) from Ginkgo biloba seed kernels in transgenic cucumber plants inhibited fungus growth when infected with Fusarium oxysporum32. Similarly, overexpressing the Cast_Gnk2-like gene in nodular embryogenic structures of holm oak significantly increased the survival period during infection with Pc37.

In intragenesis and cisgenesis procedures, genes transferred contain the same gene pool as in conventional breeding78. As Cast_Gnk2-like gene has been identified in Japanese chestnut, we would be overexpressing a gene from the same family as the cork oak, and consequently the plants obtained will not be strictly transgenic, being able to call intragenic if Cast_Gnk2-like gene was driven by a constitutive plant promoter (i.e. Ubiquitin 10 (UBQ10) promoter of Arabidopsis). In the last decade, the scientific community has showed great interest in promoting less stringent regulations for cisgenic/intragenic plants79–81. Several surveys and focus group interviews in the USA, Europe and other countries show that intragenic and cisgenic plants are more acceptable than strict transgenic plants, encouraging their use78,81. Hence, less rigorous plant transgenic regulations could reduce the costs and boost their commercial use78.

Material and methods

Plant material

For the genetic transformation, three embryogenic lines (ALM6, ALM80 and TGR3) induced from leaves derived from centenary trees were used20. These embryogenic lines have been maintained by secondary embryogenesis, with subcultures every 6 weeks, in proliferation medium (for more information see Supplementary Information 5). The embryogenic lines were cultured in a growth chamber under a 16-h photoperiod (provided by cool-white fluorescent lamps at a photon flux density of 50–60 µmol m-2 s-1) and 25 °C light/20 °C dark temperatures (standard conditions).

Transformation vector and bacterial strain

The Cast_Gnk2-like gene, which encodes a ginkbilobin-like protein, was cloned into the pK7WG2D vector82 under the CaMV35S promoter, using the Gateway system (Invitrogen, USA), as described Corredoira et al.83. The plasmid also includes the neomycin phosphotransferase gene (NPTII), and the green fluorescent protein (EGFP) reporter gene, used as a visual indicator of transformation. The NPTII gene is regulated by the promoter of the nopaline synthase (Nos) gene, whereas the EGFP gene is controlled by the promoter of the rol root loci D (rolD) gene. The vector with the three genes was called pK7WG2D-GIN (see Supplementary Information 6) and was transferred to the Agrobacterium tumefaciens strain EHA10584 using the freeze–thaw method defined by Xu and Li85. The resulting strain, named, EHA105pK7WG2D-GIN, was employed in all transformation experiments performed in the present report.

Somatic embryo transformation

For transformation, groups of 2 or 3 somatic embryos (4–7 mg; se) in the globular or torpedo stage of the three embryogenic lines were isolated 4 weeks after the last subculture and were precultured in 24 h in Petri dishes with preculture medium under standard conditions. The precultured explants were immersed in the infection medium with EHA105pK7WG2D-GIN which was prepared as described in Corredoira et al.83. After 30 min with gentle agitation, infection medium was removed by filtration and the explants were transferred to coculture medium and cocultured for 5 days in darkness at 25 °C. Bacteria were eliminated by immersion of explants for 30 min in a washing solution consisting of sterile water with carbenicillin (CB; 300 mg/L). The washing solution was removed by filtration and the explants were transferred to selection medium that consisted of proliferation medium supplemented with 300 mg/L CB and 125 mg/L of kanamycin (kan) (for more information see Supplementary Information 5). The embryos were maintained in standard conditions, with periodic subcultures every two weeks to this selection medium.

In each embryogenic line, 120 explants were used (6 Petri dishes with 10 explants and the experiment was repeated twice) resulting in a total of 360 explants. Additionally, for each line 20 explants that were not infected and that were cultured on proliferation medium without antibiotics (control positive) and with antibiotics (negative control).

Selection and maintenance of transgenic lines

After 10 weeks on selection medium, explants that showed signs of growth were transferred to proliferation medium (Fig. 1a) supplemented with a higher concentration of kan than in the selection medium (150 mg/L) and were cultivated for 4 more weeks, in order to increase the selective pressure and confirm its resistance to kan. After 14 weeks from the beginning of the transformation, the percentage of kan-resistant (defined as the percentage of initial explants that are resistant to kan) and the transformation efficiency (defined as the percentage of initial explants that showed fluorescence (GFP +) were determined. The fluorescence was observed using a Leica M205 FA magnifier (Germany), equipped with a 200-W bulb, and a specific filter for fluorescence, with 470/40 × nm excitation and 525/50 nm emission. The images were taken with a Leica DSC7000T camera (Germany). From each GFP + explant, a single somatic embryo at cotyledonary stage was isolated, and each embryo was multiplied by secondary embryogenesis in selection medium to establish the different transgenic embryogenic lines. During establishment period, GFP expression was also verified periodically to rule out the appearance of chimeras. After approximately four months, when the multiplication of GFP + lines was obtained, they were transferred to medium without antibiotic, to eliminate selective pressure and increase multiplication rates. To test the proliferation ability, the number of cotyledonary somatic embryos > 5 mm produced by explant was evaluated and compared with those of the corresponding non-transformed (control) embryogenic line.

Molecular analysis of transgenic embryogenic lines

Gene analysis presence in embryogenic lines

Total genomic DNA was isolated from 100 mg of cotyledonary somatic embryos of kan-resistant embryogenic lines and non-transformed embryogenic lines (WT) using the Kit REALPURE (Durviz, Spain) following the instructions of the manufacturer. PCR was performed to amplify a specific DNA sequence in transgenic lines corresponding to the GFP, NTPII and Cast_Gnk2-like genes (see Supplementary Information 1). A 25 µl PCR reaction was prepared containing 250–500 ng of genomic DNA, 15 µM of primers, 0.2 μL Taq DNA Polymerase (Qiagen, Germany), 1μL MgCl2, and 2.5 mM of dNTPs. Plasmid DNA was used as a positive control, whereas DNA from a non-transformed control somatic embryos and distilled water served as negative controls. PCR products were separated in 1.5% (w/v) agarose gels (Ultra Pure Agarose, Invitrogen, USA) in 1 × TBE buffer at 90 V. After electrophoresis, gels were stained with 1% de RedSafe (iNtRON Biotechnology, Korea) and visualized under the UV light.

Gene number copy analysis in embryogenic lines

The number of Cast_Gnk2-like copies inserted in the transformed lines was estimated by quantitative real-time PCR (qPCR), similarly as described in37. In brief, DNA was extracted as in the previous section, and three biological replicates per genotype were prepared. Specific oligos for the CaMV35S promoter were used in the reactions (see Supplementary Information 1). A standard curve was established as in Song et al.86, by mixing the plasmid used for transformation with non-transformed genomic DNA. Using the value of Q. suber genome size reported by Ramos et al.87, we calculated the quantity of plasmid needed to be mixed with non-transgenic DNA to simulate 1, 2, 5 and 10 copies of transgene. qPCR cycling with a subsequent step for the melting curve was performed as in Mallón et al.88. The qPCR experiments included no template controls, two technical replicates and two repetitions. The copy number in each sample was estimated making the correspondence between the CT and the mocked copies of the transgene.

Gene expression analysis in embryogenic lines

The expression of Cast_Gnk2-like gene was analyzed by qPCR. Total RNA was extracted from non-transformed and putatively transformed lines using somatic embryos at early cotyledonary stage and a protocol with CTAB (hexadecyltrimethylammonium bromide buffer) based on89,90. Three biological replicates were considered from three different embryo clumps, for both non-transformed and transformed lines. Once RNA was extracted, it was treated with the RNase-free DNase Set (Qiagen, Hilden, Germany) to eliminate possible DNA contamination. The Qiagen RNeasy Plant Mini Kit (Qiagen, Hilden, Germany) was then used for RNA Cleanup according to the manufacturer’s instructions. The stability of the extracted RNA was verified by observing a double band in gels of 1% agarose, and the concentration was quantified with a Thermo Scientific™ NanoDrop™ One Microvolume UV–Vis Spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). Total RNA (2 μg) was used as the template for reverse transcription with RevertAid H Minus Reverse Transcriptase (Thermo Fisher Scientific, Waltham, MA, USA) and primed with an oligo(dT) primer. Specific primers for Cast_Gnk2-like are described in Supplementary Information 1. 20 ng of cDNA were used per reaction in a 15 μL final volume using 7.5 μL of NZYSupreme qPCR Green Master Mix (2x) (NZYtech, Carnide, Lisboa, Portugal). A final concentration of 0.2 μM of each primer was used in a CFX96 Touch Real-Time PCR Detection System (BioRad, Hercules, CA, USA). Reactions started with a denaturation step at 95 °C for 10 min followed by 40 cycles of denaturation at 95 °C for 15 s and annealing temperature for 30 s. Each set of reactions included a no-template control and three technical replicates. Dissociation curves were used to analyze non-specific PCR products. To normalize expression data, Elongation factor 1-alpha (EF1α)91 and β-Tubulin (TUB)92 were used (Supplementary information 1). Gene expression was calculated using the Pfaffl method93.

Plant regeneration from transgenic lines

For plant regeneration from cork oak somatic, previously defined procedure94 was used. Briefly, somatic embryos at cotyledonary stage (≥ 5 mm) were isolated from 6-week cultures of non-transformed lines and the transformed lines with the best proliferation rates. Embryos were transferred to 100 ml glass jars with plastic lids with 30 ml of proliferation medium and were stored at 4 °C in darkness. After two months, the embryos were transferred to 500 ml glass jars with 70 ml of medium germination (for more information see Supplementary Information 5). After 8 weeks of culture under standard conditions the following parameters were determined: the percentage of embryos that only developed roots ≥ 5 mm, and the percentage of embryos that developed a complete plant, as well as the length of the root (mm) and the shoot (mm) and the number of leaves per plant. For each WT line and their respective transgenic lines 36 somatic embryos (6 se per jar) were cultured.

In vitro tolerance assays

Plant incubation conditions

Plants derived from germination of somatic embryos as described above were transferred from agar to plastic sterile tubes containing 1 ml of sterilized water and they were incubated in an Aralab chamber (Lisbon, Portugal), at 50% humidity (v/v), temperature of 24 °C during the day and 18 °C during the night, with a 16-h light/8-h dark photoperiod and light intensity of 150 µE. m-2 per second for all experiments. They were maintained in these conditions for a week.

Plant infection

The explants from agar tubes medium are incubated in an Aralab chamber (Lisbon, Portugal), at 50% humidity (v/v), temperature of 24 °C during the day and 18 °C during the night, with a 16-h light/8-h dark photoperiod and light intensity of 150 µE.m-2 per second for all experiments for a week, before to be transferred to new sterile culture tubes DeWit® (Sgl, Barcelona, Spain), containing 1 ml of sterile water and maintained for an additional two days for further acclimation. Phytophtora cinnamomi Rands was kindly provided by TRAGSA (Maceda nursery, Orense, Spain) and isolated by the Center for Research and Technology of Extremadura (CYCYTEX), located in Mérida, Spain.The Phytophtora cinnamomi mycelium culture was initially obtained from a fresh plate using previous describe methodology66–68. Then a zoospore stock was obtained, and saved at − 80 °C until use, following the previous optimized protocol65. He plants were divided into two groups: a control group and an inoculated group. The explants were inoculated into the sterile culture tubes DeWit® (Sgl, Barcelona, Spain), with Phytophthora cinnamomi zoospores stock at a final concentration of 107 zoospores/ml or with sterilized water on controls.

Disease evaluation symptoms

The plant disease symptoms produced by Phytophthora cinnamomi Rands were considered following the standards approved for diagnostic protocols for regulated pests by EPPO council95. Those symptoms included principally root rot brown lesions followed by necrosis produced directly by the pathogen and secondary symptoms of decline producing leave chlorosis and later necrosis and dead. Disease symptoms were followed for 4 days using a modified infection criteria described by96 as followed: 0. No symptoms, 1. Leaf chlorosis and light necrosis on roots, 2. Apparent necrosis on leaves and light necrosis on roots, 3. High necrosis on leaves and roots, 4. Decayed explants.

Fresh weight (FW) was measured at 0 and 4 dpi, to calculate the fresh weight lost (FWL) Following the formula:FW%=FW4dpi×100FWT0;FWL%=100%-FW%

Photographs of the explants were captured for assessing root necrosis. The quantification of root infection and necrosis involved measuring the total length of the infected root and the segments showing necrosis at 4 dpi using the software ImageJ. The calculations were performed using the following formula:R

Rootnecrosis(\%)=(Lengthofinfectedarea(cm)×100)Totallength(cm)

Necrotic lesions were measured by Trypan blue staining, performed as previously described97. Briefly, a Trypan Blue solution was prepared by combining 10 ml lactic acid (85% w:w), 10 ml phenol (TE balanced buffer, pH 7.5–8.0), 10 ml glycerol (≥ 99%), 10 ml distilled water, and 40 mg Trypan Blue, resulting in a final concentration of 10 mg/ml. Roots and leaves were stained for 20 min, rinsed with 100% ethanol over two days, and preserved in 60% glycerol until microscopic observation.

Callose deposition was detected by Aniline Blue staining applying an adapted procedure of the method described by98. Briefly, 0.1% aniline blue was prepared in 0.1 M of monosodium phosphate buffer (NaH2PO4·2H2O; pH = 9) and 2% glycerol (v/v). The root explants were immersed in the solution overnight in darkness. The stained samples were stored at 4 °C until microscopic observation.

Aniline blue was then detected using a DAPI/UV filter by fluorescence microscopy using a stereomicroscope (A292/21 Microscopy iScope IS.3153-PLFi/6 with Fluorescence—IS.3153- PLi/6,nEWF 10x/22, with Plan Fluarex PLFi, 4 × , 10 × , 20 × , 40 × , and 100 × oil lenses including fluorescence: Blue, Green, UV-DAPI, and Red filters, Microsercon SLU, Madrid, Spain) with a charge-coupled device (CCD) digital cooled camera (A292/21 Euromex 20 MP USB 3.0 with 1 inch CMOS sensor), to obtain digital photos. At least five photos were quantified per root area and treatment, considering taken on three different root areas of differentiation, elongation, and meristematic region and three plants per treatment. No callose was detected into apical root area.

Quantitative real-time-PCR analyses in infected plants

Tissues derived from inoculated plants were frozen using liquid nitrogen and subsequently pulverized with a mortar and pestle to obtain a powder, which was then transferred into Eppendorf tubes. Total RNA extraction was conducted in two steps. In the first step, the method with CTAB buffer already described was employed. The second step involved purifying the RNA solution from the first step using the RNeasy Plus Mini Kit (QIAGEN, Germany) in accordance with the kit's protocol. cDNA was synthesized using NZY First-Strand cDNA Synthesis Kit according to the manufacturer’s protocol.

The qPCR experiments were performed using a fluorescence green Power Sybr Green PCR Master mix (Thermofisher, USA) with reactions at a final volume of 20 µl per well. Specific primers for Cast-Gnk2-like and for Quercus suber β-Tubulin as a referencegene are described in Supplementary Information 1. The successful specificity of both oligonucleotides amplified in Quercus suber was confirmed by getting a melting curve at the last step of corresponding qPCR reactions.

The Cast_Gnk2-like gene expression in the transgenic plants inoculated or non-inoculated was quantified using the 2 − λλCT method relatively to wild-type explants and using three biological triplicates to each condition and assay. Data points were compared using t-tests (see below). Three biological replicates were used from two different individuals for both control and inoculated entire plants, in ALM6-WT and their transgenic lines (ALM6-GIN1 and ALM6-GIN3), respectively.

Statistical analysis

The SPSS 26.0 statistical package for Windows (Chicago, USA) was used to analyze the transformation and regeneration experiments. In transformation experiments, one-way analysis of variance (ANOVA I) was applied to determine the effect of genotype on genetic transformation of somatic embryos (Table 1) and the effect of genotype on the embryo production (Fig. 2) and plant regeneration (Table 3). In gene expression analysis–embryogenic lines, ANOVA I complemented by a Tukey HSD test was used for statistical significance analysis of the transgenic lines compared to the wild-type genotypes at p < 0.05 (Fig. 4). In the tolerance assay, the Stat Graphics Centurion XVI.II program (Stat Point Technologies, Inc., Warrenton, VA, United States) was used for all data analysis related to plant fresh weight lost, root necrosis and the expression of Cast_Gnk2-like gene. ANOVA I and Duncan’s mean comparison test were performed for all experiments and t-tests with a significance level of 0.05% (Fig. 5). Percentages were subjected to square-root transformation prior to analysis to normalize the data. Non-transformed data are presented in Tables and Figures. In the case of non-homogeneous variance, a non-parametric Kruskal–Wallis test was used.

Conclusions

The results reported here are preliminary, but they indicate that Cast_Gnk2-like gene has an important role in the tolerance to P. cinnamomi and provides new insights to investigate the infection process of this oomycete. Furthermore, Cast_Gnk2-like gene is a promising tool for engineering plants for pathogen tolerance. However, the mechanism of this gene and its role in cork oak tolerance requires further research. While the plants generated in this study incorporate viral promoters and reporter/marker genes, the promising outcomes achieved with GFP selection and the incorporation of a gene isolated from a taxonomically close species suggest the potential for future cultivation of cork oak plants overexpressing this gene as enhanced and well-adapted varieties.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-70272-2.

Acknowledgements

Authors thank Fátima Mosteiro and Sol Campañó for their technical assistance.

Author contributions

E.C. and T.M. conceived the study, planned the experiments, designed the constructs, and performed the cloning, genetic transformation, and plant regeneration. P.P. performed PCR analysis. S.S., G.C., R.M., and R.C. performed the transgene copy number and gene expression analysis. S.S. and M.B.-L. performed assays of tolerance in vitro and qPCR and callose assays on infected explants forming part of her PhD advised by MB-L at UPM. All authors contributed to writing of the manuscript and approved the submitted version.

Funding

Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature. The research described in the present paper was partly supported by Ministerio de Ciencia, Innovación y Universidades (MICIN, Spain) through the project PID2020-112627RB-C33 (AEI/10.13039/501100011033). SeS was funded by an international fellow from the UPM-Algerian Government from MHESR, Decree No. 14–196. SuS is being funded by an FCT contract in the project 2022.06990.PTDC (DOI: 10.54499/2022.06990.PTDC). At BioISI, work was supported by UIDB/04046/2020 (DOI: 10.54499/UIDB/04046/2020) and UIDP/04046/2020 (DOI: 10.54499/UIDP/04046/2020) Centre grants from FCT, Portugal.

Data availability

All data generated or analyzed during this study are included in this article and its supplementary information file. The images of callose do not correspond to protein blots, corresponds to aniline blue staining so the images are just provided directly by the fluorescence microscope. Because of the long size of root explants, the amplification zoom corresponds to 10 × so the original images are obtained directly from it. We provided a link file with original images obtained into this work including replicates used for quantifications of snapshots corresponding to aniline blue staining. Original callose imagenes have been included.

Competing interests

The authors declare no competing interests.

The original online version of this Article was revised: The original version of this Article contained an error in the spelling of the author Gonçalo Candeias which was incorrectly given as Gonçlo Candeias.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Change history

9/12/2024

A Correction to this paper has been published: 10.1038/s41598-024-72759-4
==== Refs
References

1. Dettori S Filigheddu MR Deplano G Molgora JE Ruiu M Sedda L Employing a spatio-temporal contingency table for the analysis of cork oak cover change in the Sa Serra region of Sardinia Sci. Rep. 2018 8 1 16946 10.1038/s41598-018-35319-1 30446680
Dettori, S. et al. Employing a spatio-temporal contingency table for the analysis of cork oak cover change in the Sa Serra region of Sardinia. Sci. Rep. 8(1), 16946. 10.1038/s41598-018-35319-1 (2018).30446680 10.1038/s41598-018-35319-1
2. Sierra-Pérez J Boschmonart-Rives J Gabarrell X Production and trade analysis in the Iberian cork sector: economic characterization of a forest industry Resour. Conserv. Recycl. 2015 98 55 66 10.1016/j.resconrec.2015.02.011
Sierra-Pérez, J., Boschmonart-Rives, J. & Gabarrell, X. Production and trade analysis in the Iberian cork sector: economic characterization of a forest industry. Resour. Conserv. Recycl. 98, 55–66. 10.1016/j.resconrec.2015.02.011 (2015).10.1016/j.resconrec.2015.02.011
3. Bugalho MN Plieninger T Aronson J Ellatifi M Crespo DG Aronson J Pereira JS Pausas JG Open woodlands: a diversity of uses (and Overuses) Cork Oak Woodlands on the Edge 2009 Island Press 33 44
Bugalho, M. N., Plieninger, T., Aronson, J., Ellatifi, M. & Crespo, D. G. Open woodlands: a diversity of uses (and Overuses). In Cork Oak Woodlands on the Edge (eds Aronson, J. et al.) 33–44 (Island Press, 2009).
4. EUFORGEN. Distribution map of cork oak (Quercus suber). www.euforgen.org (2009).
5. Gibson IJ Farries AM The mechanics of three-dimensional cellular materials Proc. R. Soc. Lond. A 1982 382 43 59 10.1098/rspa.1982.0088
Gibson, I. J. & Farries, A. M. The mechanics of three-dimensional cellular materials. Proc. R. Soc. Lond. A 382, 43–59. 10.1098/rspa.1982.0088 (1982).10.1098/rspa.1982.0088
6. Gil L Cork: A strategic material Front. Chem. 2014 2 16 10.3389/fchem.2014.00016 24790984
Gil, L. Cork: A strategic material. Front. Chem. 2, 16. 10.3389/fchem.2014.00016 (2014).24790984 10.3389/fchem.2014.00016
7. Pereira H Cork: Biology, Production and Uses 2007 Elsevier
Pereira, H. Cork: Biology, Production and Uses (Elsevier, 2007).
8. Duarte AP Bordado JC Cork – a renewable raw material: Forecast of industrial potential and development priorities Front. Mater. 2015 2 2 10.3389/fmats.2015.00002
Duarte, A. P. & Bordado, J. C. Cork – a renewable raw material: Forecast of industrial potential and development priorities. Front. Mater. 2, 2. 10.3389/fmats.2015.00002 (2015).10.3389/fmats.2015.00002
9. Amorim Cork Composites S.A. (n.d.). A prancha de surf de cortiça de McNamara. Retrieved october 24, 2022 from Amorim Cork Composites: https://amorimcorkcomposites.com/pt/inova%C3%A7%C3%A3o/casos-deestudo/a-prancha-de-surf-de-corti%C3%A7a-de-mcnamara/ (2022).
10. Sousa EMR Kadiri ZL déclenchement des perturbations physiologiques des peuplements de chêne-liège: une synthèse des relations agent/hôte IOBC/wprs Bull. 2005 28 9 16
Sousa, E. M. R. & Kadiri, Z. L. déclenchement des perturbations physiologiques des peuplements de chêne-liège: une synthèse des relations agent/hôte. IOBC/wprs Bull. 28, 9–16 (2005).
11. Tiberi R Branco M Bracalini M Croci F Panzavolta T Cork oak pests: A review of insect damage and management Ann. For. Sci. 2016 73 219 232 10.1007/s13595-015-0534-1
Tiberi, R., Branco, M., Bracalini, M., Croci, F. & Panzavolta, T. Cork oak pests: A review of insect damage and management. Ann. For. Sci. 73, 219–232. 10.1007/s13595-015-0534-1 (2016).10.1007/s13595-015-0534-1
12. Moricca S Linaldeddu BT Ginetti B Scanu B Franceschini A Ragazzi A Endemic and emerging pathogens threatening cork oak trees: Management options for conserving a unique forest ecosystem Plant Dis. 2016 100 2184 2193 10.1094/PDIS-03-16-0408-FE 30682920
Moricca, S. et al. Endemic and emerging pathogens threatening cork oak trees: Management options for conserving a unique forest ecosystem. Plant Dis. 100, 2184–2193 (2016).30682920 10.1094/PDIS-03-16-0408-FE
13. Corcobado T Miranda-Torres JJ Martín-García J Jung T Solla A Early survival of Quercus ilex subspecies from different populations after infections and co-infections by multiple Phytophthora species Plant Pathol. 2017 66 92 804 10.1111/ppa.12627
Corcobado, T., Miranda-Torres, J. J., Martín-García, J., Jung, T. & Solla, A. Early survival of Quercus ilex subspecies from different populations after infections and co-infections by multiple Phytophthora species. Plant Pathol. 66, 92–804. 10.1111/ppa.12627 (2017).10.1111/ppa.12627
14. Rodríguez-Molina MC Blanco-Santos A Palo-Núñez EJ Torres-Vila LM Torres-Álvarez E Suárez-de-la-Cámara MA Seasonal and spatial mortality patterns of holm oak seedlings in a reforested soil infected with Phytophthora cinnamomi For. Pathol. 2005 35 41 422 10.1111/j.1439-0329.2005.00423.x
Rodríguez-Molina, M. C. et al. Seasonal and spatial mortality patterns of holm oak seedlings in a reforested soil infected with Phytophthora cinnamomi. For. Pathol. 35, 41–422 (2005).10.1111/j.1439-0329.2005.00423.x
15. Vieitez AM Corredoira E Martínez MT San-José MC Sánchez C Valladares S Vidal N Ballester A Application of biotechnological tools to Quercus improvement Eur. J. For. Res. 2012 131 519 539 10.1007/s10342-011-0526-0
Vieitez, A. M. et al. Application of biotechnological tools to Quercus improvement. Eur. J. For. Res. 131, 519–539. 10.1007/s10342-011-0526-0 (2012).10.1007/s10342-011-0526-0
16. Martínez MT San José MC Arrillaga I Cano V Morcillo M Cernadas MJ Corredoira E Holm oak somatic embryogenesis: Current status and future perspectives Front. Plant Sci. 2019 10 239 10.3389/fpls.2019.00239 30967881
Martínez, M. T. et al. Holm oak somatic embryogenesis: Current status and future perspectives. Front. Plant Sci. 10, 239. 10.3389/fpls.2019.00239 (2019).30967881 10.3389/fpls.2019.00239
17. Ivana, K. Plants obtained by certain new genomic techniques, EPRS: European Parliamentary Research Service. Belgium. Retrieved from https://policycommons.net/artifacts/11321567/plants-obtained-by-certain-new-genomic techniques/12207361/ on 14 Feb 2024. CID: 20.500.12592/4mw6s5r (2024).
18. Corredoira E Costa RL Application of tissue culture in plant reproduction Forests 2021 12 342 10.3390/f12030342
Corredoira, E. & Costa, R. L. Application of tissue culture in plant reproduction. Forests 12, 342. 10.3390/f12030342 (2021).10.3390/f12030342
19. Corredoira E Merkle SA Martínez MT Toribio M Canhoto JM Correia SI Ballester A Vieitez AM Non-zygotic embryogenesis in hardwood species Crit. Rev. Plant Sci. 2019 38 29 97 10.1080/07352689.2018.1551122
Corredoira, E. et al. Non-zygotic embryogenesis in hardwood species. Crit. Rev. Plant Sci. 38, 29–97. 10.1080/07352689.2018.1551122 (2019).10.1080/07352689.2018.1551122
20. Hernández I Celestino C Toribio M Vegetative propagation of Quercus suber L. by somatic embryogenesis. I. Factors affecting the induction in leaves from mature cork oak trees Plant Cell Rep. 2003 21 759 764 10.1007/s00299-003-0605-x 12789519
Hernández, I., Celestino, C. & Toribio, M. Vegetative propagation of Quercus suber L. by somatic embryogenesis. I. Factors affecting the induction in leaves from mature cork oak trees. Plant Cell Rep. 21, 759–764. 10.1007/s00299-003-0605-x (2003).12789519 10.1007/s00299-003-0605-x
21. Hernández I Cuenca B Carneros E Alonso-Blázquez N Ruiz M Celestino C Ocaña L Alegre J Toribio M Application of plant regeneration of selected cork oak trees by somatic embryogenesis to implement multivarietal forestry for cork production Tree For. Sci. Biotechnol. 2011 5 19 26
Hernández, I. et al. Application of plant regeneration of selected cork oak trees by somatic embryogenesis to implement multivarietal forestry for cork production. Tree For. Sci. Biotechnol. 5, 19–26 (2011).
22. Alvarez R Alonso P Cortizo M Celestino C Hernández I Toribio T Ordás RJ Genetic transformation of selected mature cork oak (Quercus suber L.) Plant Cell Rep. 2004 23 218 223 10.1007/s00299-004-0810-2 15185122
Alvarez, R. et al. Genetic transformation of selected mature cork oak (Quercus suber L.). Plant Cell Rep. 23, 218–223. 10.1007/s00299-004-0810-2 (2004).15185122 10.1007/s00299-004-0810-2
23. Sánchez N Manzanera JA Pintos B Bueno MA Agrobacterium-mediated transformation of cork oak (Quercus suber L) somatic embryos New For. 2005 29 169 176 10.1007/s11056-005-0208-1
Sánchez, N., Manzanera, J. A., Pintos, B. & Bueno, M. A. Agrobacterium-mediated transformation of cork oak (Quercus suber L) somatic embryos. New For. 29, 169–176 (2005).10.1007/s11056-005-0208-1
24. Alvarez R Alvarez JM Humara JM Revilla A Ordás RJ Genetic transformation of cork oak (Quercus suber L.) for herbicide resistance Biotechnol. Lett. 2009 31 1477 1483 10.1007/s10529-009-0033-2 19543858
Alvarez, R., Alvarez, J. M., Humara, J. M., Revilla, A. & Ordás, R. J. Genetic transformation of cork oak (Quercus suber L.) for herbicide resistance. Biotechnol. Lett. 31, 1477–1483. 10.1007/s10529-009-0033-2 (2009).19543858 10.1007/s10529-009-0033-2
25. Cano V Martínez MT Couselo JL Varas E Vieitez FJ Corredoira E Efficient transformation of somatic embryos and regeneration of cork oak plantlets with a gene (CsTL1) encoding a chestnut thaumatin-like protein Int. J. Mol. Sci. 2021 22 1757 10.3390/ijms22041757 33578757
Cano, V. et al. Efficient transformation of somatic embryos and regeneration of cork oak plantlets with a gene (CsTL1) encoding a chestnut thaumatin-like protein. Int. J. Mol. Sci. 22, 1757. 10.3390/ijms22041757 (2021).33578757 10.3390/ijms22041757
26. Stuiver M Custers J Engineering disease resistance in plants Nature 2001 411 865 868 10.1038/35081200 11459071
Stuiver, M. & Custers, J. Engineering disease resistance in plants. Nature 411, 865–868. 10.1038/35081200 (2001).11459071 10.1038/35081200
27. Selitrennikoff CP Antifungal proteins Appl. Environ. Microbiol. 2001 67 2883 2894 10.1128/AEM.67.7.2883-2894.2001 11425698
Selitrennikoff, C. P. Antifungal proteins. Appl. Environ. Microbiol. 67, 2883–2894 (2001).11425698 10.1128/AEM.67.7.2883-2894.2001
28. Ali S Ganai BA Kamili AN Bhat AA Mir ZA Bhat JA Tyagi A Islam ST Mushtaq M Yadav P Rawat S Grover A Pathogenesis-related proteins and peptides as promising tools for engineering plants with multiple stress tolerance Microbiol. Res. 2018 212–213 29 37 10.1016/j.micres.2018.04.008 29853166
Ali, S. et al. Pathogenesis-related proteins and peptides as promising tools for engineering plants with multiple stress tolerance. Microbiol. Res. 212–213, 29–37. 10.1016/j.micres.2018.04.008 (2018).29853166 10.1016/j.micres.2018.04.008
29. Chiu T Poucet T Li Y The potential of plant proteins as antifungal agents for agricultural applications Synth. Syst. Biotechnol. 2022 7 4 1075 1083 10.1016/j.synbio.2022.06.009 35891944
Chiu, T., Poucet, T. & Li, Y. The potential of plant proteins as antifungal agents for agricultural applications. Synth. Syst. Biotechnol. 7(4), 1075–1083. 10.1016/j.synbio.2022.06.009 (2022).35891944 10.1016/j.synbio.2022.06.009
30. Vaattovaara A Brandt B Rajaraman S Safronov O Veidenberg A Mechanistic insights into the evolution of DUF26-containing proteins in land plants Commun. Biol. 2019 2 56 10.1038/s42003-019-0306-9 30775457
Vaattovaara, A. et al. Mechanistic insights into the evolution of DUF26-containing proteins in land plants. Commun. Biol. 2, 56. 10.1038/s42003-019-0306-9 (2019).30775457 10.1038/s42003-019-0306-9
31. Sawano Y Miyakawa T Yamazaki H Tanokura M Hatano K Purification, characterization, and molecular gene cloning of an antifungal protein from Ginkgo biloba seeds Biol. Chem. 2007 388 273 280 10.1515/BC.2007.030 17338634
Sawano, Y., Miyakawa, T., Yamazaki, H., Tanokura, M. & Hatano, K. Purification, characterization, and molecular gene cloning of an antifungal protein from Ginkgo biloba seeds. Biol. Chem. 388, 273–280 (2007).17338634 10.1515/BC.2007.030
32. Liu J Tian H Wang Y Guo A Overexpression of a novel antifungal protein gene GNK2-1 results in elevated resistance of transgenic cucumber to Fusarium oxysporum Chin. Bull. Bot. 2010 45 411 418
Liu, J., Tian, H., Wang, Y. & Guo, A. Overexpression of a novel antifungal protein gene GNK2-1 results in elevated resistance of transgenic cucumber to Fusarium oxysporum. Chin. Bull. Bot. 45, 411–418 (2010).
33. Miyakawa T Hatano K-I Miyauchi Y Suwa Y-I Sawano Y Tanokura MA Secreted protein with plant-specific cysteine-rich motif functions as a mannose-binding lectin that exhibits antifungal activity Plant Physiol. 2014 166 766 778 10.1104/pp.114.242636 25139159
Miyakawa, T. et al. Secreted protein with plant-specific cysteine-rich motif functions as a mannose-binding lectin that exhibits antifungal activity. Plant Physiol. 166, 766–778. 10.1104/pp.114.242636 (2014).25139159 10.1104/pp.114.242636
34. Serrazina S Santos C Machado H Pesquita C Vicentini R Pais MS Sebastiana M Costa R Castanea root transcriptome in response to Phytophthora cinnamomi challenge Tree Genet. Genomes 2015 11 6 10.1007/s11295-014-0829-7
Serrazina, S. et al. Castanea root transcriptome in response to Phytophthora cinnamomi challenge. Tree Genet. Genomes 11, 6. 10.1007/s11295-014-0829-7 (2015).10.1007/s11295-014-0829-7
35. Santos C Duarte S Tedesco S Fevereiro P Costa RL Expression profiling of castanea genes during resistant and susceptible interactions with the oomycete pathogen Phytophthora cinnamomi reveal possible mechanisms of immunity Front. Plant Sci. 2017 8 515 10.3389/fpls.2017.00515 28443110
Santos, C., Duarte, S., Tedesco, S., Fevereiro, P. & Costa, R. L. Expression profiling of castanea genes during resistant and susceptible interactions with the oomycete pathogen Phytophthora cinnamomi reveal possible mechanisms of immunity. Front. Plant Sci. 8, 515. 10.3389/fpls.2017.00515 (2017).28443110 10.3389/fpls.2017.00515
36. Colavolpe MB Vaz Dias F Serrazina S Malhó R Lourenço Costa R Castanea crenata Ginkbilobin-2-like recombinant protein reveals potential as an antimicrobial against Phytophthora cinnamomi, the causal agent of ink disease in European chestnut Forests 2023 14 785 10.3390/f14040785
Colavolpe, M. B., Vaz Dias, F., Serrazina, S., Malhó, R. & Lourenço Costa, R. Castanea crenata Ginkbilobin-2-like recombinant protein reveals potential as an antimicrobial against Phytophthora cinnamomi, the causal agent of ink disease in European chestnut. Forests 14, 785. 10.3390/f14040785 (2023).10.3390/f14040785
37. Serrazina S Martínez MT Cano V Malhó R Costa RL Corredoira E Genetic transformation of Quercus ilex somatic embryos with a Gnk2-like protein that reveals a putative anti-oomycete action Plants 2022 11 304 10.3390/plants11030304 35161285
Serrazina, S. et al. Genetic transformation of Quercus ilex somatic embryos with a Gnk2-like protein that reveals a putative anti-oomycete action. Plants 11, 304. 10.3390/plants11030304 (2022).35161285 10.3390/plants11030304
38. Martínez MT Cuenca B Mosteiro F Piñeiro P Pérez F Solla A Corredoira E Screening of cork oak for resistance to Phytophthora cinnamomi and micropropagation of tolerant seedlings Horticulturae 2023 9 692 10.3390/horticulturae9060692
Martínez, M. T. et al. Screening of cork oak for resistance to Phytophthora cinnamomi and micropropagation of tolerant seedlings. Horticulturae 9, 692. 10.3390/horticulturae9060692 (2023).10.3390/horticulturae9060692
39. Ali S Mir ZA Tyagi A Mehari H Meena RP Bhat JA Yadav P Papalou P Rawat S Grover A Overexpression of NPR1 in Brassica juncea confers broad spectrum resistance to fungal pathogens Front Plant Sci. 2017 8 1693 10.3389/fpls.2017.01693 29046679
Ali, S. et al. Overexpression of NPR1 in Brassica juncea confers broad spectrum resistance to fungal pathogens. Front Plant Sci. 8, 1693. 10.3389/fpls.2017.01693 (2017).29046679 10.3389/fpls.2017.01693
40. Li H Liu D He H Zhang N Ge F Chen C Molecular cloning of a 14–3-3 protein gene from Lilium regale Wilson and overexpression of this gene in tobacco increased resistance to pathogenic fungi Sci. Hort. 2014 168 9 16 10.1016/j.scienta.2013.12.034
Li, H. et al. Molecular cloning of a 14–3-3 protein gene from Lilium regale Wilson and overexpression of this gene in tobacco increased resistance to pathogenic fungi. Sci. Hort. 168, 9–16. 10.1016/j.scienta.2013.12.034 (2014).10.1016/j.scienta.2013.12.034
41. Wang Z Han Q Zi Q Lv S Qiu D Zeng H Enhanced disease resistance and drought tolerance in transgenic rice plants overexpressing protein elicitors from Magnaporthe oryzae PLoS One 2017 12 4 e0175734 10.1371/journal.pone.0175734 28419172
Wang, Z. et al. Enhanced disease resistance and drought tolerance in transgenic rice plants overexpressing protein elicitors from Magnaporthe oryzae. PLoS One 12(4), e0175734. 10.1371/journal.pone.0175734 (2017).28419172 10.1371/journal.pone.0175734
42. Nyaboga EN Njiru JM Tripathi L Factors influencing somatic embryogenesis, regeneration, and Agrobacterium-mediated transformation of cassava (Manihot esculenta Crantz) cultivar TME14 Front. Plant Sci. 2015 6 411 10.3389/fpls.2015.00411 26113851
Nyaboga, E. N., Njiru, J. M. & Tripathi, L. Factors influencing somatic embryogenesis, regeneration, and Agrobacterium-mediated transformation of cassava (Manihot esculenta Crantz) cultivar TME14. Front. Plant Sci. 6, 411. 10.3389/fpls.2015.00411 (2015).26113851 10.3389/fpls.2015.00411
43. Giri CC Shyamkumar B Anjaneyulu C Progress in tissue culture, genetic transformation and applications of biotechnology to trees: An overview Trees 2004 18 115 135 10.1007/s00468-003-0287-6
Giri, C. C., Shyamkumar, B. & Anjaneyulu, C. Progress in tissue culture, genetic transformation and applications of biotechnology to trees: An overview. Trees 18, 115–135. 10.1007/s00468-003-0287-6 (2004).10.1007/s00468-003-0287-6
44. Hazubska-Przybyl T Bojarczuk K Tree somatic embryogenesis in science and forestry Dendrobiology 2016 76 105 116 10.12657/denbio.076.010
Hazubska-Przybyl, T. & Bojarczuk, K. Tree somatic embryogenesis in science and forestry. Dendrobiology 76, 105–116 (2016).10.12657/denbio.076.010
45. Alvarez R Ordás RJ Improved genetic transformation protocol for cork oak (Quercus suber L.) Plant Cell Tissue Organ Cult. 2007 91 45 52 10.1007/s11240-007-9276-6
Alvarez, R. & Ordás, R. J. Improved genetic transformation protocol for cork oak (Quercus suber L.). Plant Cell Tissue Organ Cult. 91, 45–52 (2007).10.1007/s11240-007-9276-6
46. Polin LD Liang H Rothrock RE Nishii M Diehl DL Newhouse AE Mairn CJ Powell WA Maynard CA Agrobacterium-mediated transformation of American chestnut (Castanea dentata (Marsh) Borkh.) somatic embryos Plant Cell Tissue Organ Cult. 2006 84 69 79 10.1007/s11240-005-9002-1
Polin, L. D. et al. Agrobacterium-mediated transformation of American chestnut (Castanea dentata (Marsh) Borkh.) somatic embryos. Plant Cell Tissue Organ Cult. 84, 69–79 (2006).10.1007/s11240-005-9002-1
47. Andrade GM Nairn CJ Le HT Merkle SA Sexually mature transgenic American chestnut trees via embryogenic suspension-based transformation Plant Cell Rep. 2009 28 1385 1397 10.1007/s00299-009-0738-7 19578855
Andrade, G. M., Nairn, C. J., Le, H. T. & Merkle, S. A. Sexually mature transgenic American chestnut trees via embryogenic suspension-based transformation. Plant Cell Rep. 28, 1385–1397 (2009).19578855 10.1007/s00299-009-0738-7
48. Corredoira E San José MC Vieitez AM Allona I Aragoncillo C Ballester A Agrobacterium-mediated transformation of European chestnut somatic embryos with a Castanea sativa (Mill.) endochitinase gene New For. 2016 47 669 684 10.1007/s11056-016-9537-5
Corredoira, E. et al. Agrobacterium-mediated transformation of European chestnut somatic embryos with a Castanea sativa (Mill.) endochitinase gene. New For. 47, 669–684. 10.1007/s11056-016-9537-5 (2016).10.1007/s11056-016-9537-5
49. Corredoira E San-José MC Vieitez AM Ballester A Improving genetic transformation of European chestnut and cryopreservation of transgenic lines Plant Cell Tissue Organ Cult. 2007 91 281 288 10.1007/s11240-007-9295-3
Corredoira, E., San-José, M. C., Vieitez, A. M. & Ballester, A. Improving genetic transformation of European chestnut and cryopreservation of transgenic lines. Plant Cell Tissue Organ Cult. 91, 281–288 (2007).10.1007/s11240-007-9295-3
50. McGuigan L Fernandes P Oakes A Stewart K Powell W Transformation of american chestnut (Castanea dentata (Marsh) Borkh) using RITA® temporary immersion bioreactors and we vitro containers Forests 2020 11 11 1196 10.3390/f11111196
McGuigan, L., Fernandes, P., Oakes, A., Stewart, K. & Powell, W. Transformation of american chestnut (Castanea dentata (Marsh) Borkh) using RITA® temporary immersion bioreactors and we vitro containers. Forests 11(11), 1196. 10.3390/f11111196 (2020).10.3390/f11111196
51. Mallón R Valladares S Corredoira E Vieitez AM Vidal N Overexpression of the chestnut CsTL1 gene coding for a thaumatin-like protein in somatic embryos of Quercus robur Plant Cell Tissue Organ Cult. 2014 116 141 151 10.1007/s11240-013-0390-3
Mallón, R., Valladares, S., Corredoira, E., Vieitez, A. M. & Vidal, N. Overexpression of the chestnut CsTL1 gene coding for a thaumatin-like protein in somatic embryos of Quercus robur. Plant Cell Tissue Organ Cult. 116, 141–151. 10.1007/s11240-013-0390-3 (2014).10.1007/s11240-013-0390-3
52. Cano V Martínez MT San José MC Couselo JL Varas E Bouza-Morcillo L Toribio M Corredoira E Regeneration of transgenic plants by Agrobacterium-mediated transformation of Quercus ilex L. somatic embryos with the gene CsTL1 New For. 2020 51 1003 1021 10.1007/s11056-020-09771-9
Cano, V. et al. Regeneration of transgenic plants by Agrobacterium-mediated transformation of Quercus ilex L. somatic embryos with the gene CsTL1. New For. 51, 1003–1021. 10.1007/s11056-020-09771-9 (2020).10.1007/s11056-020-09771-9
53. Newell CA Plant transformation technology. Developments and applications Mol. Biotechnol. 2000 16 1 53 65 10.1385/MB:16:1:53 11098468
Newell, C. A. Plant transformation technology. Developments and applications. Mol. Biotechnol. 16(1), 53–65. 10.1385/MB:16:1:53 (2000).11098468 10.1385/MB:16:1:53
54. Liang H Kumar PA Nain V Powell WA Carlson JE Kole C Michler CH Abbott AG Hall TC Selection and screening strategies Principles and Development 2010 Springer-Verlag 85 143
Liang, H., Kumar, P. A., Nain, V., Powell, W. A. & Carlson, J. E. Selection and screening strategies. In Principles and Development (eds Kole, C. et al.) 85–143 (Springer-Verlag, 2010). 10.1007/978-3-642-04809-8_4.
55. De Buck S De Paepe A Depicker A Christou P Savin R Costa-Pierce BA Misztal I Whitelaw CBA Transgene expression in plants, control of Sustainable Food Production 2013 Springer 1570 1593
De Buck, S., De Paepe, A. & Depicker, A. Transgene expression in plants, control of. In Sustainable Food Production (eds Christou, P. et al.) 1570–1593 (Springer, 2013). 10.1007/978-1-4614-5797-8_412.
56. Jongedijk E Tigelaar H van Roekel JSC Bres-Vloemans SA Dekker I van den Elzen PJM Cornelissen BJC Melchers LS Synergistic activity of chitinases and β-1,3-glucanases enhances fungal resistance in transgenic tomato plants Euphytica 1995 85 173 180 10.1007/BF00023946
Jongedijk, E. et al. Synergistic activity of chitinases and β-1,3-glucanases enhances fungal resistance in transgenic tomato plants. Euphytica 85, 173–180 (1995).10.1007/BF00023946
57. Punja ZK Raharjo HT Response of transgenic cucumber and carrot plants expressing different chitinase enzymes to inoculation with fungal pathogens Plant Dis. 1996 80 999 1005 10.1094/PD-80-0999
Punja, Z. K. & Raharjo, H. T. Response of transgenic cucumber and carrot plants expressing different chitinase enzymes to inoculation with fungal pathogens. Plant Dis. 80, 999–1005 (1996).10.1094/PD-80-0999
58. Chen W Punja Z Transgenic herbicide- and disease-tolerant carrot (Daucus carota L.) plants obtained through Agrobacterium-mediated transformation Plant Cell Rep. 2002 20 929 935 10.1007/s00299-001-0419-7
Chen, W. & Punja, Z. Transgenic herbicide- and disease-tolerant carrot (Daucus carota L.) plants obtained through Agrobacterium-mediated transformation. Plant Cell Rep. 20, 929–935. 10.1007/s00299-001-0419-7 (2002).10.1007/s00299-001-0419-7
59. Fagoaga C Rodrigo I Conejero V Hinarejos C Tuset JJ Arnau J Pina JA Navarro L Peña L Increased tolerance to Phytophthora citrophthora in transgenic orange plants constitutively expressing a tomato pathogenesis related protein PR-5 Mol. Breed. 2001 7 175 185 10.1023/A:1011358005054
Fagoaga, C. et al. Increased tolerance to Phytophthora citrophthora in transgenic orange plants constitutively expressing a tomato pathogenesis related protein PR-5. Mol. Breed. 7, 175–185. 10.1023/A:1011358005054 (2001).10.1023/A:1011358005054
60. Pérez-Jiménez M Pérez-Tornero O In Vitro plant evaluation trial: Reliability test of salinity assays in Citrus plants Plants 2020 9 10 1352 10.3390/plants9101352 33066234
Pérez-Jiménez, M. & Pérez-Tornero, O. In Vitro plant evaluation trial: Reliability test of salinity assays in Citrus plants. Plants 9(10), 1352. 10.3390/plants9101352 (2020).33066234 10.3390/plants9101352
61. Pérez-Tornero O Tallón CI Porras I Navarro JM Physiological and growth changes in micropropagated Citrus macrophylla explants due to salinity J. Plant Physiol. 2009 166 1923 1933 10.1016/j.jplph.2009.06.009 19604601
Pérez-Tornero, O., Tallón, C. I., Porras, I. & Navarro, J. M. Physiological and growth changes in micropropagated Citrus macrophylla explants due to salinity. J. Plant Physiol. 166, 1923–1933. 10.1016/j.jplph.2009.06.009 (2009).19604601 10.1016/j.jplph.2009.06.009
62. Thakur M Sharma DR Sharma SK In vitro selection and regeneration of carnation (Dianthus caryophyllus L.) plants resistant to culture filtrate of Fusarium oxysporum f. sp. dianthi Plant Cell Rep. 2002 20 825 828 10.1007/s00299-001-0412-1
Thakur, M., Sharma, D. R. & Sharma, S. K. In vitro selection and regeneration of carnation (Dianthus caryophyllus L.) plants resistant to culture filtrate of Fusarium oxysporum f. sp. dianthi. Plant Cell Rep. 20, 825–828 (2002).10.1007/s00299-001-0412-1
63. Berrocal-Lobo M Molina A Ethylene response factor 1 mediates Arabidopsis resistance to the soilborne fungus Fusarium oxysporum Mol. Plant Microbe Interact. 2004 17 7 763 770 10.1094/MPMI.2004.17.7.763 15242170
Berrocal-Lobo, M. & Molina, A. Ethylene response factor 1 mediates Arabidopsis resistance to the soilborne fungus Fusarium oxysporum. Mol. Plant Microbe Interact. 17(7), 763–770 (2004).15242170 10.1094/MPMI.2004.17.7.763
64. Challagulla V Bhattarai S Midmore DJ In-vitro vs in-vivo inoculation: Screening for resistance of australian rice genotypes against blast fungus Rice Sci. 2015 22 132 137 10.1016/j.rsci.2015.05.017
Challagulla, V., Bhattarai, S. & Midmore, D. J. In-vitro vs in-vivo inoculation: Screening for resistance of australian rice genotypes against blast fungus. Rice Sci. 22, 132–137. 10.1016/j.rsci.2015.05.017 (2015).10.1016/j.rsci.2015.05.017
65. Del Castillo-González L Soudani S de la Cruz-Gómez N Manzanera JA Berrocal-Lobo M An improved method to study Phytophthora cinnamomi Rands zoospores interactions with host BMC Plant Biol. 2024 24 508 10.1186/s12870-024-05205-2 38844843
Del Castillo-González, L., Soudani, S., de la Cruz-Gómez, N., Manzanera, J. A. & Berrocal-Lobo, M. An improved method to study Phytophthora cinnamomi Rands zoospores interactions with host. BMC Plant Biol. 24, 508. 10.1186/s12870-024-05205-2 (2024).38844843 10.1186/s12870-024-05205-2
66. Chen D Zentmyer GA Production of sporangia by Phytophthora cinnamomi in axenic culture Mycologia 1970 62 397 402 10.1080/00275514.1970.12018977
Chen, D. & Zentmyer, G. A. Production of sporangia by Phytophthora cinnamomi in axenic culture. Mycologia 62, 397–402. 10.1080/00275514.1970.12018977 (1970).10.1080/00275514.1970.12018977
67. Byrt P Grant BR Some conditions governing zoospore production in axenic culturesof Phytophthora cinnamomi rands Aust. J. Bot. 1972 2 103 115 10.1071/BT9790103
Byrt, P. & Grant, B. R. Some conditions governing zoospore production in axenic culturesof Phytophthora cinnamomi rands. Aust. J. Bot. 2, 103–115. 10.1071/BT9790103 (1972).10.1071/BT9790103
68. Byrt PN Holland AA Infection of axenic Eucalyptus seedlings with Phytophthora cinnamomi zoo-spores Aust. J. Bot. 1978 2 169 176 10.1071/BT9780169
Byrt, P. N. & Holland, A. A. Infection of axenic Eucalyptus seedlings with Phytophthora cinnamomi zoo-spores. Aust. J. Bot. 2, 169–176. 10.1071/BT9780169 (1978).10.1071/BT9780169
69. Wang B Andargie M Fang R The function and biosynthesis of callose in high plants Heliyon 2022 8 4 e09248 10.1016/j.heliyon.2022.e09248 35399384
Wang, B., Andargie, M. & Fang, R. The function and biosynthesis of callose in high plants. Heliyon 8(4), e09248 (2022).35399384 10.1016/j.heliyon.2022.e09248
70. Van den Berg N Christie JB Aveling TAS Engelbrecht J Callose andb-1,3-glucanase inhibit Phytophthora cinnamomi in a resistant avocado rootstock Plant Pathol. 2018 67 1150 1160 10.1111/ppa.12819
Van den Berg, N., Christie, J. B., Aveling, T. A. S. & Engelbrecht, J. Callose andb-1,3-glucanase inhibit Phytophthora cinnamomi in a resistant avocado rootstock. Plant Pathol. 67, 1150–1160 (2018).10.1111/ppa.12819
71. Hinch JM Clarke AE Callose formation in Zea mays as a response to infection with Phytophthora cinnamomi Physiol. Plant Pathol. 1982 21 113 120 10.1016/0048-4059(82)90014-5
Hinch, J. M. & Clarke, A. E. Callose formation in Zea mays as a response to infection with Phytophthora cinnamomi. Physiol. Plant Pathol. 21, 113–120 (1982).10.1016/0048-4059(82)90014-5
72. Beffa RS Hofer R-M Thomas M Meins F Jr Decreased susceptibility to viral disease of [beta]-1,3-glucanase-deficient plants generated by antisense Transformation Plant Cell 1996 8 1001 1011 10.2307/3870211 12239410
Beffa, R. S., Hofer, R.-M., Thomas, M. & Meins, F. Jr. Decreased susceptibility to viral disease of [beta]-1,3-glucanase-deficient plants generated by antisense Transformation. Plant Cell 8, 1001–1011 (1996).12239410 10.2307/3870211
73. Tomczynska I Stumpe M Doan TG Mauch F A Phytophthora effector protein promotes symplastic cell-to-cell trafficking by physical interaction with plasmodesmata-localised callose synthases New Phytol. 2020 227 5 1467 1478 10.1111/nph.16653 32396661
Tomczynska, I., Stumpe, M., Doan, T. G. & Mauch, F. A Phytophthora effector protein promotes symplastic cell-to-cell trafficking by physical interaction with plasmodesmata-localised callose synthases. New Phytol. 227(5), 1467–1478 (2020).32396661 10.1111/nph.16653
74. Nishimura MT Stein M Hou BH Vogel JP Edwards H Somerville SC Loss of a callose synthase results in salicylic acid-dependent disease resistance Science 2003 301 969 972 10.1126/science.1086716 12920300
Nishimura, M. T. et al. Loss of a callose synthase results in salicylic acid-dependent disease resistance. Science 301, 969–972. 10.1126/science.1086716 (2003).12920300 10.1126/science.1086716
75. Cahill D Weste G Formation of callose deposits as a response to infection with Phytophthora cinnamomi Trans. Br. Mycol. Soc. 1983 80 23 29 10.1016/S0007-1536(83)80161-2
Cahill, D. & Weste, G. Formation of callose deposits as a response to infection with Phytophthora cinnamomi. Trans. Br. Mycol. Soc. 80, 23–29 (1983).10.1016/S0007-1536(83)80161-2
76. Pillai SE Kumar C Patel HK Sonti RV Overexpression of a cell wall damage induced transcription factor, OsWRKY42, leads to enhanced callose deposition and tolerance to salt stress but does not enhance tolerance to bacterial infection BMC Plant Biol. 2018 18 177 10.1186/s12870-018-1391-5 30176792
Pillai, S. E., Kumar, C., Patel, H. K. & Sonti, R. V. Overexpression of a cell wall damage induced transcription factor, OsWRKY42, leads to enhanced callose deposition and tolerance to salt stress but does not enhance tolerance to bacterial infection. BMC Plant Biol. 18, 177 (2018).30176792 10.1186/s12870-018-1391-5
77. Chen XY Kim JY Callose synthesis in higher plants Plant Signal Behav. 2009 4 489 492 10.4161/psb.4.6.8359 19816126
Chen, X. Y. & Kim, J. Y. Callose synthesis in higher plants. Plant Signal Behav. 4, 489–492. 10.4161/psb.4.6.8359 (2009).19816126 10.4161/psb.4.6.8359
78. Vasudevan SN Pooja SK Raju TJ Damini CS Cisgenics and intragenics: Boon or bane for crop improvement Front. Plant Sci. 2023 14 1275145 10.3389/fpls.2023.1275145 38089798
Vasudevan, S. N., Pooja, S. K., Raju, T. J. & Damini, C. S. Cisgenics and intragenics: Boon or bane for crop improvement. Front. Plant Sci. 14, 1275145. 10.3389/fpls.2023.1275145 (2023).38089798 10.3389/fpls.2023.1275145
79. Espinoza C Schlechter R Herrera D Torres E Serrano A Medina C Arce-Johnson P Cisgenesis and intragenesis: New tools for improving crops Biol. Res. 2013 46 4 323 331 10.4067/S0716-97602013000400003 24510134
Espinoza, C. et al. Cisgenesis and intragenesis: New tools for improving crops. Biol. Res. 46(4), 323–331. 10.4067/S0716-97602013000400003 (2013).24510134 10.4067/S0716-97602013000400003
80. EFSA Panel on Genetically Modified Organisms Scientific opinion addressing the safety assessment of plants developed through cisgenesis and intragenesis Eur. Food Saf. Auth. 2020 10 2561 33
EFSA Panel on Genetically Modified Organisms. Scientific opinion addressing the safety assessment of plants developed through cisgenesis and intragenesis. Eur. Food Saf. Auth. 10(2561), 33 (2020).
81. Dayé C Spök A Allan AC Yamaguchi T Sprink T Chaurasia A Kole C Social acceptability of cisgenic plants: Public perception, consumer preferences, and legal regulation Cisgenic crops: Safety, legal and social issues concepts and strategies in plant sciences 2023 Springer 43 75
Dayé, C., Spök, A., Allan, A. C., Yamaguchi, T. & Sprink, T. Social acceptability of cisgenic plants: Public perception, consumer preferences, and legal regulation. In Cisgenic crops: Safety, legal and social issues concepts and strategies in plant sciences (eds Chaurasia, A. & Kole, C.) 43–75 (Springer, 2023). 10.1007/978-3-031-10721-4_3.
82. Karimi M Inzé D Depicker A GATEWAY™ vectors for Agrobacterium-mediated plant transformation Trends Plant Sci. 2002 7 193 195 10.1016/s1360-1385(02)02251-3 11992820
Karimi, M., Inzé, D. & Depicker, A. GATEWAY™ vectors for Agrobacterium-mediated plant transformation. Trends Plant Sci. 7, 193–195. 10.1016/s1360-1385(02)02251-3 (2002).11992820 10.1016/s1360-1385(02)02251-3
83. Corredoira E Valladares S Vieitez AM Ballester A Wang K Chestnut, European (Castanea sativa) Methods in Molecular Biology Agrobacterium Protocols 2015 Springer- Verlag 163 176
Corredoira, E., Valladares, S., Vieitez, A. M. & Ballester, A. Chestnut, European (Castanea sativa). In Methods in Molecular Biology Agrobacterium Protocols (ed. Wang, K.) 163–176 (Springer- Verlag, 2015).
84. Hood EE Gelvin SB Melchers LS Hoekema A New Agrobacterium helper plasmids for gene transfer to plants Transgenic Res. 1993 2 208 218 10.1007/BF01977351
Hood, E. E., Gelvin, S. B., Melchers, L. S. & Hoekema, A. New Agrobacterium helper plasmids for gene transfer to plants. Transgenic Res. 2, 208–218. 10.1007/BF01977351 (1993).10.1007/BF01977351
85. Xu R Li QQ Protocol: Streamline cloning of genes into binary vectors in Agrobacterium via the Gateway-TOPO vector system Plant Methods 2008 4 4 10.1186/1746-4811-4-4 18211693
Xu, R. & Li, Q. Q. Protocol: Streamline cloning of genes into binary vectors in Agrobacterium via the Gateway-TOPO vector system. Plant Methods 4, 4. 10.1186/1746-4811-4-4 (2008).18211693 10.1186/1746-4811-4-4
86. Song P Cai C Skokut M Kosegi BD Petolino JF Quantitative real-time PCR as a screening tool for estimating transgene copy number in WHISKERS™-derived transgenic maize Plant Cell Rep. 2002 20 948 954 10.1007/s00299-001-0432-x
Song, P., Cai, C., Skokut, M., Kosegi, B. D. & Petolino, J. F. Quantitative real-time PCR as a screening tool for estimating transgene copy number in WHISKERS™-derived transgenic maize. Plant Cell Rep. 20, 948–954. 10.1007/s00299-001-0432-x (2002).10.1007/s00299-001-0432-x
87. Ramos A Usié A Barbosa P Barros PM Capote T The draft genome sequence of cork oak Sci. Data 2018 5 180069 10.1038/sdata.2018.69 29786699
Ramos, A. et al. The draft genome sequence of cork oak. Sci. Data 5, 180069. 10.1038/sdata.2018.69 (2018).29786699 10.1038/sdata.2018.69
88. Mallón R Vieitez AM Vidal N High-efficiency agrobacterium-mediated transformation in Quercus robur: Selection by use of a temporary immersion system and assessment by quantitative PCR Plant Cell Tiss Organ Cult. 2013 114 171 185 10.1007/s11240-013-0313-3
Mallón, R., Vieitez, A. M. & Vidal, N. High-efficiency agrobacterium-mediated transformation in Quercus robur: Selection by use of a temporary immersion system and assessment by quantitative PCR. Plant Cell Tiss Organ Cult. 114, 171–185. 10.1007/s11240-013-0313-3 (2013).10.1007/s11240-013-0313-3
89. Chang S Puryear J Cairney J A simple and efficient method for isolating RNA from pine trees Plant Mol. Biol. Rep. 1993 11 113 116 10.1007/BF02670468
Chang, S., Puryear, J. & Cairney, J. A simple and efficient method for isolating RNA from pine trees. Plant Mol. Biol. Rep. 11, 113–116. 10.1007/BF02670468 (1993).10.1007/BF02670468
90. Gambino G Perrone I Gribaudo I A Rapid and effective method for RNA extraction from different tissues of grapevine and other woody plants Phytochem. Anal. 2008 19 6 520 525 10.1002/pca.1078 18618437
Gambino, G., Perrone, I. & Gribaudo, I. A Rapid and effective method for RNA extraction from different tissues of grapevine and other woody plants. Phytochem. Anal. 19(6), 520–525. 10.1002/pca.1078 (2008).18618437 10.1002/pca.1078
91. Kotrade P Sehr EM Wischnitzki E Brüggemann W Comparative transcriptomics-based selection of suitable reference genes for normalization of RT-qPCR experiments in drought-stressed leaves of three European Quercus species Tree Genet. Genomes 2019 15 3 38 10.1007/s11295-019-1347-4
Kotrade, P., Sehr, E. M., Wischnitzki, E. & Brüggemann, W. Comparative transcriptomics-based selection of suitable reference genes for normalization of RT-qPCR experiments in drought-stressed leaves of three European Quercus species. Tree Genet. Genomes 15(3), 38. 10.1007/s11295-019-1347-4 (2019).10.1007/s11295-019-1347-4
92. Marum L Miguel A Ricardo CP Miguel C Reference gene selection for quantitative real-time PCR normalization in Quercus suber PLoS One 2012 7 4 e35113 10.1371/journal.pone.0035113 22529976
Marum, L., Miguel, A., Ricardo, C. P. & Miguel, C. Reference gene selection for quantitative real-time PCR normalization in Quercus suber. PLoS One 7(4), e35113. 10.1371/journal.pone.0035113 (2012).22529976 10.1371/journal.pone.0035113
93. Pfaffl MW A new mathematical model for relative quantification in real-time RT-PCR Nucleic Acids Res. 2001 29 9 45e 445 10.1093/nar/29.9.e45
Pfaffl, M. W. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 29(9), 45e–445. 10.1093/nar/29.9.e45 (2001).10.1093/nar/29.9.e45
94. Toribio M Celestino C Molinas M Jain SM Gupta PK Cork oak, Quercus suber L. Protocol for Somatic Embryogenesis in Woody Plants 2005 Springer 445 457
Toribio, M., Celestino, C. & Molinas, M. Cork oak, Quercus suber L. In Protocol for Somatic Embryogenesis in Woody Plants (eds Jain, S. M. & Gupta, P. K.) 445–457 (Springer, 2005).
95. OEPP/EPPO Bulletin. EPPO Standards Diagnostic protocols for regulated pests. 34, 155–157 (2004).
96. Soudani S Poza-Carrión C De la Cruz Gómez N González-Coloma A Andrés MF Berrocal-Lobo M Essential oils prime epigenetic and metabolomic changes in tomato defense against Fusarium oxysporum Front. Plant Sci. 2022 13 804104 10.3389/fpls.2022.80410 35422834
Soudani, S. et al. Essential oils prime epigenetic and metabolomic changes in tomato defense against Fusarium oxysporum. Front. Plant Sci. 13, 804104. 10.3389/fpls.2022.80410 (2022).35422834 10.3389/fpls.2022.80410
97. Fernández-Bautista N Domínguez-Núñez JA Moreno MMC Berrocal-Lobo M Plant tissue trypan blue staining during phytopathogen infection Bio-protocol 2016 6 24 e2078 e2078 10.21769/BioProtoc.2078
Fernández-Bautista, N., Domínguez-Núñez, J. A., Moreno, M. M. C. & Berrocal-Lobo, M. Plant tissue trypan blue staining during phytopathogen infection. Bio-protocol 6(24), e2078–e2078 (2016).10.21769/BioProtoc.2078
98. Chen JC Fang SC The long pollen tube journey and in vitro pollen germination of Phalaenopsis orchids Plant Reprod. 2016 29 179 188 10.1007/s00497-016-0280-z 27016359
Chen, J. C. & Fang, S. C. The long pollen tube journey and in vitro pollen germination of Phalaenopsis orchids. Plant Reprod. 29, 179–188 (2016).27016359 10.1007/s00497-016-0280-z
