Next Article in Journal
Common Marmoset Cell Lines and Their Applications in Biomedical Research
Next Article in Special Issue
Genome-Wide Identification and Characterization of CDPK Gene Family in Cultivated Peanut (Arachis hypogaea L.) Reveal Their Potential Roles in Response to Ca Deficiency
Previous Article in Journal
Hypometabolism, Alzheimer’s Disease, and Possible Therapeutic Targets: An Overview
Previous Article in Special Issue
Structure of Chlorella ohadii Photosystem II Reveals Protective Mechanisms against Environmental Stress
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Biotechnological Potential of the Stress Response and Plant Cell Death Regulators Proteins in the Biofuel Industry

by
Maciej Jerzy Bernacki
1,2,†,
Jakub Mielecki
1,†,
Andrzej Antczak
3,†,
Michał Drożdżek
3,†,
Damian Witoń
1,
Joanna Dąbrowska-Bronk
4,
Piotr Gawroński
1,
Paweł Burdiak
1,
Monika Marchwicka
3,
Anna Rusaczonek
5,
Katarzyna Dąbkowska-Susfał
6,
Wacław Roman Strobel
2,
Ewa J. Mellerowicz
7,
Janusz Zawadzki
3,
Magdalena Szechyńska-Hebda
8 and
Stanisław Karpiński
1,*
1
Department of Plant Genetics, Breeding and Biotechnology, Institute of Biology, Warsaw University of Life Sciences, Nowoursynowska Street 159, 02-776 Warsaw, Poland
2
Institute of Technology and Life Sciences—National Research Institute, Falenty, Al. Hrabska 3, 05-090 Raszyn, Poland
3
Institute of Wood Sciences and Furniture, Warsaw University of Life Sciences—SGGW, 02-776 Warsaw, Poland
4
Department of Plant Physiology, Institute of Biology, Warsaw University of Life Sciences, 02-776 Warsaw, Poland
5
Department of Botany, Institute of Biology, Warsaw University of Life Sciences, 02-776 Warsaw, Poland
6
Faculty of Chemical and Process Engineering, Warsaw University of Technology, 00-645 Warsaw, Poland
7
Department of Forest Genetics and Plant Physiology, Umeå Plant Science Centre, Swedish University of Agricultural Sciences, 901-83 Umeå, Sweden
8
W. Szafer Institute of Botany Polish Academy of Sciences, Lubicz 46, 31-512 Kraków, Poland
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Cells 2023, 12(16), 2018; https://doi.org/10.3390/cells12162018
Submission received: 3 July 2023 / Revised: 26 July 2023 / Accepted: 1 August 2023 / Published: 8 August 2023

Abstract

:
Production of biofuel from lignocellulosic biomass is relatively low due to the limited knowledge about natural cell wall loosening and cellulolytic processes in plants. Industrial separation of cellulose fiber mass from lignin, its saccharification and alcoholic fermentation is still cost-ineffective and environmentally unfriendly. Assuming that the green transformation is inevitable and that new sources of raw materials for biofuels are needed, we decided to study cell death—a natural process occurring in plants in the context of reducing the recalcitrance of lignocellulose for the production of second-generation bioethanol. “Members of the enzyme families responsible for lysigenous aerenchyma formation were identified during the root hypoxia stress in Arabidopsis thaliana cell death mutants. The cell death regulatory genes, LESION SIMULATING DISEASE 1 (LSD1), PHYTOALEXIN DEFICIENT 4 (PAD4) and ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1) conditionally regulate the cell wall when suppressed in transgenic aspen. During four years of growth in the field, the following effects were observed: lignin content was reduced, the cellulose fiber polymerization degree increased and the growth itself was unaffected. The wood of transgenic trees was more efficient as a substrate for saccharification, alcoholic fermentation and bioethanol production. The presented results may trigger the development of novel biotechnologies in the biofuel industry.

1. Introduction

One of the most important challenges in the 21st century is to provide reliable, renewable and environmentally friendly sources of energy. Moreover, the current geopolitical situation forces the invention of an alternative solution to traditional energy sources even more urgently. One of the answers is improvement of renewable biofuels production. Currently, bioethanol is still mostly produced as a first-generation biofuel from edible plants. As a result, it reduces food sources and thus is questionable. Although lignocellulosic biomass is widely obtained from tree plantations or agricultural waste, only ~3% of bioethanol is produced from this source as second-generation bioethanol [1,2,3,4]. The main limitation of efficient, industrial production of second-generation bioethanol is its complex composition and the structure of plant cell walls. Mostly composed of cellulose, the cell wall is highly heterogeneous as it also contains hemicelluloses, pectin and lignin. The outcome structure is determined by the function of secondary cell walls [5,6,7]. During technological processing, lignocellulosic biomass requires the separation of its components and breaking down complex polysaccharides into fermentable sugars [8,9]. In particular, lignin presents a great challenge as it affects the activity of cellulases [9,10] and also because phenolic compounds are inhibitors of alcohol fermentation in brewers’ yeast (Saccharomyces cerevisiae) [11]. Pretreatments to reduce recalcitrance increase the cost of production and negatively impact the environment with toxic wastes [12].
Plants are capable of efficient saccharification in vivo during induction of certain types of Programmed Cell Death (PCD) [13] e.g., during lysigenous aerenchyma formation in response to root hypoxia [14]. Discovered proteins witch regulated lysigenous aerenchyma formation opened a new pathway to study this form of plants adaptation to stress. We also postulated that the ability of plants to degradation of their cell wall can be an interesting target to improve second-generation bioethanol production. According to our knowledge, plant enzymes involved in the natural saccharification of the cell wall have not been explored in the context of their use as a potential source for industrial enzymes. Similarly, the role of genes encoding molecular PCD regulators has not been studied in the cell wall lignification and cellulose fiber polymerization-related processes.
For this study, we used two model plants in biotechnology, well-established Arabidopsis thaliana and Populus tremula x tremuloides which have been proposed as the best tree for genetic research and considered as Arabidopsis for forestry [15]. Although Arabidopsis thaliana is a herbaceous and annual plant, it shows the ability to form secondary xylem. It was shown that when Arabidopsis thaliana plants grown in appropriate conditions exhibit secondary growth in the hypocotyl, with the development of both a vascular and a cork cambium. Arabidopsis thaliana was previously proposed as a model for wood formation in trees [16]. We believe that our approach is valuable because it shows that the role of the tested proteins is strongly conserved in the plant kingdom, which makes it possible to transfer the knowledge obtained to other crop species.
Here we present the biotechnological potential of the conditional PCD regulators [17,18,19] i.e., LESION SIMULATING DISEASE 1 (LSD1), ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1), and PHYTOALEXIN DEFFICIENT 4 (PAD4). The very regulators were earlier recognized as involved in the lysigenous aerenchyma formation in response to root hypoxia stress in Arabidopsis thaliana [14]. Conditional PCD regulation by LSD1, EDS1 and PAD4 was demonstrated in runaway PCD during the growth of lsd1 mutant in ambient laboratory conditions and in field conditions [17,18], which affected its seed yield and productivity. Based on the results in stable transgenic aspen lines with deregulated PtEDS1, PtLSD1, and PtPAD4, we obtained four independent transgenic lines with lower lignin content and a higher cellulose polymerization degree in four-year-old aspen wood. Higher efficiency of the industrial bioethanol production from transgenic wood with the lowest lignin content was found.

2. Materials and Methods

2.1. Arabidopsis thaliana Plant Material

Arabidopsis thaliana seeds used in this study were available in our lab. Experiments were performed using lsd1, eds1, lsd1/eds1 mutants in two ecotypes, Col-0 and Ws-0.

2.2. Populus tremula × tremuloides Plant Material and Genetic Transformation

The silencing constructs for LSD1/PAD4 (XM_002317030, XM_002310499) and EDS1/PAD4 (XM_002318590, XM_002310499) genes were carried out based on the RNA interference using a binary vector pH7GWIWG2(I) under the control of the 35S promoter [20] (Figure S1). Short cDNA sequences of target genes were amplified using primers described in the supplementary data (Figure S2). Inserts of analyzed genes were introduced into the pH7GWIWG2(I) vector according to the ‘Gateway’ technology [20] and transferred to Agrobacterium tumefaciens C58. Next, some bacteria were used in the transformation of hybrid aspen plants. Transgenic hybrid aspen lines (P. tremula × tremuloides, background T89) were generated at the Umeå Plant Science Center according to the transformation procedure [21] from among 19 independent transgenic P. tremula × tremuloides lines, 4 of them with significant deregulation in expression levels of target genes were chosen for further experiments. For it is known that in Arabidopsis thaliana LSD1, EDS1 and PAD4 are co-regulated, we checked the expression level of PtLSD1, PtEDS1 and PtPAD4 in aspen leaves from field conditions. Based on these results we described the following lines: Lines 1, 2, 3 and 4.

2.3. Aspen Growth Condition

Transgenic aspen lines were propagated under in vitro conditions. Plants were grown in phytotron during a long photoperiod (16 h/8 h day/night) with a light intensity of 120 μmol photons m−2 s−1, temperature 25 °C/20 °C day/night and relative humidity 55%. Every 2 weeks, plants were transferred to a fresh Murashige and Skoog (MS) medium [22]. Well-rooted plants were transferred to the greenhouse, kept 3 days in open jars to acclimatize and then transferred to pots with soil. Transgenic and controlled aspen plants were grown in the greenhouse under natural light supported by low-pressure sodium lamps (Philips, Amsterdam, Netherlands) emitting light intensity (300 μmol photons m−2 s−1). After 9 months the transgenic and control lines were transplanted to a field and grown in the Wolica experimental field (52°8′30″ N, 21°4′12″ E) in 2017. Experiments and measurements were performed on 4-year-old transgenic and wild-type aspen plants. Trees were randomized in the experimental field to minimize the effect of placement on observed phenotype. The scheme of the field experiment is shown in the figure (Figure S3). We performed analysis on all four independent transgenic lines (with a minimum of six trees per Line) obtained from independent transformation events as we wanted to be sure that the phenotype that we observe is a consequence of targeted genetic manipulation and not a side effect of transgene integration into aspen genome.

2.4. Inducing of Lysigenous Aerenchyma Formation by Root Hypoxia Stress

In general, the experiment was conducted as described [14] above. Arabidopsis thaliana was grown in short-day photoperiod conditions (8 h of day and 16 h of night) for 8 weeks. After that, they were submerged in degassed water. The oxygen content in the water was reduced by boiling. After it cooled down, we checked the O2 content in water using an oxygen sensor DO210 (Extech Instruments, Nashua, NH, USA) before submerging the plants. An eventual evaporation was mitigated by topping up the water. After 7 days of hypoxia stress, hypocotyls and main roots material were harvested for future analysis (Figure S4A,B).

2.5. Determining to Reduce Sugar Content in Arabidopsis thaliana Hypocotyls

Reducing sugar content was determined using the Amplex Red Glucose Assay kit (Thermo Fisher Scientific, Waltham, MA, USA) accordingly to the manufacturer’s protocol.

2.6. Cellulase and Xylanase Activity Assays

Proteins from hypocotyls were extracted using a plastic homogenizer pestle and 1.5 mL eppendorf type tube in the PBS buffer on the ice at 4 °C. Protein content in samples was measured using the Bradford method. For activity analysis we used 20 μg of total protein. Cellulase and xylanase activities were measured against CMC (Sigma Aldrich, St. Louis, MO, USA) and beechwood xylan (Roth, Karlsruhe, Deutschland) as substrates, both at 1% (W:V) concentration. We measured the reducing sugar content in each sample before adding substrate and normalized results so that they could reflect how much the reducing sugars were produced from the substrate during the assay and not during the stress itself. We measured reducing sugar content after every hour, enzymatic reactions were carried out up to 3 h.

2.7. Native-PAGE Separation and Zymographic Activity Assay

We separated 20 ug of total proteins from each sample on 8% polyacrylamide gel. As we wanted to visualize the enzymatic activity, we did not use any reducing or denaturizing agents such as SDS, β-mercaptoethanol or DTT. Gels were supplemented with 1% CMC as substrate. After separation (2 h, 60–70 V) at 4 °C each gel was submerged in PBS and incubated at room temperature for 18 h. The gels were stained with Congo Red (Merck, Darmstadt, Germany) which binds to glycosidic bonds in cellulose. Fifteen min later the gels were detained using 2 M NaCl until clear zones were observed on the gel which indicate activity of cellulases (Figure S4C).

2.8. Maldi-TOF Analysis of Isolated Proteins

Bands exhibiting cellulase activity were cut off from each sample and sent to analysis on Maldi-TOF to The Environmental Mass Spectrometry Laboratory in the Department of Biophysics of the IBB PAN (Warsaw, Poland). Samples were treated with trypsin before separating and analyzing them on a mass spectrometer. The results were analyzed using the Mascot software and trimmed so there would be no more than 1% of false positives in each sample. The equipment used was sponsored in part by the Centre for Preclinical Research and Technology (CePT), a project co-sponsored by the European Regional Development Fund and Innovative Economy, The National Cohesion Strategy of Poland.

2.9. RNA Isolation and qPCR Analysis

RNA isolation, cDNA synthesis and qPCR analysis were performed as described before [23]. UBQ (LP-GCTTGAAGATGGGAGGACTCT, RP-AAATCTGCATTCCTCCACGG), PP2A (LP-ACTCTCTGCACTGTTGAGGAA, RP-ACCAATCAACCAGGTCAGTCT) and ACT (LP-TTACATGTTCACCACCACTGC RP-CTGCTCATAGTCAAGGGCAAC) were used as reference genes. For checking LSD1, EDS1 and PAD4 the following primers were used: LSD1 (LP-GACTCAAACTGTCGTTGTGGA, RP-CTGTAGTGACACCAACAACAAC, EDS1 (LP-CGACACCTCAAGAATGAAGA, RP-GTGCATTTATACCGTCTTGG), PAD4 (LP-CATCACCTTTGGCTCTCCAT, RP-AAATTTGCCACCCCATCTTT).

2.10. RNAseq Analysis

Fastq files were trimmed using TrimGalore 0.6.6 (https://github.com/FelixKrueger/TrimGalore (accessed on 20 June 2021)) with the following settings: --paired -j 4 --fastqc --quality 30 --length 50 --clip_R1 10 --clipR2 10 --three_prime_clip_R1 10 --three_prime_clip_R2 10. The quality of trimmed reads was assessed using the FastQC version (https://github.com/s-andrews/FastQC (accessed on 20 June 2021)). Reads were mapped to the Populus trichocarpa v3.0 genome [24] with annotations in version 4.1 downloaded from JGI (https://data.jgi.doe.gov/ (accessed on 20 June 2021)) using STAR aligner 2.7.6a [25] with the following settings: --outFilterScoreMinOverLread 0.3 --outFilterMatchNminOverLread 0.3 --outFilterMultimapNmax 100 --outReadsUnmapped Fastx --alignIntronMax 11,000 --quantMode GeneCounts. After mapping reads were counted using HTseq with the settings: -f bam -r pos -s no. Differentially expressed genes were identified with DEseq2 [26].

2.11. Biometric Parameters and Gas Exchanging Measurement

Trees were weighed and measured using an industrial scale and a commercially available measure. CO2 measurement was performed using LI-6400XT Portable Photosynthesis System (Li-COR Inc., Lincoln, NE, USA) just as described before [23].

2.12. Determination of Photosynthetic Pigments Content

The photosynthetic pigments measurements were determined as precisely described in the previous chapters [27]. Pigments extracted from the frozen tissue (grinded in liquid nitrogen before and stored at −80 °C) were separated on a SynergiTM 4 μm Max-RP 80Å 250 × 4.6 mm column (Phenomenex, Torrance, CA, USA) at 30 °C using Shimadzu HPLC System (Shimadzu, Kyoto, Japan). The results were expressed as peak area per μg of fresh weight.

2.13. Gene Ontology

Gene ontology analyses were performed using Panther (http://www.pantherdb.org/ (5 August 2022).

2.14. Cellulose, Hemicellulose and Lignin Content

Dry wood material was milled using a Retsch mill (SM 200, Retsch, Germany) and sieved (AS 200, Retsch, Germany). Wood powder with particle sizes from 0.43 mm to 1.02 mm was used for chemical analyses. The material was extracted in Soxhlet apparatus with a mixture of chloroform-ethanol 93:7 w/w for 10 h before analyses. Cellulose content was determined by the Kürschner–Hoffer method and holocellulose content was analyzed according to Wise [28]. Mineral substances content was determined by [29] on dust fraction (under 0.43 mm) and lignin content was determined according to TAPPI UM 250 (1985) and TAPPI T222 om-02 (2006).

2.15. The Alkaline Pre-Treatment and Enzymatic Hydrolysis of Aspen Wood

For the alkaline pre-treatment, a 2% aqueous NaOH solution was used in the proportions of 100 cm3 of sodium hydroxide per 5 g of dry aspen wood powder from the fraction of 0.43 mm to 1.02 mm. The alkaline pre-treatment was performed in an autoclave under the following conditions: 121 °C, p = 0.1 MPa, 30 min. Afterward the activated material was subjected to enzymatic hydrolysis according to the procedure described by us previously [30] with some changes specified below. The hydrolysis Cellic CTec2 enzyme (Novozymes, Krogshoejvej, Denmark) was used and activity was 148 FPU/cm3 determined by the NREL method [31]. The enzymatic hydrolysis process was carried out in the Erlenmeyer flasks and the total volume of the mixture was 100 cm3. In the process 50 cm3 of 0.1 M citrate buffer solution at pH = 4. was used and due to the use of yeast in the next step, no sodium azide was added during the hydrolysis. In the hydrolysis process, 3.325 cm3 of 25% (v/v) Cellic CTec2 enzyme solution on each sample was used in the ratio of 1 g of concentrated enzyme to 1 g of absolutely dry biomass. The hydrolysis proceeded in an incubator IKA KS 3000i control (IKA, Warsaw, Poland) with shaking of 150 rpm at 50 °C. The total time of enzymatic hydrolysis was 72 h. After the very process, glucose and xylose concentrations in the supernatant were analyzed by the HPLC method.

2.16. The Ethanol Fermentation Process

Saccharomyces cerevisiae KKP 50 yeast strain on YPD agar plate obtained from the Institute of Agricultural and Food Biotechnology culture collection (Warsaw, Poland) was used in the fermentation study. Yeast inoculum was prepared by transferring single colonies from the YPD agar plate using a sterile loop into 25 cm3 of autoclaved (15 min at 121 °C) liquid YPD medium (Millipore, Burlington, MA, USA) in 100 cm3 Erlenmeyer flask and incubating flask overnight at 30 °C with shaking of 240 rpm. Fermentation broths were prepared as follows: 20 cm3 of each enzymatic hydrolysate in 100 cm3 Erlenmeyer flasks were supplemented with 0.2 g of yeast extract (Biocorp, Issoire, France) and autoclaved for 15 min at 121 °C. Subsequently, 1 cm3 of inoculum cultured overnight was added into each flask. Batch cultures were carried out in triplicate at 30 °C with shaking of 240 rpm for 47.5 h. During cultivation 0.5 cm3 samples were taken from fermentation broths to 1.5 cm3 Eppendorf tubes, cooled in ice water and then centrifuged for 10 min.) to remove the yeast biomass. The obtained supernatants were frozen until the HPLC analysis.

2.17. The HPLC Analysis

The concentration of glucose and xylose in the enzymatic hydrolysates (before fermentation) or glucose, xylose and ethanol in the fermentation broths (after fermentation) were determined using the high-performance liquid chromatography (HPLC, Varian 635 CL System, USA) with a refractive index (RI) detector (Smartline 2300, Knauer, Germany) and thermostated at 60 °C Rezex ROA Organic acid H+ (8%) (300 × 7.8 mm) column with Security Guard Cartridge Carbo-H (4 × 3.0 mm) (Phenomenex, Torrance, CA, USA) using 0.001 N H2SO4 as the mobile phase at a flow rate of 0.4 cm3 min−1. All the HPLC analyses were executed in triplicate. The concentrations were expressed in g * dm−3 of hydrolysates or fermentation broths.

2.18. The Analysis of Cellulose Polymerization Degree

The degree of cellulose polymerization was determined using the SEC (Size Exclusion Chromatography) technique. Cellulose was isolated from aspen wood with the Kürschner-Hoffer method, dissolved in 8% lithium chloride/N, N-dimethylacetamide system and analyzed as described in detail in the previous publications [8,15]. All of the SEC analyses were executed in triplicate or six times.

2.19. Statistic

For the aspen wood properties, biometric measurement and Arabidopsis thaliana the statistical analysis was performed using the Tukey honest significant difference (HSD) test. The number of biological and technical repetitions can be found in the description of the Figures. The qPCR data statistical analysis was performed using LinReg.

3. Results

3.1. EDS1 and LSD1 Are Involved in Lysigenous Aerenchyma Formation

It was confirmed that LSD1, EDS1 and PAD4 are involved in the lysigenous aerenchyma formation [14] (Figure 1 and Figure S5A–D). Changes in the content of reducing sugars in roots and hypocotyls were determined in plants growing in controlled conditions and in plants flooded for 7 days (root hypoxia stress, see material and methods) to establish a cell wall saccharification and lysis occurring during the formation of lysigenous aerenchyma in hypocotyls. The elevated level of reducing sugars before hypoxia was only observed in the lsd1 mutant when compared to the wild-type (Ws-0) plant (Figure 1A), which is related to the fact that this mutant induces cell death in long day conditions without the action of external stress factors (Figure S4B) [32]. After the stress, the lsd1 mutant exhibited an elevated (compared to Ws-0) but stable (as compared to lsd1 plants before hypoxia) level of reducing sugars. It is related to the fact that lsd1 exhibit enhanced cell death levels in long-day conditions in comparison to the wild type [32]. In contrast, the eds1 mutant exhibited the highest level of reducing sugars, whereas the eds1/lsd1 double mutant did not differ from Ws-0 after 7 days of stress (Figure 1B) Moreover we found no difference between WT and double lsd1/eds1 mutant in both conditions (Figure 1A,B). This is caused by the actions of these two proteins which are closely related and LSD1 is antagonistic to EDS1. It was many times showed that phenotype of lsd1/eds1 is similar to the WT phenotype, for example in response to various types of biotic and abiotic stresses [17,18]. These results suggest that both LSD1 and EDS1 are regulators of saccharification and lysis occurring during the formation of lysigenous aerenchyma in hypocotyls. However, LSD1 seems to be a regulator itself, while EDS1 acts as a regulator together with PAD4, and the regulatory role of EDS1-PAD4 [33] in this process is important under stress conditions.

3.2. In Vitro Cellulolytic and Xylanolytic Activities of Enzymes Isolated from eds1 Hypocotyls

Two different substrates i.e., carboxymethyl cellulose (CMC) and beechwood xylan were added to proteins extracted from Ws-0 and eds1 hypocotyls. CMC is a substrate for β-D-glucanase hydrolyzing β-1-4 glycosidic linkages between glucopyranose units (in amorphous cellulose), while xylan is a substrate for xylanase acting on β-1-4 glycosidic linkages between xylopyranose units in xylan backbone. The Arabidopsis thaliana enzymes induced in roots and hypocotyls during hypoxia stress have the potential to break down both cellulose and xylan, however, proteins extracted from the eds1 mutant exhibited higher activity than those extracted from Ws-0 (Figure 1C,D). In contrast, enzymatic extracts from lsd1 and lsd1/eds1 double mutant exhibited lower activity against both CMC and xylan than eds1 mutant and wild-type plants (Figure S5A–D). The high reducing sugar content after the stress and stable enzymatic activity of extracts of eds1 mutant indicated that cellulase, xylanases and also cell wall loosening (CWL) proteins were intact and unaffected probably by endogenous proteolysis. Similarly, higher activity was observed in the protein extracts from eds1 in the Col-0 background examined in the zymography activity assay (Figure S4C; proteins were separated with Native-PAGE in the presence of 1% CMC as substrate). The ecotype Col-0 showed lower activity than Ws-0 (Figure S4C) and, in the same way, it was observed in the eds1 mutant. Generally, proteins with molecular weight higher than 150 kDa had enzymatic activity.

3.3. Enzymes Generating Reducing Ends Were Found in Hypocotyls and Roots Undergoing Lysigenous Aerenchyma Formation

In the protein fraction with enzymatic activity that was isolated from the zymographic activity assay, many proteins in individual genotypes were identified during the Maldi-TOF analysis (Table S1). Some of them were found in two or more genotypes, but after removing duplicates, we found a list of 563 proteins (Table S2). Approximately 6.4% (37 proteins) of all identified proteins were assigned/aligned to known or predicted proteins capable of generating reducing sugars (Figure 2A,B and Table S4). We identified enzymes from different CAZY families, with known or predicted CWL activities (homogalacturonan, cellulose, xylan, mannan, callose, xyloglucan, and hydroxyproline-rich or arabinogalactan, Tables S3 and S4). The protein ontology analysis (Figure 2C and Figure S6) using Panther (Pantherdb.org) software was performed for eds1 mutants in both ecotypes, due to exhibited saccharification potential. We found four proteins whose abundance depends on EDS1 dysfunction (gene expression inhibition) (Figure 2D). The very proteins likely representing cell wall dismantling machinery could become interesting targets as plant enzymes for future use in plant biotechnology. Most importantly, we found proteins involved in the polysaccharide catabolite processes and proteins involved in the lignin precursor catabolic processes (Table S5). However, many (120 out of 563) of identified proteins, which could be valuable in biofuels or pulp industry, did not annotate protein classes or other gene ontology classes in the Panther Arabidopsis thaliana database (Table S6). These results confirm that EDS1 is a negative regulator of CWL and saccharification during the formation of lysigenous aerenchyma.

3.4. Cell Walls of Aspen with Deregulated PCD Growing in the Field Contain Less Lignin and Longer Cellulose Fibers

It was demonstrated earlier that transgenic aspen with an individually lowered expression of PtLSD1, PtEDS1, or PtPAD4, differentially modified wood density, wood swelling and cell wall thickness [34,35] in comparison to the wild type in the field conditions. Here four transgenic lines with reduced expression of PtLSD1, PtEDS1, and PtPAD4 (lines: Line 1, Line 2, Line 3 and Line 4) (Figure 3A–C) were tested. All three genes, i.e., PtLSD1, PtEDS1, and PtPAD4 were deregulated regardless of the double transgene combination (even if the third gene was not modified in a double transgenic line), probably because these genes are conditionally co-regulated, while LSD1 and PAD4 interact with EDS1 protein [17,18,19,34,35] (Figure 3A–C). Transgenic trees exhibited significantly lowered lignin content in the wood and higher cellulose fiber polymerization degree (Figure 3D). This phenotype (feature) correlated with the expression level of LSD1, EDS1, and PAD4 (Figure S7). The relative reduction of lignin (Wt = 100%) was in the range of 10% in Line 1 to 22% in Line 4 (Figure 3D). The cellulose content was also significantly lower (up to 7% in Line 4), however, its polymerization degree was significantly increased to 23% (Figure 3D).

3.5. Strong Reduction in Lignin and Higher Cellulose Polymerization Degree in the Wood of Line 4 Allows High Fermentation Efficiency

Lignin is an inhibitor of alcoholic fermentation [36]. Therefore, we confirm that Line 4 was the most suitable for bioethanol production. Higher glucose yield after alkaline pre-treatment and enzymatic hydrolysis of aspen wood (wild-type plants and Line 4) (Figure 3E) as well as higher ethanol yield after fermentation (Figure 3F) were obtained. The high effectiveness of the alcohol fermentation process was confirmed, as almost all the glucose obtained from enzymatic hydrolysis was converted to ethanol for Line 4 wood but not for the wild type (Figure 3G).

3.6. PtLSD1, PtEDS1, and PtPAD4 Strongly Influence Cambium and Xylem Transcriptome in Line 4

Having analyzed the transcriptome of cambium and differentiating xylem tissue in Line 4 we found that it differed from the wild-type tree (Figure S8) with 2648 differentially expressed genes (DEG) (Table S7) of which 2299 were known (recognized by Panther.org, Table S8) and 342 were unknown (Table S9). For all the known genes, gene ontology analysis was performed considering their molecular function (Table S10 and Figure S9) or biological process (Table S11 and Figure S10). Many differently expressed genes encoded with putative enzymes (Figure S10). The overrepresentation of genes assigned to the regulation of the phenylpropanoid metabolic process, which is related to lignin biosynthesis [37], was found (Table S11). Some of the 63 genes (ca. 3.2%) were involved in lignin and cellulose synthesis/metabolism, cell wall biogenesis and cell wall organization (Table S12). Moreover, we compared genes encoding enzymes identified in the Maldi-TOF experiment in Arabidopsis thaliana (Table S2) with their homologs among genes deregulated in transgenic Line 4. For 19 Line 4 DEGs, its Arabidopsis homologs were identified in the Maldi-TOF experiment (Table S13). This list included two cell wall-loosening proteins and enzymes generating reducing sugars; ALPHA-L-FUCOSIDASE 2 and BETA-D-XYLOSIDASE 7, which could be interesting targets for the biofuel industry. In order to find some gene clusters, which can be particularly important in the context of aspen wood quality, we compared Line 4 DEGs with the known CAZY family protein database (Table S14) [38] to genes involved in the cell wall integrity (malectin and malectin-like domain-containing proteins (Table S15) [39]. In Line 4 DEGs 252 we found genes encoding proteins from CAZY families and 23 genes encoding malectin and malectin-like domain-containing proteins. 14 genes from Line 4 DEGs were found in both groups (Table S16 and Figure S11). These results strongly indicate that conditional PCD regulators LSD1, EDS1 and PAD4 are developmental regulators of the secondary cell wall structure in the woody plants.

3.7. Transgenic Lines Exhibited Similar Phenotype and Development Level as Wild Type Trees

Insignificant differences in tree height, stem fresh weight and stem diameter during annual growth were measured. In contrast, CO2 assimilation was changed in transgenic lines, even if differences were not found in chlorophyll content (Figure S12).

4. Discussion

Cell wall biotechnological amelioration with decreased lignin content in the wood a goal of many studies [3,40]. However, the knowledge about cell wall structure and composition regulators, plant cellulolytic and CWL enzymes naturally decomposing cell walls is limited [41,42]. Wood tissue and cell wall development require the occurrence of PCD [43]. Therefore, conditional regulators of PCD and lysigenous aerenchyma formation i.e., LSD1, EDS1, and PAD4 [17,18,19,44,45] were studied in both Arabidopsis thalian and aspen for their biotechnological potential in biofuels production.
Based on the results of this study, we report that the formation of lisogenous aerenchyma in response to root hypoxia stress is regulated by LSD1 and EDS1 and involves at least 37 CWL enzymes in the cell wall saccharification and decomposition process. We also report that conditional PCD regulators (LSD1, EDS1 and PAD4) are developmental regulators in woody plants of the secondary cell wall, since lines with the reduced expression of these genes, growing in four subsequent seasons in the field conditions, in stable transgenic aspen significantly reduced lignin content and significantly increased cellulose fiber polymerization degree. Most importantly, these results are strongly supported by significant similarities (homology) between proteins identified in the MALDI-TOF experiment on Arabidopsis thaliana and DEGs identified in the cambium of transgenic aspen Line 4 (Tables S2, S8 and S13).
Usually, the transgenic plants with modified wood composition exhibited impaired growth, especially in the field. Downregulation of cinnamate 4-hydroxylase (an enzyme involved in an early stage of phenylpropanoid synthesis pathway) by RNAi resulted in ~30% less lignin content in hybrid aspen but also reduced tree growth and wood mechanical properties [46]. Another widely used target for genetic manipulation of lignin content is 4-coumarate: coenzyme A ligase (4CL). Suppression of this enzyme leads to a 40–45% reduction in lignin content in aspen [47,48,49], however, this phenotype was confirmed in the greenhouse conditions [50]. During the experiments carried out in the field, the decrease in lignin content was accompanied by stunted growth [51,52,53]. One of the most promising field experiments carried out on poplar with RNAi of 4CL and CCoAOMT resulted in 28% less lignin content and unchanged growth parameters, however, the experiment was performed on 1-year-old plants (one vegetation season) [54]. To sum up, the reduction of lignin is considered often a defense trait against biotic stress. It is not usually studied in the context of plant growth [55,56,57] or reduced-lignin phenotype when growing in the greenhouse. It is also significantly weakened in fluctuating field conditions [17,18,58]. Here, we demonstrate that the deregulation of genes encoding these proteins rather than the lignin biosynthesis pathway itself appears much more promising in biotechnological improvements of trees growing for several subsequent seasons in the field. Most importantly, these regulators control the phenylpropanoid pathway and the genes involved in CWL and lignin biosynthesis (e.g., CAD) in Arabidopsis thaliana and aspen (GO:0045551, Tables S5, S12 and S14). Moreover, it is shown that the deregulation of the conditional PCD regulators results also in the deregulation of many genes from the CAZY family in the cambium of transgenic aspen Line 4 in the field conditions [59,60,61].
Lowering the lignin level up to 22% in transgenic aspen lines reduced their growth only slightly in comparison to the wild type (Figure S12). Additionally, a higher degree of cellulose polymerization (up to 23%) was found in all transgenic lines. The higher cellulose polymerization facilitates the formation of hydrogen bonds among neighboring cellulose fiber and thus plays a pivotal role in increasing the mechanical performance of wood, easier separation of lignin and hemicellulose from cellulose fibers, simultaneously suppressing the decomposition of cellulose during lignin removal [62]. These features were desirable to produce bioethanol from the wood of trees growing in the field for four subsequent seasons (optimal period for harvesting of fast-growing aspen plantation). Indeed, it was confirmed during experiments on a semi-industrial scale. Considering the above, the deregulation of well-known conditional PCD regulators (LSD1, EDS1 and PAD4) [19,44,45,63] gives an unexpected effect in a completely new context i.e., important for the biofuel industry.
Apart from improving wood as a raw technological material, the improvement of the enzymatic hydrolysis and fermentation process could be considered through the identification of the unknown so far plant CWL enzymes. In the industry, the most abundant enzymes are those obtained from the fungus Aspergillus niger. Many of the identified Arabidopsis thaliana proteins taking part in the lysigenous aerenchyma formation and aspen genes involved in the cell wall modifications have many unknown functions, thus they give us a wide range of opportunities to search for new biotechnological innovations for bioethanol, paper, production process or improvement of wood quality.

5. Conclusions

To sum up, we demonstrate for the first time that the deregulation of genes encoding conditional PCD regulators i.e., LSD1, EDS1 and PAD4 can be a breakthrough biotechnology for the biofuel industry. We clearly showed that LSD1 and EDS1 were involved in the plant cell wall saccharification and hypoxia stress and, moreover, we demonstrated that the deregulation in PtLSD1, PtEDS1, and PtPAD4 obtained by genetic enginery results in trees are a better feedstock for biofuel production. Although it needs to be further investigated, some of our discoveries can be applied directly to improve bioethanol production. Our hypothetical model of the role of LSD1, EDS1 and PAD4 in conditional regulation of cell death, stress responses and desired utility feature of wood (level of lignin and cellulose polymerization degree) is presented in Figure 4.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cells12162018/s1, Additional file S1. Supplementary Figure S1. Vector pH7GWIWG2(I) used for transformation of P. tremula × tremuloides. Supplementary Figure S2. cDNA sequences of LSD1, EDS1 and PAD4. Short cDNA sequences inserted into silencing vector are highlighted in yellow, primers are highlighted in grey and underlined. Supplementary Figure S3. Scheme of aspen cultivation and field experiment. Schematic presentation of experimental objects randomization in the field (a). Satellite image of the experimental field (GoogleMaps, USA), the field is marked with a red frame (b). Picture of the field experiment in the second year of the experiment (c). Supplementary Figure S4. Technical details of experiments. (a) Model of experimental setup for lysigenous aerenchyma formation in response to hypoxia stress. (b). Higher reducing sugar content in a lsd1 mutant before and after hypoxia stress was detected. It is easily explainable because the lsd1 mutant induces PCD during the 8-week growth in short photoperiod (8 h of light and 16 h of night) in controlled conditions and cannot survive the hypoxia stress. By contrast, a eds1 mutant increased its reducing sugar content during the hypoxia stress. We did not observe decrease in reducing sugar content during the hypoxia stress in any other genotype. Comparison of a 8-week-old Arabidopsis thaliana plant, Ws-0 (left) and lsd1 (right) (c). The eds1 and eds1/lsd1 mutants did not visually differ from Ws-0. Gel stained with Coomassie Brilliant Blue showing proteins of roots and Supplementary Figure S5 Glucanase and xylanase activities in roots and hypocotyls measured after hypoxia stress. Protein extracts isolated before hypoxia stress exhibited lower activity against carboxymethyl cellulose (a) and xylan (c), compared to extracts isolated after hypoxia stress (b,d). Means values (±SD) are derived from separate biological replicates (n = 3 to 6). Stars above the bars indicate statistically significant differences in comparison to the WT plants, according to the Tukey HSD at level p < 0.05 (*), p < 0.01 (**), and p < 0.001 (***). Supplementary Figure S6. Protein classes analysis by Panther.org. Pie chart presenting different protein classes taking part in hypoxia stress response and their abundance among identified proteins in tested genotypes; Col-0 (a), Ws-0 (b) and eds1 (Ws-0 background) (c) and pie chart presenting division of categories metabolite interconversion enzyme (PC00262) into subcategories (Right site). In order to obtain material for Maldi-TOF analysis, 6 separate experiments lasting one year in total were carried out. Then the Maldi-TOF analysis was performed and Data were based on 6 separate experiments (biological replicates) and two technical repetitions per sample. The results were analyzed using the Mascot software and trimmed so there would be not more than 1% of false positives in each sample. Supplementary Figure S7. Correlation of wood lignin content and deregulation of PtLSD1, PtEDS1 and PtPAD4 gene expression in leafs of tested lines. Correlation of lignin content (percent) and relative expression of LSD1 (a), EDS1 (b) and PAD4 (c). Supplementary Figure S8. RNAseq analysis of RNA isolated from cambium and differentiating xylem of wild type and transgenic Line 4 with the lowest relative lignin content (see Figure 3). The PCA analysis is showing huge differences between RNA samples isolated from cambium and differentiating xylem tissue of wood of wild type and Line 4 (a) and heat map of differently deregulated genes in wild type and Line 4 (b). The analysis was made on three independent biological repeats (n = 3). Supplementary Figure S9. Gene ontology analysis of deregulated genes in Line 4 in terms of their molecular function. Supplementary Figure S10. Gene ontology analysis of deregulated genes in Line 4 in terms of biological processes. Supplementary Figure S11. Results of Line 4 DEGs and its comparison to genes encoding proteins from CAZY families and genes encoding malectin and malectin-like domain-containing proteins. Supplementary Figure S12. Biometric parameters, CO2 assimilation, and chlorophyll content in transgenic lines. Height of the main stem (a), fresh weight of the main stem (b), diameter of the main stem (c), average annual growth of the main stem (c), CO2 assimilation (d), total chlorophyll content in leaves (f). Stars above the bars indicate statistically significant differences in comparison to the WT plants, according to Tukey HSD at level p < 0.05 (*), p < 0.01 (**), and p < 0.001 (***). Mean values (±SD) for five to eleven different biological replications for data on graphs A and B (n = 5 to 11), for three different biological replications on graph C (n = 3), for four different biological replications on graph D (n = 4), for 30 different leaves per genotype on graph E (n = 30) and from nine different technical replications of pooled leaf sample per genotype on graph F (n = 9). Additional file S2. Table S1. List of all genes encoding proteins identified in the Maldi-TOF analysis of 150 kDa protein fraction from Arabidopsis hypocotyls and roots samples of Col-0, Ws-0 and eds 1 mutants collected 7 days after start of root hypoxia stress. Table S2. List of all genes encoding proteins identified in Table S1 using Maldi-TOF analysis of 150 kD protein fraction after removing those present in two or more genotypes (563 proteins). Table S3. List of all genes encoding cell wall loosening proteins and enzymes generating reducing sugars identified in each ecotype (Col-0, Ws-0) and in eds1 mutants. Proteins identified in Maldi-TOF analysis of 150 kDa fraction from hypocotyls and roots samples collected 7 days after start of root hypoxia stress. Gene hierarchy in the Table corresponding to the number of peptides identified for a given protein in the Maldi-TOF analysis of 150 kDa fraction (Table S2). Table S4. List of all genes encoding cell wall loosening (CWL) proteins and enzymes generating reducing sugars. Proteins identified in Maldi-TOF analysis of 150 kDa fraction from Arabidopsis hypocotyls and roots samples collected 7 days after start of root hypoxia stress. Gene hierarchy in the Table corresponding to the number of peptides identified for a given protein in the Maldi-TOF analysis of 150 kDa fraction (Table S3). Table S5. Overrepresentation test of genes encoding proteins identified in Maldi-TOF analysis of 150 kDa fraction of proteins from hypocotyls and roots samples collected 7 days after the start of root hypoxia stress. Table S6. List of genes encoding proteins detected in Maldi-TOF analisis of 150 kDa fraction (Table S2) and not identified by PANTHER.ORG. Table S7. List of all differentially expressed genes (DEGs) compared to wild-type in cambium and differentiating xylem tissues isolated from 4-years-old transgenic aspen Line 4 growing in the field conditions. Table S8. List of differentially expressed genes (DEGs) detected in cambium and differentiating xylem tissues isolated from 4-years-old transgenic aspen Line 4 growing in the field conditions with PANTHER family/subfamily annotation. Table S9. List of differentially expressed genes (DEGs) detected in cambium and differentiating xylem tissues isolated from 4-years-old transgenic aspen Line 4 growing in the field conditions not identified in PANTHER.ORG. Table S10. Enrichment of molecular function analysis of differentially expressed genes (DEGs) detected in cambium and differentiating xylem tissues isolated from 4-years-old transgenic aspen Line 4 growing in the field condition using PANTHER.ORG software. Table S11. Enrichment of biological process analysis of differentially expressed genes (DEGs) detected in cambium and differentiating xylem tissues isolated from 4-years-old transgenic aspen Line 4 growing in the field condition using PNTHER.ORG software. Table S12. The list of differentially expressed genes (DEGs) detected in cambium and differentiating xylem tissues isolated from 4-years-old transgenic aspen Line 4 growing in the field condition involved in cellulose and lignin biosynthesis, cell wall construction and organization and other processes that may have contributed to better utility values. Table S13. List of differentially regulated genes in Line 4 (Table S8) and their corresponding orthologs or homologs identified in the Maldi-TOF analysis of 150 kDa protein fraction from Arabidopsis hypocotyls and roots (Table S2). Table S14. Genes encoding CAZY family proteins in poplar. Table S15. Poplar genes encoding proteins with malectin and malectin-like domains and their predicted characteristics. Table S16. Comparison of DEGs from cambium and differentiating xylem tissues in transgenic aspen Line 4 with genes coding proteins from CAZY families and proteins with malectin and malectin-like domains in poplar.

Author Contributions

S.K. gave idea for all the experiments, formulated the hypothesis and conceived the general research plan; S.K., M.S.-H., J.Z. and W.R.S. supervised the research and experiments; M.J.B. was involved in the aspen field experiments, bioinformatic and other data analysis and writing of the whole manuscript; J.M. carried out the roots hypoxia experiments, determination of reduction sugar content, cellulases activity assays and samples preparation for Maldi-TOF; M.M., A.A., M.D. and J.Z. were involved in the wood chemistry and in the cell wall analysis; D.W. was involved in the hypoxia experiment and biometric measurement of aspen in the field; J.D.-B. was involved in the plasmids preparation for aspen transformation, biometric aspen analysis and preparation of aspen cambium and xylem samples for RNAseq and in vitro aspen propagation; P.G. was involved in the bioinformatic of RNAseq analysis and participated in writing of the manuscript; P.B. was involved in the LSD1, EDS1, and PAD4 mRNA level analysis in all aspen lines; A.R. was responsible for the foliar chlorophyll content determination and in vitro aspen propagation; E.J.M. was involved in supervision of transgenic aspen lines generation in UPSC and the manuscript writing; K.D.-S. performed the alkaline pretreatment and enzymatic hydrolysis of aspen and ethanol fermentation and its efficiency analysis; M.J.B., J.M., M.S.-H., E.J.M., J.Z., W.R.S. and S.K. interpreted results and wrote the very manuscript with participation of A.A. and M.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported partially privately by Stanisław Karpiński and also in part by the “Biostrateg 2” project (BIOSTRATEG2/298241/10/NCBR/2016) granted to Stanisław Karpiński by the National Center for Research and Development and by the “OPUS 15” project (UMO-2018/29/B/NZ3/01198) granted to Stanisław Karpiński by the National Science Centre.

Institutional Review Board Statement

Does not apply.

Data Availability Statement

All data generated or analyzed during this study are included in the very published article and its Supplementary files.

Acknowledgments

We would like to thank Alex Białas for his technical help in dendrometry measurements. Also, we would like to express our gratitude to Aleksandra Lechańska for proofreading.

Conflicts of Interest

The authors declare that they have no competing interest.

Abbreviations

LSD1LESION SIMULATING DISEASE 1
PAD4PHYTOALEXIN DEFICIENT 4
EDS1ENHANCED DISEASE SUSCEPTIBILITY
PCDProgramed Cell Death
CMCCarboxymethyl cellulose
DEGDifferential expressed genes
4CL4-coumarate:coenzyme A ligase
CWLcell wall loosening
DEGdifferentially expressed genes
MSMurashige and Skoog medium

References

  1. Hoengenaert, L.; Wouters, M.; Kim, H.; De Meester, B.; Morreel, K.; Vandersyppe, S.; Pollier, J.; Desmet, S.; Goeminne, G.; Ralph, J.; et al. Overexpression of the Scopoletin Biosynthetic Pathway Enhances Lignocellulosic Biomass Processing. Sci. Adv. 2022, 8, eabo5738. [Google Scholar] [CrossRef]
  2. Zoghlami, A.; Paës, G. Lignocellulosic Biomass: Understanding Recalcitrance and Predicting Hydrolysis. Front. Chem. 2019, 7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Shukla, A.; Kumar, D.; Girdhar, M.; Kumar, A.; Goyal, A.; Malik, T.; Mohan, A. Strategies of Pretreatment of Feedstocks for Optimized Bioethanol Production: Distinct and Integrated Approaches. Biotechnol. Biofuels Bioprod. 2023, 16, 44. [Google Scholar] [CrossRef]
  4. Czekała, W.; Pulka, J.; Jasiński, T.; Szewczyk, P.; Bojarski, W.; Jasiński, J. Waste as Substrates for Agricultural Biogas Plants: A Case Study from Poland. J. Water Land Dev. 2023, 45–50. [Google Scholar] [CrossRef]
  5. Turumtay, H. Cell Wall Engineering by Heterologous Expression of Cell Wall-Degrading Enzymes for Better Conversion of Lignocellulosic Biomass into Biofuels. Bioenerg. Res. 2015, 8, 1574–1588. [Google Scholar] [CrossRef]
  6. Khan, I.; Akhtar, M.W. BioenergyProduction from Plant Biomass: Bioethanol from Concept to Reality. Nat. Prec. 2011, 1. [Google Scholar] [CrossRef] [Green Version]
  7. Alonso, D.M.; Hakim, S.H.; Zhou, S.; Won, W.; Hosseinaei, O.; Tao, J.; Garcia-Negron, V.; Motagamwala, A.H.; Mellmer, M.A.; Huang, K.; et al. Increasing the Revenue from Lignocellulosic Biomass: Maximizing Feedstock Utilization. Sci. Adv. 2017, 3, e1603301. [Google Scholar] [CrossRef] [Green Version]
  8. Sidana, A.; Yadav, S.K. Recent Developments in Lignocellulosic Biomass Pretreatment with a Focus on Eco-Friendly, Non-Conventional Methods. J. Clean. Prod. 2022, 335, 130286. [Google Scholar] [CrossRef]
  9. Antczak, A.; Szadkowski, J.; Szadkowska, D.; Zawadzki, J. Assessment of the Effectiveness of Liquid Hot Water and Steam Explosion Pretreatments of Fast-Growing Poplar (Populus trichocarpa) Wood. Wood Sci. Technol. 2022, 56, 87–109. [Google Scholar] [CrossRef]
  10. Lu, X.; Zheng, X.; Li, X.; Zhao, J. Adsorption and Mechanism of Cellulase Enzymes onto Lignin Isolated from Corn Stover Pretreated with Liquid Hot Water. Biotechnol. Biofuels 2016, 9, 118. [Google Scholar] [CrossRef] [Green Version]
  11. Fletcher, E.; Baetz, K. Multi-Faceted Systems Biology Approaches Present a Cellular Landscape of Phenolic Compound Inhibition in Saccharomyces Cerevisiae. Front. Bioeng. Biotechnol. 2020, 8, 539902. [Google Scholar] [CrossRef] [PubMed]
  12. Li, Y.; Fu, Q.; Rojas, R.; Yan, M.; Lawoko, M.; Berglund, L. Lignin-Retaining Transparent Wood. ChemSusChem 2017, 10, 3445–3451. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Escamez, S.; Tuominen, H. Contribution of Cellular Autolysis to Tissular Functions during Plant Development. Curr. Opin. Plant Biol. 2017, 35, 124–130. [Google Scholar] [CrossRef] [PubMed]
  14. Mühlenbock, P.; Plaszczyca, M.; Plaszczyca, M.; Mellerowicz, E.; Karpinski, S. Lysigenous Aerenchyma Formation in Arabidopsis Is Controlled by LESION SIMULATING DISEASE1. Plant Cell 2007, 19, 3819–3830. [Google Scholar] [CrossRef] [Green Version]
  15. Taylor, G. Populus: Arabidopsis for Forestry. Do We Need a Model Tree? Ann. Bot. 2002, 90, 681–689. [Google Scholar] [CrossRef] [Green Version]
  16. Chaffey, N.; Cholewa, E.; Regan, S.; Sundberg, B. Secondary Xylem Development in Arabidopsis: A Model for Wood Formation. Physiol. Plant 2002, 114, 594–600. [Google Scholar] [CrossRef]
  17. Wituszyńska, W.; Ślesak, I.; Vanderauwera, S.; Szechyńska-Hebda, M.; Kornaś, A.; Kelen, K.V.D.; Mühlenbock, P.; Karpińska, B.; Maćkowski, S.; Breusegem, F.V.; et al. LESION SIMULATING DISEASE1, ENHANCED DISEASE SUSCEPTIBILITY1, and PHYTOALEXIN DEFICIENT4 Conditionally Regulate Cellular Signaling Homeostasis, Photosynthesis, Water Use Efficiency, and Seed Yield in Arabidopsis. Plant Physiol. 2013, 161, 1795–1805. [Google Scholar] [CrossRef] [Green Version]
  18. Bernacki, M.J.; Czarnocka, W.; Rusaczonek, A.; Witoń, D.; Kęska, S.; Czyż, J.; Szechyńska-Hebda, M.; Karpiński, S. LSD1, EDS1 and PAD4-Dependent Conditional Correlation among Salicylic Acid, Hydrogen Peroxide, Water Use Efficiency, and Seed Yield in Arabidopsis Thaliana. Physiol. Plant. 2018, 165, 369–382. [Google Scholar] [CrossRef]
  19. Czarnocka, W.; Van Der Kelen, K.; Willems, P.; Szechyńska-Hebda, M.; Shahnejat-Bushehri, S.; Balazadeh, S.; Rusaczonek, A.; Mueller-Roeber, B.; Van Breusegem, F.; Karpiński, S. The Dual Role of LESION SIMULATING DISEASE 1 as a Condition-Dependent Scaffold Protein and Transcription Regulator. Plant Cell Environ. 2017, 40, 2644–2662. [Google Scholar] [CrossRef] [Green Version]
  20. Karimi, M.; Inzé, D.; Depicker, A. GATEWAYTM Vectors for Agrobacterium-Mediated Plant Transformation. Trends Plant Sci. 2002, 7, 193–195. [Google Scholar] [CrossRef]
  21. Nilsson, O.; Aldén, T.; Sitbon, F.; Little, C.H.A.; Chalupa, V.; Sandberg, G.; Olsson, O. Spatial Pattern of Cauliflower Mosaic Virus 35S Promoter-Luciferase Expression in Transgenic Hybrid Aspen Trees Monitored by Enzymatic Assay and Non-Destructive Imaging. Transgenic Res. 1992, 1, 209–220. [Google Scholar] [CrossRef]
  22. Murashige, T.; Skoog, F. A Revised Medium for Rapid Growth and Bio Assays with Tobacco Tissue Cultures. Physiol. Plant. 1962, 15, 473–497. [Google Scholar] [CrossRef]
  23. Bernacki, M.J.; Czarnocka, W.; Witoń, D.; Rusaczonek, A.; Szechyńska-Hebda, M.; Ślesak, I.; Dąbrowska-Bronk, J.; Karpiński, S. ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1) Affects Development, Photosynthesis, and Hormonal Homeostasis in Hybrid Aspen (Populus tremula L. × P. tremuloides). J. Plant Physiol. 2018, 226, 91–102. [Google Scholar] [CrossRef]
  24. Tuskan, G.A.; DiFazio, S.; Jansson, S.; Bohlmann, J.; Grigoriev, I.; Hellsten, U.; Putnam, N.; Ralph, S.; Rombauts, S.; Salamov, A.; et al. The Genome of Black Cottonwood, Populus Trichocarpa (Torr. & Gray). Science 2006, 313, 1596–1604. [Google Scholar] [CrossRef] [Green Version]
  25. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast Universal RNA-Seq Aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  26. Love, M.I.; Huber, W.; Anders, S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Rusaczonek, A.; Czarnocka, W.; Kacprzak, S.; Witoń, D.; Ślesak, I.; Szechyńska-Hebda, M.; Gawroński, P.; Karpiński, S. Role of Phytochromes A and B in the Regulation of Cell Death and Acclimatory Responses to UV Stress in Arabidopsis Thaliana. J. Exp. Bot. 2015, 66, erv375. [Google Scholar] [CrossRef] [Green Version]
  28. Wise, L.E.; Murphy, M.; Addieco, A.A.d’. Chlorite Holocellulose, Its Fractionation and Bearing on Summative Wood Analysis and Studies on the Hemicelluloses. Pap. Trade J. 1946, 122, 35–43. [Google Scholar]
  29. Sluiter, A. Determination of Extractives in Biomass: Laboratory Analytical Procedure (LAP); Issue Date 7/17/2005. Tech. Rep. 2008, 12. [Google Scholar]
  30. Antczak, A.; Marchwicka, M.; Szadkowski, J.; Drożdżek, M.; Gawron, J.; Radomski, A.; Zawadzki, J. Sugars Yield Obtained after Acid and Enzymatic Hydrolysis of Fast-Growing Poplar Wood Species. BioResources 2018, 13, 8629–8645. [Google Scholar] [CrossRef]
  31. Adney, B.; Baker, J. Measurement of Cellulase Activities: Laboratory Analytical Procedure (LAP); Issue Date: 08/12/1996. Tech. Rep. 2008, 11. [Google Scholar]
  32. Mateo, A.; Mühlenbock, P.; Rustérucci, C.; Chang, C.C.-C.; Miszalski, Z.; Karpinska, B.; Parker, J.E.; Mullineaux, P.M.; Karpinski, S. LESION SIMULATING DISEASE 1 Is Required for Acclimation to Conditions That Promote Excess Excitation Energy. Plant Physiol. 2004, 136, 2818–2830. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Feys, B.J.; Moisan, L.J.; Newman, M.-A.; Parker, J.E. Direct Interaction between the Arabidopsis Disease Resistance Signaling Proteins, EDS1 and PAD4. EMBO J. 2001, 20, 5400–5411. [Google Scholar] [CrossRef] [PubMed]
  34. Ślesak, I.; Szechyńska-Hebda, M.; Fedak, H.; Sidoruk, N.; Dąbrowska-Bronk, J.; Witoń, D.; Rusaczonek, A.; Antczak, A.; Drożdżek, M.; Karpińska, B.; et al. PHYTOALEXIN DEFICIENT 4 Affects Reactive Oxygen Species Metabolism, Cell Wall and Wood Properties in Hybrid Aspen (Populus tremula L.  × Tremuloides). Plant Cell Environ. 2015, 38, 1275–1284. [Google Scholar] [CrossRef] [PubMed]
  35. Szechyńska-Hebda, M.; Czarnocka, W.; Hebda, M.; Karpiński, S. PAD4, LSD1 and EDS1 Regulate Drought Tolerance, Plant Biomass Production, and Cell Wall Properties. Plant Cell Rep. 2016, 35, 527–539. [Google Scholar] [CrossRef]
  36. Pinto, A.S.S.; Brondi, M.G.; de Freitas, J.V.; Furlan, F.F.; Ribeiro, M.P.A.; Giordano, R.C.; Farinas, C.S. Mitigating the Negative Impact of Soluble and Insoluble Lignin in Biorefineries. Renew. Energy 2021, 173, 1017–1026. [Google Scholar] [CrossRef]
  37. Douglas, C.J. Phenylpropanoid Metabolism and Lignin Biosynthesis: From Weeds to Trees. Trends Plant Sci. 1996, 1, 171–178. [Google Scholar] [CrossRef]
  38. Kumar, V.; Hainaut, M.; Delhomme, N.; Mannapperuma, C.; Immerzeel, P.; Street, N.R.; Henrissat, B.; Mellerowicz, E.J. Poplar Carbohydrate-Active Enzymes: Whole-Genome Annotation and Functional Analyses Based on RNA Expression Data. Plant J. 2019, 99, 589–609. [Google Scholar] [CrossRef] [Green Version]
  39. Kumar, V.; Donev, E.N.; Barbut, F.R.; Kushwah, S.; Mannapperuma, C.; Urbancsok, J.; Mellerowicz, E.J. Genome-Wide Identification of Populus Malectin/Malectin-Like Domain-Containing Proteins and Expression Analyses Reveal Novel Candidates for Signaling and Regulation of Wood Development. Front. Plant Sci. 2020, 11, 588846. [Google Scholar] [CrossRef]
  40. Taherzadeh, M.J.; Karimi, K. Pretreatment of Lignocellulosic Wastes to Improve Ethanol and Biogas Production: A Review. Int. J. Mol. Sci. 2008, 9, 1621–1651. [Google Scholar] [CrossRef] [Green Version]
  41. Zhang, J.; Liu, Y.; Li, C.; Yin, B.; Liu, X.; Guo, X.; Zhang, C.; Liu, D.; Hwang, I.; Li, H.; et al. PtomtAPX Is an Autonomous Lignification Peroxidase during the Earliest Stage of Secondary Wall Formation in Populus Tomentosa Carr. Nat. Plants 2022, 8, 828–839. [Google Scholar] [CrossRef] [PubMed]
  42. Taylor-Teeples, M.; Lin, L.; de Lucas, M.; Turco, G.; Toal, T.W.; Gaudinier, A.; Young, N.F.; Trabucco, G.M.; Veling, M.T.; Lamothe, R.; et al. An Arabidopsis Gene Regulatory Network for Secondary Cell Wall Synthesis. Nature 2015, 517, 571–575. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Schuetz, M.; Smith, R.; Ellis, B. Xylem Tissue Specification, Patterning, and Differentiation Mechanisms. J. Exp. Bot. 2013, 64, 11–31. [Google Scholar] [CrossRef]
  44. Karpiński, S.; Szechyńska-Hebda, M.; Wituszyńska, W.; Burdiak, P. Light Acclimation, Retrograde Signalling, Cell Death and Immune Defences in Plants. Plant Cell Environ. 2013, 36, 736–744. [Google Scholar] [CrossRef]
  45. Mühlenbock, P.; Szechynska-Hebda, M.; Plaszczyca, M.; Baudo, M.; Mateo, A.; Mullineaux, P.M.; Parker, J.E.; Karpinska, B.; Karpinski, S. Chloroplast Signaling and LESION SIMULATING DISEASE1 Regulate Crosstalk between Light Acclimation and Immunity in Arabidopsis. Plant Cell 2008, 20, 2339–2356. [Google Scholar] [CrossRef] [Green Version]
  46. Bjurhager, I.; Olsson, A.-M.; Zhang, B.; Gerber, L.; Kumar, M.; Berglund, L.A.; Burgert, I.; Sundberg, B.; Salmén, L. Ultrastructure and Mechanical Properties of Populus Wood with Reduced Lignin Content Caused by Transgenic Down-Regulation of Cinnamate 4-Hydroxylase. Biomacromolecules 2010, 11, 2359–2365. [Google Scholar] [CrossRef] [PubMed]
  47. Hu, W.-J.; Harding, S.A.; Lung, J.; Popko, J.L.; Ralph, J.; Stokke, D.D.; Tsai, C.-J.; Chiang, V.L. Repression of Lignin Biosynthesis Promotes Cellulose Accumulation and Growth in Transgenic Trees. Nat. Biotechnol. 1999, 17, 808–812. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Jia, C.; Zhao, H.; Wang, H.; Xing, Z.; Du, K.; Song, Y.; Wei, J. Obtaining the Transgenic Poplars with Low Lignin Content through Down-Regulation Of4CL. Chin. Sci. Bull. 2004, 49, 905–909. [Google Scholar] [CrossRef]
  49. Li, L.; Zhou, Y.; Cheng, X.; Sun, J.; Marita, J.M.; Ralph, J.; Chiang, V.L. Combinatorial Modification of Multiple Lignin Traits in Trees through Multigene Cotransformation. Proc. Natl. Acad. Sci. USA 2003, 100, 4939–4944. [Google Scholar] [CrossRef]
  50. Sen, S.K.; Min, D.; Savithri, D.; Lu, F.; Jameel, H.; Chiang, V.; Chang, H. Field-Grown Transgenic Hybrid Poplar with Modified Lignin Biosynthesis to Improve Enzymatic Saccharification Efficiency. ACS Sustain. Chem. Eng. 2017, 5, 2407–2414. [Google Scholar] [CrossRef]
  51. Stout, A.T.; Davis, A.A.; Domec, J.-C.; Yang, C.; Shi, R.; King, J.S. Growth under Field Conditions Affects Lignin Content and Productivity in Transgenic Populus Trichocarpa with Altered Lignin Biosynthesis. Biomass Bioenergy 2014, 68, 228–239. [Google Scholar] [CrossRef] [Green Version]
  52. Leplé, J.-C.; Dauwe, R.; Morreel, K.; Storme, V.; Lapierre, C.; Pollet, B.; Naumann, A.; Kang, K.-Y.; Kim, H.; Ruel, K.; et al. Downregulation of Cinnamoyl-Coenzyme A Reductase in Poplar: Multiple-Level Phenotyping Reveals Effects on Cell Wall Polymer Metabolism and Structure. Plant Cell 2007, 19, 3669–3691. [Google Scholar] [CrossRef] [Green Version]
  53. Van Acker, R.; Leplé, J.-C.; Aerts, D.; Storme, V.; Goeminne, G.; Ivens, B.; Légée, F.; Lapierre, C.; Piens, K.; Van Montagu, M.C.E.; et al. Improved Saccharification and Ethanol Yield from Field-Grown Transgenic Poplar Deficient in Cinnamoyl-CoA Reductase. Proc. Natl. Acad. Sci. USA 2014, 111, 845–850. [Google Scholar] [CrossRef]
  54. Tian, X.; Xie, J.; Zhao, Y.; Lu, H.; Liu, S.; Qu, L.; Li, J.; Gai, Y.; Jiang, X. Sense-, Antisense- and RNAi-4CL1 Regulate Soluble Phenolic Acids, Cell Wall Components and Growth in Transgenic Populus Tomentosa Carr. Plant Physiol. Biochem. 2013, 65, 111–119. [Google Scholar] [CrossRef]
  55. Min, D.; Li, Q.; Jameel, H.; Chiang, V.; Chang, H. The Cellulase-Mediated Saccharification on Wood Derived from Transgenic Low-Lignin Lines of Black Cottonwood (Populus trichocarpa). Appl Biochem Biotechnol 2012, 168, 947–955. [Google Scholar] [CrossRef]
  56. Min, D.; Yang, C.; Chiang, V.; Jameel, H.; Chang, H. The Influence of Lignin–Carbohydrate Complexes on the Cellulase-Mediated Saccharification II: Transgenic Hybrid Poplars (Populus nigra L. and Populus maximowiczii A.). Fuel 2014, 116, 56–62. [Google Scholar] [CrossRef]
  57. Coleman, H.D.; Park, J.-Y.; Nair, R.; Chapple, C.; Mansfield, S.D. RNAi-Mediated Suppression of p-Coumaroyl-CoA 3’-Hydroxylase in Hybrid Poplar Impacts Lignin Deposition and Soluble Secondary Metabolism. Proc. Natl. Acad. Sci. USA 2008, 105, 4501–4506. [Google Scholar] [CrossRef]
  58. Wituszyńska, W.; Gałązka, K.; Rusaczonek, A.; Vanderauwera, S.; Van Breusegem, F.; Karpiński, S. Multivariable Environmental Conditions Promote Photosynthetic Adaptation Potential in Arabidopsis Thaliana. J. Plant Physiol. 2013, 170, 548–559. [Google Scholar] [CrossRef]
  59. Gao, F.; Shu, X.; Ali, M.B.; Howard, S.; Li, N.; Winterhagen, P.; Qiu, W.; Gassmann, W. A Functional EDS1 Ortholog Is Differentially Regulated in Powdery Mildew Resistant and Susceptible Grapevines and Complements an Arabidopsis Eds1 Mutant. Planta 2010, 231, 1037–1047. [Google Scholar] [CrossRef]
  60. Guo, J.; Bai, P.; Yang, Q.; Liu, F.; Wang, X.; Huang, L.; Kang, Z. Wheat Zinc Finger Protein TaLSD1, a Negative Regulator of Programmed Cell Death, Is Involved in Wheat Resistance against Stripe Rust Fungus. Plant Physiol. Biochem. 2013, 71, 164–172. [Google Scholar] [CrossRef]
  61. Guan, Q.; Ma, H.; Wang, Z.; Wang, Z.; Bu, Q.; Liu, S. A Rice LSD1-like-Type ZFP Gene OsLOL5 Enhances Saline-Alkaline Tolerance in Transgenic Arabidopsis Thaliana, Yeast and Rice. BMC Genom. 2016, 17, 142. [Google Scholar] [CrossRef] [Green Version]
  62. Fang, Z.; Li, B.; Liu, Y.; Zhu, J.; Li, G.; Hou, G.; Zhou, J.; Qiu, X. Critical Role of Degree of Polymerization of Cellulose in Super-Strong Nanocellulose Films. Matter 2020, 2, 1000–1014. [Google Scholar] [CrossRef]
  63. Dietrich, R.A.; Richberg, M.H.; Schmidt, R.; Dean, C.; Dangl, J.L. A Novel Zinc Finger Protein Is Encoded by the Arabidopsis LSD1 Gene and Functions as a Negative Regulator of Plant Cell Death. Cell 1997, 88, 685–694. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Lysigenous aerenchyma formation in Arabidopsis thaliana roots and hypocotyls in response to roots hypoxia stress monitored by reducing sugar content and in vitro activity of cellulolytic and xylanolytic enzymes. The reducing sugar content was measured before (A) and after (B) 7 days of hypoxia stress expressed as D-glucose content. Analysis of the hydrolytic activity of protein extracts against carboxymethyl cellulose (C) and xylan (D) was measured in controlled conditions and after 7-day-hypoxia stress which caused the release of reducing sugars during 3-h incubation. Means values (±SD) were calculated from a minimum of three and a maximum of six separate biological replicates (n = 3 to 6). Stars above the bars indicate statistically significant differences in comparison to the Ws-0 plants, according to Tukey HSD at levels p < 0.05 (*), p < 0.01 (**), and p < 0.001 (***).
Figure 1. Lysigenous aerenchyma formation in Arabidopsis thaliana roots and hypocotyls in response to roots hypoxia stress monitored by reducing sugar content and in vitro activity of cellulolytic and xylanolytic enzymes. The reducing sugar content was measured before (A) and after (B) 7 days of hypoxia stress expressed as D-glucose content. Analysis of the hydrolytic activity of protein extracts against carboxymethyl cellulose (C) and xylan (D) was measured in controlled conditions and after 7-day-hypoxia stress which caused the release of reducing sugars during 3-h incubation. Means values (±SD) were calculated from a minimum of three and a maximum of six separate biological replicates (n = 3 to 6). Stars above the bars indicate statistically significant differences in comparison to the Ws-0 plants, according to Tukey HSD at levels p < 0.05 (*), p < 0.01 (**), and p < 0.001 (***).
Cells 12 02018 g001
Figure 2. Identification of proteins isolated from roots and hypocotyls during lysogenic aerenchyma formation due to root hypoxia stress by the Maldi-TOF analysis. A number of enzymes generating reducing ends (A) and their presence in particular genotypes (B). Protein classes analysis by the Panther is shown as a pie chart presenting different protein classes taking part in root hypoxia stress response and their abundance among identified proteins in eds1 (Col-0) and a pie chart presenting the division of categories metabolite interconversion enzyme (PC00262) into subcategories (C). EDS1 abundance dependent generating reducing ends enzymes (D) Data are based on six separate experiments (biological replicates) and two technical repetitions per sample. The results were analyzed using the Mascot software and trimmed so there will be no more than 1% of false positives in each sample.
Figure 2. Identification of proteins isolated from roots and hypocotyls during lysogenic aerenchyma formation due to root hypoxia stress by the Maldi-TOF analysis. A number of enzymes generating reducing ends (A) and their presence in particular genotypes (B). Protein classes analysis by the Panther is shown as a pie chart presenting different protein classes taking part in root hypoxia stress response and their abundance among identified proteins in eds1 (Col-0) and a pie chart presenting the division of categories metabolite interconversion enzyme (PC00262) into subcategories (C). EDS1 abundance dependent generating reducing ends enzymes (D) Data are based on six separate experiments (biological replicates) and two technical repetitions per sample. The results were analyzed using the Mascot software and trimmed so there will be no more than 1% of false positives in each sample.
Cells 12 02018 g002
Figure 3. Integrated deregulation of PtLSD1, PtEDS1, and PtPAD4 affected lignin content and other wood chemical properties. Transgenic aspen trees (Lines 1–4) grew for four subsequent seasons in field conditions. The relative expression level of PtLSD1, PtEDS1 and PtPAD4 in leaves (AC), relative lignin, cellulose, holocellulose content and a degree of cellulose polymerization (D), glucose content after saccharification of wood (E), the amount of ethanol in the fermentation process (F) and the amount of un-fermented glucose in the pulp after fermentation (G). Mean values (±SD) are derived for three different biological and three technical replications for each transgenic Line on graphs (AC) (n = 9) and for three to six separate biological replications (n = 3 to 6) for data (DG). Stars above the bars indicate statistically significant differences in comparison to the WT plants, according to Tukey HSD at levels p < 0.05 (*), p < 0.01 (**), and p < 0.001 (***).
Figure 3. Integrated deregulation of PtLSD1, PtEDS1, and PtPAD4 affected lignin content and other wood chemical properties. Transgenic aspen trees (Lines 1–4) grew for four subsequent seasons in field conditions. The relative expression level of PtLSD1, PtEDS1 and PtPAD4 in leaves (AC), relative lignin, cellulose, holocellulose content and a degree of cellulose polymerization (D), glucose content after saccharification of wood (E), the amount of ethanol in the fermentation process (F) and the amount of un-fermented glucose in the pulp after fermentation (G). Mean values (±SD) are derived for three different biological and three technical replications for each transgenic Line on graphs (AC) (n = 9) and for three to six separate biological replications (n = 3 to 6) for data (DG). Stars above the bars indicate statistically significant differences in comparison to the WT plants, according to Tukey HSD at levels p < 0.05 (*), p < 0.01 (**), and p < 0.001 (***).
Cells 12 02018 g003
Figure 4. Hypothetical model of the role of LSD1, EDS1 and PAD4 in conditional regulation of the cell wall properties. Based on our results we propose the following model of the role of LSD1/EDS1/PAD4 trimer in conditional regulation of lignin biosynthesis and cellulose polymerization degree during plant growth in a natural environment. Both in Arabidopsis hypocotyls undergoing lisogenous aerenchyma formation proteomic analysis and from xylem/floem transcriptome analysis in transgenic Line 4, many proteins/genes from the phenylpropanoid pathway were deregulated (Tables S5 and S11). Knowing that this is one of the most important pathways in the lignin biosynthesis process we proposed this as the main reason for lower lignin content in the cell wall of transgenic aspen lines. The molecular/biochemical pathway of cellulose polymerization degree could be regulated by other unknown regulators (11%) as a response to transgenic deregulation of LSD1, EDS1 and PAD4. We hypothesize that a higher cellulose polymerization degree in transgenic Line 4 is a form of compensation for reduced lignin content.
Figure 4. Hypothetical model of the role of LSD1, EDS1 and PAD4 in conditional regulation of the cell wall properties. Based on our results we propose the following model of the role of LSD1/EDS1/PAD4 trimer in conditional regulation of lignin biosynthesis and cellulose polymerization degree during plant growth in a natural environment. Both in Arabidopsis hypocotyls undergoing lisogenous aerenchyma formation proteomic analysis and from xylem/floem transcriptome analysis in transgenic Line 4, many proteins/genes from the phenylpropanoid pathway were deregulated (Tables S5 and S11). Knowing that this is one of the most important pathways in the lignin biosynthesis process we proposed this as the main reason for lower lignin content in the cell wall of transgenic aspen lines. The molecular/biochemical pathway of cellulose polymerization degree could be regulated by other unknown regulators (11%) as a response to transgenic deregulation of LSD1, EDS1 and PAD4. We hypothesize that a higher cellulose polymerization degree in transgenic Line 4 is a form of compensation for reduced lignin content.
Cells 12 02018 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bernacki, M.J.; Mielecki, J.; Antczak, A.; Drożdżek, M.; Witoń, D.; Dąbrowska-Bronk, J.; Gawroński, P.; Burdiak, P.; Marchwicka, M.; Rusaczonek, A.; et al. Biotechnological Potential of the Stress Response and Plant Cell Death Regulators Proteins in the Biofuel Industry. Cells 2023, 12, 2018. https://doi.org/10.3390/cells12162018

AMA Style

Bernacki MJ, Mielecki J, Antczak A, Drożdżek M, Witoń D, Dąbrowska-Bronk J, Gawroński P, Burdiak P, Marchwicka M, Rusaczonek A, et al. Biotechnological Potential of the Stress Response and Plant Cell Death Regulators Proteins in the Biofuel Industry. Cells. 2023; 12(16):2018. https://doi.org/10.3390/cells12162018

Chicago/Turabian Style

Bernacki, Maciej Jerzy, Jakub Mielecki, Andrzej Antczak, Michał Drożdżek, Damian Witoń, Joanna Dąbrowska-Bronk, Piotr Gawroński, Paweł Burdiak, Monika Marchwicka, Anna Rusaczonek, and et al. 2023. "Biotechnological Potential of the Stress Response and Plant Cell Death Regulators Proteins in the Biofuel Industry" Cells 12, no. 16: 2018. https://doi.org/10.3390/cells12162018

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop